ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Systematic Review

Is MTHFD1 polymorphisms rs 2236225 (c.1958G>A) associated with the susceptibility of NSCL/P? A systematic review and meta-analysis

[version 1; peer review: 1 approved, 1 approved with reservations]
PUBLISHED 04 Jun 2015
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

Aims: To investigate the association between the methylenetetrahydrofolate dehydrogenase 1 (MTHFD1) polymorphism rs 2236225 (c.1958G>A) and susceptibility to non-syndromic cleft of the lip and/or palate (NSCL/P).
Methods: An extensive literature review has been conducted using PubMed, Web of Science, Cochrane Library, Google Scholar, the China National Knowledge Infrastructure (CNKI), and Wanfang Database for eligible researches. The terms for searching were “cleft lip OR cleft palate OR CLP OR CL/P OR oral facial cleft OR OFC” AND “methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1 OR methenyltetrahydrofolate cyclohydrolase formyltetrahydrofolate synthetase OR MTHFD1 OR MTHFD”. Two independent researchers screened, evaluated and extracted the data of included studies. The pooled odds ratios (OR) with 95% confidence intervals (95% CI) were calculated by random effects model under five gene models. Subgroup, sensitivity analysis and publication bias were also assessed.
Results: Ten case-control studies have been included in the systematic review and eight studies have been considered for the meta-analysis. Overall, the MTHFD1 polymorphism rs2236225 and the risk of NSCL/P showed pooled OR (95% CI) of 1.02 (0.86-1.21) under allelic model. A higher degree of heterogeneity was observed in Asian countries (I2 = 75.6%) compared to non-Asian countries (I2 = 48.9%). Similar consequence appeared in the subgroup of children (I2 = 78.6%) compared with that of mothers (I2 = 0.0%). There was no significant difference in the publication bias by the Begg’s funnel plot (P = 0.711) and Egger’s regression test (P = 0.746).
Conclusion: Our assessment suggested there was no significant association between the MTHFD1 polymorphism rs 2236225 (c.1958G>A) and the susceptibility to NSCL/P. Further investigations using a large sample size and a more advanced technique should be adopted to reach a more precise conclusion in the future.

Keywords

MTHFD1, Polymorphisms, NSCL/P susceptibility, Meta-analysis

Introduction

Cleft of the lip and/or palate (CL/P) is one of the most common facial malformations13 and a societal burden, affecting the patient ability to eat and speak and influencing social integration4. Non-syndromic CL/P, accounting for about 70% of CL/P, is considered closely related to genetic and environmental factors5. Recent studies suggested that using folic acid could reduce the rates of oral clefts6,7 and single nucleotide polymorphisms of some genes such as MTHFR8,9, MTR40 and MTRR involved in the metabolism of folic acid have been associated to high risk of NSCL/P8,9. Methylenetetrahydrofolate dehydrogenase 1 (MTHFD1), a key gene associated with three sequential enzymatic reactions in the metabolism of folic acid, might play a potential role in the risk of NSCL/P, especially the polymorphism rs2236225 (c.1958G>A)10. Indeed, different observations that linked the polymorphism rs2236225 to the risk of NSCL/P have been reported11,12. The suggestion of a link between rs2236225 polymorphism and susceptibility to NSCL/P might be result of the limitations in sample size, different ethnic populations and other environmental factors. Therefore, we conducted a systematic review and meta-analysis of eligible case-control studies to reveal a more precise connection between the MTHFD1 polymorphism rs2236225 and the risk of NSCL/P.

Materials and methods

Identification of studies

A systematic search based on the principle of evidence-based medicine13 was performed in PubMed, Web of Science, Cochrane Library, Google Scholar, China National Knowledge Infrastructure (CNKI) and WanFang Database. The final update was made on April 5th, 2015. In line with our knowledge background, the Medical Subject Headings (MESH) terms in PubMed and the known aliases of the genes of interests in the National Center of Biotechnology Information (NCBI), the following terms were used for searching: “cleft lip OR cleft palate OR CLP OR CL/P OR oral facial cleft OR OFC” AND “methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1 OR methenyltetrahydrofolate cyclohydrolase formyltetrahydrofolate synthetase OR MTHFD1 OR MTHFD”, which were slightly adjusted to optimize search results (Table S1; PubMed). We didn’t limit the search depending on publication types, data and language. Of course, the review of the published literature was examined carefully and manual search was conducted to avoid missing potential data. Two of the authors (Huaxiang Zhao and Mengqi Zhang) were in charge of the search independently and a third author (Jieni Zhang) conducted a random inspection.

