118 SNPs of folate-related genes and risks of spina bifida and conotruncal heart defects

Periconceptional vitamin supplementation with folic acid substantially reduces risks of women having neural tube defect-affected pregnancies [1, 2] and has been implicated in reducing risks of several other congenital anomalies, including orofacial clefts and selected heart defects [3–11]. Mechanisms underlying these reduced risks have not been elucidated, although it has been speculated that supplementation with vitamins containing folic acid restores some normal developmental function that is genetically compromised in selected infants.

Investigating genetic variation that influences cellular absorption, transport, and metabolism of folate may offer insight into this unknown developmentally protective mechanism. Indeed, numerous investigations of genes that are specifically involved with folate metabolism have yielded at least one gene, 5, 10-methylenetetrahydrofolate reductase (MTHFR), that has been associated with a modest increased risk of neural tube defects (e.g., [12–17]), and possibly heart defects [18, 19]. Observed risks with the two principal MTHFR variants, however, do not appear to account for a large proportion of the etiologic fraction of any of these defects, under the assumption that MTHFR variants have a causal role [17]. Thus, further investigation of other folate-related genes is necessary to reveal clues about mechanisms underlying the potential embryonic protective effects of folic acid supplementation.

We hypothesized that genetic susceptibility of fetal metabolism or transport of folate puts fetuses at risk for selected congenital anomalies. Using population-based data, we investigated 118 single nucleotide polymorphisms (SNPs) in 14 genes in the complex folate pathway as risk factors for spina bifida and conotruncal heart defects.

This population-based case-control study included infants with spina bifida or conotruncal heart defects diagnosed within 1 year after birth among infants and fetal deaths delivered to women residing in most California counties. Data were derived from the California Birth Defects Monitoring Program [20], a population-based active surveillance system for collecting information on infants and fetuses with congenital malformations. Diagnostic and demographic information was collected by program staff from multiple sources of medical records for all liveborn and stillborn fetuses (defined as >20 weeks gestation). Overall ascertainment for major malformations has been estimated as 97% complete [21]. Eligible were live born infants only because the source of DNA was from newborn screening cards.

Included were 259 infants with spina bifida and a random sample of 359 nonmalformed control infants born during 1983–86 and 1994–95 in selected counties in California. Also included for study were 214 infants with conotruncal heart defects, specifically d-transposition of the great arteries and tetralogy of Fallot. The random sample of 1983–86 controls for conotruncal heart defects included 220 of the overall 359. Newborn bloodspots were obtained from the State of California and their use in this study was consistent with the consent procedures at the time of sample collection. The protocol for this study was reviewed and approved by the State of California Health and Welfare Agency Committee for the Protection of Human Subjects.

Genomic DNA was extracted from dried blood spots on filter paper using the Puregene DNA Extraction Kit (Gentra, Minneapolis, MN). Prior to genotyping, genomic DNA was amplified using a commercial multiple displacement amplification (MDA) kit, GenomePhi (GE Healthcare, Piscataway, NJ). The MDA method relies on isothermal amplification using the DNA polymerase of the bacteriophage phi29 and is a recently developed technique for high performance WGA. MDA has been demonstrated to be reliable for genotyping, with the most favorable call rates, best genomic coverage, and lowest amplification bias [22]. Studies indicate no discernable difference between WGA samples with GenomiPhi kit and the original DNA templates [23, 24]. The whole genome amplification (WGA) product was then quantified using RNase P method (AppliedBiosystems, Foster City, CA). 150 ng WGA product was then used for each SNPlex assay pool which contained about 48 SNPs.

