Skip to main content
  • Research article
  • Open access
  • Published:

Replication and exploratory analysis of 24 candidate risk polymorphisms for neural tube defects



Neural tube defects (NTDs), which are among the most common congenital malformations, are influenced by environmental and genetic factors. Low maternal folate is the strongest known contributing factor, making variants in genes in the folate metabolic pathway attractive candidates for NTD risk. Multiple studies have identified nominally significant allelic associations with NTDs. We tested whether associations detected in a large Irish cohort could be replicated in an independent population.


Replication tests of 24 nominally significant NTD associations were performed in racially/ethnically matched populations. Family-based tests of fifteen nominally significant single nucleotide polymorphisms (SNPs) were repeated in a cohort of NTD trios (530 cases and their parents) from the United Kingdom, and case-control tests of nine nominally significant SNPs were repeated in a cohort (190 cases, 941 controls) from New York State (NYS). Secondary hypotheses involved evaluating the latter set of nine SNPs for NTD association using alternate case-control models and NTD groupings in white, African American and Hispanic cohorts from NYS.


Of the 24 SNPs tested for replication, ADA rs452159 and MTR rs10925260 were significantly associated with isolated NTDs. Of the secondary tests performed, ARID1A rs11247593 was associated with NTDs in whites, and ALDH1A2 rs7169289 was associated with isolated NTDs in African Americans.


We report a number of associations between SNP genotypes and neural tube defects. These associations were nominally significant before correction for multiple hypothesis testing. These corrections are highly conservative for association studies of untested hypotheses, and may be too conservative for replication studies. We therefore believe the true effect of these four nominally significant SNPs on NTD risk will be more definitively determined by further study in other populations, and eventual meta-analysis.

Peer Review reports


Neural tube defects (NTDs) are common, severe birth defects that affect the development of the spinal cord and/or brain. NTDs mainly present as spina bifida, anencephaly and encephalocele. Although the NTD rate has been reduced in areas with folic acid fortification [1]-[4], reported NTD rates include ~1 in 1000 births in the US [5],[6], 0.1-1.1 in 1000 births in Europe [7], 1–4.8 in 1000 pregnancies in China [3], 3.9-8.8 in 1000 births in India [8], and 0.2-9.6 in 1000 births in Latin America [9]. Both environmental and genetic factors are involved in the development of NTDs. In western countries, dietary folate is the strongest known environmental factor. Mothers with an NTD pregnancy have significantly lower levels of folate [10], and periconceptional folate supplementation can lower NTD risk by 50-70% [3],[11],[12]. The most firmly established genetic factor is a common polymorphism (rs1801133, c.677C > T) in the 5, 10-methylene-tetrahydrofolate reductase (MTHFR) gene. The alternate nucleotide in this single nucleotide polymorphism (SNP) changes an alanine to a valine (p.A222V), resulting in a thermolabile form of the protein [13]. The MTHFR c.677 TT genotype has been significantly associated with NTD cases in many populations (reviewed in [14],[15]). Other NTD risk SNPs have been identified in a number of populations (e.g., rs2236225 (c.1958G > A, p.R653Q) in methylenetetrahydrofolate dehydrogenase (NADP + dependent) 1, methenyltetrahydrofolate cyclohydrolase, formyltetrahydrofolate synthetase (MTHFD1) [16]-[21], rs1805087 (c.2756A > G, p.D919G) in methionine synthase (MTR) [15],[16],[22]-[29], rs1051266 (c.80A > G, p.H27R) in the reduced folate carrier (RFC1/SLC19A1) [15],[21],[30]-[33]), and rs70991108 (19 bp indel; c.86 + 59_86 + 60insACCTGGGCGGGACGCGCCA) in dihydrofolate reductase (DHFR) [34]-[36]). The majority of associations have not been subjected to robust replication.

Rigorously establishing genetic risk for any multifactorial disorder is important but inherently difficult. Genetic risk is optimally detected with a large number of racially/ethnically homogenous samples that include cases with a well-defined, highly penetrant phenotype. These are difficult criteria to meet for rare, complex disorders, meaning that statistical power is often compromised to some extent. This, in combination with the bias that significant associations are more likely to be published [37]-[40], contributes to many associations in the literature that would not survive correction for multiple tests.

Specifically in the field of genetic risk for NTDs, studies including large numbers of candidate polymorphisms and/or genes have implicated many candidate SNPs exhibiting nominal NTD associations [21],[41]-[45]. Each of the cited studies examined a minimum of 37 SNPs for NTD risk, yet only the association of rs1907362 in cubilin (CUBN) with NTDs survived correction in the original study [21]. The most likely explanation is that the reported nominally significant associations were due to chance. However, the large number of tests performed, reduced power due to limited availability of samples, and/or the population-specific effect of a risk allele can all contribute to Type II errors. Therefore replication studies are essential for the validation of any genetic risk factor for NTDs.

In one of the largest studies to date we previously reported evaluation of the common genetic variation in 82 candidate genes for NTD risk in an Irish population [45]. The inclusion of 570 NTD cases (~95% spina bifida cases), their parents and 999 controls made it possible to use both family based methods (the transmission disequilibrium test and log-linear analysis) and case–control based methods (logistic regression) to evaluate genetic risk for NTDs in cases (case effect) as well as genetic risk for an NTD pregnancy in mothers (maternal effect). In the current study, we test 24 of the resulting nominally significant association signals using 9 case–control tests and 14 family-based tests in racially/ethnically similar populations. The 9 SNPs that showed nominal significance by case–control analysis in the previously studied Irish NTD cohort were additionally evaluated for NTD risk using varying models in racially/ethnically similar and distinct populations.


Statistical methods

For replication analyses, each candidate risk allele was retested using the same association test and model for which it was nominally significant in the primary study. These tests are noted in Tables 1 and 2, and include: 1) case–control tests of logistic regression using a continuous, recessive or dominant genetic model; 2) log-linear tests of case or maternal effect using a dominant or recessive genetic model; and 3) the Spielman transmission disequilibrium test (TDT [46]).

