Insights Into the Complexity of Craniofacial Development From a Cellular Perspective

Depolarization induces calcium-dependent BMP4 release from mouse embryonic palate mesenchyme

Ion channels are essential for proper morphogenesis of the craniofacial skeleton. However, the molecular mechanisms underlying this phenomenon are unknown. Loss of the Kcnj2 potassium channel disrupts Bone Morphogenetic Protein (BMP) signaling within the developing palate. BMP signaling is essential for the correct development of several skeletal structures, including the palate, though little is known about the mechanisms that govern BMP secretion. We introduce a tool to image the release of bone morphogenetic protein 4 (BMP4) from mammalian cells. Using this tool, we show that depolarization induces BMP4 release from mouse embryonic palate mesenchyme cells in a calcium-dependent manner. We show native transient changes in intracellular calcium occur in cranial neural crest cells, the cells from which embryonic palate mesenchyme derives. Waves of transient changes in intracellular calcium suggest that these cells are electrically coupled and may temporally coordinate BMP release. These transient changes in intracellular calcium persist in palate mesenchyme cells from embryonic day (E) 9.5 to 13.5 mice. Disruption of Kcnj2 significantly decreases the amplitude of calcium transients and the ability of cells to secrete BMP. Together, these data suggest that temporal control of developmental cues is regulated by ion channels, depolarization, and changes in intracellular calcium for mammalian craniofacial morphogenesis.

SUMMARY

We show that embryonic palate mesenchyme cells undergo transient changes in intracellular calcium. Depolarization of these cells induces BMP4 release suggesting that ion channels are a node in BMP4 signaling.

Dynamic enhancer landscapes in human craniofacial development

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Exposure to the aryl hydrocarbon receptor agonist dioxin disrupts formation of the muscle, nerves, and vasculature in the developing jaw

Human exposure to environmental pollutants can disrupt embryonic development and impact juvenile and adult health outcomes by adversely affecting cell and organ function. Notwithstanding, environmental contamination continues to increase due to industrial development, insufficient regulations, and the mobilization of pollutants as a result of extreme weather events. Dioxins are a class of structurally related persistent organic pollutants that are highly toxic, carcinogenic, and teratogenic. 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most potent dioxin compound and has been shown to induce toxic effects in developing organisms by activating the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor targeted by multiple persistent organic pollutants. Contaminant-induced AHR activation results in malformations of the craniofacial cartilages and neurocranium; however, the mechanisms mediating these phenotypes are not well understood. In this study, we utilized the optically transparent zebrafish model to elucidate novel cellular targets and potential transcriptional targets underlying TCDD-induced craniofacial malformations. To this end, we exposed zebrafish embryos at 4 h post fertilization to TCDD and employed a mixed-methods approach utilizing immunohistochemistry staining, transgenic reporter lines, fixed and in vivo confocal imaging, and timelapse microscopy to determine the targets mediating TCDD-induced craniofacial phenotypes. Our data indicate that embryonic TCDD exposure reduced jaw and pharyngeal arch Sox10+ chondrocytes and Tcf21+ pharyngeal mesoderm progenitors. Exposure to TCDD correspondingly led to a reduction in collagen type II deposition in Sox10+ domains. Embryonic TCDD exposure impaired development of tissues derived from or guided by Tcf21+ progenitors, namely: nerves, muscle, and vasculature. Specifically, TCDD exposure disrupted development of the hyoid and mandibular arch muscles, decreased neural innervation of the jaw, resulted in compression of cranial nerves V and VII, and led to jaw vasculature malformations. Collectively, these findings reveal novel structural targets and potential transcriptional targets of TCDD-induced toxicity, showcasing how contaminant exposures lead to congenital craniofacial malformations.

Mutation in the Bone Morphogenetic Protein signaling pathway sensitize zebrafish and humans to ethanol-induced jaw malformations

