Few associations of candidate genes with nonsyndromic orofacial clefts in the population of Lithuania

Nonsyndromic orofacial clefting (NS-OFC) is a common complex multifactorial trait with a considerable genetic component and a number of candidate genes suggested by various approaches. Twenty biallelic and microsatellite DNA markers in the strong candidate loci TGFA, TGFB3, GABRB3, RARA, and BCL3 were analysed for allelic association with the NS-OFC phenotype in 112 nuclear families (proband + both parents) from Lithuania by using the transmission disequilibrium test (TDT). Associations were found between the TGFA gene marker rs2166975 and nonsyndromic cleft palate only (CPO) phenotype (p = 0.045, df 1) as well as between the D2S292 marker and the cleft lip with or without cleft palate (CL/CP) phenotype in allele-wise TDT (P = 0.005, df 9) and genotype-wise TDT (P = 0.021, df 24). A weak association (P = 0.085, df 3) of the BCL3 marker (BCL3 gene) with the risk of CPO was also found. Thus our initial results support the contribution of allelic variation in the TGFA locus to the aetiology of CL/CP in the population of Lithuania but they do not point to TGFA as a major causal gene. Different roles of the TGFA and BCL3 genes in the susceptibility to NS-OFC phenotypes are suggested.     

Linkage disequilibrium between GABRB3 gene and nonsyndromic familial cleft lip with or without cleft palate

The malformation of nonsyndromic cleft lip with or without cleft palate (CL/P) is a common congenital disease that affects approximately 1/1000 newborns in Caucasian populations. Genetic studies indicate that CL/P has the characteristics of a complex genetic trait. Linkage analysis and mouse-model knockout studies have suggested several candidate genes mapping in different chromosome regions for CL/P malformation. On these grounds, we have investigated, by linkage disequilibrium (LD) and parametric and nonparametric linkage analyses, five different candidate genes, including those for the beta3 subunit of the gamma-aminobutyric acid receptor (GABRB3), glutamic acid decarboxylase 1 (GAD1), retinoic acid receptor alpha (RARA), transforming growth factor beta3 (TGFB3), and msh ( Drosophila) homeobox homolog 1 (MSX1). Interestingly, a significant LD between GABRB3 and CL/P was obtained ( P-value=0.008 in the allele-wise analysis for multiallelic markers), suggesting that the GABRB3 gene is involved in this congenital disease. This new finding in humans is in agreement with previously reported data obtained with the murine model. Indeed, mouse studies indicate a role for gamma-aminobutyric acid (GABA) and its receptor in normal palate development. Exclusion of the GAD1 gene, which encodes the GABA-producing enzyme, in CL/P pathogenesis was obtained in our study. Moreover, we were unable to confirm the involvement of the MSX1 gene in nonsyndromic CL/P. Modest evidence of LD between marker alleles and CL/P was found at the RARA and TGFB3 loci suggesting a minor role for these genes in our family set of nonsyndromic CL/P.     

Association of genetic variation of the transforming growth factor-alpha gene with cleft lip and palate.

Complex segregation analysis of pedigrees having nonsyndromic cleft lip with or without cleft palate (CL/P) (Chung et al. 1986; Marazita et al. 1986) has shown that a major-locus model best explains the observed recurrence of CL/P in Caucasian families. To identify this major gene, we compared the frequencies of 12 RFLPs at five loci-epidermal growth factor, transforming growth factor-alpha, epidermal growth factor receptor, glucocorticoid receptor, and estrogen receptor-in both a group of 80 subjects with nonsyndromic CL/P and 102 controls. These candidate genes were selected because studies in rodents had suggested their possible involvement in palatogenesis. A significant association was observed between two RFLPs at the transforming-growth-factor-alpha (TGFA) locus and the occurrence of clefting (P = .0047 and P = .0052). This suggests that either the TGFA gene itself or DNA sequences in an adjacent region contribute to the development of a portion of cases of CL/P in humans and provides an opportunity to begin to examine the molecular events underlying lip and palate formation.     

Cleft lip with or without cleft palate: associations with transforming growth factor alpha and retinoic acid receptor loci.

