Association of Tagging SNPs in the MTHFR Gene with Risk of Type 2 Diabetes Mellitus and Serum Homocysteine Levels in a Chinese Population
Diabetes is a global public health crisis, and the prevalence is increasing rapidly. Folate supplementation is proved to be effective in reducing the risk of diabetes or improving its symptoms. Methylenetetrahydrofolate reductase is an important enzyme involved in folate metabolism. The aim of this study is to examine whether polymorphisms in the MTHFR gene are associated with risk of type 2 diabetes mellitus (T2DM) and fasting total serum homocysteine (tHcy) levels. We genotyped nine tagging SNPs in the MTHFR gene in a case-control study, including 595 T2DM cases and 681 healthy controls in China. We found that C allele of rs9651118 had significant decreased risk of T2DM (adjusted odds ratio (OR) = 0.69, 95% confidence interval (CI): 0.55–0.87, ) compared with T allele. Haplotype analysis also showed that MTHFR CTCCGA haplotype (rs12121543-rs13306553-rs9651118-rs1801133-rs2274976-rs1801131) had significant reduced risk of T2DM (adjusted OR = 0.71, 95% CI: 0.58–0.87, ) compared with CTTTGA haplotype. Besides, the MTHFR rs1801133 was significantly associated with serum levels of tHcy in healthy controls (). These associations were still significant after Bonferroni corrections (). These findings suggest that variants in the MTHFR gene may influence the risk of T2DM and tHcy levels.
Diabetes mellitus (DM) is one of the primary chronic diseases that causes premature death and disability. It was estimated to be the third most challenging disease that is threatening public health just after malignant tumors and cardiocerebral vascular diseases in the world . In recent years, the prevalence of DM was rising rapidly worldwide, particularly in developing countries . In China, age-standardized incidence rate of DM was 2.4% in 1994 but raised to 9.7% in 2007-2008 [3, 4]. The importance of conducting prevention and early detection of DM has been increasingly emphasized.
Type 2 DM (T2DM) is the most common form of diabetes, which accounts for about 90% of all diagnosed DM cases. It is a multifactorial disease that is due to a combination of environmental and genetic risk factors, including dietary habit, obesity, physical inactivity, smoking, alcohol consumption, and genetic variants . Studies have confirmed the relationship between homocysteinemia and diabetic complications, such as diabetic nephropathy, cardiovascular disease, and hypertension [6–8]. It is also reported that a high concentration of homocysteine is associated with the development of T2DM [9, 10]. MTHFR (methylenetetrahydrofolate reductase), a folate-dependent enzyme, plays an important role in the conversion of homocysteine to methionine by converting 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate. Studies have found that specific genetic polymorphisms in the MTHFR gene could lead to change of MTHFR enzyme activity [11, 12]. Previous epidemiologic studies have also shown that genetic polymorphisms in MTHFR gene may be related to cancers, ischemic stroke, hypertension, and so forth [13–15]. But the association of MTHFR genotypes and susceptibility to T2DM in Chinese population has not been fully studied.
To help clarify whether the MTHFR variants are associated with susceptibility of T2DM and serum total homocysteine (tHcy) level, we examined nine tagging SNPs in the MTHFR gene (rs12121543, rs13306561, rs13306553, rs9651118, rs1801133, rs2274976, rs4846048, rs1801131, and rs17037396) in a case-control study in China.
2. Material and Methods
2.1. Study Participants
In this study, we consecutively enrolled 595 T2DM patients in the first affiliated hospital of Jilin university, Changchun, China, aged from 35 to 80, between June 2011 and June 2013. A case of T2DM was considered to be confirmed if the participant’s reported glucose level met the well-established diagnostic criteria recommended by American diabetes association: with either the fasting plasma glucose (FPG) of ≥7.0 mmol/L or the 2-h value in the 75-g oral glucose tolerance test (OGTT) ≥11.1 mmol/L . A total of 655 healthy control subjects were selected during the same period and from the same hospital and were frequency matched to the cases by age (5-year age groups) and gender. All controls were individuals free of diabetes that are determined by medical history or clinical examinations. At enrollment, demographic and medical histories, including age, gender, body mass index, smoking habits, and hypertension status were collected from each subject by a trained interviewer using a structured questionnaire. Blood samples were collected after a 12-hour overnight fast and then separated into serum, red blood cells, and buffy coat. Written informed consent was obtained from all enrolled participants and this study was approved by the Ethics Committee of Jilin University.
2.2. SNP Selection
Tag SNPs were selected by searching Han Chinese data from the HapMap project (http://www.hapmap.org/) using the Tagger program. The following criteria were used to identify tagging SNPs: (a) SNPs located in the gene or within the 2-kb region flanking the gene, (b) a minor allele frequency ≥0.1, and (c) other unselected SNPs could be captured by one of the tagging SNPs with a linkage disequilibrium of . As a result, a total of 9 tag SNPs were identified.
2.3. Laboratory Tests
Serum levels of total cholesterol (TC), low density lipoprotein-cholesterol (LDL-C), and high density lipoprotein-cholesterol (HDL-C) were determined in cases and controls by commercial kits from the Nanjing Jiancheng Bio-company (Nanjing, China). Serum tHcy levels in healthy controls were measured by cycling enzymatic method (Nanjing Jiancheng Bio-company, Nanjing, China) and were analyzed with an automatic biochemical analyzer (AU 2700 Olympus, First Chemical Ltd, Tokyo, Japan).
Genomic DNA was extracted from peripheral blood leucocytes using DNA Extraction Kit (Qiagen, Hilden, Germany). Genotyping was performed on the ABI PRISM 7900HT Sequence Detection System (Applied Biosystems, Foster City, CA, USA), using the TaqMan assay. The genotyping call rate was >95%, and the completion rate was >99%. The quality and potential misclassification of the genotyping were assessed by regenotyping 10% of duplicate DNA samples that were randomly selected from the whole population and placed within the same reaction plates used for the study subjects. The concordance rate for the quality control samples was 100%.
2.4. Statistical Analysis
We used SAS software (version 9.3; SAS Institute, Inc.) for the statistical analyses. statistics and the -test were used to evaluate case-control differences in the distribution of risk factors. Variables were tested for normality with Shapiro-Wilk statistics. Skewed data, including age, BMI, TC, LDL-C, HDL-C, and tHcy, were log transformed and expressed as medians and interquartile ranges. A two tailed value of 0.05 was considered statistically significant. The odds ratios (ORs) and 95% confidence intervals (CIs) for the associations between the SNPs and T2DM risk were estimated by unconditional logistic regression. Hardy-Weinberg equilibrium for genotypic distribution and linkage disequilibrium between loci were assessed by HaploView version 4.0 (Daly Lab at the Broad Institute, Cambridge, MA, USA) . Associations between haplotypes (>1% frequency) and the risk of T2DM were evaluated by computing OR and 95% CI using HAPSTAT, assuming an additive model, using the most common haplotype as the referent category . Both univariate ANOVA and multivariate ANCOVA analyses adjusting for age, sex, BMI, and hypertension were performed to determine the effects of the MTHFR polymorphisms on serum tHcy levels in healthy controls. values (from the additive model or dominant model) of less than 0.0056 were deemed to be statistically significant after Bonferroni corrections.
Characteristics of the study subjects are shown in Table 1. Cases and controls were evenly matched by age and gender. Cases were more probably to have higher BMI (24.6 kg/m2 versus 23.1 kg/m2) and hypertension (43.9% versus 24.5%). Besides, cases have significant lower levels of serum HDL-C and higher levels of serum total cholesterol and LDL-C than that in controls.
The associations of MTHFR variants and risk of T2DM are presented in Table 2. The genotype distributions of these nine polymorphisms showed no deviation from the expected Hardy-Weinberg equilibrium among controls (). Of these SNPs, C allele of rs9651118 conferred significant lowered risk of T2DM (crude OR = 0.71, 95% CI: 0.57–0.88, ) compared with T allele, and rs1801131 AC genotype conferred significant increased risk of T2DM compared with AA genotype (crude OR = 1.29, 95% CI: 1.02–1.64, ). These associations remained significant after adjustment for other risk factors (age, sex, BMI, and hypertension). But only rs9651118 was still significant after Bonferroni correction in dominant model (adjusted OR = 0.69, 95% CI: 0.55–0.87, ). None of the other SNPs examined was associated with the risk of T2DM.
Six SNPs in the MTHFR gene (rs12121543, rs13306553, rs9651118, rs1801133, rs2274976, and rs1801131) were in linkage disequilibrium with ranging from 0.79 to 0.99 and ranging from 0.04 to 0.86. Subjects carrying the MTHFR CTCCGA haplotype had a significant reduced risk of T2DM (OR = 0.73, 95% CI: 0.60–0.88, ) compared with those carrying the CTTTGA haplotype (Table 3). This association remained significant after adjustment for other risk factors and after Bonferroni correction.
Finally, we investigated the associations between the MTHFR SNPs and serum tHcy levels in the control population. Carriers of the mutant alleles of rs1801133 were significantly associated with increased serum level of tHcy [C/T: 12.8(10.9–14.4) μmol/L; T/T: 11.9(9.5–16.6) μmol/L] compared with carriers of the CC genotype [11.3(9.3–13.7) μmol/L]. Rs9651118 T/C and C/C genotypes [12.1(10.1–14.1) μmol/L; 11.3(9.3–13.4) μmol/L] conferred significant decreased serum level of tHcy in controls compared with TT genotype [12.5(9.8–14.6) μmol/L]. These associations were also significant in multivariate ANCOVA analyses. But after Bonferroni correction only rs1801133 was significantly associated with serum level of tHcy. None of the other studied SNPs was associated with serum tHcy level (Table 4).
In this molecular epidemiologic study, polymorphisms in the MTHFR gene were fully studied for their association with susceptibility to T2DM and serum tHcy levels. We demonstrated that one genetic mutation (rs9651118) in MTHFR gene was significantly associated with risk of T2DM, and rs1801133 was significantly associated with serum tHcy levels in healthy population.
MTHFR gene is located on chromosome 1p36.3 in humans. To date, over 40 point genetic mutations in MTHFR gene have been identified, of which A1298C (rs1801131) and C677T (rs1801133) showed the most clinical significance. The MTHFR rs1801133 is a missense mutation that results in approximately 70% and 35% reduction of normal MTHFR enzyme activity in TT and CT genotype carriers, respectively . MTHFR activity has an adverse effect on serum tHcy level. MTHFR rs1801133 TT genotype was also associated with elevated tHcy levels, predominantly in individuals who have a low plasma folate level (<15.4 nmol/L), compared to the normal genotype . In our study, although T allele of rs1801133 has nonsignificant raised risk of T2DM (OR = 1.22, 95%CI: 0.97–1.55), we found that the mutant alleles were associated with significantly higher level of tHcy in healthy controls. Besides, we observed that the average tHcy level in TT genotype carriers was 13.2 μmol/L which was higher than that in CT genotype carriers (13.0 μmol/L). These results were in accordance with previous studies.
A1298C (rs1801131) is a common mutation in the MTHFR gene that results in the conversion of adenine to cytosine, which results in reduction of MTHFR enzyme activity . Few studies have evaluated the relationship between the MTHFR A1298C polymorphism and susceptibility of diabetes, and the results remain inclusive. Previous study suggested that MTHFR A1298C polymorphism is a risk factor for T2DM in Egyptian patients, while other failed to find any association in Taiwanese or Moroccan [21–23]. A recent Meta-analysis found that MTHFR A1298C has a significant association with diabetes in Asian population under dominant model but not in Caucasians . The association between A1298C and serum homocysteine level was also controversial. Kumar et al. found that subjects with MTHFR rs1801131 CC genotype had significant higher level of homocysteine (16.3 μmol/L) compared with AA genotype (14.4 μmol/L) in Indian population, but this association was not found in Scotch [25, 26]. In our study, although we found AC genotype had relatively higher risk of T2DM, the significance no longer exists after Bonferroni correction. In addition, we did not find its relation with homocysteine level.
Rs9651118 is located in the intron region in MTHFR gene with unknown genetic function. It was once reported that C allele of rs9651118 was associated with reduced lung cancer risk in never smokers, and a nominally significant association with schizophrenia in the form of haplotypes (rs1801133, rs17421511, rs17037396, and rs9651118) in Japanese population [27, 28]. Rs9651118 TT genotype was reported to confer significant elevated tHcy level compared with CC genotype in a population-based CoLaus study . In our research, rs9651118 was associated with risk of T2DM, and we also found that carriers of TT genotype had nonsignificant higher level of tHcy relative to CC genotype carriers, which was in accordance with the previous study.
In conclusion, the present study suggests that MTHFR rs9651118 was associated with T2DM risk and MTHFR rs1801133 may affect serum tHcy levels.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Y. Bi, T. Wang, M. Xu et al., “Advanced research on risk factors of type 2 diabetes,” Diabetes/Metabolism Research and Reviews, vol. 28, supplement 2, pp. 32–39, 2012.View at: Google Scholar
S. Mao, W. Xiang, S. Huang, and A. Zhang, “Association between homocysteine status and the risk of nephropathy in type 2 diabetes mellitus,” Clinica Chimica Acta, vol. 431, pp. 206–210, 2014.View at: Google Scholar
E. K. Hoogeveen, P. J. Kostense, P. J. Beks et al., “Hyperhomocysteinemia is associated with an increased risk of cardiovascular disease, especially in non-insulin-dependent diabetes mellitus: a population-based study,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 18, no. 1, pp. 133–138, 1998.View at: Publisher Site | Google Scholar
T. Huang, J. Ren, and J. Huang, “Association of homocysteine with type 2 diabetes: a meta-analysis implementing Mendelian randomization approach,” BMC Genomics, vol. 14, article 867, 2013.View at: Google Scholar
M. O. Ebesunun and E. O. Obajobi, “Elevated plasma homocysteine in type 2 diabetes mellitus: a risk factor for cardiovascular diseases,” The Pan African medical journal, vol. 12, p. 48, 2012.View at: Google Scholar
K. Shahzad, A. Hai, A. Ahmed, N. Kizilbash, and J. Alruwaili, “A structured-based model for the decreased activity of Ala222Val and Glu429Ala Methylenetetrahydrofolate Reductase (MTHFR) mutants,” Bioinformation, vol. 9, no. 18, pp. 929–936, 2013.View at: Google Scholar
P. Li and C. Qin, “Methylenetetrahydrofolate reductase (MTHFR) gene polymorphisms and susceptibility to ischemic stroke: a meta-analysis,” Gene, vol. 535, no. 2, pp. 359–364, 2014.View at: Google Scholar
B. Yang, S. Fan, and X. Zhi, “Associations of MTHFR gene polymorphisms with hypertension and hypertension in pregnancy: a meta-analysis from 114 studies with 15411 cases and 21970 controls,” PLoS ONE, vol. 9, no. 2, Article ID e87497, 2014.View at: Google Scholar
G. Friedman, N. Goldschmidt, Y. Friedlander et al., “A common mutation A1298C in human methylenetetrahydrofolate reductase gene: association with plasma total homocysteine and folate concentrations,” Journal of Nutrition, vol. 129, no. 9, pp. 1656–1661, 1999.View at: Google Scholar
N. R. AbdRaboh, S. Badr, and S. Ali, “Prevalence of methylenetetrahydrofolate reductase C677T and A1298C polymorphisms in Egyptian patients with type 2 diabetes mellitus,” Egyptian Journal of Medical Human Genetics, vol. 14, no. 1, pp. 87–93, 2013.View at: Google Scholar
H. Benrahma, O. Abidi, L. Melouk et al., “Association of the C677T polymorphism in the human methylenetetrahydrofolate reductase (MTHFR) gene with the genetic predisposition for type 2 diabetes mellitus in a Moroccan population,” Genetic Testing and Molecular Biomarkers, vol. 16, no. 5, pp. 383–387, 2012.View at: Publisher Site | Google Scholar
Y. Yan, H. Liang, S. Yang et al., “Methylenetetrahydrofolate reductase A1298C polymorphism and diabetes risk: evidence from a meta-analysis,” Renal Failure, vol. 15, pp. 1–5, 2014.View at: Google Scholar
S. Narayanan, J. McConnell, J. Little et al., “Associations between two common variants C677T and A1298C in the methylenetetrahydrofolate reductase gene and of measures of folate metabolism and DNA stability (strand breaks, misincorporated uracil, and DNA methylation status) in human lymphocytes in vivo,” Cancer Epidemiology Biomarkers and Prevention, vol. 13, no. 9, pp. 1436–1443, 2004.View at: Google Scholar
A. Yoshimi, B. Aleksic, Y. Kawamura et al., “Gene-wide association study between the methylenetetrahydrofolate reductase gene (MTHFR) and schizophrenia in the Japanese population, with an updated meta-analysis on currently available data,” Schizophrenia Research, vol. 124, no. 1–3, pp. 216–222, 2010.View at: Publisher Site | Google Scholar