A non-syndromic orofacial cleft risk locus links tRNA splicing defects to neural crest cell pathologies

Orofacial clefts are the most common form of congenital craniofacial malformation worldwide. The etiology of these birth defects is multifactorial, involving genetic and environmental factors. However, in most cases, the underlying causes remain unexplained, precluding a molecular understanding of disease mechanisms. Here, we integrated genome-wide association data, targeted resequencing of case and control cohorts, tissue- and cell-type-specific epigenomic profiling, and genome architecture analyses to molecularly dissect a genomic locus associated with an increased risk of non-syndromic orofacial cleft. We found that common and rare risk variants associated with orofacial cleft intersect with an enhancer (e2p24.2) that is active in human embryonic craniofacial tissue. We mapped e2p24.2 long-range interactions to a topologically associated domain harboring MYCN, DDX1, and CYRIA. We found that MYCN and DDX1, but not CYRIA, are required during craniofacial development in chicken embryos. We investigated the role of DDX1, a key component of the tRNA splicing complex, in cranial neural crest cells (cNCCs). The loss of DDX1 in cNCCs resulted in the accumulation of unspliced tRNA fragments, depletion of mature intron-containing tRNAs, and ribosome stalling at codons decoded by these tRNAs. This was accompanied by defects in both global protein synthesis and cNCC migration. We further showed that the induction of tRNA fragments is sufficient to disrupt craniofacial development. Together, these results uncovered a molecular mechanism in which impaired tRNA splicing affects cNCCs and craniofacial development and positioned MYCN, DDX1, and tRNA processing defects as risk factors in the pathogenesis of orofacial clefts.

Gene-by-environment interactions involving maternal exposures with orofacial cleft risk in Filipinos

Maternal exposures are known to influence the risk of isolated cleft lip with or without cleft palate (CL/P) - a common and highly heritable birth defect with a multifactorial etiology. To identify new CL/P risk loci, we conducted a genome-wide gene-environment interaction (GEI) analysis of CL/P on a sample of 540 cases and 260 controls recruited from the Philippines, incorporating the interaction effects of genetic variants with maternal smoking and vitamin use. As GEI analyses are typically low in power and the results can be difficult to interpret, we used multiple testing frameworks to evaluate potential GEI effects: 1 degree-of-freedom (1df) GxE test, the 3df joint test, and the two-step EDGE approach. While we did not detect any genome-wide significant interactions, we detected 12 suggestive GEI with smoking and 25 suggestive GEI with vitamin use between all testing frameworks. Several of these loci showed biological plausibility. Notable interactions with smoking include loci near FEZF1 , TWIST2, and NET1. While FEZF1 is involved in early neuronal development, TWIST2 and NET1 regulate epithelial-mesenchymal transition which is required for proper lip and palate fusion. Interactions with vitamins encompass CECR2 - a chromatin remodeling protein required for neural tube closure-and FURIN, a critical protease during early embryogenesis that activates various growth factor and extracellular-matrix protein. The activity of both proteins is influenced by folic acid. Our findings highlight the critical role of maternal exposures in identifying genes associated with structural birth defects such as CL/P and provide new paths to explore for CL/P genetics.

Combining genetic and single-cell expression data reveals cell types and novel candidate genes for orofacial clefting

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is one of the most common birth defects and has a multifactorial etiology. To date, over 45 loci harboring common risk variants have been identified. However, the effector genes at these loci, and the cell types that are affected by risk alleles, remain largely unknown. To address this, we combined genetic data from an nsCL/P genome-wide association study (GWAS) with single-cell RNA sequencing data obtained from the heads of unaffected human embryos. Using the recently developed single-cell disease relevance score (scDRS) approach, we identified two major cell types involved in nsCL/P development, namely the epithelium and the HAND2+ pharyngeal arches (PA). Combining scDRS with co-expression networks and differential gene expression analysis, we prioritized nsCL/P candidate genes, some of which were additionally supported by GWAS data (e.g., CTNND1, PRTG, RPL35A, RAB11FIP1, KRT19). Our results suggest that specific epithelial and PA sub-cell types are involved in nsCL/P development, and harbor a substantial fraction of the genetic risk for nsCL/P.

