Next Article in Journal
Risk of Ischemic Heart Disease Associated with Primary Dysmenorrhea: A Population-Based Retrospective Cohort Study
Next Article in Special Issue
A Paradigm Shift in Understanding the Pathological Basis of Autism Spectrum Disorder: From the Womb to the Tomb
Previous Article in Journal
Identification of Patients with Potential Atrial Fibrillation during Sinus Rhythm Using Isolated P Wave Characteristics from 12-Lead ECGs
Previous Article in Special Issue
Understanding Challenging Behaviors in Autism Spectrum Disorder: A Multi-Component, Interdisciplinary Model
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JPM
Volume 12
Issue 10
10.3390/jpm12101609
jpm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Folate Pathway Gene Single Nucleotide Polymorphisms and Neural Tube Defects: A Systematic Review and Meta-Analysis

by
Ahmad K. Almekkawi
1,2,*,
Marwa W. AlJardali
3,
Hicham M. Daadaa
4,
Alison L. Lane
1,2,
Ashley R. Worner
1,2,
Mohammad A. Karim
1,2,
Adrienne C. Scheck
1,2 and
Richard E. Frye
5
1
Phoenix Children’s Hospital, 1919 E. Thomas Rd, Ambulatory Building, Phoenix, AZ 85016, USA
2
College of Medicine, Department of Child Health, University of Arizona, 475 N. 5th Street, Phoenix, AZ 85004, USA
3
Rose-Marie Chagoury School of Medicine, Lebanese American University Gilbert, Byblos 1102 2801, Lebanon
4
Department of Oncology, St James University Hospital, Beckett St., Harehills, Leeds LS9 7TF, UK
5
Rossignol Medical Center, Phoenix, AZ 85050, USA
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2022, 12(10), 1609; https://doi.org/10.3390/jpm12101609
Submission received: 4 August 2022 / Revised: 5 September 2022 / Accepted: 22 September 2022 / Published: 29 September 2022
(This article belongs to the Special Issue A Personalized Medicine Approach to the Diagnosis and Management of Autism Spectrum Disorder: Beyond Genetic Syndromes)
Browse Figures
Review Reports Versions Notes

Abstract

:
Neural tube defects (NTDs) are congenital abnormalities in the central nervous system. The exact etiology of NTDs is still not determined, but several genetic and epigenetic factors have been studied. Folate supplementation during gestation is recommended to reduce the risk of NTDs. In this review we examine single nucleotide polymorphisms (SNPs) of the genes in the folate pathway associated with NTD. We reviewed the literature for all papers discussing both NTDs and SNPs in the folate pathway. Data were represented through five different genetic models. Quality assessment was performed using the Newcastle–Ottawa Scale (NOS) and Cohen’s Kappa inter-rater coefficient assessed author agreement. Fifty-nine papers were included. SNPs in MTHFR, MTRR, RFC genes were found to be highly associated with NTD risk. NOS showed that high quality papers were selected, and Kappa Q-test was 0.86. Our combined results support the notion that SNPs significantly influence NTDs across the population, particularly in Asian ethnicity. Additional high-quality research from diverse ethnicities is needed and meta-regression analysis based on a range of criteria may provide a more complete understanding of the role of folate metabolism in NTDs.

1. Introduction

Folates are part of the B-group family of vitamins consisting of an aromatic pteridine ring bound to p-aminobenzoic acid and a glutamate residue through a methylene bond. Folate is an essential vitamin; it cannot be synthesized by human metabolism and needs to be derived entirely from the diet. Folic acid is the synthetic source of folate found in fortified foods and many vitamin supplements [1]. Folates are very important in human health and have been implicated in many developmental disorders. Adequate folate during pregnancy is necessary to prevent congenital defects such as neural tube defects (NTDs) as well as neurodevelopmental disorders such as autism spectrum disorder (ASD) [1]. Folate is essential during childhood and has been implicated in ASD and neurological disorders such as cerebral folate transport deficiency (CFD) [2]. Folate is also implicated in adult health, being linked to dementia [3].

1.1. Folate Absorption and Transport

Upon ingestion, folate is bound to folate-binding protein (FBP). Folyl poly-γ-glutamate carboxypeptidase (FGCP), which is encoded by the folate hydrolase 1 (FOLH1) gene, is an apical brush border transmembrane enzyme found in the small intestine which releases folate from the bonded proteins in the intestinal lumen. Two folate carriers located on the apical surface of the enterocyte, the proton-coupled folate transporter (PCFT) and the reduced folate carrier (RFC), transport folate from the intestinal lumen into the enterocytes [4]. The multidrug resistance associated protein 3 (MRP3), located on the basolateral membrane of the enterocyte, transports folate out of the enterocyte where it is transported to the liver via the hepatic portal vein. Folate is transported into hepatocytes primarily by the PCFT and RFC. Once in the hepatocyte folate is either metabolized, secreted into the bile [5] or converted into a different form of folate and released into the bloodstream.
In addition to the RFC and PCFT transmembrane folate carriers, there are several receptor mediated folate transporters. Specifically, three glycoproteins known as folate receptor alpha (FRα), beta (FRβ) and gamma (FRγ) mediate cellular folate uptake [1].

1.2. Folate Interconversion: Reduction of Folate and the Folate Cycle

The human body’s metabolic systems use the reduced form of folate. Folic acid, the synthetic form of folate, is an oxidized form of folate. Oxidized folates, like folic acid, need to be reduced to tetrahydrofolate (THF) through dihydrofolate reductase (DHF). THF is the entry into the folate cycle. The folate cycle supports two important metabolic pathways, the production of purines and the methylation cycle. Through the folate cycle, THF is converted to 5,10-methylene–THF by the pyridoxal phosphate dependent enzyme serine hydroxymethyltransferase (SHMT) using L-serine and producing glycine in the process. 5,10-methylene–THF interconverts to 5-formyltetrahydrofolate, commonly known as leucovorin or folinic acid, which is the precursor for purine nucleotides. 5,10-methylene–THF is converted into 5-methyltetrahydrofolate (5-methyl-THF) through the enzyme methylenetetrahydrofolate reductase (MTHFR). 5-methyl-THF is the folate form that drives the methylation cycle. Methionine synthase (MTR), a vitamin B12-dependent enzyme, converts 5-methyl-THF back to THF to complete the folate cycle, and in the process produces methionine from homocysteine (HCY).

1.3. Methylation Cycle

Methyl donors are essential for regulation of gene expression and enzyme activity. Methionine produces S-adenosyl methionine (SAM), the body’s major methyl donor, by methionine adenosyltransferase (MAT) through an adenosine triphosphate (ATP) dependent reaction. Donating a methyl group results in the production of S-adenosylhomocysteine (SAH). SAH is then converted in a reversible reaction to HCY and adenosine by S-adenosylhomocysteine hydrolase (SAHH) through a reversible reaction. In addition to the reaction using 5-methyl-THF, methionine can also be produced from HCY using betaine-homocysteine methyltransferase (BHMT), an enzyme exclusively found in the liver and the kidney. BHMT converts betaine and homocysteine to dimethylglycine and methionine. Importantly, the methylation cycle is connected to transsulfuration reactions which produce glutathione, the major antioxidant of the body. HCY produces cysteine, the rate limiting precursor to GSH production, through two pyridoxal phosphate-dependent reactions catalyzed by cystathionine β-synthase (CBS) and γ-cystathionase (CTH).

1.4. Neurodevelopmental Disorders and Neural Tube Defects: Links to Folate Metabolism

