Journal of Diabetes Research

Journal of Diabetes Research / 2012 / Article
Special Issue

Pathophysiological Insights into Cardiovascular Health in Metabolic Syndrome

View this Special Issue

Research Article | Open Access

Volume 2012 |Article ID 960318 | 7 pages |

Common Variants of Homocysteine Metabolism Pathway Genes and Risk of Type 2 Diabetes and Related Traits in Indians

Academic Editor: Jun Ren
Received27 May 2011
Revised26 Jul 2011
Accepted01 Aug 2011
Published25 Sep 2011


Hyperhomocysteinemia, a risk factor for cardiovascular disorder, obesity, and type 2 diabetes, is prevalent among Indians who are at high risk of these metabolic disorders. We evaluated association of common variants of genes involved in homocysteine metabolism or its levels with type 2 diabetes, obesity, and related traits in North Indians. We genotyped 90 variants in initial phase (2.115 subjects) and replicated top signals in an independent sample set (2.085 subjects). The variant MTHFR-rs1801133 was the top signal for association with type 2 diabetes ( (95%   –0.92), ) and was also associated with 2 h postload plasma glucose ( ), high-density lipoprotein cholesterol ( ), and total cholesterol ( ) in control subjects. These associations were neither replicated nor significant after meta-analysis. Studies involving a larger study population and different ethnic groups are required before ruling out the role of these important candidate genes in type 2 diabetes, obesity, and related traits.

1. Introduction

Homocysteine, a thiol containing amino acid, has emerged as a determinant for many complex metabolic disorders in the past few years. Hyperhomocysteinemia is a well-known independent risk factor for cardiovascular diseases [1, 2] and is also associated with diabetic complications [36], obesity, [7, 8] and metabolic syndrome [9]. Hyperhomocysteinemia may lead to any of the above metabolic disorders by elicitation of oxidative stress [10], systemic inflammation, [11] and/or endothelial dysfunction [12]. These factors are known to promote insulin resistance and β-cell dysfunction, two important underlying causes for type 2 diabetes [13]. Indians are at high risk for these metabolic disorders, as reflected by high prevalence of type 2 diabetes patients and metabolic syndrome [14, 15]. The observed trend of hyperhomocysteinemia among Indians [16] indicates the possibility of aberrant homocysteine metabolism with increased risk of metabolic disorders in Indian population. Hence, it becomes important to examine the association of variants of genes involved in homocysteine metabolism or modulating plasma homocysteine level with type 2 diabetes and related traits in Indians.

Previous studies suggested contribution of variants in homocysteine metabolism pathway genes in susceptibility to obesity, type 2 diabetes, or related traits [1719]. Methylene tetrahydrofolate reductase (MTHFR) is shown to act synergistically with angiotensin-I-converting enzyme (ACE) to modulate type 2 diabetes risk [17]. The C677T polymorphism (rs1801133) of MTHFR is the most studied genetic variation and is associated with hyperhomocysteinemia [20]. Many studies have also reported association of C677T with type 2 diabetes and related complications [2126]. Variants of MTHFR and other homocysteine metabolism pathway genes like MTR and MTRR have also been shown to be associated with obesity [27]. However, little insight is available in this regard from limited reports showing association of variants in homocysteine metabolism pathway genes with plasma homocysteine levels in Indians [28, 29]. Here, we assessed association of 90 SNPs in 18 genes, from homocysteine metabolism pathway or those associated with homocysteine levels, with type 2 diabetes, obesity, and related traits in 4.200 North Indians.

2. Materials and Methods

2.1. Subjects

The study was performed in two phases: evaluation of 90 SNPs for association with type 2 diabetes and related traits in the first phase and replication of top signal in the second phase. All study participants were unrelated, urban dwellers of Indo-European ethnicity from North India. All nondiabetic controls were recruited through “diabetes awareness camps” conducted in Delhi and National Capital Region, while type 2 diabetes patients were recruited from Endocrinology Clinic, All India Institute of Medical Sciences, New Delhi. Study population for the first phase consisted of 2,115 subjects comprising of 1.073 type 2 diabetes patients (from clinics) and 1.042 nondiabetic controls. Replication phase comprised of an independent sample set of 2.085 subjects including 1.047 type 2 diabetes patients and 1.308 nondiabetic controls. Subjects were characterized as type 2 diabetes patients and controls based on inclusion and exclusion criteria described previously [30]. Briefly, type 2 diabetes was defined by the WHO criteria [31], while criteria for non-diabetic controls were as follows: age ≥40 years, HbA1c <6.0%, fasting plasma glucose <6.11 mmol/L, and no family history of diabetes.

Anthropometric and biochemical estimations were performed as described earlier [30]. The measured variables included height, weight, waist circumferences, hip circumferences, systolic and diastolic blood pressure, body mass index (BMI) and waist to hip ratio (WHR), glycosylated haemoglobin 1c (HbA1c), fasting plasma glucose (FPG), 2h postload plasma glucose (2h PPG), fasting plasma insulin (FPI), C-reactive protein (CRP), C-peptide, total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low density lipoprotein cholesterol (LDL-C), triglyceride (TG), blood urea, and creatinine. Subjects with BMI < 23 kg/m2 were classified as normal weight (NW), whereas subjects with BMI ≥ 23 kg/m2 were considered as overweight/obese (OW/OB), according to WHO Expert Committee recommendations for Asians [32]. Homeostasis model assessments of insulin resistance (HOMA-IR), β-cell function (HOMA-B), and insulin sensitivity (HOMA-S) were obtained using formulae as reported previously [33, 34].