Table 1 case–control replication analyses in the New York State NTD cohort
Table 2 Family-based replication analyses in the combined United Kingdom NTD cohort

Secondary hypotheses involved more broadly testing the candidate SNPs for NTD risk in other populations and/or NTD classifications. Three models (continuous, recessive and dominant) of logistic regression were performed to generate genotype relative risks (GRRs) and 95% confidence intervals (CI).

Study populations

Family-based replication was performed in a United Kingdom cohort consisting of trios of NTD cases and their parents. Exclusion of the “Other NTDs” class, which may include a small number of families with multiple defects or spina bifida occulta, did not substantially change the results. This cohort consists of three groups (Table 3). First, the 400 case families recruited with the assistance of the UK Association for Spina Bifida and Hydrocephalus (ASBAH) (England and Wales) are fully described elsewhere [47]. Second, 131 case families were recruited from Northern Ireland. Ethical approval of the use of these samples was granted by the UK Multi-Centre Research Ethics Committee (University of Newcastle, UK, and University of Northern Ireland, Coleraine) and the Institutional Review Board at the National Human Genome Research Institute (Bethesda, MD, USA).

Table 3 Number of NTD family members in the UK family-based cohort

Replication of findings from case–control analyses was performed in a NYS case–control cohort. NTD cases from the entire state born 1998–2005 were identified by their inclusion in the NYS Congenital Malformations Registry (NYSCMR). Matched controls were selected as a random sample of non-malformed control infants born 1998–2005 from NYS Newborn Screening Program records. Demographic data was obtained by matching to NYS birth certificates. Four controls matched for maternal race/ethnicity were selected for each NTD case. Due to low numbers, the Asian and “other” racial/ethnic subgroups were excluded, and subjects coded as “Hispanic, white” and “Hispanic, other” were combined into a single Hispanic group for analysis. Case diagnoses allowed classification of the NTDs into six subgroups based on NTD subtype (anencephaly, encephalocele, spina bifida) and whether the NTD was isolated or part of multiple defects (Table 4). To perform the replication analyses, only isolated NTDs of all three subtypes among the NHWs were considered in order to most closely match the composition of the NTD cases in the original Irish cohort. Secondary hypotheses involved analyses of: 1) all NTD cases and all controls, stratified by race/ethnicity; and 2) isolated spina bifida cases and corresponding matched controls. Approval for the use of the de-identified samples was granted by the New York State Department of Health Institutional Review Board and the NIH Office for Human Research Protections.

Table 4 Number of NTD cases and controls by race/ethnicity in the NY state cohort


MTHFD1 rs2236225 was genotyped in the UK cohort as a restriction fragment length polymorphism (PCR-RFLP) using MspI as previously described [17],[48]. Genotypes were obtained for 91.6% of NTD fathers, 92.5% of NTD cases and 93.6% of NTD mothers. Concordance was 100% for repeated (n = 83) and for re-plated (the testing of a second sample from the DNA source) (n = 92) samples using an independent assay based on detection of allele-specific primer extension using matrix-assisted laser desorption/ionization – time of flight (MALDI-TOF) mass spectrometry (Sequenom, San Diego, CA, USA). The remaining SNPs in Table 2 were also genotyped in the UK cohort using the Sequenom platform. Two independent assays of folate hydrolase (FOLH1) rs16906205 failed, so the proxy FOLH1 rs11040291 (r2 = 1 in Hapmap CEU) was genotyped and reported instead. Two independent assays to genotype methionine adenosyltransferase II, beta (MAT2B) rs17535909, were attempted but both failed, and there was no proxy for this singleton SNP. For this set of 13 SNPs, the average call rates were ≥96.8% for each family group (NTD cases, NTD mothers or NTD fathers). Re-plated and re-genotyped samples covered >18% of the cohort with 99.2% genotype concordance for this set of 13 SNPs. The 14 SNPs typed in the UK cohort exhibited non-Mendelian inheritance in <1% of families. These SNPs were also in Hardy Weinberg Equilibrium (HWE, p > 0.01) for each family group. Genotypes for families exhibiting non-Mendelian inheritance and other discordant genotypes were excluded from analysis.

For NYS samples, DNA was extracted from one 3-mm archived dried blood spot specimen [49] and whole-genome amplified using a primer extension pre-amplification method, as described previously [50]. SNPs were genotyped by KBiosciences (Herts, UK) using KASPar chemistry. Eight SNPs were genotyped in duplicate using independently whole-genome amplified DNA aliquots with 100% concordance in genotype calls. FOLH1 rs383028 was genotyped using genomic DNA because data from amplified DNA did not pass quality control criteria. The average call rate for 9 SNPs (Table 1) was 99.9% for both cases and controls. Replated samples covered 6.5% of the cohort with genotype concordance of 100%. No SNPs deviated from HWE (p > 0.01) in any case or control group for any race/ethnicity.


The primary aim of this study was to perform replication analyses of the nominally significant NTD associations identified in a recent study in an Irish population [45]. The secondary aim was to test a subset of these candidate SNPs for association using alternate risk models and populations.

Replication analyses

Replication criteria

The replication strategy was designed to retest nominally significant NTD-associated SNPs in racially/ethnically matched populations using the same association tests and genetic models that previously yielded the lowest p-values among 1441 SNPs in 82 candidate genes tested in an Irish population [45]. case–control association tests were performed in a cohort of 190 isolated NTD cases and 941 controls from non-Hispanic white (NHW) mothers from NYS, and family-based tests were performed in NTD trios (n = 530) consisting of NTD cases and their parents from the United Kingdom, including centers in Northern Ireland, England and Wales.