Neurulation, the process of closure of the neural tube, occurs 21 and 28 days after conception. Interruption of this process can cause neural tube defects (NTDs), resulting in congenital deformities of the brain and spinal cord [6]. NTDs include encephalocele, anencephaly, and spina bifida, which are among the most common and serious birth abnormalities and associated with a high mortality rate [7]. The prevalence of NTDs is 1 per 1000 live births across the world, representing 20 to 25% of birth abnormalities in China [8]. Previous studies have shown that both environmental and genetic variables play a significant role in the pathogenesis of NTDs.
Folate was found to be severely insufficient in the diets of women who had offspring with NTDs, according to ground-breaking research by Smithells et al. [9]. Furthermore, it was found that folic acid fortification could prevent NTDs [10,11,12]. Berry et al. [13] found that folic acid supplements, at doses ranging from 360 mcg to 4000 mcg, lowered the incidence of NTDs by 50–75%. Folic acid supplementation was found to reduce both the initial occurrence and recurrence of NTDs [14,15].
To participate in the methylation cycle and contribute to the metabolism of folate, folic acid must first transform into the naturally bioactive form known as THF. A deficiency in neural tube development may be brought on by inhibition in the folate metabolism pathway. The genes in the folate pathway are being extensively studied for this reason, and it has been shown that proteins encoded by certain important genes, including MTHFR, RFC, and methionine synthase reductase (MTRR), work together with MTR to transfer the methyl group to homocysteine [16,17]. When the first study on single-nucleotide polymorphisms (SNPs) appeared in 1995, the relationship between genetic variance and NTDs was discovered [18]. Since then, numerous studies have demonstrated the highest correlation between NTDs and aberrant gene mutations that alter folate metabolism [14,19,20,21]. Due to their crucial roles in the folate metabolism pathway, SNPs C677T and A1298C in MTHFR, A2756G in MTR, A66G in MTRR, and A80G in RFC have received the most attention and may account for a sizeable fraction of the risk of having NTDs [18]. Growing data from epidemiological case–control studies have shown that genetic factors may account for up to 70% of NTD prevalence [22].
It’s interesting that despite numerous research studies examining the relationship between NTDs and SNPs, no consensus has been achieved. Folate metabolism may be inhibited by folate pathway gene SNPs [23,24,25,26,27]. However, numerous follow-up investigations failed to confirm the connection [28,29,30,31]. To investigate the relationship between important SNPs in the main folate pathway genes and NTDs, we carried out this systematic review and meta-analysis.

2. Materials and Methods

To ensure the rigor of this meta-analysis, we followed the Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) guidelines [32].

2.1. Search Strategy and Identification of Relevant Studies

We searched three databases, PubMed, EMBase and Web for published papers, from June 1996 to June 2022, which examined any SNPs relating to the folate pathway associated with NTDs. The search strategy was based on English keywords, “MTHFR”, “MTRR”, “RFC”, “folate pathway” “polymorphism” or “SNP” and “NTDs or Neural Tube Defects or Spina Bifida”. References for the articles and the collected studies were reviewed.
The following inclusion criteria had to be met: (1) case–control research and cohort study design; (2) data on some, any, or all SNPs in the folate pathway; (3) presentation of data appropriate for the measurement of odds ratios (ORs); (4) consistent description of NTDs. Animal research, maternal studies, case reports and studies without original data were omitted. Duplicated studies were removed.

2.2. Data Extraction

All data were analyzed independently by two reviewers (A.K.A. & R.E.F.). The following information was extracted from the qualifying literature: year of publication, first author’s name, country, ethnicity, genotyping process, source of control, and case matching control variables. Allele or genotype counts were extracted or estimated from the published data for the re-calculation of crude ORs and their 95% confidence intervals (95% CIs) for the assessment of the association of SNPs in the folate pathway with NTDs.

2.3. Statistical Analysis

Only SNPs with 4 or more reports that included at least 5 homozygous cases for the risk allele were included in the analysis. ORs with 95% CIs were used to evaluate the strength of the link between the SNP and the risk of NTDs. The pooled OR values were determined using the Z-test, and p < 0.05 was regarded as statistically significant. All five genetic models [allele (C vs. A), homozygote (CC vs. AA), heterozygote (CA vs. AA), dominant (CC + CA vs. AA), and recessive (CC vs. CA + AA)] were used to compute the aggregated estimate of the OR and corresponding 95% CI for the association. Egger’s exact test was used to assess heterogenicity between papers. The pooled ORs and 95% CIs were calculated using a fixed-effects model (the Mantel-Haenszel technique) when there was no substantial heterogeneity (p > 0.1 or I2 50%) between trials; otherwise, a random-effects model (the DerSimonian-Laird method) was chosen. We conducted subgroup analyses by ethnicity to examine the impact of between-study heterogeneity on our pooled data. By using Pearson’s chi-square test, Hardy–Weinberg equilibrium (HWE) in controls was assessed, with p < 0.05 indicating a departure from HWE. By excluding each selected study one at a time and reanalyzing the remaining papers, sensitivity analysis was used to determine the stability of the meta-analysis results. Additionally, the sensitivity analysis was carried out by removing studies that violated the HWE. The Egger’s linear regression test and the Egger’s quantitative test were used to assess any potential publication bias. A significant publication bias was deemed to exist when the Egger’s test plot was asymmetric, and the p-value for the Egger’s test was less than 0.05.

2.4. Quality Assessment

Analysis of the quality of the final selected papers was evaluated using the Newcastle Ottawa Scale [33]. The scale consists of nine questions that rate the outcome, comparability, and choice of paper. Using a scale of 0 to 10, two authors (A.L.L. & A.R.W.), independently evaluated the literature. To gauge inter-rater reliability between reviewers, the Cohen’s kappa coefficient was determined.

3. Results

3.1. Characteristics of Eligible Studies

The selection and identification of literature are shown in Figure S1. The search resulted in 1284 pertinent articles. After the initial screening, 540 of these papers were disqualified based on their titles and abstracts. Then, full-text papers were found and evaluated in accordance with the inclusion criteria. Of the remaining 111 full-text articles, 52 studies were excluded because they were irrelevant, provided data for other SNPs or only parental SNPs, analyzed siblings, twins, and trios of parents, provided invalid data, or overlapped with other publications. 59 case–control papers in total were chosen for this meta-analysis (Figure S1).

3.2. Meta-Analysis Results

Study characteristics of the papers included in the meta-analysis for association between the polymorphism and NTDs are listed in Table S1. Pooled data showed that the MTHFR C677T SNP was significantly associated with an increased risk of NTD in all of the five genetic models, i.e., allele (T vs. C: OR = 0.73, 95% CI 0.67–0.79, p ≤ 0.0001), recessive (TT vs. TC + CC: OR = 0.70, 95% CI 0.62–0.78, p ≤ 0.0001), homozygote (TT vs. CC: OR = 0.55, 95% CI 0.46–0.64, p ≤ 0.00001), heterozygote (TT vs. TC: OR = 0.76, 95% CI 0.67–0.86, p ≤ 0.01), dominant (TT + TC vs. CC: OR = 0.64, 95% CI 0.56–0.74, p ≤ 0.0001) (Figure 1). Egger’s test shows no bias (p > 0.10) and the funnel plot shows no heterogeneity (Figure S2).
For MTHFR A1298C, none of the models were significant (Figure 2). Egger’s test showed no bias (p > 0.10) but funnel plots showed heterogeneity (Figure S3).
MTR A2756G SNP has been associated with NTDs. Our pooled analysis showed no association between MTR and NTDs. For example, the recessive model was not significant (GG vs. GA + AA, OR = 0.99, 95% CI 0.84–1.17, p > 0.10) (Figure 3). The Egger’s test showed no bias (p = 0.75) and the funnel plot showed no heterogeneity (Figure S4).
For the MTRR A66G SNP, the dominant model was significant (GG + GA vs. AA, OR = 1.27, 95% CI 1.05–1.53, p = 0.01) (Figure 4). The Egger’s test showed no bias (p = 0.73) but the funnel plot show heterogeneity (Figure S5).
MTHFD1 was examined in 6 studies. There was no association between MTHFD1 and NTDs in any of the 5 different allele models. Egger’s test show no bias (p = 0.49) and funnel plots showed little heterogeneity for the dominate model (Figure S6).
RFC was associated with NTDs in 3 of the 5 models: allele (G vs. A: OR = 0.77, 95% CI 0.66–0.90, p < 0.01; Figure 5A), dominant (GG + GA vs. AA: OR = 0.62, 95% CI 0.49–0.80, p < 0.01; Figure 5B), and homozygote (GG vs. AA: OR = 0.62, 95% CI 0.49–0.80, p = 0.01; Figure 5C). It was not significant for the heterozygote (GG vs. GA: OR = 0.95, 95% CI 0.73–1.25, p > 0.10; Figure 5D) or recessive (GG vs. GA + AA: OR = 0.82, 95% CI 0.64–1.06, p > 0.10) models. Funnel plot showed no heterogeneity and Egger’s test show no bias for the allele model (p = 0.53; Figure S7).
Subgroup analysis was conducted using all 5 genetic models to examine differences between ethnicities and SNPs in the SNPs where associations were found in the overall analysis (see Table 1). Ethnic differences were found in the association between the SNPs and NTDs. The association between MTHFR C677T and NTDs was significant (p < 0.01) for Asian and Caucasian ethnicities in 3 of the 5 genetic models. Only Asian ethnicity was significant in the recessive and heterozygous models for MTHFR. The RFC association with NTDs was only seen in the Asian ethnicity and it was significant in 4 of the 5 models. None of the models were significant for Caucasian ethnicity.