The first association study of cleft lip with or without cleft palate (CL/P), with candidate genes, found an association with the transforming growth-factor alpha (TGFA) locus. This finding has since been replicated, in whole or in part, in three independent studies. Here we extend our original analysis of the TGFA TaqI RFLP to two other TGFA RFLPs and seven other RFLPs at five candidate genes in 117 nonsyndromic cases of CL/P and 113 controls. The other candidate genes were the retinoic acid receptor (RARA), the bcl-2 oncogene, and the homeobox genes 2F, 2G, and EN2. Significant associations with the TGFA TaqI and BamHI RFLPs were confirmed, although associations of clefting with previously reported haplotypes did not reach significance. Of particular interest, in view of the known teratogenic role of retinoic acid, was a significant association with the RARA PstI RFLP (P = .016; not corrected for multiple testing). The effect on risk of the A2 allele appears to be additive, and although the A2A2 homozygote only has an odds ratio of about 2 and recurrence risk to first-degree relatives (lambda 1) of 1.06, because it is so common it may account for as much as a third of the attributable risk of clefting. There is no evidence of interaction between the TGFA and RARA polymorphisms on risk, and jointly they appear to account for almost half the attributable risk of clefting.     

Medical sequencing of candidate genes for nonsyndromic cleft lip and palate

Nonsyndromic or isolated cleft lip with or without cleft palate (CL/P) occurs in wide geographic distribution with an average birth prevalence of 1/700. We used direct sequencing as an approach to study candidate genes for CL/P. We report here the results of sequencing on 20 candidate genes for clefts in 184 cases with CL/P selected with an emphasis on severity and positive family history. Genes were selected based on expression patterns, animal models, and/or role in known human clefting syndromes. For seven genes with identified coding mutations that are potentially etiologic, we performed linkage disequilibrium studies as well in 501 family triads (affected child/mother/father). The recently reported MSX1 P147Q mutation was also studied in an additional 1,098 cleft cases. Selected missense mutations were screened in 1,064 controls from unrelated individuals on the Centre d'Étude du Polymorphisme Humain (CEPH) diversity cell line panel. Our aggregate data suggest that point mutations in these candidate genes are likely to contribute to 6% of isolated clefts, particularly those with more severe phenotypes (bilateral cleft of the lip with cleft palate). Additional cases, possibly due to microdeletions or isodisomy, were also detected and may contribute to clefts as well. Sequence analysis alone suggests that point mutations in FOXE1, GLI2, JAG2, LHX8, MSX1, MSX2, SATB2, SKI, SPRY2, and TBX10 may be rare causes of isolated cleft lip with or without cleft palate, and the linkage disequilibrium data support a larger, as yet unspecified, role for variants in or near MSX2, JAG2, and SKI. This study also illustrates the need to test large numbers of controls to distinguish rare polymorphic variants and prioritize functional studies for rare point mutations.     

Transforming growth factor-alpha (TGFA): genomic structure, boundary sequences, and mutation analysis in nonsyndromic cleft lip/palate and cleft palate only.

Transforming growth factor-alpha (TGFA) has been proposed as a candidate gene in the etiology of nonsyndromic cleft lip with or without cleft palate (NS-CL/P) and of nonsyndromic cleft palate only (NS-CPO). Biologic support for a role of TGFA arises from its presence at high levels in the epithelial tissue of the medial edge of the palatal shelves at the time of shelf fusion in mice. Genetic support for the role of TGFA in clefting comes from the reported association of TGFA alleles with human NS-CPO and NS-CL/P. In this study we report the sequence and structure of human genomic TGFA and the search for causal TGFA mutations in 250 individuals with NS-CL/P or NS-CPO by conformational analysis of the coding sequence, splice junctions, and a portion of the 3' untranslated region strongly homologous between human and mouse. We confirm that human TGFA is composed of six exons and here report several new sequence substitutions and their frequencies. Five variants in conserved segments may represent rare causes for clefting in humans and provide support for the role of TGFA in facial morphogenesis.     

Aryl hydrocarbon receptor nuclear translocator 2 (ARNT2): structure,gene mapping,polymorphisms, and candidate evaluation for human orofacial clefts

