Craniofacial abnormalities often involve sutures, the growth centers of the skull. To characterize the organization and processes governing their development, we profile the murine frontal suture, a model for sutural growth and fusion, at the tissue- and single-cell level on embryonic days (E)16.5 and E18.5. For the wild-type suture, bulk RNA sequencing (RNA-seq) analysis identifies mesenchyme-, osteogenic front-, and stage-enriched genes and biological processes, as well as alternative splicing events modifying the extracellular matrix. Single-cell RNA-seq analysis distinguishes multiple subpopulations, of which five define a mesenchyme-osteoblast differentiation trajectory and show variation along the anteroposterior axis. Similar analyses of in vivo mouse models of impaired frontal suturogenesis in Saethre-Chotzen and Apert syndromes, Twist1+/- and Fgfr2+/S252W, demonstrate distinct transcriptional changes involving angiogenesis and ribogenesis, respectively. Co-expression network analysis reveals gene expression modules from which we validate key driver genes regulating osteoblast differentiation. Our study provides a global approach to gain insights into suturogenesis.

Craniofacial development depends on formation and maintenance of sutures between bones of the skull. In sutures, growth occurs at osteogenic fronts along the edge of each bone, and suture mesenchyme separates adjacent bones. Here, we perform single-cell RNA-seq analysis of the embryonic, wild type murine coronal suture to define its population structure. Seven populations at E16.5 and nine at E18.5 comprise the suture mesenchyme, osteogenic cells, and associated populations. Expression ofHhip, an inhibitor of hedgehog signaling, marks a mesenchymal population distinct from those of other neurocranial sutures. Tracing of the neonatalHhip-expressing population shows that descendant cells persist in the coronal suture and contribute to calvarial bone growth. InHhip^−/−coronal sutures at E18.5, the osteogenic fronts are closely apposed and the suture mesenchyme is depleted with increased hedgehog signaling compared to those of the wild type. Collectively, these data demonstrate thatHhipis required for normal coronal suture development.

Midface dysgenesis is a feature of over 200 genetic conditions, in which upper airway anomalies frequently cause respiratory distress, but its etiology is poorly understood. Mouse models of Apert and Crouzon craniosynostosis syndromes exhibit midface dysgenesis similar to the human conditions. They carry activating mutations of Fgfr2, which is expressed in multiple craniofacial tissues during development. Magnetic resonance microscopy of three mouse models of Apert and Crouzon syndromes revealed decreased nasal passage volume in all models at birth. Histological analysis suggested overgrowth of the nasal cartilage in the two Apert syndrome mouse models. We used tissue-specific gene expression and transcriptome analysis to further dissect the structural, cellular, and molecular alterations underlying midface and upper airway dysgenesis in Apert Fgfr2^+/S252W mutants. Cartilage thickened progressively during embryogenesis due to increased chondrocyte proliferation in the presence of Fgf2. Oral-epithelium expression of mutant Fgfr2, which resulted in a distinctive nasal septal fusion defect, and premature facial suture fusion contributed to the overall dysmorphology. Midface dysgenesis in Fgfr2-craniosynostosis is a complex phenotype arising from the combined effects of aberrant signaling in multiple craniofacial tissues.

The phenotype currently accepted as Pierre Robin syndrome/sequence/anomalad/complex (PR) is characterized by mandibular dysmorphology, glossoptosis, respiratory obstruction, and in some cases, cleft palate. A causative sequence of developmental events is hypothesized for PR, but few clear causal relationships between discovered genetic variants, dysregulated gene expression, precise cellular processes, pathogenesis, and PR-associated anomalies are documented. This review presents the current understanding of PR phenotypes, the proposed pathogenetic processes underlying them, select genes associated with PR, and available animal models that could be used to better understand the genetic basis and phenotypic variation of PR.