3.3. Quality Assessment and Inter-Rater Agreement

Newcastle-Ottowa Scale was conducted to assess quality of the literature by 2 authors (A.L.L. & A.R.W.). Results were represented as a scale of 1–10. Majority of the papers selected had a score of 7 and above which represents good quality (Table S2). Cohen’s Kappa Inter-rater agreement coefficient calculated was 0.82, which represents agreeability between raters.

4. Discussion

NTDs are the most prevalent congenital disorder, second only to cardiovascular disorders [69]. The effectiveness of folic acid-containing vitamin supplementation during pregnancy for significantly lowering the risk of NTDs in offspring was been proven over the last three decades [27]. The production and maintenance of DNA as well as the methylation of DNA, RNA, and protein are only a few of the vital biological pathways that depend on folate [70]. Due to this, genes associated with folate have been studied closely [71]. To our knowledge, this work is the first to combine research on SNPs in folate metabolism associated with NTDs. In a comprehensive meta-analysis of case–control studies, our findings show a substantial association between MTHFR C677T and NTDs. Further, the ensuing cumulative meta-analysis provide strong support for the association. Our comprehensive meta-analysis includes research on the MTHFR variant A1298C, the MTRR variants A66G, the MTR variant A2756G, and the RFC variant A80G.
This study identifies MTHFR C677T, MTRR A66G and RFC A80G SNPs as risk factors for NTDs. Our findings are summarized in Figure 6. One caveat is that RFC A80G was only significant for Asian ethnicity but not Caucasian ethnicity, potentially suggesting an underappreciated effect of ethnicity for SNPs and NTDs.
Birth abnormalities can result from gene mutations in essential enzymes that control folate metabolism and transportation. Several mutations can severely affect MTHFR activity, reducing folate status and potentially explaining a quarter of the frequency of NTDs [72,73]. TT mutation at nucleotide (nt) 677 can diminish the MTHFR activity by more than 65%, as was initially shown by Frosst et al. [74]. A less potent version of this effect was observed in A1298C [75]. The results of this comprehensive meta-analysis suggested that MTHFR C677T may be a risk factor for NTDs.
Loss-of-function mutations in MTR and mutations of the MTR chaperone MTRR may influence homocysteine levels resulting in severe disease phenotypes [76,77]. Our analysis potentially supports an association of MTRR A66G as a risk factor for NTDs. The absence of a previously identified association may be due to several factors. First, previous studies used small sample sizes, raising the possibility of missing an effect due to sample heterogenicity. Second, the severity of NTD subtypes vary from livebirth to stillbirth. Therefore, if studies only include livebirths (i.e., less severe cases), the impact of genetic variations on risk of NTDs may be understated. Third, the results could be impacted by potential gene–gene and gene–environment interactions.
The association of RCF SNPs with NTDs is still unclear. Relton et al. [28] found that the 80A allele, not the 80G allele, increased the risk of NTDs, in contrast to most studies. We ruled out the chance that they might have reported the reverse strand allele. Furthermore, the heterogeneity was dramatically decreased when this study was excluded, demonstrating that it was the primary cause of the heterogeneity. Thus, since our analysis included the Relton study, excluding it should produce a more significant effect. This could be the explanation for the fact that the RFC SNP was only associated with NTDs for Asian ethnicity since the Relton study included a Caucasian population. Clearly further studies are needed to examine this association between the RFC and NTDs, particularly with respect to its variation across ethnicities.
MTHFD1, also known as “C1-THF synthase,” is a trifunctional cytoplasmic enzyme that is nicotinamide adenine dinucleotide phosphate (NADP) dependent. Additionally, the functional SNP MTHFD1 G1958A has been investigated in several populations, and earlier research has shown that MTHFD1 genetic variations can cause a variety of disorders [78,79,80,81,82,83,84,85,86]. According to Carroll et al. [87], both case (genotype and allele frequencies) and maternal (allele frequencies only) groups have highly significant associations between NTD risk and the MTHFD1 SNP G1958A (R653Q). The idea of maternal folate metabolism increasing the risk of NTDs is also supported by Green et al. [88] where it was found that both genetic and epigenetic mechanisms contributed to the increased risk of NTDs and overall birth defect risk in the maternal lineage as more genotypes were linked to poorer folate metabolism in maternal relatives as compared to paternal relatives [89]. While the other known variants in the MTHFD2 gene are located within intronic regions that are unlikely to contain regulatory sequences and are not obvious candidates for a disease association, Parle-McDermott et al. [88] confirmed that some SNPs on MTHFD1 may have undergone evolutionary selection. However, our paper found no significant association between MTHFD1 and NTDs in all 5 genetic models.
A good technique for combining data from all relevant research to increase statistical power is meta-analysis. However, a few concerns with this meta-analysis need to be considered. First, the number of mixed population studies that were considered for subgroup analysis was limited, with the bulk of the subjects being Asians. Additionally, following stratification, the sample sizes for Caucasian ethnicity was reduced, which may have limited the statistical power of our findings. Second, we were unable to locate any articles written in languages other English and Chinese, which could have influenced our results. Third, some studies with unfavorable findings tend not to be published (i.e., publication bias), which could bias the available data. Fourth, NTDs encompass a wide spectrum of malformations, including encephalocele and spinal bifida. Additionally, the severity and prognosis of NTDs vary along a continuum, from anencephalic stillbirth to meningocele spina bifida with mild symptoms. Therefore, if studies only include live births and less severe instances, the effects of SNPs on the development of NTD types may be understated. Finally, epigenetic variables may be linked to the occurrence of NTDs. Due to the lack of data related to gene–gene and gene–environment interactions in the formation of NTDs these potential contributors to NTDs could not be studied.
Despite these caveats we believe our meta-analysis is robust in its findings. First, sensitivity analysis showed that none of the included studies independently changed the meta-analysis results. Indeed, removing the HWE-violating studies had no effect on the outcomes of our pooled data. This suggests that the pooled results are sufficiently reliable and stable. Second, this is the most thorough meta-analysis to include data from earlier research on the connection of the SNPs in the folate pathway with risk of NTDs and spina bifida. In conclusion, our combined findings are consistent with the idea that SNPs are significantly associated with NTDs across the population, notably in the Asian ethnicity.

5. Conclusions

The current study, which performed a cumulative meta-analysis of 42 studies of MTHFR C677T that provided results with a relatively narrow 95% confidence range, significantly supports the association between the MTHFR C677T SNPs and NTD risk. We found a possible link between RFC A80G and the risk of NTDs but this was only confirmed for Asian ethnicity, potentially because of the contradictory results in studies performed in the Caucasian population. Significant results were found for the MTRR SNP but for only one model. In contrast, we were unable to identify an association for the remaining SNPs. Future work will benefit from the incorporation of additional high-quality studies from various ethnicities when they become available. Additionally, when further studies start to incorporate more environmental variables, meta-regression analysis based on a variety of additional parameters of interest may help reduce heterogenicity between studies.
Most importantly, we have come to appreciate the differences between folic acid, which is an oxidized folate compound with limited bioavailability. Indeed, we are understanding that reduced folate compounds such as folinic acid, leucovorin and methyltetrahydrofolate may be able to bypass blockage in the folate pathway in order to prevent and treat disease [1]. Thus, a better understanding of the variations in the folate pathway will lead to a more personalized medicine approach to prenatal care that may reduce the occurrence and severity of NTD and related neurodevelopmental disorders.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jpm12101609/s1, Table S1: Table representing eligible study characteristics. Table S2: Newcastle Ottawa Scale describing the quality of the literature included. Results were represented as a scale from 0–10 and conducted by 2 independent authors. High quality described with a score of greater than 7 was seen in most of the literature included in the study. Figure S1: Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA) flowchart. Figure S2: Funnel plot and Egger’s Test for bias and heterogeneity. Figure S3: Funnel plot and Egger’s Test for bias and heterogeneity. Figure S4: Funnel plot and Egger’s Test for bias and heterogeneity. Figure S5: Funnel plot and Egger’s Test for bias and heterogeneity. Figure S6: Funnel plot and Egger’s Test for bias and heterogeneity. Figure S7: Funnel plot and Egger’s Test for bias and heterogeneity.