Retinoic Acid Upregulates METTL14 Expression and the m6A Modification Level to Inhibit the Proliferation of Embryonic Palate Mesenchymal Cells in Cleft Palate Mice

Cleft palate only (CPO) is one of the most common craniofacial birth defects. Environmental factors can induce cleft palate by affecting epigenetic modifications such as DNA methylation, histone acetylation, and non-coding RNA. However, there are few reports focusing on the RNA modifications. In this study, all-trans retinoic acid (atRA) was used to simulate environmental factors to induce a C57BL/6J fetal mouse cleft palate model. Techniques such as dot blotting and immunofluorescence were used to find the changes in m6A modification when cleft palate occurs. RNA-seq and KEGG analysis were used to screen for significantly differentially expressed pathways downstream. Primary mouse embryonic palate mesenchymal (MEPM) cells were successfully isolated and used for in vitro experimental verification. We found that an increased m6A methylation level was correlated with suppressed cell proliferation in the palatine process mesenchyme of cleft palate mice. This change is due to the abnormally high expression of m6A methyltransferase METTL14. When using siRNAs and the m6A methyltransferase complex inhibitor SAH to interfere with the expression or function of METTL14, the teratogenic effect of atRA on primary cells was partially alleviated. In conclusion, METTL14 regulates palatal mesenchymal cell proliferation and cycle-related protein expression relies on m6A methylation modification, affecting the occurrence of cleft palate.

Genome-wide meta-analyses identify five new risk loci for nonsyndromic orofacial clefts in the Chinese Han population

BACKGROUND

Nonsyndromic orofacial clefts (NSOFCs) are the most common craniofacial birth malformations in humans and are generally classified as nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Genome-wide association studies (GWASs) of NSOFCs have demonstrated multiple risk loci and candidate genes; however, published risk factors are able to explain only a small fraction of the observed NSOFCs heritability.

METHODS

Here, we performed GWASs of 1615 NSCPO cases and 2340 controls, and then conducted genome-wide meta-analyses of NSOFCs, totaling 6812 NSCL/P cases, 2614 NSCPO cases, and 19,165 controls from the Chinese Han population.

RESULTS

We identify 47 risk loci with genome-wide pmeta -value <5.0 × 10-8 , 5 risk loci (1p32.1, 3p14.1, 3p14.3, 3p21.31, and 13q22.1) of which are new. All of the 47 susceptibility loci conjointly account for 44.12% of the NSOFCs' heritability in the Chinese Han population.

CONCLUSION

Our results improve the comprehending of genetic susceptibility to NSOFCs and provide new views into the genetic etiology of craniofacial anomalies.

Identification of putative regulatory single-nucleotide variants in NTN1 gene associated with NSCL/P

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common polygenetic disease. Although genome-wide association studies (GWAS) identified NTN1 gene as a high-priority candidate of NSCL/P, the comprehensive genetic architecture of NTN1 weren't yet known. Thus, this study aimed to determine full-scale genetic variants of NTN1 for NSCL/P in Chinese Han people. Initially, targeted sequencing of NTN1 gene was performed on 159 NSCL/P patients to identify susceptible single nucleotide polymorphisms (SNPs) associated with NSCL/P. Then, association analysis and burden analysis were separately used to validate the common variants and rare variants identified among large size of samples (1608 NSCL/P cases and 2255 controls). Additionally, NSCL/P subtype association analysis was applied to elucidate the etiology discrepancy of non-syndromic cleft lip with palate (NSCLP) and non-syndromic cleft lip only (NSCLO). Lastly, bioinformatics analysis was performed to annotate and prioritize candidate variants. We found 15 NSCL/P-associated SNPs including rs4791774 (P = 1.10E-08, OR = 1.467, 95% CI: 1.286~1.673) and rs9788972 (P = 1.28E-07, OR = 1.398, 95% CI : 1.235~1.584) originally detected by previous GWASs in Chinese Han ancestry. Four NSCLO risk-associated SNPs and eight specific NSCLP associated SNPs were found. Three SNPs (rs4791331, rs4791774 and rs9900753) were predicted to locate at regulatory region of NTN1. Our study validated the association between NTN1 gene and pathogenesis of NSCL/P and reinforced the hypothesis that NSCLP have a different etiology from NSCLO. We also identified three putative regulatory SNPs in NTN1 gene.