The study was approved by ethics committee of participating institutions and was conducted in accordance with the principles of Helsinki Declaration. Written informed consent was obtained from each study subject.

2.2. Genotyping

We selected 79 SNPs from 15 genes which are involved in homocysteine metabolism pathway and 11 SNPs from 2 genes (NOS3 and ACE) reported to be associated with homocysteine levels [17, 35] (see Table  1 in Supplementary Material available online at: 10.1155/2012/960318). SNPs in functional regions, with minor allele frequency (MAF) >0.05, polymorphic in at least two populations and tag SNPs were preferred. Genotyping in the first phase was performed using Illumina Golden Gate assay. Stringent quality control (QC) criteria were applied as provided previously [36]. Also, 7.2% (147 samples) of total samples were genotyped in duplicate, and an error rate < 0.01 was figured. Of the 78 QC passed SNPs, those with MAF <0.05 or deviating from Hardy-Weinberg equilibrium (HWE) in control subjects were also excluded. Consequently, a total of 67 SNPs were considered for further association analyses.

In replication phase, we genotyped the top signal of first phase, rs1801133 (MTHFR), using iPLEX on a MassARRAY System (Sequenom, San Deigo, Calif, USA). The average genotyping success rate was 97% with >99.8% consistency rate in genotype calls in 5% duplicates.

2.3. Statistical Analysis

