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

Single-cell transcriptomic resolution of osteogenesis during craniofacial morphogenesis

Craniofacial morphogenesis depends on complex cell fate decisions during the differentiation of post-migratory cranial neural crest cells. Molecular mechanisms of cell differentiation of mesenchymal cells to developing bones, cartilage, teeth, tongue, and other craniofacial tissues are still poorly understood. We performed single-cell transcriptomic analysis of craniofacial mesenchymal cells derived from cranial NCCs in mouse embryo. Using FACS sorting of Wnt1-Cre2 progeny, we carefully mapped the cell heterogeneity in the craniofacial region during the initial stages of cartilage and bone formation. Transcriptomic data and in vivo validations identified molecular determinants of major cell populations involved in the development of lower and upper jaw, teeth, tongue, dermis, or periocular mesenchyme. Single-cell transcriptomic analysis of Meis2-deficient mice revealed critical gene expression differences, including increased osteogenic and cell adhesion markers. This leads to affected mesenchymal cell differentiation and increased ossification, resulting in impaired bone, cartilage, and tongue formation.

Mesenchymal Wnts are required for morphogenetic movements of calvarial osteoblasts during apical expansion

Apical expansion of calvarial osteoblast progenitors from the cranial mesenchyme (CM) above the eye is integral to calvarial growth and enclosure of the brain. The cellular behaviors and signals underlying the morphogenetic process of calvarial expansion are unknown. Time-lapse light-sheet imaging of mouse embryos revealed calvarial progenitors intercalate in 3D in the CM above the eye, and exhibit protrusive and crawling activity more apically. CM cells express non-canonical Wnt/planar cell polarity (PCP) core components and calvarial osteoblasts are bidirectionally polarized. We found non-canonical ligand Wnt5a-/- mutants have less dynamic cell rearrangements and protrusive activity. Loss of CM-restricted Wntless (CM-Wls), a gene required for secretion of all Wnt ligands, led to diminished apical expansion of Osx+ calvarial osteoblasts in the frontal bone primordia in a non-cell autonomous manner without perturbing proliferation or survival. Calvarial osteoblast polarization, progressive cell elongation and enrichment for actin along the baso-apical axis were dependent on CM-Wnts. Thus, CM-Wnts regulate cellular behaviors during calvarial morphogenesis for efficient apical expansion of calvarial osteoblasts. These findings also offer potential insights into the etiologies of calvarial dysplasias.

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.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Ezh2-dependent methylation in oral epithelia promotes secondary palatogenesis

BACKGROUND

In addition to genomic risk variants and environmental influences, increasing evidence suggests epigenetic modifications are important for orofacial development and their alterations can contribute to orofacial clefts. Ezh2 encodes a core catalytic component of the Polycomb repressive complex responsible for addition of methyl marks to Histone H3 as a mechanism of repressing target genes. The role of Ezh2 in orofacial clefts remains unknown.

AIMS

To investigate the epithelial role of Ezh2-dependent methylation in secondary palatogenesis.

METHODS

We used conditional gene-targeting methods to ablate Ezh2 in the surface ectoderm-derived oral epithelium of mouse embryos. We then performed single-cell RNA sequencing combined with immunofluorescence and RT-qPCR to investigate gene expression in conditional mutant palate. We also employed double knockout analyses of Ezh1 and Ezh2 to address if they have synergistic roles in palatogenesis.

RESULTS

We found that conditional inactivation of Ezh2 in oral epithelia results in partially penetrant cleft palate. Double knockout analyses revealed that another family member Ezh1 is dispensable in orofacial development, and it does not have synergistic roles with Ezh2 in palatogenesis. Histochemistry and single-cell RNA-seq analyses revealed dysregulation of cell cycle regulators in the palatal epithelia of Ezh2 mutant mouse embryos disrupts palatogenesis.

CONCLUSION

Ezh2-dependent histone H3K27 methylation represses expression of cell cycle regulator Cdkn1a and promotes proliferation in the epithelium of the developing palatal shelves. Loss of this regulation may perturb movement of the palatal shelves, causing a delay in palate elevation which may result in failure of the secondary palate to close altogether.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underly facial shape variation, yet how those in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest even during the late migratory phase results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both Tfap2 members dysregulated numerous midface GRN components involved in midface fusion, patterning, and differentiation. Notably, Alx1/3/4 (Alx) transcript levels are reduced, while ChIP-seq analyses suggest TFAP2 directly and positively regulates Alx gene expression. TFAP2 and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish further implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a mutant zebrafish present abnormal alx3 expression patterns, and the two genes display a genetic interaction in this species. Together, these data demonstrate a critical role for TFAP2 in regulating vertebrate midfacial development in part through ALX transcription factor gene expression.

Single cell sequencing of the mouse anterior palate reveals mesenchymal heterogeneity

BACKGROUND

Cleft palate is one of the most prevalent birth defects. Mice are useful for studying palate development because of their morphological and genetic similarities to humans. In mice, palate development occurs between embryonic days (E)11.5 to 15.5. Single cell transcriptional profiles of palate cell populations have been a valuable resource for the craniofacial research community, but we lack a single cell transcriptional profile for anterior palate at E13.5, at the transition from proliferation to shelf elevation.

RESULTS

A detailed single cell RNA sequencing analysis reveals heterogeneity in expression profiles of the cell populations of the E13.5 anterior palate. Hybridization chain reaction RNA fluorescent in situ hybridization (HCR RNA FISH) reveals epithelial populations segregate into layers. Mesenchymal populations spatially segregate into four domains. One of these mesenchymal populations expresses ligands and receptors distinct from the rest of the mesenchyme, suggesting that these cells have a unique function. RNA velocity analysis shows two terminal cell states that contribute to either the proximal or distal palatal regions emerge from a single progenitor pool.

CONCLUSION

This single cell resolution expression data and detailed analysis from E13.5 anterior palate provides a powerful resource for mechanistic insight into secondary palate morphogenesis for the craniofacial research community.

The role of Evi/Wntless in exporting Wnt proteins

Intercellular communication by Wnt proteins governs many essential processes during development, tissue homeostasis and disease in all metazoans. Many context-dependent effects are initiated in the Wnt-producing cells and depend on the export of lipidated Wnt proteins. Although much focus has been on understanding intracellular Wnt signal transduction, the cellular machinery responsible for Wnt secretion became better understood only recently. After lipid modification by the acyl-transferase Porcupine, Wnt proteins bind their dedicated cargo protein Evi/Wntless for transport and secretion. Evi/Wntless and Porcupine are conserved transmembrane proteins, and their 3D structures were recently determined. In this Review, we summarise studies and structural data highlighting how Wnts are transported from the ER to the plasma membrane, and the role of SNX3-retromer during the recycling of its cargo receptor Evi/Wntless. We also describe the regulation of Wnt export through a post-translational mechanism and review the importance of Wnt secretion for organ development and cancer, and as a future biomarker.