Identification of a new familial case of 3q29 deletion syndrome associated with cleft lip and palate via whole-exome sequencing

Many unbalanced large copy number variants reviewed in the paper are associated with syndromic orofacial clefts, including a 1.6 Mb deletion on chromosome 3q29. The current report presents a new family with this recurrent deletion identified via whole-exome sequencing and confirmed by array comparative genomic hybridization. The proband exhibited a more severe clinical phenotype than his affected mother, comprising right-sided cleft lip/alveolus and cleft palate, advanced dental caries, heart defect, hypospadias, psychomotor, and speech delay, and an intellectual disability. Data analysis from the 3q29 registry revealed that the 3q29 deletion increases the risk of clefting by nearly 30-fold. No additional rare and pathogenic nucleotide variants were identified that could explain the clefting phenotype and observed intrafamilial phenotypic heterogeneity. These data suggest that the 3q29 deletion may be the primary risk factor for clefting, with additional genomic variants located outside the coding sequences, methylation changes, or environmental exposure serving as modifiers of this risk. Additional studies, including whole-genome sequencing or methylation analyses, should be performed to identify genetic factors underlying the phenotypic variation associated with the recurrent 3q29 deletion.

Decoding the Human Face: Progress and Challenges in Understanding the Genetics of Craniofacial Morphology

Variations in the form of the human face, which plays a role in our individual identities and societal interactions, have fascinated scientists and artists alike. Here, we review our current understanding of the genetics underlying variation in craniofacial morphology and disease-associated dysmorphology, synthesizing decades of progress on Mendelian syndromes in addition to more recent results from genome-wide association studies of human facial shape and disease risk. We also discuss the various approaches used to phenotype and quantify facial shape, which are of particular importance due to the complex, multipartite nature of the craniofacial form. We close by discussing how experimental studies have contributed and will further contribute to our understanding of human genetic variation and then proposing future directions and applications for the field.

Identification of novel susceptibility genes for non-syndromic cleft lip with or without cleft palate using NGS-based multigene panel testing

For non-syndromic cleft lip with or without cleft palate (ns-CL/P), the proportion of heritability explained by the known risk loci is estimated to be about 30% and is captured mainly by common variants identified in genome-wide association studies. To contribute to the explanation of the "missing heritability" problem for orofacial clefts, a candidate gene approach was taken to investigate the potential role of rare and private variants in the ns-CL/P risk. Using the next-generation sequencing technology, the coding sequence of a set of 423 candidate genes was analysed in 135 patients from the Polish population. After stringent multistage filtering, 37 rare coding and splicing variants of 28 genes were identified. 35% of these genetic alternations that may play a role of genetic modifiers influencing an individual's risk were detected in genes not previously associated with the ns-CL/P susceptibility, including COL11A1, COL17A1, DLX1, EFTUD2, FGF4, FGF8, FLNB, JAG1, NOTCH2, SHH, WNT5A and WNT9A. Significant enrichment of rare alleles in ns-CL/P patients compared with controls was also demonstrated for ARHGAP29, CHD7, COL17A1, FGF12, GAD1 and SATB2. In addition, analysis of panoramic radiographs of patients with identified predisposing variants may support the hypothesis of a common genetic link between orofacial clefts and dental abnormalities. In conclusion, our study has confirmed that rare coding variants might contribute to the genetic architecture of ns-CL/P. Since only single predisposing variants were identified in novel cleft susceptibility genes, future research will be required to confirm and fully understand their role in the aetiology of ns-CL/P.