Genotype analyses were performed using SNPlex assays (AppliedBiosystems, Foster City, CA). SNP markers were selected using the SNPBrowser™ program (version 3.0) provided by AppliedBiosystems Inc. This program allowed selection of SNP markers from the HapMap database. For each target gene, tagging SNPs were selected based on the pairwise r² > = 0.8. SNPs with minor allele frequencies lower than 10% in Caucasians were excluded. All validated non-synonymous SNPs were included. Successful rates for SNPlex assays were >96% for 75 SNPs, from 90% to 96% for 32 SNPs, from 70% to 90% for 7 SNPs. 15 SNPs suffered from more than 30% failure rates. In a subsequent effort to fill in the missing genotyping data and obtain higher call rate, we performed TaqMan SNP assays (Appliedbiosystems, Foster City, CA) for 22 of these SNPs on an ABI 7900 Genetic Analyzer.

Genotypes among controls were analyzed to verify that their distributions fit Hardy-Weinberg expectations. Genotypes for each SNP were statistically consistent with Hardy-Weinberg expectations. Odds ratios and 95% confidence intervals (CI) were used to estimate risks. These measures were calculated using SAS software (version 9.1). Information on maternal race/ethnicity was obtained for case and control infants from California birth certificates. Logistic regression was used to compute risk estimates adjusted for maternal race/ethnicity (white Hispanic; white nonHispanic, and other). Analyses estimated defect risks (spina bifida or conotruncal heart defects) for each SNP assuming a recessive model, i.e., homozygous variant genotype compared to homozygous reference genotype and heterozygous variant genotype compared to homozygous reference genotype. In addition to single SNP-at-a-time analyses, we explored haplotype block analyses. Haplotype analyses were performed using Haploview version 3.32. Identified blocks were assessed with odds ratios.

We examined risks for each of the 118 SNPs and for each of the two birth defect outcome (Additional file 1). Few odds ratios (ORs) revealed sizable departures from 1.0. Given the large number of comparisons (n = 472) we expected more ORs to be substantially different from 1.0 by chance. With respect to spina bifida, we observed ORs with confidence intervals that did not include 1.0 for the following SNPs (heterozygous or homozygous) relative to the reference genotype: BHMT (rs3733890) OR = 1.8 (1.1–3.1), CBS (rs2851391) OR = 2.0 (1.2–3.1); CBS (rs234713) OR = 2.9 (1.3–6.7); MTHFD1 (rs2236224) OR = 1.7 (1.1–2.7); MTHFD1 (hcv11462908) OR = 0.2 (0–0.9); MTHFD2 (rs702465) OR = 0.6 (0.4–0.9); MTHFD2 (rs7571842) OR = 0.6 (0.4–0.9); MTHFR (rs1801133) OR = 2.0 (1.2–3.1); MTRR (rs162036) OR = 3.0 (1.5–5.9); MTRR (rs10380) OR = 3.4 (1.6–7.1); MTRR (rs1801394) OR = 0.7 (0.5–0.9); MTRR (rs9332) OR = 2.7 (1.3–5.3); TYMS (rs2847149) OR = 2.2 (1.4–3.5); TYMS (rs1001761) OR = 2.4 (1.5–3.8); and TYMS (rs502396) OR = 2.1 (1.3–3.3). Each gene involving multiple SNP associations was investigated for linkage disequilibrium. Modest to strong evidence for linkage disequilibrium was observed for SNPs in each gene, i.e., D' ranged from 0.44 to 1.0 with all p values < 10^-4. With respect to conotruncal heart defects, we did not observe any OR with a confidence interval that did not include 1.0.

We did not observe evidence to indicate that risk patterns were confounded by race/ethnicity groupings, i.e., observed ORs were not substantially altered after adjusting for maternal race/ethnicity (not shown, available from authors upon request).

Haplotype analyses were stratified by race/ethnic background (Hispanic white and nonHispanic white). We observed evidence of a nonrandom haplotype association with TYMS for spina bifida and conotruncal heart defects among nonHispanic whites. Lack of evidence for other haplotypes that were observed overall was likely the result of smaller sample sizes from stratification.

In this California population we found only modest evidence that polymorphisms in 14 folate-related genes contributed to risk of spina bifida. SNPs contributing risks were in BHMT, CBS, MTHFD1, MTHFD2, MTHFR, MTRR, and TYMS. Haplotype association analyses further identified TYMS and MTHFR as potential contributors to spina bifida risk. In general, however, most of these folate-related genes showed little evidence for a gene-only effect on risk of spina bifida, and even less, on risks of conotruncal heart defects.

