Non-syndromic cleft palate: Association analysis on three gene polymorphisms of the folate pathway in Asian and Italian populations

Abstract

Periconceptional folic acid supplementation can reduce the risk of inborn malformations, including orofacial clefts. Polymorphisms of MTHFR, TCN2, and CBS folate-related genes seem to modulate the risk of cleft lip with or without cleft palate (CL/P) in some populations. CL/P and cleft palate only (CPO) are different malformations that share several features and possibly etiological causes. In the present investigation, we conducted a family-based, candidate gene association study of non-syndromic CPO. Three single nucleotide polymorphisms, namely, rs1801133 of MTHFR, rs1801198 of TCN2, and rs4920037 of CBS, were investigated in a sample that included 129 Italian and 65 Asian families. No evidence of association between the three genotyped polymorphisms and CPO was found in the Italian and Asian cases, indeed the transmission disequilibrium test did not detect any asymmetry of transmission of alleles. This investigation, although with some limitation, further supports that CL/P and CPO diverge in their genetic background.

Keywords

Asian and Italian populations cleft palate folic acid

Introduction

Folate metabolism plays a critical role in embryonic development. The relationship between periconceptional folic acid supplementation and reduced risk for neural tube defects is well established.¹ Neural tube and oral cavity are embryologically related. In fact, the neural crest cells that contribute to the orofacial and tooth morphology originate from the dorsolateral regions of the neural tube.

Evidence supporting a preventive effect for folic acid against birth defects have been reported, but contrasting data have been obtained about clefting.² A meta-analysis carried out on 12 case–control studies on folic acid-containing supplement consumption during pregnancy and risk for oral clefts evidenced a combined relative risk of 0.77 and 0.80 for cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO), respectively.³ Non-syndromic orofacial clefts are complex traits with a strong genetic component that interacts with environmental factors. Folate has long been considered one of these factors. Thus, besides periconceptional folic acid dietary intake and supplementation, a number of genetic variants in folate and methionine metabolism-related genes have been considered as potential susceptibility factors.^4,5 Our studies on Italian population provided evidence that maternal genotype at C677T MTHFR gene, as well as specific alleles of one-carbon metabolic pathway genes, can influence the risk for non-syndromic CL/P (nsCL/P).^6–9 On the contrary, the same three polymorphisms (rs1801133 of MTHFR, rs1801198 of TCN2, and rs4920037 of CBS) investigated in independent samples of different ancestries showed discordant results. In particular, 677T allele at MTHFR gene seemed to have a protective effect at the child level, reducing the risk of cleft in an Asian population sample.¹⁰

nsCL/P and non-syndromic CPO (nsCPO) are considered distinct but related malformations.^11,12 It is still a matter of debate whether CL/P and CPO share etiological factors. The aim of this study was to verify whether MTHFR, TCN2, and CBS alleles, which were previously found associated with the increased risk of nsCL/P, could be involved in CPO etiology in the Italian population and in a sample of Asian population from Iran, Tibet, and Bangladesh.

Materials and methods

Sample study

This family-based association study was performed with two different cohorts. The first one included 129 unrelated CPO patients and their parents with Italian ancestry, recruited by the Department of Medical and Surgical Sciences of S. Orsola-Malpighi Polyclinic, Bologna University. The second one that included 65 unrelated Asian probands and their parents was composed of 44 families of Fars ethnicity from the central region of Iran, 11 families from Yushu Tibetan Autonomous Prefecture, and 10 from Bangladesh. We have implemented systematic guidelines to collect samples uniformly from different countries. A clinician ascertained the non-syndromic status of all the probands and that their mothers had neither smoked nor used clefting drugs such as phenytoin, warfarin, and ethanol during pregnancy. In the interview, part of the mothers reported folic acid supplementation in late pregnancy, but all of them excluded that supplementation occurred in periconceptional, or in the early weeks of pregnancy. The study was approved by the ethics committees of the Centers involved in the study, and it complied with the Declaration of Helsinki’s Ethical Principles for Medical Research Involving Human Subjects. Written informed consent was obtained from all patients and parents before study entry. Among the Italian patients enrolled in the study, 120 (93%) were considered sporadic or non-familial cases, while for the remaining 9 probands (7%) at least one more (second, third, or fourth degree) relative was reported as affected by cleft of the palate. Among Asian cases, 56 (86%) were sporadic, while 9 (14%) familial.