Author Contributions

A.K.A. was responsible for paper conceptualization, determining article eligibility, paper design, data collection, statistical analysis and writing the manuscript with revision. M.W.A. was responsible for second stage eligibility, first draft design, and quality assessment. H.M.D. was responsible for drafting the introduction of the first draft and reviewing the manuscript. A.L.L. and A.R.W. were responsible for completing the quality assessment and reviewing the manuscript. M.A.K. and A.C.S. reviewed the manuscript. R.E.F. is the principal investigator in this paper he is responsible for manuscript conceptualization, revision, editing and finalization. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Frye, R.E.; Slattery, J.C.; Quadros, E.V. Folate metabolism abnormalities in autism: Potential biomarkers. Biomark. Med. 2017, 11, 687–699. [Google Scholar] [CrossRef] [PubMed]
  2. Rossignol, D.A.; Frye, R.E. Cerebral Folate Deficiency, Folate Receptor Alpha Autoantibodies and Leucovorin (Folinic Acid) Treatment in Autism Spectrum Disorders: A Systematic Review and Meta-Analysis. J. Pers. Med. 2021, 11, 1141. [Google Scholar] [CrossRef] [PubMed]
  3. Gil Martínez, V.; Avedillo Salas, A.; Santander Ballestín, S. Vitamin Supplementation and Dementia: A Systematic Review. Nutrients 2022, 14, 1033. [Google Scholar] [CrossRef] [PubMed]
  4. Visentin, M.; Diop-Bove, N.; Zhao, R.; Goldman, I.D. The intestinal absorption of folates. Annu. Rev. Physiol. 2014, 76, 251–274. [Google Scholar] [CrossRef] [PubMed]
  5. Boyer, J.L.; Soroka, C.J. Bile formation and secretion: An update. J. Hepatol. 2021, 75, 190–201. [Google Scholar] [CrossRef]
  6. Christianson, A.; Howson, C.P.; Modell, B. March of Dimes: Global Report on Birth Defect, the Hidden Toll of Dying and Disabled Children; March of Dimes Birth Defects Foundation: New York, NY, USA, 2006; pp. 10–16. [Google Scholar]
  7. Pitkin, R.M. Folate and neural tube defects. Am. J. Clin. Nutr. 2007, 85, 285S–288S. [Google Scholar] [CrossRef]
  8. Copp, A.J.; Greene, N.D.; Murdoch, J.N. The genetic basis of mammalian neurulation. Nat. Rev. Genet. 2003, 4, 784–793. [Google Scholar] [CrossRef]
  9. Smithells, R.W.; Sheppard, S.; Schorah, C.J. Vitamin dificiencies and neural tube defects. Arch. Dis. Child. 1976, 51, 944–950. [Google Scholar] [CrossRef]
  10. Wolff, T.; Witkop, C.T.; Miller, T.; Syed, S.B. Folic acid supplementation for the prevention of neural tube defects: An update of the evidence for the US Preventive Services Task Force. Ann. Intern. Med. 2009, 150, 632–639. [Google Scholar] [CrossRef]
  11. Laurence, K. Prevention of neural tube defects by improvement in maternal diet and preconceptional folic acid supplementation. Prog. Clin. Biol. Res. 1985, 163, 383–388. [Google Scholar]
  12. Christiansen, M.; Garne, E. Prevention of neural tube defects with periconceptional folic acid supplementation in Europe. Ugeskr. Laeger 2005, 167, 2875–2876. [Google Scholar]
  13. Berry, R.J.; Li, Z.; Erickson, J.D.; Li, S.; Moore, C.A.; Wang, H.; Mulinare, J.; Zhao, P.; Wong, L.-Y.C.; Gindler, J. Prevention of neural-tube defects with folic acid in China. N. Engl. J. Med. 1999, 341, 1485–1490. [Google Scholar] [CrossRef]
  14. Kibar, Z.; Capra, V.; Gros, P. Toward understanding the genetic basis of neural tube defects. Clin. Genet. 2007, 71, 295–310. [Google Scholar] [CrossRef] [PubMed]
  15. De-Regil, L.M.; Peña-Rosas, J.P.; Fernández-Gaxiola, A.C.; Rayco-Solon, P. Effects and safety of periconceptional oral folate supplementation for preventing birth defects. Cochrane Database Syst Rev. 2015, 12. [Google Scholar] [CrossRef]
  16. Bassuk, A.G.; Kibar, Z. Genetic basis of neural tube defects. Semin. Pediatric Neurol. 2009, 16, 101–110. [Google Scholar] [CrossRef] [PubMed]
  17. Ulrich, C.M.; Robien, K.; McLeod, H.L. Cancer pharmacogenetics: Polymorphisms, pathways and beyond. Nat. Rev. Cancer 2003, 3, 912–920. [Google Scholar] [CrossRef]
  18. Molloy, A.M.; Brody, L.C.; Mills, J.L.; Scott, J.M.; Kirke, P.N. The search for genetic polymorphisms in the homocysteine/folate pathway that contribute to the etiology of human neural tube defects. Birth Defects Res. Part A Clin. Mol. Teratol. 2009, 85, 285–294. [Google Scholar] [CrossRef]
  19. Richter, B.; Stegmann, K.; Röper, B.; Böddeker, I.; Ngo, E.T.; Koch, M.C. Interaction of folate and homocysteine pathway genotypes evaluated in susceptibility to neural tube defects (NTD) in a German population. J. Hum. Genet. 2001, 46, 105–109. [Google Scholar] [CrossRef] [PubMed]
  20. Collin, S.M.; Metcalfe, C.; Zuccolo, L.; Lewis, S.J.; Chen, L.; Cox, A.; Davis, M.; Lane, J.; Donovan, J.; Smith, G.D. Association of Folate-Pathway Gene Polymorphisms with the Risk of Prostate Cancer: A Population-Based Nested Case-Control Study, Systematic Review, and Meta-analysisFolate-Pathway Gene Polymorphisms and Prostate Cancer. Cancer Epidemiol. Biomark. Prev. 2009, 18, 2528–2539. [Google Scholar] [CrossRef]
  21. Van der Linden, I.J.; Afman, L.A.; Heil, S.G.; Blom, H.J. Genetic variation in genes of folate metabolism and neural-tube defect risk. Proc. Nutr. Soc. 2006, 65, 204–215. [Google Scholar] [CrossRef]
  22. Greene, N.D.; Copp, A.J. Neural tube defects. Annu. Rev. Neurosci. 2014, 37, 221. [Google Scholar] [CrossRef] [PubMed]
  23. Boduroglu, K.; Alanay, Y.; Alikasifoglu, M.; Aktas, D.; Tunçbilek, E. Analysis of MTHFR 1298A> C in addition to MTHFR 677C> T polymorphism as a risk factor for neural tube defects in the Turkish population. Turk. J. Pediatrics 2005, 47, 327. [Google Scholar]
  24. Kirke, P.N.; Mills, J.L.; Molloy, A.M.; Brody, L.C.; O’Leary, V.B.; Daly, L.; Murray, S.; Conley, M.; Mayne, P.D.; Smith, O. Impact of the MTHFR C677T polymorphism on risk of neural tube defects: Case-control study. BMJ 2004, 328, 1535–1536. [Google Scholar] [CrossRef] [Green Version]
  25. De Marco, P.; Calevo, M.G.; Moroni, A.; Arata, L.; Merello, E.; Finnell, R.H.; 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. [Google Scholar] [CrossRef]
  26. Parle-McDermott, A.; Mills, J.L.; Kirke, P.N.; O’Leary, V.B.; Swanson, D.A.; Pangilinan, F.; Conley, M.; Molloy, A.M.; Cox, C.; Scott, J.M. Analysis of the MTHFR 1298A→ C and 677C→ T polymorphisms as risk factors for neural tube defects. J. Hum. Genet. 2003, 48, 190–193. [Google Scholar] [CrossRef] [PubMed]
  27. Shaw, G.M.; Rozen, R.; Finnell, R.H.; Wasserman, C.R.; Lammer, E.J. Maternal vitamin use, genetic variation of infant methylenetetrahydrofolate reducatase, and risk for spina bifida. Am. J. Epidemiol. 1998, 148, 30–37. [Google Scholar] [CrossRef]
  28. Relton, C.; Wilding, C.; Pearce, M.; Laffling, A.; Jonas, P.; Lynch, S.; Tawn, E.; 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. [Google Scholar] [CrossRef]
  29. Godbole, K.; Gayathri, P.; Ghule, S.; Sasirekha, B.V.; Kanitkar-Damle, A.; Memane, N.; Suresh, S.; Sheth, J.; Chandak, G.R.; Yajnik, C.S. Maternal one-carbon metabolism, MTHFR and TCN2 genotypes and neural tube defects in India. Birth Defects Res. Part A: Clin. Mol. Teratol. 2011, 91, 848–856. [Google Scholar] [CrossRef]
  30. De Franchis, R.; Buoninconti, A.; Mandato, C.; Pepe, A.; Sperandeo, M.P.; Del Gado, R.; Capra, V.; Salvaggio, E.; Andria, G.; Mastroiacovo, P. The C677T mutation of the 5, 10-methylenetetrahydrofolate reductase gene is a moderate risk factor for spina bifida in Italy. J. Med. Genet. 1998, 35, 1009–1013. [Google Scholar] [CrossRef] [PubMed]
  31. Mornet, E.; Muller, F.; Lenvoisé-Furet, A.; Delezoide, A.-L.; Col, J.-Y.; Simon-Bouy, B.; Serre, J.-L. Screening of the C677T mutation on the methylenetetrahydrofolate reductase gene in French patients with neural tube defects. Hum. Genet. 1997, 100, 512–514. [Google Scholar] [CrossRef] [PubMed]
  32. Moher, D.; Altman, D.G.; Liberati, A.; Tetzlaff, J. PRISMA statement. Epidemiology 2011, 22, 128. [Google Scholar] [CrossRef] [PubMed]
  33. Stang, A. Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur. J. Epidemiol. 2010, 25, 603–605. [Google Scholar] [CrossRef] [PubMed]
  34. Van der Put, N.; Thomas, C.; Eskes, T.; Trijbels, F.; Steegers-Theunissen, R.; Mariman, E.; De Graaf-Hess, A.; Smeitink, J.; Blom, H. Altered folate and vitamin B12 metabolism in families with spina bifida offspring. QJM Mon. J. Assoc. Physicians 1997, 90, 505–510. [Google Scholar] [CrossRef] [Green Version]
  35. Garcia-Fragoso, L.; García-García, I.; de la Vega, A.; Renta, J.; Cadilla, C.L. Presence of the 5, 10-methylenetetrahydrofolate reductase C677T mutation in Puerto Rican patients with neural tube defects. J. Child Neurol. 2002, 17, 30–32. [Google Scholar] [CrossRef] [PubMed]
  36. Perez, A.B.A.; D’Almeida, V.; Vergani, N.; de Oliveira, A.C.; de Lima, F.T.; Brunoni, D. Methylenetetrahydrofolate reductase (MTHFR): Incidence of mutations C677T and A1298C in Brazilian population and its correlation with plasma homocysteine levels in spina bifida. Am. J. Med. Genet. Part A 2003, 119, 20–25. [Google Scholar] [CrossRef] [PubMed]
  37. Sadewa, A.H.; Sutomo, R.; Istiadjid, M.; Nishiyama, K.; Shirakawa, T.; Matsuo, M.; Nishio, H. C677T mutation in the MTHFR gene was not found in patients with frontoethmoidal encephalocele in East Java, Indonesia. Pediatrics Int. 2004, 46, 409–414. [Google Scholar] [CrossRef] [PubMed]
  38. Gao, C.M.; Tang, J.H.; Cao, H.X.; Ding, J.H.; Wu, J.Z.; Wang, J.; Liu, Y.T.; Li, S.P.; Su, P.; Matsuo, K.; et al. MTHFR polymorphisms, dietary folate intake and breast cancer risk in Chinese women. J. Hum. Genet. 2009, 54, 414–418. [Google Scholar] [CrossRef]
  39. Eser, B.; Cosar, M.; Eser, O.; Erdogan, M.O.; Aslan, A.; Yildiz, H.; Boyaci, G.; Buyukbas, S.; Solak, M. 677C>T and 1298A>C polymorphisms of methylenetetrahydropholate reductase gene and biochemical parameters in Turkish population with spina bifida occulta. Turk. Neurosurg. 2010, 20, 9–15. [Google Scholar] [PubMed]
  40. Harisha, P.N.; Devi, B.I.; Christopher, R.; Kruthika-Vinod, T.P. Impact of 5, 10-methylenetetrahydrofolate reductase gene polymorphism on neural tube defects. J. Neurosurg. Pediatrics 2010, 6, 364–367. [Google Scholar] [CrossRef] [PubMed]
  41. Liu, J.; Qi, J.; Yu, X.; Zhu, J.; Zhang, L.; Ning, Q.; Luo, X. Investigations of single nucleotide polymorphisms in folate pathway genes in Chinese families with neural tube defects. J. Neurol. Sci. 2014, 337, 61–66. [Google Scholar] [CrossRef] [PubMed]
  42. Nasri, K.; Midani, F.; Kallel, A.; Ben Jemaa, N.; Aloui, M.; Boulares, M.; Lassoued, M.; Ben Halima, M.; Ben Wafi, S.; Soussi, M.; et al. Association of MTHFR C677T, MTHFR A1298C, and MTRR A66G Polymorphisms with Neural Tube Defects in Tunisian Parents. Pathobiology 2019, 86, 190–200. [Google Scholar] [CrossRef] [PubMed]
  43. Zhang, Q.; Bai, B.-L.; Liu, X.-Z.; Miao, C.-Y.; Li, H. Association of folate metabolism genes MTRR and MTHFR with complex congenital abnormalities among Chinese population in Shanxi Province, China. Zhongguo Dang Dai Er Ke Za Zhi. 2014, 16, 840–845. [Google Scholar] [PubMed]
  44. Yu, J.; Chen, B.; Zhang, G.; Fu, S.; Li, P. The 677 C→T Mutation in the Methylenetetrahydrofolate Reductase (MTHFR) Gene in Five Chinese Ethnic Groups. Hum. Hered. 2000, 50, 268–270. [Google Scholar] [CrossRef] [PubMed]
  45. Pardo, R.; Suazo, J.; Castillo, S.; Vargas, M.; Zalavari, A.; Santos, J.L.; Blanco, R.; Rotter, K.; Solar, M.; Tapia, E. Methylenetetrahydrofolate reductase polymorphisms as risk factors for myelomeningocele. Rev. Med. Chil. 2014, 142, 587–592. [Google Scholar] [CrossRef]
  46. Wang, Y.; Liu, Y.; Ji, W.; Qin, H.; Wu, H.; Xu, D.; Turtuohut, T.; Wang, Z. Variants in MTHFR gene and neural tube defects susceptibility in China. Metab. Brain Dis. 2015, 30, 1017–1026. [Google Scholar] [CrossRef] [PubMed]
  47. Dutta, H.K.; Borbora, D.; Baruah, M.; Narain, K. Evidence of gene-gene interactions between MTHFD1 and MTHFR in relation to anterior encephalocele susceptibility in Northeast India. Birth Defects Res. 2017, 109, 432–444. [Google Scholar] [CrossRef]
  48. Fang, Y.; Zhang, R.; Zhi, X.; Zhao, L.; Cao, L.; Wang, Y.; Cai, C. Association of main folate metabolic pathway gene polymorphisms with neural tube defects in Han population of Northern China. Child’s Nerv. Syst. 2018, 34, 725–729. [Google Scholar] [CrossRef]
  49. Stegmann, K.; Zietgler, A.; Ngo, E.T.; Kohlschmidt, N.; Schröter, B.; Ermert, A.; Koch, M.C. Linkage disequilibrium of MTHFR genotypes 677C/T-1298A/C in the German population and association studies in probands with neural tube defects(NTD). Am. J. Med Genet. 1999, 87, 23–29. [Google Scholar] [CrossRef]
  50. Barber, R.; Shalat, S.; Hendricks, K.; Joggerst, B.; Larsen, R.; Suarez, L.; Finnell, R. Investigation of Folate Pathway Gene Polymorphisms and the Incidence of Neural Tube Defects in a Texas Hispanic Population. Mol. Genet. Metab. 2000, 70, 45–52. [Google Scholar] [CrossRef]
  51. A Volcik, K.; Blanton, S.H.; Tyerman, G.H.; Jong, S.T.; Rott, E.J.; Palge, T.Z.; Romaine, N.K.; Northrup, H. Methylenetetrahydrofolate reductase and spina bifida: Evaluation of level of defect and maternal genotypic risk in Hispanics. Am. J. Med Genet. 2000, 95, 21–27. [Google Scholar] [CrossRef]
  52. Cunha, A.L.; Hirata, M.H.; Kim, C.; Guerra-Shinohara, E.M.; Nonoyama, K.; Hirata, R.D.C. Metabolic effects of C677T and A1298C mutations at the MTHFR gene in Brazilian children with neural tube defects. Clin. Chim. Acta. 2002, 318, 139–143. [Google Scholar] [CrossRef]
  53. Félix, T.M.; Leistner, S.; Giugliani, R. Metabolic effects and the methylenetetrahydrofolate reductase (MTHFR) polymorphism associated with neural tube defects in southern Brazil. Birth Defects Res. Part A: Clin. Mol. Teratol. 2004, 70, 459–463. [Google Scholar] [CrossRef] [PubMed]
  54. Gos, M.; Sliwerska, E.; Szpecht-Potocka, A. Mutation incidence in folate metabolism genes and regulatory genes in Polish families with neural tube defects. J. Appl. Genet. 2004, 45, 363–368. [Google Scholar] [PubMed]
  55. Grandone, E.; Corrao, A.M.; Colaizzo, D.; Vecchione, G.; Di Girgenti, C.; Paladini, D.; Sardella, L.; Pellegrino, M.; Zelante, L.; Martinelli, P.; et al. Homocysteine metabolism in families from southern Italy with neural tube defects: Role of genetic and nutritional determinants. Prenat. Diagn. 2006, 26, 1–5. [Google Scholar] [CrossRef]
  56. Gonzalez-Herrera, L.; Castillo-Zapata, I.; Garcia-Escalante, G.; Pinto-Escalante, D. A1298C polymorphism of the MTHFR gene and neural tube defects in the state of Yucatan, Mexico. Birth Defects Res A Clin Mol Teratol. 2007, 79, 622–626. [Google Scholar] [CrossRef]
  57. Behunova, J.; Klimcakova, L.; Zavadilikova, E.; Potocekova, D.; Sykora, P.; Podracka, L. Methylenetetrahydrofolate reductase gene polymorphisms and neural tube defects epidemiology in the Slovak population. Birth Defects Res. Part A: Clin. Mol. Teratol. 2010, 88, 695–700. [Google Scholar] [CrossRef]
  58. Selvi, N.; Ergun, R.; Gube, O.; Yİldİrİ, U.; Tetİk, A.; Gunduz, C.; Ulman, I.; Eroglu, Z. Genetic Polymorphism of Methylenetetrahydrofolate Reductase. Middle-East Journal of Scientific Research. 2010, 6, 93–98. [Google Scholar]
  59. Yildiz, S.H.; Erdogan, M.O.; Solak, M.; Eser, O.; Terzi, E.A.; Eser, B.; Kocabaş, V.; Aslan, A. Lack of association between the methylenetetrahydropholate reductase gene A1298C polymorphism and neural tube defects in a Turkish study group. Genet Mol Res. 2016, 15. [Google Scholar] [CrossRef]
  60. Morrison, K.; Papapetrou, C.; Hol, F.A.; Mariman, E.C.M.; Lynch, S.A.; Burn, J.; Edwards, Y.H. Susceptibility to spina bifida; an association study of five candidate genes. Ann. Hum. Genet. 1998, 62, 379–396. [Google Scholar] [CrossRef]
  61. Christensen, B.; Arbour, L.; Tran, P.; Leclerc, D.; Sabbaghian, N.; Platt, R.; Gilfix, B.M.; Rosenblatt, D.S.; Gravel, R.A.; Forbes, P.; et al. 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. [Google Scholar] [CrossRef]
  62. Johanning, G.L.; Tamura, T.; Johnston, K.E.; Wenstrom, K.D. Comorbidity of 5,10-methylenetetrahydrofolate reductase and methionine synthase gene polymorphisms and risk for neural tube defects. J. Med Genet. 2000, 37, 949–951. [Google Scholar] [CrossRef]
  63. O’Leary, V.B.; Mills, J.L.; Parle-McDermott, A.; Pangilinan, F.; Molloy, A.M.; Cox, C.; Weiler, A.; Conley, M.; Kirke, P.N.; Scott, J.M.; et al. Screening for new MTHFR polymorphisms and NTD risk. Am. J. Med. Genet. Part A. 2005, 138, 99–106. [Google Scholar] [CrossRef] [PubMed]
  64. Doudney, K.; Grinham, J.; Whittaker, J.; Lynch, S.; Thompson, D.; Moore, G.; Copp, A.L.; Greene, N.; Stanier, P. Evaluation of folate metabolism gene polymorphisms as risk factors for open and closed neural tube defects. Am. J. Med Genet. Part A. 2009, 149A, 1585–1589. [Google Scholar] [CrossRef] [PubMed]
  65. Wilson, A.; Platt, R.; Wu, Q.; Leclerc, D.; Christensen, B.; Yang, H.; Gravel, R.A.; Rozen, R. A Common Variant in Methionine Synthase Reductase Combined with Low Cobalamin (Vitamin B12) Increases Risk for Spina Bifida. Mol. Genet. Metab. 1999, 67, 317–323. [Google Scholar] [CrossRef]
  66. Pietrzyk, J.J.; Bik-Multanowski, M.; Salnak, M.; Twardowska, M. Polymorphisms of the 5,10-methylenetetrahydrofolate and the methionine synthase reductase genes as independent risk factors for spina bifida. J. Appl. Genet. 2003, 44, 111–113. [Google Scholar] [PubMed]
  67. Pei, L.; Zhu, H.; Retn, A.; Li, Z.; Hao, L.; Finnell, R.; Li, Z. Reduced folate carrier gene is a risk factor for neural tube defects in a Chinese population. Birth Defects Res. Part A: Clin. Mol. Teratol. 2005, 73, 430–433. [Google Scholar] [CrossRef] [PubMed]
  68. Cao, L.; Wang, Y.; Zhang, R.; Dong, L.; Cui, H.; Fang, Y.; Zhao, L.; Shi, O.; Cai, C. Association of neural tube defects with gene polymorphisms in one-carbon metabolic pathway. Child’s Nerv. Syst. 2018, 34, 277–284. [Google Scholar] [CrossRef]
  69. Kotova, N.; Maichuk, V.; Fedorenko, O. The differentiated approach to prevention of neural tube defects in children. Georgian Med. News. 2018, 274, 52–59. [Google Scholar]
  70. Crider, K.S.; Yang, T.P.; Berry, R.J.; Bailey, L.B. Folate and DNA methylation: A review of molecular mechanisms and the evidence for folate’s role. Adv. Nutr. 2012, 3, 21–38. [Google Scholar] [CrossRef]
  71. Zheng, J.; Lu, X.; Liu, H.; Zhao, P.; Li, K.; Li, L. MTHFD1 polymorphism as maternal risk for neural tube defects: A meta-analysis. Neurol. Sci. 2015, 36, 607–616. [Google Scholar] [CrossRef]
  72. Friso, S.; Choi, S.-W.; Girelli, D.; Mason, J.B.; Dolnikowski, G.G.; Bagley, P.J.; Olivieri, O.; Jacques, P.F.; Rosenberg, I.H.; Corrocher, R. A common mutation in the 5, 10-methylenetetrahydrofolate reductase gene affects genomic DNA methylation through an interaction with folate status. Proc. Natl. Acad. Sci. USA 2002, 99, 5606–5611. [Google Scholar] [CrossRef] [PubMed]
  73. Kondo, A.; Kamihira, O.; Ozawa, H. Neural tube defects: Prevalence, etiology and prevention. Int. J. Urol. 2009, 16, 49–57. [Google Scholar] [CrossRef] [PubMed]
  74. Frosst, P.; Blom, H.; Milos, R.; Goyette, P.; Sheppard, C.A.; Matthews, R.; Boers, G.; Den Heijer, M.; Kluijtmans, L.; Van Den Heuve, L. A candidate genetic risk factor for vascular disease: A common mutation in methylenetetrahydrofolate reductase. Nat. Genet. 1995, 10, 111–113. [Google Scholar] [CrossRef]
  75. Weisberg, I.; Tran, P.; Christensen, B.; Sibani, S.; Rozen, R. A second genetic polymorphism in methylenetetrahydrofolate reductase (MTHFR) associated with decreased enzyme activity. Mol. Genet. Metab. 1998, 64, 169–172. [Google Scholar] [CrossRef]
  76. Li, Y.N.; Gulati, S.; Baker, P.J.; Brody, L.C.; Banerjee, R.; Kruger, W.D. Cloning, mapping and RNA analysis of the human methionine synthase gene. Hum. Mol. Genet. 1996, 5, 1851–1858. [Google Scholar] [CrossRef] [Green Version]
  77. Selzer, R.R.; Rosenblatt, D.S.; Laxova, R.; Hogan, K. Adverse effect of nitrous oxide in a child with 5, 10-methylenetetrahydrofolate reductase deficiency. N. Engl. J. Med. 2003, 349, 45–50. [Google Scholar] [CrossRef]
  78. Sutherland, H.G.; Hermile, H.; Sanche, R.; Menon, S.; Lea, R.A.; Haupt, L.M.; Griffiths, L.R. Association study of mthfd 1 coding polymorphisms r 134 k and r 653 q with migraine susceptibility. Headache J. Head Face Pain. 2014, 54, 1506–1514. [Google Scholar] [CrossRef] [PubMed]
  79. Christensen, K.E.; Dahhou, M.; Kramer, M.S.; Rozen, R. The MTHFD1 1958G>A variant is associated with elevated C-reactive protein and body mass index in Canadian women from a premature birth cohort. Mol. Genet. Metab. 2014, 111, 390–392. [Google Scholar] [CrossRef] [PubMed]
  80. Cui, Y.; Jing, Y.; Sun, Z. Lack of association between MTHFD1 G401A polymorphism and ovarian cancer susceptibility. Tumor Biol. 2014, 35, 3385–3389. [Google Scholar] [CrossRef]
  81. Da Silva, L.M.; Galbiatti, A.L.; Ruiz, M.T.; Raposo, L.S.; Maniglia, J.V.; Pavarino, E.C.; Goloni-Bertollo, E.M. MTHFD1 G1958A, BHMT G742A, TC2 C776G and TC2 A67G polymorphisms and head and neck squamous cell carcinoma risk. Mol. Biol. Rep. 2012, 39, 887–893. [Google Scholar] [CrossRef]
  82. García-González, I.; Flores-Ocampo, J.; Alejos-Mex, R.; Herrera-Sánchez, L.; González-Herrera, L. G894T (NOS3) and G1958A (MTHFD1) gene polymorphisms and risk of ischemic heart disease in Yucatan, Mexico. Clin. Investig. Arterioscler. Publ. Of. Soc. Esp. Arterioscler. 2014, 27, 64–73. [Google Scholar]
  83. Jiang, J.; Zhang, Y.; Wei, L.; Sun, Z.; Liu, Z. Association between MTHFD1 G1958A polymorphism and neural tube defects susceptibility: A meta-analysis. PLoS ONE 2014, 9, e101169. [Google Scholar] [CrossRef]
  84. Keller, M.D.; Ganesh, J.; Heltzer, M.; Paessler, M.; Bergqvist, A.C.; Baluarte, H.J.; Watkins, D.; Rosenblatt, D.S.; Orange, J.S. Severe combined immunodeficiency resulting from mutations in MTHFD1. Pediatrics 2013, 131, e629–e634. [Google Scholar] [CrossRef] [PubMed]
  85. Lorenc, A.; Seremak-Mrozikiewicz, A.; Barlik, M.; Wolski, H.; Drews, K. The role of 401a>G polymorphism of methylenetetrahydrofolate dehydrogenase gene (MTHFD1) in fetal hypotrophy. Ginekol. Pol. 2014, 85, 494–499. [Google Scholar] [CrossRef]
  86. Murthy, J.; Gurramkonda, V.B.; Lakkakula, B.V. Significant association of MTHFD1 1958G>A single nucleotide polymorphism with nonsyndromic cleft lip and palate in Indian population. Med. Oral Patol. Oral Y Cir. Bucal 2014, 19, e616–e621. [Google Scholar] [CrossRef]
  87. Carroll, N.; Pangilinan, F.; Molloy, A.M.; Troendle, J.; Mills, J.L.; Kirke, P.N.; Brody, L.C.; Scott, J.M.; Parle-McDermott, A. Analysis of the MTHFD1 promoter and risk of neural tube defects. Hum. Genet. 2009, 125, 247–256. [Google Scholar] [CrossRef] [Green Version]
  88. Parle-McDermott, A.; Pangilinan, F.; O’Brien, K.K.; Mills, J.L.; Magee, A.M.; Troendle, J.; Sutton, M.; Scott, J.M.; Kirke, P.N.; Molloy, A.M.; et al. A common variant in MTHFD1L is associated with neural tube defects and mRNA splicing efficiency. Hum. Mutat. 2009, 30, 1650–1656. [Google Scholar] [CrossRef]
  89. Green, R.F.; Byrne, J.; Crider, K.S.; Gallagher, M.; Koontz, D.; Berry, R.J. Folate-related gene variants in Irish families affected by neural tube defects. Front. Genet. 2013, 4, 223. [Google Scholar] [CrossRef] [Green Version]