Allele-specific transcription factor binding in a cellular model of orofacial clefting

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is a frequent congenital malformation with multifactorial etiology. While recent genome-wide association studies (GWAS) have identified several nsCL/P risk loci, the functional effects of the associated non-coding variants are largely unknown. Furthermore, additional risk loci remain undetected due to lack of power. As genetic variants might alter binding of transcription factors (TF), we here hypothesized that the integration of data from TF binding sites, expression analyses and nsCL/P GWAS might help to (i) identify functionally relevant variants at GWAS loci, and (ii) highlight novel risk variants that have been previously undetected. Analysing the craniofacial TF TFAP2A in human embryonic palatal mesenchyme (HEPM) cells, we identified 2845 TFAP2A ChIP-seq peaks, several of which were located near nsCL/P candidate genes (e.g. MSX1 and SPRY2). Comparison with independent data suggest that 802 of them might be specific to craniofacial development, and genes near these peaks are enriched in processes relevant to nsCL/P. Integration with nsCL/P GWAS data, however, did not show robust evidence for co-localization of common nsCL/P risk variants with TFAP2A ChIP-seq peaks. This data set represents a new resource for the analyses of craniofacial processes, and similar approaches with additional cell lines and TFs could be applied to generate further insights into nsCL/P etiology.

Mycn deficiency underlies the development of orofacial clefts in mice and humans

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is the most common subphenotype of non-syndromic orofacial clefts (NSOFCs) arising from genetic and/or environmental perturbations during embryonic development. We previously identified 2p24.2 as a risk locus associated with NSCL/P in the Chinese Han population, and MYCN is a candidate risk gene in this region. To understand the potential function of MYCN in craniofacial development, we generated Wnt1-Cre;Mycnflox/flox mice that exhibited cleft palate, microglossia, and micrognathia, resembling the Pierre Robin sequence (PRS) in humans. Further analyses indicated that the cleft palate was secondary to the delayed elevation of palatal shelves caused by micrognathia. The micrognathia resulted from impaired chondrogenic differentiation in Merkel's cartilage, which limited tongue development, leading to microglossia. In terms of mechanism, Mycn deficiency in cranial neural crest cells (CNCCs) downregulated Sox9 expression by inhibiting Wnt5a in a CNCC-derived chondrogenic lineage in Merkel's cartilage. To investigate whether MYCN deficiency contributed to NSCL/P, we performed direct sequencing targeting all exons and exon-intron boundaries of MYCN in 104 multiplex families with mendelian NSCL/P and identified a novel pathogenic variant in MYCN. Taken together, our data indicate that ablation of Mycn in mouse CNCCs could resemble PRS by suppressing the Wnt5a-Sox9 signaling pathway in Merkel's cartilage and that mutations in MYCN may be novel potential causes of NSCL/P.

MiRNA-149 as a Candidate for Facial Clefting and Neural Crest Cell Migration

Nonsyndromic cleft lip with or without palate (nsCL/P) ranks among the most common human birth defects and has a multifactorial etiology. Human neural crest cells (hNCC) make a substantial contribution to the formation of facial bone and cartilage and are a key cell type in terms of nsCL/P etiology. Based on increasing evidence for the role of noncoding regulatory mechanisms in nsCL/P, we investigated the role of hNCC-expressed microRNAs (miRNA) in cleft development. First, we conducted a systematic analysis of miRNAs expressed in human-induced pluripotent stem cell-derived hNCC using Affymetrix microarrays on cell lines established from 4 unaffected donors. These analyses identified 152 candidate miRNAs. Based on the hypothesis that candidate miRNA loci harbor genetic variation associated with nsCL/P risk, the genomic locations of these candidates were cross-referenced with data from a previous genome-wide association study of nsCL/P. Associated variants were reanalyzed in independent nsCL/P study populations. Jointly, the results suggest that miR-149 is implicated in nsCL/P etiology. Second, functional follow-up included in vitro overexpression and inhibition of miR-149 in hNCC and subsequent analyses at the molecular and phenotypic level. Using 3'RNA-Seq, we identified 604 differentially expressed (DE) genes in hNCC overexpressing miR-149 compared with untreated cells. These included TLR4 and JUNB, which are established targets of miR-149, and NOG, BMP4, and PAX6, which are reported nsCL/P candidate genes. Pathway analyses revealed that DE genes were enriched in pathways including regulation of cartilage development and NCC differentiation. At the cellular level, distinct hNCC migration patterns were observed in response to miR-149 overexpression. Our data suggest that miR-149 is involved in the etiology of nsCL/P via its role in hNCC migration.

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