Genotyping

Three polymorphisms of one-carbon metabolic pathway genes were selected for this study, because they can influence the risk for nsCL/P (Table 1).^6–9 Polymorphisms were typed using the TaqMan SNP Genotyping Assay (Assay-On-Demand ID: C_1202883_20 for MTHFR rs1801133, C_325467_10 for TCN2 rs1801198, and C_1605440_1 for CBS rs4920037) on a 7500 Sequence Detection System (Life Technologies, Foster City, CA, USA) following the manufacturer’s protocol.

Table 1.

Characteristics of investigated SNPs.

Gene	dbSNP ID^a	Genomic position^b	Location	Alleles^c	Residue change	MAF
MTHFR	rs1801133	chr1:11796321	Exon 5	C/T	A [Ala] → V [Val]	0.38
TCN2	rs1801198	chr22:30615623	Exon 6	C/G	P [Pro] → R [Arg]	0.41
CBS	rs4920037	chr21:43061781	Intron 11	G/A	−	0.19

SNP: single nucleotide polymorphism; MAF: minor allele frequency.

NCBI-SNP database accession numbers.

UCSC Genome Browser on Human December 2013 (GRCh38/hg38) Assembly.

Major allele first.

Statistical analysis

Power calculation indicated that the sample size of the study provided a power >80% to detect linkage disequilibrium with mutations that increase two times the relative risk of CPO. The Hardy–Weinberg equilibrium for genotype distribution, in both proband and parent groups, was examined using Pearson’s χ² test with a significant threshold of P < 0.05.

Association analysis was performed on a Microsoft Windows platform by the PLINK v1.07 software.¹³ Basically, the family-based association analysis was performed by the transmission disequilibrium test (TDT) that examined allele transmission from heterozygous parents to the affected probands.¹⁴ The odds ratio (OR) was used to evaluate the level of association; this corresponds to the ratio M1/M2, where M1 denotes the number of transmissions of the reference allele from heterozygous parents to affected offspring and M2 the number of times that the allele was not transmitted.¹⁵

Family-based association was analyzed with the whole sample to maximize the power of the study. Then, the Italian and Asian samples were analyzed separately in order to avoid that genetic heterogeneity, or different population structure, mask association signals.

Results

Three candidate gene polymorphisms were investigated for association with nsCPO. Genotype frequencies in patients and parents were distributed according to the Hardy–Weinberg equilibrium law, and no Mendelian errors were detected.

The family-based association study was performed by the TDT. The analysis revealed that alleles of candidate genes were transmitted from heterozygous parents to CPO patients according to Mendelian expectation, indicating no association with the malformation (Table 2). This result was obtained with the entire data set and then confirmed with the Italian and Asian populations, singly.

Table 2.

Family-based association analysis of folate pathway polymorphisms in nsCPO.

Sample	Gene	SNP ID	Alleles^a	T^b	U^c	TDT P value	OR (95% CI)
Italian	MTHFR	rs1801133	C/T	57	58	0.93	0.98 (0.68–1.42)
Italian	TCN2	rs1801198	C/G	61	62	0.93	0.98 (0.69–1.40)
Italian	CBS	rs4920037	G/A	39	33	0.48	1.18 (0.74–1.88)
Asian	MTHFR	rs1801133	C/T	14	12	0.69	1.17 (0.54–2.52)
Asian	TCN2	rs1801198	C/G	22	20	0.76	1.10 (0.60–2.02)
Asian	CBS	rs4920037	G/A	14	14	1.00	1.00 (0.48–2.10)
Whole	MTHFR	rs1801133	C/T	71	70	0.93	1.01 (0.73–1.41)
Whole	TCN2	rs1801198	C/G	83	82	0.94	1.01 (0.75–1.37)
Whole	CBS	rs4920037	G/A	53	47	0.55	1.13 (0.76–1.67)

