Risk allelic load in Th2 and Th3 cytokines genes as biomarker of susceptibility to HPV-16 positive cervical cancer: a case control study
BMC Cancer volume 16, Article number: 330 (2016)
Alterations in the host cellular immune response allow persistent infections with High-Risk Human Papillomavirus (HR-HPV) and development of premalignant cervical lesions and cervical cancer (CC). Variations of immunosuppressive cytokine levels in cervix are associated with the natural history of CC. To assess the potential role of genetic host immunity and cytokines serum levels in the risk of developing CC, we conducted a case–control study paired by age.
Peripheral blood samples from patients with CC (n = 200) and hospital controls (n = 200), were used to evaluate nine biallelic SNPs of six cytokine genes of the adaptive immune system by allelic discrimination and cytokines serum levels by ELISA.
After analyzing the SNP association by multivariate logistic regression adjusted by age, CC history and smoking history, three Th2 cytokines (IL-4, IL-6 and IL-10) and one Th3 (TGFB1) cytokine were significantly associated with CC. Individuals with at least one copy of the following risk alleles: T of SNP (−590C > T IL-4), C of SNP (−573G > C IL-6), A of SNP (−592C > A IL-10), T of SNP (−819C > T IL-10) and T of SNP (−509C > T TGFB1), had an adjusted odds ratio (OR) of 2.08 (95 % CI 1.475–2.934, p = 0.0001), an OR of 1.70 (95 % CI 1.208–2.404, p = 0.002), an OR of 1.87 (95 % CI 1.332–2.630, p = 0.0001), an OR of 1.67 (95 % CI 1.192–2.353, p = 0.003) and an OR of 1.91 (95 % CI 1.354–2.701, p = 0.0001), respectively, for CC. The burden of carrying two or more of these risk alleles was found to have an additive effect on the risk of CC (p trend = 0.0001). Finally, the serum levels of Th2 and Th3 cytokines were higher in CC cases than the controls; whereas IFNG levels, a Th1 cytokine, were higher in controls than CC cases.
The significant associations of five SNPs with CC indicate that these polymorphisms are potential candidates for predicting the risk of development of CC, representing a risk allelic load for CC and can be used as a biomarker of susceptibility to this disease.
Cervical cancer (CC) ranks third as a cause of death among women worldwide, with an estimated overall mortality rate of 15 per 100,000 women . CC was the second most common cause of death in Mexican women in 2011 (10.4 %) . The immune system plays key roles during HPV-associated carcinogenesis, as HPV clearance is determined by specific immunological reactions . Thus, CC seems to be due in part to a failure of the immune system, which is unable to eliminate persistent HPV infections and virus-transformed cells .
A local predominance of T helper 2 (Th2) cytokine profile expression (IL-4, IL-6, IL-10) and an impaired cellular immune response induced by immune suppressor cytokines, such as IL-10 and Transforming Growth Factor Beta 1 (TGFB1) in association with a diminished Th1 profile, has been demonstrated in patients with CC [5, 6]. Cytokines are key players in the immune system and several polymorphisms identified in their gene promoters may be responsible for variations in expression levels observed between individuals in different diseases .
Single nucleotide polymorphisms (SNPs) in IL-4, IL-6, IL-10, TGFB1, Tumor Necrosis Factor-alpha (TNF) and Interferon-gamma (IFNG) genes that cause variations in host immune response may contribute to CC risk. Since the boom of the genomic era, there are have been numerous reports of polymorphisms associated with CC risk; however, very few of the studied polymorphisms in this neoplasm have shown consistent associations across studies and populations. A large number of SNPs have been identified in IL-4, IL-6, IL-10, TGFB1, TNF e IFNG and many of them have been studied regarding their association with CC in different populations, particularly -592C > A, −819C > T and -1082A > G of IL-10, −509C > T TGFB1 and -308G > A of TNF .
In this study, we analyzed the association between -590C > T (IL-4), −573G > C (IL-6), −592C > A (IL-10), −819C > T (IL-10), −1082A > G (IL-10), −509C > T (TGFB1), −800G > A (TGFB1), −308G > A (TNF) and -1615C > T (IFNG) gene promoter polymorphisms and the risk for CC in Mexican women. These SNPs were selected based on two criteria: validated SNPs for frequency or for utilization in the HAPmap Project and SNPs in the promoter region with potential role in transcriptional regulation of each cytokine evaluated. Furthermore, we determined whether there are variations in serum protein levels of these cytokines according to diagnostic and allelic variation of each SNP.
Study design and population
A clinic-based case–control study matched by age (±5 years) 1:1, (n = 200 per group) was conducted at the National Cancer Institute’s (INCan) Gynecology Service in Mexico City between September 2010 and December 2011. Cases with cervical squamous cell carcinoma confirmed by two pathologists were included, and women with a negative Papanicolaou study and a normal colposcopy were selected as controls and matched to the cases by age (±5 years). All participants were HPV 16-positive, since this was an inclusion criterion. The analyzed population was Mexican mestizo with a history of three previous generations being born in Mexico and with a residency period of > 1 year in the study area. To eliminate any selection bias, the authors screened CC cases and controls to ensure that they had never been diagnosed with chronic inflammatory diseases. The Bioethics and Research Committees at INCan (reference INCan/CC/326/10CB/609) and at the Instituto Nacional de Salud Pública (INSP, reference CI814 and CI1289 No. 1624) approved the study, which was carried out according to the Helsinki Declaration. All participants signed an informed consent to participate in this study, agreeing also to the subsequent publication of the results. Each subject was interviewed regarding lifestyle and socio-demographic and hormonal factors known to be associated with an increased risk of CC.
Specimen collection and sample processing
Cervical epithelial cell scrapings from the controls and fresh cell biopsies from the women diagnosed with CC were used for this study. Peripheral blood mononuclear cells (PBMC) from all subjects were obtained by Ficoll-hypaque density gradients (Hystopaque, Sigma Chemical Co). Genomic DNA was extracted from PBMC using the Genomic DNA Purification Kit (Fermentas Life Sciences, Lithuania) and from cervical epithelial scrapings and biopsies previously digested with proteinase K. DNA concentration and purity were evaluated with a Thermo Scientific NanoDropTM 1000 Spectrophotometer (260/280) and the integrity of the DNA was determined by electrophoresis in agarose gels at 0.8 %.
Cervical epithelial scrapings and biopsy specimens were tested for HPV by PCR amplification, using consensus primers MY09/MY11, LIC1/LIC2, and GP5/GP6. GAPDH (250 bp) was used as an internal control for DNA quality; SiHa was used as a positive control, and deionized H2O as a negative control. The purified DNA band was sequenced using the Sanger method. After analyzing the sequences by BLAST, only HPV 16-positive women were included in this study. Serum protein levels of IL-4, IL-6, IL-10, TGFB1, TNF, and IFNG were determined using the Human high sensitivity ELISA kit (Abcam, Cambridge UK). All assays were done by duplicate and the final concentration (picograms per millilitre [pg/ml]) and TGFB1 (nanograms per millilitre [ng/ml]) of each cytokine corresponded to the average of the duplicate readings.
The SNPs selection criteria were as follows: 1) Validated SNPs for frequency or for utilization in the HAPmap Project; 2) SNPs in the promoter region with potential role in the transcriptional regulation of each cytokine evaluated (the Ensemble program was used for this selection); 3) SNPs in IL-4, IL-6, IL-10, TGFB1, TNF and IFNG promoter, located in the binding sites of the transcription factors that potentially influence the transcriptional activity, reported in the following database: SNPper.URL: http//snpper.chip.org. Nine SNPs in the promoter region were genotyped using PCR with TaqMan fluorogenic probes in the ViiA™ 7 (Applied Biosystems, Foster City, CA, USA), (−590C > T (IL-4) rs2243250, −573G > C (IL-6) rs1800796, −592C > A (IL-10) rs1800872, −819C > T (IL-10) rs1800871, −1082A > G (IL-10) rs1800896, −509C > T (TGFB1) rs1800469, −800G > A (TGFB1) rs1800468, −308G > A (TNF) rs1800629 and -1615C > T (IFNG) rs2069705). All tests were performed in duplicate. The alleles were assigned using the ViiA™ 7 software (Applied Biosystems). A call rate of 0.99 for controls and CC was used for quality control. When the call rate was less than 0.99, DNA was reextracted and a new genotyping was performed (only six samples required this procedure).
Reproductive factors and sexual lifestyle were evaluated by multinomial logistic regression models adjusted for age. The Hardy-Weinberg equilibrium for each SNP was assessed using the allelic frequencies of the control group. Genotype-specific and allele risks were estimated as odds ratios (OR) with associated 95 % confidence intervals (CI). All p-values were based on a two-sided hypothesis test using logistic regression analyses in the three inheritance models and adjusting by potential confounders (age, CC history and smoking history). The Bonferroni method was used to correct the multiple comparisons (α = 0.05/9 = 0.0055). The effect of having one or more risk-associated alleles with CC was evaluated using multiple logistic regressions. All possible 2-way interactions among SNPs and between SNPs and serum level of IL-4, IL-6, IL-10, TGFB1, TNF e IFNG were tested by multivariate logistic models. Tertiles for the serum levels of each cytokine according to the respective observed distribution in the HPV-positive control group were created to evaluate their association with CC by multinomial logistic regression models. The mean differences of the serum levels of cytokines between CC and controls stratified by genotypes were assessed by linear regression models adjusted by age, CC history and smoking history.
Stratification by ethnicity that could potentially confound genetic analyses was avoided by confining our analysis to the Mexican mestizo population with two previous generations born in Mexico. The power calculation for each of the analyzed SNPs was evaluated taking into account the value of n fixed for the study and the minor allele frequency for each SNP in the NCL (P1), with an expected OR of 1.2, 1.5, 2, 2.5 and 3. The established alpha value was 0.05. The P2 calculation was performed using the following formula: P1*OR/(1-P1) + (P1*OR). The statistic power obtained in the regression analysis for each SNP was as follows: −590 C > T of IL-4 (0.99); −573 G > C of IL-6 (0.92); −592 C > A (0.98), −819 C > T (0.91) and −1082 A > G (0.14) of IL-10; −509 C > T of TGFB1 (0.98); −308 G > A of TNF (0.28) and −1615 C > T of IFNG (0.99). All the statistical analyses were performed in the STATA program, version 13.0 (StataCorp, Collage Station, TX, EUA).
The analyses confirmed that known reproductive and sexual lifestyle risk factors for CC were statistically different between study groups. As expected, we found a significant positive association of CC diagnosis with age at first intercourse, parity, number of lifetime sexual partners, cancer family history and smoking history. The study did not find the history of previous Sexually Transmitted Diseases (STDs) and use of contraceptive methods to be risk factors for CC in this population (Table 1). Given the best score in the goodness of fit tests performed for all the logistic models evaluated, we carried out further multinomial logistic regression analysis, adjusting for cancer family history, smoking history and age. We observed no significant deviations from HWE among controls in any of the genotyped SNPs. The SNP -800G > A (TGFB1) rs1800468, was found not to be polymorphic with a minor allele frequency (MAF) of <1 % (Table 2).
We found a highly significant positive association with CC for minor allele homozygotes of -590C > T (IL-4), −573G > C (IL-6), −592C > A (IL-10), −819C > T (IL-10) and -509C > T (TGFB1). Odds ratios and p trend were OR = 4.36, 95 % CI 2.166–8.803 (p = 0.0001); OR = 3.2, 95 % CI 1.572–6.535 (p = 0.002); OR = 3.28, 95 % CI 1.660–6.508 (p = 0.0001); OR = 2.6 95 % CI 1.330–5.116 (p = 0.004); and OR = 3.73, 95 % CI 1.842–7.572 (p = 0.0001), respectively. For -1615C > T (IFNG) minor allele homozygotes, we also found a significant negative association with CC (OR = 0.11, 95 % CI (0.038–0.351)), accompanied with a significant negative trend (p = 0.0001).
In order to explore whether the evaluated SNPs could act together to increase CC risk, allele load was calculated (Table 3). We indeed found that having two or more risk alleles had a statistically significant association with CC, compared with having no risk alleles, p trend = 0.0001. Levels in serum of IL-4, IL-6 and IL-10 were much higher in patients with CC than in HPV-positive control patients (p = 0.00001), (Fig. 1). In contrast, IFNG and TNF serum levels were lower than those observed in HPV-positive control patients (p = 0.00001), (Fig. 2). Serum levels of TGFB1 were significantly higher in comparison with the other cytokines and the difference between control and CC patients was statistically significant (p < 0.00001) (Fig. 3).
We stratified the SNPs of interest by polymorphism genotypes to explore whether they had a relationship with levels of serum cytokines seen in CC patients. Only polymorphisms -592C > A and -819C > T of IL-10, −308G > A (TNF) and -509C > T (TGFB1) showed statistically significant differences in serum cytokines levels between minor allele homozygous and ancestral allele homozygous in CC patients (Table 4). When we evaluated the interaction between each of the SNPs: −590C > T (IL-4), −573G > C (IL-6), −592C > A, −819C > T and -1082A > G (IL-10), −509C > T (TGFB1), −308G > A (TNF) and -1615C > T (IFNG) and respective serum levels, as well as the SNP-SNP interaction, no statistically significant interaction on CC was found (data not shown).
The main findings of this study were a significant positive association between CC with the SNP’s: −590C > T (IL-4), −573G > C (IL-6), −592C > A (IL10), −819C > T (IL10) and -509C > T (TGFB1) and their serum levels. Although the effect of cytokine gene polymorphisms on CC have been reported in different populations, our results show for the first time a risk allelic load in Th2 and Th3 cytokines genes as a biomarker of susceptibility to CC.
HR-HPV persistent infection is a necessary but not sufficient cause for CC ; a combination of several risk factors is required for the disease to develop. The association found in this study for the variables early-age of first sexual intercourse, high number of sexual partners and multiparity, confirms the known association reported in previous epidemiological studies .
Genetic predisposing factors may influence the likelihood of, sensitivity to or persistence of HPV infection, as well as the rate of tumor development . Various immune-suppressive states associated with a reduction in cellular immune responses are associated with increasingly severe cervical dysplasia and increased viral shedding, suggesting that a T-helper bias that opposes Th1 responses (ie, Th2 bias) might predispose to chronic infection [12–15].
To our knowledge, this is the first study to show that the IL-4 SNP rs2243250 is significantly associated with the risk of CC. IL-4 is an anti-inflammatory cytokine and the T allele of the SNP -590C > T has been associated with increased transcriptional activity in vitro . The exact mechanism for this increment is not known. However, since the SNP is located within 5′UTR of the gene, it may be possible that alterations in this gen could be influencing in its transcription and/or mRNA stabilization . The IL-4 -590 SNP is a transition (C → T) that has been associated with oral cancer and is a suitable genetic marker for screening for this condition . Increased plasma concentrations of IL-4 can influence the immune status of an affected individual through several mechanisms and result in many phenotypes beyond the scope of the immune system .
Likewise, our results showed that the genotype C/C of the IL-6 SNP -573G > C (rs1800796, previously denoted as -572G > C), was significantly associated with the risk of CC when compared with the heterozygous genotype GC. The IL-6 -573 SNP is a transversion (G → C) that was reported as a genetic risk factor of lung cancer risk in Singaporean Chinese non-smoking females . IL-6 is a multifunctional cytokine that can regulate immune and inflammatory responses. In CC patients, high expression of IL-6 correlate with a promoting effect in tumor cell growth by autocrine and/or paracrine processes . Accumulated data of over-expression of IL-6 mRNA and protein in CC cells have showed that the IL-6 protein played important roles in CC development [12, 21, 22]. High levels of IL-6 are associated with the severity of disease and may promote tumor angiogenesis and cancer, which could be genetically determined at an individual level . Clinical studies in CC patients report that elevated IL-6 serum levels are associated with poor prognosis . Variations in the IL-6 gene have been found to be associated with the plasma levels of the protein  and functional relevance has been ascribed to several IL-6 variants located in the promoter region, including -573G > C . However, in this study the concentrations of angiogenic cytokines, such as IL-6, do not vary with genotype. This result is similar to other study which explored whether polymorphisms in angiogenic cytokine genes may affect the levels of cytokines and which reported no variation .
With respect to IL-10 cytokine gene polymorphisms it has been reported that several important polymorphic sites in the IL-10 gene, including three in the promoter region (−1082 A > G, −819 C > T, −592 C > A) may influence the transcription of IL-10 messenger RNA and the expression of IL-10 in vitro . Several studies have investigated the possible role of IL-10 cytokine gene polymorphisms in CC . The current study has shown that individuals homozygous for the A-allele of the −592 SNP are at three time’s greater odds of having CC as compared to controls. The −592 SNP rs1800872 is a transversion (C → A) that is located in the IL-10 promoter in a region with negative enhancer activity and is associated with loss of this activity , between putative consensus binding sequences for Sp1 and a sequence with similarity to that recognized by members of the ETS family proteins. Steinke et al., have demonstrated that the C to A nucleotide exchange results in increased IL-10 gene promoter activity and that the C variant of this SNP is a repressor element . A previous study carried out in 695 CC patients, 115 family-based patients and 586 unrelated controls, in Caucasian population, revealed an increased risk CC for individuals heterozygous for the A-allele of this SNP . Like our study, two recent meta-analyses reported this SNP as a risk factor for developing CC, especially for Asians [28, 31].
Only one other study has investigated the IL-10 SNP -819C > T rs1800871 and CC. In this study the CT + TT genotype combination and T allele, was slightly higher in 150 CC cases as compared with 162 age- and ethnically-matched cervical cytology negative healthy controls. However, contrary to our investigation, the study did not show statistical significance . At this time, a series of seven molecular epidemiological studies and meta-analyses have investigated the association between the IL-10 -819 SNP, a transition (C → T) and the susceptibility to different cancer types among different populations . Nevertheless, the results from these studies are inconsistent. Most studies investigating this SNP have focused on gastric cancer and some have found a reduced risk of gastric cancer among Asians but not among Caucasians . The mechanism of the influence of IL-10-819 SNP on carcinogenesis is still unknown .
The IL-10 -1082 SNP rs1800896 is a transition (A → G) and is the most extensively studied polymorphism in the IL-10 gene in cancer susceptibility . Studies link this SNP to high/low IL-10 producer status. Functional analyses have shown that the −1082 region contains a putative ETS-like transcription factor-binding site and that nuclear factors from a monocyte cell line bind to this region. Transient transfection studies in an Epstein-Barr virus-transformed B cell line indicated that the −1082 A allele confers a two fold increase in transcriptional activity of the IL-10 promoter compared to the G allele . This polymorphism has been associated with CC risk in populations from Zimbabwe and Japan [30, 37], while studies in The Netherlands, South Africa, Hungary and China, similarly to our study did not find any association . In the Japanese population, IL-10 -1082 genotypes corresponding to high production of interleukin (GA and GG) were significantly associated with CC severity . The discrepancy of reported results among studies may be explained by the frequency of the G allele in healthy controls [30, 35, 37]. A previously meta-analyses reported in the subgroup analysis by ethnicity of the association of IL-10 -1082A allele with decreased CC susceptibility among whites only (A vs G: OR, 0.39; 95 % CI, 0.32–0.47) . As expected, the G allele frequency obtained in this study (0.25) was similar to that reported for Italy (0.37)  and Argentina (0.32) , showing no relation with CC risk.
On the other hand, the common TGFB1 promoter SNP -509C-T (rs1800469) was associated with CC in our study, for all models of inheritance and analysis alleles. This SNP is a transition (C → T) linked to a nearly twofold difference in plasma levels among individuals and with risk, progression, and outcome of numerous diseases . This association with the pathogenesis of certain TGFB1-related diseases can be due to transcriptional suppression by AP1 binding to -509C .
Only one other study has investigated the SNP -509C > T rs1800469 and CC. In this study 150 CC patients with -509TT had marginal low risk for stage I (p = 0.04, OR = 0.95, 95 % CI = 0.91–0.99). However, −509TT genotype of TGFB1 was associated with increased risk of stage II of cancer (p = 0.07, OR = 3.13, 95 % CI = 0.87–11.14). In gene-environment interaction, carriers of TGFB1 -509TT genotype with tobacco usage were at higher risk of CC (OR = 3.67, 95 % CI = 0.38–35.1) . The association found in this study is according to previous report where TGF-beta family members play a key role in cellular growth, proliferation, differentiation and angiogenesis of several cancer types  and a marginal protection for early stage 1B but risk for stage II of CC of TGFB1 -509 T allele has been reported .
Our results indicated that TNF −308 SNP rs1800629, a transition (G → A), had no effect on CC risk in the overall analysis. This polymorphism is the most studied polymorphism in terms of association with CC and severity with conflicting results. Four recent meta-analysis, pooled data from eight, twelve, fifteen and eighteen case–control studies respectively, have suggested that TNF-308G > A polymorphism is associated with an increased risk of CC [43–46]. The analysis stratified by ethnicity, showed that there was a significant association of this polymorphism and increased risk of CC in Asians  and in Caucasian and African populations . In our study, the lack of association between TNF −308 G > A polymorphism and CC was consistent with previous reports from African [48, 49] and Caucasian ethnicity [39, 50, 51]. Differences in ethnic compositions; genetic background, sizes of samples, population stratification, or selection of symptoms of disease between analyzed populations could explain these controversial results .
Additionally, our results indicate that the variant heterozygote CT and the homozygous TT of SNP -1615C > T of IFNG was associated with a decreased risk of CC. To our knowledge, this is the first study to show that the IFNG SNP rs2069705 was significantly associated with decreased risk of CC. Studies so far have only investigated this SNP in connection with the risk of breast cancer and SLE [53, 54]. Based on the evidence that the T allele is associated with a higher level of IFNG than the C allele the association of this SNP with this type of cancer is predictable, given that at breast cancer diagnosis, inflammatory cytokines are detected at elevated levels . In contrast, a protective effect for the T allele in CC could be related to an increased expression of IFNG and play a role in protection against persistent HPV infection associated to CC development .
Nevertheless, a consensus seems to be emerging in that the combined effect of several cytokine SNPs may play a more crucial role in development disease. Thus, in order to explore whether the polymorphisms we examined could act together to increase risk for CC, we combined alleles and tested the effect on risk for the combinations. Individuals carrying two or more of the risk-associated alleles (−590 (IL-4)/-573 (IL-6)/-592 (IL-10)/-819 (IL-10)/-509 (TGFB1) were at increased risk for CC, compared to those with no risk-associated alleles. Although our study did not have the power to examine multiplicative interactions of specific genotypes, the combination of risk-associated alleles produced notably high odds ratios. This exploratory approach supported the idea that combinations of risk-associated alleles have an additive effect on the risk of CC.
In addition, we have found increased serum level of IL-4, IL-6, IL-10 and TGFB1 cytokines among our CC patients. Levels of IFNG were lower in the group of women with CC compared to controls. This is in accordance with previous works indicating a shift towards a T helper 2-type (Th2) cytokine pattern in CC  and that increased serum levels of IL-6, IL-10, and TGFB1 are associated with progression of CC during different stages, all of which have the potential to be angiogenic amplifiers .
Given that cytokines contribute importantly to each step of cancer development and progression, and that deregulated levels of cytokines and cytokine receptors can be detected in cancer patients locally and systemically , this study demonstrated an anti-inflammatory genetic profile, which is associated with increased susceptibility to cervical cancer. However, no statistically significant interactions were obtained between the polymorphisms cytokines of interest and respective serum levels.
IL-4 could possibly be exerting inhibitory effects on the expression and release of proinflammatory cytokines, IL-6 could be inducing prometastatic genes that subsequently lead to proliferation, and prolonged survival of cancer cells, given that the IL-6 serum levels are associated with poor prognosis . Increased serum IL-10 levels could facilitate development of tumors by suppressing the expression of MHC class I and II antigens and preventing tumor antigen presentation to CD8-cytotic T lymphocytes . In addition to its immunosuppressive effects, IL10 can induce transcription of the HPV16 E7 oncoprotein . Additionally, in one of our previous studies, systemic IL-10 mRNA expression level and the IL-10 protein level in serum are significantly higher in cervical lesions compared to women without lesions. We postulate that IL-10 and TGFB1 induce immune system evasion through an immunosuppressive state in the environment of the cervix in women infected with HPV .
Additionally, we are aware of a few limitations in our study, some of which cannot be overcame. First, −based in the patients and the control group randomly from the same hospital, we cannot completely rule out the inherent selection bias. Nevertheless, in order to minimize potential biases, we have made a rigorous epidemiological design of study subjects and more statistical adjustments for known risk factors. Second, our sample is medium sized, which may weaken the statistical power of this study. Under the current sample size, we have 99 %, 98 %, 92 % and 91 % power at a 0.05 significance level to detect an OR of 1.2 or higher for −590 IL-4 and −1615 IFNG, for −592 IL-10 and −509 TGFB1, for −573 IL-6 and for −819 IL-10 polymorphisms, respectively. However, well-designed large, prospective studies with detailed information about HPV infection are required to validate our findings and confirm this hypothesis.
In conclusion, our results suggest an association of three SNP’s Th2 cytokines (IL-4, IL-6 and IL-10) and Th3 (TGFB1) cytokine with CC. Evaluating the association of the burden of risk alleles −590 (IL-4)/-573 (IL-6)/-592 (IL-10)/-819 (IL-10)/-509 (TGFB1) with CC, we found that carrying three or more risk alleles has an additive effect on the risk of CC. We demonstrate a significant difference in the level of serum proteins IL-4, IL-6, IL-10 and TGFB1 in cases compared to controls (p < 0.001). Thus, women with CC have an imbalance in the pattern of cytokines Th1/Th2 given systemically by an increase in type II cytokines (IL-4, IL-6 and IL-10, suppressor of the cellular immune response) and reduced in type I cytokine (INFG), allowing HPV persistence and development of lesions and CC. Finally, Significant associations found for SNP’s −590 (IL-4)/- 573 (IL-6)/-592 (IL-10)/-819 (IL-10)/-509 (TGFB1) with CC, indicate that these polymorphisms are potential candidates for predicting the risk of development of CC, representing a risk allelic load for CC and can be used as a biomarker of susceptibility to this disease. These markers could be useful for early prevention programs and be clinically relevant as immunological biomarkers in patients with cervical lesions.
CC: cervical cancer; CI: confidence interval; ETS family proteins: E26 transformationspecific family proteins; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; HPV: human papilloma virus; HR-HPV: high-risk human papillomavirus; IFNG: Interferon-gamma; IL-10: interleukin 10; IL-4: interleukin 4; IL-6: interleukin 6; MAF: minor allele frequency; mRNA: messenger RNA; NCL: non-cervical lesions; OR: odds ratio; PBMC: peripheral blood mononuclear cells; PCR: polymerase chain reaction; SNP: single nucleotide polymorphism; STD: sexually transmitted disease; TGFB1: transforming growth factor beta 1; Th: T helper cells; TNF: tumor necrosis factor.
Arbyn M, Castellsagué X, de Sanjosé S, Bruni L, Saraiya M, Bray F, et al. Worldwide burden of cervical cancer in 2008. Ann Oncol. 2011;22(12):2675–86.
INEGI. Estadísticas de mortalidad. Cubos dinámicos y CONAPO 2012. Proyecciones de la población de México 2012–2050. 2012.
Insinga RP, Perez G, Wheeler CM, Koutsky LA, Garland SM, Leodolter S, et al. Incident cervical HPV infections in young women: transition probabilities for CIN and infection clearance. Cancer Epidemiol Biomarkers Prev. 2011;20:287–96.
Conesa-Zamora P. Immune responses against virus and tumor in cervical carcinogenesis: treatment strategies for avoiding the HPV-induced immune escape. Gynecol Oncol. 2013;131(2):480–8.
Bermúdez-Morales VH, Gutierrez LX, Alcocer-Gonzalez JM, Burguete A, Madrid-Marina V. Correlation between IL-10 gene expression and HPV infection in cervical cancer: a mechanism for immune response escape. Cancer Invest. 2008;26(10):1037–43.
Torres-Poveda K, Burguete-García AI, Cruz M, Martínez-Nava GA, Bahena-Román M, Ortíz-Flores E, et al. The SNP at −592 of human IL-10 gene is associated with serum IL-10 levels and increased risk for human papillomavirus cervical lesion development. Infect Agent Cancer. 2012;7(1):32–40.
Basavaraju U, Shebl FM, Palmer AJ, Berry S, Hold GL, El-Omar EM, et al. Cytokine gene polymorphisms, cytokine levels and the risk of colorectal neoplasia in a screened population of Northeast Scotland. Eur J Cancer Prev. 2015;24(4):296–304.
Zhang X, Zhang L, Tian C, Yang L, Wang Z. Genetic variants and risk of cervical cancer: epidemiological evidence, meta-analysis and research review. BJOG. 2014;121(6):664–74.
Muñoz N, Castellsagué X, de González AB, Gissmann L. Chapter 1: HPV in the etiology of human cancer. Vaccine. 2006;24 Suppl 3:S3/1–10.
Louie KS, de Sanjose S, Diaz M, Castellsagué X, Herrero R, Meijer CJ, et al. Early age at first sexual intercourse and early pregnancy are risk factors for cervical cancer in developing countries. Br J Cancer. 2009;100(7):1191–7.
Magnusson PK, Lichtenstein P, Gyllensten UB. Heritability of cervical tumours. Int J Cancer. 2000;88(5):698–701.
Alcocer-González JM, Berumen J, Taméz-Guerra R, Bermúdez-Morales V, Peralta-Zaragoza O, Hernández-Pando R, et al. In vivo expression of immunosuppressive cytokines in human papillomavirus-transformed cervical cancer cells. Viral Immunol. 2006;19(3):481–91.
Díaz-Benítez CE, Navarro-Fuentes KR, Flores-Sosa JA, Juárez-Díaz J, Uribe-Salas FJ, Román-Basaure E, et al. CD3zeta expression and T cell proliferation are inhibited by TGF-beta1 and IL-10 in cervical cancer patients. J Clin Immunol. 2009;29(4):532–44.
Piersma SJ. Immunosuppressive tumor microenvironment in cervical cancer patients. Cancer Microenviron. 2011;4(3):361–75.
Torres-Poveda K, Bahena-Román M, Madrid-González C, Burguete-García AI, Bermúdez-Morales VH, Peralta-Zaragoza O, et al. Role of IL-10 and TGF-β1 in local immunosuppression in HPV-associated cervical neoplasia. World J Clin Oncol. 2014;5(4):753–63.
Shamran HA, Hamza SJ, Yaseen NY, Al-Juboory AA, Taub DD, Price RL, et al. Impact of single nucleotide polymorphism in IL-4, IL-4R genes and systemic concentration of IL-4 on the incidence of glioma in Iraqi patients. Int J Med Sci. 2014;11(11):1147–53.
Tsai MH, Chen WC, Tsai CH, Hang LW, Tsai FJ. Interleukin-4 gene, but not the interleukin-1 beta gene polymorphism, is associated with oral cancer. J Clin Lab Anal. 2005;19(3):93–8.
Zheng Z, Li X, Li Z, Ma XC. IL-4 -590C/T polymorphism and susceptibility to liver disease: a meta-analysis and meta-regression. DNA Cell Biol. 2013;32(8):443–50.
Lim WY, Chen Y, Ali SM, Chuah KL, Eng P, Leong SS, et al. Polymorphisms in inflammatory pathway genes, host factors and lung cancer risk in Chinese female never-smokers. Carcinogenesis. 2011;32(4):522–9.
Castro FA, Haimila K, Sareneva I, Schmitt M, Lorenzo J, Kunkel N, et al. Association of HLA-DRB1, interleukin-6 and cyclin D1 polymorphisms with cervical cancer in the Swedish population—A candidate gene approach. Int J Cancer. 2009;125(8):1851–8.
Veerapaneni P, Kirma N, Nair HB, Hammes LS, Hall KL, Tekmal RR. Elevated aromatase expression correlates with cervical carcinoma progression. Gynecol Oncol. 2009;114:496–500.
Wei LH, Kuo ML, Chen CA, Cheng WF, Cheng SP, Hsieh FJ, et al. Interleukin-6 in cervical cancer: the relationship with vascular endothelial growth factor. Gynecol Oncol. 2001;82:49–56.
Shi TY, Zhu ML, He J, Wang MY, Li QX, Zhou XY, et al. Polymorphisms of the Interleukin 6 gene contribute to cervical cancer susceptibility in Eastern Chinese women. Hum Genet. 2013;132:301–12.
Grimm C, Watrowski R, Baumühlner K, Natter C, Tong D, Wolf A, et al. Genetic variations of interleukin-1 and −6 genes and risk of cervical intraepithelial neoplasia. Gynecol Oncol. 2011;121(3):537–41.
Bennermo M, Held C, Stemme S, Ericsson CG, Silveira A, Green F, et al. Genetic predisposition of the interleukin-6 response to inflammation: implications for a variety of major diseases? Clin Chem. 2004;50(11):2136–40.
Kong SY, Lee HL, Eom HS, Park WS, Yun T, Kim HJ, et al. Reference intervals for circulating angiogenic cytokines. Clin Chem Lab Med. 2008;46(4):545–50.
Kingo K, Rätsep R, Kõks S, Karelson M, Silm H, Vasar E, et al. Influence of genetic polymorphisms on interleukin-10 mRNA expression and psoriasis susceptibility. J Dermatol Sci. 2005;37:111–3.
Ni J, Ye Y, Teng F, Wu Q. Interleukin 10 polymorphisms and cervical cancer risk: a meta-analysis. Int J Gynecol Cancer. 2013;23(1):126–33.
Steinke JW, Barekzi E, Hagman J, Borish L. Functional Analysis of −571 IL-10 promoter polymorphism reveals a repressor element controlled by Sp1. J Immunol. 2004;173:3215–22.
Stanczuk GA, Sibanda EN, Perrey C, Chirara M, Pravica V, Hutchinson IV, et al. Cancer of the uterine cervix may be significantly associated with a gene polymorphism coding for increased IL-10 production. Int J Cancer. 2001;94(6):792–4.
Ding Q, Shi Y, Fan B, Fan Z, Ding L, Li F, et al. The Interleukin-10 Promoter Polymorphism rs1800872 (−592C.A), Contributes to Cancer Susceptibility: Meta-Analysis of 16 785 Cases and 19 713 Controls. PLoS ONE. 2013;8(2):e57246.
Singh H, Jain M, Sachan R, Mittal B. Association of TNFA (−308G > A) and IL-10 (−819C > T) promoter polymorphisms with risk of cervical cancer. Int J Gynecol Cancer. 2009;19(7):1190–4.
Yu Z, Liu Q, Huang C, Wu M, Li G. The interleukin 10–819C/T polymorphism and cancer risk: a HuGE review and meta-analysis of 73 studies including 15,942 cases and 22,336 controls. OMICS. 2013;17(4):200–14.
Xue H, Lin B, An J, Zhu Y, Huang G. Interleukin-10-819 promoter polymorphism in association with gastric cancer risk. BMC Cancer. 2012;12:102.
Wang J, Ding Q, Shi Y, Cao Q, Qin C, Zhu J, et al. The interleukin-10-1082 promoter polymorphism and cancer risk: a meta-analysis. Mutagenesis. 2012;27(3):305–12.
Rees L, Wood N, Gillespie K, Lai K, Gaston K, Mathieson P. The interleukin-10-1082 G/A polymorphism: allele frequency in different populations and functional significance. Cell Mol Life Sci. 2002;59(3):560–9.
Matsumoto K, Oki A, Satoh T, Okada S, Minaguchi T, Onuki M, et al. Interleukin-10–1082 gene polymorphism and susceptibility to cervical cancer among Japanese women. Jpn J Clin Oncol. 2010;40:1113–6.
Poli F, Nocco A, Berra S, Scalamogna M, Taioli E, Longhi E, et al. Allele frequencies of polymorphisms of TNFA, IL-6, IL-10 and IFNG in an Italian Caucasian population. Eur J Immunogenet. 2002;29:237–40.
Barbisan G, Pérez LO, Contreras A, Golijow CD. TNF-α and IL-10 promoter polymorphisms, HPV infection, and cervical cancer risk. Tumor Biol. 2012;33:1549–56.
Shah R, Hurley CK, Posch PE. A molecular mechanism for the differential regulation of TGF-beta1 expression due to the common SNP -509C-T (c. -1347C > T). Hum Genet. 2006;120(4):461–9.
Singh H, Jain M, Mittal B. Role of TGF-beta1 (−509C > T) promoter polymorphism in susceptibility to cervical cancer. Oncol Res. 2009;18(1):41–5.
Ramos-Flores C, Romero-Gutiérrez T, Delgado-Enciso I, Maldonado GE, Plascencia VM, Vazquez-Vuelvas OF, et al. Polymorphisms in the genes related to angiogenesis are associated with uterine cervical cancer. Int J Gynecol Cancer. 2013;23(7):1198–204.
Liu L, Yang X, Chen X, Kan T, Shen Y, Chen Z, et al. Association between TNF-α polymorphisms and cervical cancer risk: a meta-analysis. Mol Biol Rep. 2012;39:2683–8.
Ding B, Fu S, Wang M, Yue C, Wang W, Zhou D, et al. Tumor necrosis factor α -308 G > A polymorphisms and cervical cancer risk: a meta-analysis. Int J Gynecol Cancer. 2012;22(2):213–9.
Pan F, Tian J, Ji CS, He YF, Han XH, Wang Y, et al. Association of TNF-α-308 and −238 Polymorphisms with Risk of Cervical Cancer: A Meta-analysis. Asian Pac J Cancer Prev. 2012;13(11):5777–83.
Jin Y. Association of Single Nucleotide Polymorphisms in Tumor Necrosis Factor-Alpha with Cervical Cancer Susceptibility. Cell Biochem Biophys. 2015;71(1):77–84.
Zhang H-L, Zhang YJ. A systemic assessment of the association between tumor necrosis factor alpha 308G/A polymorphism and risk of cervical cancer. Tumor Biol. 2013;34(3):1659–65.
Stanczuk GA, Sibanda EN, Tswana SA, Bergstrom S. Polymorphism at the −308- promoter position of the tumor necrosis factor-alpha (TNF-alpha) gene and cervical cancer. Int J Gynecol Cancer. 2003;13:14853.
Govan VA, Constant D, Hoffman M, Williamson AL. The allelic distribution of −308 Tumor Necrosis Factor-alpha gene polymorphism in South african women with cervical cancer and control women. BMC Cancer. 2006;6:24.
Ivansson EL, Magnusson JJ, Magnusson PK, Erlich HA, Gyllensten UB. MHC loci affecting cervical cancer risk: distinguishing the effects of HLA-DQB1 and non-HLA genes TNF, LTA, TAP1 and TAP2. Genes Immun. 2008;9(7):613–23.
Calhoun ES, McGovern RM, Janney CA, Cerhan JR, Iturria SJ, Smith DI, et al. Host genetic polymorphism analysis in cervical cancer. Clin Chem. 2002;48(8):1218–24.
Baena A, Leung JY, Sullivan AD, Landires I, Vasquez-Luna N, Quiñones-Berrocal J, et al. TNF-alpha promoter single nucleotide polymorphisms are markers of human ancestry. Genes Immun. 2002;3:482–7.
Kim K, Cho SK, Sestak A, Namjou B, Kang C, Bae SC. Interferon-gamma gene polymorphisms associated with susceptibility to systemic lupus erythematosus. Ann Rheum Dis. 2010;69(6):1247–50.
Li CJ, Dai Y, Fu YJ, Tian JM, Li JL, Lu HJ, et al. Correlations of IFN-γ genetic polymorphisms with susceptibility to breast cancer: a meta-analysis. Tumour Biol. 2014;35(7):6867–77.
Song SH, Lee JK, Lee NW, Saw HS, Kang JS, Lee KW. Interferon-gamma (IFN-gamma): a possible prognostic marker for clearance of high-risk human papillomavirus (HPV). Gynecol Oncol. 2008;108(3):543–8.
Lages EL, Belo AV, Andrade SP, Rocha MÂ, de Freitas GF, Lamaita RM, et al. Analysis of systemic inflammatory response in the carcinogenic process of uterine cervical neoplasia. Biomed Pharmacother. 2011;65(7):496–9.
Ali KS, Ali HY, Jubrael JM. Concentration levels of IL-10 and TNFα cytokines in patients with human papilloma virus (HPV) DNA+ and DNA− cervical lesions. J Immunotoxicol. 2012;9(2):168–72.
Chechlinska M, Kowalska M, Kaminska J. Cytokines as potential tumour markers. Expert Opin Med Diagn. 2008;2(6):691–711.
Arany I, Grattendick KG, Tyring SK. Interleukin-10 induces transcription of the early promoter of human papillomavirus type 16 (HPV16) through the 5′-segment of the upstream regulatory region (URR). Antiviral Res. 2002;55(2):331–9.
The authors wish to thank Juraj Lord for manuscript assistance. The authors wish to thank all CAPASAM members and the administrative personnel of INSP. This work was submitted in partial fulfillment of the requirements for the PhD. degree of Torres-Poveda KJ from the Doctoral Program in Public Health of the School of Public Health of Mexico. Torres-Poveda thanks the Ministry of Foreign Affairs of Mexico for the PhD fellowship.
This work was supported by the Instituto Nacional de Salud Pública, Mexico and grants from Consejo Nacional de Ciencia y Tecnología (MX) CONACyT-FONSEC SSA/IMSS/ISSSTE-2008-C01-87701, CONACyT-Fondo E0013 APOYO COMPLEMENTARIO CATEDRAS-2014-C01-245520, CONACyT-FONSEC SSA/IMSS/ISSSTE-2014-C01-234149 and CONACyT-Fondo APOYOS COMPLEMENTARIOS PARA LA ADQUISICIÓN DE EQUIPO CIENTÍFICO 2013–205707, Mexico.
Availability of data and materials
SNPs studied in this investigation were validated SNPs for frequency or for utilization in the HAPmap Project and SNPs reported in public SNPper.URL: http//snpper.chip.org. Likewise, the SNP_ID of each SNP are included in the Genotyping section within of the Methods section.
TPK participated in the design of the study, the collection of clinical materials, carried out the genotyping and cytokines analysis in serum, the statistical analysis and drafted the manuscript. BGAI participated in the design of the study and the statistical analysis. BRM participated in the cytokines analysis in serum. MMR and ZDMA participated in the processing of samples. LEG, DRK and CD participated in the selection of study population, gynecological sampler and patient care. PZO and BMVH participated in the analysis of the results. GCA and MMV participated in the analysis of the results and drafted the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was conducted according to the principles expressed in the Declaration of Helsinki. It was approved by the Research, Ethics and Biosafety Committees at INSP (CI814 and CI1289 No. 1624) and INCan (INCan/CC/326/10CB/609). Written informed consent was obtained from all participants.
About this article
Cite this article
Torres-Poveda, K., Burguete-García, A.I., Bahena-Román, M. et al. Risk allelic load in Th2 and Th3 cytokines genes as biomarker of susceptibility to HPV-16 positive cervical cancer: a case control study. BMC Cancer 16, 330 (2016). https://doi.org/10.1186/s12885-016-2364-4