Interaction between ELMO1 gene polymorphisms and environment factors on susceptibility to diabetic nephropathy in Chinese Han population
Diabetology & Metabolic Syndrome volume 11, Article number: 97 (2019)
The association of diabetic nephropathy (DN) risk with single nucleotide polymorphisms (SNPs) within Engulfment and Cell Motility 1 (ELMO1) gene and gene–environment synergistic effect have not been extensively examined in, therefore, the purpose of this study is to explore the association between multiple SNPs in ELMO1 gene, and the relationship between gene–environment synergy effect and the risk of DN.
Genotyping for 4 SNPs was performed with polymerase chain reaction (PCR) and following restriction fragment length polymorphism (RFLP) methods. Hardy–Weinberg balance of the control group was tested by SNPstats (online software: http://bioinfo.iconologia.net/snpstats). The best combination of four SNPs of ELMO1 gene and environmental factors was screened by GMDR model. Logistic regression was used to calculating the OR values between different genotypes of ELMO1 gene and DN.
The rs741301-G allele and the rs10255208-GG genotype were associated with an increased risk of DN risk, adjusted ORs (95% CI) were 1.75 (1.19–2.28) and 1.41 (1.06–1.92), respectively, both p-values were < 0.001. We also found that the others SNPs-rs1345365 and rs7782979 were not significantly associated with susceptibility to DN. GMDR model found a significant gene–alcohol drinking interaction combination (p = 0.0107), but no significant gene–hypertension interaction combinations. Alcohol drinkers with rs741301-AG/GG genotype also have the highest DN risk, compared to never drinkers with rs741301-AA genotype, OR (95% CI) 3.52 (1.93–4.98).
The rs741301-G allele and the rs10255208-GG genotype, gene–environment interaction between rs741301 and alcohol drinking were all associated with increased DN risk.
Diabetic nephropathy (DN) is one of the most destructive microvascular complications of all the complications of diabetes, and it is the main cause of end-stage renal disease (ESRD), which accounts for a large proportion of the morbidity and mortality of diabetes . DN is also an important risk factor for the development of cardiovascular disease and chronic kidney disease  and one of the main causes of death . The pathogenesis of DN was multifactorial and was also not been well known. Previous studies have reported that some genetic and environmental risk factors may lead to its progress, including hyperglycemia, hypertension and hyperlipidemia, which were common modifiable risk factors. In addition, some unmodifiable risk factors should also draw public attentions, such as genetic factors [4, 5].
Genome-Wide Association Studies (GWAS) and meta-analysis [6, 7] have reported some genes related with susceptibility to DN in different populations. In 2005, a Japanese GWAS study  confirmed that Engulfment and Cell Motility 1 (ELMO1) gene is a new candidate gene. In this study, they evaluated more than 80,000 loci and found that one SNP locus in intron 18 of ELMO1 gene was closely related to DN. Then, some studies also reported that this gene was associated with kidney disease attributed to T2DM in different populations [9,10,11], including Chinese populations . The detailed mechanism for association between ELMO1 gene polymorphisms and DN susceptibility was not well known. However, a previous study  showed that mutations in elmo1 gene can lead to the disorder of extracellular matrix (ECM) metabolism, which leads to the accumulation of ECM, thickening of renal tubules and glomerular basement membrane, thus increasing the risk of DN. However, to date, the association between DN risk and SNPs within ELMO1 gene have not been extensively examined in Chinese populations. Strong evidences [13, 14] have suggested the important role for the combination of environmental exposures and genetic factors for development or progression of the DN. However, till now, no study investigated the ELMO1 gene–environment interaction. Therefore, the present study aimed to evaluate the impact of four SNPs in the ELMO1 gene, and its interaction with environmental risk factors on susceptibility to DN.
Materials and methods
There were a total of 1325 participants including 660 T2DM patients with DN and 665 T2DM patients without DN. Those participants with a fasting glucose ≥ 126 mg/dl (7.0 mmol/l), or a 2 h postprandial blood glucose ≥ 200 mg/dl (11.0 mmol/l), or if hypoglycemic therapy (oral agents or insulin) had been started were considered as T2DM patients . The diagnostic criteria of DN was made according to World Health Organization 1999 . Those T2DM patients were considered as cases with DN, including those with persistent urine albuminuria more than 300 mg/l in consecutive twice measurements, or without renal failure: serum creatinine more than 1 mg/dl. Those T2DM patients without DN were included in the control group. All subjects in the study received detailed clinical and biochemical examinations. Those T2DM patients with end-stage renal disease (ESRD) were excluded from both case and control group. All the subjects were Han people, and there was no genetic and blood relationship between any two subjects. Before participating in the study, all subjects understood the content of the study and signed the informed consent. All participants will fill in a self-designed questionnaire, including general demographic information, lifestyle risk factors such as smoking and drinking, physical indicators measurement, the history of major diseases and other data. Smokers are those who report smoking at least once a day for a year or more. Drinking is the sum of milliliters of alcohol extracted weekly from wine, beer and spirits. Hypertension patients were those whose SBP was equal or more than 140 mmHg and/or DBP was equal or more than 90 mmHg and/or use of antihypertensive medication.
Genomic DNA extraction and genotyping
We selected the SNPs according to two methods as following: firstly, we selected SNPs from ELMO1 gene family; secondly, we selected the SNPs within ELMO1 gene, which have been reported in previous studies, but no consistent results were obtained on the relation between SNPs and DN risk. According to the standard instructions, we will take 3 ml blood samples from all participants, which are processed by EDTA and stored in a refrigerator at – 20 °C for DNA genomic DNA extraction. Genotyping was performed with a polymerase chain reaction (PCR) and following restriction fragment length polymorphism (RFLP). The description, primers and enzyme for the four SNPs were shown in Additional file 1: Table S1, including rs741301, rs1345365, rs10255208 and rs7782979. A 25 μl reaction mixture including 1.25 μl SNP Genotyping Assays (20×), 12.5 μl Genotyping Master Mix (2×), 20 ng DNA. PCR cycling conditions consisted of an initial denaturation at 96 °C for 7 min, followed by 35 cycles of 96 °C for 20 s, 57 °C for 30 s, and 72 °C for 45 s, ending with a final elongation step at 72 °C for 5 min. For quality control, genotyping was performed with blinding to the grouping of participants, and both controls and cases were randomly selected 10%, which were genotyped for two times by different staffs, and the reproducibility should be 100%.
In current study, the Hardy–Weinberg equilibrium (HWE) test and comparison of distribution of alleles and genotypes between case and control groups were performed with Chi-squared test. The mean ± standard deviation (SDS) was used to represent the continuous variables of normal distribution, and Student t test was used to compare the case group and the control group. Logistic regression model was used to calculate the statistical relationship between four SNPs and DN risk. The best interaction combinations associated with DN among the 4 SNPs in ELMO1 gene was determined with generalized multifactor dimensionality reduction (GMDR), a sign or permutation test (providing empirical p values) used for predicting accuracy was employed to measure the significance of an identified model. Statistically significance was determined when p-values was less than 0.05.
A description for DN patients and normal controls regarding demographic and general or clinical characteristics is shown in Table 1. A total of 1325 participants including 660 T2DM patients with DN and 665 T2DM patients without DN (controls). The average age for all participants was 66.0 ± 12.8 years. There was no significant difference found in several parameters, including gender, smoking rates, age, duration of diabetes, FPG and BMI between the two groups (all p-values were more than 0.05). In contrast, the percentages of participants who consumed alcohol, hypertension patients, means of HbA1c, creatinine, urea and urine albumin/creatinine ratios (ACR) were higher in the DN patients than that in controls.
The genotype frequencies in the control of the current study were all distributed accordingly to HWE (Table 2). The frequency of the rs741301-G allele was 30.7% in DN patients and 20.2% in normal controls, in addition, the frequency of the rs10255208-G allele was 28.3% in DN patients and 20.4% in normal controls, which was also indicating a significantly statistical difference. The rs741301-G allele and the rs10255208-GG genotype were associated with an increased risk of DN risk, adjusted ORs (95% CI) were 1.75 (1.19–2.28) and 1.41 (1.06–1.92), respectively. We also found that the others SNPs-rs1345365 and rs7782979 were not correlated with DN susceptibility significantly.
The gene–hypertension or alcohol drinking interaction test were determined by GMDR model (Table 3). We found a significant gene–alcohol drinking interaction combination, but no significant gene–hypertension interaction combinations, after adjusting for age, gender, BMI, smoking and alcohol drinking covariates. A two-locus including rs741301 and alcohol drinking was significantly in the GMDR model. In order to obtain the odds ratios and 95% CI for the joint effects of gene–alcohol drinking on DN, we conducted stratified analysis for interaction effect using logistic regression. We found that alcohol drinkers with rs741301-AG/GG genotype also have the highest DN risk, compared to never drinkers with rs741301-AA genotype, OR (95% CI) 3.52 (1.93–4.98) (Fig. 1).
In the current study, we found that rs741301-G allele and the rs10255208-GG genotype were associated with an increased risk of DN risk respectively. We also found that the others SNPs-rs1345365 and rs7782979 were not significantly associated with susceptibility to DN. DN is the most common cause of ESRD and the most common microvascular complication of diabetes mellitus . However, its detailed pathogenesis was not well known yet, and both environmental and genetic factors were reported involving in the occurrence and development of DN susceptibility [18, 19]. ELMO1 is a protein composed of 720 amino acids, encoded by elmo1 gene and located on chromosome 7p14.2-14.1. It may be a new and powerful candidate gene of DN, which has influence on the movement and phagocytosis of apoptotic cells . Animal studies in mice  suggest that ELMO1 protein plays an important role in the pathogenesis of DN and proteinuria. Over the last decade, some studies were performed to test the relationship between ELMO1 gene SNPs and DN risk, but these studies concluded inconsistent results, and few studies were performed in Chinese populations. Mehrabzadeh et al.  investigated the relationship between ELMO1-rs741301 and the risk of DN in Iranian population, and considered that ELMO1-rs741301 is an important candidate gene for DN susceptibility. Bodhini et al.  reported a significant association of the ELMO1 rs741301 SNPs with DN in south Indians. Wu et al.  also verified the association of rs741301 with DN, and rs1345365 was not associated with DN risk, however the sample size for this study was very small, just 200 unrelated Chinese subjects were enrolled. Some studies including GWAS and replication studies were performed in Japanese type 2 diabetic patients , replication studies in American Indian study  and GoKinD collection study , African American cohorts . The key role of ELMO1 as a DN susceptibility gene has been verified. Although many studies have confirmed that elmo1 is a candidate gene for DN, the risk locus and risk allele are not consistent in all populations . Different perspectives also existed in some studies. Yahya et al.  suggested that ELMO1-rs74130 was not associated with susceptibility to DN. Hanson et al.  indicated that ELMO1-rs1345365, but not rs741301 associated with DN risk. In terms of rs10255208, which was less studies previously, but this SNP has been reported association with T2DM in Tunisian Arabs .
The pathogenesis and process of DN are very complex. The factors of DN susceptibility mainly include genetic factors, environmental factors and the synergistic effect of genetic and environmental factors [28, 29]. Some studies have showed that the interaction between gene and environment factors were associated with the DN risk. This study is the first to confirm the synergistic effect of elmo1 gene polymorphism and environmental factors on DN susceptibility. We found a significant gene–alcohol drinking interaction combination, but no significant gene–hypertension interaction combinations. Alcohol drinkers with rs741301-AG/GG genotype also have the highest DN risk, compared to never drinkers with rs741301-AA genotype. The detailed mechanism for ELMO1 gene-DN association was not well known. However, a previous study  showed that mutations in elmo1 gene can lead to the disorder of ECM metabolism, the accumulation of which could thicken of renal tubules and glomerular basement membrane, thus increasing the risk of DN. Another study  suggested that ELMO1 also could plays an important role in the development of DN by increasing oxidative stress (OS), the level of elmo1 gene expression in diabetic mice was consistent with the degree of renal fibrosis and urinary albumin excretion. There was also a positive correlation between ROS and ELMO1 expression.
The limitations of this study were: Firstly, more environment factors should be included in the interaction analysis; secondly, the study population were all Chinese Han, which may not present on behalf of the Chinese populations, because there were 56 races in China. Secondly, just four SNPs were selected for genotyping, more SNPs should be included, and gene–gene interaction should be investigated in the future.
In conclusion, we found that the rs741301-G allele and the rs10255208-GG genotype, gene–environment interaction between rs741301 and alcohol drinking were all associated with increased DN risk. This study not only investigated whether the ELMO1 gene was related to DN susceptibility, but also investigated whether the interaction effect existed between this gene and some significant environmental factors, such as hypertension, alcohol drinking. Previous studies merely investigated this interaction effect between ELMO1 gene and alcohol drinking, especially in Chinese Han populations.
Availability of data and materials
Zelmanovitz TGF, Balthazar AS, Thomazelli FCS, Jorge D, Matos JD, Canani LH. Diabetic nephropathy. Diabetol Metab Syndr. 2009;1:1–17.
Gnudi L, Coward RJM, Long DA. Diabetic nephropathy: perspective on novel molecular mechanisms. Trends Endocrinol Metab. 2016;27(11):820–30.
Kos I, Prkacin I. Diabetic nephropathy as a cause of chronic kidney disease. Acta Med Croatica. 2014;68:375–81.
Jenum AK, Holme I, Graf-Iversen S, Birkeland KI. Ethnicity and sex are strong determinants of diabetes in an urban Western society: implications for prevention. Diabetologia. 2005;48:435–9.
Temelkova-Kurktschiev T, Stefanov T. Lifestyle and genetics in obesity and type 2 diabetes. Exp Clin Endocrinol Diabetes. 2012;120:1–6.
McKnight AJ, Duffy S, Maxwell AP. Genetics of diabetic nephropathy: a long road of discovery. Curr Diab Rep. 2015;15(7):1–11.
Mooyaart AL, Valk EJ, van Es LA, Bruijn JA, de Heer E, Freedman BI, Dekkers OM, Baelde HJ. Genetic associations in diabetic nephropathy: a meta-analysis. Diabetologia. 2011;54:544–53.
Shimazaki A, Kawamura Y, Kanazawa A, Sekine A, Saito S, Tsunoda T, Koya D, Babazono T, Tanaka Y, Matsuda M, Kawai K, Iiizumi T, Imanishi M, Shinosaki T, Yanagimoto T, Ikeda M, Omachi S, Kashiwagi A, Kaku K, Iwamoto Y, Kawamori R, Kikkawa R, Nakajima M, Nakamura Y, Maeda S. Genetic variations in the gene encoding ELMO1 are associated with susceptibility to diabetic nephropathy. Diabetes. 2005;54(4):1171–8.
Pezzolesi MG, Katavetin P, Kure M, Poznik GD, Skupien J, Mychaleckyj JC, Rich SS, Warram JH, Krolewski AS. Confirmation of genetic associations at ELMO1 in the GoKinD collection supports its role as a susceptibility gene in diabetic nephropathy. Diabetes. 2009;58(11):2698–702.
Hanson RL, Millis MP, Young NJ, Kobes S, Nelson RG, Knowler WC, DiStefano JK. ELMO1 variants and susceptibility to diabetic nephropathy in American Indians. Mol Genet Metab. 2010;101(4):383–90.
Wu HY, Wang Y, Chen M, Zhang X, Wang D, Pan Y, Li L, Liu D, Dai XM. Association of ELMO1 gene polymorphisms with diabetic nephropathy in Chinese population. J Endocrinol Investig. 2013;36(5):298–302.
Shimazaki A, Tanaka Y, Shinosaki T, Ikeda M, Watada H, Hirose T, Kawamori R, Maeda S. ELMO1 increases expression of extracellular matrix proteins and inhibits cell adhesion to ECMs. Kidney Int. 2006;70:1769–76.
Pettitt DJ, Saad MF, Bennett PH, Nelson RG, Knowler WC. Familial predisposition to renal disease in two generations of Pima Indians with type 2 (non-insulin-dependent) diabetes mellitus. Diabetologia. 1990;33(7):438–43.
Quinn M, Angelico MC, Warram JH, Krolewski AS. Familial factors determine the development of diabetic nephropathy in patients with IDDM. Diabetologia. 1996;39(8):940–5.
American Diabetes Association. Diagnosis and classification of diabetes mellitus. Diabetes Care. 2012;35(Suppl. 1):S64–71.
World Health Organization. Part 1: diagnosis and classification of diabetes mellitus: report of a WHO consultation. In: Alwan A, King H, editors. Definition, diagnosis and classification of diabetes mellitus and its complications. Geneva: Word Health Department of Noncommunicable Disease Surveillance; 1999.
Betz B, Conway BR. Recent advances in animal models of diabetic nephropathy. Nephron Exp Nephrol. 2014;126(4):191–5.
Sidaway P. Diabetes: epigenetic changes lead to impaired wound healing in patients with T2DM. Nat Rev Endocrinol. 2015;11(2):65.
Keating ST, El-Osta A. Epigenetic changes in diabetes. Clin Genet. 2013;84(1):1–10.
van Ham TJ, Kokel D, Peterson RT. Apoptotic cells are cleared by directional migration and elmo1-dependent macrophage engulfment. Curr Biol. 2012;22(9):830–6.
Hathaway CK, Chang AS, Grant R, Kim HS, Madden VJ, Bagnell CR Jr, Jennette JC, Smithies O, Kakoki M. High Elmo1 expression aggravates and low Elmo1expression prevents diabetic nephropathy. Proc Natl Acad Sci USA. 2016;113(8):2218–22.
Mehrabzadeh M, Pasalar P, Karimi M, Abdollahi M, Daneshpour M, Asadolahpour E, Razi F. Association between ELMO1 gene polymorphisms and diabetic nephropathy in an Iranian population. J Diabetes Metab Disord. 2016;15:43.
Bodhini D, Chidambaram M, Liju S, Revathi B, Laasya D, Sathish N, Kanthimathi S, Ghosh S, Anjana RM, Mohan V, Radha V. Association of rs11643718 SLC12A3 and rs741301 ELMO1 variants with diabetic nephropathy in south Indian population. Ann Hum Genet. 2016;80(6):336–41.
Leak T, Perlegas P, Smith S, Keene KL, Hicks P, Langefeld C, et al. Variants in intron 13 of the ELMO1 gene are associated with diabetic nephropathy in African Americans. Ann Hum Genet. 2009;73(2):152–9.
Chang Y-C, Chang EY-C, Chuang L-M. Recent progress in the genetics of diabetic microvascular complications. World J Diabetes. 2015;6(5):715.
Yahya MJ, Ismail PB, Nordin NB, Akim ABM, Yusuf WSBM, Adam NLB, Yusoff MJ. Association of CCL2, CCR5, ELMO1, and IL8 polymorphism with diabetic nephropathy in Malaysian type 2 diabetic patients. Int J Chronic Dis. 2019;2019:2053015.
Turki A, Mzoughi S, Mtitaoui N, Khairallah M, Marmouch H, Hammami S, Mahjoub T, Almawi WY. Gender differences in the association of ELMO1 genetic variants with type 2 diabetes in Tunisian Arabs. J Endocrinol Investig. 2018;41(3):285–91.
Cai W, Li J, Xu JX, Liu Y, Zhang W, Xiao JR, Zhu LY, Liu JY. Association of 2184AG polymorphism in the RAGE gene with diabetic nephropathy in chinese patients with type 2 diabetes. J Diabetes Res. 2015;2015:310237.
Zhang Y, Jia N, Hu F, Fan N, Guo X, Du H, Mei C, Gao C. Association of single-nucleotide polymorphisms in the RAGE gene and its gene–environment interactions with diabetic nephropathy in Chinese patients with type 2 diabetes. Oncotarget. 2017;8(57):96885–92.
We appreciate the cooperation of the families and individuals who cooperated in this study.
Ethics approval and consent to participate
We have obtained consent to publish from the participant to report individual patient data.
Consent for publication
We have obtained consent to publish from the participant.
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
Hou, Y., Gao, Y., Zhang, Y. et al. Interaction between ELMO1 gene polymorphisms and environment factors on susceptibility to diabetic nephropathy in Chinese Han population. Diabetol Metab Syndr 11, 97 (2019). https://doi.org/10.1186/s13098-019-0492-0
- Diabetic nephropathy
- Single nucleotide polymorphisms
- Engulfment and Cell Motility 1