The 14 genes studied here have been implicated in the complex metabolic cycle involving folate (e.g., [25–27]). To our knowledge, this study contained the largest number of SNPs in folate-related genes interrogated as risk factors for human spina bifida or conotruncal heart defects. Previous studies have included some of the SNPs examined here. For example, Boyles and colleagues [28] studied 28 SNPs in 11 folate-related genes and found that only BHMT (rs3733890) was associated with increased spina bifida risk. This BHMT association is consistent with our findings that showed an odds ratio of 1.8 (1.1–3.1).

Many studies have explored MTHFR 677 (rs1801133) polymorphism. A range of risks, including no-effect, has been reported for this SNP relative to spina bifida. Botto and Yang [15] in a meta-analysis demonstrated a pooled odds ratio of 1.8 for spina bifida among infants homozygous for 677T. A few studies have also explored this 677 SNP in MTHFR as a risk factor for selected congenital heart defects, with most investigations finding no or little association [18, 19, 29–31]. We did observe a 2-fold increased risk of spina bifida associated with this SNP for homozygous infants. Further, haplotype analyses showed some association for the MTHFR gene as well.

Methionine synthase (MTR) is a vitamin B₁₂ dependent enzyme that is essential for the remethylation of homocysteine to methionine. The enzyme is required by cells for the essential accumulation of folate [32]. One particular SNP (A2756G; rs1805087) has been considerably investigated, with increased risks of NTDs reported in some studies [33–35], but not in others [36, 37]. We did not find an increased risk for spina bifida or conotruncal heart defects associated with this SNP or any other SNP of MTR.

Cystathione beta synthase (CBS) is critical to the degradation of homocysteine to cysteine. Regulation of this pyridoxal phosphate-dependent enzyme catalyzes the hydroxyl group of serine with the thiolate of homocysteine [38]. The polymorphism in the CBS gene that has received the most study is a 68 bp insertion (844ins68), with predominantly no associations observed for NTDs [27]. This polymorphism was not investigated in the current study. We did observe, however, two CBS SNPs (rs2851391 and rs234713) that showed increased risks for spina bifida. Boyles et al [28], albeit using a different study design than ours, observed that these two SNPs were not differentially transmitted from parents of infants with spina bifida.

MTRR gene polymorphisms (particularly rs1801394) have been investigated as a risk factor for both spina bifida and congenital heart defects. Polymorphisms in MTRR could alter homocysteine levels because methionine synthase reductase participates in maintaining the vitamin B₁₂-dependent conversion of homocysteine to methionine [32]. The most frequently studied MTRR polymorphism has been the 66A>G (rs1801394). This polymorphism in infants was associated with a 2.6-fold increased risk of spina bifida in an earlier study by us [33], it was associated with increased risk for spina bifida in another study only when vitamin B₁₂ levels were low [39], or in combination with MTHFR CC genotype [35]. The polymorphism in mothers of infants with neural tube defects has been associated with increased risk in one study [40], but not in another study [41]. Recent work from the Netherlands has shown a lack of association between this polymorphism and risk for conotruncal heart defects [42] as well as no increased risks for a broader phenotypic group of heart defects [43]. In this study, the 66A>G polymorphism was not associated with increased risks for either spina bifida or conotruncal heart defects. We did observe, however, approximately 3-fold elevated risks for spina bifida associated with three other MTRR SNPs (rs162036, rs10380, and rs9332). The significance of these observations will have to be explored in future studies.