Inclusion and exclusion criteria

Researches included in our systematic review and meta-analysis meet the following criteria: (1) evaluating the association between the NSCL/P and MTHFD1 polymorphism rs2236225, (2) focusing on humans, (3) case-control studies. Exclusion criteria were: (1) no association between NSCL/P and MTFHD1, (2) not focusing on humans but animal models or in vitro studies, (3) duplication of previous researches, (4) not original literature such as reviews, meta-analyses, comments and editorials.

Data collection

Data from eligible studies were extracted by two independent researchers (Huaxiang Zhao and Mengqi Zhang) in accordance with the inclusion and exclusion criteria. In case of any discrepancies, the third chief author (Feng Chen) would make a further investigation or bring it into a group-discussion. A special table was used for collecting information from the selected articles and the following entries were recorded: authors (year), country, location of geography, subjects, methods for genotyping, sample size of cases/controls, descriptions of samples rolled in the study, P for HWE (Hardy-Weinberg equilibrium) of control group, whether included in the meta-analysis or not.

Methodological quality assessment

A methodological quality assessment adapted from previous studies1416 was carried on included studies (Table S2). Cases, source of controls, sample sizes and Hardy-Weinberg equilibrium (HWE) were considered as important aspects in this systematic review.

Statistical analysis

The PRISMA checklist (Supplementary material S3) was used as a protocol in our meta-analysis17. Odd ratios (ORs) and 95% confidence intervals (CIs) were calculated to estimate the association between the susceptibility to NSCL/P and MTHFD1. Five genetic models were used in the process of pooling the OR and 95% CIs: allelic comparison (A versus G), heterozygote model (AG versus GG), homozygote model (AA versus GG), dominant model (AA + AG versus GG), recessive model (AA versus AG + GG). The significance of the pooled effects was determined by Z-test with P value less than 0.05. The Q-statistic and the I2 test were used to evaluated; P < 0.05 in Q statistic or I2 > 50%18,19, would indicate a significant heterogeneity. When P > 0.05 in Q statistic or I2 < 50%, the fixed pooling model (Mantel-Haenszel) was conducted; if not, the random pooling model (M-H heterogeneity) was used. We also carried subgroup analyses in which different subjects (mothers or children), location of geography (non-Asian countries or Asian countries) were considered potential source of heterogeneity. A sensitivity analysis was conducted by omitting each study in turn to evaluate the single study’s influence on the overall estimation. We used Begg’s funnel plot and Egger’s linear regression test to find out the publication bias of the included studies2022. The studies with disequilibrium of HWE among control group were added into a supplementary meta-analysis as described previously23. Meanwhile, as for the studies included but not carried into the meta-analysis, to achieve a qualitative analysis we adopted a method described by others24. Results were considered significant when P < 0.05. Stata 12.0 (Stata Corp, College Station, TX, USA) was used for the analysis.

Results

Data retrieval

A total of 251 articles resulted from the search described above (PubMed: 86, Web of Science: 8, Google Scholar: 135, Cochrane Library: 0, CNKI: 18, Wanfang: 4). After being imported into EndNote X6 software (Thomson Corporation, Stamford), a screening process was conducted among 102 articles– that is, duplicates were removed using the ‘Discard Duplicates’ function as well as by handwork. Following paper selection by two independent researchers, 15 studies were then thoroughly reviewed. Of these, five studies were excluded, among which two had no control groups25,26, one no relation to MTFHD127, and the other two presented data previously published28,29. Finally, 10 studies that met the criteria were included in the systematic review (Table 1)1012,3036 and mathematic data from eight studies were used for reference to carry out the meta-analysis1012,3133,35,36. The selection process is shown in Figure 1.

Table 1. Characteristics of studies included in the systematic review and meta-analysis.

