Molecular mechanisms of midfacial developmental defects

Single-cell transcriptome and chromatin accessibility mapping of upper lip and primary palate fusion

Cleft lip and/or primary palate (CL/P) represent a prevalent congenital malformation, the aetiology of which is highly intricate. Although it is generally accepted that the condition arises from failed fusion between the upper lip and primary palate, the precise mechanism underlying this fusion process remains enigmatic. In this study, we utilized transposase-accessible chromatin sequencing (scATAC-seq) and single-cell RNA sequencing (scRNA-seq) to interrogate lambdoidal junction tissue derived from C57BL/6J mouse embryos at critical stages of embryogenesis (10.5, 11.5 and 12.5 embryonic days). We successfully identified distinct subgroups of mesenchymal and ectodermal cells involved in the fusion process and characterized their unique transcriptional profiles. Furthermore, we conducted cell differentiation trajectory analysis, revealing a dynamic repertoire of genes that are sequentially activated or repressed during pseudotime, facilitating the transition of relevant cell types. Additionally, we employed scATAC data to identify key genes associated with the fusion process and demonstrated differential chromatin accessibility across major cell types. Finally, we constructed a dynamic intercellular communication network and predicted upstream transcriptional regulators of critical genes involved in important signalling pathways. Our findings provide a valuable resource for future studies on upper lip and primary palate development, as well as congenital defects.

DeepFace: Deep-learning-based framework to contextualize orofacial-cleft-related variants during human embryonic craniofacial development

Orofacial clefts (OFCs) are among the most common human congenital birth defects. Previous multiethnic studies have identified dozens of associated loci for both cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP). Although several nearby genes have been highlighted, the "casual" variants are largely unknown. Here, we developed DeepFace, a convolutional neural network model, to assess the functional impact of variants by SNP activity difference (SAD) scores. The DeepFace model is trained with 204 epigenomic assays from crucial human embryonic craniofacial developmental stages of post-conception week (pcw) 4 to pcw 10. The Pearson correlation coefficient between the predicted and actual values for 12 epigenetic features achieved a median range of 0.50-0.83. Specifically, our model revealed that SNPs significantly associated with OFCs tended to exhibit higher SAD scores across various variant categories compared to less related groups, indicating a context-specific impact of OFC-related SNPs. Notably, we identified six SNPs with a significant linear relationship to SAD scores throughout developmental progression, suggesting that these SNPs could play a temporal regulatory role. Furthermore, our cell-type specificity analysis pinpointed the trophoblast cell as having the highest enrichment of risk signals associated with OFCs. Overall, DeepFace can harness distal regulatory signals from extensive epigenomic assays, offering new perspectives for prioritizing OFC variants using contextualized functional genomic features. We expect DeepFace to be instrumental in accessing and predicting the regulatory roles of variants associated with OFCs, and the model can be extended to study other complex diseases or traits.

Dyrk1a is required for craniofacial development in Xenopus laevis

Loss of function variations in the dual specificity tyrosine-phosphorylation-regulated kinase 1 A (DYRK1A) gene are associated with craniofacial malformations in humans. Here we characterized the effects of deficient DYRK1A in craniofacial development using a developmental model, Xenopus laevis. Dyrk1a mRNA and protein were expressed throughout the developing head and both were enriched in the branchial arches which contribute to the face and jaw. Consistently, reduced Dyrk1a function, using dyrk1a morpholinos and pharmacological inhibitors, resulted in orofacial malformations including hypotelorism, altered mouth shape, slanted eyes, and narrower face accompanied by smaller jaw cartilage and muscle. Inhibition of Dyrk1a function resulted in misexpression of key craniofacial regulators including transcription factors and members of the retinoic acid signaling pathway. Two such regulators, sox9 and pax3 are required for neural crest development and their decreased expression corresponds with smaller neural crest domains within the branchial arches. Finally, we determined that the smaller size of the faces, jaw elements and neural crest domains in embryos deficient in Dyrk1a could be explained by increased cell death and decreased proliferation. This study is the first to provide insight into why craniofacial birth defects might arise in humans with variants of DYRK1A.

Genetic mechanism regulating diversity in the placement of eyes on the head of animals

Despite the conservation of genetic machinery involved in eye development, there is a strong diversity in the placement of eyes on the head of animals. Morphogen gradients of signaling molecules are vital to patterning cues. During Drosophila eye development, Wingless (Wg), a ligand of Wnt/Wg signaling, is expressed anterolaterally to form a morphogen gradient to determine the eye- versus head-specific cell fate. The underlying mechanisms that regulate this process are yet to be fully understood. We characterized defective proventriculus (dve) (Drosophila ortholog of human SATB1), a K50 homeodomain transcription factor, as a dorsal eye gene, which regulates Wg signaling to determine eye versus head fate. Across Drosophila species, Dve is expressed in the dorsal head vertex region where it regulates wg transcription. Second, Dve suppresses eye fate by down-regulating retinal determination genes. Third, the dve-expressing dorsal head vertex region is important for Wg-mediated inhibition of retinal cell fate, as eliminating the Dve-expressing cells or preventing Wg transport from these dve-expressing cells leads to a dramatic expansion of the eye field. Together, these findings suggest that Dve regulates Wg expression in the dorsal head vertex, which is critical for determining eye versus head fate. Gain-of-function of SATB1 exhibits an eye fate suppression phenotype similar to Dve. Our data demonstrate a conserved role for Dve/SATB1 in the positioning of eyes on the head and the interocular distance by regulating Wg. This study provides evidence that dysregulation of the Wg morphogen gradient results in developmental defects such as hypertelorism in humans where disproportionate interocular distance and facial anomalies are reported.

A BMP-controlled metabolic/epigenetic signaling cascade directs midfacial morphogenesis

Craniofacial anomalies, especially midline facial defects, are among the most common birth defects in patients and are associated with increased mortality or require lifelong treatment. During mammalian embryogenesis, specific instructions arising at genetic, signaling, and metabolic levels are important for stem cell behaviors and fate determination, but how these functionally relevant mechanisms are coordinated to regulate craniofacial morphogenesis remain unknown. Here, we report that bone morphogenetic protein (BMP) signaling in cranial neural crest cells (CNCCs) is critical for glycolytic lactate production and subsequent epigenetic histone lactylation, thereby dictating craniofacial morphogenesis. Elevated BMP signaling in CNCCs through constitutively activated ACVR1 (ca-ACVR1) suppressed glycolytic activity and blocked lactate production via a p53-dependent process that resulted in severe midline facial defects. By modulating epigenetic remodeling, BMP signaling-dependent lactate generation drove histone lactylation levels to alter essential genes of Pdgfra, thus regulating CNCC behavior in vitro as well as in vivo. These findings define an axis wherein BMP signaling controls a metabolic/epigenetic cascade to direct craniofacial morphogenesis, thus providing a conceptual framework for understanding the interaction between genetic and metabolic cues operative during embryonic development. These findings indicate potential preventive strategies of congenital craniofacial birth defects via modulating metabolic-driven histone lactylation.

Single-cell multiomics decodes regulatory programs for mouse secondary palate development

Perturbations in gene regulation during palatogenesis can lead to cleft palate, which is among the most common congenital birth defects. Here, we perform single-cell multiome sequencing and profile chromatin accessibility and gene expression simultaneously within the same cells (n = 36,154) isolated from mouse secondary palate across embryonic days (E) 12.5, E13.5, E14.0, and E14.5. We construct five trajectories representing continuous differentiation of cranial neural crest-derived multipotent cells into distinct lineages. By linking open chromatin signals to gene expression changes, we characterize the underlying lineage-determining transcription factors. In silico perturbation analysis identifies transcription factors SHOX2 and MEOX2 as important regulators of the development of the anterior and posterior palate, respectively. In conclusion, our study charts epigenetic and transcriptional dynamics in palatogenesis, serving as a valuable resource for further cleft palate research.

Dyrk1a is required for craniofacial development in Xenopus laevis