Extending the allelic spectrum at noncoding risk loci of orofacial clefting

Genome-wide association studies (GWAS) have generated unprecedented insights into the genetic etiology of orofacial clefting (OFC). The moderate effect sizes of associated noncoding risk variants and limited access to disease-relevant tissue represent considerable challenges for biological interpretation of genetic findings. As rare variants with stronger effect sizes are likely to also contribute to OFC, an alternative approach to delineate pathogenic mechanisms is to identify private mutations and/or an increased burden of rare variants in associated regions. This report describes a framework for targeted resequencing at selected noncoding risk loci contributing to nonsyndromic cleft lip with/without cleft palate (nsCL/P), the most frequent OFC subtype. Based on GWAS data, we selected three risk loci and identified candidate regulatory regions (CRRs) through the integration of credible SNP information, epigenetic data from relevant cells/tissues, and conservation scores. The CRRs (total 57 kb) were resequenced in a multiethnic study population (1061 patients; 1591 controls), using single-molecule molecular inversion probe technology. Combining evidence from in silico variant annotation, pedigree- and burden analyses, we identified 16 likely deleterious rare variants that represent new candidates for functional studies in nsCL/P. Our framework is scalable and represents a promising approach to the investigation of additional congenital malformations with multifactorial etiology.

The Intersection of the Genetic Architectures of Orofacial Clefts and Normal Facial Variation

Unaffected relatives of individuals with non-syndromic cleft lip with or without cleft palate (NSCL/P) show distinctive facial features. The presence of this facial endophenotype is potentially an expression of underlying genetic susceptibility to NSCL/P in the larger unselected population. To explore this hypothesis, we first partitioned the face into 63 partially overlapping regions representing global-to-local facial morphology and then defined endophenotypic traits by contrasting the 3D facial images from 264 unaffected parents of individuals with NSCL/P versus 3,171 controls. We observed distinct facial features between parents and controls across 59 global-to-local facial segments at nominal significance (p ≤ 0.05) and 52 segments at Bonferroni corrected significance (p < 1.2 × 10-3), respectively. Next, we quantified these distinct facial features as univariate traits in another dataset of 8,246 unaffected European individuals and performed a genome-wide association study. We identified 29 independent genetic loci that were associated (p < 5 × 10-8) with at least one of the tested endophenotypic traits, and nine genetic loci also passed the study-wide threshold (p < 8.47 × 10-10). Of the 29 loci, 22 were in proximity of loci previously associated with normal facial variation, 18 were near genes that show strong evidence in orofacial clefting (OFC), and another 10 showed some evidence in OFC. Additionally, polygenic risk scores for NSCL/P showed associations with the endophenotypic traits. This study thus supports the hypothesis of a shared genetic architecture of normal facial development and OFC.

A clinical and multi‑omics study of Van der Woude syndrome in three generations of a Chinese family

Previous studies have suggested that pathogenic variants in interferon regulatoryse factor 6 (IRF6) can account for almost 70% of familial Van der Woude Syndrome (VWS) cases. However, gene modifiers that account for the phenotypic variability of IRF6 in the context of VWS remain poorly characterized. The aim of this study was to report a family with VWS with variable expressivity and to identify the genetic cause. A 4-month-old boy initially presented with cleft palate and bilateral lower lip pits. Examination of his family history identified similar, albeit milder, clinical features in another four family members, including bilateral lower lip pits and/or hypodontia. Peripheral blood samples of eight members in this three-generation family were subsequently collected, and whole-exome sequencing was performed to detect pathogenic variants. A heterozygous missense IRF6 variant with a c.1198C>T change in exon 9 (resulting in an R400W change at the amino acid level) was detected in five affected subjects, but not in the other three unaffected subjects. Moreover, subsequent structural analysis was indicative of damaged stability to the structure in the mutant IRF protein. Whole-transcriptome sequencing, expression analysis and Gene Ontology enrichment analysis were conducted on two groups of patients with phenotypic diversity from the same family. These analyses identified significant differentially expressed genes and enriched pathways in these two groups. Altogether, these findings provide insight into the mechanism underlying the variable expressivity of VWS.

Evidence for craniofacial enhancer variation underlying nonsyndromic cleft lip and palate

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect for which only ~ 20% of the underlying genetic variation has been identified. Variants in noncoding regions have been increasingly suggested to contribute to the missing heritability. In this study, we investigated whether variation in craniofacial enhancers contributes to NSCLP. Candidate enhancers were identified using VISTA Enhancer Browser and previous publications. Prioritization was based on patterning defects in knockout mice, deletion/duplication of craniofacial genes in animal models and results of whole exome/whole genome sequencing studies. This resulted in 20 craniofacial enhancers to be investigated. Custom amplicon-based sequencing probes were designed and used for sequencing 380 NSCLP probands (from multiplex and simplex families of non-Hispanic white (NHW) and Hispanic ethnicities) using Illumina MiSeq. The frequencies of identified variants were compared to ethnically matched European (CEU) and Los Angeles Mexican (MXL) control genomes and used for association analyses. Variants in mm427/MSX1 and hs1582/SPRY1 showed genome-wide significant association with NSCLP (p ≤ 6.4 × 10^-11). In silico analysis showed that these enhancer variants may disrupt important transcription factor binding sites. Haplotypes involving these enhancers and also mm435/ABCA4 were significantly associated with NSCLP, especially in NHW (p ≤ 6.3 × 10^-7). Importantly, groupwise burden analysis showed several enhancer combinations significantly over-represented in NSCLP individuals, revealing novel NSCLP pathways and supporting a polygenic inheritance model. Our findings support the role of craniofacial enhancer sequence variation in the etiology of NSCLP.

A Retrospective Mixed Longitudinal Study of Tooth Formation in Children With Clefts

OBJECTIVE

To test for systematic age changes in cleft children based on dental age.

DESIGN

Retrospective case-control longitudinal study.

SETTING

One orthodontic solo practice.

PATIENTS

Nonsyndromic, complete cleft lip and palateCLP cases, either unilateral or bilateral (102 children; 370 radiographs), between 4 and 16 years of age.

INTERVENTIONS

Children were treated with a team approach, but only orthodontic radiographs were studied.

MAIN OUTCOME MEASURE

The principal outcome measure was dental age of the cleft cases compared to a sex-specific sample of phenotypically normal children (1107 children), from the same geographical region. Multiple panoramic radiographs taken during the course of orthodontic treatment were examined to track patterns of dental age as children matured. Analysis used linear mixed models primarily testing for sex, cleft type (unilateral, bilateral), and hypodontia differences. Initial expectation was that cleft children would exhibit delayed dental ages from postnatal stressors and would become more deviant with maturity.

RESULTS

In childhood (4-6 years), both sexes were significantly delayed (P < .001), but dental age normalized around 8 to 10 years. Boys experienced faster maturation thereafter than girls (P < .001). Only trivial differences occurred between unilateral CLP and bilateral CLP samples.Hypodontia further depressed maturation rates (P < .001). Dental age improved in a decidedly curvilinear fashion (P < .001), with greater change at earlier ages.

CONCLUSIONS

This report agrees with other contemporary studies, showing childhood catch-up. Older studies observed that clefting caused significant delays that worsened with growth. This potential "seachange" suggests better recovery and quicker normalization of children with clefts, perhaps due to improved management.

High-Resolution Epigenomic Atlas of Human Embryonic Craniofacial Development