The top 25 groups of SNPs sharing high linkage disequilibrium (LD; D’ > 0.9) with the lowest p-values for any test were selected for replication (52 SNPs total). Without access to NTD mothers and corresponding controls several loci could not be retested. This included 9 independent mother-control signals (17 SNPs) in adenosine deaminase (ADA), alcohol dehydrogenase 1 family, member A2 (ALDH1A2), catechol-O-methyltransferase (COMT), CUBN, MTHFD1, methylenetetrahydrofolate dehydrogenase (NADP + dependent) 1-like (MTHFD1L), brachyury (T), and transcription factor AP-2 alpha (TFAP2A). Of the remaining 35 SNPs, the SNP with the lowest observed p-value was selected for replication testing whenever SNPs with a minimum p-value for the same test and model shared high LD (D’ > 0.9). This reduced the number of tests to 24 (9 case–control tests and 15 family-based tests). Results of these tests are shown in Tables 1 and 2.

Replication of associations detected in case–control analyses

We used the NTD case samples from NYS to replicate case effects previously observed in the Irish NTD cohort. Each SNP was tested for association with NTDs by logistic regression with the same genetic model used for the original observation. No SNP was observed to be significantly associated with isolated NTDs in white cases from NYS (Table 1). This lack of replication was accompanied by the corresponding genotype relative risk (GRR) values indicating an inconsistent effect of the candidate risk allele to that observed in the original study for eight of the nine SNPs.

Replication of associations detected in family-based analyses

Replication analyses of previously observed case or maternal effects detected by log-linear analysis were repeated with the same genetic model in a combined cohort of NTD triads from the United Kingdom. Two SNPs showed nominally significant association with NTDs: ADA rs452159 in a dominant model (GRR = 1.86 [1.01-3.40], p = 0.045), and MTR rs10925260 in a recessive model (GRR = 0.71 [0.53-0.94], p = 0.019). In both cases, the magnitude of the effect and the specific risk allele was similar to the original observation in the Irish NTD cohort (Table 2). If applied, the significance of these results would not withstand Bonferroni correction for multiple tests when considering all SNPs tested in the current study. In contrast to the case–control analyses, the direction of effect for the GRRs of these family-based associations largely agreed between the initial and replication studies (10 of 14), regardless of significance.

Secondary hypotheses - exploratory analyses in new populations

Applying analyses that yielded nominal associations from the initial study

In addition to replication, these candidate SNPs were tested for association in other NTD populations using various models. The nine SNPs selected for case–control replication in the NYS cohort were first tested in African American and Hispanic cases with an isolated NTD and controls using the same tests and models for which each SNP had been originally observed to be nominally associated in the Irish population (Table 1, 2 tests/SNP). Of the nine SNPs examined in each of the two racial/ethnic groups, only one, ALDH1A2 rs7169289, was found to be nominally associated with isolated NTDs -- in African Americans in a continuous model (GRR = 0.57 [0.34-0.98], p = 0.041). This “protective” effect is in contrast to the risk effect seen for this SNP in the Irish cohort using a dominant model.

case–control analyses in all NTD cases vs. isolated Spina Bifida cases

These SNPs were also tested by performing logistic regression using three genetic models (continuous, dominant and recessive) in all NTD cases and controls in each of the three racial/ethnic groups (9 tests/SNP). Of these nine SNPs tested in nine models (N = 81 tests), AT rich interactive domain 1A (ARID1A) rs11247593 was the only SNP found to be nominally associated with NTDs -- in non-Hispanic whites in a dominant model (GRR = 0.58 [0.35-0.97], p = 0.037). This protective effect is in contrast to the risk effect (GRR > 1) seen in the Irish cohort using a continuous model (Table 1).

Lastly, the same tests were performed in a restricted subset of isolated spina bifida cases and controls in the three racial/ethnic groups (9 tests/SNP). These 81 tests generated two significant findings: ALDH1A2 rs7169289 was nominally associated with isolated spina bifida cases in the African American population in continuous (GRR = 0.46 [0.25-0.87], p = 0.017) and dominant (GRR = 0.47 [0.24-0.94], p = 0.033) models. Similarly, a significant risk effect was observed in the original Irish population, but not the NYS NHW population (Table 1).

These results are summarized in Table 5. These results are nominally significant, though would not withstand Bonferroni correction using the total number of tests as the correction factor.

Table 5 Nominally significant associations in the New York state cohort - secondary hypotheses


Our study addressed two questions. First, are SNPs that were nominally significant upon testing for NTD association in an Irish population also associated with NTDs in a similar, independent population? Second, are the SNPs that were nominally significant upon testing for NTD association in an Irish population by logistic regression also associated with NTDs in different racial/ethnic populations when using a broader range of association tests, genetic models and NTD case groupings?

Our criteria for replication were stringent, and only performed for a SNP when the same association test and genetic model could be applied. Because we did not have samples from mothers of NYS NTD cases, we were unable to test for a maternal risk for nine independent mother-control signals in seven genes. Nine SNPs were re-tested for a case effect by logistic regression in a white NTD sample from NYS, and log-linear analyses were used to re-test ten SNPs for a case effect and five SNPs for a maternal effect in a UK NTD sample. Of 14 SNPs previously observed to be associated with NTDs by family-based analysis in an Irish population, two showed nominally significant NTD association in trios from the UK (Table 2). ADA rs452159 falls in the first intron, and is at the border of a D’ block encompassing the first exon and intron of the gene. MTR rs10925260 is in intron 23 of the gene, and is part of a large block of D’ LD encompassing the entire gene. Due to their strong D’ linkage, these nominally significant associations may reflect a direct signal or that of a causative SNP linked to the tested SNP(s).