Loss of function mutations in the dual specificity tyrosine-phosphorylation-regulated kinase 1A (DYRK1A) gene are associated with craniofacial malformations in humans. Here we characterized the effects of deficient DYRK1A in craniofacial development using a developmental model, Xenopus laevis . Dyrk1a mRNA and protein was expressed throughout the developing head and was enriched in the branchial arches which contribute to the face and jaw. Consistently, reduced Dyrk1a function, using dyrk1a morpholinos and pharmacological inhibitors, resulted in orofacial malformations including hypotelorism, altered mouth shape, slanted eyes, and narrower face accompanied by smaller jaw cartilage and muscle. Inhibition of Dyrk1a function resulted in misexpression of key craniofacial regulators including transcription factors and members of the retinoic acid signaling pathway. Two such regulators, sox9 and pax3 are required for neural crest development and their decreased expression corresponds with smaller neural crest domains within the branchial arches. Finally, we determined that the smaller size of the faces, jaw elements and neural crest domains in embryos deficient in Dyrk1a could be explained by increased cell death and decreased proliferation. This study is the first to provide insight into why craniofacial birth defects might arise in humans with DYRK1A mutations.

Sagittal Growth Restriction of the Midface Following Isolated Cleft Lip Repair: A Systematic Review and Meta-Analysis

Midface hypoplasia (MFH) is a long-term sequela of cleft lip and palate repair, and is poorly understood. No study has examined the aggregate data on sagittal growth restriction of the midface following repair of the lip, but not palate, in these patients.

A systematic review of 3780 articles was performed. Twenty-four studies met inclusion criteria and 11 reported cephalometric measurements amenable to meta-analysis. Patients with Veau class I-III palatal clefts were included so long as they had undergone only lip repair. Groups were compared against both noncleft and unrepaired controls.

Cephalometrics were reported for 326 patients (31.3% female). Noncleft controls had an average SNA angle of 81.25° ± 3.12°. The only patients demonstrating hypoplastic SNA angles were those with unilateral CLP with isolated lip repair (77.4° ± 4.22°). Patients with repaired CL had SNA angles similar to noncleft controls (81.4° ± 4.02°). Patients with unrepaired CLP and CL tended toward more protruding maxillae, with SNA angles of 83.3° ± 4.04° and 87.9° ± 3.11°, respectively. Notably, when comparing SNA angles between groups, patients with CLP with isolated lip repair had significantly more hypoplastic angles compared to those with repaired CL (P < .0001). Patients with CLP with isolated lip repair were also more hypoplastic than noncleft controls (P < .0001). In contrast, there was no significant difference between the SNA of patients with repaired CL and controls (P = .648).

We found that cleft lip repair only appeared to contribute to MFH in the setting of concurrent cleft palate pathology, suggesting that scarring from lip repair itself is unlikely to be the predominant driver of MFH development. However, studies generally suffered from inadequate reporting of timing, technique, follow-up time, and cleft severity.

A non-coding insertional mutation of Grhl2 causes gene over-expression and multiple structural anomalies including cleft palate, spina bifida and encephalocele

Orofacial clefts, including cleft lip and palate (CL/P) and neural tube defects (NTDs) are among the most common congenital anomalies, but knowledge of the genetic basis of these conditions remains incomplete. The extent to which genetic risk factors are shared between CL/P, NTDs and related anomalies is also unclear. While identification of causative genes has largely focused on coding and loss of function mutations, it is hypothesized that regulatory mutations account for a portion of the unidentified heritability. We found that excess expression of Grainyhead-like 2 (Grhl2) causes not only spinal NTDs in Axial defects (Axd) mice but also multiple additional defects affecting the cranial region. These include orofacial clefts comprising midline cleft lip and palate and abnormalities of the craniofacial bones and frontal and/or basal encephalocele, in which brain tissue herniates through the cranium or into the nasal cavity. To investigate the causative mutation in the Grhl2Axd strain, whole genome sequencing identified an approximately 4 kb LTR retrotransposon insertion that disrupts the non-coding regulatory region, lying approximately 300 base pairs upstream of the 5' UTR. This insertion also lies within a predicted long non-coding RNA, oriented on the reverse strand, which like Grhl2 is over-expressed in Axd (Grhl2Axd) homozygous mutant embryos. Initial analysis of the GRHL2 upstream region in individuals with NTDs or cleft palate revealed rare or novel variants in a small number of cases. We hypothesize that mutations affecting the regulation of GRHL2 may contribute to craniofacial anomalies and NTDs in humans.

Analysis of candidate genes for cleft lip ± cleft palate using murine single-cell expression data

Introduction: Cleft lip ± cleft palate (CL/P) is one of the most common birth defects. Although research has identified multiple genetic risk loci for different types of CL/P (i.e., syndromic or non-syndromic forms), determining the respective causal genes and understanding the relevant functional networks remain challenging. The recent introduction of single-cell RNA sequencing (scRNA-seq) has provided novel opportunities to study gene expression patterns at cellular resolution. The aims of our study were to: (i) aggregate available scRNA-seq data from embryonic mice and provide this as a resource for the craniofacial community; and (ii) demonstrate the value of these data in terms of the investigation of the gene expression patterns of CL/P candidate genes. Methods and Results: First, two published scRNA-seq data sets from embryonic mice were re-processed, i.e., data representing the murine time period of craniofacial development: (i) facial data from embryonic day (E) E11.5; and (ii) whole embryo data from E9.5-E13.5 from the Mouse Organogenesis Cell Atlas (MOCA). Marker gene expression analyses demonstrated that at E11.5, the facial data were a high-resolution representation of the MOCA data. Using CL/P candidate gene lists, distinct groups of genes with specific expression patterns were identified. Among others we identified that a co-expression network including Irf6, Grhl3 and Tfap2a in the periderm, while it was limited to Irf6 and Tfap2a in palatal epithelia, cells of the ectodermal surface, and basal cells at the fusion zone. The analyses also demonstrated that additional CL/P candidate genes (e.g., Tpm1, Arid3b, Ctnnd1, and Wnt3) were exclusively expressed in Irf6+ facial epithelial cells (i.e., as opposed to Irf6- epithelial cells). The MOCA data set was finally used to investigate differences in expression profiles for candidate genes underlying different types of CL/P. These analyses showed that syndromic CL/P genes (syCL/P) were expressed in significantly more cell types than non-syndromic CL/P candidate genes (nsCL/P). Discussion: The present study illustrates how scRNA-seq data can empower research on craniofacial development and disease.

Spatiotemporal MicroRNA-Gene Expression Network Related to Orofacial Clefts