nsCPO: non-syndromic cleft palate only; SNP: single nucleotide polymorphism; TDT: transmission disequilibrium test; OD: odds ratio; CI: confidence interval.

Major allele first.

Number of times the minor allele was transmitted from heterozygous parents.

Number of times the minor allele was untransmitted from heterozygous parents.

Discussion

The orofacial shape is determined by several facial prominences that grow, converge, and differentiate during the second–third months of embryo development. The processes that lead to the orofacial morphology are regulated by an intricate network of inductive signals from ectoderm, mesoderm, endoderm, and neural crest in each prominence.¹⁶ A failure of these crucial steps could originate a malformation. An important observation is the awareness that mothers who increase their intake of folic acid can reduce the risk of having a child with a birth defect, such as neural tube defect or cleft palate.¹⁷ How this could happen is still a matter of study. The link between folate and craniofacial development has also been established in several vertebrate models, in which researchers found a correlation between folate deficiency and/or mutations in folate pathway genes and the presence of orofacial defects.^18–20 Recently, Wahl et al.²¹ tried to elucidate the mechanism of protection acted by folic acid during palate development, in a Xenopus model. The authors demonstrated that folic acid supplementation can reduce orofacial median clefts by enhancing cell survival or boosting proliferation.

The increasing mass of evidence supporting a critical role of folate in orofacial development and the finding of association between one-carbon metabolic pathway genes and nsCL/P collected so far pose the foundation for an nsCPO investigation. This study aimed to verify the involvement of polymorphisms of three folate-related genes—MTHFR, TCN2, and CBS—in nsCPO etiology. Nevertheless, it was not possible to investigate a potential gene–environment interaction by taking into account of folic acid intake during the periconceptional period. Indeed, dietary intake was not evaluated, while folic acid supplementation in the periconceptional period was excluded. Some mothers reported of folic acid intake during pregnancy, but supplementation began in the second trimester, too late to influence the risk of orofacial cleft in offspring.

In a previous attempt, Mills et al.²² did not find a significant association with the variant at C677T MTHFR polymorphism in CPO cases, but instead they found support for a maternal effect. In the same investigation, they excluded a role for C667G TCN2 in their sample study. The CBS rs2124459 minor allele was seen associated with a reduced risk of CPO, in both mothers and children of French and Belgian nsCPO cohorts.²³ Since the risk of nsCPO associated with specific variants in genes crucial for the folic acid metabolism was not extensively studied and it differs among ethnic groups, we investigated our cohorts of Italian and Asian CPO cases for three emblematic polymorphisms linked to folate, vitamin B12, and methionine pathways.

Our results did not support the hypothesis of an involvement of the three investigated genes, at least for the three polymorphisms considered. This could be attributable to the small size of the sample study or to its fragmentation. However, our data confirm the difficulty to find common genetic causes for both CL/P and CPO.

Footnotes

Acknowledgements

We are indebted to the families participating in the study for their invaluable contribution, as well as to all the personnel involved in clinical data and specimen collection.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported in part by a grant from the Interethnos-Interplast Association, Italy and in part by the grant number 189087 from deputy for research, Isfahan University of Medical Sciences, Isfahan, Iran.

ORCID iDs

Francesco Carinci

Luca Scapoli

References

Wolff

Witkop

Miller

, et al. (2009) Folic acid supplementation for the prevention of neural tube defects: An update of the evidence for the U.S. Preventive Services Task Force. Annals of Internal Medicine 150: 632–639.

De-Regil

Pena-Rosas

Fernandez-Gaxiola

, et al. (2015) Effects and safety of periconceptional oral folate supplementation for preventing birth defects. The Cochrane Database of Systematic Reviews 12: CD007950.

Badovinac

Werler

Williams

, et al. (2007) Folic acid-containing supplement consumption during pregnancy and risk for oral clefts: A meta-analysis. Birth Defects Research Part A, Clinical and Molecular Teratology 79(1): 8–15.

Marini

Yang

Asrani

, et al. (2016) Sequence variation in folate pathway genes and risks of human cleft lip with or without cleft palate. American Journal of Medical Genetics Part A 170: 2777–2787.