Although none of the 9 SNPs found to be nominally significant in Irish NTDs were replicated in isolated NYS white NTD cases, two were nominally associated with NTDs under different conditions. ARID1A rs11247593 did not replicate by logistic regression using a continuous model in isolated NTDs (Table 1), but was significant when all white NTD cases and controls were tested using a dominant model (Table 5). ARID1A rs11247593 is intronic and is part of a large D’ block extending over the entire gene, so the signal may be due to the tested SNP or a linked causative SNP. In addition, ALDH1A2 rs7169289 was significantly associated with isolated NTDs and isolated spina bifida cases in the African American population (Table 5). These results most likely represent a single association signal as the genetic models (continuous, dominant) and NTD sets (isolated NTDs, isolated SB) for these analyses overlap. The estimated effects observed in African American NTD cases from NYS (Table 5; OR = 0.46-0.57, p = 0.017-0.041) are similar to the original effect observed in an Irish population (Table 1; OR = 0.67 [0.52-0.86], p = 0.0016). ALDH1A2 rs7169289 is just downstream of the gene. It falls between two blocks of D’ LD, and the association signal may be due to this SNP or a linked SNP in either block.

As was found in the original study, none of the observed associations in the NYS and UK sample sets would have survived correction for multiple tests. One interpretation is that these associations are indeed due to chance and that the tested SNPs do not contribute to NTD risk. It is not clear, however, that Bonferroni correction is appropriate for replication studies. All the tests performed in this study were clearly based on individual, a priori hypotheses generated for our initial study and supported by previous data. Additionally, factors that can contribute to Type II error should be considered, such as population stratification and genetic heterogeneity.

The most important factor may be limited sample size, which compromises the power to detect true associations. Compared to the number of NTD cases (n ~ 570) in the original Irish cohort, a limited number of white NYS isolated NTD cases (n ~ 190) were available in the current study. With ~530 NTD trios, however, the UK NTD cohort appears comparable, but this is before taking into consideration the power required to replicate an association. Due to the “winner’s curse,” or the tendency of an initial observation to overestimate the effect size or significance of a true association, replication studies generally require larger study sample numbers to detect the original effect [51],[52]. This may explain why the well known NTD risk allele of MTHFR rs1801133 (T) is not associated with NTDs by TDT (p = 0.742) in our large UK NTD cohort. In fact, it is nominally significant as a protective factor in a recessive model by log-linear analysis (GRR = 0.663 [0.441-0.997], p = 0.049). Considering that this association has been replicated in many studies in other populations supporting its role in NTD risk, it would be surprising if the MTHFR rs1801133 TT genotype does not contribute to NTDs in the UK population. A lack of replication power and chance seem the most likely explanations for this failure to replicate, which must be taken into consideration for the other SNPs tested for replication in this study. Authentic validation of any genetic association requires accumulation of evidence over time, involving multiple studies in independent populations.

Lack of adequate power is a pervasive problem in the field of NTD genetics. One recent review estimated that ~1000 cases would be required to attain 80% power to detect an odds ratio of 2 or under, yet approximately one quarter of published NTD association studies used fewer than 100 cases [53]. This problem is compounded when attempting to evaluate candidate risk SNPs in less studied populations. While we were able to examine candidate risk SNPs of interest in African American and Hispanic populations from NYS, far fewer cases were available compared to whites (Table 4). Although some population-based studies have included African American cases in aggregate analyses of genetic risk for NTDs [54]-[56], and there are studies of NTD risk in African cohorts [57], this is the first report of candidate NTD SNPs evaluated in an African American case–control cohort. The nominal association observed for ALDH1A2 rs7169289 may be real, but needs confirmation.

Identifying potential NTD risk SNPs has proven much easier than validating them. A single replication study is not definitive, and the limited numbers in existing NTD cohorts may contribute to subsequent underpowered studies failing to confirm reported associations. By combining data from multiple published studies, meta-analyses increase power and confidence in whether a SNP truly contributes to NTD risk. Meta-analyses involving aggregate data from several hundred to several thousand NTD cases and controls have confirmed MTHFR rs1801133 (c.677C > T) as a maternal [58] and case risk factor for NTDs [14],[15],[58],[59]. Although comparatively fewer data were available, meta-analyses show 5-methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR) rs1801394 (c/66A > G) is a maternal risk factor for NTDs [60],[61], while MTR rs1805087 (c/2756A > G) does not contribute to maternal or case risk [15],[22],[60],[62]. However, meta-analysis requires a large amount of data including full genotype information. These data are absent in publications for the majority of NTD risk SNPs that show nominal significance. It is therefore essential to publish all association study data in a format that allows future meta-analyses to be performed (i.e., genotype counts or a way to unambiguously determine them). As such, we have reported the genotype data for all 24 SNPs tested in our study (Additional file 1: Table S1 and Additional file 2: Table S2).


In summary, this study represents both the effort to validate nominally significant NTD associations via replication, and to test a subset of these SNPs in alternate NTD groupings and racial/ethnic populations. The NTD associations we observed with four SNPs would not have survived correction for multiple tests. However, we believe it would be overly severe to apply Bonferroni correction to ADA rs452159 and MTR rs10925260, which were nominally significant by strict replication analysis, as well as ARID1A rs11247593, which was found to be associated by a different model and expanded definition of NTDs. Additionally, to our knowledge, this is the first report of a genetic association with NTDs in African Americans and so the nominal significance of ALDH1A2 rs7169289 in the small NYS NTD cohort warrants follow-up. Further testing in independent populations and meta-analyses are needed to clarify the role of these and a multitude of nominally significant SNPs associated with NTDs.

Additional files



Linkage disequilibrium


Neural tube defects


Single nucleotide polymorphism


Transmission disequilibrium test


New York State


Non-Hispanic white