Craniofacial structures change dynamically in morphology during development through the coordinated regulation of various cellular molecules. However, it remains unclear how these complex mechanisms are regulated in a spatiotemporal manner. Here we applied natural cubic splines to model gene and microRNA (miRNA) expression from embryonic day (E) 10.5 to E14.5 in the proximal and distal regions of the maxillary processes to identify spatiotemporal patterns of gene and miRNA expression, followed by constructing corresponding regulatory networks. Three major groups of differentially expressed genes (DEGs) were identified, including 3,927 temporal, 314 spatial, and 494 spatiotemporal DEGs. Unsupervised clustering further resolved these spatiotemporal DEGs into 8 clusters with distinct expression patterns. Interestingly, we found 2 clusters of differentially expressed miRNAs: 1 had 80 miRNAs monotonically decreasing and the other had 97 increasing across developmental stages. To evaluate the phenotypic relevance of these DEGs during craniofacial development, we integrated data from the CleftGeneDB database and constructed the regulatory networks of genes related to orofacial clefts. Our analysis revealed 2 hub miRNAs, mmu-miR-325-3p and mmu-miR-384-5p, that repressed cleft-related genes Adamts3, Runx2, Fgfr2, Acvr1, and Edn2, while their expression increased over time. On the contrary, 2 hub miRNAs, mmu-miR-218-5p and mmu-miR-338-5p, repressed cleft-related genes Pbx2, Ermp1, Snai1, Tbx2, and Bmi1, while their expression decreased over time. Our experiments indicated that these miRNA mimics significantly inhibited cell proliferation in mouse embryonic palatal mesenchymal (MEPM) cells and O9-1 cells through the regulation of genes associated with cleft palate and validated the role of our regulatory networks in orofacial clefts. To facilitate interactive exploration of these data, we developed a user-friendly web tool to visualize the gene and miRNA expression patterns across developmental stages, as well as the regulatory networks (https://fyan.shinyapps.io/facebase_shiny/). Taken together, our results provide a valuable resource that serves as a reference map for future research in craniofacial development.

FaceBase: A Community-Driven Hub for Data-Intensive Research

The FaceBase Consortium, funded by the National Institute of Dental and Craniofacial Research of the National Institutes of Health, was established in 2009 with the recognition that dental and craniofacial research are increasingly data-intensive disciplines. Data sharing is critical for the validation and reproducibility of results as well as to enable reuse of data. In service of these goals, data ought to be FAIR: Findable, Accessible, Interoperable, and Reusable. The FaceBase data repository and educational resources exemplify the FAIR principles and support a broad user community including researchers in craniofacial development, molecular genetics, and genomics. FaceBase demonstrates that a model in which researchers "self-curate" their data can be successful and scalable. We present the results of the first 2.5 y of FaceBase's operations as an open community and summarize the data sets published during this period. We then describe a research highlight from work on the identification of regulatory networks and noncoding RNAs involved in cleft lip with/without cleft palate that both used and in turn contributed new findings to publicly available FaceBase resources. Collectively, FaceBase serves as a dynamic and continuously evolving resource to facilitate data-intensive research, enhance data reproducibility, and perform deep phenotyping across multiple species in dental and craniofacial research.

Single-cell transcriptomic signatures and gene regulatory networks modulated by Wls in mammalian midline facial formation and clefts

Formation of highly unique and complex facial structures is controlled by genetic programs that are responsible for the precise coordination of three-dimensional tissue morphogenesis. However, the underlying mechanisms governing these processes remain poorly understood. We combined mouse genetic and genomic approaches to define the mechanisms underlying normal and defective midfacial morphogenesis. Conditional inactivation of the Wnt secretion protein Wls in Pax3-expressing lineage cells disrupted frontonasal primordial patterning, cell survival and directional outgrowth, resulting in altered facial structures, including midfacial hypoplasia and midline facial clefts. Single-cell RNA sequencing revealed unique transcriptomic atlases of mesenchymal subpopulations in the midfacial primordia, which are disrupted in the conditional Wls mutants. Differentially expressed genes and cis-regulatory sequence analyses uncovered that Wls modulates and integrates a core gene regulatory network, consisting of key midfacial regulatory transcription factors (including Msx1, Pax3 and Pax7) and their downstream targets (including Wnt, Shh, Tgfβ and retinoic acid signaling components), in a mesenchymal subpopulation of the medial nasal prominences that is responsible for midline facial formation and fusion. These results reveal fundamental mechanisms underlying mammalian midfacial morphogenesis and related defects at single-cell resolution.

Evaluation of the multiple tissue factors in bone of primary osteoplasty and rhinoplasty in patients affected by cleft lip palate

Clefts of the lip and/or palate (CLP) are visible disruptions of standard facial structure. The aim of our study was to determine a relative number and appearance of the tissue factors in bone of patients with CLP during first time plastic alveolar osteoplasty or rhinoplasty. Immunohistochemistry was performed with matrix metalloproteinase-8 (MMP-8), matrix metalloproteinase-9 (MMP-9), osteopontin (OPN), osteocalcin (OC), Runt-related transcription factor 2 (Runx2), beta-defensin-2 (βdef-2), beta-defensin-3 (βdef-3), interleukin-1 alpha (IL-1α), and interleukin-10 (IL-10). The bone formation was observed by Masson-trichrome (Masson) staining. For the quantification of structures, the semi-quantitative census method was used. Spearman rank order correlation coefficient and Mann-Whitney U test were used for the statistical analysis. A significantly higher number of OPN positive osteocytes was observed in the CLP group when compared to the control group (p=0.002). The number of OC positive osteocytes (p=0.000) and βdef-2 positive osteocytes (p=0.003) was significantly lower in the CLP group in comparison to the control group. Strong, positive correlations between IL-10 and OC (rs=0.608; p=0.002), IL-1α and MMP-9 (rs=0.666; p=0.000), OPN and MMP-8 (rs=0.620; p=0.002) were detected in the CLP group. A tendency for the increased appearance of MMP-8, MMP-9 positive osteocytes of the patients with CLP, suggests elevated tissue remodelling properties. Increased appearance of OPN positive osteocytes in bone of the patients with CLP shows increased bone homeostasis based on seriously decreased mineralization, which may be a possible compensatory reaction to decreased quality of postsurgical bone.

Prenatal alcohol exposure contributes to negative pregnancy outcomes by altering fetal vascular dynamics and the placental transcriptome

Background

Prenatal alcohol exposure (PAE) has been shown to alter fetal blood flow in utero and is also associated with placental insufficiency and intrauterine growth restriction (IUGR), suggesting an underlying connection between perturbed circulation and pregnancy outcomes.

Methods

Timed‐pregnant C57/BL6NHsd mice, bred in‐house, were exposed by gavage on gestational day 10 (GD10) to ethanol (3 g/kg) or purified water, as a control. Pulse‐wave Doppler ultrasound measurements for umbilical arteries and ascending aorta were obtained post‐gavage (GD12, GD14, GD18) on 2 fetuses/litter. RNA from the non‐decidual (labyrinthine and junctional zone) portion of placentas was isolated and processed for RNA‐seq and subsequent bioinformatic analyses, and the association between transcriptomic changes and fetal phenotypes assessed.

Results

Exposure to ethanol in pregnant mice on GD10 attenuates umbilical cord blood flow transiently during gestation, and is associated with indices of IUGR, specifically decreased fetal weight and morphometric indices of cranial growth. Moreover, RNA‐seq of the fetal portion of the placenta demonstrated that this single exposure has lasting transcriptomic changes, including upregulation of Tet3, which is associated with spontaneous abortion. Weighted gene co‐expression network analysis (WGCNA) identified erythrocyte differentiation and homeostasis as important pathways associated with improved umbilical cord blood flow as gestation progresses. WGCNA also identified sensory perception of chemical stimulus/odorant and receptor activity as important pathways associated with cranial growth.

Conclusion

Our data suggest that PAE perturbs the expression of placental genes relevant for placental hematopoiesis and environmental sensing, resulting in transient impairment of umbilical cord blood flow and, subsequently, IUGR.

Toward clinical and molecular dissection of frontonasal dysplasia with facial skin polyps: From Pai syndrome to differential diagnosis through a series of 27 patients

Unique or multiple congenital facial skin polyps are features of several rare syndromes, from the most well-known Pai syndrome (PS), to the less recognized oculoauriculofrontonasal syndrome (OAFNS), encephalocraniocutaneous lipomatosis (ECCL), or Sakoda complex (SC). We set up a research project aiming to identify the molecular bases of PS. We reviewed 27 individuals presenting with a syndromic frontonasal polyp and initially referred for PS. Based on strict clinical classification criteria, we could confirm only nine (33%) typical and two (7%) atypical PS individuals. The remaining ones were either OAFNS (11/27-41%) or presenting with an overlapping syndrome (5/27-19%). Because of the phenotypic overlap between these entities, OAFNS, ECCL, and SC can be either considered as differential diagnosis of PS or part of the same spectrum. Exome and/or genome sequencing from blood DNA in 12 patients and from affected tissue in one patient failed to identify any replication in candidate genes. Taken together, our data suggest that conventional approaches routinely utilized for the identification of molecular etiologies responsible for Mendelian disorders are inconclusive. Future studies on affected tissues and multiomics studies will thus be required in order to address either the contribution of mosaic or noncoding variation in these diseases.

Analysis of Gene-Environment Interactions Related to Developmental Disorders

Various genetic and environmental factors are associated with developmental disorders (DDs). It has been suggested that interaction between genetic and environmental factors (G × E) is involved in the etiology of DDs. There are two major approaches to analyze the interaction: genome-wide and candidate gene-based approaches. In this mini-review, we demonstrate how these approaches can be applied to reveal the G × E related to DDs focusing on zebrafish and mouse models. We also discuss novel approaches to analyze the G × E associated with DDs.

Osteogenic and angiogenic profiles of the palatal process of the maxilla and the palatal process of the palatine bone

Hard palate consists anteriorly of the palatal process of the maxilla (ppmx) and posteriorly of the palatal process of the palatine (ppp). Currently, palatal osteogenesis is receiving increasing attention. This is the first study to provide an overview of the osteogenesis process of the mouse hard palate. We found that the period in which avascular mesenchymal condensation becomes a vascularized bone structure corresponds to embryonic day (E) 14.5 to E16.5 in the hard palate. The ppmx and ppp differ remarkably in morphology and molecular respects during osteogenesis. Osteoclasts in the ppmx and ppp are heterogeneous. There was a multinucleated giant osteoclast on the bone surface at the lateral-nasal side of the ppmx, while osteoclasts in the ppp were more abundant and adjacent to blood vessels but were smaller and had fewer nuclei. In addition, bone remodeling in the hard palate was asymmetric and exclusively occurred on the nasal side of the hard palate at E18.5. During angiogenesis, CD31-positive endothelial cells were initially localized in the surrounding of palatal mesenchymal condensation and then invaded the condensation in a sprouting fashion. At the transcriptome level, we found 78 differentially expressed genes related to osteogenesis and angiogenesis between the ppmx and ppp. Fifty-five related genes were up/downregulated from E14.5 to E16.5. Here, we described the morphogenesis and the heterogeneity in the osteogenic and angiogenic genes profiles of the ppmx and ppp, which are significant for subsequent studies of normal and abnormal subjects.

Distinct and redundant roles for zebrafish her genes during mineralization and craniofacial patterning

The Notch pathway is a cell-cell communication system which is critical for many developmental processes, including craniofacial development. Notch receptor activation induces expression of several well-known canonical targets including those encoded by the hes and her genes in mammals and zebrafish, respectively. The function of these genes, individually and in combination, during craniofacial development is not well understood. Here, we used zebrafish genetics to investigate her9 and her6 gene function during craniofacial development. We found that her9 is required for osteoblasts to efficiently mineralize bone, while cartilage is largely unaffected. Strikingly, gene expression studies in her9 mutants indicate that although progenitor cells differentiate into osteoblasts at the appropriate time and place, they fail to efficiently lay down mineralized matrix. This mineralization role of her9 is likely independent of Notch activation. In contrast, her9 also functions redundantly with her6 downstream of Jagged1b-induced Notch activation during dorsoventral craniofacial patterning. These studies disentangle distinct and redundant her gene functions during craniofacial development, including an unexpected, Notch independent, requirement during bone mineralization.

The Role of Sonic Hedgehog in Human Holoprosencephaly and Short-Rib Polydactyly Syndromes

The Hedgehog (HH) signalling pathway is one of the major pathways controlling cell differentiation and proliferation during human development. This pathway is complex, with HH function influenced by inhibitors, promotors, interactions with other signalling pathways, and non-genetic and cellular factors. Many aspects of this pathway are not yet clarified. The main features of Sonic Hedgehog (SHH) signalling are discussed in relation to its function in human development. The possible role of SHH will be considered using examples of holoprosencephaly and short-rib polydactyly (SRP) syndromes. In these syndromes, there is wide variability in phenotype even with the same genetic mutation, so that other factors must influence the outcome. SHH mutations were the first identified genetic causes of holoprosencephaly, but many other genes and environmental factors can cause malformations in the holoprosencephaly spectrum. Many patients with SRP have genetic defects affecting primary cilia, structures found on most mammalian cells which are thought to be necessary for canonical HH signal transduction. Although SHH signalling is affected in both these genetic conditions, there is little overlap in phenotype. Possible explanations will be canvassed, using data from published human and animal studies. Implications for the understanding of SHH signalling in humans will be discussed.

Inhibition of Cyclooxygenase-2 Alters Craniofacial Patterning during Early Embryonic Development of Chick

A recent study from our lab revealed that the inhibition of cyclooxygenase-2 (COX-2) exclusively reduces the level of PGE2 (Prostaglandin E2) among prostanoids and hampers the normal development of several structures, strikingly the cranial vault, in chick embryos. In order to unearth the mechanism behind the deviant development of cranial features, the expression pattern of various factors that are known to influence cranial neural crest cell (CNCC) migration was checked in chick embryos after inhibiting COX-2 activity using etoricoxib. The compromised level of cell adhesion molecules and their upstream regulators, namely CDH1 (E-cadherin), CDH2 (N-cadherin), MSX1 (Msh homeobox 1), and TGF-β (Transforming growth factor beta), observed in the etoricoxib-treated embryos indicate that COX-2, through its downstream effector PGE2, regulates the expression of these factors perhaps to aid the migration of CNCCs. The histological features and levels of FoxD3 (Forkhead box D3), as well as PCNA (Proliferating cell nuclear antigen), further consolidate the role of COX-2 in the migration and survival of CNCCs in developing embryos. The results of the current study indicate that COX-2 plays a pivotal role in orchestrating craniofacial structures perhaps by modulating CNCC proliferation and migration during the embryonic development of chicks.

Defining the critical period of hedgehog pathway inhibitor-induced cranial base dysplasia in mice

BACKGROUND

The hedgehog signaling pathway is critical for developmental patterning of the limb, craniofacial and axial skeleton. Disruption of this pathway in mice leads to a series of structural malformations, but the exact role and critical period of the Hh pathway in the early development of the cranial base have been rarely described.

RESULTS

Embryos exposed to vismodegib from E7.5, E9.5, and E10.5 had a higher percentage of cranial base fenestra. The peak incidence of hypoplasia in sphenoid winglets and severe craniosynostosis in cranial base synchondroses was observed when vismodegib was administered between E9.5 and E10.5. Cranial base craniosynostosis results from accelerating terminal differentiation of chondrocytes and premature osteogenesis.

CONCLUSIONS

We define the critical periods for the induction of cranial base deformity by vismodegib administration at a meticulous temporal resolution. Our findings suggest that the Hh pathway may play a vital role in the early development of the cranial base. This research also establishes a novel and easy-to-establish mouse model of synostosis in the cranial base using a commercially available pathway-selective inhibitor.

Adjunctive orthodontic therapy for prosthetic rehabilitation in a growing child with Axenfeld-Rieger syndrome: A case report

AIM

The congenital oligodontia impeding the development of the alveolar process resulting in disproportionate jaw growth has been previously reported. This case report describes the interdisciplinary management of an 11-year-old girl with Axenfeld-Rieger syndrome exhibiting oligodontia and maxillary hypoplasia.

METHODS AND RESULTS

An adjunctive orthodontic therapy was performed by taking advantage of bone-anchored maxillary protraction (BAMP) therapy using miniplates and 24-hour traction by intermaxillary Class III elastics. After 6 months of active treatment, the maxilla advanced by approximately 5 mm and upper lip relation improved by 3 mm without any significant changes in vertical relations. All the changes were maintained at 15th-month follow-up. A temporary removable prosthesis was given for immediate esthetics, and the definitive management is discussed.

CONCLUSION

The BAMP therapy can be a befitting alternative in cases exhibiting complex presentation involving skeletal and dentoalveolar components. An appreciable profile improvement without any dentoalveolar side effects can be achieved with BAMP therapy.

Phenytoin Inhibits Cell Proliferation through microRNA-196a-5p in Mouse Lip Mesenchymal Cells

Cleft lip (CL) is one of the most common birth defects. It is caused by either genetic mutations or environmental factors. Recent studies suggest that environmental factors influence the expression of noncoding RNAs [e.g., microRNA (miRNA)], which can regulate the expression of genes crucial for cellular functions. In this study, we examined which miRNAs are associated with CL. Among 10 candidate miRNAs (miR-98-3p, miR-101a-3p, miR-101b-3p, miR-141-3p, miR-144-3p, miR-181a-5p, miR-196a-5p, miR-196b-5p, miR-200a-3p, and miR-710) identified through our bioinformatic analysis of CL-associated genes, overexpression of miR-181a-5p, miR-196a-5p, miR-196b-5p, and miR-710 inhibited cell proliferation through suppression of genes associated with CL in cultured mouse embryonic lip mesenchymal cells (MELM cells) and O9-1 cells, a mouse cranial neural crest cell line. In addition, we found that phenytoin, an inducer of CL, decreased cell proliferation through miR-196a-5p induction. Notably, treatment with a specific inhibitor for miR-196a-5p restored cell proliferation through normalization of expression of CL-associated genes in the cells treated with phenytoin. Taken together, our results suggest that phenytoin induces CL through miR-196a-5p induction, which suppresses the expression of CL-associated genes.

Neural crest-specific deletion of Bmp7 leads to midfacial hypoplasia, nasal airway obstruction, and disordered breathing modelling Obstructive Sleep Apnea

Pediatric obstructive sleep apnea (OSA), a relatively common sleep-related breathing disorder (SRBD) affecting approximately 1-5% of children, is often caused by anatomical obstruction and/or collapse of the nasal and/or pharyngeal airways. The resulting sleep disruption and intermittent hypoxia lead to various systemic morbidities. Predicting the development of OSA from craniofacial features alone is currently not possible and a controversy remains if upper airway obstruction facilitates reduced midfacial growth or vice-versa. Currently, there is no rodent model that recapitulates both the development of craniofacial abnormalities and upper airway obstruction to address these questions. Here, we describe that mice with a neural crest-specific deletion of Bmp7 (Bmp7ncko) present with shorter, more acute angled cranial base, midfacial hypoplasia, nasal septum deviation, turbinate swelling and branching defects, and nasal airway obstruction. Interestingly, several of these craniofacial features develop after birth during periods of rapid midfacial growth and precede the development of an upper airway obstruction. We identified that in this rodent model, no single feature appeared to predict upper airway obstruction, but the sum of those features resulted in a reduced breathing frequency, apneas and overall reduced oxygen consumption. Metabolomics analysis of serum from peripheral blood identified increased levels of hydroxyproline, a metabolite upregulated under hypoxic conditions. As this model recapitulates many features observed in OSA, it offers unique opportunities for studying how upper airway obstruction affects breathing physiology and leads to systemic morbidities.

A developmental stage-specific network approach for studying dynamic co-regulation of transcription factors and microRNAs during craniofacial development

Craniofacial development is regulated through dynamic and complex mechanisms that involve various signaling cascades and gene regulations. Disruption of such regulations can result in craniofacial birth defects. Here, we propose the first developmental stage-specific network approach by integrating two crucial regulators, transcription factors (TFs) and microRNAs (miRNAs), to study their co-regulation during craniofacial development. Specifically, we used TFs, miRNAs and non-TF genes to form feed-forward loops (FFLs) using genomic data covering mouse embryonic days E10.5 to E14.5. We identified key novel regulators (TFs Foxm1, Hif1a, Zbtb16, Myog, Myod1 and Tcf7, and miRNAs miR-340-5p and miR-129-5p) and target genes (Col1a1, Sgms2 and Slc8a3) expression of which changed in a developmental stage-dependent manner. We found that the Wnt-FoxO-Hippo pathway (from E10.5 to E11.5), tissue remodeling (from E12.5 to E13.5) and miR-129-5p-mediated Col1a1 regulation (from E10.5 to E14.5) might play crucial roles in craniofacial development. Enrichment analyses further suggested their functions. Our experiments validated the regulatory roles of miR-340-5p and Foxm1 in the Wnt-FoxO-Hippo subnetwork, as well as the role of miR-129-5p in the miR-129-5p-Col1a1 subnetwork. Thus, our study helps understand the comprehensive regulatory mechanisms for craniofacial development.

Non-syndromic Cleft Palate: An Overview on Human Genetic and Environmental Risk Factors

The epithelial and mesenchymal cells involved in early embryonic facial development are guided by complex regulatory mechanisms. Any factor perturbing the growth, approach and fusion of the frontonasal and maxillary processes could result in orofacial clefts that represent the most common craniofacial malformations in humans. The rarest and, probably for this reason, the least studied form of cleft involves only the secondary palate, which is posterior to the incisive foramen. The etiology of cleft palate only is multifactorial and involves both genetic and environmental risk factors. The intention of this review is to give the reader an overview of the efforts made by researchers to shed light on the underlying causes of this birth defect. Most of the scientific papers suggesting potential environmental and genetic causes of non-syndromic cleft palate are summarized in this review, including genome-wide association and gene–environment interaction studies.

Critical microRNAs and regulatory motifs in cleft palate identified by a conserved miRNA-TF-gene network approach in humans and mice

Cleft palate (CP) is the second most common congenital birth defect. The etiology of CP is complicated, with involvement of various genetic and environmental factors. To investigate the gene regulatory mechanisms, we designed a powerful regulatory analytical approach to identify the conserved regulatory networks in humans and mice, from which we identified critical microRNAs (miRNAs), target genes and regulatory motifs (miRNA-TF-gene) related to CP. Using our manually curated genes and miRNAs with evidence in CP in humans and mice, we constructed miRNA and transcription factor (TF) co-regulation networks for both humans and mice. A consensus regulatory loop (miR17/miR20a-FOXE1-PDGFRA) and eight miRNAs (miR-140, miR-17, miR-18a, miR-19a, miR-19b, miR-20a, miR-451a and miR-92a) were discovered in both humans and mice. The role of miR-140, which had the strongest association with CP, was investigated in both human and mouse palate cells. The overexpression of miR-140-5p, but not miR-140-3p, significantly inhibited cell proliferation. We further examined whether miR-140 overexpression could suppress the expression of its predicted target genes (BMP2, FGF9, PAX9 and PDGFRA). Our results indicated that miR-140-5p overexpression suppressed the expression of BMP2 and FGF9 in cultured human palate cells and Fgf9 and Pdgfra in cultured mouse palate cells. In summary, our conserved miRNA-TF-gene regulatory network approach is effective in detecting consensus miRNAs, motifs, and regulatory mechanisms in human and mouse CP.

Msx1 deficiency interacts with hypoxia and induces a morphogenetic regulation during mouse lip development

Nonsyndromic clefts of the lip and palate are common birth defects resulting from gene-gene and gene-environment interactions. Mutations in human MSX1 have been linked to orofacial clefting and we show here that Msx1 deficiency causes a growth defect of the medial nasal process (Mnp) in mouse embryos. Although this defect alone does not disrupt lip formation, Msx1-deficient embryos develop a cleft lip when the mother is transiently exposed to reduced oxygen levels or to phenytoin, a drug known to cause embryonic hypoxia. In the absence of interacting environmental factors, the Mnp growth defect caused by Msx1 deficiency is modified by a Pax9-dependent 'morphogenetic regulation', which modulates Mnp shape, rescues lip formation and involves a localized abrogation of Bmp4-mediated repression of Pax9 Analyses of GWAS data revealed a genome-wide significant association of a Gene Ontology morphogenesis term (including assigned roles for MSX1, MSX2, PAX9, BMP4 and GREM1) specifically for nonsyndromic cleft lip with cleft palate. Our data indicate that MSX1 mutations could increase the risk for cleft lip formation by interacting with an impaired morphogenetic regulation that adjusts Mnp shape, or through interactions that inhibit Mnp growth.

Development and gene expression of C57BL/6 mouse embryo palate shelves in rotary organ culture

The aim of the present study was to improve methods for the suspension culture of mouse palatal shelves by comparing the expression of platelet-derived growth factor receptor (PDGFR)-α in palatal shelves in vivo, to that in vitro. The palatal shelves of C57BL/6 mouse embryos were obtained on gestation days (GDs) 13.5, 14.5, 15.0 and 15.5 for in vivo experiments. The palatal shelves were removed and observed under a stereomicroscope to investigate palatal development. For in vitro experiments, the palatal shelves were dissected under a stereomicroscope on GD 13.5 and then subjected to rotary culture for 0, 24, 36 or 48 h. The expression of PDGFR-α at different time points was detected by immunohistochemical staining and western blot analysis. Both methods of analysis displayed PDGFR-α expression in mesenchymal and epithelial cells at GD 13.5, 14.5, 15.0 and 15.5, in vivo and in vitro. The level of PDGFR-α expression peaked on GD 14.5. The expression of PDGFR-α in palatal shelves in in vitro rotary culture was consistent with that in vivo. Therefore, the novel technique of palatal rotary organ culture presented in the current study could provide a good model for studying the mechanism of pathological palatal fusion in vitro. Additionally, the present study further confirmed that PDGFR-α gene expression was associated with the development of palatal shelves.

MicroRNA-124-3p suppresses mouse lip mesenchymal cell proliferation through the regulation of genes associated with cleft lip in the mouse

Background

Cleft lip (CL), one of the most common congenital birth defects, shows considerable geographic and ethnic variation, with contribution of both genetic and environmental factors. Mouse genetic studies have identified several CL-associated genes. However, it remains elusive how these CL-associated genes are regulated and involved in CL. Environmental factors may regulate these genes at the post-transcriptional level through the regulation of non-coding microRNAs (miRNAs). In this study, we sought to identify miRNAs associated with CL in mice.

Results

Through a systematic literature review and a Mouse Genome Informatics (MGI) database search, we identified 55 genes that were associated with CL in mice. Subsequent bioinformatic analysis of these genes predicted that a total of 33 miRNAs target multiple CL-associated genes, with 20 CL-associated genes being potentially regulated by multiple miRNAs. To experimentally validate miRNA function in cell proliferation, we conducted cell proliferation/viability assays for the selected five candidate miRNAs (miR-124-3p, let-7a-5p, let-7b-5p, let-7c-5p, and let-7d-5p). Overexpression of miR-124-3p, but not of the others, inhibited cell proliferation through suppression of CL-associated genes in cultured mouse embryonic lip mesenchymal cells (MELM cells) isolated from the developing mouse lip region. By contrast, miR-124-3p knockdown had no effect on MELM cell proliferation. This miRNA-gene regulatory mechanism was mostly conserved in O9–1 cells, an established cranial neural crest cell line. Expression of miR-124-3p was low in the maxillary processes at E10.5, when lip mesenchymal cells proliferate, whereas it was greatly increased at later developmental stages, suggesting that miR-124-3p expression is suppressed during the proliferation phase in normal palate development.

Conclusions

Our findings indicate that upregulated miR-124-3p inhibits cell proliferation in cultured lip cells through suppression of CL-associated genes. These results will have a significant impact, not only on our knowledge about lip morphogenesis, but also on the development of clinical approaches for the diagnosis and prevention of CL.

Nubp2 is required for cranial neural crest survival in the mouse

The N-ethyl-N-nitrosourea (ENU) ←forward genetic screen is a useful tool for the unbiased discovery of novel mechanisms regulating developmental processes. We recovered the dorothy mutation in such a screen designed to recover recessive mutations affecting craniofacial development in the mouse. Dorothy embryos die prenatally and exhibit many striking phenotypes commonly associated with ciliopathies, including a severe midfacial clefting phenotype. We used exome sequencing to discover a missense mutation in nucleotide binding protein 2 (Nubp2) to be causative. This finding was confirmed by a complementation assay with the dorothy allele and an independent Nubp2 null allele (Nubp2^null). We demonstrated that Nubp2 is indispensable for embryogenesis. NUBP2 is implicated in both the cytosolic iron/sulfur cluster assembly pathway and negative regulation of ciliogenesis. Conditional ablation of Nubp2 in the neural crest lineage with Wnt1-cre recapitulates the dorothy craniofacial phenotype. Using this model, we found that the proportion of ciliated cells in the craniofacial mesenchyme was unchanged, and that markers of the SHH, FGF, and BMP signaling pathways are unaltered. Finally, we show evidence that the phenotype results from a marked increase in apoptosis within the craniofacial mesenchyme.

A combination insecticide at sub-lethal dose debilitated the expression pattern of crucial signalling molecules that facilitate craniofacial patterning in domestic chick Gallus domesticus

Pesticides despite being agents that protect the plants and humans from noxious pests, are infamous for their potential to cause detrimental health issues in nontargeted species. In order to ascertain the latter, a set of experiments were conducted by exposing early chick embryos to a widely used combination insecticide (Ci, 50% chlorpyrifos and 5% cypermethrin). The results revealed a myriad of congenital defects pertaining to craniofacial development such as anophthalmia, microphthalmia, exencephaly as well as deformed beak and cranial structures. These teratological manifestations could be attributed to the Ci induced alteration in the titre of major regulators of neurulation and ossification. Therefore, the mRNA and/or the protein level expression pattern of genes which are reported to be involved in the craniofacial development were studied at selected time points of embryonic development. The analysis of the result showed that there have been significant alternations in the expression patterns of the signalling molecules such as SHH, WNTs, CDH1, CDH2, L1CAM, PAX6, HOX, PCNA, GLI3, BMP7, FGF8, GLIs, SOX9, RUNX2, DLX5, COL10A1, CASPASE3 etc. on embryonic days 2, 4 and/or 10. Concurrently, on day 10, whole-mount skeletal staining and biochemical estimation of hydroxyproline were carried out in the cranial tissues of the embryos. The overall result of the current study indicates that exposure to Ci during early development impede the crucial regulatory signals that orchestrate the morphogenesis of cranial neural crest cells thereby hindering the normal progression of neural tube and endochondral ossification which collectively lead to craniofacial dysmorphism in domestic chicks.

The molecular anatomy of mammalian upper lip and primary palate fusion at single cell resolution

The mammalian lip and primary palate form when coordinated growth and morphogenesis bring the nasal and maxillary processes into contact, and the epithelia co-mingle, remodel and clear from the fusion site to allow mesenchyme continuity. Although several genes required for fusion have been identified, an integrated molecular and cellular description of the overall process is lacking. Here, we employ single cell RNA sequencing of the developing mouse face to identify ectodermal, mesenchymal and endothelial populations associated with patterning and fusion of the facial prominences. This analysis indicates that key cell populations at the fusion site exist within the periderm, basal epithelial cells and adjacent mesenchyme. We describe the expression profiles that make each population unique, and the signals that potentially integrate their behaviour. Overall, these data provide a comprehensive high-resolution description of the various cell populations participating in fusion of the lip and primary palate, as well as formation of the nasolacrimal groove, and they furnish a powerful resource for those investigating the molecular genetics of facial development and facial clefting that can be mined for crucial mechanistic information concerning this prevalent human birth defect.

MicroRNA-655-3p and microRNA-497-5p inhibit cell proliferation in cultured human lip cells through the regulation of genes related to human cleft lip

Background

The etiology of cleft lip with or without cleft palate (CL/P), a common congenital birth defect, is complex and involves the contribution of genetic and environmental factors. Although many candidate genes have been identified, the regulation and interaction of these genes in CL/P remain unclear. In addition, the contribution of microRNAs (miRNAs), non-coding RNAs that regulate the expression of multiple genes, to the etiology of CL/P is largely unknown.

Methods

To identify the signatures of causative biological pathways for human CL/P, we conducted a systematic literature review for human CL/P candidate genes and subsequent bioinformatics analyses. Functional enrichment analyses of the candidate CL/P genes were conducted using the pathway databases GO and KEGG. The miRNA-mediated post-transcriptional regulation of the CL/P candidate genes was analyzed with miRanda, PITA, and TargetScan, and miRTarbase. Genotype-phenotype association analysis was conducted using GWAS. The functional significance of the candidate miRNAs was evaluated experimentally in cell proliferation and target gene regulation assays in human lip fibroblasts.

Results

Through an extensive search of the main biomedical databases, we mined 177 genes with mutations or association/linkage reported in individuals with CL/P, and considered them as candidate genes for human CL/P. The genotype-phenotype association study revealed that mutations in 12 genes (ABCA4, ADAM3A, FOXE1, IRF6, MSX2, MTHFR, NTN1, PAX7, TP63, TPM1, VAX1, and WNT9B) were significantly associated with CL/P. In addition, our bioinformatics analysis predicted 16 microRNAs (miRNAs) to be post-transcriptional regulators of CL/P genes. To validate the bioinformatics results, the top six candidate miRNAs (miR-124-3p, miR-369-3p, miR-374a-5p, miR-374b-5p, miR-497-5p, and miR-655-3p) were evaluated by cell proliferation/survival assays and miRNA-gene regulation assays in cultured human lip fibroblasts. We found that miR-497-5p and miR-655-3p significantly suppressed cell proliferation in these cells. Furthermore, the expression of the predicted miRNA-target genes was significantly downregulated by either miR-497-5p or miR-655-3p mimic.

Conclusion

Expression of miR-497-5p and miR-655-3p suppresses cell proliferation through the regulation of human CL/P-candidate genes. This study provides insights into the role of miRNAs in the etiology of CL/P and suggests possible strategies for the diagnosis of CL/P.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0535-2) contains supplementary material, which is available to authorized users.

