Association of a Polymorphism in the Gene Encoding Methylenetetrahydrofolate Dehydrogenase 1 (MTHFD1) 1958G>A with Orofacial Cleft

ABSTRACT

Objective:

To evaluate the possible association of a polymorphism in the gene encoding methylenetetrahydrofolate dehydrogenase 1 (MTHFD1), 1958G>A, with the susceptibility to orofacial cleft in an Indonesian population.

Material and Methods:

A total of 200 stored secondary biological samples from 30 cases of orofacial cleft and 170 unaffected controls were analyzed to determine the polymorphism status at base 1958. The analysis was conducted using the PCR-restriction fragment length polymorphism technique after digestion with the Msp1 restriction enzyme. The samples were then subjected to agarose gel electrophoresis to investigate the presence or absence of the following fragments: genotype GG, 196, 86 and 40 base pairs (bp); genotype AA, 282 and 28 bp and genotype AG, 282, 196, 86, 40 and 28 bp. The test groups were compared using the Chi-square test.

Results:

The wild-type allele containing 1958G, as well as the genotype GG, were significantly more common in the control group than in the orofacial cleft group.

Conclusion:

The MTHFD1 1958G>A polymorphism was significantly associated with orofacial cleft susceptibility in the tested Indonesian population.

Keywords:
Cleft Palate; Cleft Lip; Genetic Variation; Polymorphism, Genetic

Introduction

Orofacial cleft is a type of congenital malformation of the oral cavity and palate [¹]. Previous epidemiologic studies have suggested that the overall prevalence of orofacial cleft is approximately 2 cases per 1000 births. In Indonesia, the corresponding prevalence is approximately 1 case per 1000 live births in 24-59-month-old children [²]. Non-syndromic orofacial cleft is a multifactorial condition attributed to both genetic and environmental factors, including maternal smoking, alcohol consumption and folate deficiency during the peri-conceptual period [³,⁴]. Particularly, folic acid deficiency has a known impact not only on the manifestation of orofacial cleft, but also on other congenital malformations such as neural tube defects.

In addition to environmental factors, several genes have been identified as potential contributors to the developmental embryogenesis of orofacial cleft. These genes include MTHFD1, a gene that encodes methylenetetrahydrofolate dehydrogenase 1 (MTHFD1) and is thus involved in folate metabolism [⁵,⁶]. MTHFD1, which has a chromosomal location of 14q23.3, plays an important role in the generation of one-carbon metabolic derivatives of tetrahydrofolate, including 5,10-methylenetetrahydrofolate cyclohydrolase, 10-formyl tetrahydrofolate and 5,10-methylenetetrahydrofolate dehydrogenase [⁵].

One previous study revealed the associations of MTHFD1 variants with homocysteine and folic acid levels. Moreover, one common MTHFD1 polymorphism, 1958G>A, can lead to the amino acid substitution of arginine for glutamine, which may reduce the activity and stability of the resulting enzyme and increase the risk of deviations such as orofacial cleft [⁷]. This polymorphism may also be associated with other craniofacial structural abnormalities, given the expression of MTHFD1 within localized regions in the brain, neural tube, limb buds and tail bud at all stages of embryogenesis [⁸]. We thereby hypothesized that the MTHFD1 1958G>A polymorphism may be associated with susceptibility to orofacial cleft, and aimed to explore this hypothesis in an Indonesian population.

Material and Methods

Study Design

The study included secondary biological samples of DNA from peripheral blood of 30 cases of orofacial cleft and 170 unaffected controls using the techniques by previous authors [⁹-¹¹], which were made available from the Oral Biology Laboratory of the Faculty of Dentistry, Universitas Indonesia.

Laboratory Procedures

For DNA isolation, blood samples were inserted into Falcon tubes containing red blood cell anticoagulant (1:3), and the tube was repeatedly inverted for 10 min. The tubes were then centrifuged for 15 min at 1500 rpm at room temperature, supernatant was discarded, and this was repeated 3 times. After adding 2 mL of cell lysis solution (CLS) to the pellet and mixing by inverting, the solution was incubated in a water bath at 37°C for 60 min. Protein precipitation (PP) solution of 1.3 mL was added, vortexed and centrifuged at 40°C/3000 rpm for 5 min. The supernatant was poured into a new tube containing 2.3 mL of cold isopropanol, and the tube was inverted repeatedly until DNA in it appeared white. Then, 1.3 mL of cold 70% ethanol was added to precipitate DNA that was dried by inverting the tube. Washing was repeated two times, and the tube was dried in tilted position. Then, 0.3 mL of Tris-EDTA buffer (TE) was added (adjusted by the amount of DNA), incubated in a water bath at 37°C for 2 h, and finally the solution was transferred into Eppendorf tubes and stored at −20°C.