BACKGROUND: Nonsyndromic orofacial clefts have an estimated incidence of 1/1000 live births. Population genetic and embryologic studies suggest that cleft palate only (CPO) may be a distinct clinical entity from cleft lip with or without cleft palate (CL/P). Both CPO and CL/P are thought to be multifactorial in etiology, with evidence indicating that genetic, environmental, and developmental determinants may all play a role. The ARNT2 gene localizes to a conserved linkage group on mouse chromosome 7 that is syntenic with human chromosome 15q23-25. This chromosomal region was previously identified as a teratogen-induced clefting susceptibility locus in a genome-wide scan of AXB and BXA recombinant inbred mice. Arnt2 is expressed in the first branchial arch in mice. The teratogen 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) acts through the aryl hydrocarbon receptor (Ahr) pathway to produce dose-dependent CPO and thymic wasting in mice exposed in utero. Arnt2 and Ahr proteins dimerize in vitro. TCDD exposure is also associated with orofacial clefting in children of parents involved in agricultural work. METHODS: To determine whether ARNT2 influences human craniofacial development, we identified the human ARNT2 gene and conducted genomic structural analysis. Mutational screening was performed in infants with nonsyndromic CPO or CL/P who were identified by the Iowa Birth Defects Registry. RESULTS: A common amino acid polymorphism was detected but, no obvious disease-causing mutations were detected by SSCP analysis. The microsatellite marker, GATA89D04 (D15S823) was identified within intron 11 of the human ARNT2 gene, and linkage disequilibrium of nonsyndromic CPO and CL/P parent-infant trios was conducted. CONCLUSIONS: No association was demonstrated with CPO (n = 45) and CL/P (n = 37). Teratology 66:85-90, 2002.     

Genetic variation in transforming growth factor alpha: possible association of BamHI polymorphism with bilateral sporadic cleft lip and palate

Non-syndromic cleft lip with or without cleft palate (CL/P) is one of the most common birth defects affecting 1/1000 Caucasians. Genetic factors are thought to contribute to the development of this disorder. A significant association between two restriction fragment length polymorphisms, the TGF TaqI 2.7-kb allele and the TGF BamHI 40-kb allele, at the transforming growth factor alpha (TGF) locus and the occurrence of clefting has previously been reported. A total of 98 Caucasian patients of Alsacian ancestry was recruited from our registry of congenital malformations. These patients had isolated CL/P but no other anomalies. In addition 57 patients with cleft palate, but without cleft lip, were studied. A control group comprised 99 unrelated healthy Caucasians of the same Alsacian ancestry. TaqI and BamHI identify two-allele polymorphisms. The TGFA Taq and BamHI alleles showed no significant association with the presence of clefting, the only exception being that the BamHI 10.0-kb allele was significantly more frequent in patients with bilateral CL/P.     

Nonsyndromic cleft lip with or without cleft palate: evidence of linkage to BCL3 in 17 multigenerational families.

Nonsyndromic cleft lip with or without cleft palate (CL/P) is a common craniofacial developmental defect. Recent segregation analyses have suggested that major genes play a role in the etiology of CL/P. Linkage to 22 candidate genes was tested in 11 multigenerational families with CL/P, and 21 of these candidates were excluded. APOC2, 19q13.1, which is linked to the proto-oncogene BCL3, gave suggestive evidence for linkage to CL/P. The study was expanded to include a total of 39 multigenerational CL/P families. Linkage was tested in all families, using an anonymous marker, D19S178, and intragenic markers in BCL3 and APOC2. Linkage was tested under two models, autosomal dominant with reduced penetrance and affecteds only. Homogeneity testing on the two-point data gave evidence of heterogeneity at APOC2 under the affecteds-only model. Both models showed evidence of heterogeneity, with 43% of families linked at zero recombination to BCL3 when marker data from BCL3 and APOC2 were included. A maximum multipoint LOD score of 7.00 at BCL3 was found among the 17 families that had posterior probabilities > = 50% in favor of linkage. The transmission disequilibrium test provided additional evidence for linkage with the 3 allele of BCL3 more often transmitted to affected children. These results suggest that BCL3, or a nearby gene, plays a role in the etiology of CL/P in some families.     