Fetal Alcohol Spectrum Disorders (FASD) describe ethanol-induced developmental defects including craniofacial malformations. While ethanol-sensitive genetic mutations contribute to facial malformations, the impacted cellular mechanisms remain unknown. Bmp signaling is a key regulator of epithelial morphogenesis driving facial development, providing a possible ethanol-sensitive mechanism. We found that zebrafish mutants for Bmp signaling components are ethanol-sensitive and affect anterior pharyngeal endoderm shape and gene expression, indicating ethanol-induced malformations of the anterior pharyngeal endoderm cause facial malformations. Integrating FASD patient data, we provide the first evidence that variants in the human Bmp receptor gene BMPR1B associate with ethanol-related differences in jaw volume. Our results show that ethanol exposure disrupts proper morphogenesis of, and tissue interactions between, facial epithelia that mirror overall viscerocranial shape changes and are predictive for Bmp-ethanol associations in human jaw development. Our data provide a mechanistic paradigm linking ethanol to disrupted epithelial cell behaviors that underlie facial defects in FASD.

RA-induced prominence-specific response resulted in distinctive regulation of Wnt and osteogenesis

Proper retinoic acid (RA) signaling is essential for normal craniofacial development. Both excessive RA and RA deficiency in early embryonic stage may lead to a variety of craniofacial malformations, for example, cleft palate, which have been investigated extensively. Dysregulated Wnt and Shh signaling were shown to underlie the pathogenesis of RA-induced craniofacial defects. In our present study, we showed a spatiotemporal-specific effect of RA signaling in regulating early development of facial prominences. Although inhibited Wnt activities was observed in E12.5/E13.5 mouse palatal shelves, early exposure of excessive RA induced Wnt signaling and Wnt-related gene expression in E11.5/E12.5 mouse embryonic frontonasal/maxillary processes. A conserved regulatory network of miR-484-Fzd5 was identified to play critical roles in RA-regulated craniofacial development using RNA-seq. In addition, subsequent osteogenic/chondrogenic differentiation were differentially regulated in discrete mouse embryonic facial prominences in response to early RA induction, demonstrated using both in vitro and in vivo analyses.

Disruption of fos causes craniofacial anomalies in developing zebrafish

Craniofacial development is a complex and tightly regulated process and disruptions can lead to structural birth defects, the most common being nonsyndromic cleft lip and palate (NSCLP). Previously, we identified FOS as a candidate regulator of NSCLP through family-based association studies, yet its specific contributions to oral and palatal formation are poorly understood. This study investigated the role of fos during zebrafish craniofacial development through genetic disruption and knockdown approaches. Fos was expressed in the periderm, olfactory epithelium and other cell populations in the head. Genetic perturbation of fos produced an abnormal craniofacial phenotype with a hypoplastic oral cavity that showed significant changes in midface dimensions by quantitative facial morphometric analysis. Loss and knockdown of fos caused increased cell apoptosis in the head, followed by a significant reduction in cranial neural crest cells (CNCCs) populating the upper and lower jaws. These changes resulted in abnormalities of cartilage, bone and pharyngeal teeth formation. Periderm cells surrounding the oral cavity showed altered morphology and a subset of cells in the upper and lower lip showed disrupted Wnt/β-catenin activation, consistent with modified inductive interactions between mesenchymal and epithelial cells. Taken together, these findings demonstrate that perturbation of fos has detrimental effects on oral epithelial and CNCC-derived tissues suggesting that it plays a critical role in zebrafish craniofacial development and a potential role in NSCLP.

Anatomical network modules of the human central nervous-craniofacial skeleton system

Anatomical network analysis (AnNA) is a systems biological framework based on network theory that enables anatomical structural analysis by incorporating modularity to model structural complexity. The human brain and facial structures exhibit close structural and functional relationships, suggestive of a co-evolved anatomical network. The present study aimed to analyze the human head as a modular entity that comprises the central nervous system, including the brain, spinal cord, and craniofacial skeleton. An AnNA model was built using 39 anatomical nodes from the brain, spinal cord, and craniofacial skeleton. The linkages were identified using peripheral nerve supply and direct contact between structures. The Spinglass algorithm in the igraph software was applied to construct a network and identify the modules of the central nervous system-craniofacial skeleton anatomical network. Two modules were identified. These comprised an anterior module, which included the forebrain, anterior cranial base, and upper-middle face, and a posterior module, which included the midbrain, hindbrain, mandible, and posterior cranium. These findings may reflect the genetic and signaling networks that drive the mosaic central nervous system and craniofacial development and offer important systems biology perspectives for developmental disorders of craniofacial structures.

A novel de novo TP63 mutation in whole-exome sequencing of a Syrian family with Oral cleft and ectrodactyly

