- Research article
- Open Access
Analysis for interaction between interleukin-35 genes polymorphisms and risk factors on susceptibility to coronary heart disease in the Chinese Han population
BMC Cardiovascular Disorders volume 21, Article number: 6 (2021)
The relationship between IL-35 genes polymorphism and susceptibility to coronary heart disease has not been tested in the largest Han population in China. The aim of this study was to explore the effect of single nucleotide polymorphisms (SNPs) of interleukin-35 (IL-35) genes and its relationship with environment on the risk of coronary heart disease (CHD).
We performed Hardy–Weinberg equilibrium test on the control group. The relationship between the four SNPs of IL-35 genes and the risk of coronary heart disease was studied by multivariate logistic regression. The best interaction was identified with generalized multifactor dimensionality reduction (GMDR). Logistic regression was used for investigation on association between four SNPs and CHD risk.
Logistic regression analysis showed that the C allele of rs428253 and the G allele of rs2243115 were independently correlated with increased risk of CHD, and adjusted ORs (95% CI) were 1.91 (1.28–2.64) and 1.80 (1.30–2.23), respectively. However, there was no significant association between CHD and rs4740 or rs568408. GMDR model indicated a best model for CHD risk consisted of rs428253 and current smoking, which scored 10/10 for both the sign test and cross-validation consistency (p = 0.010). Therefore, this overall multi-dimensional model had the highest cross-validation consistency, regardless of how the data were divided. This provided an evidence of gene–environment interaction effects. We also found that current smokers with rs428253-GC/CC genotype have the highest CHD risk, compared to never smokers with rs428253-GG genotype, OR (95% CI) = 3.04 (1.71–4.41), after adjustment for age, gender, hypertension, T2DM and alcohol consumption status.
In this study, the C allele of rs428253 and the G allele of rs2243115, and the interaction rs428253 and current smoking were correlated with increased risk of CHD.
In recent years, coronary heart disease (CHD) has become the main cause of the incidence rate and mortality rate in developed China and developing countries . In recent years, the incidence and mortality rate of CHD in China has increased rapidly, and the age of onset is younger. In young people (below 40 years old), about 700,000 people die from CHD every year [2, 3]. The etiology of CHD is very complex, and risk of CHD is not only affected by conventional environmental factors, but also by genetic variation . Previous study has reported many environmental risk factors related to coronary heart disease, including blood lipid concentration, blood pressure, smoking, diabetes and so on . In addition, inflammation also was reported associations with atherosclerosis development, which was an important risk factor for CHD .
Previous studies have reported that some inflammatory cytokine gene polymorphisms are associated with the risk of CHD, including cytokines in the IL-12 family [7, 8]. IL-35 is an anti-inflammatory cytokine and a member of IL-12 family. It is a heterodimer composed of p35 (IL-12A) and EBI3 subunits [9, 10]. Previous studies have reported that this genetic variation is related to the susceptibility of Alzheimer's disease and Graves' disease [11, 12], while EBI3 genetic variation may affect the risk of allergic rhinitis (AR) and tuberculosis in Chinese [13, 14]. Recently, only two studies [15, 16] have also studied the relationship between IL-35 genes polymorphism and CHD susceptibility. As far as we know, the relationship between IL-35 genes polymorphism and susceptibility to coronary heart disease has not been tested in the largest Han population in China. Lin et al.  performed a case–control study for Chinese population, but these participants were all Chinese Zhuang. In addition, CHD is a multifactorial disease caused by both genetic and environmental factors  and gene–environment interactions . The aim of this study was to evaluate the influence of SNPs within IL-35 gene and their interaction with environment on susceptibility to CHD.
The population of current study was composed of 921 CHD patients and 926 age- and gender-matched controls. All CHD patients were recruited from the First Naval Hospital of Southern Theater Command. CHD was defined according to the World Health Organization criteria . Participants with the following diseases will be excluded from the study cohort: heart related diseases, autoimmune diseases, chronic inflammatory diseases. The control participants were all the subjects who came to our hospital for routine occupational physical examination and voluntarily participated in this study. And those with CHD, hypertension and history of CHD were excluded from the control group. All the subjects were Han people, and there was no genetic relationship between them. Informed consent was obtained from each participant. In the study, we collected general demographic data and physical examination data of the subjects by questionnaire survey and physical measurement. The questionnaire includes general demographic information, life-style, smoking and alcohol consumption. This study has been approved by ethics committee of the First Naval Hospital of Southern Theater Command.
Extraction of genomic DNA and genotyping
We selected 4 SNPs according to the following criteria from dbSNP algorithm (http://www.ncbi.nlm.nih.gov/projects/SNP): firstly, the MAF > 5% in the database; secondly, the relationship between SNPs and CHD was not verified in the previous studies. Genomic DNA from whole blood containing EDTA was isolated strictly following the instructions of the manufacturer. Four SNPs including rs428253, rs4740, rs2243115 and rs568408 were selected. Genotyping for rs2243115 and rs568408 was performed using polymerase chain reaction (PCR) and following restriction fragment length polymorphism (RFLP). Genotyping for rs428253 and rs4740 using TaqMan genotyping assays on an ABI Prism 7900HT Fast Real-Time PCR System according to the manufacturer’s instructions (Applied Biosystems, Foster City, CA, USA). The primers and assays used for genotyping were listed in Table 1.
Hardy Weinberg equilibrium (HWE) was tested using SNPstats (http://bioinfo.iconcologia.net/SNPstats). Chi-square test was used to compare the distribution of alleles, genotypes and categorical variables among groups. The mean ± standard deviation (SD) was used to represent the continuous variables of normal distribution, and Student t test was used to compare the differences between the two groups. Generalized multifactor dimensionality reduction (GMDR)  was used to determine the optimal interaction combination of four SNPs of IL-35 genes and environmental risk factors. Logistic regression was used to analyze the relationship between four SNPs and the risk of CHD. When p value is less than 0.05, the statistical significance is significant.
In this study, a total of 1847 participants were enrolled, including 921 patients with CHD and 926 normal controls. Table 2 showed the different demographic characteristics in both case and control group. The mean age of all participants is 60.9 ± 10.8 years old. There are no statistically significant differences between case and control group in males, BMI and age (p > 0.05). Compared with control group, the mean of C-reactive protein (CRP) levels, the rates for hypertension and T2DM of case group are significantly different (all p < 0.05). In addition, the case group had higher alcohol consumption and smoking rates than control group.
All genotypes of the four SNPs follow by HWE distribution in the control group (all p values in HWE testing are greater than 0.05). The C allele frequency of rs428253 gene is 20.5% in the normal control group and 30.2% in patients with CHD. The minor (G allele) allele frequencies of rs2243115 gene are 26.3% in patients with CHD and 17.8% in normal controls. Logistic regression analysis showed that both the C allele of rs428253 and the G allele of rs2243115 were independently related to the increased risk of CHD. The adjusted OR (95% CI) was 1.91 (1.28–2.64) and 1.80 (1.30–2.23), respectively. However, there was no significant correlation between CHD and rs4740 or rs568408 (Table 3).
Further investigation on gene–gene and gene–environment interaction was tested using the GMDR methods, and the cross-validation consistency and testing accuracy were calculated. With covariate adjustments, the combination consist of rs428253 and current smoking is the best model, the cross-validation consistency of which was 10/10, and testing balanced accuracy of which was 0.632 (p = 0.010, Table 4). The fore-mentioned results provided an evidence for interaction effects between rs428253 and current smoking. But all gene–gene interaction combinations obtained from GMDR test have no statistical significance. We also conducted stratified analysis for interaction effect between rs428253 and current smoking using logistic regression. We found that current smokers with rs428253-GC/CC genotype have the highest CHD risk, compared to never smokers with rs428253-GG genotype, OR (95% CI) = 3.04 (1.71–4.41) (Fig. 1).
The current study indicated that both the C allele of rs428253 and the G allele of rs2243115 were correlated with increased risk of CHD, but rs4740 or rs568408 were not statistically associated with CHD risk. The rs428253 belongs to EBI3 gene, which located on chromosome 19p13.3. The rs2243115 belongs to IL-12A subunit, which is located on chromosome 3q25.33. Previously, several studies have investigated the association between the two SNPs and inflammatory diseases [11, 12]. So far, there are few studies about the effect of IL-35 polymorphism on cardiovascular disease (CVD). Zhang et al.  confirmed that EBI3-rs428253 polymorphism is related to the chronic rhinosinusitis and allergic rhinitis risk reduction. So far, only two studies [15, 16] have confirmed the relationship between IL-35 genes and CHD risk, but the conclusions of the two studies are inconsistent. A study  for Mexico populations suggested a relationship of the EBI3 SNPs with IL-35 levels, and EBI3-rs428253 and IL-12a-rs2243115 gene polymorphisms play an important role in the mechanism of CHD risk reduction. However, a recent case–control study in the Chinese population  suggested a statistical correlation between the EBI3 rs428253 mutation genotype and the risk of CHD, but there was no statistical significance between the mutation genotype of IL-12A rs2243115 and the risk of CHD. The study also found that there was no significant difference in the level of IL-35 between different genotypes in the healthy control group. Although two case–control studies have been performed previously, but the two studies concluded inconsistent results, and just one study was performed for Chinese Zhuang population, but no study focused on Chinese Han population, which was the largest race in China. So, the different race may lead to different results. The biological mechanism of the association between IL-35 genes and CHD susceptibility is not well established, previous studies suggested that EBI3-rs428253 may be involved in the modification of LEF1 binding site [15, 21], and play an important role in granulocyte proliferation and differentiation . SNP-rs428253 in EBI3 gene is related to the occurrence and development of CHD, which may be caused by regulating β—Catenin pathway and Treg pathway rather than by influencing the production of IL-35, because in the IL-35 levels was not different among different genotype of rs428253 in controls of study by Lin et al. . In this study, we found that rs4740 and rs568408 was not associated with CHD risk, which was consistent with study by Lin et al. . Although Posadas-Sánchez et al.  found that IL-35 genes rs4740 and rs568408 were not associated with CHD risk in Mexicans, but the two SNPs genotypes were correlated with IL-35 levels in healthy Mexicans.
The pathogenesis of CHD is very complex, and is not only affected by genetic factors and environmental factors independently, but also by the synergistic effect between them [23, 24]. We were all known that alcohol consumption and smoking were two main modifiable risk factors for CHD [25, 26]. In this study, the smoking and alcohol consumption rates were higher in cases than controls, which indicated that smoking and alcohol consumption were two risk factors for CHD. So, we performed a GMDR analysis for gene–smoking or alcohol consumption interaction, and the results suggested a significant interaction between rs428253 and current smoking. Current smokers with rs428253—GC/CC genotype have the highest CHD risk, compared to never smokers with rs428253—GG genotype. Environmental factors can cause phenotypic differences by influencing gene expression regulation. Therefore, the study of IL-35 gene–environment interaction is helpful for us to better understand the occurrence of CHD.
There were several limitations in this study. Firstly, the sample size is not large enough, so this is only a preliminary study of the polymorphism of this locus, the results obtained from current study need to be verified in the study with larger sample size and in different populations. Secondly, we just selected four SNPs in this study. In the future, we will investigate the relationship between multiple loci of IL-35 gene and CHD, in order to better understand the mechanism of CHD from the perspective of genetics. Thirdly, the G allele frequency of rs2243115 was higher than that in gene database, so the selection bias may exist in this study. Lastly, we did not obtain the information on lipid gram data, such as total cholesterol, LDL or perhaps Lp (a) values, so these lipid values could not be adjusted in the regression models. IL-35 levels were not measured in the study subjects, so we could not observe the difference of IL-35 levels between cases and controls.
The C allele of rs428253 and the G allele of rs2243115 are correlated with increased risk of CHD. We also find a significant interaction between rs428253 and current smoking associated with CHD risk. In addition, we also found the interaction between rs428253 and current smoking, which added more detailed mechanism for relationship between gene, environmental risk factors and CHD susceptibility.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
Coronary heart disease
Single nucleotide polymorphisms
Generalized multifactor dimensionality reduction
Polymerase chain reaction—based restriction fragment length polymorphism
Heusch G, Libby P, Gersh B, Yellon D, Bohm M, Lopaschuk G, Opie L. Cardiovascular remodelling in coronary artery disease and heart failure. Lancet. 2014;383(9932):1933–43.
Wang F, Xu CQ, He Q, Cai JP, Li XC, Wang D, Xiong X, Liao YH, Zeng QT, Yang YZ, Cheng X, Li C, Yang R. Genome-wide association identifies a susceptibility locus for coronary artery disease in the Chinese Han population. Nat Genet. 2011;43(4):345–9.
Yang G, Wang Y, Zeng Y, Gao GF, Liang X, Zhou M, Wan X, Yu S, Jiang Y, Naghavi M, Vos T, Wang H, Lopez AD, Murray CJ. Murray, rapid health transition in China, 1990–2010: findings from the Global Burden of Disease Study 2010. Lancet. 2013;381(9882):1987–2015.
Pranavchand R, Reddy BM. Current status of understanding of the genetic etiology of coronary heart disease. J Postgrad Med. 2013;59:30–41.
Batty GD, Shipley M, Smith GD, Kivimaki M. Long term risk factors for coronary heart disease and stroke: influence of duration of follow-up over four decades of mortality surveillance. Eur J Prev Cardiol. 2015;22(9):1139–45.
Cichoń N, Lach D, Dziedzic A, Bijak M, Saluk J. The inflammatory processes in atherogenesis. Pol Merkur Lekarski. 2017;42:125–8.
Posadas-Sánchez R, Pérez-Hernández N, Rodríguez-Pérez JM, Coral-Vázquez RM, Roque-Ramírez B, Llorente L, Lima G, Flores-Dominguez C, Villarreal-Molina T, Posadas-Romero C, Vargas-Alarcón G. Interleukin-27 polymorphisms are associated with premature coronary artery disease and metabolic parameters in the Mexican population: the genetics of atherosclerotic disease (GEA) Mexican study. Oncotarget. 2017;8:64459–70.
Tabrez S, Ali M, Jabir NR, Firoz CK, Ashraf GM, Hindawi S, Damanhouri GA, Nabil AM. A putative association of interleukin-10 promoter polymorphisms with cardiovascular disease. IUBMB Life. 2017;69:522–7.
Huang Y, Lin YZ, Shi Y, Ji QW. IL-35: a potential target for the treatment of atherosclerosis. Pharmazie. 2013;68:793–5.
Vignali DA, Kuchroo VK. IL-12 family cytokines: immunological playmakers. Nat Immunol. 2012;13:722–8.
Guo T, Yang S, Liu N, Wang S, Cui B, Ning G. Association study of interleukin-12A gene polymorphisms with Graves’ disease in two Chinese populations. Clin Endocrinol (Oxf). 2011;74:125–9.
Zhu XC, Tan L, Jiang T, Tan MS, Zhang W, Yu JT. Association of IL-12A and IL-12B polymorphisms with Alzheimer’s disease susceptibility in a Han Chinese population. J Neuroimmunol. 2014;274:180–4.
Zheng R, Liu H, Song P, Feng Y, Qin L, Huang X, Chen J, Yang H, Liu Z, Cui Z, Hu Z, Ge B. Epstein-Barr virus-induced gene 3 (EBI3) polymorphisms and expression are associated with susceptibility to pulmonary tuberculosis. Tuberculosis (Edinb). 2015;95:497–504.
Zhang Y, Duan S, Wei X, Zhao Y, Zhao L, Zhang L. Association between polymorphisms in FOXP3 and EBI3 genes and the risk for development of allergic rhinitis in Chinese subjects. Hum Immunol. 2012;73:939–45.
Posadas-Sánchez R, Pérez-Hernández N, Angeles-Martínez J, López-Bautista F, Villarreal-Molina T, Rodríguez-Pérez JM, Posadas-Romero C, Vargas-Alarcón G. Interleukin 35 polymorphisms are associated with decreased risk of premature coronary artery disease, metabolic parameters, and IL-35 levels: The Genetics of Atherosclerotic Disease (GEA) Study. Mediators Inflamm. 2017;2017:6012795.
Lin Y, Xue Y, Huang X, Lu J, Yang Z, Ye J, Zhang S, Liu L, Liu Y, Shi Y. Association between interleukin-35 polymorphisms and coronary heart disease in the Chinese Zhuang population: a case–control study. Coron Artery Dis. 2018;29(5):423–8.
Song C, Chang Z, Magnusson PK, Ingelsson E, Pedersen NL. Genetic factors may play a prominent role in the development of coronary heart disease dependent on important environmental factors. J Intern Med. 2014;275:631–9.
World Health Organization. Nomenclature and criteria for diagnosis of ischemic heart disease. Report of the joint international society and federation of cardiology/world health organization task force on standardization of clinical nomenclature. Circulation. 1979; 59:607–609.
Lou XY, Chen GB, Yan L, Ma JZ, Zhu J, Elston RC, Li MD. A generalized combinatorial approach for detecting gene-by-gene and gene-by-environment interactions with application to nicotine dependence. Am J Hum Genet. 2007;80:1125–37.
Zhang Y, Wang C, Zhao Y, Zhang L. Some polymorphisms in Epstein–Barr virus-induced gene 3 modify the risk for chronic rhinosinusitis. Am J Rhinol Allergy. 2013;27(2):91–7.
Skokowa J, Klimiankou M, Klimenkova O, Lan D, Gupta K, Hussein K, Carrizosa E, Kusnetsova I, Li Z, Sustmann C, Ganser A, Zeidler C, Kreipe HH, Burkhardt J, Grosschedl R, Welte K. Interactions among HCLS1, HAX1 and LEF-1 proteins are essential for G-CSF-triggered granulopoiesis. Nat Med. 2012;18:1550–9.
Sharma A, Hu XT, Napier TC, Al-Harthi L. Methamphetamine and HIV-1 Tat down regulate β-catenin signaling: implications for methampetamine abuse and HIV-1 co-morbidity. J Neuroimmune Pharmacol. 2011;6:597–607.
Song Y, Shan Z, Luo C, Kang C, Yang Y, He P, Li S, Chen L, Jiang X, Liu L. Glutathione S-transferase T1 (GSTT1) null polymorphism, smoking, and their interaction in coronary heart disease: a comprehensive meta-analysis. Heart Lung Circ. 2017;26(4):362–70.
Han Y, Dorajoo R, Ke T, Ayala B, Chang X, Khor CC, van Dam RM, Yuan JM, Koh WP, Liu J, Goh DY, Friedlander Y. Interaction effects between Paraoxonase 1 variants and cigarette smoking on risk of coronary heart disease in a Singaporean Chinese population. Atherosclerosis. 2015;240(1):40–5.
Yahagi K, Joner M, Virmani R. Impact of smoking on coronary heart disease: is there a smoker’s paradox? Coron Artery Dis. 2015;26(6):466–8.
Rehm J, Sempos CT, Trevisan M. Alcohol and cardiovascular disease–more than one paradox to consider. Average volume of alcohol consumption, patterns of drinking and risk of coronary heart disease—a review. J Cardiovasc Risk. 2003;10(1):15–20.
We appreciate the cooperation of the partners and staffs cooperated in this study.
Ethics approval and consent to participate
This study has been approved by ethics committee of the First Naval Hospital of Southern Theater Command. The written informed consent was obtained from all participants.
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.
About this article
Cite this article
Li, H., Liu, YX., Huang, JY. et al. Analysis for interaction between interleukin-35 genes polymorphisms and risk factors on susceptibility to coronary heart disease in the Chinese Han population. BMC Cardiovasc Disord 21, 6 (2021). https://doi.org/10.1186/s12872-020-01811-8
- Coronary heart disease
- Single nucleotide polymorphisms