Association of IFT88 gene variants with nonsyndromic cleft lip with or without cleft palate

BACKGROUND

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect with multifactorial etiology. Genetic studies have identified numerous gene variants in association with NSCLP. IFT88 (intraflagellar transport 88) has been suggested to play a major role in craniofacial development, as Ift88 mutant mice exhibit cleft palate and mutations in IFT88 were identified in individuals with NSCLP.

OBJECTIVE

To investigate the association of IFT88 single nucleotide gene variants (SNVs) with NSCLP in a large family data set consisting of non-Hispanic white (NHW) and Hispanic families.

METHODS

Nine SNVs in/nearby IFT88 were genotyped in 482 NHW families and 301 Hispanic NSCLP families. Genotyping was performed using TaqMan® chemistry. Single- and pairwise-SNV association analyses were performed for all families stratified by ethnicity and family history of NSCLP using the family-based association test (FBAT), and association in the presence of linkage (APL). Bonferroni correction was used to adjust for multiple testing and p values ≤.0055 were considered statistically significant.

RESULTS

Significant association was found between IFT88 rs9509311 and rs2497490 and NSCLP in NHW all families (p = .004 and .005, respectively), while nominal associations were found for rs7998361 and rs9509307 (p < .05). Pairwise association analyses also showed nominal associations between NSCLP in both NHW and Hispanic data sets (p < .05). No association was found between individual variants in IFT88 and NSCLP in Hispanics.