The following primers were used to amplify a 310-bp PCR product (Figure 1a) containing the target SNP: forward, 5´-CCTGGTTTCCACAGGGCACTC-3´ and reverse, 5´-CCACGTGGGGGCAGAGGCCGGAATACCGG-3´. The PCR amplification protocol comprised an initial denaturation step at 94°C for 7 minutes, followed by 36 cycles of denaturation at 94°C for 60 seconds, annealing at 60°C for 60 seconds and extension at 72°C for 60 seconds.

Figure 1
Visualization of a) the 310-bp PCR product in all lanes except M, which contains a 50-bp reference ladder, and b) MTHFD1 1958G>A polymorphism genotypes AA (lane 1), AG (lane 2) and GG (lane 3).

The PCR products were subject to genotyping via the restriction fragment length polymorphism (RFLP) method after digestion with the restriction enzyme Msp1. After a 16-hour incubation in a water bath at 37°C, the RFLP products were incubated in a thermoblock at 65°C for 20 minutes for enzyme inactivation. The samples were then subjected to agarose gel electrophoresis to investigate the presence or absence of the following fragments: genotype GG, 196, 86 and 40 base pairs (bp); genotype AA, 282 and 28 bp and genotype AG, 282, 196, 86, 40 and 28 bp (Figure 1b).

Data Analysis

The test groups were compared using the Chi-square test. SPSS software (IBM Corp., Armonk, NY, USA) was used for the statistical analysis. A p-value of <0.05 was considered to indicate a statistically significant difference.

Ethical Clearance

The protocol was subject to the Ethical approval no.17/Ethical Exempted/FKGUI/XII/2015, protocol No. 01560615. All samples originated from native Indonesian patients who resided mainly on the island of Java.

Results

Table 1 presents the distributions of genotypes and alleles of the MTHFD1 1958G>A polymorphism in the orofacial cleft and control groups. The GG genotype (wild type) and G allele were clearly dominant in the control group, but appeared at significantly reduced frequencies in the case group. The genotype distribution of Table 1 for the case group is consistent with the Hardy-Weinberg equilibrium (HWE) (p=0.543), but the same is not true for the control group (p=<0.001).

Table 1
Distributions of genotypes and alleles associated with the MTHFD1 1958G>A polymorphism in test groups with and without orofacial cleft.

	Orofacial Cleft
Variables	Yes	No	p-value
	N (%)	N (%)
Genotype
GG	2 (6.7)	78 (45.9)	0.001
AA	14 (46.7)	43 (25.3)
GA	14 (46.7)	49 (28.8)
Total	30 (100.0)	170 (100.0)
Allele
G	18 (30.0)	205 (60.3)	0.001
A	42 (70.0)	135 (39.7)
Total	60 (100.0)	340 (100.0)

Discussion

Our study of samples from an Indonesian population suggested a significant association of the MTHFD1 1958G>A (rs2236225) polymorphism with orofacial cleft (p=0.001), consistent with the results from previous studies conducted in India and China [⁶,¹²]. In contrast to studies of Asian populations, most previous studies involving European (Norwegian, Italian, and Polish) and South American (Brazilian) populations indicated no significant association of this polymorphism with the risk of orofacial cleft [¹³-¹⁷], although one study also suggested such an association in a European (Irish) population [¹⁸]. More recently, a review including meta-analyses has casting doubts on the reported association when including studies on Asian populations [¹⁹]. Another study has confirmed the lacking association in a wider range of European populations from 11 countries [²⁰]. In other words, although the presumed mechanism underlying the effect of this polymorphism on folate metabolism appears reasonable, the associated evidence does not consistently support this association. The possible effects of ethnic origin, as well as the type of oral cleft (lip, palate or both), may warrant further study [⁴]. The conclusions have been generally consistent for a number of the reviewed studies that have included all these types of nonsyndromic cleft cases, and a wide range of populations with different ethnic origins [¹⁹,²⁰]. Nevertheless, the overall coverage may be instructive but by no means complete: for example, the coverage on African and South-East Asian populations remains rather poor. Unfortunately, the number of samples from patients with orofacial cleft available for the present study was too small for a meaningful subgroup analysis.

While the genotype distribution for the tested polymorphism in the case group from 30 samples of orofacial cleft patients was consistent with HWE, the sample was relatively small. With 170 samples, the control group was considerably larger, but not consistent with HWE-. In comparison, in the other Asian (Indian and Chinese) studies that indicated a significant association of the same polymorphism and orofacial cleft, the genotypes of the control group were consistent with HWE [⁶,¹²,¹⁹]. Therefore, it appears reasonable to suggest an extension of the present study to a larger case group, for example classified according to the type of cleft, and to extend the comparison to another control group, both from the same ethnic origin.