Genotype relative risk


  1. Canfield MA, Honein MA, Yuskiv N, Xing J, Mai CT, Collins JS, Devine O, Petrini J, Ramadhani TA, Hobbs CA, Kirby RS: National estimates and race/ethnic-specific variation of selected birth defects in the United States, 1999−2001. Birth Defects Res A Clin Mol Teratol. 2006, 76: 747-756. 10.1002/bdra.20294.

    Article  CAS  PubMed  Google Scholar 

  2. Botto LD, Lisi A, Bower C, Canfield MA, Dattani N, De Vigan C, De Walle H, Erickson DJ, Halliday J, Irgens LM, Lowry RB, McDonnell R, Metneki J, Poetzsch S, Ritvanen A, Robert-Gnansia E, Siffel C, Stoll C, Mastroiacovo P: Trends of selected malformations in relation to folic acid recommendations and fortification: an international assessment. Birth Defects Res A Clin Mol Teratol. 2006, 76: 693-705. 10.1002/bdra.20307.

    Article  CAS  PubMed  Google Scholar 

  3. Berry RJ, Li Z, Erickson JD, Li S, Moore CA, Wang H, Mulinare J, Zhao P, Wong LY, Gindler J, Hong SX, Correa A: Prevention of neural-tube defects with folic acid in China. China-U.S. Collaborative Project for Neural Tube Defect Prevention. N Engl J Med. 1999, 341: 1485-1490. 10.1056/NEJM199911113412001.

    Article  CAS  PubMed  Google Scholar 

  4. Blencowe H, Cousens S, Modell B, Lawn J: Folic acid to reduce neonatal mortality from neural tube disorders. Int J Epidemiol. 2010, 39 (Suppl 1): i110-i121. 10.1093/ije/dyq028.

    Article  PubMed  PubMed Central  Google Scholar 

  5. A Report from the International Clearinghouse for Birth Defects Monitoring Systems. 1991, Elsevier Science Publishers, Amsterdam

  6. Busby A, Abramsky L, Dolk H, Armstrong B: Preventing neural tube defects in Europe: population based study. BMJ. 2005, 330: 574-575. 10.1136/bmj.330.7491.574.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Botto LD, Lisi A, Robert-Gnansia E, Erickson JD, Vollset SE, Mastroiacovo P, Botting B, Cocchi G, de Vigan C, de Walle H, Feijoo M, Irgens LM, McDonnell B, Merlob P, Ritvanen A, Scarano G, Siffel C, Metneki J, Stoll C, Smithells R, Goujard J: International retrospective cohort study of neural tube defects in relation to folic acid recommendations: are the recommendations working?BMJ 2005, 330:571.,

  8. Cherian A, Seena S, Bullock RK, Antony AC: Incidence of neural tube defects in the least-developed area of India: a population-based study. Lancet. 2005, 366: 930-931. 10.1016/S0140-6736(05)67319-9.

    Article  PubMed  Google Scholar 

  9. Rosenthal J, Casas J, Taren D, Alverson CJ, Flores A, Frias J: Neural tube defects in Latin America and the impact of fortification: a literature review. Public Health Nutr. 2014, 17 (3): 537-550. 10.1017/S1368980013000256.

    Article  PubMed  Google Scholar 

  10. Kirke PN, Molloy AM, Daly LE, Burke H, Weir DG, Scott JM: Maternal plasma folate and vitamin B12 are independent risk factors for neural tube defects. Q J Med. 1993, 86: 703-708.

    CAS  PubMed  Google Scholar 

  11. Prevention of neural tube defects: results of the Medical Research Council Vitamin Study. Lancet. 1991, 338: 131-137. 10.1016/0140-6736(91)90133-A.

  12. Czeizel AE, Dudas I: Prevention of the first occurrence of neural-tube defects by periconceptional vitamin supplementation. N Engl J Med. 1992, 327: 1832-1835. 10.1056/NEJM199212243272602.

    Article  CAS  PubMed  Google Scholar 

  13. Frosst P, Blom HJ, Milos R, Goyette P, Sheppard CA, Matthews RG, Boers GJ, den Heijer M, Kluijtmans LA, van den Heuvel LP: A candidate genetic risk factor for vascular disease: a common mutation in methylenetetrahydrofolate reductase. Nat Genet. 1995, 10: 111-113. 10.1038/ng0595-111.

    Article  CAS  PubMed  Google Scholar 

  14. Botto LD, Yang Q: 5,10-Methylenetetrahydrofolate reductase gene variants and congenital anomalies: a HuGE review. Am J Epidemiol. 2000, 151: 862-877. 10.1093/oxfordjournals.aje.a010290.

    Article  CAS  PubMed  Google Scholar 

  15. Zhang T, Lou J, Zhong R, Wu J, Zou L, Sun Y, Lu X, Liu L, Miao X, Xiong G: Genetic variants in the folate pathway and the risk of neural tube defects: a meta-analysis of the published literature.PLoS One 2013, 8:e59570.,

  16. Doudney K, Grinham J, Whittaker J, Lynch SA, Thompson D, Moore GE, Copp AJ, Greene ND, Stanier P: Evaluation of folate metabolism gene polymorphisms as risk factors for open and closed neural tube defects. Am J Med Genet A. 2009, 149A: 1585-1589. 10.1002/ajmg.a.32937.

    Article  CAS  PubMed  Google Scholar 

  17. Brody LC, Conley M, Cox C, Kirke PN, McKeever MP, Mills JL, Molloy AM, O’Leary VB, Parle-McDermott A, Scott JM, Swanson DA: A polymorphism, R653Q, in the trifunctional enzyme methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase/formyltetrahydrofolate synthetase is a maternal genetic risk factor for neural tube defects: report of the Birth Defects Research Group. Am J Hum Genet. 2002, 71: 1207-1215. 10.1086/344213.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. De Marco P, Merello E, Calevo MG, Mascelli S, Raso A, Cama A, Capra V: Evaluation of a methylenetetrahydrofolate-dehydrogenase 1958G > A polymorphism for neural tube defect risk. J Hum Genet. 2006, 51: 98-103. 10.1007/s10038-005-0329-6.

    Article  PubMed  Google Scholar 

  19. Parle-McDermott A, Kirke PN, Mills JL, Molloy AM, Cox C, O’Leary VB, Pangilinan F, Conley M, Cleary L, Brody LC, Scott JM: Confirmation of the R653Q polymorphism of the trifunctional C1-synthase enzyme as a maternal risk for neural tube defects in the Irish population. Eur J Hum Genet. 2006, 14: 768-772. 10.1038/sj.ejhg.5201603.

    Article  CAS  PubMed  Google Scholar 

  20. van der Linden IJ, Heil SG, Kouwenberg IC, den Heijer M, Blom HJ: The methylenetetrahydrofolate dehydrogenase (MTHFD1) 1958G > A variant is not associated with spina bifida risk in the Dutch population. Clin Genet. 2007, 72: 599-600. 10.1111/j.1399-0004.2007.00904.x.

    Article  CAS  PubMed  Google Scholar 

  21. Franke B, Vermeulen SH, Steegers-Theunissen RP, Coenen MJ, Schijvenaars MM, Scheffer H, den Heijer M, Blom HJ: An association study of 45 folate-related genes in spina bifida: Involvement of cubilin (CUBN) and tRNA aspartic acid methyltransferase 1 (TRDMT1). Birth Defects Res A Clin Mol Teratol. 2009, 85: 216-226. 10.1002/bdra.20556.

    Article  CAS  PubMed  Google Scholar 

  22. Yang M, Yang L, Qi L, Guo Y, Lin X, Zhang Y, Du Y: Association between the methionine synthase A2756G polymorphism and neural tube defect risk: a meta-analysis. Gene. 2013, 520: 7-13. 10.1016/j.gene.2013.02.005.

    Article  CAS  PubMed  Google Scholar 

  23. Morrison K, Papapetrou C, Hol FA, Mariman EC, Lynch SA, Burn J, Edwards YH: Susceptibility to spina bifida; an association study of five candidate genes. Ann Hum Genet. 1998, 62: 379-396. 10.1046/j.1469-1809.1998.6250379.x.

    Article  CAS  PubMed  Google Scholar 

  24. Christensen B, Arbour L, Tran P, Leclerc D, Sabbaghian N, Platt R, Gilfix BM, Rosenblatt DS, Gravel RA, Forbes P, Rozen R: Genetic polymorphisms in methylenetetrahydrofolate reductase and methionine synthase, folate levels in red blood cells, and risk of neural tube defects. Am J Med Genet. 1999, 84: 151-157. 10.1002/(SICI)1096-8628(19990521)84:2<151::AID-AJMG12>3.0.CO;2-T.

    Article  CAS  PubMed  Google Scholar 

  25. Akar N, Akar E, Deda G, Arsan S: Spina bifida and common mutations at the homocysteine metabolism pathway. Clin Genet. 2000, 57: 230-231. 10.1034/j.1399-0004.2000.570310.x.

    Article  CAS  PubMed  Google Scholar 

  26. Johanning GL, Tamura T, Johnston KE, Wenstrom KD: Comorbidity of 5,10-methylenetetrahydrofolate reductase and methionine synthase gene polymorphisms and risk for neural tube defects. J Med Genet. 2000, 37: 949-951. 10.1136/jmg.37.12.949.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. De Marco P, Calevo MG, Moroni A, Arata L, Merello E, Finnell RH, Zhu H, Andreussi L, Cama A, Capra V: Study of MTHFR and MS polymorphisms as risk factors for NTD in the Italian population. J Hum Genet. 2002, 47: 319-324. 10.1007/s100380200043.

    Article  CAS  PubMed  Google Scholar 

  28. Zhu H, Wicker NJ, Shaw GM, Lammer EJ, Hendricks K, Suarez L, Canfield M, Finnell RH: Homocysteine remethylation enzyme polymorphisms and increased risks for neural tube defects. Mol Genet Metab. 2003, 78: 216-221. 10.1016/S1096-7192(03)00008-8.

    Article  CAS  PubMed  Google Scholar 

  29. Brouns R, Ursem N, Lindemans J, Hop W, Pluijm S, Steegers E, Steegers-Theunissen R: Polymorphisms in genes related to folate and cobalamin metabolism and the associations with complex birth defects. Prenat Diagn. 2008, 28: 485-493. 10.1002/pd.2006.

    Article  CAS  PubMed  Google Scholar 

  30. De Marco P, Calevo MG, Moroni A, Merello E, Raso A, Finnell RH, Zhu H, Andreussi L, Cama A, Capra V: Reduced folate carrier polymorphism (80A– > G) and neural tube defects. Eur J Hum Genet. 2003, 11: 245-252. 10.1038/sj.ejhg.5200946.

    Article  CAS  PubMed  Google Scholar 

  31. Wang HG, Wang JL, Zhang J, Zhao LX, Zhai GX, Xiang YZ, Chang P: Reduced folate carrier A80G polymorphism and susceptibility to neural tube defects: a meta-analysis. Gene. 2012, 510: 180-184. 10.1016/j.gene.2012.02.020.

    Article  CAS  PubMed  Google Scholar 

  32. Relton CL, Wilding CS, Pearce MS, Laffling AJ, Jonas PA, Lynch SA, Tawn EJ, Burn J: Gene-gene interaction in folate-related genes and risk of neural tube defects in a UK population. J Med Genet. 2004, 41: 256-260. 10.1136/jmg.2003.010694.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Shaw GM, Lammer EJ, Zhu H, Baker MW, Neri E, Finnell RH: Maternal periconceptional vitamin use, genetic variation of infant reduced folate carrier (A80G), and risk of spina bifida. Am J Med Genet. 2002, 108: 1-6. 10.1002/ajmg.10195.

    Article  PubMed  Google Scholar 

  34. Johnson WG, Stenroos ES, Spychala JR, Chatkupt S, Ming SX, Buyske S: New 19 bp deletion polymorphism in intron-1 of dihydrofolate reductase (DHFR): a risk factor for spina bifida acting in mothers during pregnancy?. Am J Med Genet A. 2004, 124A: 339-345. 10.1002/ajmg.a.20505.

    Article  PubMed  Google Scholar 

  35. Parle-McDermott A, Pangilinan F, Mills JL, Kirke PN, Gibney ER, Troendle J, O’Leary VB, Molloy AM, Conley M, Scott JM, Brody LC: The 19-bp deletion polymorphism in intron-1 of dihydrofolate reductase (DHFR) may decrease rather than increase risk for spina bifida in the Irish population. Am J Med Genet A. 2007, 143A: 1174-1180. 10.1002/ajmg.a.31725.

    Article  CAS  PubMed  Google Scholar 

  36. van der Linden IJ, Nguyen U, Heil SG, Franke B, Vloet S, Gellekink H, den Heijer M, Blom HJ: Variation and expression of dihydrofolate reductase (DHFR) in relation to spina bifida. Mol Genet Metab. 2007, 91: 98-103. 10.1016/j.ymgme.2007.01.009.

    Article  CAS  PubMed  Google Scholar 

  37. Valachis A, Mauri D, Neophytou C, Polyzos NP, Tsali L, Garras A, Papanikolau EG: Translational medicine and reliability of single-nucleotide polymorphism studies: can we believe in SNP reports or not?. Int J Med Sci. 2011, 8: 492-500. 10.7150/ijms.8.492.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Munafo MR, Clark TG, Flint J: Assessing publication bias in genetic association studies: evidence from a recent meta-analysis. Psychiatry Res. 2004, 129: 39-44. 10.1016/j.psychres.2004.06.011.

    Article  PubMed  Google Scholar 

  39. Dwan K, Altman DG, Arnaiz JA, Bloom J, Chan AW, Cronin E, Decullier E, Easterbrook PJ, Von Elm E, Gamble C, Ghersi D, Ioannidis JP, Simes J, Williamson PR: Systematic review of the empirical evidence of study publication bias and outcome reporting bias.PLoS One 2008, 3:e3081.,

  40. Kyzas PA, Denaxa-Kyza D, Ioannidis JP: Almost all articles on cancer prognostic markers report statistically significant results. Eur J Cancer. 2007, 43: 2559-2579. 10.1016/j.ejca.2007.08.030.

    Article  PubMed  Google Scholar 

  41. Shaw GM, Lu W, Zhu H, Yang W, Briggs FB, Carmichael SL, Barcellos LF, Lammer EJ, Finnell RH: 118 SNPs of folate-related genes and risks of spina bifida and conotruncal heart defects.BMC Med Genet 2009, 10:49.,

  42. Lu W, Guzman AR, Yang W, Chapa CJ, Shaw GM, Greene RM, Pisano MM, Lammer EJ, Finnell RH, Zhu H: Genes encoding critical transcriptional activators for murine neural tube development and human spina bifida: a case–control study.BMC Med Genet 2010, 11:141.,

  43. Martinez CA, Northrup H, Lin JI, Morrison AC, Fletcher JM, Tyerman GH, Au KS: Genetic association study of putative functional single nucleotide polymorphisms of genes in folate metabolism and spina bifida. Am J Obstet Gynecol. 2009, 201 (394): e311-e391.

    Google Scholar 

  44. Carter TC, Pangilinan F, Troendle JF, Molloy AM, VanderMeer J, Mitchell A, Kirke PN, Conley MR, Shane B, Scott JM, Brody LC, Mills JL: Evaluation of 64 candidate single nucleotide polymorphisms as risk factors for neural tube defects in a large Irish study population. Am J Med Genet A. 2011, 155A: 14-21. 10.1002/ajmg.a.33755.

    Article  PubMed  Google Scholar 

  45. Pangilinan F, Molloy AM, Mills JL, Troendle JF, Parle-McDermott A, Signore C, O’Leary VB, Chines P, Seay JM, Geiler-Samerotte K, Mitchell A, VanderMeer JE, Krebs KM, Sanchez A, Cornman-Homonoff J, Stone N, Conley M, Kirke PN, Shane B, Scott JM, Brody LC: Evaluation of common genetic variants in 82 candidate genes as risk factors for neural tube defects.BMC Med Genet 2012, 13:62.,

  46. Spielman RS, McGinnis RE, Ewens WJ: Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM). Am J Hum Genet. 1993, 52: 506-516.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Pangilinan F, Mitchell A, VanderMeer J, Molloy AM, Troendle J, Conley M, Kirke PN, Sutton M, Sequeira JM, Quadros EV, Scott JM, Mills JL, Brody LC: Transcobalamin II receptor polymorphisms are associated with increased risk for neural tube defects. J Med Genet. 2010, 47: 677-685. 10.1136/jmg.2009.073775.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Hol FA, van der Put NM, Geurds MP, Heil SG, Trijbels FJ, Hamel BC, Mariman EC, Blom HJ: Molecular genetic analysis of the gene encoding the trifunctional enzyme MTHFD (methylenetetrahydrofolate-dehydrogenase, methenyltetrahydrofolate-cyclohydrolase, formyltetrahydrofolate synthetase) in patients with neural tube defects. Clin Genet. 1998, 53: 119-125. 10.1111/j.1399-0004.1998.tb02658.x.

    Article  CAS  PubMed  Google Scholar 

  49. Saavedra-Matiz CA, Isabelle JT, Biski CK, Duva SJ, Sweeney ML, Parker AL, Young AJ, Diantonio LL, Krein LM, Nichols MJ, Caggana M: Cost-effective and scalable DNA extraction method from dried blood spots. Clin Chem. 2013, 59: 1045-1051. 10.1373/clinchem.2012.198945.

    Article  CAS  PubMed  Google Scholar 

  50. Browne ML, Carter TC, Kay DM, Kuehn D, Brody LC, Romitti PA, Liu A, Caggana M, Druschel CM, Mills JL: Evaluation of genes involved in limb development, angiogenesis, and coagulation as risk factors for congenital limb deficiencies. Am J Med Genet A. 2012, 158A: 2463-2472. 10.1002/ajmg.a.35565.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Gorroochurn P, Hodge SE, Heiman GA, Durner M, Greenberg DA: Non-replication of association studies: “pseudo-failures” to replicate?. Genet Med. 2007, 9: 325-331. 10.1097/GIM.0b013e3180676d79.

    Article  PubMed  Google Scholar 

  52. Zollner S, Pritchard JK: Overcoming the winner’s curse: estimating penetrance parameters from case–control data. Am J Hum Genet. 2007, 80: 605-615. 10.1086/512821.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Au KS, Ashley-Koch A, Northrup H: Epidemiologic and genetic aspects of spina bifida and other neural tube defects. Dev Disabil Res Rev. 2010, 16: 6-15. 10.1002/ddrr.93.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Shaw GM, Rozen R, Finnell RH, Wasserman CR, Lammer EJ: Maternal vitamin use, genetic variation of infant methylenetetrahydrofolate reductase, and risk for spina bifida. Am J Epidemiol. 1998, 148: 30-37. 10.1093/oxfordjournals.aje.a009555.

    Article  CAS  PubMed  Google Scholar 

  55. Volcik KA, Shaw GM, Zhu H, Lammer EJ, Laurent C, Finnell RH: Associations between polymorphisms within the thymidylate synthase gene and spina bifida. Birth Defects Res A Clin Mol Teratol. 2003, 67: 924-928. 10.1002/bdra.10029.

    Article  CAS  PubMed  Google Scholar 

  56. Lupo PJ, Canfield MA, Chapa C, Lu W, Agopian AJ, Mitchell LE, Shaw GM, Waller DK, Olshan AF, Finnell RH, Zhu H: Diabetes and obesity-related genes and the risk of neural tube defects in the national birth defects prevention study. Am J Epidemiol. 2012, 176: 1101-1109. 10.1093/aje/kws190.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Ubbink JB, Christianson A, Bester MJ, Van Allen MI, Venter PA, Delport R, Blom HJ, van der Merwe A, Potgieter H, Vermaak WJ: Folate status, homocysteine metabolism, and methylene tetrahydrofolate reductase genotype in rural South African blacks with a history of pregnancy complicated by neural tube defects. Metabolism. 1999, 48: 269-274. 10.1016/S0026-0495(99)90046-X.

    Article  CAS  PubMed  Google Scholar 

  58. Yan L, Zhao L, Long Y, Zou P, Ji G, Gu A, Zhao P: Association of the maternal MTHFR C677T polymorphism with susceptibility to neural tube defects in offsprings: evidence from 25 case–control studies.PLoS One 2012, 7:e41689.,

  59. van der Put NM, Eskes TK, Blom HJ: Is the common 677C– > T mutation in the methylenetetrahydrofolate reductase gene a risk factor for neural tube defects? A meta-analysis. QJM. 1997, 90: 111-115. 10.1093/qjmed/90.2.111.

    Article  CAS  PubMed  Google Scholar 

  60. Ouyang S, Li Y, Liu Z, Chang H, Wu J: Association between MTR A2756G and MTRR A66G polymorphisms and maternal risk for neural tube defects: a meta-analysis. Gene. 2013, 515: 308-312. 10.1016/j.gene.2012.11.070.

    Article  CAS  PubMed  Google Scholar 

  61. van der Linden IJ, den Heijer M, Afman LA, Gellekink H, Vermeulen SH, Kluijtmans LA, Blom HJ: The methionine synthase reductase 66A > G polymorphism is a maternal risk factor for spina bifida. J Mol Med (Berl). 2006, 84: 1047-1054. 10.1007/s00109-006-0093-x.

    Article  Google Scholar 

  62. Ouyang S, Liu Z, Li Y, Wu J: Meta-analyses on the association of MTR A2756G and MTRR A66G polymorphisms with neural tube defect risks in Caucasian children. J Matern Fetal Neonatal Med. 2013, 26 (12): 1166-1170. 10.3109/14767058.2013.777699.

    Article  CAS  PubMed  Google Scholar 