With respect to MTHFD1 and MTHFD2, two studies have demonstrated an association with one polymorphism (rs 2236225) in MTHFD1 and NTD risk. One study showed a 1.5-fold increase in risk of an NTD-affected pregnancy in Irish women who were homozygous AA [44], a finding that confirmed an earlier increased risk that was identified in Irish women. Another study showed a similar risk for Italian women as well as a 1.9-fold risk for infants with the AA genotype to have spina bifida [45]. For this particular SNP, we observed a similar magnitude of risk (OR = 1.6) for infants with the homozygous genotype, but the estimate was relatively imprecise. We did observe a modestly elevated spina bifida risk for individuals who were homozygous for another MTHFD1 SNP (rs2236224) and modestly lowered risks for three others (hcv11462908, rs702465, and rs7571842). These observations will need to be replicated in future studies.

Polymorphisms in the DHFR gene have not been well-studied for their role in risks of birth defects. Three studies have investigated a 19-bp deletion with mixed results [46–48]. That particular polymorphism was not interrogated in the current study.

Our analyses did not show associations with SNPs in RFC1. Previous investigations of this gene have focused on a particular SNP, rs1051266, and have found mixed results [37, 41, 49–53]. This particular SNP was not analyzed here as a result of too many samples failing to be genotyped for this SNP using the SNPlex platform.

Recent studies have focused on the importance of TYMS in the folate metabolic pathway, including associations between TYMS polymorphisms and folate levels [54–56]. This folate-dependent enzyme catalyzes the reductive methylation of deoxyuridylate (dUMP) to thymidylate (dTMP), thereby playing a central role in DNA synthesis and repair by serving as the primary intracellular source of dTMP [54, 57–59]. We previously [56] observed a 4-fold increased risk of spina bifida in nonHispanic white infants who had a polymorphism for a 28 bp insertion in the promoter region. This observation, however, was not replicated in a population from the northern UK [55]. This particular polymorphism was not interrogated in the current study. Three of the five TYMS SNPs (rs284179, rs1001761, and rs502396) investigated here showed elevated risks for spina bifida for both heterozygote or homozygote individuals. This finding and the corresponding haplotype finding (Table 3) will be important to explore in future studies.

The strengths of this study were: 1) it investigated the potential effects of a large number of folate pathway SNPs, as well as investigated haplotype associations; 2) it had population-based ascertainment of two case phenotypes and controls; and 3) it included cases and controls born before the US food supply was fortified with folic acid, thus we would expect a sizable proportion of cases to have been folate-responsive.

Conversely, our study was limited in its effect estimation owing to small sample sizes for some comparisons. For example, our study had 80% power to detect risks of 2.5 or more associated with genotypes that were observed in at least 4% of controls. Another potential limitation is the lack of information on maternal folate status. Our working hypothesis is that transient elevation in maternal serum folate from supplementation or dietary intake could prevent birth defects by overcoming metabolic inefficiencies or transport-related issues. Absence of information on low folate status would make it more difficult to find putative genotypes. It is also possible that the protective effect of folic acid relates to correction of a maternal metabolic defect, rather than the fetus. Our study was limited to infant genotype information. Thus, we were unable to investigate the potential effects of maternal genotype. As with any study that seeks to explore associations with a large number of genotypes, findings are subject to chance owing to multiple comparisons. As noted above, we conducted 472 analytic comparisons and thus expected more "statistically significant" findings to arise by chance alone. Further, our findings may have been influenced by uncontrolled confounding by population stratification undetectable in analyses stratified or adjusted by race/ethnicity [60, 61]. Lastly, the selected SNPs represent only a fraction of the potential variation of the studied genes. Thus, full gene coverage was not achieved even though a large number of SNPs was studied.

Despite compelling evidence that folate intake by women in early pregnancy substantially reduces risks of selected birth defects, the underlying mechanisms have not been elucidated. Our study attempted to determine genetic mechanisms responsible for folic acid's preventive effects. Our observations do not implicate a particular folate transport or metabolism gene to be strongly associated with risks for spina bifida or conotruncal defects. Although we explored a sizable number of polymorphic areas in these genes, we clearly did not capture all the genetic variation. Thus, these genes may continue to be candidates for further inquiry. Alternatively, the preventive role of folate may be via other biological mechanisms such as methylation of nonfolate-related genes that participate in the closure of the neural tube or the development of the heart.