CONCLUSIONS

Our results suggest that variation in IFT88 may contribute to NSCLP risk, particularly in multiplex families from a non-Hispanic white population.

Wnt signaling in orofacial clefts: crosstalk, pathogenesis and models

Diverse signaling cues and attendant proteins work together during organogenesis, including craniofacial development. Lip and palate formation starts as early as the fourth week of gestation in humans or embryonic day 9.5 in mice. Disruptions in these early events may cause serious consequences, such as orofacial clefts, mainly cleft lip and/or cleft palate. Morphogenetic Wnt signaling, along with other signaling pathways and transcription regulation mechanisms, plays crucial roles during embryonic development, yet the signaling mechanisms and interactions in lip and palate formation and fusion remain poorly understood. Various Wnt signaling and related genes have been associated with orofacial clefts. This Review discusses the role of Wnt signaling and its crosstalk with cell adhesion molecules, transcription factors, epigenetic regulators and other morphogenetic signaling pathways, including the Bmp, Fgf, Tgfβ, Shh and retinoic acid pathways, in orofacial clefts in humans and animal models, which may provide a better understanding of these disorders and could be applied towards prevention and treatments.

Transcriptome analysis of Xenopus orofacial tissues deficient in retinoic acid receptor function

BACKGROUND

Development of the face and mouth is orchestrated by a large number of transcription factors, signaling pathways and epigenetic regulators. While we know many of these regulators, our understanding of how they interact with each other and implement changes in gene expression during orofacial development is still in its infancy. Therefore, this study focuses on uncovering potential cooperation between transcriptional regulators and one important signaling pathway, retinoic acid, during development of the midface.

RESULTS

Transcriptome analyses was performed on facial tissues deficient for retinoic acid receptor function at two time points in development; early (35 hpf) just after the neural crest migrates and facial tissues are specified and later (60 hpf) when the mouth has formed and facial structures begin to differentiate. Functional and network analyses revealed that retinoic acid signaling could cooperate with novel epigenetic factors and calcium-NFAT signaling during early orofacial development. At the later stage, retinoic acid may work with WNT and BMP and regulate homeobox containing transcription factors. Finally, there is an overlap in genes dysregulated in Xenopus embryos with median clefts with human genes associated with similar orofacial defects.

CONCLUSIONS

This study uncovers novel signaling pathways required for orofacial development as well as pathways that could interact with retinoic acid signaling during the formation of the face. We show that frog faces are an important tool for studying orofacial development and birth defects.

Choanal atresia and stenosis: Development and diseases of the nasal cavity

Proper craniofacial development in vertebrates depends on growth and fusion of the facial processes during embryogenesis. Failure of any step in this process could lead to craniofacial anomalies such as facial clefting, which has been well studied with regard to its molecular etiology and cellular pathogenesis. Nasal cavity invagination is also a critical event in proper craniofacial development, and is required for the formation of a functional nasal cavity and airway. The nasal cavity must connect the nasopharynx with the primitive choanae to complete an airway from the nostril to the nasopharynx. In contrast to orofacial clefts, defects in nasal cavity and airway formation, such as choanal atresia (CA), in which the connection between the nasal airway and nasopharynx is physically blocked, have largely been understudied. This is also true for a narrowed connection between the nasal cavity and the nasopharynx, which is known as choanal stenosis (CS). CA occurs in approximately 1 in 5,000 live births, and can present in isolation but typically arises as part of a syndrome. Despite the fact that CA and CS usually require immediate intervention, and substantially affect the quality of life of affected individuals, the etiology and pathogenesis of CA and CS have remained elusive. In this review I focus on the process of nasal cavity development with respect to forming a functional airway and discuss the cellular behavior and molecular networks governing this process. Additionally, the etiology of human CA is discussed using examples of disorders which involve CA or CS. This article is categorized under: Signaling Pathways > Cell Fate Signaling Comparative Development and Evolution > Model Systems Birth Defects > Craniofacial and Nervous System Anomalies.

Mouse models for the study of cranial base growth and anomalies

The cranial base is a central and integral component of the cranioskeleton, yet little is known about its growth. Despite the dissimilarities between human and murine cranioskeletal form, mouse models are proving instrumental in studying craniofacial growth. The objectives of this review are to summarize recent findings from numerous mouse models that display growth defects in one or more cranial base synchondroses, with accompanying changes in chondrocyte cellular zones. Many of these models also display altered growth of the cranial vault and/or the facial region. FGFR, PTHrP, Ihh, BMP and Wnt/β-catenin, as well as components of primary cilia, are the major genes and signalling pathways identified in cranial base synchondroses. Together, these models are helping to uncover specific genetic influences and signalling pathways operational at the cranial base synchondroses. Many of these genes are in common with those of importance in the cranial vault and the facial skeleton, emphasizing the molecular integration of growth between the cranial base and other cranial regions. Selected models are also being utilized in testing therapeutic agents to correct defective craniofacial and cranial base growth.

The Association between Maxillary Sinus Dimensions and Midface Parameters during Human Postnatal Growth

Objective

The aim of the study based on CT images was to assess the age-related changes in maxillary sinus diameters in relation to diameters of the facial skeleton.

Materials and Methods

The retrospective analysis of CT images of the head of 170 patients aged 0-18 years (85 females and 85 males) was performed. Specific orientation points (zy, zm, pr, ns, n, and P) were identified in every patient and the following distances were measured: zy-zy, maximum facial width; zm-zm, midfacial width; n-pr, upper facial height; ns-pr, alveolar facial height; and ns-P, distance not indicated in craniometry.

Results

The maxillary sinuses of every patient were bilaterally measured in three planes. Three diameters were obtained: maximum transverse (horizontal) diameter called MSW, maximum vertical diameter called MSH, and maximum anteroposterior diameter (length) called MSL. In females, the correlation of MSW, MSH, and MSL and zy-zy, as well as n-pr distances, is very strong. Moreover, the significant correlation was found between all measurements of maxillary sinus and ns-pr as well as ns-P distances in females. The correlation between MSL and all measurements of midface as well as MSH and MSW and all measurements except ns-P is stronger in females than in males. In males, all measurements of maxillary sinus correlate with ns-P distance very strongly.

Conclusions

The statistical analysis (correlation and determination coefficient) showed that all measurements of maxillary sinuses correlate with midface dimensions.

Neural crest and the patterning of vertebrate craniofacial muscles

Patterning of craniofacial muscles overtly begins with the activation of lineage-specific markers at precise, evolutionarily conserved locations within prechordal, lateral, and both unsegmented and somitic paraxial mesoderm populations. Although these initial programming events occur without influence of neural crest cells, the subsequent movements and differentiation stages of most head muscles are neural crest-dependent. Incorporating both descriptive and experimental studies, this review examines each stage of myogenesis up through the formation of attachments to their skeletal partners. We present the similarities among developing muscle groups, including comparisons with trunk myogenesis, but emphasize the morphogenetic processes that are unique to each group and sometimes subsets of muscles within a group. These groups include branchial (pharyngeal) arches, which encompass both those with clear homologues in all vertebrate classes and those unique to one, for example, mammalian facial muscles, and also extraocular, laryngeal, tongue, and neck muscles. The presence of several distinct processes underlying neural crest:myoblast/myocyte interactions and behaviors is not surprising, given the wide range of both quantitative and qualitative variations in craniofacial muscle organization achieved during vertebrate evolution.

Role of Matrix Gla protein in midface development: Recent advances

Craniofacial development is a delicate process that involves complex interactions among cells of multiple developmental origins, their migration, proliferation, and differentiation. Tissue morphogenesis of the craniofacial skeleton depends on genetic and environmental factors, and on specific signaling pathways, which are still not well understood. Developmental defects of the midface caused by the absence, delays, or premature fusion of nasal and maxillary prominences vary in severity; leading to clefts, hypoplasias, and midline expansion. In the current review, we focus on the importance of the chondrocranium in craniofacial growth and how its impaired development leads to midface hypoplasia. More importantly, we reported how Matrix Gla protein (MGP), a potent inhibitor of extracellular matrix mineralization, facilitates midface development by preventing ectopic calcification of the nasal septum. In fact, MGP may act as a common link in multiple developmental pathologies all showing midface hypoplasia caused by abnormal cartilage calcification. This brief review discusses the gap in knowledge in the field, raises pertinent questions, which remain unanswered, and sheds light on the future research directions.

Genes and microRNAs associated with mouse cleft palate: A systematic review and bioinformatics analysis

Cleft palate (CP) is the most prevalent craniofacial deformity, with ethnic and geographic variation in prevalence in humans. Mice have been used as an animal model to study the cause(s) of CP by several approaches, including genetic and chemical-induced approaches. Mouse genetic approaches revealed that significant amounts of genes are involved in the CP pathology. The aim of this study was to identify common features of CP-associated genes and to explore the roles of microRNAs (miRNAs) as important post-transcriptional regulators that may be involved in the regulation of CP genes. To generate an accurate list of genes associated with CP, we first conducted systematic literature searches through main databases such as Medline, Embase, and PubMed, as well as other sources such as Scopus and Mouse Genome Informatics. We found that 195 mouse strains with single-gene mutations and 140 mouse strains with compound-gene mutations were reported to have CP. The CP genes were categorized by functions and pathways using the Kyoto Encyclopedia of Genes and Genomes and Gene Ontology annotations, highlighting the contribution of cellular metabolism to CP. A total of 18 miRNAs were involved in the regulation of multiple CP genes. Human genotype-phenotype analysis revealed that variants in five human homologous CP genes (IRF6, FOXE1, VAX1, WNT9B, and GAD1) significantly contributed to the human CP phenotype. Thus, our results suggest that cellular metabolism and miRNAs play an important role in the regulation of genetic pathways and networks crucial for palatal formation.

Mllt10 knockout mouse model reveals critical role of Af10-dependent H3K79 methylation in midfacial development

Epigenetic regulation is required to ensure the precise spatial and temporal pattern of gene expression that is necessary for embryonic development. Although the roles of some epigenetic modifications in embryonic development have been investigated in depth, the role of methylation at lysine 79 (H3K79me) is poorly understood. Dot1L, a unique methyltransferase for H3K79, forms complexes with distinct sets of co-factors. To further understand the role of H3K79me in embryogenesis, we generated a mouse knockout of Mllt10, the gene encoding Af10, one Dot1L complex co-factor. We find homozygous Mllt10 knockout mutants (Mllt10-KO) exhibit midline facial cleft. The midfacial defects of Mllt10-KO embryos correspond to hyperterolism and are associated with reduced proliferation of mesenchyme in developing nasal processes and adjacent tissue. We demonstrate that H3K79me level is significantly decreased in nasal processes of Mllt10-KO embryos. Importantly, we find that expression of AP2α, a gene critical for midfacial development, is directly regulated by Af10-dependent H3K79me, and expression AP2α is reduced specifically in nasal processes of Mllt10-KO embryos. Suppression of H3K79me completely mimicked the Mllt10-KO phenotype. Together these data are the first to demonstrate that Af10-dependent H3K79me is essential for development of nasal processes and adjacent tissues, and consequent midfacial formation.

Systems biology of facial development: contributions of ectoderm and mesenchyme

The rapid increase in gene-centric biological knowledge coupled with analytic approaches for genomewide data integration provides an opportunity to develop systems-level understanding of facial development. Experimental analyses have demonstrated the importance of signaling between the surface ectoderm and the underlying mesenchyme are coordinating facial patterning. However, current transcriptome data from the developing vertebrate face is dominated by the mesenchymal component, and the contributions of the ectoderm are not easily identified. We have generated transcriptome datasets from critical periods of mouse face formation that enable gene expression to be analyzed with respect to time, prominence, and tissue layer. Notably, by separating the ectoderm and mesenchyme we considerably improved the sensitivity compared to data obtained from whole prominences, with more genes detected over a wider dynamic range. From these data we generated a detailed description of ectoderm-specific developmental programs, including pan-ectodermal programs, prominence- specific programs and their temporal dynamics. The genes and pathways represented in these programs provide mechanistic insights into several aspects of ectodermal development. We also used these data to identify co-expression modules specific to facial development. We then used 14 co-expression modules enriched for genes involved in orofacial clefts to make specific mechanistic predictions about genes involved in tongue specification, in nasal process patterning and in jaw development. Our multidimensional gene expression dataset is a unique resource for systems analysis of the developing face; our co-expression modules are a resource for predicting functions of poorly annotated genes, or for predicting roles for genes that have yet to be studied in the context of facial development; and our analytic approaches provide a paradigm for analysis of other complex developmental programs.