Download references


These studies would not be possible without the participation of the affected families, and their recruitment by the United Kingdom Association of Spina Bifida and Hydrocephalus. The authors acknowledge research support from the Intramural Research programs of the National Human Genome Research Institute, the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD contract number: N01-DK-7-3431), and the Health Research Board, Ireland. We thank Dr. Charlotte Druschel, Director of the Congenital Malformations Registry, for identifying and reviewing cases for this investigation, and Sandra D. Richardson at the NYSCMR for subject selection and data management. We also thank April J. Atkins, Adam C. Gearhart, and Robert J. Sicko at the Wadsworth Center, New York State Department of Health for providing laboratory and technical assistance.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Lawrence C Brody.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

FP, AMM, JLM, PNK, BS, JMS and LCB formulated the study design. FP drafted the manuscript. MS managed subject recruitment and data collection for the UK subjects. AMM managed the UK DNA samples and their selection. DK, MB, MC and ECM managed and selected the NYS DNA samples, and contracted their genotyping. FP, AP-M, and HA performed genotyping. JFT contributed to study design and performed the statistical analyses. All authors read, edited and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S1.: Genotype counts for candidate SNPs in the NTD cases, fathers and mothers from the United Kingdom (England, Wales, and Northern Ireland). (XLSX 46 KB)


Additional file 2: Table S2.: Genotype counts for candidate SNPs in the New York state NTD cases and controls by ethnicity and each analyzed NTD subtype. (XLSX 49 KB)

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pangilinan, F., Molloy, A.M., Mills, J.L. et al. Replication and exploratory analysis of 24 candidate risk polymorphisms for neural tube defects. BMC Med Genet 15, 102 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: