Experimental Diabetes ResearchVolume 2012 (2012), Article ID 960318, 7 pagesdoi:10.1155/2012/960318

Research Article

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

Ganesh Chauhan,¹ Ismeet Kaur,¹ Rubina Tabassum,¹ Om Prakash Dwivedi,¹ Saurabh Ghosh,² Nikhil Tandon,³ and Dwaipayan Bharadwaj¹

¹Genomics and Molecular Medicine Unit, CSIR-Institute of Genomics and Integrative Biology, New Delhi 110 007, India²Human Genetics Unit, Indian Statistical Institute, Kolkata 700 108, India³Department of Endocrinology, All India Institute of Medical Sciences, New Delhi 110 029, India

Received 27 May 2011; Revised 26 July 2011; Accepted 1 August 2011

Academic Editor: Jun Ren

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.

Abstract.

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 [3–6], 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 [17–19]. 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 [21–26]. 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/m² were classified as normal weight (NW), whereas subjects with BMI ≥ 23 kg/m² 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 χ ² 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 (http://hydra.usc.edu/gxe/). Statistical analyses were performed using SPSS version 17.0 (SPSS, Chicago, Ill), Stata 10.1 (Stata Corporation, College Station, Tex, United States), and PLINK (http://pngu.mgh.harvard.edu/purcell/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.

Table 1: Anthropometric and clinical characteristics of the study populations.

Characteristics
Phase 1
Phase 2

Type 2 diabetes patients
Control Subjects
Type 2 diabetes Patients
Control 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/m²)

Men
23.80 (22.00–26.00)
23.10 (20.10–5.70)
24.82 (22.68–27.76)
24.68 (22.35–27.39)

Women
26.70 (24.20–29.20)
25.00 (21.10–8.50)
27.39 (24.56–30.63)
26.37 (23.13–29.33)

WHR

Men
1.00 (0.97–1.03)
0.97 (0.92–1.0)
0.98 (0.95–1.03)
0.97 (0.93–1.01)

Women
1.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-IR
5.20 (2.30–9.60)
1.20 (0.60–2.00)
4.30 (2.30–9.30)
1.60 (1.00–2.40)

HOMA-B
60.25 (25.3–154.51)
73.40 (40.70–138.60)
61.66 (25.61–131.86)
105.27 (65.13–167.32)

HOMA-S
0.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).

Table 2: Variants showing association (uncorrected P-value) with type 2 diabetes and obesity.

Trait
Gene
SNP
Allele (major/minor)
MAF (affected/unaffected)
OR (95% CI)
P

Type 2 diabetes (total study population: 1019 cases versus 1006 controls)†
MTHFR
rs1801133
C/T
0.18/0.21
0.78 (0.67–0.92)
0.003

Type 2 diabetes (Normal weight subjects: 290 cases versus 436 controls)†
MTHFR
rs1801133
C/T
0.18/0.21
0.62 (0.46–0.85)
0.002

MTHFR
rs9651118
C/T
0.26/0.25
1.31 (1.01–1.68)
0.04

CHDH
rs4563403
G/A
0.10/0.11
0.66 (0.45–0.95)
0.02

CBS
rs706208
T/C
0.41/0.39
1.30 (1.03–1.63)
0.03

MTHFD1L
rs1555179
C/T
0.26/0.23
1.36 (1.04–1.78)
0.03

Type 2 diabetes (Obese/over-weight subjects: 691cases versus 562 controls)†    
MTHFD1
rs1076991
A/G
0.40/0.44
0.81 (0.69–0.95)
0.01

TCN2
rs1801198
G/C
0.42/0.44
0.85 (0.72–0.99)
0.04

Obesity (562 Obese/over-weight versus 436 normal weights)‡
CHDH
rs4563403
G/A
0.09/0.13
0.69 (0.52–0.92)
0.01

TCN2
rs1801198
G/C
0.46/0.41
1.24 (1.04–1.48)
0.02

MTR
rs16834521
T/C
0.31/0.35
0.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.

Table 3: Association of rs1801133 with all traits in initial phase, replication phase, and meta-analysis.

Trait
Initial phase
Replication phase
Meta-analysis

OR/β(95% CI)
P
OR/β (95% CI)
P

Q
I

Type2 diabetes (all) †
0.78 (0.67–0.92)
0.003
1.01 (0.87–1.19)
0.86
0.05
0.38
0.90
0.89
0.024
80.44