Analysis of candidate genes on chromosome 2 in oral cleft case-parent trios from three populations

Isolated oral clefts, including cleft lip with/without cleft palate (CL/P) and cleft palate (CP), have a complex and heterogeneous etiology. Case-parent trios from three populations were used to study genes spanning chromosome 2, where single nucleotide polymorphic (SNP) markers were analyzed individually and as haplotypes. Case-parent trios from three populations (74 from Maryland, 64 from Singapore and 95 from Taiwan) were genotyped for 962 SNPs in 104 genes on chromosome 2, including two well-recognized candidate genes: TGFA and SATB2. Individual SNPs and haplotypes (in sliding windows of 2–5 SNPs) were used to test for linkage and disequilibrium separately in CL/P and CP trios. A novel candidate gene (ZNF533) showed consistent evidence of linkage and disequilibrium in all three populations for both CL/P and CP. SNPs in key regions of ZNF533 showed considerable variability in estimated genotypic odds ratios and their significance, suggesting allelic heterogeneity. Haplotype frequencies for regions of ZNF533 were estimated and used to partition genetic variance into among-and within-population components. Wright’s fixation index, a measure of genetic diversity, showed little difference between Singapore and Taiwan compared with Maryland. The tensin-1 gene (TNS1) also showed evidence of linkage and disequilibrium among both CL/P and CP trios in all three populations, albeit at a lower level of significance. Additional genes (VAX2, GLI2, ZHFX1B on 2p; WNT6–WNT10A and COL4A3–COL4A4 on 2q) showed consistent evidence of linkage and disequilibrium only among CL/P trios in all three populations, and TGFA showed significant evidence in two of three populations.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Velopharyngeal variations in relatives of cleft-affected individuals

Eighteen families each with two or more cleft lip and palate patients were studied by speech cephalometry for evidence of velopharyngeal inadequacy (VPI). With this total of 56 persons, three groups were recognized: 1) patients with cleft lip only (N = 7), 2) unaffected sibs of CL(P) probands and the unaffected parents with the positive clefting history on their side of the family (N = 33), and 3) unaffected parents with negative CL(P) history to their side of the family (N = 16). The latter served as controls. The velopharyngeal mechanism in function was evaluated by voicing the fricative/S/. The results showed no significant differences in the length of either the resting soft palate or pharyngeal depth among the three groups. Even though a significant (P less than 0.01) increase in soft palate length while voicing /S/ was found in group 2 relatives compared to group 3 controls, the failure to find differences in either resting palate length or pharyngeal depth coupled with a failure to demonstrate VPI in group 2 subjects by speech testing leaves the value of this observation uncertain.     

Distinct functions for Bmp signaling in lip and palate fusion in mice

Previous work suggested that cleft lip with or without cleft palate (CL/P) is genetically distinct from isolated cleft secondary palate (CP). Mutations in the Bmp target gene Msx1 in families with both forms of orofacial clefting has implicated Bmp signaling in both pathways. To dissect the function of Bmp signaling in orofacial clefting, we conditionally inactivated the type 1 Bmp receptor Bmpr1a in the facial primordia, using the Nestin cre transgenic line. Nestin cre; Bmpr1a mutants had completely penetrant, bilateral CL/P with arrested tooth formation. The cleft secondary palate of Nestin cre; Bmpr1a mutant embryos was associated with diminished cell proliferation in maxillary process mesenchyme and defective anterior posterior patterning. By contrast, we observed elevated apoptosis in the fusing region of the Nestin cre; Bmpr1a mutant medial nasal process. Moreover, conditional inactivation of the Bmp4 gene using the Nestin cre transgenic line resulted in isolated cleft lip. Our data uncover a Bmp4-Bmpr1a genetic pathway that functions in lip fusion, and reveal that Bmp signaling has distinct roles in lip and palate fusion.     

Definition of Critical Periods for Hedgehog Pathway Antagonist-Induced Holoprosencephaly,Cleft Lip,and Cleft Palate