No.Authors (year)CountryGeographical
location
SubjectsMethods for
genotyping
Sample size of
case/control
group (just for
the patients)
Descriptions of
samples from study
participants
P for
HWE* of
control
group
Whether
included
in meta-
analysis
or not
casecontrol
1Mostowska et al.
(2006)
PolandEuropeMothersPCR-RFLPγ12282The case samples
came from healthy
mothers of NSCL/P
children, while the
control group includes
samples from healthy
mothers of children
without NSCL/P. There
was no difference
between the two
groups in their age,
habit of smoking.
NMψYes
2Boyles et al.
(2008)
NorwayEuropeMothers
and
children
MALDI-TOF
MSξ
573763377 cases were CL/P
and 196 cases CPO.
Most mothers in
the case group use
supplemental folate
during the pregnancy.
NMψNo
3Mills et al.
(2008)
IrelandEuropeMothers,
fathers
and
children
PCR-RFLPγ10301000536 were CLP
consisted of 494 cases
with isolated defects
23 with one additional
defect, 18 with multiple
defects, and one with
Pierre Robin. 426 cases
with CPO consisted of
321 isolated defects,
15 with one additional
defect, 21 with multiple
defects, and 69 with
Pierre Robin Sequence.
0.03Yes
4Bufalino et al.
(2010)
BrazilSouth AmericaMothersPCR-RFLPγ106184Mothers who smoke,
drink and use anti-
hypertensives and
drugs that could
potentially impair the
function of folic acids
were not included in
this study.
0.66Yes
5Mostowska et al.
(2010)
PolandEuropeChildrenPCR-RFLPγ174176The patients with clefts
palate only (CPO) were
excluded because the
researchers thought
the pathogenesis of
NSCL/P and the CPO
was different.
0.11Yes
6Li et al. (2013)ChinaAsianChildrenPCR-RFLPγ187157The patients in the case
group consisted of 126
boys and 61 girls.
0.89Yes
7Yuan (2013)ChinaAsianMothers,
fathers
and
children
PCR-RFLPγ15015068 CLO and 82 CLP
were enrolled in the
case group.
0.92Yes
8Zhao et al.
(2013)
ChinaAsianChildrenPCR-RFLPγ294126There were 191 CLP
and 103 CPO in the
patients group.
0.08Yes
9de Aquino et al.
(2013)
BrazilSouth AmericaMothers,
fathers
and
children
Real-Time
PCR
181478Patients with clefts
palate only (CPO) were
excluded. 65 clefts lip
only (CLO) and 116
clefts lip and palate
(CLP) were included in
this study, consisting
of 101 males and 80
females.
NMψNo
10Murthy et al.
(2014)
IndiaAsianChildrenPCR-RFLPγ142141There were 123 CLP
and 19 CPO in the case
group.
0.94Yes

HWE*: Hardy-Weinberg equilibrium.

NMψ: Not mentioned in the study.

PCR-RFLPγ: PCR-restriction fragment length polymorphism.

MALDI-TOF MSξ: matrix-assisted laser desorption/ionization time-of-flight mass spectrometry.

7382c295-0b6c-4c1a-9201-70ddb38f024c_figure1.gif

Figure 1. Flow chart showing study selection in the systematic and meta-analysis.

Study characteristics

Eventually, all 10 studies containing 6216 samples (2959 cases and 3257 controls) were analyzed in our review. The characteristics of every study can be seen in Table 1. To summarize briefly, there were four studies from European groups, four from Asian groups and two from South American groups, among which two studies focused on the genotype of patients’ mothers only, four on children’ s genotype only and four on both of them. PCR-restriction fragment length polymorphism (PCR-RFLP) was the major method of genotyping, while other techniques had been used as well.

Association between MTHFD1 polymorphism rs2236225 (c.1958G>A) and NSCL/P susceptibility

The association between MTHFD1 polymorphism rs2236225 (c.1958G>A) and NSCL/P susceptibility was analyzed through a meta-analysis and qualitative analysis. In the meta-analysis, since significant heterogeneity had been identified by Q-test and I2 statistic in every genetic model, the random effect models were used. Overall, a significant association was not found in any genetic model (A versus G: OR = 1.02, 95% CI 0.86–1.21, PH = 0.010, Figure 2; AG versus GG: OR = 0.97, 95% CI 0.75–1.26, PH = 0.019, Figure 3A; AA versus GG: OR = 1.07, 95% CI 0.70–1.65, PH = 0.005, Figure 3B; AA + AG versus GG: OR = 1.00, 95% CI 0.76–1.31, PH = 0.006, Figure 3C; AA versus AG + GG: OR = 1.05, 95% CI 0.71–1.53, PH = 0.014, Figure 3D). On the other hand, no association was found in the genotypes of children, mothers or fathers in the qualitative analysis30,34.

7382c295-0b6c-4c1a-9201-70ddb38f024c_figure2.gif

Figure 2. Forest plot of allelic comparison of MTHFD1 polymorphism rs2236225 (c.1958G>A) for overall comparison (A versus G).