Genotype distributions were analyzed for deviation from HWE using χ2 test. Association of SNPs with type 2 diabetes and obesity was determined by logistic regression assuming log additive model. Linear regression was performed under log additive model for association analysis of SNPs with quantitative traits in control subjects. The analyses were adjusted for sex, age, and BMI as appropriate. Odds ratios (ORs), β, and 95% confidence intervals (CI,) were calculated with respect to minor allele. To assess difference in allele frequency of the two study populations, equality of proportions Z-test was applied. Meta-analysis under fixed and random effect models was done by merging summary statistics of two phases. After correcting for multiple testing (Bonferroni correction), a P value was considered significant (58 independent SNPs and 26 traits). Power of study was calculated assuming 10% disease risk at under log additive model using Quanto ( Statistical analyses were performed using SPSS version 17.0 (SPSS, Chicago, Ill), Stata 10.1 (Stata Corporation, College Station, Tex, United States), and PLINK ( unless specified otherwise.

3. Results

Clinical characters of the study subjects (phase 1 and phase 2) are provided in Table 1. The study had >82.9% power to detect association of variants with allele frequency of 0.20 and genetic effect size of 1.25. Allele frequency of controls and cases of the initial and replication phase for SNP rs1801133 did not differ significantly. The allelic and genotypic distribution of SNPs for type 2 diabetes patients and non-diabetic control subjects recruited in the first phase are provided in Supplementary Table  1, whereas details for normal weight (NW) and overweight/obese (OW/OB) subjects are presented in Supplementary Table  2.

CharacteristicsPhase 1Phase 2
Type 2 diabetes patientsControl SubjectsType 2 diabetes PatientsControl Subjects

N (Men/Women)1019 (592/427)1006 (606/400)1047 (619/428)1038 (516/522)
Age (years)53 (45–62)50 (44–60)55 (49–62)54 (45–64)
BMI (Kg/m2)
Men23.80 (22.00–26.00)23.10 (20.10–5.70)24.82 (22.68–27.76)24.68 (22.35–27.39)
Women26.70 (24.20–29.20)25.00 (21.10–8.50)27.39 (24.56–30.63)26.37 (23.13–29.33)
Men1.00 (0.97–1.03)0.97 (0.92–1.0)0.98 (0.95–1.03)0.97 (0.93–1.01)
Women1.00 (0.97–1.03)0.86 (0.82–0.92)0.93 (0.87–0.97)0.86 (0.82–0.90)
Systolic BP (mmHg)130 (130–140)120 (112–132)130 (122–140)130 (120–140)
Diastolic BP (mmHg)80 (78–90)80 (70–88)82 (80–90)80 (78–90)
HbA1c (%)7.80 (6.50–9.40)5.20 (4.90–5.60)8.20 (6.90–9.60)5.65 (5.33–5.89)
FPG (mmol/L)7.90 (6.40–10.30)4.90 (4.50–5.30)7.79 (6.21–10.27)4.87 (4.43–5.27)
2h PPG (mmol/L)5.57 (4.83–6.30)5.63 (4.88–6.22)
FPI (pmol/L)82.80 (42.0–166.80)32.40 (17.40–57.60)74.40 (26.00–96.00)43.80 (28.20–63.60)
HOMA-IR5.20 (2.30–9.60)1.20 (0.60–2.00)4.30 (2.30–9.30)1.60 (1.00–2.40)
HOMA-B60.25 (25.3–154.51)73.40 (40.70–138.60)61.66 (25.61–131.86)105.27 (65.13–167.32)
HOMA-S0.20 (0.10–0.46)0.86 (0.49–1.69)0.23 (0.11–0.44)0.64 (0.42–1.03)
C-peptide (nmol/L)0.89 (0.56–1.36)0.53 (0.36–0.73)1.05 (0.69–1.60)0.66 (0.50–0.86)
CRP (mg/L)2.20 (0.90–4.70)1.30 (0.60–3.00)1.86 (0.90–3.44)1.61 (0.90–3.04)
TC (mmol/L)4.20 (3.50–5.00)4.40 (3.7–5.10)4.64 (3.86–5.42)4.91 (4.22–5.52)
LDL-C (mmol/L)2.57 (1.99–3.36)2.79 (2.33–3.41)2.73 (2.10–3.42)3.01 (2.49–3.51)
HDL-C (mmol/L)1.03 (0.88–1.22)1.06 (0.88–1.28)1.11 (0.94–1.34)1.24 (1.06–1.46)
TG (mmol/L)1.60 (1.10–2.20)1.30 (1.00–1.80)1.43 (0.98–2.13)1.22 (0.86–1.64)

Values provided are median (interquartile range).
N: number of individuals, BMI: body mass index, WHR: waist to hip ratio, BP: blood pressure, HbA1c: glycosylated haemoglobin 1c, FPG: fasting plasma glucose, 2h PPG: 2h postload plasma glucose, FPI: fasting plasma insulin, HOMA-IR: homeostasis model assessments of insulin resistance, HOMA-B: homeostasis model assessments for β-cell function, HOMA-S: homeostasis model assessments for insulin sensitivity, CRP: C-reactive protein, TC: total cholesterol, LDL-C: low-density lipoprotein cholesterol, HDL-C: high-density lipoprotein-cholesterol, TG: triglyceride.

Association analysis results of all the SNPs with type 2 diabetes are presented in Supplementary Table 3. MTHFR variant rs1801133 showed strongest association with type 2 diabetes in initial phase (OR=0.78 (95% CI 0.67–0.92), P = 0.003) (Table 2). However, this could not be replicated in the second phase (OR=1.01 (95% CI 0.87–1.19), P = 0.86). Meta-analysis also failed to detect any association of rs1801133 with type 2 diabetes ( random effect , fixed effect ) (Table 3).

OR (95% CI)P

Type 2 diabetes (total study population: 1019 cases versus 1006 controls)MTHFRrs1801133C/T0.18/0.210.78 (0.67–0.92)0.003

Type 2 diabetes (Normal weight subjects: 290 cases versus 436 controls)MTHFRrs1801133C/T0.18/0.210.62 (0.46–0.85)0.002
MTHFRrs9651118C/T0.26/0.251.31 (1.01–1.68)0.04
CHDHrs4563403G/A0.10/0.110.66 (0.45–0.95)0.02
CBSrs706208T/C0.41/0.391.30 (1.03–1.63)0.03
MTHFD1Lrs1555179C/T0.26/0.231.36 (1.04–1.78)0.03

Type 2 diabetes (Obese/over-weight subjects: 691cases versus 562 controls)MTHFD1rs1076991A/G0.40/0.440.81 (0.69–0.95)0.01
TCN2rs1801198G/C0.42/0.440.85 (0.72–0.99)0.04

Obesity (562 Obese/over-weight versus 436 normal weights)CHDHrs4563403G/A0.09/0.130.69 (0.52–0.92)0.01
TCN2rs1801198G/C0.46/0.411.24 (1.04–1.48)0.02
MTRrs16834521T/C0.31/0.350.82 (0.68–0.99)0.04

Analyses adjusted for sex, age, and BMI. Analysis performed only in control subjects and adjusted for sex and age, MAF: minor allele frequency, OR: odds ratio, CI: confidence interval.

Trait Initial phase Replication phaseMeta-analysis
OR/β(95% CI)POR/β (95% CI)P QI

Type2 diabetes (all)0.78 (0.67–0.92)0.0031.01 (0.87–1.19)0.860.050.380.900.890.02480.44
Type2 diabetes (NW)0.62 (0.46–0.85)0.0021.11(0.80–1.54)0.530.070.520.810.830.01284.30
Type2 diabetes (OW/OB)0.86 (0.71–1.04)0.130.98 (0.82–1.17)0.850.250.250.930.930.320.25
Obesity (NW versus OW/OB)1.07 (0.86–1.33)0.561.11 (0.86–1.43)0.420.330.331.091.090.820
BMI0.32 ((−0.19)–(−0.84))0.220.25 (−0.26–0.76)0.340.
Weight1.41 (0.02–2.80)0.050.43 (−0.95–1.81)0.540.070.070.920.920.330
WC0.47 (−0.76–1.70)0.450.28 (−1.12–1.68)0.700.410.410.390.390.840
WHR−0.003 (−0.01–0.005)0.46−0.005 (−0.01–0.003)−0.004−0.0040.730
FPG0.55 (−0.66–1.77)0.37-0.83 (−2.14–0.47)0.210.840.86−0.09−0.120.1356.98
2h PPG−5.06 ((−9.72)–(−0.40))0.03−1.13 (−4.84–2.58)0.550.070.15−2.65−2.830.2040.15
HbA1c0.02 (−0.03–0.07)0.41−0.03 (−0.07–0.01)0.150.580.79−0.009−0.0070.1259.73
FPI0.43 (−0.44–1.30)0.33−0.77 ((−1.42)–(−0.11))−0.33−0.200.0378.49
HOMA-IR0.11 (−0.09–0.31)0.30−0.20 ((−0.36)–(−0.04))−0.08−0.050.0281.77
HOMA-B70.25 (−34.99–175.50)0.19−652.8 ((−1165)–(−140.1))0.010.440.4941.02−245.940.006886.36
HOMA-S−0.09 (−0.26–0.07)0.270.62 (0.13–1.11)0.010.820.53−
C-peptide0.01 (−0.10–0.12)0.86−0.07 (−0.17–0.04)0.220.440.44−0.03−0.030.320
CRP−0.16 (−0.38–0.05)0.14−0.12 (−0.34–0.1)−0.14−0.140.790
HDL-C2.14 (0.68–3.59)0.004−0.05 (−1.31–1.21)0.940.070.360.891.010.0379.83
LDL-C3.25 (−0.32–6.81)0.072.61 (−1.03–6.26)
TC5.61 (1.28–9.94)0.012.02 (−2.43–6.47)0.370.010.033.863.850.2622.03
TG7.08 (−0.88–15.05)0.08−3.73 (−10.26–2.81)0.260.810.790.621.430.0476.37
Systolic BP1.53 (−0.41–3.48)0.12−0.77 (−2.58–1.03)0.400.660.760.290.350.0965.50
Diastolic BP0.77 (−0.30–1.85)0.16−0.74 (−1.82–0.35)0.180.950.980.020.020.0573.47
Creatinine0.02 (−0.004–0.04)0.11−0.01 (−0.05–0.02)0.500.330.680.010.010.1650.00
Urea0.28 (−0.63–1.20)0.54−0.12 (−1.14–0.9)0.820.760.760.100.100.570
Uric acid0.37 (−0.02–0.76)0.060.05 (−0.1–0.19)0.550.220.310.090.150.1356.68

Analyses adjusted for sex, age, and BMI.
Analysis performed only in control subjects and adjusted for sex and age.
OR: odds ratio, CI: confidence interval, : value for fixed effect, : value for random effect ORf: odds ratio for fixed effect, β f: beta value for fixed effect, ORr: odds ratio for random effect, : beta value for random effect, Q: value for Cochrane's Q statistic, and I: heterogeneity index (0–100).

Our previous studies have indicated influence of obesity status on association of genetic variants with type 2 diabetes in North Indians [30, 37]. Thus, we stratified our subjects into NW and OW/OB groups and assessed association with type 2 diabetes in these subgroups. Among NW subjects, we observed nominal association of MTHFR-rs1801133, MTHFR-rs9651118, CHDH-rs4563403, CBS-rs706208, and MTHFD1L-rs1555179 with type 2 diabetes (P value range= 0.002–0.04) (Table 2). On the other hand, in OW/OB individuals, MTHFD1-rs1076991 ( (95% CI 0.69–0.95) ) and TCN2-rs1801198 ( (95% CI 0.72–0.99), ) showed nominal association with type 2 diabetes. In the second phase, association of MTHFR-rs1801133 with type 2 diabetes in NW subjects could not be replicated either in the second phase or upon Meta-analysis (Table 3).

Further, we investigated association of these SNPs with obesity status (OW/OB versus NW) in non-diabetic control subjects (Supplementary Table  4). Only a nominal association of CHDH-rs4563403 ( (95% CI 0.52–0.92), ), TCN2-rs1801198 ( (95% CI 1.04–1.48), ), and MTR-rs16834521 ( (95% CI 0.68–0.99), ) was observed in the first phase (Table 2). The top signal MTHFR-rs1801133 was not associated with obesity in both first and second phases or after meta-analysis (Tables 2 and 3).

We also assessed association of variants with quantitative traits related to type 2 diabetes and obesity. The first phase control subjects showed only nominal association with the quantitative traits as shown in Supplementary Table  5. Most significant association was observed at AMD1-rs1007274 with LDL-C levels ( (95% CI (−9.61)–(−3.1)), ). MTHFR variant rs1801133 was nominally associated with 2h PPG ( (95% CI (−9.72)–(−0.40)), ), HDL-C ( (95% CI 0.68–3.59), ) and TC ( (95% CI 1.28–9.94), ) in the first phase (Table 3). In the second phase, while none of the above associations were replicated, new associations were observed with FPI ( (95% CI (−1.42)–(−0.11)), ), HOMA-IR ( (95% CI (−0.36)–(−0.04)), ), HOMA-B ( (95% CI (−1165)–(−140.1)), ) and HOMA-S ( (95% CI 0.13–1.11), ). Meta-analysis revealed association of MTHFR-rs1801133 with LDL-C ( , , ) and TC ( , , ). The association analysis results of MTHFR-rs1801133 with all traits (dichotomous and continuous) in the initial and replication phases and after meta-analysis are presented in Table 3. However, none of the associations observed in the present study were significant after correcting for multiple testing.

4. Discussion

Homocysteine plays an important role in cell metabolism as it is involved as an essential intermediate in the transfer of activated methyl groups in the activated methyl cycle. This cycle is responsible for global and gene promoter-specific DNA methylation, an important factor in regulating gene expression [3840]. The biological relevance of homocysteine metabolism and its association with metabolic disorders make it an important candidate pathway for type 2 diabetes. To the best of our knowledge, the present study is first to comprehensively evaluate association of variants of homocysteine metabolism pathway genes as well as genes associated with homocysteine levels with type 2 diabetes, obesity, and related traits.

The C677T (rs1801133) variant of MTHFR is an established variant for plasma homocysteine levels [20, 41] and has also been reported to be associated with type 2 diabetes, its complications, and related traits like LDL-C levels [42, 43]. Association of MTHFR-rs1801133 with type 2 diabetes in different population has been inconsistent [13, 17, 44]. The association studies of MTHFR-rs1801133 with type 2 diabetes and LDL-C levels have also been limited to studies with small sample sizes. In the present study, MTHFR-rs1801133 showed association with type 2 diabetes and 2h PPG, HDL-C, and TC in initial phase but could not be replicated in the second phase. Interestingly, MTHFR-rs1801133 showed nominal association with LDL-C levels in meta-analysis as also suggested in a previous study though this association was not significant after correcting for multiple testing. The effect of MTHFR-rs1801133 variant on type 2 diabetes is most likely to be modulated by influencing levels of homocysteine. However, in the absence of levels of homocysteine in the present study population, this aspect could not be probed.

We also observed association of other investigated variants with quantitative traits. Earlier reports have suggested variants from MTR to be associated with obesity [27], similar to this, we also observed variants of MTR (rs16834521) to be associated with BMI. Variants from this gene were also associated with TC (rs1805087 and rs10737812) and HbA1c (rs1805087) levels. Interestingly, for the first time, we found association of variant rs1007274 of AMD1 with plasma levels of LDL-C, CRP and TC. Though these associations could not reach the threshold after multiple testing corrections, they provide important leads for follow-up studies in larger study populations.

The present study is the largest among the studies investigating association of MTHFR-rs1801133 and was sufficiently powered to detect the previously reported effect sizes. However, we were unable to confirm the associations in Indian population. Differences in genetic architecture and LD pattern due to ethnic differences [17, 21, 44] could be one of the possible reasons for the observed inconsistency in the association of MTHFR-rs1801133 with type 2 diabetes. Previous studies by our group also observed inconsistency in association of variants of TNF-LTA [45] and TNFRSF1B [46] in Indians as compared to other ethnic groups indicating difference in genetic structure.

Population stratification may also give rise to spurious associations in a case-control study. However, cases and controls in our study populations were recruited from geographical locations that form a homogenous cluster as suggested by Indian genome variation consortium [47]. The multidimensional scaling (MDS) plot drawn from the information of 608 independent variants genotyped in the same set of samples indicates genetic homogeneity of the samples (Supplementary Figure  1). Study subjects of present study are also a part of an ongoing genome-wide association study of type 2 diabetes (unpublished data) and have been found to belong to single cluster after performing MDS using genome-wide SNP data. Hence, the subjects of two study populations were genetically similar, but phenotypic difference in terms of biochemical traits was observed between the two study groups which could have influenced association results. It is also known that dietary habits can influence levels of homocysteine and in Indians; individuals with vegetarian diet have been shown to have higher levels of homocysteine [48]. We observe difference in the number of individuals on vegetarian and nonvegetarian diets among both cases and controls of phase1 and phase2 (Supplementary Table  6). We believe this difference in status of vegetarian/nonvegetarian diet would have led to the difference in the association results of the two phases, but additional information on the frequency and kind of non-vegetarian diet along with other dietary information would have best answered this question.

5. Conclusions

In conclusion, variants from genes involved in homocysteine metabolism showed suggestive evidence for association with type 2 diabetes, obesity, and other related traits. However, these associations were not confirmed upon replication. Thus, studies involving a larger study population and different ethnic groups are required before completely ruling out the role of these important candidate genes in type 2 diabetes, obesity, and related traits.


The authors express their gratitude to all the study participants. They are thankful to Mr. A. K. Sharma for his help in sample collection. This study was supported by “Diabetes mellitus-New drug discovery R&D, molecular mechanisms and genetic & epidemiological factors” (NWP0032-4) funded by Council of Scientific and Industrial Research (CSIR), Government of India. G. Chauhan and I. Kaur contributed equally to this work.

Supplementary Materials

The supplementary materials provided contains allelic and genotypic distributions of SNPs in type 2 diabetes patients and non-diabetic control subjects (Supplementary Table 1); allelic and genotypic distribution of SNPs in normal-weight and over-weight/obese subjects (Supplementary Table 2); association results of SNPs with type 2 diabetes (Supplementary Table 3); association results of SNPs with obesity (Supplementary Table 4); association results of SNPs with quantitative traits (Supplementary Table 5); proportion of individuals on a vegetarian and non-vegetarian diet in cases and controls of the two phases (Supplementary Table 6) and multidimensional scaling plot for first phase subjects (Supplementary Figure 1).

  1. Supplementary Materials


  1. L. L. Humphrey, R. Fu, K. Rogers, M. Freeman, and M. Helfand, “Homocysteine level and coronary heart disease incidence: a systematic review and meta-analysis,” Mayo Clinic Proceedings, vol. 83, no. 11, pp. 1203–1212, 2008. View at: Publisher Site | Google Scholar
  2. M. Soinio, J. Marniemi, M. Laakso, S. Lehto, and T. Rönnemaa, “Elevated plasma homocysteine level is an independent predictor of coronary heart disease events in patients with type 2 diabetes mellitus,” Annals of Internal Medicine, vol. 140, no. 2, pp. 94–100, 2004. View at: Google Scholar
  3. J. B. Meigs, P. F. Jacques, J. Selhub et al., “Fasting plasma homocysteine levels in the insulin resistance syndrome: the framingham offspring study,” Diabetes Care, vol. 24, no. 8, pp. 1403–1410, 2001. View at: Google Scholar
  4. L. Brazionis, K. Rowley, C. Itsiopoulos, C. A. Harper, and K. O'Dea, “Homocysteine and diabetic retinopathy,” Diabetes Care, vol. 31, no. 1, pp. 50–56, 2008. View at: Publisher Site | Google Scholar
  5. H. C. Looker, A. Fagot-Campagna, E. W. Gunter et al., “Homocysteine as a risk factor for nephropathy and retinopathy in Type 2 diabetes,” Diabetologia, vol. 46, no. 6, pp. 766–772, 2003. View at: Publisher Site | Google Scholar
  6. M. Buysschaert, A. S. Dramais, P. E. Wallemacq, and M. P. Hermans, “Hyperhomocysteinemia in type 2 diabetes: relationship to macroangiopathy, nephropathy, and insulin resistance,” Diabetes Care, vol. 23, no. 12, pp. 1816–1822, 2000. View at: Google Scholar
  7. P. F. Jacques, A. G. Bostom, P. W. F. Wilson, S. Rich, I. H. Rosenberg, and J. Selhub, “Determinants of plasma total homocysteine concentration in the framingham offspring cohort,” American Journal of Clinical Nutrition, vol. 73, no. 3, pp. 613–621, 2001. View at: Google Scholar
  8. R. Martos, M. Valle, R. Morales, R. Cañete, M. I. Gavilan, and V. Sánchez-Margalet, “Hyperhomocysteinemia correlates with insulin resistance and low-grade systemic inflammation in obese prepubertal children,” Metabolism, vol. 55, no. 1, pp. 72–77, 2006. View at: Publisher Site | Google Scholar
  9. G. Ntaios, C. Savopoulos, S. Chatzopoulos, D. Mikhailidis, and A. Hatzitolios, “Iatrogenic hyperhomocysteinemia in patients with metabolic syndrome: a systematic review and metaanalysis,” Atherosclerosis, vol. 214, no. 1, pp. 11–19, 2011. View at: Publisher Site | Google Scholar
  10. N. Weiss, S. J. Heydrick, O. Postea, C. Keller, J. F. Keaney Jr., and J. Loscalzo, “Influence of hyperhomocysteinemia on the cellular redox state—impact on homocysteine-induced endothelial dysfunction,” Clinical Chemistry and Laboratory Medicine, vol. 41, no. 11, pp. 1455–1461, 2003. View at: Publisher Site | Google Scholar
  11. M. A. Hofmann, E. Lalla, Y. Lu et al., “Hyperhomocysteinemia enhances vascular inflammation and accelerates atherosclerosis in a murine model,” Journal of Clinical Investigation, vol. 107, no. 6, pp. 675–683, 2001. View at: Google Scholar
  12. J. S. Stamler, J. A. Osborne, O. Jaraki et al., “Adverse vascular effects of homocysteine are modulated by endothelium- derived relaxing factor and related oxides of nitrogen,” Journal of Clinical Investigation, vol. 91, no. 1, pp. 308–318, 1993. View at: Google Scholar
  13. Y. Song, N. R. Cook, C. M. Albert, M. Van Denburgh, and J. E. Manson, “Effect of homocysteine-lowering treatment with folic acid and B vitamins on risk of type 2 diabetes in women: a randomized, controlled trial,” Diabetes, vol. 58, no. 8, pp. 1921–1928, 2009. View at: Publisher Site | Google Scholar
  14. S. Wild, G. Roglic, A. Green, R. Sicree, and H. King, “Global prevalence of diabetes: estimates for the year 2000 and projections for 2030,” Diabetes Care, vol. 27, no. 10, pp. 1047–1053, 2004. View at: Publisher Site | Google Scholar
  15. R. Gupta, P. C. Deedwania, A. Gupta, S. Rastogi, R. B. Panwar, and K. Kothari, “Prevalence of metabolic syndrome in an Indian urban population,” International Journal of Cardiology, vol. 97, no. 2, pp. 257–261, 2004. View at: Publisher Site | Google Scholar
  16. N. K. Vikram, R. M. Pandev, R. Sharma, and A. Misra, “Hyperhomocysteinemia in healthy Asian Indians,” American Journal of Hematology, vol. 72, no. 2, pp. 151–152, 2003. View at: Publisher Site | Google Scholar
  17. S. Mehri, N. Koubaa, A. Nakbi et al., “Relationship between genetic polymorphisms of angiotensin-converting enzyme and methylenetetrahydrofolate reductase as risk factors for type 2 diabetes in Tunisian patients,” Clinical Biochemistry, vol. 43, no. 3, pp. 259–266, 2010. View at: Publisher Site | Google Scholar
  18. F. I. V. Errera, M. E. R. Silva, E. Yeh et al., “Effect of polymorphisms of the MTHFR and APOE genes on susceptibility to diabetes and severity of diabetic retinopathy in Brazilian patients,” Brazilian Journal of Medical and Biological Research, vol. 39, no. 7, pp. 883–888, 2006. View at: Publisher Site | Google Scholar
  19. J. T. Bazzaz, M. Shojapoor, H. Nazem et al., “Methylenetetrahydrofolate reductase gene polymorphism in diabetes and obesity,” Molecular Biology Reports, vol. 37, no. 1, pp. 105–109, 2010. View at: Publisher Site | Google Scholar
  20. P. Frosst, H. J. Blom, R. Milos et al., “A candidate genetic risk factor for vascular disease: a common mutation in methylenetetrahydrofolate reductase,” Nature Genetics, vol. 10, no. 1, pp. 111–113, 1995. View at: Google Scholar
  21. P. Beneš, K. Kaňková, J. Mužík et al., “Methylenetetrahydrofolate reductase polymorphism, type II diabetes mellitus, coronary artery disease, and essential hypertension in the Czech population,” Molecular Genetics and Metabolism, vol. 73, no. 2, pp. 188–195, 2001. View at: Publisher Site | Google Scholar
  22. A. R. Chen, H. G. Zhang, Z. P. Wang et al., “C-reactive protein, vitamin B12 and C677T polymorphism of N-5,10-methylenetetrahydrofolate reductase gene are related to insulin resistance and risk factors for metabolic syndrome in Chinese population,” Clinical and Investigative Medicine, vol. 33, no. 5, pp. E290–E297, 2010. View at: Google Scholar
  23. R. Nemr, R. A. Salman, L. H. Jawad, E. A. Juma, S. H. Keleshian, and W. Y. Almawi, “Differential contribution of MTHFR C677T variant to the risk of diabetic nephropathy in Lebanese and Bahraini Arabs,” Clinical Chemistry and Laboratory Medicine, vol. 48, no. 8, pp. 1091–1094, 2010. View at: Publisher Site | Google Scholar
  24. M. Rahimi, A. Hasanvand, Z. Rahimi et al., “Synergistic effects of the MTHFR C677T and A1298C polymorphisms on the increased risk of micro- and macro-albuminuria and progression of diabetic nephropathy among Iranians with type 2 diabetes mellitus,” Clinical Biochemistry, vol. 43, no. 16-17, pp. 1333–1339, 2010. View at: Publisher Site | Google Scholar
  25. G. T. Russo, A. Di Benedetto, D. Magazzù et al., “Mild hyperhomocysteinemia, C677T polymorphism on methylenetetrahydrofolate reductase gene and the risk of macroangiopathy in type 2 diabetes: a prospective study,” Acta Diabetologica, vol. 48, no. 2, pp. 95–101, 2011. View at: Publisher Site | Google Scholar
  26. J. Z. Sun, Y. Xu, H. Lu, and Y. Zhu, “Polymorphism of the methylenetetrahydrofolate reductase gene association with homocysteine and ischemic stroke in type 2 diabetes,” Neurology India, vol. 57, no. 5, pp. 589–593, 2009. View at: Publisher Site | Google Scholar
  27. I. Terruzzi, P. Senesi, I. Fermo, G. Lattuada, and L. Luzi, “Are genetic variants of the methyl group metabolism enzymes risk factors predisposing to obesity?” Journal of Endocrinological Investigation, vol. 30, no. 9, pp. 747–753, 2007. View at: Google Scholar
  28. J. Kumar, G. Garg, A. Kumar et al., “Single nucleotide polymorphisms in homocysteine metabolism pathway genes: association of CHDH A119C and MTHFR C677T with hyperhomocysteinemia,” Circulation, vol. 2, no. 6, pp. 599–606, 2009. View at: Publisher Site | Google Scholar
  29. U. K. Misra, J. Kalita, A. K. Srivastava, and S. Agarwal, “MTHFR gene polymorphism and its relationship with plasma homocysteine and folate in a North Indian population,” Biochemical Genetics, vol. 48, no. 3-4, pp. 229–235, 2010. View at: Publisher Site | Google Scholar
  30. R. Tabassum, S. Chavali, O. P. Dwivedi, N. Tandon, and D. Bharadwaj, “Genetic variants of FOXA2: risk of type 2 diabetes and effect on metabolic traits in North Indians,” Journal of Human Genetics, vol. 53, no. 11-12, pp. 957–965, 2008. View at: Publisher Site | Google Scholar
  31. J. R. Gavin, K. G. M. M. Alberti, M. B. Davidson et al., “Report of the expert committee on the diagnosis and classification of diabetes mellitus,” Diabetes Care, vol. 26, no. 1, pp. S5–S20, 2003. View at: Google Scholar
  32. C. Barba, T. Cavalli-Sforza, J. Cutter et al., “Appropriate body-mass index for Asian populations and its implications for policy and intervention strategies,” The Lancet, vol. 363, no. 9403, pp. 157–163, 2004. View at: Publisher Site | Google Scholar
  33. D. R. Matthews, J. P. Hosker, A. S. Rudenski, B. A. Naylor, D. F. Treacher, and R. C. Turner, “Homeostasis model assessment: insulin resistance and β-cell function from fasting plasma glucose and insulin concentrations in man,” Diabetologia, vol. 28, no. 7, pp. 412–419, 1985. View at: Google Scholar
  34. J. C. Levy, D. R. Matthews, and M. P. Hermans, “Correct homeostasis model assessment (HOMA) evaluation uses the computer program,” Diabetes Care, vol. 21, no. 12, pp. 2191–2192, 1998. View at: Google Scholar
  35. K. S. Brown, L. A. J. Kluijtmans, I. S. Young et al., “Genetic evidence that nitric oxide modulates homocysteine: the NOS3 894TT genotype is a risk factor for hyperhomocystenemia,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 23, no. 6, pp. 1014–1020, 2003. View at: Publisher Site | Google Scholar
  36. G. Chauhan, C. J. Spurgeon, R. Tabassum et al., “Impact of common variants of PPARG, KCNJ11, TCF7L2, SLC30A8, HHEX, CDKN2A, IGF2BP2, and CDKAL1 on the risk of type 2 diabetes in 5,164 Indians,” Diabetes, vol. 59, no. 8, pp. 2068–2074, 2010. View at: Publisher Site | Google Scholar
  37. R. Tabassum, A. Mahajan, G. Chauhan et al., “Evaluation of DOK5 as a susceptibility gene for type 2 diabetes and obesity in North Indian population,” BMC Medical Genetics, vol. 11, article 35, 2010. View at: Publisher Site | Google Scholar
  38. P. Yi, S. Melnyk, M. Pogribna, I. P. Pogribny, R. J. Hine, and S. J. James, “Increase in plasma homocysteine associated with parallel increases in plasma S-adenosylhomocysteine and lymphocyte DNA hypomethylation,” Journal of Biological Chemistry, vol. 275, no. 38, pp. 29318–29323, 2000. View at: Publisher Site | Google Scholar
  39. S. Bleich, B. Lenz, M. Ziegenbein et al., “Epigenetic DNA hypermethylation of the HERP gene promoter induces down-regulation of its mRNA expression in patients with alcohol dependence,” Alcoholism Research, vol. 30, no. 4, pp. 587–591, 2006. View at: Publisher Site | Google Scholar
  40. R. Thaler, M. Agsten, S. Spitzer et al., “Homocysteine suppresses the expression of the collagen cross-linker lysyl oxidase involving IL-6, Fli1, and epigenetic DNA methylation,” Journal of Biological Chemistry, vol. 286, no. 7, pp. 5578–5588, 2011. View at: Publisher Site | Google Scholar
  41. G. Paré, D. I. Chasman, A. N. Parker et al., “Novel associations of CPS1, MUT, NOX4, and DPEP1 with plasma Homocysteine in a healthy population a genome-wide evaluation of 13 974 participants in the women's genome health study,” Circulation, vol. 2, no. 2, pp. 142–150, 2009. View at: Publisher Site | Google Scholar
  42. L. Zhang, R. X. Yin, W. Y. Liu et al., “Association of methylenetetrahydrofolate reductase C677T polymorphism and serum lipid levels in the Guangxi Bai Ku Yao and Han populations,” Lipids in Health and Disease, vol. 9, article 123, 2010. View at: Publisher Site | Google Scholar
  43. R. Kawamoto, K. Kohara, Y. Tabara et al., “An association of 5,10-methylenetetrahydrofolate reductase (MTHFR) gene polymorphism and common carotid atherosclerosis,” Journal of Human Genetics, vol. 46, no. 9, pp. 506–510, 2001. View at: Publisher Site | Google Scholar
  44. T. Helfenstein, F. A. H. Fonseca, W. G. M. Relvas et al., “Prevalence of myocardial infarction is related to hyperhomocysteinemia but not influenced by C677T methylenetetrahydrofolate reductase and A2756G methionine synthase polymorphisms in diabetic and non-diabetic subjects,” Clinica Chimica Acta, vol. 355, no. 1-2, pp. 165–172, 2005. View at: Publisher Site | Google Scholar
  45. A. Mahajan, R. Tabassum, S. Chavali et al., “Obesity-dependent association of TNF-LTA locus with type 2 diabetes in North Indians,” Journal of Molecular Medicine, vol. 88, no. 5, pp. 515–522, 2010. View at: Publisher Site | Google Scholar
  46. R. Tabassum, S. Chavali, A. Mahajan et al., “Association analysis of TNFRSF1B polymorphisms with type 2 diabetes and its related traits in North India,” Genomic Medicine, vol. 2, no. 3-4, pp. 93–100, 2008. View at: Publisher Site | Google Scholar
  47. S. K. Brahmachari, P. P. Majumder, M. Mukerji et al., “Genetic landscape of the people of India: a canvas for disease gene exploration,” Journal of Genetics, vol. 87, no. 1, pp. 3–20, 2008. View at: Publisher Site | Google Scholar
  48. J. Kumar, S. K. Das, P. Sharma, G. Karthikeyan, L. Ramakrishnan, and S. Sengupta, “Homocysteine levels are associated with MTHFR A1298C polymorphism in Indian population,” Journal of Human Genetics, vol. 50, no. 12, pp. 655–663, 2005. View at: Publisher Site | Google Scholar

Copyright © 2012 Ganesh Chauhan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

1398 Views | 909 Downloads | 9 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.