Jpm 12 01609 g001 550
Figure 1. Homozygous model for MTHFR C677T. Van Der Put [34], De Franchis [30], Bodurofulo [23], Garcia-Fragoso [35], Perez [36], Relton [28], Sadewa [37], Kirke [24], Gao [38], Eser [39], Harisha [40], Wei [41], Liu [42], Qin [43], Yu [44], Pardo [45], Wang [46], Dutta [47], Fang [48].
Figure 1. Homozygous model for MTHFR C677T. Van Der Put [34], De Franchis [30], Bodurofulo [23], Garcia-Fragoso [35], Perez [36], Relton [28], Sadewa [37], Kirke [24], Gao [38], Eser [39], Harisha [40], Wei [41], Liu [42], Qin [43], Yu [44], Pardo [45], Wang [46], Dutta [47], Fang [48].
Jpm 12 01609 g001
Jpm 12 01609 g002 550
Figure 2. Homozygous model for MTHFR A1298C. Van der Put [34], Stegmann [49], Barber [50], Volcik [51], Richter [19], Cunha [52], De Marco [25], Parle-McDermott [26], Perez [36], Felix [53], Gos [54], Relton [28], Sadewa [37], Grandone [55], Gonzalez-herrera [56], Behunova [57], Eser [39], Selvi [58], Pardo [45], Wang [46], Yildiz [59], Dutta [47].
Figure 2. Homozygous model for MTHFR A1298C. Van der Put [34], Stegmann [49], Barber [50], Volcik [51], Richter [19], Cunha [52], De Marco [25], Parle-McDermott [26], Perez [36], Felix [53], Gos [54], Relton [28], Sadewa [37], Grandone [55], Gonzalez-herrera [56], Behunova [57], Eser [39], Selvi [58], Pardo [45], Wang [46], Yildiz [59], Dutta [47].
Jpm 12 01609 g002
Jpm 12 01609 g003 550
Figure 3. Recessive mode for MTR A2756G. Van Der Put [34], Morrison [60], Christensen [61], Shaw [27], Johanning [62], De Marco [25], Gos [54], O’Leary [63], Doudney [64].
Figure 3. Recessive mode for MTR A2756G. Van Der Put [34], Morrison [60], Christensen [61], Shaw [27], Johanning [62], De Marco [25], Gos [54], O’Leary [63], Doudney [64].
Jpm 12 01609 g003
Jpm 12 01609 g004 550
Figure 4. Association between MTRR A66G. Wilson [65], Pietrzyk [66], Gos [54], Relton [28], O’Leary [63], Van Der Linden [21], Doudney [64].
Figure 4. Association between MTRR A66G. Wilson [65], Pietrzyk [66], Gos [54], Relton [28], O’Leary [63], Van Der Linden [21], Doudney [64].
Jpm 12 01609 g004
Jpm 12 01609 g005 550
Figure 5. (A) Allele model for RFC; Marco [25], Relton [28], Pei [67], Cao [68](B) Dominant model for RFC; (C) Homozygote model for RFC; (D) Heterozygous model for RFC.
Figure 5. (A) Allele model for RFC; Marco [25], Relton [28], Pei [67], Cao [68](B) Dominant model for RFC; (C) Homozygote model for RFC; (D) Heterozygous model for RFC.
Jpm 12 01609 g005
Jpm 12 01609 g006 550
Figure 6. The folate pathway with significant compounds (red boxes), enzymes (blue circles) and transporters (green triangles) outlined. Single nucleotide polymorphisms studied in this meta-analysis are provided in boxes below the enzymes and transporters with the ones found to be significantly associated with NTD risk highlighted in red. MTHFR: Methylenetetrahydrofolate reductase; MTHFD: Methylenetetrahydrofolate dehydrogenase; MTR: Methionine synthase; MTRR: Methionine synthase reductase; THF: Tetrahydrofolate; TYMS: thymidylate synthase; DHFR: Dihydrofolate Reductase; DHF: Dihydrofolate; RFC: Reduced Folate Carrier; FOLH1: Folate Hydrolase 1; SHMT1: Serine Hydroxymethyltransferase; BHMT: Betaine—homocysteine S-methyltransferase; SAH: S-adenosylhomocysteine; SAM, S-adenosylmethionine; AHCY, S-adeonsyl-L-homocysteine hydrolase; MAT, methionine adenosyltransferase; Hcy: Homocysteine; CBS: Cystathione ß-synthase; CTH: Cystathionase; GS: Glutathione synthase; PCFR: Protein Coupled Folate Receptor; MT: Methyltransferase.
Figure 6. The folate pathway with significant compounds (red boxes), enzymes (blue circles) and transporters (green triangles) outlined. Single nucleotide polymorphisms studied in this meta-analysis are provided in boxes below the enzymes and transporters with the ones found to be significantly associated with NTD risk highlighted in red. MTHFR: Methylenetetrahydrofolate reductase; MTHFD: Methylenetetrahydrofolate dehydrogenase; MTR: Methionine synthase; MTRR: Methionine synthase reductase; THF: Tetrahydrofolate; TYMS: thymidylate synthase; DHFR: Dihydrofolate Reductase; DHF: Dihydrofolate; RFC: Reduced Folate Carrier; FOLH1: Folate Hydrolase 1; SHMT1: Serine Hydroxymethyltransferase; BHMT: Betaine—homocysteine S-methyltransferase; SAH: S-adenosylhomocysteine; SAM, S-adenosylmethionine; AHCY, S-adeonsyl-L-homocysteine hydrolase; MAT, methionine adenosyltransferase; Hcy: Homocysteine; CBS: Cystathione ß-synthase; CTH: Cystathionase; GS: Glutathione synthase; PCFR: Protein Coupled Folate Receptor; MT: Methyltransferase.
Jpm 12 01609 g006
Table
Table 1. Subgroup analysis using 5 genetic models. Model Code: R = Random, F = Fixed.
Table 1. Subgroup analysis using 5 genetic models. Model Code: R = Random, F = Fixed.
MTHFR C677T
ModelEthnicityNumber of Studies Test of
Association		 Test of
Heterogeneity		 Bias
OR95% CIp-ValueModelp-ValueI^2p
Allele contrast
(A vs. a)		Overall190.722 [0.6337; 0.8227]0.000001Random0.00220.550.75
Asian100.629 [0.5021; 0.7878]0.00005Random0.00400.630.33
Caucasian70.7676 [0.6462; 0.9118]0.003Random0.08250.460.47
Mixed20.9273 [0.6956; 1.2363]0.61Fixed0.5556
Recessive model
(AA vs. Aa + aa)		Overall190.682 [0.5759; 0.8077]0.000009Random0.03020.420.51
Asian100.5676 [0.4623; 0.6970]0.00000006Fixed0.23460.230.74
Caucasian70.7817 [0.6053; 1.0094]0.06Random0.05170.520.72
Mixed20.8206 [0.5697; 1.1819]0.29Fixed0.3387
Dominant model (AA + Aa vs. aa)Overall180.6441 [0.5219; 0.7948]0.00004Random0.01770.460.87
Asian90.5954 [0.4205; 0.8430]0.003Random0.00690.620.23
Caucasian70.5727 [0.4648; 0.7057]0.0000001Fixed0.57670.000.28
Mixed21.3088 [0.6597; 2.5965]0.44Fixed0.8181
Homozygote
(AA vs. aa)		Overall180.5519 [0.4380; 0.6955]0.0000005Random0.05390.380.89
Asian90.4627 [0.3204; 0.6681]0.00004Random0.07740.440.25
Caucasian70.5443 [0.4334; 0.6837]0.0000002Fixed0.20660.290.52
Mixed21.1424 [0.5639; 2.3145]0.71Fixed0.9638
Heterozygous
(AA vs. Aa)		Overall190.7328 [0.6141; 0.8745]0.0006Random0.04420.390.46
Asian100.6034 [0.4824; 0.7546]0.00001Fixed0.31250.140.80
Caucasian70.8635 [0.6633; 1.1240]0.28Random0.06750.490.93
Mixed20.7676 [0.5231; 1.1265]0.18Fixed0.26620.19
RFC A80G
Allele contrast
(A vs. a)		Overall40.7729 [0.6607; 0.9043]0.001Fixed0.32300.140.53
Asian20.6468 [0.5053; 0.8280]0.0005Fixed0.72030.00
Caucasian20.8721 [0.7117; 1.0688]0.19Fixed0.98750.00
Recessive model (AA vs. Aa + aa)Overall40.8223 [0.6382; 1.0596]0.13Fixed0.46880.000.15
Asian20.6321 [0.4062; 0.9838]0.04Fixed0.61880.00
Caucasian20.9352 [0.6864; 1.2743]0.67Fixed0.60640.00
Dominant model (AA + Aa vs. aa)Overall40.6248 [0.4858; 0.8036]0.0002Fixed0.57970.000.76
Asian20.5304 [0.3676; 0.7653]0.0007Fixed0.87620.00
Caucasian20.7229 [0.5115; 1.0217]0.07Fixed0.48320.00
Homozygote (AA vs. aa)Overall40.6028 [0.4410; 0.8238]0.002Fixed0.51590.000.55
Asian20.4536 [0.2756; 0.7467]0.002Fixed0.65430.00
Caucasian20.7246 [0.4852; 1.0821]0.12Fixed0.87980.00
Heterozygous (AA vs. Aa)Overall40.9533 [0.7268; 1.2504]0.73Fixed0.67690.000.03
Asian20.8028 [0.5006; 1.2876]0.36Fixed0.63840.00
Caucasian21.0374 [0.7447; 1.4451]0.83Fixed0.46050.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Almekkawi, A.K.; AlJardali, M.W.; Daadaa, H.M.; Lane, A.L.; Worner, A.R.; Karim, M.A.; Scheck, A.C.; Frye, R.E. Folate Pathway Gene Single Nucleotide Polymorphisms and Neural Tube Defects: A Systematic Review and Meta-Analysis. J. Pers. Med. 2022, 12, 1609. https://doi.org/10.3390/jpm12101609