7382c295-0b6c-4c1a-9201-70ddb38f024c_figure3.gif

Figure 3. Forest plot of heterozygote, homozygote, dominant and recessive model comparison of MTHFD1 polymorphisms rs2236225 (c.1958G>A) for overall comparison.

(A) Heterozygote model, AG versus GG. (B) Homozygote model, AA versus GG. (C) Dominant model, AA + AG versus GG. (D) Recessive model, AA versus AG + GG.

Next we conducted the subgroup analysis using allelic A versus G model according to the location of geography and subjects (mothers or children). It turned out that there was no significant difference between Asian (OR = 1.03, 95% CI 0.75–1.40, PH = 0.003) or non-Asian population (OR = 1.06, 95% CI 0.86–1.30, PH = 0.118). However, a higher degree of heterogeneity was observed in the Asian countries compared to non-Asian countries (Figure 4A). A similar result was observed in the subgroup analysis between mothers and children. The heterogeneity was much higher in the children group (OR = 0.99, 95% CI 0.72–1.36, PH = 0.001) than in the mothers’ group (OR = 1.11, 95% CI 0.98–1.27, PH = 0.630), while no significant difference was observed in both groups (Figure 4B).

7382c295-0b6c-4c1a-9201-70ddb38f024c_figure4.gif

Figure 4.

Subgroup analysis by locations of geography (A) and subjects (B) under allelic comparison of MTHFD1 polymorphism rs2236225 (c.1958G>A).

Sensitivity analysis and publication bias

To access the influence of each individual study on the pooled ORs, a sensitivity analysis was performed by omitting each study at a time. The results of sensitivity suggests that no individual study affects the pooled ORs of the associations between MTHFD1 polymorphism rs2236225 (c.1958G>A) and NSCL/P risk under allelic model (Figure 5).

7382c295-0b6c-4c1a-9201-70ddb38f024c_figure5.gif

Figure 5. Sensitivity analysis of the association between MTHFD1 polymorphism rs2236225 (c.1958G>A) and susceptibility to NSCL/P under allelic model (A versus G).

We used the Begg’s funnel plot and Egger’s regression test (both used the allelic A versus G model) to estimate the publication bias. Our results indicate that there is no significant publication bias both in the symmetry of Begg’s funnel plot (P = 0.711, Figure 6) and Egger’s regression test (P = 0.746).

7382c295-0b6c-4c1a-9201-70ddb38f024c_figure6.gif

Figure 6. Begg’s funnel plot of the association between MTHFD1 polymorphism rs2236225 (c.1958G>A) and the susceptibility to NSCL/P under allelic model (A versus G).

Discussion

CL/P is one of the most common facial malformations, affecting approximately 1.7/1000 people around the world with ethnic and geographic variation1. Although CL/P is not considered one of the major causes of infant mortality, individuals affected by CL/P it may have difficulties in feeding, speaking, difficult social integration4. Approximately 70% of CL/P cases are considered to be non-syndromic37,38, and their susceptibility has been linked to the expression of various candidate genes through twin studies, familial clustering studies and genome-wide studies39.

Recent studies suggest that using folic acid could reduce the rates of oral clefts6,7. Some genes involved in the metabolism of folic acid such as MTHFR8,9, MTR40, and MTRR41 have been identified. MTHFD1, a crucial gene associated with three sequential enzymatic reactions among 5,10-methylenetetrahydrofolate, 5,10-methenyltetrahydrofolate, 10-formyltetrahydrofolate, tetrahydrofolate, might play a potential role in NSCL/P10. However, controversial results about the MTHFD1 polymorphism rs2236225 (c.1958G>A) have been reported in different articles10,12.

In this systematic review, 10 independent case-control studies were included (eight studies for meta-analysis and two studies qualitatively analyzed) containing 6216 samples (2959 cases and 3257 controls). All the eligible studies of meta-analysis and qualitative analysis showed no significant association of MTHFD1 rs2236225 to the risk of NSCL/P, whether in the whole analysis of five model (A versus G, AG versus GG, AA versus GG, AA + AG versus GG, AA versus AG + GG) or in the subgroup of subjects (mothers or children) and the location of geography (non-Asian countries or Asian countries). Meanwhile, high heterogeneity was observed, which might be the reason for the genetic drift and natural selection among different ethnic groups42. Also, small sample size of different studies might be a possible reason for the disparate results. Our findings suggest that the MTHFD1 polymorphism rs2236225 (c.1958G>A) might not be an appropriate biomarker in predicting the susceptibility of an individual to NSCL/P.

