IL1A polymorphisms is a risk factor for colorectal cancer in Chinese Han population: a case control study
BMC Cancer volume 19, Article number: 181 (2019)
Colorectal cancer (CRC) is one of the most common cancers worldwide, and genetic variations exert distinct roles in its pathogenesis. Single nucleotide polymorphisms (SNPs) in interleukin 1 alpha (IL1A) were reported to be correlated to the susceptibility of diverse cancers. The aim of this study was to assess the association of IL1A SNPs with the risk of colorectal cancer in a Chinese Han population.
To evaluate the correlation between IL1A polymorphisms and CRC risk, Agena MassARRAY platform was used for genotype determination among 248 CRC patients and 463 controls. The relationships between IL1A variants and CRC susceptibility were examined by logistic regression analysis. Stratified analysis was conducted for the association detection in males and females. Haplotype construction and analysis were applied to evaluate the potential relationship between the genetic block and the risk of CRC. SNP functional exploration was performed with available bioinformatics datasets.
After adjusting for age and gender, the “AA” genotype of rs2856838 exhibited a risk association with colorectal cancer in the recessive model (adjusted OR = 1.98, 95% CI: 1.05–3.72, p = 0.036). With stratified analysis, the recessive models of rs3783550 (OR = 2.17, 95% CI: 1.03–4.60, p = 0.043), rs2856838 (OR = 2.58, 95% CI: 1.13–5.87, p = 0.024), rs1609682 (OR = 2.20, 95% CI: 1.04–4.65, p = 0.040), and rs3783521 (OR = 2.13, 95% CI: 1.01–4.49, p = 0.048) revealed significant relationships between these variants and an increased CRC risk only in females. Bioinformatics analysis also revealed the putative functions of the selected SNPs.
This study demonstrated that rs2856838 could influence the susceptibility to CRC in Chinese Han population from northwest China. IL1A variants rs3783550, rs2856838, rs1609682, and rs3783521 were associated with CRC risk only in females.
Colorectal cancer (CRC) is one of the most common cancers worldwide, and it is the third leading cause of cancer-related deaths in both genders. In the past two decades, an increasing incidence of CRC was found in China as lifestyles and eating habits have become progressively changed [16, 33]. The etiology of CRC are multifactorial, and epidemiological evidence suggests that smoking cigarettes, unhealthy eating habits and lack of exercises are three well-known causal factors for CRC development [5, 25]. It has been confirmed that long-term exposure to these adverse habits may eventually increase the individual CRC risk, however, genetic aberrations, such as single-nucleotide polymorphisms (SNPs), play more important roles in the pathogenesis of CRC [2, 8, 23, 27, 28].
Accumulating evidence has demonstrated that most cancers develop from long-time chronic inflammation, including CRC, which could be attributed to two common inflammatory bowel diseases, namely, Crohn’s disease and ulcerative colitis [7, 18]. Inflammatory and tumour cells could produce cytokines and chemokines that facilitate tumour promotion and progression . Interleukin 1 (IL-1) is a family of pleiotropic cytokine with multiple functions in human innate inflammatory and immune responses . It has also been reported that IL-1 proteins are involved in tumour angiogenesis, invasion, proliferation, and metastasis . Together with IL-1β, IL-1α, encoded by IL1A (interleukin 1 alpha), is defined as an “alarm cytokine” that belongs to IL-1 cluster and plays dual roles in malignant tumour progression . Membrane-associated IL-1α could elevate immunogenicity of the tumour cell and exert positive effect on anti-tumour immune surveillance and tumour regression. However, secretable form of IL-1α in the microenvironment of tumour cells has been proved to facilitate tumour invasiveness and angiogenesis .
Furthermore, it is widely recognized that the polymorphisms of IL-1A is associated with diverse diseases. Gao et al. has reported an insertion/deletion (ins/del) polymorphism (rs3783553, TTCA) in IL1A that may contribute to hepatocellular cancer susceptibility, and revealed the regulatory function of the variation on IL1α expression by disrupting the binding sites for miR-122 and miR-378 . This functional polymorphism has also been demonstrated in gastric cancer and cervical carcinoma [13, 20, 31]. Additionally, the relationships between IL1A polymorphisms and the risk of papillary thyroid cancer, nasopharyngeal cancer and epithelial ovarian cancer have been reported as well [11, 29, 32].
A recent work has elucidated a decreased tumour expression of IL1A in colorectal adenocarcinoma, which indicated the potential role of IL1A in the etiology of CRC. However, few studies have examined the associations of IL1A polymorphisms with the risk of CRC. In our study, we investigated the effects of five IL1A variants (rs3783550, rs3783546, rs2856838, rs1609682, and rs3783521) on the susceptibility to CRC, which is supposed to provide more evidence for IL1A in CRC pathogenesis and contribute to early CRC risk estimation among the individuals of Chinese Han ancestry.
The current research involved a total of 248 CRC patients (143 males and 105 females) and 463 controls (265 males and 198 females) with unrelated Chinese Han ancestries. All CRC cases were diagnosed and confirmed by two independent pathological examinations. As for the eligible CRC cases selection, the individuals without other types of cancer, familial adenomatous polyposis, hereditary nonpolyposis colorectal cancer or undergone previous radiotherapy and chemotherapy were included. With regard to healthy controls, the subjects who had suffered from chronic or severe endocrine and metabolic diseases, digestive system disorders and malignant diseases were excluded from this study. The controls were polyp free individuals at recruitment. Moreover, the individuals with family colorectal disease and cancer history were excluded from control group as well.
By reading the previous publications of IL1A polymorphisms, we selected SNPs which could impact cancer risk, and searched the genetic data of them in dbSNP database (https://www.ncbi.nlm.nih.gov/snp/) and 1000 Genomes database (http://www.internationalgenome.org/). Only the SNPs whose minor allele frequency (MAF) beyond 5% in Asian populations were included in this study in order to achieve adequate statistical power. Finally, five candidate variations rs3783550, rs3783546, rs2856838, rs1609682, and rs3783521 in IL1A gene were selected for genotyping. Five millilitres venous blood was collected from the subjects when they were recruited in this study. Genomic DNA was extracted from the blood with the Whole Blood Genomic DNA Purification Kit (GoldMag, Xi’an, China). The PCR primers used in multiplexed PCR assay were designed by Agena Bioscience Assay Design Suite V2.0 software and are showed in Additional file 1: Table S1. In order to improve the PCR efficiency and ensure that the amplification primers, extension primers and extension products could be differed by their molecular weight, a tag 10-base sequence (5′-ACGTTGGATG-3′) was added to the 5′ end of each amplification primer. SNP genotyping was carried out with the usage of the Agena Nanodispenser (Agena Bioscience, San Diego, USA) and MassARRAY iPLEX platform (Agena Bioscience, San Diego, USA) . The procedure for SNP genotyping was described as follows. First, a first round PCR was performed to increase and concentrate the target genomic fragments containing the polymorphisms. Second, using the products obtained from the last step as the templates, only one ‘mass-modified’ nucleotide was added to the polymorphic locus at the end of the extension primer fragment in extension reaction. Third, the analyte mixture was dropped to a SpectroCHIP Array by Agena Bioscience Nanodispenser and further processed by the matrix-assisted laser desorption/ionization—time of flight (MALDI-TOF) mass spectrometry. With this method, the analytes could be differed according to their flight time in the system and the nucleotide at the each SNP locus was identified. Finally, the genotyping results were managed and presented in the Agena Bioscience TYPER version 4.0 software [10, 24].
SPSS 18.0 software (SPSS, Chicago, IL, USA) was employed for primary statistical analysis. The heterozygosity for the five SNPs in control cohort was assessed by the Fisher’s exact test to detect their compliance with Hardy-Weinberg equilibrium (HWE). Allele and genotype frequencies of each SNP in CRC cases and healthy controls were calculated and evaluated for the difference with Pearson’s χ2 test. Logistic regression analysis with adjustment for age and gender was applied to the odds ratio (OR) and 95% confidence interval (CI) assessment . Adjustment process is usually implemented in the statistical analysis to eliminate the influence of confounding factors, such as age and gender. In all statistical analysis, p ≤ 0.05 was regarded as statistical difference. The associations of the five candidate variants with the CRC risk were examined in codominant, dominant and recessive genetic models using SNPstats online software (https://www.snpstats.net/start.htm). Additionally, we used the Haploview software package, version 4.2 (http://analysis.bio-x.cn/myAnalysis.php) for linkage disequilibrium (LD) construction and haplotype evaluation . Using the genotyping results, we first constructed the LD pattern among genetic variants. The status of LD was analyzed using two parameters: r2 and D’. And D’ = 1 is defined as complete linkage disequilibrium. The allele combination of the variants in one LD block was considered as haplotype. We selected the haplotypes with frequency greater than 5% to evaluate whether the haplotypes were associated with CRC risk. These associations were detected using SHEsis software with available OR (95% CI) and p value.
Bioinformatics prediction and expression quantitative trait loci (eQTL) analysis
HaploReg v4.1 database (https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php) was applied for SNP functional annotation using bioinformatics methods. Expression quantitative trait loci (eQTL) analysis were carried out with the data available on the Genotype-Tissue Expression (GTEx) project (https://gtexportal.org/home/). The impact of each candidate SNP on IL1A expression was determined in several human tissues and p-values were provided for statistical significance assessment.
The demographic information of the enrolled CRC patients and healthy controls are listed in Additional file 1: Table S2. After applicable statistical evaluation, there was no significant difference in gender observed between patients and controls. The mean age [± standard deviation (SD)] of the patient group was 58.69 ± 12.77 years at diagnosis and of the controls was 50.65 ± 11.79 years at recruitment.
SNPs and the risk of CRC
Basic information and allele frequencies of the IL1A polymorphisms are presented in Table 1. None of the five SNPs was excluded because of the deviation from HWE at the 5% p level in the control population. The minor allele of each SNP was assumed a risk factor compared to the wild-type allele. Results of the genetic model analysis are showed in Table 2. The frequency of the homozygous “A/A” genotype of the IL1A SNP rs2856838 differed significantly between patients and controls (8.9% vs. 5.2%, p < 0.05). An association between allele “A/A” of rs2856838 and an increased risk of CRC was observed under the recessive model (adjusted OR = 1.98, 95% CI: 1.05–3.72, p = 0.036) after adjustments for age and gender.
In accordance with the stratified analysis by gender, we found that the IL1A variants rs3783550, rs2856838, rs1609682 and rs3783521 did not show any significant association with CRC risk in males, whereas the statistical significance was detected in females (Table 3). The frequency of homozygous “T/T” genotype of IL1A SNPs rs3783550 and rs1609682 differed significantly between patients and controls in females (rs3783550: 15.2% vs. 8.1%; rs1609682: 15.4% vs. 8.1%, p < 0.05), respectively. Moreover, rs3783550 and rs1609682 were associated with an increased risk of CRC based on the results of the recessive model (rs3783550: OR = 2.17, 95% CI: 1.03–4.60, p = 0.043; rs1609682: OR = 2.20, 95% CI: 1.04–4.65, p = 0.040). The frequency of homozygous “A/A” genotype of rs2856838 differed significantly between patients and controls in females (13.3% vs. 6.1%, p < 0.05) as well. The recessive genetic model also provided the evidence for the role of rs2856838 in CRC susceptibility among Chinese Han women (OR = 2.58, 95% CI: 1.13–5.87, p = 0.024). Additionally, rs3783521 homozygous genotype “G/G” distributed differently between the case and control group (15.2% vs. 8.1%, p < 0.05). Further genetic model analysis revealed that SNP rs3783521 was linked to an increased risk of CRC in females according to the recessive model (OR = 2.13, 95% CI: 1.01–4.49, p = 0.048). However, no statistically significant association was detected between SNP rs3783546 and CRC risk in Chinese Han women.
CRC and haplotypes on chromosome 2q13
Finally, five IL1A polymorphisms (rs3783550, rs3783546, rs2856838, rs1609682, rs3783521) were mapped to a 10-kb LD block and displayed three haplotypes with frequencies of more than 5% in our subjects (Table 4). In Fig. 1, the red squares of the IL1A LD block represented significant linkage between the five SNPs. The LD degree is displayed by standard color schemes with bright red for very strong LD (D’ = 1). According to the haplotype analysis, there were no statistically significant differences appearing between patients and controls among any of the IL1A haplotype frequencies.
SNP functional evaluation
In order to evaluate the possible function harbored by the five selected variants, we performed a bioinformatics analysis using HaploReg v4.1 database. All the variations were successfully predicted as regulatory SNPs with different biological functions (Additional file 1: Table S3). Furthermore, eQTL analysis uncovered remarkable relationships between these SNPs and IL1A expression in skin, testis or pituitary, which indicated the effects of different genotypes on gene expression (Additional file 1: Table S4).
In this study, we genotyped five polymorphisms of IL1A and evaluated their correlations with the risk of CRC in a Chinese Han population. Our results first showed that four SNPs (rs3783546, rs2856838, rs1609682, rs3783521) were associated with CRC susceptibility in Chinese Han females.
Interleukin-1 alpha (IL-1α) is a major member of the interleukin-1 (IL-1) family, which has been well established as pro-inflammatory cytokine cluster involved in tumorigenesis. Several studies have showed that IL-1α contributes to tumour invasion, tumour proliferation, tumour metastases and interaction between the host immune system and malignant tumour cells [9, 19]. Kasza has reported that IL-1α serves as a pleiotropic cytokine in human cancer progression as well as various immune and inflammatory responses, and the secreted IL-1α is an active form which induces the tissue damage and tumour growth together with other regulatory factors . Accordingly, it is biologically reasonable that functional IL1A polymorphisms may play causal roles in the development of cancer.
Increasing evidence has showed that IL1A polymorphisms are associated with several cancers . Fei Wang et al. has found that four IL1A SNPs, rs3783550, rs3783546, rs1609682, and rs3783521, were associated with a significantly increased risk of renal cell carcinoma . Xiaofeng Zeng et al. has demonstrated the relationship between the IL1A rs3783553 polymorphism and a significantly decreased risk of gastric cancer (OR = 0.48, 95% CI:0.26–0.90, p = 0.02) . Abazis-Stamboulieh D et al. has proved that the IL1A rs1800587 polymorphism was related to a higher risk of multiple myeloma development . Our data showed that rs2856838 was significantly associated with CRC susceptibility (p = 0.036). After stratification by gender, we found that the recessive model of rs3783550, rs2856838, rs1609682, and rs3783521 were associated with increased CRC risk in females. Moreover, it should be noted that the incidence rate of CRC differs significantly according to gender . Considering the predicted functions of the selected SNPs in our study and their influences on gene expression provided by database, we speculated that SNPs might affect oncogenic processes by altering gene expression and this process could be influenced by individual gender background, thus leading to different outcomes on CRC risk between males and females. Significant IL1A variants rs3783550, rs2856838, rs1609682 were located in the intron, and rs3783521 were in promoter. As both the intron and promoter play important roles in gene expression, the alterations of the SNPs might also influence gene via modulating the functions of these regions. In addition, gender may reflect the different levels of exposure to risk factors related to diet, occupational exposure, or lifestyle. However, due to the limitations of the collected data, it was only possible to study the effect of gender rather than the complete impacts of other factors contributing to the gene environment. Therefore, it is necessary to carry out a complete analysis on the interaction of these factors with genetic polymorphisms in CRC susceptibility to verify the results of this study.
Several limitations should be considered in this research. Firstly, the influence of the selected SNPs on IL1A gene expression did not verify in our population and the underlying mechanism upon this regulatory process still need to be further elucidated. Further research will be design to investigate the influence of the promising SNPs on intron and promoter functions in gene regulation, and the gender will be considered as an important factor in future research. Secondly, although variations of age and gender were considered during our analysis, some living habits were not evaluated in this study, such as alcohol consumption and cigarette smoking, which might facilitate the progression of CRC.
Our study first provided objective evidence that the IL1A genetic polymorphism rs2856838 might contribute to the risk of CRC in Chinese Han population. In addition, rs3783550, rs1609682, rs3783521, and rs2856838 were susceptible pathogenic variants for CRC among females. These results could suggest new targets for CRC risk assessment, prevention and prognosis.
Expression quantitative trait loci
Genotype-Tissue expression; SD: Standard deviation
Minor allele frequency
- MALDI-TOF mass spectrometry:
Matrix-assisted laser desorption/ionization—time of flight mass spectrometry
Single nucleotide polymorphisms
Abazis-Stamboulieh D, Oikonomou P, Papadoulis N, Panayiotidis P, Vrakidou E, Tsezou A. Association of interleukin-1A, interleukin-1B and interleukin-1 receptor antagonist gene polymorphisms with multiple myeloma. Leukemia & Lymphoma. 2007;48:2196–203.
An J, Zhao J, Zhang X, Ding R, Geng T, Feng T, Jin T. Impact of multiple alcohol dehydrogenase gene polymorphisms on risk of laryngeal, esophageal, gastric and colorectal cancers in Chinese Han population. Am J Cancer Res. 2015;5:2508–15.
Apte RN, Voronov E. Is interleukin-1 a good or bad ‘guy’ in tumor immunobiology and immunotherapy? Immunol Rev. 2008;222:222.
Bland JM, Altman DG. The odds ratio Bmj. 2000;320:1468.
Botteri E, Iodice S, Bagnardi V, Raimondi S, Lowenfels AB, Maisonneuve P. Smoking and colorectal Cancer: a meta-analysis Jama. J Am Med Assoc. 2008;102:996–7.
Burada F, Dumitrescu T, Nicoli R, Ciurea ME, Angelescu C, Mixich F, Ioana M. IL-1RN +2018T>C polymorphism is correlated with colorectal cancer. Mol Biol Rep. 2013;40:2851–7.
Coussens LM, Werb Z. Inflammation and cancer. Nature. 2002;420:860–7.
Duan X, et al. Genetic association of PLCE1, C11orf92-C11orf93, and NOC3L with colorectal cancer risk in the Han population. Tumor Biol. 2014;35:1813–7.
Elaraj DM, et al. The role of interleukin 1 in growth and metastasis of human cancer xenografts clinical Cancer research An official journal of the American association for. Cancer Res. 2006;12:1088.
Gabriel S, Ziaugra L, Tabbaa D (2009) SNP genotyping using the Sequenom MassARRAY iPLEX platform Current Protocols in Human Genetics Chapter 2:Unit 2.12.
Gao L, et al. Association between a functional insertion/deletion polymorphism in IL1A gene and risk of papillary thyroid carcinoma. Tumor Biol. 2014;35:3861.
Gao Y, et al. An insertion/deletion polymorphism at miRNA-122-binding site in the interleukin-1alpha 3′ untranslated region confers risk for hepatocellular carcinoma. Carcinogenesis. 2009;30:2064–9.
Huang J, Ni S, Li D, He Y. An insertion/deletion polymorphism at miRNA-122 binding site in the IL1A is associated with a reduced risk of cervical squamous cell carcinoma. Genetic testing and molecular biomarkers. 2015;19:331–4. https://doi.org/10.1089/gtmb.2015.0015.
Karjalainen J, et al. The IL1A genotype is associated with nasal polyposis in asthmatic adults. Allergy. 2003;58:393–6.
Kasza A. IL-1 and EGF regulate expression of genes important in inflammation and cancer. Cytokine. 2013;62:22–33.
Li M, Gu J. Changing patterns of colorectal cancer in China over a period of 20 years. World J Gastroenterol. 2005;11:4685–8.
Martinelli M, Scapoli L, Cura F, Rodia MT, Ugolini G, Montroni I, Solmi R. Colorectal cancer susceptibility: apparent gender-related modulation by ABCB1 gene polymorphisms. J Biomed Sci. 2014;21:89.
Moss SF, Blaser MJ. Mechanisms of disease: inflammation and the origins of cancer. Nat Clin Pract Oncol. 2005;2:90.
Pantschenko AG, et al. The interleukin-1 family of cytokines and receptors in human breast cancer: implications for tumor progression. Int J Oncol. 2003;23:269.
Pu Y, et al. Association of an insertion/deletion polymorphism in IL1A 3′-UTR with risk for cervical carcinoma in Chinese Han women. Hum Immunol. 2014;75:740–4.
Rider P, Carmi Y, Voronov E, Apte RN. Interleukin-1alpha. Semin Immunol. 2013;25:430–8. https://doi.org/10.1016/j.smim.2013.10.005.
Sakamoto K, et al. Inhibitor of kappaB kinase beta regulates gastric carcinogenesis via interleukin-1alpha expression. Gastroenterology. 2010;139:226.
Su Q, Wang Y, Zhao J, Ma C, Wu T, Jin T, Xu J. Polymorphisms of PRLHR and HSPA12A and risk of gastric and colorectal cancer in the Chinese Han population. BMC Gastroenterol. 2015;15. https://doi.org/10.1186/s12876-015-0336-9.
Thomas RK, et al. High-throughput oncogene mutation profiling in human cancer. Nat Genet. 2007;39:347–51.
Tsoi KK, Pau CY, Wu WK, Chan FK, Griffiths S, Sung JJ. Cigarette smoking and the risk of colorectal cancer: a meta-analysis of prospective cohort studies. Clinical Gastroenterology & Hepatology the Official Clinical Practice Journal of the American Gastroenterological Association. 2009;7:682.
Wang F, et al. IL1 genes polymorphism and the risk of renal cell carcinoma in Chinese Han population. Oncotarget. 2017;8:56021–9.
Wang N, et al. Multiple genetic variants are associated with colorectal cancer risk in the Han Chinese population. European journal of cancer prevention : the official journal of the European Cancer Prevention Organisation (ECP). 2015;24:1–5. https://doi.org/10.1097/cej.0000000000000012.
Yang H, Gao Y, Feng T, Jin TB, Kang LL, Chen C. Meta-analysis of the rs4779584 polymorphism and colorectal cancer risk. PLoS One. 2014;9:e89736.
Yang ZH, Dai Q, Zhong L, Zhang X, Guo QX, Li SN. Association of IL-1 polymorphisms and IL-1 serum levels with susceptibility to nasopharyngeal carcinoma. Mol Carcinog. 2011;50:208–14.
Yong Y, Lin H. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005;15:97–8.
Zeng XF, Li J, Li SB. A functional polymorphism in IL-1A gene is associated with a reduced risk of gastric cancer. Tumour biology the journal of the international society for Oncodevelopmental biology & Medicine. 2014;35:265–8.
Zhang Z, et al. A polymorphism at miRNA-122-binding site in the IL-1α 3'UTR is associated with risk of epithelial ovarian cancer. Familial Cancer. 2014;13:595.
Zhong R, et al. Genetic variations in the TGFβ signaling pathway, smoking and risk of colorectal cancer in a Chinese population. Carcinogenesis. 2013;34:936–42.
We are grateful to the patients and control individuals for providing blood samples. We also appreciate all editors and reviewers for their patience and suggestions to this work.
Availability of data and materials
The data analyzed during the current research can be available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Our research was approved by the Ethics Committee of the Second Affiliated Hospital of Xi’an Jiaotong University. All procedures performed in this study were in accordance with the ethical standards of the Ethics Committee of the Second Affiliated Hospital of Xi’an Jiaotong University and with the 1964 Helsinki declaration and its later amendments. Informed consents were obtained from all individual participants included in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primers used for identification of the IL1A polymorphisms. Table S2. Distributions of age and gender in CRC patients and controls. Table S3. Functional annotation of the selected variants provided by HaploReg 4.1. Table S4. eQTL analysis for the five IL1A SNPs provided by GTEx database. (DOCX 24 kb)
About this article
Cite this article
Ji, H., Lu, L., Huang, J. et al. IL1A polymorphisms is a risk factor for colorectal cancer in Chinese Han population: a case control study. BMC Cancer 19, 181 (2019). https://doi.org/10.1186/s12885-019-5395-9
- Single-nucleotide polymorphisms
- Stratification analysis
- Colorectal cancer
- Chinese Han population