Other polymorphisms of MTHFD1 than 1958G>A have been previously considered for association with oral cleft susceptibility. In particular, the polymorphisms 31136A>G (rs2357694) and 18913T>C (rs8006686) of MTHFD1 have shown no significant association so far [¹²], and the same was the case for the polymorphisms rs1950902, rs10813 and rs17857382 in another study [¹⁶]. In parallel, other genes encoding enzymes of the folate metabolism, such as MTHFR, MTR and MTRR, and their polymorphic variants, have been addressed for significance in previous studies, often in studies together with MTHFD1 [¹³,¹⁴,¹⁶,¹⁷]. Variants with indicated association to oral cleft include e.g. MTHFR rs2274976 [¹⁶,¹⁷] and with more controversial results, rs1801133 (C677T) of the same gene [¹⁷,²¹,²²].

Although individual polymorphisms may exhibit little or weak association with oral cleft, it may still exert influence in combination with other specific polymorphisms. Such influence has been indicated also for MTHFD1 [¹⁶,¹⁹,²⁰].

For studies on the genetic component in acquiring oral cleft, further applied options include sampling of twins or triads for wider comparisons, and looking into details of e.g. lifestyle and environmental features during pregnancy, to explore confounding factors [¹²,¹⁶]. Such details were not available for the source material of the present study, and it is suggested to expand it first by simply increasing the number of comparable samples.

Conclusion

Our study suggests a significant association of the MTHFD1 1958G>A polymorphism with susceptibility to orofacial cleft in an Indonesian population. Nonetheless, the contradictory results produced by other studies may warrant further analyses of the effects of this polymorphism on other aspects, such as the type of cleft.

• Financial Support
  University of Indonesia (Grant No. NKB-4188/UN2.RST/HKP.05.00/2020).
• Data Availability
  The data used to support the findings of this study can be made available upon request to the corresponding author.

References

• ^[1]
  Allam E, Windsor LJ, Stone C, Allam E, Windsor L, Stone C. Cleft lip and palate: etiology, epidemiology, preventive and intervention strategies. Anat Physiol 2014; 4:150. https://doi.org/10.4172/2161-0940.1000150
  » https://doi.org/10.4172/2161-0940.1000150
• ^[2]
  KEMENKESRI. Riset kesehatan dasar. Jakarta: Badan Penelitian dan Pengembangan Kesehatan Departemen Kesehatan Republik Indonesia; 2013. p. 228. [In Indonesian].
• ^[3]
  Allagh KP, Shamanna BR, Murthy GV, Ness AR, Doyle P, Neogi SB, et al. Birth prevalence of neural tube defects and orofacial clefts in India: A systematic review and meta-analysis. Plos One 2015; 10(3):e0118961. https://doi.org/10.1371/journal.pone.0118961
  » https://doi.org/10.1371/journal.pone.0118961
• ^[4]
  Lesmana S, Auerkari EI. Genes contributing in cleft lip and cleft palate: A literature review. J Int Dent Med Res 2017; 9:441-8.
• ^[5]
  Neagos D, Cretu R, Tutulan-Cunita A, Stoian V, Bohiltea LC. Methylenetetrahydrofolate dehydrogenase (MTHFD) enzyme polymorphism as a maternal risk factor for trisomy 21: a clinical study. J Med Life 2010; 3(4):454-7.
• ^[6]
  Murthy J, Gurramkonda VB, Lakkakula BV. Significant association of MTHFD1 1958G<A single nucleotide polymorphism with nonsyndromic cleft lip and palate in Indian population. Med Oral Patol Oral Cir Bucal 2014; 19(6):e616-21. https://doi.org/10.4317/medoral.19796
  » https://doi.org/10.4317/medoral.19796
• ^[7]
  Wehby GL, Murray JC. Folic acid and orofacial clefts: a review of the evidence. Oral Dis 2010; 16(1):11-9. https://doi.org/10.1111/j.1601-0825.2009.01587.x
  » https://doi.org/10.1111/j.1601-0825.2009.01587.x
• ^[8]
  Momb J, Lewandowski JP, Bryant JD, Fitch R, Surman DR, Vokes SA, et al. Deletion of Mthfd1l causes embryonic lethality and neural tube and craniofacial defects in mice. Proc Natl Acad Sci U S A 2013; 110(2):549-54. https://doi.org/10.1073/pnas.1211199110
  » https://doi.org/10.1073/pnas.1211199110
• ^[9]
  Tania SF, Kim KR, Suhartono AW, Latief EI Auerkari. Distribution of E-cadherin 1 (CDH1) promoter methylation in patients with orofacial cleft. J Phys Conf Ser 2018; 1073(3):032064. https://doi.org/10.1088/1742-6596/1073/3/032064
  » https://doi.org/10.1088/1742-6596/1073/3/032064
