CTLA4 tagging polymorphisms and risk of colorectal cancer: a case–control study involving 2,306 subjects

Abstract

Background: CTLA4 is a candidate gene which has been implicated in the development of colorectal cancer (CRC).

Patients and Methods: To determine the important role of CTLA-4 polymorphisms on risk of CRC, we genotyped four CTLA-4 tagging polymorphisms and calculated crude/adjusted ORs with their 95% CIs. We recruited 1,003 sporadic CRC cases and 1,303 controls.

Results: The findings suggested that CTLA-4 rs231775 G>A polymorphism increased the risk of CRC (homozygote model: adjusted OR=1.40, 95% CI=1.05–1.87,P=0.022; dominant model: adjusted OR=1.19, 95% CI=1.00–1.41,P=0.047; and recessive model: adjusted OR=1.38, 95% CI=1.05–1.82, P=0.021). In a stratified analysis by site of tumor, this association was also found in colon cancer. We also found that CTLA-4 rs231775 GA/AA genotypes might be associated with an increased risk of CRC in Zhenjiang cohort. In addition, we found the CTLA-4 rs16840252 C>T polymorphism was associated with the risk of colon cancer. Haplotype comparison analysis showed that CTLA-4 G_rs3087243C_rs16840252C_rs733618 A_rs231775, G_rs3087243C_rs16840252T_rs733618A_rs231775, and other haplotypes increased the risk of CRC (P<0.001, <0.001, and 0.002, respectively).

Conclusion: This study evidences an association of CTLA-4 tagging polymorphisms and haplotypes with CRC risk. Additional well-designed studies with large sample sizes are required to confirm our findings.

Introduction

Colorectal cancer (CRC) is one of the most common malignancies and is becoming the fifth leading cause of cancer death in China.¹ It is reported that the incidence rate and CRC-related mortality are increasing rapidly worldwide.^2,3 These may be attributed to certain lifestyles and environmental factors, including physically inactive, overweight, smoking, and drinking. Recently, Katsidzira et al⁴ reported that diabetes mellitus, previous schistosomiasis, approximation to a western lifestyle, and family history were the predominant associations with CRC. Besides these unhealthy lifestyle and environmental risk factors, genetic factors may also affect the development of CRC. A previous study demonstrated that genetic risk factors may contribute to ~35% etiology of CRC cases.⁵ Up to now, the inherited factor of CRC remains controversial. Recently, a number of investigations have been devoted to exploring the potential molecular mechanism of CRC carcinogenesis, and inherited factors have been considered to play a vital role in the occurrence and development of CRC.^6,7 Investigations of these conceivable inherited factors may enrich our view on the etiology of CRC.

The cytotoxic T-lymphocyte antigen-4 (CTLA-4) gene, also named as a cluster of differentiation 152 (CD152), is one of the costimulatory molecule genes involved in immune response. CTLA-4 is transiently expressed on some activated T cells.⁸ The expression of CTLA-4 inhibits cytokine production of T cells and then provides a negative signal to T cells.^9,10 The structure of CTLA-4 shares some homologies with CD28 and binds to B7.1 and B7.2 ligands competitively. However, CTLA-4 has a higher binding affinity with B7 molecules compared with CD28. Through interaction of CTLA-4 with B7 molecules, T-cell proliferation, activation, and cytokine production can be inhibited.^11–13 The CTLA-4 gene is located in chromosome 2q33, which belongs to several immune regulatory gene regions. Since CTLA-4 acts as a vital regulatory factor for some immune responses, any genetic variation in CTLA-4 gene may influence normal immune function and then alter the risk of cancer. Hence, exploring the impact of these genetic variations in CTLA-4 gene could determine their relationship with cancer susceptibility. CTLA-4 is polymorphic and contains more than 100 single-nucleotide polymorphisms (SNPs). Among them, some SNPs (eg, rs3087243 G>A, rs16840252 C>T, rs4553808 T>C, rs5742909 C>T, rs733618 T>C, and rs231775 G>A polymorphisms) in CTLA-4 gene were extensively studied and were reported to be correlated with risk of human malignancy.^14,15 Some case–control studies explored the relationship between CTLA-4 polymorphisms and CRC;^16,17 however, the sample sizes were limited and the results remained conflicting.