BACKGROUND

Oral clefts and ectrodactyly are common, heterogeneous birth defects. We performed whole-exome sequencing (WES) analysis in a Syrian family. The proband presented with both orofacial clefting and ectrodactyly but not ectodermal dysplasia as typically seen in ectrodactyly, ectodermal dysplasia, and cleft lip/palate syndrome-3. A paternal uncle with only an oral cleft was deceased and unavailable for analysis.

METHODS

Variant annotation, Mendelian inconsistencies, and novel variants in known cleft genes were examined. Candidate variants were validated using Sanger sequencing, and pathogenicity assessed by knocking out the tp63 gene in zebrafish to evaluate its role during zebrafish development.

RESULTS

Twenty-eight candidate de novo events were identified, one of which is in a known oral cleft and ectrodactyly gene, TP63 (c.956G > T, p.Arg319Leu), and confirmed by Sanger sequencing.

CONCLUSION

TP63 mutations are associated with multiple autosomal dominant orofacial clefting and limb malformation disorders. The p.Arg319Leu mutation seen in this patient is de novo but also novel. Two known mutations in the same codon (c.956G > A, p.(Arg319His; rs121908839, c.955C > T), p.Arg319Cys) cause ectrodactyly, providing evidence that mutating this codon is deleterious. While this TP63 mutation is the best candidate for the patient's clinical presentation, whether it is responsible for the entire phenotype is unclear. Generation and characterization of tp63 knockout zebrafish showed necrosis and rupture of the head at 3 days post-fertilization (dpf). The embryonic phenotype could not be rescued by injection of zebrafish or human messenger RNA (mRNA). Further functional analysis is needed to determine what proportion of the phenotype is due to this mutation.

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Zebrafish crocc2 mutants exhibit divergent craniofacial shape, misregulated variability, and aberrant cartilage morphogenesis

BACKGROUND

Phenotypic variation is of paramount importance in development, evolution, and human health; however, the molecular mechanisms that influence organ shape and shape variability are not well understood. During craniofacial development, the behavior of skeletal precursors is regulated by both biochemical and environmental inputs, and the primary cilia play critical roles in transducing both types of signals. Here, we examine a gene that encodes a key constituent of the ciliary rootlets, crocc2, and its role in cartilage morphogenesis in larval zebrafish.

RESULTS

Geometric morphometric analysis of crocc2 mutants revealed altered craniofacial shapes and expanded variation. At the cellular level, we observed altered chondrocyte shapes and planar cell polarity across multiple stages in crocc2 mutants. Notably, cellular defects were specific to areas that experience direct mechanical input. Cartilage cell number, apoptosis, and bone patterning were not affected in crocc2 mutants.

CONCLUSIONS

Whereas "regulatory" genes are widely implicated in patterning the craniofacial skeleton, genes that encode "structural" aspects of the cell are increasingly implicated in shaping the face. Our results add crocc2 to this list, and demonstrate that it affects craniofacial geometry and canalizes phenotypic variation. We propose that it does so via mechanosensing, possibly through the ciliary rootlet. If true, this would implicate a new organelle in skeletal development and evolution.

Cell Type- and Tissue-specific Enhancers in Craniofacial Development

The genetic basis of craniofacial birth defects and general variation in human facial shape remains poorly understood. Distant-acting transcriptional enhancers are a major category of non-coding genome function and have been shown to control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development1-3. However, a lack of accurate maps of the genomic location and cell type-specific in vivo activities of all craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combined histone modification and chromatin accessibility profiling from different stages of human craniofacial development with single-cell analyses of the developing mouse face to create a comprehensive catalogue of the regulatory landscape of facial development at tissue- and single cell-resolution. In total, we identified approximately 14,000 enhancers across seven developmental stages from weeks 4 through 8 of human embryonic face development. We used transgenic mouse reporter assays to determine the in vivo activity patterns of human face enhancers predicted from these data. Across 16 in vivo validated human enhancers, we observed a rich diversity of craniofacial subregions in which these enhancers are active in vivo. To annotate the cell type specificities of human-mouse conserved enhancers, we performed single-cell RNA-seq and single-nucleus ATAC-seq of mouse craniofacial tissues from embryonic days e11.5 to e15.5. By integrating these data across species, we find that the majority (56%) of human craniofacial enhancers are functionally conserved in mice, providing cell type- and embryonic stage-resolved predictions of their in vivo activity profiles. Using retrospective analysis of known craniofacial enhancers in combination with single cell-resolved transgenic reporter assays, we demonstrate the utility of these data for predicting the in vivo cell type specificity of enhancers. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Augmentation of bone morphogenetic protein signaling in cranial neural crest cells in mice deforms skull base due to premature fusion of intersphenoidal synchondrosis

