- Research
- Open access
- Published:
Association of the interaction between interleukin-1β gene polymorphism and smoking status with the diabetic nephropathy risk in a Chinese Han population
Diabetology & Metabolic Syndrome volume 17, Article number: 101 (2025)
Abstract
Objectives
we aimed to evaluate the association of interleukin-1β (IL-1β) gene single nucleotide polymorphisms (SNPs) and its interaction with smoking status on diabetic nephropathy (DN) risk in a Chinese Han population.
Methods
The Hardy-Weinberg equilibrium (HWE) was tested by using SNPStats (https://www.snpstats.net/start.htm), which was also used for testing the relationship between four SNPs and DN risk and haplotype analysis. The SNP- SNP and gene- smoking interaction were verified by using generalized multifactor dimensionality reduction (GMDR) model.
Results
Logistic regression suggested that the DN risks of participants with rs16944- G allele were significantly higher than those with AA genotype, adjusted OR (95%CI) = 1.62 (1.24–2.01) for AG versus AA, 1.41 (0.75–2.12) for GG versus AA. Additionally, we also found that participants with rs3917356- T allele had an obviously higher DN risk than those with CC genotype, adjusted OR (95%CI) = 1.75 (1.34–2.19) for CT versus CC, 1.87 (1.23–2.54) for TT versus CC. GMDR model found a significant two-locus model (P = 0.011) including rs16944 and smoking. Compared with non- smokers with rs16944- AA genotype, smokers with rs1225404 AG or GG genotype had the highest DN risk after covariates adjustment, OR (95%CI) was 3.04 (1.98–4.12). We also found a haplotype containing rs1143634- T and rs3917356- T was associated with higher DN risk.
Conclusions
we found that the rs16944- G and rs3917356- T allele, interaction between rs16944 and smoking, haplotype containing rs1143634- T and rs3917356- T were all associated with increased DN risk.
Introduction
Diabetes nephropathy (DN) is one of the major complications of type 2 diabetes (T2DM) and a major risk factor for end-stage renal disease (ESRD) or chronic kidney disease (CKD) [1, 2]. The main clinical features of DN include renal failure, hypertension, massive proteinuria and persistent proteinuria [3]. In the past decade, the prevalence and incidence rate of DN among Chinese residents have increased sharply, affecting the health of about 150 million people [4]. Importantly, with the increase of incidence rate of type 2 diabetes in the Chinese population, the incidence rate of diabetes nephropathy continues to increase [5]. There are many factors that affect the risk of DN, including age, gender, dyslipidemia, smoking, race, and genetic background [6]. Previous studies [7, 8] have reported that candidate genes have an obvious impact on the genetic susceptibility of DN. Some studies [9, 10] suggest that genetic variations encoding inflammatory cytokines may have a significant impact on susceptibility to DN.
The interleukin-1 (IL-1) gene family produces three cytokines, including IL-1α, IL-1β and IL-1Ra. These molecules play crucial roles in various inflammatory responses by influencing antigen recognition and the function of lymphocytes [11]. IL-1β could bind to IL-1 receptor 1, leading to fever and immune activation, and is therefore known as one of the most prominent mediators of inflammation [11]. Previously, several studies have reported that IL-1β gene single nucleotide polymorphisms (SNPs) were associated with diabetic retinopathy [12]. Studies also suggested a correlation between IL-1β gene SNPs and DN risk in different populations [13, 14]. However, to date, just one study [15] has been performed in Chinese population. In addition, smoking has been suggested as an important environmental risk factor for DN [16, 17]. Several studies [18, 19] also indicated that gene- smoking interaction could influence DN susceptibility. But to date, no study was performed focused on the interaction effect between IL-1β gene and smoking. Therefore, we performed this case- control study to evaluate the impact of IL-1β gene polymorphism, its interaction with smoking status and potential haplotype among four SNPs on the risk of DN based on a Chinese Han population.
Materials and methods
Study population
All participants were recruited continuously from the affiliated Hospital of Nantong University between July 2017 and March 2023. A total of 828 subjects were enrolled. The average age of all participants was 51.7 ± 11.7 years, involving 412 cases (T2DM patients with DN) and 416 controls (T2DM patients without DN). T2DM patients were diagnosed according to the criteria recommended by World Health Origination [20]. DN patients were diagnosed based on a consensus established by the American Diabetes Association and the National Kidney Foundation [21, 22]. Those T2DM patients with severe system diseases, including cardiovascular disease (CVD), cancer, renal disease and historical exposure to radiation were excluded. The controls were matched to cases by sex, age (± 3 years) and nationality. The affiliated Hospital of Nantong University’s ethics committee approved all study protocols for this research.
SNP selection and genotyping methods
Four SNPs of the IL-1β gene were searched in the SNP database of the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/snp/), according to the following criteria: first, those SNPs of the IL-1β gene which were not detailed studied previously; second, the minor allele frequency (MAF) of SNPs was greater than 2%. At last, 4 SNPs of IL-1β gene were selected, including: rs3917356, rs16944, rs1143634 and rs1143627. All collected EDTA-processed whole blood samples were stored in a cryogenic refrigerator. We then extracted DNA from a 2 ml whole blood sample using the Genomic DNA Extraction Kit (Tiangen, Beijing, China), adhering to the provided manufacturer’s guidelines. Polymerase chain reactions (PCRs) were performed for genotyping using specific primers for four SNPs of IL-1β gene (Table 1).
Statistical analysis
This research utilized R4.4.1 software for all statistical evaluations. We computed the means and standard deviations of continuous variables that followed a normal distribution and compared them via the student t-test. For categorical variables, percentages were calculated and then compared using the χ2 test. The Hardy-Weinberg equilibrium (HWE) was tested by using SNPStats (https://www.snpstats.net/start.htm), which was also used for investigation the relationship between four SNPs and DN risk and haplotype analysis. The SNP-SNP and gene-smoking interactions were confirmed through the application of the generalized multifactor dimensionality reduction (GMDR) model [23].
Results
This research involved 828 participants, with an average age of 51.7 ± 11.7 years. The participants consist of 412 cases (T2DM patients with DN) and 416 controls (T2DM patients without DN). Table 2 displays the general and clinical data of all participants grouped by cases or controls. The distribution for gender, alcohol consumption and the means of age, BMI, duration of diabetes, TG, HDL-C, LDL-C was not significantly different between cases and controls. More smokers were found in the cases than in the controls. The means of FPG, SBP, DBP, TC, urea and eGFR were significantly higher in cases than controls.
All genotypes in control group conformed to Hardy-Weinberg equilibrium (all P values > 0.05). The frequencies of rs16944- allele (G) and genotype (AG or/ and GG) and rs3917356 allele (T) and genotype (CT or/ and TT) in T2DM patients with DN group were significantly higher than that in T2DM patients without DN group. Logistic regression suggested that the DN risk of participants with rs16944- G allele was obviously higher than those with AA genotype, adjusted OR (95%CI) = 1.62 (1.24–2.01) for AG versus AA, 1.41 (0.75–2.12) for GG versus AA. In addition, we also found that DN risk was higher in participants with rs3917356- T allele than those participants with CC genotype, adjusted OR (95%CI) = 1.75 (1.34–2.19) for CT versus CC, 1.87 (1.23–2.54) for TT versus CC (Table 3).
We employed the GMDR model to assess the effect of interactions between SNP-SNP and gene-smoking within four SNPs of the IL-1β gene, along with smoking status on the risk of DN (Table 4). No statistically significant SNP- SNP interaction was obtained among IL-1β gene SNPs. However, we found a significant two-locus model (p = 0.011) including rs16944 and smoking, which suggesting a potential gene–environment interaction between rs16944 and smoking on DN risk. After adjusting for covariates, smokers who possess the rs1225404 AG or GG genotype, in comparison to non-smokers with the rs16944-AA genotype, exhibited the greatest risk of DN, with an odds ratio (95% CI) of 3.04 (1.98–4.12) (Table 5).
Pairwise LD method calculated the D’ values among four SNPs of the IL-1β gene. The results showed a heavy LD between rs1143634 and rs3917356. Therefore, we performed haplotype analysis for rs1143634 and rs3917356 by using SHEsis software. We found that rs1143634- C and rs3917356- C haplotype frequency was the highest in both cases and control groups (0.4811 for T2DM with DN patients, 0.5689 for T2DM without DN patients). The haplotype containing rs1143634- T and rs3917356- T allele was associated with higher DN risk after covariate adjustment (Table 6).
Discussion
In this study, we evaluated the effect of IL-1β gene SNPs on DN risk. We found that increased DN risk was associated with both rs16944- G and rs3917356- T allele. Nevertheless, we did not find any significant correlation of rs1143634 and rs1143627 with DN susceptibility, after covariates adjustment. The polymorphism rs16944 in the IL-1β gene is linked to changes in both the transcription [24] of this gene and the production levels of the IL-1β protein [25]. Previously, an animal study [26] found a significant upregulation of IL-1β levels in diabetic animals in comparison with normal rats. Yuan et al. [15] conducted a study observing a cohort in China, which indicated that IL-1β is significantly implicated in the advancement of T2DM. They also found that its concentration escalates as DN worsens. Therefore, the level of IL-1β in human plasma can serve as an important indicator for determining the occurrence and development of DN, provide reference for early diagnosis and precise treatment of DN. Bi et al. [27] suggested that IL-1β was positively correlated with insulin resistance (IR) in non-diabetic hemodialysis patients, which means that IL-1β may be involved in the pathogenesis of IR. The previous case-control study in South Korean diabetes population suggested that there was a significant correlation between the T allele of rs16944 and the risk of DN [28]. However, the other two studies performed in white patients with diabetes have contradictory results with this study, they found no significant association between rs16944 and DN risk [29, 30]. Stefanidis et al. [14] conducted a case-control study at the University Hospital of Aachen and demonstrated that the rs16944 polymorphism in the IL1β gene was statistically associated with the onset of type 2 diabetic nephropathy and resulted in increased IL-1β synthesis. Hameed et al. [31] confirmed that there was a statistical correlation between the promoter polymorphism of cytokine IL-1 β gene and the regulation of transcription level and the susceptibility to nephropathy in diabetic patients. Wu et al. [32] performed a comprehensive meta-analysis and revealed a significant relationship between rs16944 polymorphism and the increased risk of developing DN. Previously, no study focused on relationship between rs3917356 and DN risk was performed. This study firstly verified the impact of rs3917356 and DN risk in Chinese population. We also found that rs1143634- T and rs3917356- T haplotype was associated with increased DN risk. The exact mechanism of IL-1β gene SNPs mediated progression of DN is still unclear. However, IL-1β gene polymorphisms could influence IL-1β levels, IR and which was associated with DN susceptibility [15, 27].
Numerous factors, including genetics, environmental influences and the interaction between them can influence susceptibility to DN. It is widely recognized that smoking significantly elevates the risk of developing DN [18, 19]. The exact mechanism of smoking mediated progression of DN is still unclear. Potential mechanisms include persistent sympathetic overactivity, atherosclerosis, oxidative stress and hyperlipidemia [33]. There may be a synergistic effect between smoking and hyperglycemia, which aggravates kidney deterioration through continuous sympathetic nerve activity, atherosclerosis, oxidative stress and hyperlipidemia. Our research revealed a higher prevalence of smoking among patients with diabetic nephropathy. Additionally, a notable interaction between the genetic marker rs16944 and tobacco use was observed. Maybe rs16944 and smoking can affect the risk of DN through different mechanisms. When both exist, they will significantly increase the risk of individual DN.
Although we obtained some novel results, several certain limitations existed in this study. Firstly, only four SNPs of IL-1 β gene were included in this study and more SNPs should be included in future research to explore the synergistic effect of multiple SNPs on the risk of DN. Secondly, this study was conducted in a single hospital of Nantong city, Jiangsu province in southeast of China, selection bias may be present to some extent. Thirdly, smoking status was self-reported, which may lead to misclassification and bias. Lastly, this study was a case-control study, which did not allow for causal inferences. So, the results obtained from this study should be verified by future studies explored the functional mechanisms underlying the observed genetic effects.
In conclusion, our research table confirms that there is a significant statistical correlation between IL-1 β gene rs16944-G and rs3917356-T alleles and DN risk in the Chinese population. In addition, there is a significant synergistic effect between rs16944 and smoking on the risk of DN. Although this study did not find an association between rs1143634 and DN risk, we found that the haplotype rs1143634-T and rs3917556-T was associated with increased DN risk.
Data availability
No datasets were generated or analysed during the current study.
Abbreviations
- OR:
-
Odds ratios
- CI:
-
Confidence interval
- IL-1β:
-
Interleukin-1β
- DN:
-
Diabetic nephropathy
- GMDR:
-
Generalized multifactor dimensionality reduction
- HWE:
-
Hardy-weinberg equilibrium
- T2DM:
-
Type 2 diabetes mellitus
- SNPs:
-
Single nucleotide polymorphism
- ESRD:
-
End-stage renal disease
- CKD:
-
Chronic kidney disease
References
Jadawji C, Crasto W, Gillies C, Kar D, Davies MJ, Khunti K, Seidu S. Prevalence and progression of diabetic nephropathy in South Asian, white European and African Caribbean people with type 2 diabetes: A systematic review and meta-analysis. Diabetes Obes Metab. 2019;21(3):658–73.
Demir Y, Ceylan H, Turkes C, Beydemir S. Molecular Docking and Inhibition studies of vulpinic, carnosic and Usnic acids on polyol pathway enzymes. J Biomol Struct Dyn. 2022;40(22):12008–21.
Merker L, Ebert T, Guthoff M, Isermann B. Nephropathy in diabetes. Exp Clin Endocrinol Diabetes. 2023;131(1–02):61–5.
Li Y, Ning Y, Shen B, Shi Y, Song N, Fang Y, Ding X. Temporal trends in prevalence and mortality for chronic kidney disease in China from 1990 to 2019: an analysis of the global burden of disease study 2019. Clin Kidney J. 2022;16(2):312–21.
Sun H, Saeedi P, Karuranga S, Pinkepank M, Ogurtsova K, Duncan BB, Stein C, Basit A, Chan JCN, Mbanya JC, Pavkov ME, Ramachandaran A, Wild SH, James S, Herman WH, Zhang P, Bommer C, Kuo S, Boyko EJ, Magliano DJ. IDF diabetes atlas: global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res Clin Pract. 2022;183:109119.
Natesan V, Kim SJ. Diabetic Nephropathy - a review of risk factors, progression, mechanism, and dietary management. Biomol Ther (Seoul). 2021;29(4):365–72.
Brennan E, McEvoy C, Sadlier D, Godson C, Martin F. The genetics of diabetic nephropathy. Genes. 2013;4(4):596–619.
Lee SH, Lee TW, Ihm CG, Kim MJ, Woo JT, Chung JH. Genetics of diabetic nephropathy in type 2 DM: candidate gene analysis for the pathogenic role of inflammation. Nephrology (Carlton). 2005;10 Suppl: S32-6.
El Helaly RM, Elzehery RR, El-Emam OA, El Domiaty HA, Elbohy WR, Aboelenin HM, Salem NA. Genetic association between interleukin-10 gene rs1518111 and rs3021094 polymorphisms and risk of type 1 diabetes and diabetic nephropathy in Egyptian children and adolescents. Pediatr Diabetes. 2021;22(4):567–76.
Završnik M, Letonja J, Makuc J, Šeruga M, Cilenšek I, Petrovič D. Interleukin-4 (IL4) -590 C/T (rs2243250) gene polymorphism is not associated with diabetic nephropathy (DN) in Caucasians with type 2 diabetes mellitus (T2DM). Bosn J Basic Med Sci. 2018;18(4):347–51.
Yazdi AS, Ghoreschi K. The Interleukin-1 family. Adv Exp Med Biol. 2016;941:21–9.
Lin N, Lu H, Cheng X, Zhao Y, Wan Q, Luo Y, Miao Y, Bai X, Liu D, Wang C. Association between the interleukin-1B polymorphism at rs16944 T > C and diabetic retinopathy. Int J Immunogenet. 2023;50(1):34–40.
Buraczynska M, Ksiazek K, Wacinski P, Zaluska W. Interleukin-1β gene (IL1B) polymorphism and risk of developing diabetic nephropathy. Immunol Invest. 2019;48(6):577–84.
Stefanidis I, Kreuer K, Dardiotis E, Arampatzis S, Eleftheriadis T, Hadjigeorgiou GM, Zintzaras E, Mertens PR. Association between the interleukin-1β gene (IL1B) C-511T polymorphism and the risk of diabetic nephropathy in type 2 diabetes: a candidate-gene association study. DNA Cell Biol. 2014;33(7):463–8.
Yuan Y, Li L, Wang X, Zhang P, Wang J, Xiao Y. Correlation between plasma NLRP3, IL-1β, and IL-18 and diabetic nephropathy in patients with type 2 diabetes. Altern Ther Health Med. 2023;29(4):52–6.
Han Q, Wang S, Zhang J, Zhang R, Guo R, Wang Y, Li H, Xu H, Liu F. The association between cigarette smoking and diabetic nephropathy in Chinese male patients. Acta Diabetol. 2018;55(11):1131–41.
Jaimes EA, Zhou MS, Siddiqui M, Rezonzew G, Tian R, Seshan SV, Muwonge AN, Wong NJ, Azeloglu EU, Fornoni A, Merscher S, Raij L. Nicotine, smoking, podocytes, and diabetic nephropathy. Am J Physiol Ren Physiol. 2021;320(3):F442–53.
Ma L, Jiang Y, Kong X, Liu Q, Zhao H, Zhao T, Cao Y, Li P. Interaction of MTHFR C677T polymorphism with smoking in susceptibility to diabetic nephropathy in Chinese men with type 2 diabetes. J Hum Genet. 2019;64(1):23–8.
Xue P, Cao H, Ma Z, Zhou Y, Wang N. Transcription factor 7-like 2 gene- smoking interaction on the risk of diabetic nephropathy in Chinese Han population. Genes Environ. 2021;43(1):26.
Benhalima K, Van Crombrugge P, Moyson C, Verhaeghe J, Vandeginste S, Verlaenen H, Vercammen C, Maes T, Dufraimont E, De Block C, Jacquemyn Y, Mekahli F, De Clippel K, Van Den Bruel A, Loccufier A, Laenen A, Minschart C, Devlieger R, Mathieu C. Risk factor screening for gestational diabetes mellitus based on the 2013 WHO criteria. Eur J Endocrinol. 2019;180(6):353–63.
Gross JL, de Azevedo MJ, Silveiro SP, Canani LH, Caramori ML, Zelmanovitz T. Diabetic nephropathy: diagnosis, prevention, and treatment. Diabetes Care. 2005;28(1):164–76.
American Diabetes Association. 11. Microvascular complications and foot care: standards of medical care in diabetes-2020. Diabetes Care. 2020;43(Suppl 1):S135–S151.
Kim Y, Park T. Robust Gene-Gene interaction analysis in genome wide association studies. PLoS ONE. 2015;10(8):e0135016.
Rogus J, Beck JD, Offenbacher S, Huttner K, Iacoviello L, Latella MC, de Gaetano M, Wang HY, Kornman KS, Duff GW. IL1B gene promoter haplotype pairs predict clinical levels of interleukin-1beta and C-reactive protein. Hum Genet. 2008;123(4):387–98.
Hall SK, Perregaux DG, Gabel CA, Woodworth T, Durham LK, Huizinga TW, Breedveld FC, Seymour AB. Correlation of polymorphic variation in the promoter region of the interleukin-1 beta gene with secretion of interleukin-1 beta protein. Arthritis Rheum. 2004;50:1976–83.
D’Amico AG, Maugeri G, Rasà DM, Bucolo C, Saccone S, Federico C, Cavallaro S. D’Agata V. Modulation of IL-1β and VEGF expression in rat diabetic retinopathy after PACAP administration. Peptides. 2017;97:64–9.
Bi X, Ai H, Wu Q, Fan Q, Ding F, Hu C, Ding W. Insulin resistance is associated with Interleukin 1β (IL-1β) in Non-Diabetic Hemodialysis patients. Med Sci Monit. 2018;24:897–902.
Lee SH, Ihm CG, Sohn SD, Lee TW, Kim MJ, Koh G, Oh SJ, Woo JT, Kim SW, Kim JW, Kim YS, Lee BC, Kim SD, Cho BS, Lee HJ, Chung JH. Polymorphisms in interleukin-1 beta and Interleukin-1 receptor antagonist genes are associated with kidney failure in Korean patients with type 2 diabetes mellitus. Am J Nephrol. 2004;24(4):410–4.
Tarnow L, Pociot F, Hansen PM, Rossing P, Nielsen FS, Hansen BV, Parving HH. Polymorphisms in the interleukin-1 gene cluster do not contribute to the genetic susceptibility of diabetic nephropathy in Caucasian patients with IDDM. Diabetes. 1997;46(6):1075–6.
Loughrey BV, Maxwell AP, Fogarty DG, Middleton D, Harron JC, Patterson CC, Darke C, Savage DA. An interluekin 1B allele, which correlates with a high secretor phenotype, is associated with diabetic nephropathy. Cytokine. 1998;10(12):984–8.
Hameed I, Masoodi SR, Malik PA, Mir SA, Ghazanfar K, Ganai BA. Genetic variations in key inflammatory cytokines exacerbates the risk of diabetic nephropathy by influencing the gene expression. Gene. 2018;661:51–9.
Wu J, Jiang C, Hua Y, Liu X, You C. Association between polymorphisms of cytokine genes and diabetic nephropathy: A comprehensive systematic review and meta-analysis. Int J Clin Pract. 2021;75(11):e14634.
Chakkarwar V. Smoking in diabetic nephropathy: sparks in the fuel tank?? World J Diabetes. 2012;3(12):186.
Acknowledgements
None.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
All authors have read and approved the manuscript; TX performed data curation, formal analysis, methodology, writing and editing the manuscript. ZT: performed methodology and review the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Each participant understood the process of the study and signed a written informed consent before the start of the study. All study protocols of the current study were approved by the ethics committee of affiliated Hospital of Nantong University. All methods were performed in accordance with the Declaration of Helsinki.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Xie, T., Tang, Z. Association of the interaction between interleukin-1β gene polymorphism and smoking status with the diabetic nephropathy risk in a Chinese Han population. Diabetol Metab Syndr 17, 101 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13098-025-01667-y
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13098-025-01667-y