The evasion of immune surveillance and the production of immunosuppressive cytokines are two of the most common immune defects identified to be correlated with CRC. CTLA-4 is a candidate gene which has been implicated in immune response.¹⁸ Moreover, due to the emerging role of CTLA-4 as an immune checkpoint molecule, anti-CTLA-4 antibody has been tested recently in the treatment of CRC patients.¹⁹ Previous case–control studies, conducted in diverse population to assess the relationship of CRC with CTLA-4 polymorphisms, have generated conflicting findings. Hence, we undertook a study to determine whether CTLA-4 variations could cause a predisposition toward CRC. In this study, we analyzed the tagging SNPs of CTLA-4 (rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A polymorphisms) and identified whether these SNPs confer susceptibility to CRC in an Eastern Chinese Han population.

Patients and methods

Study population and patient selection

The study population consisted of 2,306 subjects (1,003 diagnosed as sporadic CRC and 1,303 controls) between October 2014 and August 2017. Among them, 757 CRC patients and 680 controls were enrolled from the Affiliated Union Hospital of Fujian Medical University (Fuzhou, China), and 246 CRC patients and 623 controls were enrolled from the Affiliated People's Hospital of Jiangsu University (Zhenjiang, China). CRC was confirmed via pathology. The age of CRC cases ranged from 21 to 90 years old (mean age at diagnosis was 61.10±12.17 years). The major exclusion criteria were autoimmune disorders, hereditary nonpolyposis CRC, and history of another malignancy. CRC cases who had received neoadjuvant chemoradiotherapy were also excluded. The age of controls ranged from 21 to 87 years old (mean age at sampling was 61.40±9.61 years). In this study, the cancer-free controls included 1,303 healthy volunteers who participated in a routine examination in hospitals mentioned above. The primary information of the participants was collected by a pre-structured questionnaire. The definitions of "ever smokers" and "ever drinkers" are described in our previous study.²⁰ In addition, according to the criterion for overweight and obesity, a body mass index (BMI) of 24 was used as the cutoff point in Chinese adults.^21,22 Each participant was informed about the present study and signed a standard informed consent form. The Ethical Committee of Fujian Medical University and Jiangsu University approved the protocols of the study (No KY-2013–11 and No 2012-00-18, respectively).

Data collection

All participants were personally questioned by two experienced doctors. The questionnaire included the primary information about demographics (eg, age, sex), smoking, drinking, height, and weight (Table 1). The clinical and pathological information of CRC cases was collected from their medical records.

Table 1 Distribution of selected demographic variables and risk factors in CRC cases and controls Notes: Italic values are statistically significant (P<0.05). aTwo-sided chi-squared test and Student's t-test. Abbreviations: CRC, colorectal cancer; SD, standard deviation; BMI, body mass index; AJCC, American Joint Committee on Cancer.

Selection of tagging SNPs