Type2 diabetes (NW) †
0.62 (0.46–0.85)
0.002
1.11(0.80–1.54)
0.53
0.07
0.52
0.81
0.83
0.012
84.30

Type2 diabetes (OW/OB) †
0.86 (0.71–1.04)
0.13
0.98 (0.82–1.17)
0.85
0.25
0.25
0.93
0.93
0.32
0.25

Obesity (NW versus OW/OB) ‡
1.07 (0.86–1.33)
0.56
1.11 (0.86–1.43)
0.42
0.33
0.33
1.09
1.09
0.82
0

BMI‡
0.32 ((−0.19)–(−0.84))
0.22
0.25 (−0.26–0.76)
0.34
0.12
0.12
0.29
0.29
0.84
0

Weight‡
1.41 (0.02–2.80)
0.05
0.43 (−0.95–1.81)
0.54
0.07
0.07
0.92
0.92
0.33
0

WC‡
0.47 (−0.76–1.70)
0.45
0.28 (−1.12–1.68)
0.70
0.41
0.41
0.39
0.39
0.84
0

WHR‡
−0.003 (−0.01–0.005)
0.46
−0.005 (−0.01–0.003)
0.24
0.18
0.18
−0.004
−0.004
0.73
0

FPG†
0.55 (−0.66–1.77)
0.37
-0.83 (−2.14–0.47)
0.21
0.84
0.86
−0.09
−0.12
0.13
56.98

2h PPG†
−5.06 ((−9.72)–(−0.40))
0.03
−1.13 (−4.84–2.58)
0.55
0.07
0.15
−2.65
−2.83
0.20
40.15

HbA1c†
0.02 (−0.03–0.07)
0.41
−0.03 (−0.07–0.01)
0.15
0.58
0.79
−0.009
−0.007
0.12
59.73

FPI†
0.43 (−0.44–1.30)
0.33
−0.77 ((−1.42)–(−0.11))
0.02
0.21
0.73
−0.33
−0.20
0.03
78.49

HOMA-IR†
0.11 (−0.09–0.31)
0.30
−0.20 ((−0.36)–(−0.04))
0.01
0.21
0.73
−0.08
−0.05
0.02
81.77

HOMA-B†
70.25 (−34.99–175.50)
0.19
−652.8 ((−1165)–(−140.1))
0.01
0.44
0.49
41.02
−245.94
0.0068
86.36

HOMA-S†
−0.09 (−0.26–0.07)
0.27
0.62 (0.13–1.11)
0.01
0.82
0.53
−0.02
0.23
0.0067
86.41

C-peptide†
0.01 (−0.10–0.12)
0.86
−0.07 (−0.17–0.04)
0.22
0.44
0.44
−0.03
−0.03
0.32
0

CRP†
−0.16 (−0.38–0.05)
0.14
−0.12 (−0.34–0.1)
0.28
0.07
0.07
−0.14
−0.14
0.79
0

HDL-C†
2.14 (0.68–3.59)
0.004
−0.05 (−1.31–1.21)
0.94
0.07
0.36
0.89
1.01
0.03
79.83

LDL-C†
3.25 (−0.32–6.81)
0.07
2.61 (−1.03–6.26)
0.16
0.02
0.02
2.94
2.94
0.81
0

TC†
5.61 (1.28–9.94)
0.01
2.02 (−2.43–6.47)
0.37
0.01
0.03
3.86
3.85
0.26
22.03

TG†
7.08 (−0.88–15.05)
0.08
−3.73 (−10.26–2.81)
0.26
0.81
0.79
0.62
1.43
0.04
76.37

Systolic BP†
1.53 (−0.41–3.48)
0.12
−0.77 (−2.58–1.03)
0.40
0.66
0.76
0.29
0.35
0.09
65.50

Diastolic BP†
0.77 (−0.30–1.85)
0.16
−0.74 (−1.82–0.35)
0.18
0.95
0.98
0.02
0.02
0.05
73.47

Creatinine†
0.02 (−0.004–0.04)
0.11
−0.01 (−0.05–0.02)
0.50
0.33
0.68
0.01
0.01
0.16
50.00

Urea†
0.28 (−0.63–1.20)
0.54
−0.12 (−1.14–0.9)
0.82
0.76
0.76
0.10
0.10
0.57
0

Uric acid†
0.37 (−0.02–0.76)
0.06
0.05 (−0.1–0.19)
0.55
0.22
0.31
0.09
0.15
0.13
56.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 [38–40]. 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.

Acknowledgments

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.

References

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.