Craniofacial anomalies (CFAs) are a diverse group of disorders affecting the shapes of the face and the head. Malformation of the cranial base in humans leads CFAs, such as midfacial hypoplasia and craniosynostosis. These patients have significant burdens associated with breathing, speaking, and chewing. Invasive surgical intervention is the current primary option to correct these structural deficiencies. Understanding molecular cellular mechanism for craniofacial development would provide novel therapeutic options for CFAs. In this study, we found that enhanced bone morphogenetic protein (BMP) signaling in cranial neural crest cells (NCCs) (P0-Cre;caBmpr1a mice) causes premature fusion of intersphenoid synchondrosis (ISS) resulting in leading to short snouts and hypertelorism. Histological analyses revealed reduction of proliferation and higher cell death in ISS at postnatal day 3. We demonstrated to prevent the premature fusion of ISS in P0-Cre;caBmpr1a mice by injecting a p53 inhibitor Pifithrin-α to the pregnant mother from E15.5 to E18.5, resulting in rescue from short snouts and hypertelorism. We further demonstrated to prevent premature fusion of cranial sutures in P0-Cre;caBmpr1a mice by injecting Pifithrin-α through E8.5 to E18.5. These results suggested that enhanced BMP-p53-induced cell death in cranial NCCs causes premature fusion of ISS and sutures in time-dependent manner.

Evolution and dysfunction of human cognitive and social traits: A transcriptional regulation perspective

Evolutionary changes in brain and craniofacial development have endowed humans with unique cognitive and social skills, but also predisposed us to debilitating disorders in which these traits are disrupted. What are the developmental genetic underpinnings that connect the adaptive evolution of our cognition and sociality with the persistence of mental disorders with severe negative fitness effects? We argue that loss of function of genes involved in transcriptional regulation represents a crucial link between the evolution and dysfunction of human cognitive and social traits. The argument is based on the haploinsufficiency of many transcriptional regulator genes, which makes them particularly sensitive to loss-of-function mutations. We discuss how human brain and craniofacial traits evolved through partial loss of function (i.e. reduced expression) of these genes, a perspective compatible with the idea of human self-domestication. Moreover, we explain why selection against loss-of-function variants supports the view that mutation-selection-drift, rather than balancing selection, underlies the persistence of psychiatric disorders. Finally, we discuss testable predictions.

Spatiotemporal single-cell regulatory atlas reveals neural crest lineage diversification and cellular function during tooth morphogenesis

Cranial neural crest cells are an evolutionary innovation of vertebrates for craniofacial development and function, yet the mechanisms that govern the cell fate decisions of postmigratory cranial neural crest cells remain largely unknown. Using the mouse molar as a model, we perform single-cell transcriptome profiling to interrogate the cell fate diversification of postmigratory cranial neural crest cells. We reveal the landscape of transcriptional heterogeneity and define the specific cellular domains during the progression of cranial neural crest cell-derived dental lineage diversification, and find that each domain makes a specific contribution to distinct molar mesenchymal tissues. Furthermore, IGF signaling-mediated cell-cell interaction between the cellular domains highlights the pivotal role of autonomous regulation of the dental mesenchyme. Importantly, we reveal cell-type-specific gene regulatory networks in the dental mesenchyme and show that Foxp4 is indispensable for the differentiation of periodontal ligament. Our single-cell atlas provides comprehensive mechanistic insight into the cell fate diversification process of the cranial neural crest cell-derived odontogenic populations.

Marsupials and Multi-Omics: Establishing New Comparative Models of Neural Crest Patterning and Craniofacial Development