Jiang

Yin

Zhao

, et al. (2015) Lack of association between MTHFR, MTR, MTRR, and TCN2 genes and nonsyndromic CL±P in a Chinese population: Case-control study and meta-analysis. The Cleft Palate-Craniofacial Journal 52(5): 579–587.

Martinelli

Scapoli

Pezzetti

, et al. (2001) C677T variant form at the MTHFR gene and CL/P: A risk factor for mothers? American Journal of Medical Genetics 98(4): 357–360.

Pezzetti

Martinelli

Scapoli

, et al. (2004) Maternal MTHFR variant forms increase the risk in offspring of isolated nonsyndromic cleft lip with or without cleft palate. Human Mutation 24: 104–105.

Martinelli

Scapoli

Palmieri

, et al. (2006) Study of four genes belonging to the folate pathway: Transcobalamin 2 is involved in the onset of non-syndromic cleft lip with or without cleft palate. Human Mutation 27: 294.

Martinelli

Masiero

Carinci

, et al. (2011) New evidence for the role of cystathionine beta-synthase in non-syndromic cleft lip with or without cleft palate. European Journal of Oral Sciences 119: 193–197.

10.

Martinelli

Girardi

Cura

, et al. (2016) Non-syndromic cleft lip with or without cleft palate in Asian populations: Association analysis on three gene polymorphisms of the folate pathway. Archives of Oral Biology 61: 79–82.

11.

Fogh-Andersen

(1942) Inheritance of Harelip and Cleft Palate. Copenhagen: Busck.

12.

Christensen

Fogh-Andersen

(1994) Etiological subgroups in non-syndromic isolated cleft palate. A genetic-epidemiological study of 52 Danish birth cohorts. Clinical Genetics 46(5): 329–335.

13.

Purcell

Neale

Todd-Brown

, et al. (2007) PLINK: A tool set for whole-genome association and population-based linkage analyses. American Journal of Human Genetics 81(3): 559–575.

14.

Spielman

McGinnis

Ewens

(1993) Transmission test for linkage disequilibrium: The insulin gene region and insulin-dependent diabetes mellitus (IDDM). American Journal of Human Genetics 52(3): 506–516.

15.

Kazeem

Farrall

(2005) Integrating case-control and TDT studies. Annals of Human Genetics 69(Pt. 3): 329–335.

16.

Szabo-Rogers

Smithers

Yakob

, et al. (2010) New directions in craniofacial morphogenesis. Developmental Biology 341(1): 84–94.

17.

Wilcox

Lie

Solvoll

, et al. (2007) Folic acid supplements and risk of facial clefts: National population based case-control study. BMJ 334(7591): 464.

18.

Tang

Santillano

Wlodarczyk

, et al. (2005) Role of Folbp1 in the regional regulation of apoptosis and cell proliferation in the developing neural tube and craniofacies. American Journal of Medical Genetics Part C, Seminars in Medical Genetics 135C(1): 48–58.

19.

Momb

Lewandowski

Bryant

, et al. (2013) Deletion of Mthfd1l causes embryonic lethality and neural tube and craniofacial defects in mice. Proceedings of the National Academy of Sciences of the United States of America 110(2): 549–554.

20.

Maldonado

Lopez

Herrera

, et al. (2018) Craniofacial structure alterations of foetuses from folic acid deficient pregnant mice. Annals of Anatomy 218: 59–68.

21.

Wahl

Kennedy

Wyatt

, et al. (2015) The role of folate metabolism in orofacial development and clefting. Developmental Biology 405(1): 108–122.

22.

Mills

Molloy

Parle-McDermott

, et al. (2008) Folate-related gene polymorphisms as risk factors for cleft lip and cleft palate. Birth Defects Research Part A, Clinical and Molecular Teratology 82(9): 636–643.

23.

Goffinet

Oussalah

Gueant-Rodriguez

, et al. (2016) Cystathionine beta-synthase genetic variant rs2124459 is associated with a reduced risk of cleft palate in French and Belgian populations. Journal of Medical Genetics 53(12): 828–834.