Defects in patterning during human embryonic development frequently result in craniofacial abnormalities. The gene regulatory programs that build the craniofacial complex are likely controlled by information located between genes and within intronic sequences. However, systematic identification of regulatory sequences important for forming the human face has not been performed. Here, we describe comprehensive epigenomic annotations from human embryonic craniofacial tissues and systematic comparisons with multiple tissues and cell types. We identified thousands of tissue-specific craniofacial regulatory sequences and likely causal regions for rare craniofacial abnormalities. We demonstrate significant enrichment of common variants associated with orofacial clefting in enhancers active early in embryonic development, while those associated with normal facial variation are enriched near the end of the embryonic period. These data are provided in easily accessible formats for both craniofacial researchers and clinicians to aid future experimental design and interpretation of noncoding variation in those affected by craniofacial abnormalities.

Disruption of FOXF2 as a Likely Cause of Absent Uvula in an Egyptian Family

This study investigated the genetic basis of an unusual autosomal dominant phenotype characterized by familial absent uvula, with a short posterior border of the soft palate, abnormal tonsillar pillars, and velopharyngeal insufficiency. Cytogenetic analysis and single-nucleotide polymorphism-based linkage analysis were investigated in a 4-generation family with 8 affected individuals. Whole exome sequencing data were overlaid, and segregation analysis identified a single missense variant, p.Q433P in the FOXF2 transcription factor, that fully segregated with the phenotype. This was found to be in linkage disequilibrium with a small 6p25.3 tandem duplication affecting FOXC1 and GMDS. Notably, the copy number imbalances of this region are commonly associated with pathologies that are not present in this family. Bioinformatic predictions with luciferase reporter studies of the FOXF2 missense variant indicated a negative impact, affecting both protein stability and transcriptional activation. Foxf 2 is expressed in the posterior mouse palate, and knockout animals develop an overt cleft palate. Since mice naturally lack the structural equivalent of the uvula, we demonstrated FOXF2 expression in the developing human uvula. Decipher also records 2 individuals with hypoplastic or bifid uvulae with copy number variants affecting FOXF2. Nevertheless, given cosegregation with the 6p25.3 duplications, we cannot rule out a combined effect of these gains and the missense variant on FOXF2 function, which may account for the rare palate phenotype observed.

Nonsyndromic cleft palate: An association study at GWAS candidate loci in a multiethnic sample

BACKGROUND

Nonsyndromic cleft palate only (nsCPO) is a common and multifactorial form of orofacial clefting. In contrast to successes achieved for the other common form of orofacial clefting, that is, nonsyndromic cleft lip with/without cleft palate (nsCL/P), genome wide association studies (GWAS) of nsCPO have identified only one genome wide significant locus. Aim of the present study was to investigate whether common variants contribute to nsCPO and, if so, to identify novel risk loci.

METHODS

We genotyped 33 SNPs at 27 candidate loci from 2 previously published nsCPO GWAS in an independent multiethnic sample. It included: (i) a family-based sample of European ancestry (n = 212); and (ii) two case/control samples of Central European (n = 94/339) and Arabian ancestry (n = 38/231), respectively. A separate association analysis was performed for each genotyped dataset, and meta-analyses were performed.

RESULTS

After association analysis and meta-analyses, none of the 33 SNPs showed genome-wide significance. Two variants showed nominally significant association in the imputed GWAS dataset and exhibited a further decrease in p-value in a European and an overall meta-analysis including imputed GWAS data, respectively (rs395572: P_MetaEU  = 3.16 × 10^-4 ; rs6809420: P_MetaAll  = 2.80 × 10^-4 ).

CONCLUSION

Our findings suggest that there is a limited contribution of common variants to nsCPO. However, the individual effect sizes might be too small for detection of further associations in the present sample sizes. Rare variants may play a more substantial role in nsCPO than in nsCL/P, for which GWAS of smaller sample sizes have identified genome-wide significant loci. Whole-exome/genome sequencing studies of nsCPO are now warranted.