The tagging SNPs across the entire region of CTLA-4 gene (16.2 kbp spanning from 203862788–203878960 in chromosome 2 [upstream and downstream of gene extending 5 kb, respectively]) were selected from the Chinese Han in Beijing cohort via the HapMap Project (http://hapmap.ncbi.nlm.nih.gov/index.html.en) and analyzed with Haploview 4.2 software using a pairwise linkage disequilibrium (LD) r² threshold of 0.8 between SNPs (with a minimum LD of r²>0.8). SNPs with a Hardy–Weinberg equilibrium (HWE) P≥0.05, minor allele frequency (MAF) ≥0.05, and call rate ≥95% in the CHB cohort were included.²³ The detailed information of the selected four SNPs is summarized in Table 2.

Table 2 Primary information for CTLA-4 polymorphisms Abbreviations: SNP, single-nucleotide polymorphism; MAF, minor allele frequency; HWE, Hardy–Weinberg equilibrium.

DNA extraction and genotyping

Ethylenediamine tetraacetic acid (EDTA)-anticoagulated intravenous blood was collected after an overnight fast. The genomic DNA was isolated using the Promega DNA Blood Mini Kit (Promega, Madison, WI, USA).

The genotyping of the CTLA-4 rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A polymorphisms were performed by a custom-by-design 48-Plex SNPscan Kit (Genesky Biotechnologies Inc., Shanghai, China) as described in previous studies.^24,25 This 48-Plex SNPscan Kit was based on double ligation and multiplex fluorescence PCR.²⁶ For quality control, 92 (4%) samples were randomly selected and were tested again by the same genotyping method. The accordance ratio was 100%.

Statistical analysis

Statistical analysis was performed using SAS version 9.4 software package for Windows (SAS Institute, Cary, NC, USA), and a P<0.05 (two-tailed) was considered for level of significance. The continuous variables were expressed as mean±SD. We used the Student's t-test to check the differences for normally distributed continuous variables between CRC cases and controls. We used the chi-squared test to determine the differences in demographic variables, risk factors (smoking, BMI, and drinking), and the frequencies of genotype between CRC cases and controls. An internet-based calculator program (http://ihg.gsf.de/cgi-bin/hw/hwa1.pl) was harnessed to examine the deviation of HWE.¹⁴ Multivariate logistic regression was used to analyze the associations between CTLA-4 rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A SNPs and risk of CRC. The relationships between CTLA-4 rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A polymorphisms genotypes and risk of CRC were assessed by crude/adjusted ORs and the corresponding 95% CIs when appropriate. The relationships were assessed in additive, homozygote, dominant, and recessive models. We used a Bonferroni correction to adjust for multiple testing.^27,28

Results

Demographic characteristics

In our study, 1,003 CRC patients were included. Of them, 620 were males (61.81%) and 383 were females (38.19%). The mean age and SD were 61.10±12.17 years. The primary tumor site was the colon in 431 (42.97%) patients and the rectum in 572 (57.03%) patients. For the control group, we recruited 1,303 non-cancer controls, 801 males (61.47%) and 502 females (38.53%). Their age mean±SD was 61.40±9.61 years. All participants were Chinese Han population. The differences of age and sex between CRC and control groups were not statistically significant (P≥0.05) (Table 1). As summarized in Table 1, significant differences were found on alcohol consumption, smoking status, and BMI between the cases and the controls (P=0.002, <0.001, and <0.001, respectively). The primary information for CTLA-4 rs733618 T>C, rs3087243 G>A, rs16840252 C>T, and rs231775 G>A SNPs is shown in Table 2. For these SNPs, the successful ratio was more than 98.50%. In controls, MAF of CTLA-4 tagging SNPs was very close to the MAF data for Chinese (Table 2). Table 2 shows the genotype frequencies of CTLA-4 tagging SNPs polymorphisms were all in HWE.

Association of CTLA-4 rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A polymorphisms with CRC in overall analysis

The frequencies of CTLA-4 rs733618 TT, TC, and CC genotypes were 35.31%, 47.35%, and 17.35% in the CRC group and 35.23%, 47.15%, and 17.62% in the control group, respectively (Table 3). The frequencies of CTLA-4 rs16840252 CC, CT, and TT genotypes were 75.71%, 22.76%, and 1.53% in the CRC group and 77.38%, 21.77%, and 0.85% in controls, respectively (Table 3). There was no statistically significant difference in genotype distribution of CTLA-4 rs16840252 C>T and rs733618 T>C polymorphisms among CRC patients and controls. The frequencies of CTLA-4 rs3087243 GG, GA, and AA genotypes were 65.00%, 30.20%, and 4.80% in the CRC group and 65.38%, 31.38%, and 3.23% in controls, respectively (Table 3). The CTLA-4 rs3087243 AA genotype was associated with a borderline statistically increased risk of CRC, compared with CTLA-4 rs3087243 GG/GA genotypes (crude OR=1.51, 95% CI=0.99–2.31, P=0.058). When adjusted for age, sex, BMI, smoking, and drinking, a borderline statistically increased risk of CRC was also found (crude OR=1.52, 95% CI=0.99–2.34, P=0.058; Table 4).

Table 3 The frequencies of CTLA-4 rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A polymorphisms in CRC patients and controls Abbreviation: CRC, colorectal cancer.

Table 4 Overall and stratified analyses of CTLA-4 rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A polymorphisms with CRC by region Note: aAdjusted for age, sex, BMI, smoking status, and alcohol use in a logistic regression model. Abbreviations: CRC, colorectal cancer; BMI, body mass index; OR, odds ratio; CI, confidence interval.

The frequencies of CTLA-4 rs231775 GG, GA, and AA genotypes were 42.59%, 45.25%, and 12.16% in the CRC group and 47.77%, 43.31%, and 8.85% in the control group, respectively (Table 3). When compared with the CTLA-4 rs231775 GG genotype, CTLA-4 rs231775 AA and GA/AA genotypes significantly increased the risk of CRC (homozygote model: crude OR=1.47, 95% CI=1.10–1.95, P=0.008; and dominant model: crude OR=1.23, 95% CI=1.05–1.46, P=0.014). When compared with CTLA-4 rs231775 GG/GA genotypes, the CTLA-4 rs231775 AA genotype also increased the risk of CRC (crude OR=1.43, 95% CI=1.09–1.87, P=0.011). When adjusting for age, sex, BMI, smoking, and drinking, the results were not essentially changed (homozygote model: adjusted OR=1.40, 95% CI=1.05–1.87, P=0.022; dominant model: adjusted OR=1.19, 95% CI=1.00–1.41, P=0.047; and recessive model: adjusted OR=1.38, 95% CI= 1.05–1.82, P=0.021; Table 4).

Association of CTLA-4 rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A polymorphisms with CRC in a stratification group by primary site of tumor

To determine whether the effect of CTLA-4 tagging SNPs was modified by the primary site of tumor, we performed a stratified analysis. For CTLA-4 rs16840252 C>T, stratified analysis revealed this polymorphism was associated with an increased risk of colon cancer (homozygote model: adjusted OR=2.51, 95% CI=1.04–6.03, P=0.040 and recessive model: adjusted OR=2.54, 95% CI=1.06–6.09, P=0.037; Table 4). For the CTLA-4 rs231775 G>A polymorphism, we found that CTLA-4 rs231775 AA genotypes might be associated with an increased risk of colon cancer (homozygote model: adjusted OR=1.61, 95% CI=1.12–2.30, P=0.010 and recessive model: adjusted OR=1.59, 95% CI=1.13–2.23, P=0.009; Table 4). The results of other genetic comparisons are summarized in Table 4.

Association of CTLA-4 rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A polymorphisms with CRC in a stratification group by geographical area

In this study, CRC patients and controls were enrolled from two different geographical areas (Fuzhou and Zhenjiang, China). We performed a stratified analysis according to geographical area. Compared with CTLA-4 rs231775 GG, we found CTLA-4 rs231775 GA/AA genotypes might be associated with an increased risk of CRC in the Zhenjiang cohort (adjusted OR=1.38, 95% CI=1.01–1.88, P=0.041; Table 5). In addition, the stratified analysis revealed the CTLA-4 rs231775 G>A polymorphism also had a tendency of increased risk to CRC in the Fuzhou cohort (recessive model: adjusted OR=1.40, 95% CI=0.99–1.98, P=0.061; Table 5).

Table 5 Logistic regression analyses of associations between CTLA-4 polymorphisms and risk of CRC in two cohorts Note: aAdjusted for age, sex, BMI, smoking status, and alcohol use. Abbreviations: CRC, colorectal cancer; OR, odds ratio; CI, confidence interval; BMI, body mass index.

SNP haplotypes

Using an expectation–maximization algorithm (SHESIS program; Bio-X Inc., Shanghai, China, http://analysis.bio-x.cn/myAnalysis.php),²⁹ we constructed seven haplotypes (Table 6). Haplotype comparison analysis suggested that CTLA4 G_rs3087243C_rs16840252C_rs733618A_rs231775, G_rs3087243C_rs16840252T_rs733618A_rs231775, and other haplotypes significantly increased the risk of CRC (P<0.001, <0.001, and 0.002, respectively, Table 6).

Table 6 CTLA-4 haplotype frequencies (%) in cases and controls and risk of CRC Abbreviations: CRC, colorectal cancer; OR, odds ratio; CI, confidence interval.

Discussion

The individual's susceptibility to CRC may be diverse, even with the same environmental exposure. Host genetic predisposition may lead to these differences. In recent years, several case–control studies have been performed to test the hypothesis that some functional variants in CTLA-4 and other immune checkpoint molecules such as HLA-G may influence the risk and the treatment of CRC.^30–36 Garziera et al³⁶ reported that HLA-G 3′UTR polymorphisms might significantly affect the development of CRC. However, the association between CRC susceptibility and CTLA-4 polymorphisms remain conflicting. In addition, a comprehensive assessment was lacking. The aim of the present study was to identify the association between CTLA-4 tagging polymorphisms (rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A) and susceptibility of CRC in a case–control study. Genotyping of 1,003 CRC cases and 1,303 cancer-free controls was carried out in an Eastern Chinese Han population. Our findings demonstrated that CTLA-4 rs231775 G>A polymorphism might be associated with the risk of CRC. In a stratified analysis by primary site of tumor, the association was also found in colon cancer. When a subgroup analysis was performed by cohort, we also found CTLA-4 rs231775 GA/AA genotypes might be associated with an increased risk of CRC in the Zhenjiang cohort. Additionally, we found the CTLA-4 rs16840252 C>T polymorphism was associated with a risk of colon cancer. Haplotype comparison analysis showed that CTLA4 G_rs3087243C_rs16840252C_rs733618A_rs231775, G_rs3087243C_rs16840252T_rs733618A_rs231775, and other haplotypes increased the risk of CRC. Based on these primary findings, we found CTLA-4 tagging polymorphisms and haplotypes might influence the susceptibility to developing CRC.

Several case–control studies focused on the association between CTLA-4 rs16840252 C>T polymorphism and risk of cancer.^37,38 The observed results indicated that the CTLA-4 rs16840252 C>T polymorphism might not confer a risk to cancer. However, CTLA-4 rs16840252 C>T located in the promoter region of the CTLA4 gene. HapMap data suggest that CTLA-4 rs16840252 C>T and rs4553808 C>T (−1,661 C>T) are in complete LD.³⁷ Interestingly, Idris et al³⁹ also reported that strong LD was found between CTLA-4 rs16840252 C>T and rs5742909 C>T (−318 C>T) across all LD structures in an Asian population. In the presence of these functional SNPs on the same LD block, it could be that the predisposing allele of CTLA-4 rs4553808 C>T or rs5742909 C>T polymorphism is in LD with the protective allele of rs16840252 C>T. Ligers et al⁴⁰ found that individuals carrying thymine at position −318 of the CTLA4 promoter (CTLA-4 rs5742909 C>T) showed significantly increased expression, both of CTLA-4 mRNA in non-stimulated cells and of cell-surface CTLA-4 after cellular stimulation. Recently, several meta-analyses indicated that CTLA-4 rs4553808 T>C and rs5742909 C>T polymorphisms were associated with the risk of cancer, especially in Asians.^33,41,42 Since CTLA-4 rs16840252 C>T, rs5742909 C>T, and rs4553808 T>C are in strong LD, the function of CTLA-4 rs16840252 C>T could be influenced by CTLA-4 rs4553808 C>T or rs5742909 C>T. To the best of our knowledge, this case–control study was the first investigation to assess the association between CTLA-4 rs16840252 C>T genotype and CRC risk. Our findings indicated the CTLA-4 rs16840252 C>T polymorphism represented a risk factor for colon cancer. Our findings are supported by those pool-analyses mentioned above.

The CTLA-4 rs231775 G>A polymorphism was the most frequently explored and was established as a functional SNP of the CTLA-4 gene.^43,44 The CTLA-4 rs231775 G>A (c.49 G>A) SNP causes p.17Ala >17 Thr change in the leading sequence of CTLA-4 receptor.^40,44 Previous studies have demonstrated that the CTLA-4 rs231775 G allele has a lower mRNA efficiency and downregulates CTLA-4 protein more than the CTLA-4 rs231775 A allele.⁴⁵ Therefore, individuals who carry the CTLA-4 rs231775 AA genotype have lower T-cell proliferation and immune response than those with the CTLA-4 rs231775 GG genotype.⁴⁴ Sun et al⁴³ also found that the p.17Ala >17 Thr substitution in CTLA-4 amino acid residue caused by the c.49 G>A SNP significantly increased the interaction of the CTLA-4 receptor with its ligand B7.1, and recombinant CTLA-4-17Thr had a higher inhibitory effect to T-cell proliferation and immune response compared with CTLA-4-17Ala. These primary studies suggested that p.17Ala>17 Thr change in CTLA-4 may lead to a significant effect of T-cell proliferation and activation. A previous study demonstrated that donor CTLA-4 rs231775 genotype modulates the immune response to minor histocompatibility antigen mismatches.⁴⁶ The CTLA-4 rs231775 genotype was also considered as a genetic determinant in autoimmune Addison's disease.⁴⁷ Recently, a number of case–control studies focused on the relationship between CTLA-4 rs231775 G>A SNP and the risk of cancer, and results of subsequent meta-analyses evidenced that the CTLA-4 rs231775 G>A polymorphism was a risk factor for multiple cancer, especially in Asian populations.^32–35 Three pooled-analysis studies also suggested that this polymorphism was associated with the development of CRC.^15–17 Although these findings tried to suggest an association between CTLA-4 rs231775 G>A polymorphism and CRC, the number of included studies and participants were limited. Thus, we conducted this case–control study with larger sample sizes to explore whether the CTLA-4 rs231775 G>A polymorphism was a risk factor for CRC. As demonstrated in Table 4, we found that this polymorphism was associated with an increased risk of CRC. We also studied the association of CTLA-4 rs231775 G>A polymorphism with CRC in different subgroups. Similar findings were also found when the Bonferroni correction was applied. The association was also significant in the colon cancer subgroup (AA vs GG: OR=1.61; 95% CI=1.12–2.30; P=0.010 and AA vs GG/GA: OR=1.59; 95% CI=1.13–2.23; P=0.009; Table 4). Results of the present study were in accordance with results of those meta-analyses and functional studies mentioned above.

CTLA-4 rs3087243 G>A, rs16840252 C>T, rs733618 T>C, and rs231775 G>A variants may not be inherited randomly. As summarized in Table 6, we found the frequency of G_rs3087243C_rs16840252C_rs733618A_rs231775, G_rs3087243C_rs16840252T_rs733618A_rs231775, and other haplotypes was significantly increased in CRC patients. We first reported the association of these CTLA-4 haplotypes with CRC susceptibility. A previous study suggested that the CTLA4 G_rs3087243C_rs16840252C_rs733618A_rs231775 haplotype significantly increased the risk of gastric cardia adenocarcinoma,¹⁴ which was similar to our findings. However, these CTLA4 haplotypes only influenced a very minor fraction (less than 2%) of the CRC patients.

Of note, we focused on the relationship of CTLA-4 tagging SNPs with CRC risk in an Eastern Chinese Han population. In addition, the sample size of our study was larger than before. Finally, the MAF in our controls was very similar to the data for Chinese in the database (Table 2).

Although there were some merits in our study, some limitations should also be addressed. First, this study was designed as a hospital-based investigation; the CRC patients and controls were recruited from hospitals in Eastern China and might not well represent the whole Eastern Chinese Han population. Second, the recruited CRC cases were moderate in stratified analyses. In the future, these findings should be verified in well-designed studies with a larger sample size. Third, because of the limited sample size of CRC patients and absence of a validation cohort, the power of the present study may be insufficient, especially in stratified analyses. Fourth, for insufficient samples, a replicated study was not conducted. Fifth, due to lack of other information, we did not carry out a further evaluation of potential interaction, such as dietary habit, family history, hormone level, intake of vitamins, other environmental factors, and lifestyles. In considering the complexity of CRC etiology, the gene–environment interaction should not be ignored. Finally, in our case–control study, we investigated four tagging SNPs in the CTLA-4 gene and did not focus on other functional SNPs. In the future, a fine-mapping study is needed to further identify any potential association.

Conclusion

In summary, the findings of our case–control study evidence that CTLA-4 rs16840252 C>T and rs231775 G>A SNPs are correlated with genetic susceptibility for development of CRC in an Eastern Chinese Han population. Additionally, this study first highlights that CTLA-4 rs16840252 C>T polymorphism increases the susceptibility of CRC. Furthermore, findings are consistent with the biological functions of tagging SNPs in the CTLA-4 gene and validate the hypothesis that CTLA-4 tagging polymorphisms, which alter CTLA-4 mRNA and/or protein expression, may influence normal immune functions and lead to an increased risk of CRC.

Acknowledgments

We appreciate all subjects who participated in this study. We wish to thank Dr Yan Liu (Genesky Biotechnologies Inc., Shanghai, China) for technical support. This study was supported by the Natural Science Foundation of Universities and colleges of Jiangsu Province (Grant No 16KJB310002), Senior Talents Scientific Research Foundation of Jiangsu University (Grant No 16JDG066), Critical Patented Project of the Science and Technology Bureau of Fujian Province (Grant No 2013YZ0002-2), the Special Program for the Development of Strategic Emerging Industries of Fujian Province (Grant No 13YZ0201), the Natural Science Foundation of Fujian Province (Grant No 2015J01435, 2017J01259), the Fujian Provincial Health and Family Planning Research Talent Training Program (Grant No 2015-CX-7, 2018-ZQN-13, 2016-1-11, 2018-1-1), the Joint Funds for the Innovation of Science and Technology, Fujian Province (Grant No 2017Y9077), and the National Clinical Key Specialty Construction Program.

Disclosure

The authors report no conflicts of interest in this work.

References

1. Chen W, Zheng R, Baade PD, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–132.

2. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.

3. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61(2):69–90.

4. Katsidzira L, Gangaidzo IT, Makunike-Mutasa R, et al. A case–control study of risk factors for colorectal cancer in an African population. Eur J Cancer Prev. Epub 2018 Apr 11:1.

5. Markowitz SD, Bertagnolli MM. Molecular origins of cancer: molecular basis of colorectal cancer. N Engl J Med. 2009;361(25):2449–2460.

6. Lichtenstein P, Holm NV, Verkasalo PK, et al. Environmental and heritable factors in the causation of cancer – analyses of cohorts of twins from Sweden, Denmark, and Finland. N Engl J Med. 2000;343(2):78–85.

7. Hemminki K, Chen B. Familial risk for colorectal cancers are mainly due to heritable causes. Cancer Epidemiol Biomarkers Prev. 2004;13(7):1253–1256.

8. Hodi FS, Mihm MC, Soiffer RJ, et al. Biologic activity of cytotoxic T lymphocyte-associated antigen 4 antibody blockade in previously vaccinated metastatic melanoma and ovarian carcinoma patients. Proc Natl Acad Sci U S A. 2003;100(8):4712–4717.

9. Walunas TL, Lenschow DJ, Bakker CY, et al. CTLA-4 can function as a negative regulator of T cell activation. Immunity. 1994;1(5):405–413.

10. van der Merwe PA, Bodian DL, Daenke S, Linsley P, Davis SJ. CD80 (B7-1) binds both CD28 and CTLA-4 with a low affinity and very fast kinetics. J Exp Med. 1997;185(3):393–404.

11. Keilholz U. CTLA-4: negative regulator of the immune response and a target for cancer therapy. J Immunother. 2008;31(5):431–439.

12. Engelhardt JJ, Sullivan TJ, Allison JP. CTLA-4 overexpression inhibits T cell responses through a CD28-B7-dependent mechanism. J Immunol. 2006;177(2):1052–1061.

13. Teft WA, Kirchhof MG, Madrenas J. A molecular perspective of CTLA-4 function. Annu Rev Immunol. 2006;24:65–97.

14. Tang W, Wang Y, Chen S, et al. Investigation of cytotoxic T-lymphocyte antigen 4 polymorphisms in gastric cardia adenocarcinoma. Scand J Immunol. 2016;83(3):212–218.

15. Wang L, Jing F, Su D, et al. Association between CTLA-4 rs231775 polymorphism and risk of colorectal cancer: a meta analysis. Int J Clin Exp Med. 2015;8(1):650–657.

16. Wang Y, Wang X, Zhao R. The association of CTLA-4 A49G polymorphism with colorectal cancer risk in a Chinese Han population. Int J Immunogenet. 2015;42(2):93–99.

17. He L, Deng T, Luo HS. Association between cytotoxic T-lymphocyte antigen-4 +49A/G polymorphism and colorectal cancer risk: a meta-analysis. Int J Clin Exp Med. 2015;8(3):3752–3760.

18. Martin M, Schneider H, Azouz A, Rudd CE. Cytotoxic T lymphocyte antigen 4 and CD28 modulate cell surface raft expression in their regulation of T cell function. J Exp Med. 2001;194(11):1675–1682.

19. Chung KY, Gore I, Fong L, et al. Phase II study of the anti-cytotoxic T-lymphocyte-associated antigen 4 monoclonal antibody, tremelimumab, in patients with refractory metastatic colorectal cancer. J Clin Oncol. 2010;28(21):3485–3490.

20. Tang W, Zhang S, Qiu H, et al. Genetic variations in MTHFR and esophageal squamous cell carcinoma susceptibility in Chinese Han population. Med Oncol. 2014;31(5):915.

21. Zhai Y, Zhao WH, Chen CM. [Verification on the cut-offs of waist circumference for defining central obesity in Chinese elderly and tall adults]. Zhonghua Liu Xing Bing Xue Za Zhi. 2010;31(6):621–625.

22. Zhang X, Zhang S, Li Y, et al. Association of obesity and atrial fibrillation among middle-aged and elderly Chinese. Int J Obes. 2009;33(11):1318–1325.

23. Carlson CS, Eberle MA, Kruglyak L, Nickerson DA. Mapping complex disease loci in whole-genome association studies. Nature. 2004;429(6990):446–452.

24. Zheng L, Yin J, Wang L, et al. Interleukin 1B rs16944 G>A polymorphism was associated with a decreased risk of esophageal cancer in a Chinese population. Clin Biochem. 2013;46(15):1469–1473.

25. Yin J, Wang L, Shi Y, et al. Interleukin 17A rs4711998 A>G polymorphism was associated with a decreased risk of esophageal cancer in a Chinese population. Dis Esophagus. 2014;27(1):87–92.

26. Yin J, Wang X, Wei J, et al. Interleukin 12 B rs3212227 T>G polymorphism was associated with an increased risk of gastric cardiac adenocarcinoma in a Chinese population. Dis Esophagus. 2015;28(3):291–298.

27. Bland JM, Altman DG. Multiple significance tests: the Bonferroni method. BMJ. 1995;310(6973):170.

28. Lesack K, Naugler C. An open-source software program for performing Bonferroni and related corrections for multiple comparisons. J Pathol Inform. 2011;2:52.

29. Shi YY, He L, Shesis HL. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005;15(2):97–98.

30. Garziera M, Virdone S, de Mattia E, et al. HLA-G 3′UTR polymorphisms predict drug-induced G3-4 toxicity related to folinic acid/ 5-fluorouracil/oxaliplatin (FOLFOX4) chemotherapy in non-metastatic colorectal cancer. Int J Mol Sci. 2017;18(7):1366.

31. Garziera M, Bidoli E, Cecchin E, et al. HLA-G 3′UTR polymorphisms impact the prognosis of stage II–III CRC patients in fluoropyrimidine-based treatment. PLoS One. 2015;10(12):e0144000.

32. Zheng J, Yu X, Jiang L, et al. Association between the cytotoxic T-lymphocyte antigen 4 +49G > A polymorphism and cancer risk: a meta-analysis. BMC Cancer. 2010;10:522.

33. Geng R, Song F, Yang X, et al. Association between cytotoxic T lymphocyte antigen-4 +49A/G, −1722T/C, and −1661A/G polymorphisms and cancer risk: a meta-analysis. Tumour Biol. 2014;35(4):3627–3639.

34. Zhang Y, Zhang J, Deng Y, et al. Polymorphisms in the cytotoxic T-lymphocyte antigen 4 gene and cancer risk: a meta-analysis. Cancer. 2011;117(18):4312–4324.

35. Wang L, Jiang Z, Qiu H, Tang W, Duan T, Wang L. Associations between CTLA-4 +49 A/G (rs231775) polymorphism and cancer risk: a meta-analysis based on 52 case-control studies. Int J Clin Exp Med. 2015;8(5):6835–6851.

36. Garziera M, Catamo E, Crovella S, et al. Association of the HLA-G 3′UTR polymorphisms with colorectal cancer in Italy: a first insight. Int J Immunogenet. 2016;43(1):32–39.

37. Welsh MM, Applebaum KM, Spencer SK, Perry AE, Karagas MR, Nelson HH. CTLA4 variants, UV-induced tolerance, and risk of non-melanoma skin cancer. Cancer Res. 2009;69(15):6158–6163.

38. Bouwhuis MG, Gast A, Figl A, et al. Polymorphisms in the CD28/CTLA4/ICOS genes: role in malignant melanoma susceptibility and prognosis? Cancer Immunol Immunother. 2010;59(2):303–312.

39. Idris ZM, Miswan N, Muhi J, Mohd TA, Kun JF, Noordin R. Association of CTLA4 gene polymorphisms with lymphatic filariasis in an East Malaysian population. Hum Immunol. 2011;72(7):607–612.

40. Ligers A, Teleshova N, Masterman T, Huang WX, Hillert J. CTLA-4 gene expression is influenced by promoter and exon 1 polymorphisms. Genes Immun. 2001;2(3):145–152.

41. Yan Q, Chen P, Lu A, Zhao P, Gu A. Association between CTLA-4 60G/A and −1661A/G polymorphisms and the risk of cancers: a meta-analysis. PLoS One. 2013;8(12):e83710.

42. Tang W. Relationship between cytotoxic T-lymphocyte antigen 4 (CTLA-4) rs5742909 C>T polymorphism and cancer risk: a meta-analysis based on thirty case-control studies. Int J Clin Exp Med. 2016;9(8):15191–15203.

43. Sun T, Zhou Y, Yang M, et al. Functional genetic variations in cytotoxic T-lymphocyte antigen 4 and susceptibility to multiple types of cancer. Cancer Res. 2008;68(17):7025–7034.

44. Sun T, Hu Z, Shen H, Lin D. Genetic polymorphisms in cytotoxic T-lymphocyte antigen 4 and cancer: the dialectical nature of subtle human immune dysregulation. Cancer Res. 2009;69(15):6011–6014.

45. Chistiakov DA, Savost’anov KV, Turakulov RI, Efremov IA, Demurov LM. Genetic analysis and functional evaluation of the C/T(−318) and A/G(−1661) polymorphisms of the CTLA-4 gene in patients affected with Graves’ disease. Clin Immunol. 2006;118(2–3):233–242.

46. Gallardo D, Bosch-Vizcaya A, Rodríguez-Romanos R, et al. Donor CTLA-4 genotype modulates the immune response to minor histocompatibility antigen mismatches. Biol Blood Marrow Transplant. 2017;23(12):2042–2047.

47. Wolff AS, Mitchell AL, Cordell HJ, et al. CTLA-4 as a genetic determinant in autoimmune Addison’s disease. Genes Immun. 2015;16(6):430–436.