AMA Style

Almekkawi AK, AlJardali MW, Daadaa HM, Lane AL, Worner AR, Karim MA, Scheck AC, Frye RE. Folate Pathway Gene Single Nucleotide Polymorphisms and Neural Tube Defects: A Systematic Review and Meta-Analysis. Journal of Personalized Medicine. 2022; 12(10):1609. https://doi.org/10.3390/jpm12101609

Chicago/Turabian Style

Almekkawi, Ahmad K., Marwa W. AlJardali, Hicham M. Daadaa, Alison L. Lane, Ashley R. Worner, Mohammad A. Karim, Adrienne C. Scheck, and Richard E. Frye. 2022. "Folate Pathway Gene Single Nucleotide Polymorphisms and Neural Tube Defects: A Systematic Review and Meta-Analysis" Journal of Personalized Medicine 12, no. 10: 1609. https://doi.org/10.3390/jpm12101609

APA Style

Almekkawi, A. K., AlJardali, M. W., Daadaa, H. M., Lane, A. L., Worner, A. R., Karim, M. A., Scheck, A. C., & Frye, R. E. (2022). Folate Pathway Gene Single Nucleotide Polymorphisms and Neural Tube Defects: A Systematic Review and Meta-Analysis. Journal of Personalized Medicine, 12(10), 1609. https://doi.org/10.3390/jpm12101609

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop