PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of bmcmdBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Medical Genetics
 
BMC Med Genet. 2011; 12: 150.
Published online Nov 21, 2011. doi:  10.1186/1471-2350-12-150
PMCID: PMC3266217
Folate network genetic variation, plasma homocysteine, and global genomic methylation content: a genetic association study
Susan M Wernimont,1 Andrew G Clark,2 Patrick J Stover,1 Martin T Wells,3 Augusto A Litonjua,4 Scott T Weiss,4 J Michael Gaziano,5 Katherine L Tucker,6 Andrea Baccarelli,7,8 Joel Schwartz,7 Valentina Bollati,8 and Patricia A Cassanocorresponding author9
1Division of Nutritional Sciences, Cornell University, Ithaca, NY, USA
2Department of Molecular Biology & Genetics, Cornell University, Ithaca, NY, USA
3Department of Biological Statistics & Computational Biology, Cornell, Ithaca, NY, USA
4Channing Laboratory, Brigham and Women's Hospital, and Harvard Medical School, Boston, MA, USA
5Veterans Administration (VA) Normative Aging Study, VA Boston Healthcare System, and Division of Aging, Brigham & Women's Hospital, Boston, MA, USA
6Department of Health Sciences, Northeastern University, Boston, MA, USA
7Departments of Environmental Health and Epidemiology, Harvard University, Boston, MA, USA
8Center of Molecular and Genetic Epidemiology, Department of Environmental and Occupational Health, Università degli Studi di Milano and IRCCS Fondazione Ca' Granda Ospedale Maggiore Policlinico, Milan, Italy
9209 Savage Hall, Division of Nutritional Sciences, Cornell University, Ithaca, NY, USA
corresponding authorCorresponding author.
Susan M Wernimont: smw38/at/cornell.edu; Andrew G Clark: ac347/at/cornell.edu; Patrick J Stover: pjs13/at/cornell.edu; Martin T Wells: mtw1/at/cornell.edu; Augusto A Litonjua: ALITONJUA/at/PARTNERS.ORG; Scott T Weiss: scott.weiss/at/channing.harvard.edu; J Michael Gaziano: jmgaziano/at/partners.org; Katherine L Tucker: KL.Tucker/at/neu.edu; Andrea Baccarelli: abaccare/at/hsph.harvard.edu; Joel Schwartz: JSCHWRTZ/at/hsph.harvard.edu; Valentina Bollati: abaccare/at/hsph.harvard.edu; Patricia A Cassano: pac6/at/cornell.edu
Received May 27, 2011; Accepted November 21, 2011.
Background
Sequence variants in genes functioning in folate-mediated one-carbon metabolism are hypothesized to lead to changes in levels of homocysteine and DNA methylation, which, in turn, are associated with risk of cardiovascular disease.
Methods
330 SNPs in 52 genes were studied in relation to plasma homocysteine and global genomic DNA methylation. SNPs were selected based on functional effects and gene coverage, and assays were completed on the Illumina Goldengate platform. Age-, smoking-, and nutrient-adjusted genotype--phenotype associations were estimated in regression models.
Results
Using a nominal P ≤ 0.005 threshold for statistical significance, 20 SNPs were associated with plasma homocysteine, 8 with Alu methylation, and 1 with LINE-1 methylation. Using a more stringent false discovery rate threshold, SNPs in FTCD, SLC19A1, and SLC19A3 genes remained associated with plasma homocysteine. Gene by vitamin B-6 interactions were identified for both Alu and LINE-1 methylation, and epistatic interactions with the MTHFR rs1801133 SNP were identified for the plasma homocysteine phenotype. Pleiotropy involving the MTHFD1L and SARDH genes for both plasma homocysteine and Alu methylation phenotypes was identified.
Conclusions
No single gene was associated with all three phenotypes, and the set of the most statistically significant SNPs predictive of homocysteine or Alu or LINE-1 methylation was unique to each phenotype. Genetic variation in folate-mediated one-carbon metabolism, other than the well-known effects of the MTHFR c.665C>T (known as c.677 C>T, rs1801133, p.Ala222Val), is predictive of cardiovascular disease biomarkers.
Folate and other B vitamins play key roles in biologic processes important to health, including DNA synthesis and the generation of cellular methylation potential. Folate status is influenced by both dietary intake and variation in genes encoding folate-related enzymes, and altered folate status due to nutritional or genetic perturbations is associated with adverse outcomes, including birth defects, cardiovascular disease (CVD), and cancer [1].
Elevated plasma homocysteine, a sulfur-containing amino acid by-product of folate metabolism, is a marker of disturbed folate-mediated one-carbon metabolism, and is associated with an increased risk of CVD [2-5]. Homocysteine levels are modulated by nutrition, particularly folate and vitamin B-12 [6], and by genetic variants, including a well-studied SNP in the methylenetetrahydrofolate reductase gene MTHFR c.665C>T (known as c.677 C>T, rs1801133, p.Ala222Val)[7].
The association of homocysteine with CVD is hypothesized to be mediated, in part, by changes in DNA methylation [8]. Folate-mediated one-carbon metabolism is linked to DNA methylation status through regulation of S-adenosylmethionine, the universal methyl donor, and through the activity of enzymes involved in methylation reactions [9,10].
LINE-1 and Alu elements are abundant, transposable elements whose methylation status has been shown to be highly correlated with genome-wide DNA methylation in some studies [11,12]. Atherosclerosis is characterized by global DNA hypomethylation and transposable element methylation levels are associated with heart disease, stroke, and total mortality; reduced LINE-1 methylation was associated with an increased incidence of ischemic heart disease and stroke in the Normative Aging Study (NAS) [13]. These findings contribute to interest in global genomic DNA methylation as a potential biomarker of CVD risk.
Most previous work investigating variation in genes contributing to folate-mediated one-carbon metabolism in relation to homocysteine and genomic methylation phenotypes focused on a small number of candidate genes; however, other enzymes and genes may also be important; thus this study represents both first report and replication efforts. To investigate the genetic and nutritional predictors of homocysteine and methylation phenotypes, this candidate gene study examined variation across the network of genes representing folate-mediated one-carbon metabolism in relation to homocysteine and methylation outcomes. 330 single nucleotide polymorphisms (SNPs) in 52 genes with a role in folate-mediated one-carbon metabolism were studied. The set of genes, the SNP markers, and the nutrients examined in this study were selected to represent the full functional variation of the folate-mediated one carbon metabolic pathway.
Study population
The Veterans' Administration (VA) established the NAS in 1961. 2,280 men aged 21-81 years (mean age of 42 y at study entry) were enrolled in the study on the basis of health criteria; details have been described elsewhere [14,15]. The analyses described herein focus on non-Hispanic white males using data from the subset of men (~ 700) with measurements of homocysteine and global genomic DNA methylation (Alu and LINE-1). This study complied with the Helsinki Declaration and was approved by the following: Brigham and Women's Hospital Human Subjects committee, VA R&D committee, Harvard School of Public Health, Cornell University Committee on Human Subjects.
DNA extraction, SNP selection and genotyping
Genomic DNA was extracted from stored frozen buffy coat of 7 ml whole blood using the QIAamp DNA Blood Kit (QIAGEN, Valencia, CA). The REPLI-g whole genome amplification kit (QIAGEN) was used to amplify genomic DNA when quantity was insufficient for genotyping.
52 genes that contribute to folate-mediated one-carbon metabolism were identified (Additional file 1). SNP selection encompassed 2 kb on either side of the gene to include promoter and/or regulatory region variants; a total of 384 SNPs were selected. 384 SNPs were submitted to the Center for Inherited Disease Research at the Johns Hopkins University for genotyping via an Illumina GoldenGate custom genotyping panel. Genotype frequencies in controls were compared with those expected in Hardy-Weinberg equilibrium (HWE). Of the 384 SNPs originally submitted, 54 were ultimately excluded, leaving 330 SNPs available for analysis (Additional file 2).
Extensive previously collected data on study participants includes physical measurements, lifestyle factors, and blood assays. Plasma folate, vitamin B-6 (as pyridoxal-5'-phosphate; PLP) and vitamin B-12 were assayed as previously described [16]. Plasma total homocysteine was assayed in the same unselected subset of stored blood samples as plasma folate, vitamin B-6, and vitamin B-12 [16]. The analysis of transposon DNA methylation was reported in prior publications [17,18].
Restricted maximum likelihood and ordinary least squares regression models evaluated the relation between SNPs and the plasma homocysteine and global DNA methylation phenotypes; maximum likelihood regression was used to evaluate epistatic interactions with the dummy-coded MTHFR SNP. Previous work in this cohort demonstrated no population substructure [19], thus no adjustments were made. All regression models were adjusted for age, smoking status, and nutrient residuals (variation in nutrient not predicted by SNP), and an extended model also adjusted for the MTHFR rs1801133 variant (coded as recessive to account for the pattern of association using the fewest model terms). For the homocysteine phenotype, further models tested the interaction of each genotype with the rs1801133 SNP. For all phenotypes, further models tested the interaction of each genotype with the nutrients.
For main effects, regression coefficients with a nominal P ≤ 0.005 were reported, and a False Discovery Rate (FDR) multiple testing correction [20] was applied, with an FDR-adjusted P value significance threshold of 0.05; final models were conditional on a first step that selected the best genetic model for each SNP, thus the FDR is conditional on this first step. For interactions, a less stringent FDR-adjusted P value significance threshold of 0.20 was used. For gene-nutrient interactions, regression coefficients with a nominal P ≤ 0.02 were reported, given few results reached the FDR threshold.
To assess effect modification, product terms between the SNP and the nutrient biomarker residual were included in models. Interactions were captured in a single model term; significance of the interaction was assessed by the P value for the interaction term. Interactions with MTHFR rs1801133, which was dummy-coded, were assessed with the likelihood ratio test (LRT). All statistical analyses were conducted with SAS v. 9.2 (SAS, Cary, NC).
Additional details on methodology are provided in online materials (Additional file 3).
Measurements of the homocysteine phenotype, the Alu element methylation phenotype, and the LINE-1 methylation phenotype were available for 760, 628 and 621 participants, respectively. All had genotype data, 533 men had data on all three phenotypes; each analysis included the maximum number possible. The phenotype groups had similar frequencies for the MTHFR rs1801133 TT genotype, but differed by age and hence differed slightly on age-related variables (Table (Table1).1). The MTHFR rs1801133 TT genotype prevalence in the largest group, the plasma homocysteine group, was 12.2%, similar to the frequency reported in a large North American sample [7].
Table 1
Table 1
Characteristics of Normative Aging Study participants, 1961-2001, with measurements on three phenotypes.
Age and current smoking status were associated with homocysteine (P ≤ 0.001), age was associated with Alu (P ≤ 0.005), and current smoking was associated with LINE-1 (P = 0.055). Folate, vitamin B-6, and vitamin B-12 were associated with homocysteine (P ≤ 0.005), vitamin B-6 was associated with Alu (P ≤ 0.05), and these biomarkers had little or no association with LINE-1. Models exploring the SNP--phenotype association were adjusted for age, smoking, and nutrient residuals. Adjusting for age and smoking made little difference to the coefficients for each SNP. The set of SNPs comprising the most significant associations was nearly identical with or without adjusting for nutrient residuals. Further adjustment for the MTHFR rs1801133 variant made little or no difference to the SNP regression coefficients. The most statistically significant SNPs for each phenotype were relatively common (MAF ≥13%), and the set of most significant SNPs was unique to each phenotype (Tables (Tables2,2, ,3,3, and and44 and Figure Figure11).
Table 2
Table 2
The most statistically significant associations (P ≤ 0.005) between single nucleotide polymorphisms and the plasma homocysteine phenotype a, e
Table 3
Table 3
The most statistically significant associations (P ≤ 0.005) between single nucleotide polymorphisms and the Alu methylation phenotype a, b, d, f
Table 4
Table 4
The most statistically significant association (P ≤ 0.005) between single nucleotide polymorphisms and the LINE-1 methylation phenotype a, b, d, e, f
Figure 1
Figure 1
Manhattan plot. Folate-related SNPs as predictors of plasma total homocysteine and global genomic DNA methylation phenotypes. Models adjusted for age, smoking status, and folate, vitamin B-6, and vitamin B-12 residuals. Horizontal lines represent nominal (more ...)
Total plasma homocysteine phenotype
Of the 20 SNPs with a nominal P ≤ 0.005, five were also significant at the FDR threshold (P ≤ 0.05) (Table (Table2).2). These 5 SNPs comprise 3 genes: formiminotransferase cyclodeaminase (FTCD; 1 SNP, intronic), solute carrier family 19 (folate transporter), member 1 (SLC19A1, 3 SNPs, representing coding nonsynonymous, 5' region, and intronic variants), and solute carrier family 19, member 3 (SLC19A3, 1 SNP, intronic). Genetic variation in all 5 SNPs was positively associated with plasma homocysteine levels, and effects were similar in direction and magnitude (variant genotypes associated with a 4.9-7.2% higher plasma total homocysteine vs. the referent genotype). In each case, the association of the genotype with homocysteine was partially mediated by nutrients; when plasma folate and vitamin B-6 or B-12 biomarkers were added to the models, the regression coefficients were reduced by 29% for FTCD rs2277820, by 43% for SLC19A1 rs1051266, rs1131596, and rs4819130, and by 34% for SLC19A3 rs13007334 (data not shown). A model containing a nonredundant set of 3 of the top 5 FDR-significant SNPs (FTCD rs2277820, SLC19A3 rs13007334, SLC19A1 rs1051266) explained 3.6% of the variation in plasma homocysteine beyond that explained by age, smoking, and folate, B-6, and B-12 residuals (data not shown); the set of 3 SNPs was statistically significant (LRT = 17.6, P = 0.0005, 3 degrees of freedom, df), and the coefficients for each SNP were similar to coefficients from single SNP models. Considering the MTHFR genotype in more detail, the TT genotype group (vs. CC) had elevated homocysteine (nominal P = 0.0052), but the CT genotype had no association with homocysteine (nominal P = 0.8107); thus, the MTHFR genotype did not pass preset FDR thresholds.
In models investigating interactions between each SNP and MTHFR rs1801133, 4 interaction terms were below the FDR threshold (FDR-adjusted P value ≤ 0.2) for the homocysteine phenotype (Additional file 4). No SNP--nutrient (folate, B-6, or B-12) interaction coefficients reached FDR-significance (FDR-adjusted P value ≤ 0.2; Additional file 5). The MTHFR--folate interaction did not reach preset statistical thresholds (pnominal = 0.0578), but the pattern of interaction supported a greater association of MTHFR TT genotype with homocysteine conditional on lower folate status.
Global genomic DNA methylation phenotype: Alu elements
In analyses of the Alu element methylation phenotype, 8 SNPs were statistically significant with a nominal P ≤ 0.005; however, none were statistically significant at the FDR threshold (FDR-adjusted P value ≤ 0.05) (Table (Table3).3). There was little or no mediation of the association by nutrients or plasma homocysteine levels (data not shown). There were no SNP--nutrient interactions with folate or B-12 that reached FDR thresholds for statistical significance (FDR-adjusted P ≤ 0.2) (Additional file 6). Three SNPs had an FDR-significant interaction with plasma vitamin B-6 (Additional file 6); these interactions involved 3 intronic SNPs in 2 genes, aminomethyltransferase (AMT, rs1464567 and rs1464566) and DNA (cytosine-5-)-methyltransferase 3 beta (DNMT3B, rs1883729). Comparing men with the AMT rs1464567 CC/CG genotype to the GG genotype, the mean Alu element methylation was 0.4 SD higher at low B-6, 0.1 SD higher at median B-6, and 0.4 SD lower at high B-6. Comparing men with the AMT rs1464566 GG/GA genotype to the AA genotype, the mean Alu element methylation was 0.4 SD higher at low B-6, 0.1 SD higher at median B-6, and 0.3 SD lower at high B-6. Comparing men with the DNMT3B rs1883729 AA genotype to the AG/GG genotype, the mean Alu element methylation was 0.1 SD lower at low B-6, 0.3 SD higher at median B-6, and 0.8 SD higher at high B-6.
Global genomic DNA methylation phenotype: LINE-1 elements
No SNP main effect associations reached the FDR-significance threshold for LINE-1 methylation (FDR-adjusted P ≤ 0.05; Table Table4).4). There were no SNP--nutrient interactions for folate or B-12 that reached FDR-significance levels (FDR-adjusted P ≤ 0.2) (Additional file 7). An interaction of plasma B-6 with 1 SNP was significant at the FDR threshold of P ≤ 0.2 (rs17080689, an intronic SNP in methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1-like, MTHFD1L) (Additional file 7), suggesting that the relation of the SNP to LINE-1 methylation varied according to plasma levels of vitamin B-6. Comparing participants with the MTHFD1L rs17080689 CA genotype to the CC/AA genotype, mean LINE-1 element methylation was 0.6 SD higher at low B-6, 0.2 SD higher at median B-6, and 0.4 SD lower at high B-6.
We investigated sequence variation in a network of candidate genes involved in one-carbon metabolism in relation to plasma total homocysteine and two measures of global genomic DNA methylation (Alu, LINE-1).
Genes involved in absorption and transport had the most statistically significant associations with the homocysteine phenotype; about 30-40% of the association was mediated through plasma folate and vitamin B-6 and B-12 levels. For the Alu-element methylation phenotype, the top hits were in genes involved in mitochondrial metabolism, nuclear metabolism, and methylation/homocysteine metabolism. For the LINE-1 methylation phenotype, the top SNP was in a gene in the methylation/homocysteine pathway. There was no evidence that nutrient biomarkers mediated the association of SNPs with the methylation phenotypes.
The set of genes represented in the top hits was unique to each phenotype, although pleiotropy was identified for plasma homocysteine and Alu element methylation involving the MTHFD1L and sarcosine dehydrogenase (SARDH) genes.
Plasma total homocysteine phenotype
SLC19A1. There were FDR-significant associations between 3 SNPs in the SLC19A1 gene and plasma total homocysteine; the direction and magnitude of association were similar. Thus, each copy of the coding nonsynonymous rs1051266 A allele, the 5'region rs1131596 C allele, and the intronic rs4819130 C allele was associated with about a 5.0% increase in plasma homocysteine. HapMap plots indicate high LD across the SLC19A1 gene, thus the three SNPs may represent a single effect. The SLC19A1 gene encodes a transporter involved in folate and thiamine uptake and may play a role in intracellular folate distribution [21]. Transporter expression may be regulated by folate status [21]. About half of the association of these three SLC19A1 SNPs with homocysteine was mediated by plasma folate and vitamins B-6/B-12. The nonsynonymous SLC19A1 rs1051266 SNP was previously associated with blood folate levels [22,23], and risk of intracranial aneurysm [24], but not with homocysteine [23,25] or abdominal aortic aneurysm [25]. The 5' region SLC19A1 rs1131596 SNP was associated with reduced RBC folate levels in coronary artery disease patients and decreased SLC19A1 protein expression [26,27]. Genetic variation in SLC19A1 may influence homocysteine levels, mediated by changes in nutrient biomarkers.
FTCD. The intronic FTCD rs2277820 SNP was associated with plasma total homocysteine. The CT genotype group was 7.2% higher on plasma total homocysteine vs. the CC/TT group. FTCD encodes a Golgi-associated enzyme involved in the production of 5,10-methenyl-tetrahydrofolate (THF) [1]. Based on HapMap LD patterns the association with the intronic rs2277820 SNP may proxy variation elsewhere in the gene. Mutations in FTCD are associated with inherited disorders of folate metabolism [28]. 29% of the association between rs2277820 and homocysteine was mediated through plasma folate and vitamins B-6/B-12.
SLC19A3. An FDR-significant association was identified between the intronic rs13007334 SNP in SLC19A3 and plasma total homocysteine. The CT genotype group was 6.9% higher on plasma total homocysteine vs. the CC/TT group. The SLC19A3 gene belongs to the folate transporter family and encodes a thiamine transporter [21]. Although SLC19A3 is not known to transport folate or vitamins B-6/B-12, 34% of the SNP--homocysteine association was mediated by these nutrients. No prior reports link SLC19A3 to biochemical or disease phenotypes, and a biological basis for the link to thiamine metabolism could not be identified.
The variability in homocysteine explained by the model containing the set of the 3 most significant nonredundant SNP hits was 3.6%, a small proportion of the estimated > 50% heritability in homocysteine [29,30], and similar to the proportion explained by age and smoking together.
There were four FDR-significant interactions between studied SNPs and MTHFR rs1801133 (Additional file 4); the most statistically significant was for the ALDH1L1 rs2305230 SNP. In participants with the ALDH1L1 rs2305230 AA genotype, men with 1 copy of the MTHFR rs1801133 T allele had plasma homocysteine 64% higher than men with no copies. However, among participants with the ALDH1L1 rs2305230 AC/CC genotype, men with 1 copy of the MTHFR rs1801133 T allele had plasma homocysteine 2.1% lower than men with no copies.
There were no FDR-significant interactions between studied SNPs and plasma folate, vitamin B-6, or vitamin B-12 for the plasma homocysteine phenotype. The null results may be due to an overly conservative FDR significance threshold, network compensation for genetic and nutritional stresses, or inadequate power to evaluate interactions involving low MAF SNPs; also, the folate status for men in the NAS was relatively high in comparison to national averages as reported in Pfeiffer et al [31], and SNP--nutrient interactions may be attenuated in this range of folate status. The MTHFR rs1801133 SNP, which is expected to interact with folate in predicting the homocysteine phenotype, had a nonsignificant interaction in these data (nominal Pinteraction = 0.0578), but the association of MTHFR with homocysteine was stronger at lower concentrations of plasma folate (data not shown).
A cluster of SNP--vitamin B-6 interactions was noted for variants in the CBS gene, but the P values for these interaction terms were about 0.1 and did not reach thresholds set prior to the analysis. These findings suggest that interactions between vitamin B-6 and genetic variants in the SHMT1 and CBS genes may only be evident with very low vitamin B-6 status, which is consistent with previous work [32,33]. A systematic review of literature published prior to August, 2009 revealed only one report of a statistically significant interaction between genetic variation in SHMT1 (rs1979277) and B-6 [34].
Global genomic DNA methylation phenotype (Alu elements)
There were no FDR-significant main effect associations for the Alu element methylation outcome. None of the SNP--folate or SNP--vitamin B-12 interaction terms reached FDR significance thresholds. Given that the Alu phenotype was measured after the introduction of mandatory folate fortification in the U.S., findings may be limited. Three FDR-significant SNP--vitamin B-6 interactions were identified, including two intronic SNPs in the AMT gene (rs1464567 and rs1464566) and one intronic SNP in the DNMT3B gene (rs1883729). The AMT gene encodes an enzyme that functions in the vitamin B-6-dependent mitochondrial glycine cleavage system [35]. B-6 interactions involving SNPs in GLDC were among the top nominally significant hits for the homocysteine and Alu methylation phenotypes, but did not reach FDR-significance. The DNMT3B gene encodes a DNA methyltransferase enzyme that is localized to the nucleus, developmentally regulated, and functions to establish de novo methylation patterns [36,37]; DNMT3B expression is associated with cancer [36-38]. Although cell culture studies have not supported Alu elements as DNMT3B targets [36,39] in both in vitro and in vivo models, DNMT3b protein levels were down-regulated by B vitamin deficiency (deficiency of folate, B-6, and B-12 together), de novo methylation was suppressed both in vitro and in vivo under conditions of B vitamin deficiency[40], and S-adenosylmethionine levels were markedly decreased in response to lowered B-6 concentrations in culture medium [32] consistent with the direction of association observed here.
Global genomic DNA methylation phenotype (LINE-1 elements)
There were no FDR-significant associations observed for the LINE-1 methylation phenotype. There were no FDR-significant interactions between SNPs and folate or vitamin B-12; the measurement of LINE-1 in Normative Aging Study men took place after the introduction of mandatory folate fortification in the U.S., and limited variation may have limited findings. A single SNP--vitamin B-6 interaction was significant at the FDR threshold for the intronic rs17080689 in the MTHFD1L gene. The MTHFD1L gene product functions downstream from the vitamin B-6-dependent glycine cleavage system [41] and intronic variation in MTHFD1L was previously associated with CVD [42].
Strengths of the present study include investigation of a large cohort with homocysteine data collected prior to the introduction of mandatory folate-fortification in the U.S. Also, SNP selection for the genotyping assay reflected functional, LD, and physical coverage of genes. Using a systematic approach, we identified the best genetic models for each SNP, then tested single SNPs, the interaction of each SNP with MTHFR rs1801133, and the interaction of each SNP with folate, vitamin B-6 and vitamin B-12. Findings were corrected for multiple comparisons and those surpassing the FDR threshold were discussed in more detail. Weaknesses of the study include the fact that methylation (but not homocysteine) measures were collected after the introduction of mandatory folate fortification in the U.S., which may have limited variation in B-vitamin status. Information on additional nutrients such as choline would have allowed a more complete evaluation of gene-nutrient interactions. S-adenosylhomocysteine and/or the ratio of S-adenosylmethionine to S-adenosylhomocysteine are likely to be more sensitive indicators of vascular disease risk than homocysteine [43,44], but were not measured. The methylation phenotypes studied here are believed to be an adequate proxy of genome-wide DNA methylation. Gene-specific methylation data was not available; furthermore, because the folate-mediated one-carbon network functions to generate cellular methylation potential and thus contributes to numerous methylation reactions, it may be more appropriate to evaluate folate network genetic variation in relation to global measures of methylation. Finally, due to genotyping failure, some key variants could not be analyzed, for example rs6922269 in MTHFD1L [42], although proxies were selected purposefully to address this limitation.
The most significant hits for the homocysteine and methylation outcomes reflected genes involved in the generation of one-carbon units, including SLC19A1 and FTCD. Because a unique set of genes was identified for each phenotype, and because the top hits could not be predicted on the basis of hypothesized impact on cellular methylation potential, this work suggests that not all folate effects are mediated through the ratio of S-adenosylmethionine to S-adenosylhomocysteine. Thus, beyond the well-described MTHFR rs1801133 SNP, polymorphisms in other genes make important contributions to homocysteine and global genomic DNA methylation phenotypes. Furthermore, some associations are sensitive to nutritional status of B vitamins. Future work should continue to include a broad evaluation of one-carbon network genetic and nutritional variation in unfortified or pre-fortification populations and extend these findings for CVD biomarkers to an investigation of CVD phenotypes.
% 5-meC: percentage of methylated cytosines; 3': 3' region; 5': 5' region; A: Additive; AHCY: Adenosylhomocysteinase; AHCYL1: Adenosylhomocysteinase-like 1; AHCYL2: Adenosylhomocysteinase-like 2, KIAA0828; ALDH1L1: Aldehyde dehydrogenase 1 family: member L1; AMT: Aminomethyltransferase; ATIC: 5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase/IMP cyclohydrolase; BHMT: Betaine-homocysteine S-methyltransferase; CBS: Cystathionine-beta-synthase; CTH: Cystathionase (cystathionine gamma-lyase); CELF1: CUGBP: Elav-like family member 1; CEPH: Centre d'Etude du Polymorphisme Humain; CV: coefficient of variation; CN: Coding nonsynonymous; CS: Coding synonymous; CVD: cardiovascular disease; df: degrees of freedom; DHFR: Dihydrofolate reductase; DMGDH: Dimethylglycine dehydrogenase; DNMT1: DNA (cytosine-5-)-methyltransferase 1; DNMT3A: DNA (cytosine-5-)-methyltransferase 3 alpha; DNMT3B: DNA (cytosine-5-)-methyltransferase 3 beta; D: Dominant; FDR: False Discovery Rate; FOLH1: Folate hydrolase (prostate-specific membrane antigen) 1; FOLR1: Folate receptor 1 (adult); FOLR2: Folate receptor 2 (fetal); FOLR3: Folate receptor 3 (gamma); FPGS: Folylpolyglutamate synthase; FTCD: Formiminotransferase cyclodeaminase; FTH1: Ferritin: heavy polypeptide 1; GART: Phosphoribosylglycinamide formyltransferase: phosphoribosylglycinamide synthetase: phosphoribosylaminoimidazole synthetase; GCSH: Glycine cleavage system protein H (aminomethyl carrier); GGH: Gamma-glutamyl hydrolase (conjugase: folylpolygammaglutamyl hydrolase); GLDC: Glycine dehydrogenase (decarboxylating); GNMT: Glycine N-methyltransferase; HSPA8: Heat shock 70 kDa protein 8; HWE: Hardy-Weinberg equilibrium; I: Intronic; LRT: likelihood ratio test; LD: linkage disequilibrium; MAF: minor allele frequency; MARS: Methionyl-tRNA synthetase; MAT1A: Methionine adenosyltransferase I: alpha; MAT2A: Methionine adenosyltransferase II: alpha; MAT2B: Methionine adenosyltransferase II: beta; M.E.: Main effect; MTHFD1: Methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1: methenyltetrahydrofolate cyclohydrolase: formyltetrahydrofolate synthetase; MTHFD1L: Methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1-like; MTHFD2: Methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 2: methenyltetrahydrofolate cyclohydrolase; MTHFR: Methylenetetrahydrofolate reductase (NADPH); MTHFS: 5,10-methenyltetrahydrofolate synthetase (5-formyltetrahydrofolate cyclo-ligase); MTR: 5-methyltetrahydrofolate-homocysteine methyltransferase; MTRR: 5-methyltetrahydrofolate-homocysteine methyltransferase reductase; NAS: Normative Aging Study; O: Overdominant; PLP: pyridoxal-5'-phosphate; R: Recessive; SARDH: Sarcosine dehydrogenase; SHMT1: Serine hydroxymethyltransferase 1 (soluble); SHMT2: Serine hydroxymethyltransferase 2 (mitochondrial); SLC19A1: Solute carrier family 19 (folate transporter): member 1; SLC19A2: Solute carrier family 19 (thiamine transporter): member 2; SLC19A3: Solute carrier family 19: member 3; SLC25A32: Solute carrier family 25: member 32; SLC46A1: Solute carrier family 46 (folate transporter): member 1; SNP: single nucleotide polymorphism; TCN1: Transcobalamin I (vitamin B-12 binding protein: R binder family); TCN2: Transcobalamin II; THF: tetrahydrofolate; TYMS: Thymidylate synthetase; UBE2I: Ubiquitin-conjugating enzyme E2I (UBC9 homolog: yeast); UBE2N: Ubiquitin-conjugating enzyme E2N (UBC13 homolog: yeast); VA: Veterans' Administration.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
SMW, PAC, AAL, STW, JMG, KLT, AB, and JS designed research; SMW, PJS, AGC, MTW, VB and PAC conducted research; SMW and PAC analyzed data; SMW, PAC, AGC, PJS, and MTW wrote the paper, and PAC had primary responsibility for all work and final content. All authors read and approved the final manuscript.
Pre-publication history
The pre-publication history for this paper can be accessed here:
Additional file 1
52 genes in the folate-mediated one-carbon pathway.
Additional file 2
330 folate-related SNPs assayed in men in the Normative Aging Study.
Additional file 3
Supplemental methods.
Additional file 4
Epistatic interactions with the MTHFR rs1801133 SNP and plasma homocysteine. The most statistically significant associations (FDR-adjusted Likelihood Ratio Test P ≤ 0.2) for SNP by MTHFR rs1801133 interactions in relation to the plasma homocysteine phenotype for men in the Normative Aging Study.
Additional file 5
Gene-nutrient interactions and plasma homocysteine. The most statistically significant associations (P ≤ 0.02) for SNP by nutrient interactions in relation to the plasma homocysteine phenotype for men in the Normative Aging Study.
Additional file 6
Gene-nutrient interactions and Alu element methylation. The most statistically significant associations (P ≤ 0.02) for SNP by nutrient interactions in relation to the global genomic DNA methylation phenotype (Alu elements) for men in the Normative Aging Study.
Additional file 7
Gene-nutrient interactions and LINE-1 element methylation. The most statistically significant associations (P ≤ 0.02) for SNP by nutrient interactions in relation to the global genomic DNA methylation phenotype (LINE-1 elements) for men in the Normative Aging Study
Acknowledgements
The authors wish to thank Ms. Ellen M. Smith, Cornell '09, for her contributions to preliminary analyses for this work. The project was supported in part by the Center for Vertebrate Genomics (PAC), the Presidents Council for Cornell Women (PAC), and the Bronfenbrenner Life Course Center (PAC), all at Cornell University; by U.S. Department of Agriculture Cooperative State Research, Education, and Extension Service grant 2003-34324-13135 (subproject to PAC); The Veterans Administration (VA) Normative Aging Study is supported by the Department of Veterans Affairs Cooperative Studies Program/ERIC and is a research component of the Massachusetts Veterans Epidemiology and Information Center (MAVERIC); by National Institute of Environmental Health Sciences grant ES014663 (JS); and by T32DK007158 (SMW), National Institute of Diabetes, Digestive and Kidney Diseases. Dr. Baccarelli receives salary support from New Investigator funding from the Harvard School of Public Health-National Institute of Environmental Health Sciences Center for Environmental Health (ES000002). Dr. Schwartz receives support through National Institute of Environmental Health Sciences grant R01ES015172. Genotyping services were provided (to PAC) by the Johns Hopkins University under federal contract number (N01-HV-48195) from the National Heart, Lung, and Blood Institute. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Diabetes And Digestive And Kidney Diseases, the National Heart, Lung, and Blood Institute, or the National Institutes of Health. The funding bodies played no role in study design, the collection, analysis, and interpretation of data, the writing of the manuscript, or in the decision to submit the manuscript for publication.
  • Fox JT, Stover PJ. Folate-mediated one-carbon metabolism. Vitam Horm. 2008;79:1–44. [PubMed]
  • Selhub J. The many facets of hyperhomocysteinemia: studies from the Framingham cohorts. The Journal of nutrition. 2006;136:1726S–1730S. [PubMed]
  • The Homocysteine Studies Collaboration. Homocysteine and risk of ischemic heart disease and stroke: a meta-analysis. JAMA. 2002;288:2015–2022. doi: 10.1001/jama.288.16.2015. [PubMed] [Cross Ref]
  • Lewis SJ, Ebrahim S, Davey Smith G. Meta-analysis of MTHFR 677C->T polymorphism and coronary heart disease: does totality of evidence support causal role for homocysteine and preventive potential of folate? Br Med J. 2005;331:1053. doi: 10.1136/bmj.38611.658947.55. [PMC free article] [PubMed] [Cross Ref]
  • Wald DS, Law M, Morris JK. Homocysteine and cardiovascular disease: evidence on causality from a meta-analysis. Br Med J. 2002;325:1202. doi: 10.1136/bmj.325.7374.1202. [PMC free article] [PubMed] [Cross Ref]
  • Homocysteine Lowering Trialists' Collaboration. Dose-dependent effects of folic acid on blood concentrations of homocysteine: a meta-analysis of the randomized trials. Am J Clin Nutr. 2005;82:806–812. [PubMed]
  • Klerk M, Verhoef P, Clarke R, Blom HJ, Kok FJ, Schouten EG. MTHFR 677C-->T polymorphism and risk of coronary heart disease: a meta-analysis. JAMA. 2002;288:2023–2031. doi: 10.1001/jama.288.16.2023. [PubMed] [Cross Ref]
  • Jamaluddin MS, Yang X, Wang H. Hyperhomocysteinemia, DNA methylation and vascular disease. Clin Chem Lab Med. 2007;45:1660–1666. doi: 10.1515/CCLM.2007.350. [PubMed] [Cross Ref]
  • Ulrey CL, Liu L, Andrews LG, Tollefsbol TO. The impact of metabolism on DNA methylation. Hum Mol Genet. 2005;14(Spec No 1):R139–R147. [PubMed]
  • Pogribny IP, Beland FA. DNA hypomethylation in the origin and pathogenesis of human diseases. Cell Mol Life Sci. 2009;66:2249–2261. doi: 10.1007/s00018-009-0015-5. [PubMed] [Cross Ref]
  • Cordaux R, Batzer MA. The impact of retrotransposons on human genome evolution. Nat Rev Genet. 2009;10:691–703. doi: 10.1038/nrg2640. [PMC free article] [PubMed] [Cross Ref]
  • Yang AS, Estecio MR, Doshi K, Kondo Y, Tajara EH, Issa JP. A simple method for estimating global DNA methylation using bisulfite PCR of repetitive DNA elements. Nucleic acids research. 2004;32:e38. doi: 10.1093/nar/gnh032. [PMC free article] [PubMed] [Cross Ref]
  • Baccarelli A, Wright R, Bollati V, Litonjua A, Zanobetti A, Tarantini L, Ischemic Heart Disease and Stroke in Relation to Blood DNA Methylation. Epidemiology (Cambridge, Mass) 2010. [PubMed]
  • Bell B, Rose CL, Damon A. The Veterans Administration longitudinal study of healthy aging. The Gerontologist. 1966;6:179–184. doi: 10.1093/geront/6.4.179. [PubMed] [Cross Ref]
  • Damon A, Seltzer CC, Stoudt HW, Bell B. Age and physique in health white veterans at Boston. Journal of gerontology. 1972;27:202–208. [PubMed]
  • Tucker KL, Qiao N, Scott T, Rosenberg I, Spiro A III. High homocysteine and low B vitamins predict cognitive decline in aging men: the Veterans Affairs Normative Aging Study. Am J Clin Nutr. 2005;82:627–635. [PubMed]
  • Baccarelli A, Wright RO, Bollati V, Tarantini L, Litonjua AA, Suh HH. et al. Rapid DNA methylation changes after exposure to traffic particles. Am J Respir Crit Care Med. 2009;179:572–578. doi: 10.1164/rccm.200807-1097OC. [PMC free article] [PubMed] [Cross Ref]
  • Bollati V, Schwartz J, Wright R, Litonjua A, Tarantini L, Suh H. et al. Decline in genomic DNA methylation through aging in a cohort of elderly subjects. Mech Ageing Dev. 2009;130:234–239. doi: 10.1016/j.mad.2008.12.003. [PMC free article] [PubMed] [Cross Ref]
  • Wilker EH, Alexeeff SE, Poon A, Litonjua AA, Sparrow D, Vokonas PS. et al. Candidate genes for respiratory disease associated with markers of inflammation and endothelial dysfunction in elderly men. Atherosclerosis. 2009;206:480–485. doi: 10.1016/j.atherosclerosis.2009.03.004. [PMC free article] [PubMed] [Cross Ref]
  • Benjamini Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society, Series B (Methodological) 1995;57:289–300.
  • Ganapathy V, Smith SB, Prasad PD. SLC19: the folate/thiamine transporter family. Pflugers Arch. 2004;447:641–646. doi: 10.1007/s00424-003-1068-1. [PubMed] [Cross Ref]
  • Stanislawska-Sachadyn A, Mitchell LE, Woodside JV, Buckley PT, Kealey C, Young IS. et al. The reduced folate carrier (SLC19A1) c.80G>A polymorphism is associated with red cell folate concentrations among women. Ann Hum Genet. 2009;73:484–491. doi: 10.1111/j.1469-1809.2009.00529.x. [PMC free article] [PubMed] [Cross Ref]
  • Fredriksen A, Meyer K, Ueland PM, Vollset SE, Grotmol T, Schneede J. Large-scale population-based metabolic phenotyping of thirteen genetic polymorphisms related to one-carbon metabolism. Hum Mutat. 2007;28:856–865. doi: 10.1002/humu.20522. [PubMed] [Cross Ref]
  • Semmler A, Linnebank M, Krex D, Gotz A, Moskau S, Ziegler A. et al. Polymorphisms of homocysteine metabolism are associated with intracranial aneurysms. Cerebrovasc Dis. 2008;26:425–429. doi: 10.1159/000155638. [PubMed] [Cross Ref]
  • Giusti B, Saracini C, Bolli P, Magi A, Sestini I, Sticchi E. et al. Genetic analysis of 56 polymorphisms in 17 genes involved in methionine metabolism in patients with abdominal aortic aneurysm. J Med Genet. 2008;45:721–730. doi: 10.1136/jmg.2008.057851. [PubMed] [Cross Ref]
  • Chatzikyriakidou A, Vakalis KV, Kolaitis N, Kolios G, Naka KK, Michalis LK. et al. Distinct association of SLC19A1 polymorphism -43T>C with red cell folate levels and of MTHFR polymorphism 677C>T with plasma folate levels. Clin Biochem. 2008;41:174–176. doi: 10.1016/j.clinbiochem.2007.11.006. [PubMed] [Cross Ref]
  • Chatzikyriakidou A, Georgiou I, Voulgari PV, Papadopoulos CG, Tzavaras T, Drosos AA. Transcription regulatory polymorphism -43T>C in the 5'-flanking region of SLC19A1 gene could affect rheumatoid arthritis patient response to methotrexate therapy. Rheumatol Int. 2007;27:1057–1061. doi: 10.1007/s00296-007-0339-0. [PubMed] [Cross Ref]
  • Hilton JF, Christensen KE, Watkins D, Raby BA, Renaud Y, de la Luna S. et al. The molecular basis of glutamate formiminotransferase deficiency. Human mutation. 2003;22:67–73. doi: 10.1002/humu.10236. [PubMed] [Cross Ref]
  • Siva A, De Lange M, Clayton D, Monteith S, Spector T, Brown MJ. The heritability of plasma homocysteine, and the influence of genetic variation in the homocysteine methylation pathway. QJM. 2007;100:495–499. doi: 10.1093/qjmed/hcm054. [PubMed] [Cross Ref]
  • Nilsson SE, Read S, Berg S, Johansson B. Heritabilities for fifteen routine biochemical values: findings in 215 Swedish twin pairs 82 years of age or older. Scand J Clin Lab Invest. 2009;69:562–569. doi: 10.1080/00365510902814646. [PubMed] [Cross Ref]
  • Pfeiffer CM, Johnson CL, Jain RB, Yetley EA, Picciano MF, Rader JI. et al. Trends in blood folate and vitamin B-12 concentrations in the United States, 1988 2004. Am J Clin Nutr. 2007;86:718–727. [PubMed]
  • Perry C, Yu S, Chen J, Matharu KS, Stover PJ. Effect of vitamin B6 availability on serine hydroxymethyltransferase in MCF-7 cells. Arch Biochem Biophys. 2007;462:21–27. doi: 10.1016/j.abb.2007.04.005. [PMC free article] [PubMed] [Cross Ref]
  • Taoka S, West M, Banerjee R. Characterization of the heme and pyridoxal phosphate cofactors of human cystathionine beta-synthase reveals nonequivalent active sites. Biochemistry. 1999;38:7406. doi: 10.1021/bi995077i. [PubMed] [Cross Ref]
  • Koushik A, Kraft P, Fuchs CS, Hankinson SE, Willett WC, Giovannucci EL. et al. Nonsynonymous polymorphisms in genes in the one-carbon metabolism pathway and associations with colorectal cancer. Cancer Epidemiol Biomarkers Prev. 2006;15:2408–2417. doi: 10.1158/1055-9965.EPI-06-0624. [PubMed] [Cross Ref]
  • Lamers Y, Williamson J, Gilbert LR, Stacpoole PW, Gregory JF III. Glycine turnover and decarboxylation rate quantified in healthy men and women using primed, constant infusions of [1,2-(13)C2]glycine and [(2)H3]leucine. The Journal of nutrition. 2007;137:2647–2652. [PubMed]
  • Weisenberger DJ, Velicescu M, Cheng JC, Gonzales FA, Liang G, Jones PA. Role of the DNA methyltransferase variant DNMT3b3 in DNA methylation. Mol Cancer Res. 2004;2:62–72. [PubMed]
  • Turek-Plewa J, Jagodzinski PP. The role of mammalian DNA methyltransferases in the regulation of gene expression. Cell Mol Biol Lett. 2005;10:631–647. [PubMed]
  • Robertson KD, Uzvolgyi E, Liang G, Talmadge C, Sumegi J, Gonzales FA. et al. The human DNA methyltransferases (DNMTs) 1, 3a and 3b: coordinate mRNA expression in normal tissues and overexpression in tumors. Nucleic acids research. 1999;27:2291–2298. doi: 10.1093/nar/27.11.2291. [PMC free article] [PubMed] [Cross Ref]
  • Choi SH, Heo K, Byun HM, An W, Lu W, Yang AS. Identification of preferential target sites for human DNA methyltransferases. Nucleic acids research. 2010. [PMC free article] [PubMed]
  • Fuso A, Nicolia V, Cavallaro RA, Scarpa S. DNA methylase and demethylase activities are modulated by one-carbon metabolism in Alzheimer's disease models. J Nutr Biochem. 2010. [PubMed]
  • Christensen KE, MacKenzie RE. Mitochondrial methylenetetrahydrofolate dehydrogenase, methenyltetrahydrofolate cyclohydrolase, and formyltetrahydrofolate synthetases. Vitam Horm. 2008;79:393–410. [PubMed]
  • Wellcome Trust Case Control Consortium. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007;447:661–678. doi: 10.1038/nature05911. [PMC free article] [PubMed] [Cross Ref]
  • Kerins DM, Koury MJ, Capdevila A, Rana S, Wagner C. Plasma S-adenosylhomocysteine is a more sensitive indicator of cardiovascular disease than plasma homocysteine. Am J Clin Nutr. 2001;74:723–729. [PubMed]
  • Wagner C, Koury MJ. S-Adenosylhomocysteine: a better indicator of vascular disease than homocysteine? Am J Clin Nutr. 2007;86:1581–1585. [PubMed]
Articles from BMC Medical Genetics are provided here courtesy of
BioMed Central