Studies across vertebrates have revealed significant insights into the processes that drive craniofacial morphogenesis, yet we still know little about how distinct facial morphologies are patterned during development. Studies largely point to evolution in GRNs of cranial progenitor cell types such as neural crest cells, as the major driver underlying adaptive cranial shapes. However, this hypothesis requires further validation, particularly within suitable models amenable to manipulation. By utilizing comparative models between related species, we can begin to disentangle complex developmental systems and identify the origin of species-specific patterning. Mammals present excellent evolutionary examples to scrutinize how these differences arise, as sister clades of eutherians and marsupials possess suitable divergence times, conserved cranial anatomies, modular evolutionary patterns, and distinct developmental heterochrony in their NCC behaviours and craniofacial patterning. In this review, I lend perspectives into the current state of mammalian craniofacial biology and discuss the importance of establishing a new marsupial model, the fat-tailed dunnart, for comparative research. Through detailed comparisons with the mouse, we can begin to decipher mammalian conserved, and species-specific processes and their contribution to craniofacial patterning and shape disparity. Recent advances in single-cell multi-omics allow high-resolution investigations into the cellular and molecular basis of key developmental processes. As such, I discuss how comparative evolutionary application of these tools can provide detailed insights into complex cellular behaviours and expression dynamics underlying adaptive craniofacial evolution. Though in its infancy, the field of "comparative evo-devo-omics" presents unparalleled opportunities to precisely uncover how phenotypic differences arise during development.

Association of soluble epoxide hydrolase 2 gene with the risk of non-syndromic cleft lip with or without cleft palate in western Han Chinese population

OBJECTIVES

This study aimed to explore the associations between soluble epoxide hydrolase 2 gene (EPHX2) variants and non-syndromic cleft lip with or without cleft palate (NSCL/P) in Chinese Han population.

METHODS

We recruited 159 NSCL/P cases from Chinese Han population and carried out targeted resequencing using the whole genome sequencing data of 542 healthy Chinese individuals from Novegene internal database as controls. We classified EPHX2 variants as common or rare according to their minor allele frequency and performed an association analysis for common variations and a burden analysis for rare variations.

RESULTS

The lowest P-value in NSCL/P was observed at rs57699806 (P=0.000 13, OR=2.849 and 95% CI: 1.691-4.800), followed by rs4732723 (P=0.006 50, OR=0.662 and 95%CI: 0.491-0.892), rs7829267 (P=0.009 20, OR=1.496 and 95%CI: 1.117-2.005), rs721619 (P=0.011 00, OR=1.474 and 95%CI: 1.098-1.980), and rs7816586 (P=0.040 00, OR=1.310 and 95%CI: 1.015-1.691). The odds ratios suggested the C allele at rs4732723 as a protective factor for NSCL/P and the reference alleles at other single nucleotide polymorphisms (SNPs) as the risk factors for NSCL/P. Burden analysis showed no statistical significance (P>0.05).

CONCLUSIONS

Through targeted resequencing, this study identified five SNPs named rs57699806, rs4732723, rs7829267, rs721619, and rs7816586 around the region of EPHX2 gene associated with NSCL/P in Chinese Han population. Four SNPs of rs57699806, rs4732723, rs7829267, and rs7816586 were first identified.

Spatiotemporal Gene Expression Regions along the Anterior-Posterior Axis in Mouse Embryos before and after Palatal Elevation

The mammalian secondary palate is formed through complex developmental processes: growth, elevation, and fusion. Although it is known that the palatal elevation pattern changes along the anterior-posterior axis, it is unclear what molecules are expressed and whether their locations change before and after elevation. We examined the expression regions of molecules associated with palatal shelf elevation (Pax9, Osr2, and Tgfβ3) and tissue deformation (F-actin, E-cadherin, and Ki67) using immunohistochemistry and RT-PCR in mouse embryos at E13.5 (before elevation) and E14.5 (after elevation). Pax9 was expressed at significantly higher levels in the lingual/nasal region in the anterior and middle parts, as well as in the buccal/oral region in the posterior part at E13.5. At E14.5, Pax9 was expressed at significantly higher levels in both the lingual/nasal and buccal/oral regions in the anterior and middle parts and the buccal/oral regions in the posterior part. Osr2 was expressed at significantly higher levels in the buccal/oral region in all parts at E13.5 and was more strongly expressed at E13.5 than at E14.5 in all regions. No spatiotemporal changes were found in the other molecules. These results suggested that Pax9 and Osr2 are critical molecules leading to differences in the elevation pattern in palatogenesis.