Some limitations of this systematic review and meta-analysis should be noted. Firstly, the choice of retrospective studies has its own limitations, as we may encounter selection bias and influence the results of our analysis43. However, a bigger size of cohort study cannot be conducted easily because of the relatively low morbidity44. Secondly, only 10 studies were included in our review, a small sample size that might not provide sufficient evidence to estimate the connections between the MTHFD1 polymorphisms and the risk of NSCL/P.

NSCL/P is also associated with gene-gene and gene-environment interactions45. Although no correlation was observed between MTFHD1 polymorphism rs2236225 (c.1958G>A) and the risk of NSCL/P, in view of MTFHD1 gene’s key role in folic acid metabolism, we cannot draw a definite conclusion that there is no association between MTFHD1 and NSCL/P’s susceptibility. The use of larger sample size studies, different techniques and considering gene-gene or gene-environment interactions should be explored in future investigations.

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 04 Jun 2015
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Zhao H, Zhang J, Zhang M et al. Is MTHFD1 polymorphisms rs 2236225 (c.1958G>A) associated with the susceptibility of NSCL/P? A systematic review and meta-analysis [version 1; peer review: 1 approved, 1 approved with reservations] F1000Research 2015, 4:142 (https://doi.org/10.12688/f1000research.6425.1)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 1
VERSION 1
PUBLISHED 04 Jun 2015
Views
21
Cite
Reviewer Report 06 Jul 2015
Jose Suazo, Institute for Research in Dental Sciences, Faculty of Dentistry, University of Chile, Santiago, Chile 
Approved with Reservations
VIEWS 21
In general terms, this article contributes to unraveling the complex genetic architecture of NSCL/P. This type of systematic reviews are always relevant due to the fact that they pooled several singles trials, and therefore, the sample sizes have an important ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Suazo J. Reviewer Report For: Is MTHFD1 polymorphisms rs 2236225 (c.1958G>A) associated with the susceptibility of NSCL/P? A systematic review and meta-analysis [version 1; peer review: 1 approved, 1 approved with reservations]. F1000Research 2015, 4:142 (https://doi.org/10.5256/f1000research.6892.r9328)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 06 Jan 2016
    Feng Chen, Laboratory Center, Peking University School and Hospital of Stomatology, Peking, 100081, China
    06 Jan 2016
    Author Response
    Thank you so much for your precious advice. We have studied your comments carefully and revised the manuscript. The point to point responds are listed as following:
    1. We have corrected all
    ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 06 Jan 2016
    Feng Chen, Laboratory Center, Peking University School and Hospital of Stomatology, Peking, 100081, China
    06 Jan 2016
    Author Response
    Thank you so much for your precious advice. We have studied your comments carefully and revised the manuscript. The point to point responds are listed as following:
    1. We have corrected all
    ... Continue reading
Views
25
Cite
Reviewer Report 17 Jun 2015
Jingtan Su, Center for Craniofacial Molecular Biology, University of Southern California, Los Angeles, CA, USA 
Approved
VIEWS 25
Non-syndromic cleft of the lip and/or palate (NSCL/P) is considered closely related to genetic and environmental factors. As a key gene in the metabolism of folic acid which is associated to high risk of NSCL/P, methylenetetrahydrofolate dehydrogenase (MTHFD1) polymorphisms rs ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Su J. Reviewer Report For: Is MTHFD1 polymorphisms rs 2236225 (c.1958G>A) associated with the susceptibility of NSCL/P? A systematic review and meta-analysis [version 1; peer review: 1 approved, 1 approved with reservations]. F1000Research 2015, 4:142 (https://doi.org/10.5256/f1000research.6892.r8906)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 06 Jan 2016
    Feng Chen, Laboratory Center, Peking University School and Hospital of Stomatology, Peking, 100081, China
    06 Jan 2016
    Author Response
    Thank you. We have studied your valuable comments and revised the manuscript according to your suggestions. The point to point responds are listed as following:
    1. We have corrected the error, thank
    ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 06 Jan 2016
    Feng Chen, Laboratory Center, Peking University School and Hospital of Stomatology, Peking, 100081, China
    06 Jan 2016
    Author Response
    Thank you. We have studied your valuable comments and revised the manuscript according to your suggestions. The point to point responds are listed as following:
    1. We have corrected the error, thank
    ... Continue reading

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 04 Jun 2015
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.