The Hedgehog (Hh) signaling pathway mediates multiple spatiotemporally-specific aspects of brain and face development. Genetic and chemical disruptions of the pathway are known to result in an array of structural malformations, including holoprosencephaly (HPE), clefts of the lip with or without cleft palate (CL/P), and clefts of the secondary palate only (CPO). Here, we examined patterns of dysmorphology caused by acute, stage-specific Hh signaling inhibition. Timed-pregnant wildtype C57BL/6J mice were administered a single dose of the potent pathway antagonist vismodegib at discrete time points between gestational day (GD) 7.0 and 10.0, an interval approximately corresponding to the 15^th to 24^th days of human gestation. The resultant pattern of facial and brain dysmorphology was dependent upon stage of exposure. Insult between GD7.0 and GD8.25 resulted in HPE, with peak incidence following exposure at GD7.5. Unilateral clefts of the lip extending into the primary palate were also observed, with peak incidence following exposure at GD8.875. Insult between GD9.0 and GD10.0 resulted in CPO and forelimb abnormalities. We have previously demonstrated that Hh antagonist-induced cleft lip results from deficiency of the medial nasal process and show here that CPO is associated with reduced growth of the maxillary-derived palatal shelves. By defining the critical periods for the induction of HPE, CL/P, and CPO with fine temporal resolution, these results provide a mechanism by which Hh pathway disruption can result in “non-syndromic” orofacial clefting, or HPE with or without co-occurring clefts. This study also establishes a novel and tractable mouse model of human craniofacial malformations using a single dose of a commercially available and pathway-specific drug.     

Identification of susceptibility loci for nonsyndromic cleft lip with or without cleft palate in a two stage genome scan of affected sib-pairs

Nonsyndromic cleft lip with or without cleft palate (CL/P) is a complex disorder of multigenic origin involving between two and ten loci. Linkage and association studies of CL/P have implicated a number of candidate genes and regions but have often proved difficult to replicate. Here, we report the findings from a two-stage genome-wide scan of 92 affected sib-pairs to identify susceptibility loci to CL/P. An initial set of 400 microsatellite markers was used, with an average spacing of 10 cM throughout the genome. Eleven regions on eight chromosomes were found to have a P-value smaller than 0.05. These eight chromosomes were then further mapped with a second set of markers to increase the average map density to 5 cM. In seven out of eleven areas densely mapped, significance was markedly increased by decreasing the marker interval. Excessive allele sharing was found at 1p (NPL=2.35, P=0.009, MLS=1.51), 2p (NPL=1.77, P= 0.04, MLS=0.66), 6p (NPL=2.35, P=0.009, MLS=1.34), 8q (NPL=2.15, P=0.015, MLS= 1.51) 11 cen (NPL=2.70, P=0.003, MLS=2.10), 12q (NPL=2.08, P=0.02, MLS= 1.5), 16p (NPL=2.1, P=0.018, MLS=0.97) and Xcen-q (NPL=2.40, P=0.008, MLS=2.68). Although none reached the level required for significant susceptibility loci, two of these areas have previously been implicated in CL/P, viz. 2p13, an area harbouring the TGFA gene, and 6p23-24. We also demonstrate highly suggestive linkage to a susceptibility locus for nonsyndromic clefting on the X chromosome. Further studies are currently underway to replicate these findings in a larger cohort of affected sib-pairs.     

Evidence, from family studies, for linkage disequilibrium between TGFA and a gene for nonsyndromic cleft lip with or without cleft palate.

The inheritance of alleles of the transforming growth factor alpha (TGFA) locus has been studied in families affected with cleft lip with or without cleft palate (CL/P), by using the transmission/disequilibrium test described by Spielman and colleagues. Only heterozygous parents with an affected child can be included in this test, but within such families a significantly greater frequency of C2 alleles were transmitted to affected children than would be expected by chance. There was no evidence that the total number of C2 alleles transmitted to affected and unaffected children differed significantly from random segregation. These data provide evidence from within families that a gene for susceptibility to CL/P is in significant linkage disequilibrium with the C2 allele of the TGFA locus.     

Segregation analysis of cleft lip with or without cleft palate in the First Nations (Amerindian) people of British Columbia and review of isolated cleft palate etiologies

BACKGROUND: The First Nations (Amerindian) population of British Columbia, Canada, has the highest reported birth prevalence in the world of cleft lip with or without cleft palate (CL/P) at nearly 3 per 1000 births. In addition, a substantial proportion of cleft palate only (CPO) cases in this population has been reported to be X‐linked. The aims of this study were to perform complex segregation analysis to investigate the mode of inheritance of CL/P in the First Nations people of British Columbia and to review the etiology of the CPO cases. METHODS: All First Nations children born in British Columbia between 1952 and 1971 with an orofacial cleft were included in the study. Multiple sources of ascertainment were used, so that nearly 100% of live births were identified and included during this time. No stillbirths were found but would likely have been ascertained. Extended pedigrees were constructed from these probands and examination of immediate family members, e.g., parents and siblings, was done wherever possible. Complex segregation analysis included all family members. In addition, a CPO case review was conducted. RESULTS: Complex segregation analysis supports the hypothesis that the most likely mode of inheritance of CL/P in this population is a mixed model; that is, an autosomal major gene with polygenic component. The review of 26 CPO cases showed that a substantial proportion are syndromic. Birth Defects Research (Part A), 2009. © 2009 Wiley‐Liss, Inc.     

Correlation of susceptibility to 6-aminonicotinamide and hydrocortisone-induced cleft palate

Our previous genome-wide Quantitative Trait Locus (QTL) mapping study using mouse A/J by C57BL/6J recombinant inbred (RI) lines suggested several chromosomal regions contain genes influencing susceptibility to phenytoin (PT)-induced cleft lip with or without cleft palate [CL(P)] and 6-aminonicotinamide (6-AN)-induced isolated cleft palate (CP). Importantly, the same chromosomal regions but different RI parental strain alleles were sometimes implicated in susceptibility to these different kinds of orofacial clefting. Here we report the susceptibility to hydrocortisone (HC)-induced CP in these RI lines. We treated pregnant females with HC and studied the incidence of CP in day 17 fetuses. RI lines showed highly correlated responses to HC and 6-AN. The A/J parental line and five RI lines showed very high levels of clefting in response to both of these teratogens. The C57BL/6J parental line and five other RI lines exhibited low incidence of CP for these teratogens. In contrast, there was no significant correlation between incidence of PT-induced CL(P) and HC-induced CP.     

Cleft lip and palate: no evidence of linkage to transforming growth factor alpha.

Recently, an association of a specific transforming growth factor alpha (TGFA) haplotype has been reported to occur in patients with nonsyndromic cleft lip and palate (CLP) who had a positive family history of the disorder. We tested linkage of CLP to TGFA in 12 multiplex families with clefting. TGFA haplotypes were constructed on the basis of three polymorphic restriction sites. No haplotype was associated with CLP, and discordant inheritance of TGFA haplotypes was observed in a single informative family. The log odds (LOD) score was -2.1 at theta = .05. These findings indicate that CLP and TGFA are not tightly linked in these families.     

Evidence for an association between nonsyndromic cleft lip with or without cleft palate and a gene located on the long arm of chromosome 4.

Recent studies suggest that the familial aggregation of nonsyndromic cleft lip with or without cleft palate (CL +/- P) is likely to be attributable to the effects of several susceptibility loci, acting in a multiplicative fashion. Two potential CL +/- P susceptibility loci (CSL), transforming growth factor alpha (TGFA) and retinoic acid receptor (RARA), have been identified through association studies. In addition, recent evidence of linkage between CL +/- P and two markers (D4S175 and D4S192) in the region 4q25-4q31.3 raised the possibility that a CSL, with a larger effect than either TGFA or RARA, may reside within this region of the human genome. The present analyses were undertaken to determine whether D4S175 or D4S192 is significantly associated with CL +/- P in a sample of unrelated patients that have previously provided evidence of associations between CL +/- P and both TGFA and RARA. The results of these analyses provide further, tentative, evidence for the presence of a CSL locus on the long arm of chromosome 4 and help to refine the location of this locus in the region of D4S175 and D4S192.