• ^[10]
  Jessica C, Alwadris TT, Prasetyo SR, Puspitawati R, Auerkari EI. Association of interleukin 8-251 A/T gene polymorphism with periodontitis in Indonesia. J Phys: Conf Ser 2018;1025(1):012071. https://doi.org/10.1088/1742-6596/1025/1/012071
  » https://doi.org/10.1088/1742-6596/1025/1/012071
• ^[11]
  Auerkari EI, Suryandari DA, Umami SS, Kusdhany LS, Siregar TWA, Rahardjo TBW, Talbot C, Hogervorst E. Gene promoter polymorphism of RUNX2 and risk of osteoporosis in postmenopausal Indonesian women. SAGE Open Med 2014; 2:2050312114531571. https://doi.org/10.1177/2050312114531571
  » https://doi.org/10.1177/2050312114531571
• ^[12]
  Wu J, Chen Y, Pei J, Pan J. MTHFD1 gene polymorphisms as risk factors involved in orofacial cleft: an independent case-control study and a meta-analysis. Int J Clin Exp Med 2015; 8(5):7737-45.
• ^[13]
  Mostowska A, Hozyasz KK, Jagodzinski PP. Maternal MTR genotype contributes to the risk of non-syndromic cleft lip and palate in the Polish population. Clin Genet 2006; 69(6): 512-7.
• ^[14]
  Boyles AL, Wilcox AJ, Taylor JA, Meyer K, Fredriksen Å, Ueland PM, et al. Folate and one-carbon metabolism gene polymorphisms and their associations with oral facial clefts. Am J Med Genet A 2008; 146A(4): 440-9.
• ^[15]
  Palmieri A, Masiero E, Martinelli M, Scapoli L, Pezzetti F, Caramelli E, et al. The MTHFD1 gene is not involved in cleft lip with or without palate onset among the Italian population. Hum Genet 2008; 72(3):297-9. https://doi.org/10.1111/j.1469-1809.2007.00428.x
  » https://doi.org/10.1111/j.1469-1809.2007.00428.x
• ^[16]
  Bufalino A, Paranaiba LMR, De Aquino SN, Martelli-Junior H, Swerts MSO, Coletta RD. Maternal polymorphisms in folic acid metabolic genes are associated with nonsyndromic cleft lip and/or palate in the Brazilian population. Birth Defects Res A 2010; 88:980-6.
• ^[17]
  De Aquino SN, Hoshi R, Bakordakis E, Pucciarelli MGR, Messetti AC, Moreira H, et al. MTHFR rs2274976 polymorphism is a risk marker for nonsyndromic cleft lip with or without cleft palate in the Brazilian population. Birth Defects Res A 2014; 100:30-5.
• ^[18]
  Mills JL, Molloy AM, Parle-McDermott A, Troendle JF, Brody LC, Conley MR, et al. Folate-related gene polymorphisms as risk factors for cleft lip and cleft palate. Birth Defects Res A 2008; 82:636-43.
• ^[19]
  Zhao H, Zhang J, Zhang M, Deng F, Zheng L, Zheng H, et al. Is MTHFD1 polymorphism rs2236225 (c. 1958G>A) associated with the susceptibility of NSCL/P? A systematic review and meta-analysis. F1000Research 2015, 4:142. https://doi.org/10.12688/f1000research.6425.2
  » https://doi.org/10.12688/f1000research.6425.2
• ^[20]
  Mossey PA, Little J, Steegers-Theunissen R, Molloy A, Peterlin B, Shaw WC, et al. Genetic interactions in nonsyndromic orofacial clefts in Europe - EUROCRAN study. Cleft Palate Craniofac J 2017; 54(6):623-30. https://doi.org/10.1597/16-037
  » https://doi.org/10.1597/16-037
• ^[21]
  Luo YL, Cheng YL, Ye P, Wang W, Xiao H, Chen Q. Association between MTHFR polymorphisms and orofacial clefts risk: a meta-analysis. Birth Defects Res A 2012; 94:237-44.
• ^[22]
  Butali A, Little J, Chevrier C, Cordier S, Steegers-Theunissen R, Jugessur S, et al. Folic acid supplementation use and the MTHFR C667T polymorphism in orofacial clefts etiology: an individual participant data pooled-analysis. Birth Defects Res A 2013; 97:509-14.

Edited by

Academic Editor: Alessandro Leite Cavalcanti

Publication Dates

• Publication in this collection
  02 Apr 2021
• Date of issue
  2021

History

• Received
  19 Jan 2020
• Reviewed
  28 Aug 2020
• Accepted
  14 Oct 2020