Susan M. Motch Perrine

Sequence variations in fibroblast growth factor receptor 2 (FGFR2) are associated with craniosynostosis syndromes, disrupting craniofacial development, yet early mandibular alterations remain poorly characterized. This study examined Meckel's cartilage (MC) and mandibular morphology in mouse models of Apert ( Fgfr2 ^+/P253R ) and Crouzon ( Fgfr2c ^C342Y/+ ) syndromes using Euclidean distance matrix analysis and finite element modeling at embryonic days 14.5 and 16.5. Apert mice exhibited significant MC shortening and elevated, widespread strain in the growth center area correlating with mandibular shortening and asymmetry. In contrast, Crouzon mice showed minor MC changes and localized strain peaks with larger mandibles. These biomechanical perturbations in Meckel's cartilage are associated with distinct mandibular phenotypes, consistent with a bidirectional relationship in which altered strain environments may both influence and reflect mandibular dysmorphogenesis during early development. The findings enhance our understanding of the biomechanical basis of mandibular anomalies associated with FGFR2 mutations, informing future investigations into craniofacial morphogenesis.

Craniofacial development depends on the formation of fibrous joints, or sutures, between skull bones. Premature fusion of sutures, or craniosynostosis, is a common human pathology. Ectopic Hedgehog (HH) signaling is one cause of craniosynostosis. Hhip encodes an inhibitor of HH ligands, and we previously identified coronal suture dysgenesis in embryonic Hhip−/- mice, in which suture mesenchyme was depleted between closely opposed but unfused osteogenic fronts at E18.5. Here, we report that the lambdoid suture fuses in Hhip−/- mice by E18.5. RNA-seq analysis of the Hhip−/- coronal and lambdoid sutures show that HH target gene expression, including Pthlh, is upregulated. Paradoxically, expression of Ihh is downregulated. We hypothesized that PTHLH, a negative regulator of Ihh expression, may reduce HH signaling to promote coronal suture patency and prevent fusion of the Hhip−/- coronal suture. We generated Hhip−/-;Pthlh−/- embryos and found that coronal sutures are fusing by E18.5. Our results reveal a previously undescribed role for Pthlh in suture development and demonstrate suture-specific roles for HH inhibitors in maintaining suture patency.

Osteogenic differentiation is essential for bone development, metabolism, and repair; however, the underlying regulatory relationships among genes remain poorly understood. To elucidate the transcriptomic changes and identify novel regulatory genes involved in osteogenic differentiation, we differentiated mesenchymal stem cells (MSCs) derived from 20 human iPSC lines into preosteoblasts (preOBs) and osteoblasts (OBs). We then performed transcriptome profiling of MSCs, preOBs and OBs. The iPSC-derived MSCs and OBs showed similar transcriptome profiles to those of primary human MSCs and OBs, respectively. Differential gene expression analysis revealed global changes in the transcriptomes from MSCs to preOBs, and then to OBs, including the differential expression of 840 genes encoding transcription factors (TFs). TF regulatory network analysis uncovered a network comprising 451 TFs, organized into five interactive modules. Multiscale embedded gene co-expression network analysis (MEGENA) identified gene co-expression modules and key network regulators (KNRs). From these analyses, KLF16 emerged as an important TF in osteogenic differentiation. We demonstrate that overexpression of Klf16 in vitro inhibited osteogenic differentiation and mineralization, while Klf16+/- mice exhibited increased bone mineral density, trabecular number, and cortical bone area. Our study underscores the complexity of osteogenic differentiation and identifies novel regulatory genes such as KLF16, which plays an inhibitory role in osteogenic differentiation both in vitro and in vivo.

Osteogenic differentiation is essential for bone development, metabolism, and repair; however, the underlying regulatory relationships among genes remain poorly understood. To elucidate the transcriptomic changes and identify novel regulatory genes involved in osteogenic differentiation, we differentiated mesenchymal stem cells (MSCs) derived from 20 human iPSC lines into preosteoblasts (preOBs) and osteoblasts (OBs). We then performed transcriptome profiling of MSCs, preOBs and OBs. The iPSC-derived MSCs and OBs showed similar transcriptome profiles to those of primary human MSCs and OBs, respectively. Differential gene expression analysis revealed global changes in the transcriptomes from MSCs to preOBs, and then to OBs, including the differential expression of 840 genes encoding transcription factors (TFs). TF regulatory network analysis uncovered a network comprising 451 TFs, organized into five interactive modules. Multiscale embedded gene co-expression network analysis (MEGENA) identified gene co-expression modules and key network regulators (KNRs). From these analyses, KLF16 emerged as an important TF in osteogenic differentiation. We demonstrate that overexpression of Klf16 in vitro inhibited osteogenic differentiation and mineralization, while Klf16 ^+/- mice exhibited increased bone mineral density, trabecular number, and cortical bone area. Our study underscores the complexity of osteogenic differentiation and identifies novel regulatory genes such as KLF16 , which plays an inhibitory role in osteogenic differentiation both in vitro and in vivo.

Achondroplasia, the most common chondrodysplasia in humans, is caused by one of two gain of function mutations localized in the transmembrane domain of fibroblast growth factor receptor 3 (FGFR3) leading to constitutive activation of FGFR3 and subsequent growth plate cartilage and bone defects. Phenotypic features of achondroplasia include macrocephaly with frontal bossing, midface hypoplasia, disproportionate shortening of the extremities, brachydactyly with trident configuration of the hand, and bowed legs. The condition is defined primarily on postnatal effects on bone and cartilage, and embryonic development of tissues in affected individuals is not well studied. Using the Fgfr3^Y367C/+ mouse model of achondroplasia, we investigated the developing chondrocranium and Meckel's cartilage (MC) at embryonic days (E)14.5 and E16.5. Sparse hand annotations of chondrocranial and MC cartilages visualized in phosphotungstic acid enhanced three‐dimensional (3D) micro‐computed tomography (microCT) images were used to train our automatic deep learning‐based 3D segmentation model and produce 3D isosurfaces of the chondrocranium and MC. Using 3D coordinates of landmarks measured on the 3D isosurfaces, we quantified differences in the chondrocranium and MC of Fgfr3^Y367C/+ mice relative to those of their unaffected littermates. Statistically significant differences in morphology and growth of the chondrocranium and MC were found, indicating direct effects of this Fgfr3 mutation on embryonic cranial and pharyngeal cartilages, which in turn can secondarily affect cranial dermal bone development. Our results support the suggestion that early therapeutic intervention during cartilage formation may lessen the effects of this condition.

Most bone develops either by intramembranous ossification where bone forms within a soft connective tissue, or by endochondral ossification by way of a cartilage anlagen or model. Bones of the skull can form endochondrally or intramembranously or represent a combination of the two types of ossification. Contrary to the classical definition of intramembranous ossification, we have previously described a tight temporo‐spatial relationship between cranial cartilages and dermal bone formation and proposed a mechanistic relationship between chondrocranial cartilage and dermal bone. Here, we further investigate this relationship through an analysis of how cells organize to form cranial cartilages and dermal bone. Using Wnt1‐Cre2 and Mesp1‐Cre transgenic mice, we determine the derivation of cells that comprise cranial cartilages from either cranial neural crest (CNC) or paraxial mesoderm (PM). We confirm a previously determined CNC‐PM boundary that runs through the hypophyseal fenestra in the cartilaginous braincase floor and identify four additional CNC‐PM boundaries in the chondrocranial lateral wall, including a boundary that runs along the basal and apical ends of the hypochiasmatic cartilage. Based on the knowledge that as osteoblasts differentiate from CNC‐ and PM‐derived mesenchyme, the differentiating cells express the transcription factor genes RUNX2 and osterix (OSX), we created a new transgenic mouse line called R2Tom. R2Tom mice carry a tdTomato reporter gene joined with an evolutionarily well‐conserved enhancer sequence of RUNX2. R2Tom mice crossed with Osx‐GFP mice yield R2Tom;Osx‐GFP double transgenic mice in which various stages of osteoblasts and their precursors are detected with different fluorescent reporters. We use the R2Tom;Osx‐GFP mice, new data on the cell derivation of cranial cartilages, histology, immunohistochemistry, and detailed morphological observations combined with data from other investigators to summarize the differentiation of cranial mesenchyme as it forms condensations that become chondrocranial cartilages and associated dermal bones of the lateral cranial wall. These data advance our previous findings of a tendency of cranial cartilage and dermal bone development to vary jointly in a coordinated manner, promoting a role for cranial cartilages in intramembranous bone formation.

The cranial endo and dermal skeletons, which comprise the vertebrate skull, evolved independently over 470 million years ago and form separately during embryogenesis. In mammals, much of the cartilaginous chondrocranium is transient, undergoing endochondral ossification or disappearing, so its role in skull morphogenesis is not well studied and it remains an enigmatic structure. We provide complete 3D reconstructions of the laboratory mouse chondrocranium from embryonic day (E) 13.5 through E17.5 using a novel methodology of uncertainty-guided segmentation of phosphotungstic enhanced 3D micro-computed tomography images with sparse annotation. We evaluate the embryonic mouse chondrocranium and dermatocranium in 3D, and delineate the effects of a Fgfr2 variant on embryonic chondrocranial cartilages and on their association with forming dermal bones using the Fgfr2c^C342Y/+ Crouzon syndrome mouse. We show that the dermatocranium develops outside of and in shapes that conform to the chondrocranium. Results reveal direct effects of the Fgfr2 variant on embryonic cartilage, on chondrocranium morphology, and on the association between chondrocranium and dermatocranium development. Histologically, we observe a trend of relatively more chondrocytes, larger chondrocytes, and/or more matrix in the Fgfr2c^C342Y/+ embryos at all timepoints before the chondrocranium begins to disintegrate at E16.5. The chondrocrania and forming dermatocrania of Fgfr2c^C342Y/+ embryos are relatively large, but a contrasting trend begins at E16.5 and continues into early postnatal (P0 and P2) timepoints, with the skulls of older Fgfr2c^C342Y/+ mice reduced in most dimensions compared to Fgfr2c^+/+ littermates. Our findings have implications for the study and treatment of human craniofacial disease, for understanding the impact of chondrocranial morphology on skull growth, and potentially on the evolution of skull morphology.

The Fgfr2c^C342Y/+ Crouzon syndrome mouse model carries a cysteine to tyrosine substitution at amino acid position 342 (Cys342Tyr; C342Y) in the fibroblast growth factor receptor 2 (Fgfr2) gene equivalent to a FGFR2 mutation commonly associated with Crouzon and Pfeiffer syndromes in humans. The Fgfr2c C342Y mutation results in constitutive activation of the receptor and is associated with upregulation of osteogenic differentiation. Fgfr2c^C342Y/+ Crouzon syndrome mice show premature closure of the coronal suture and other craniofacial anomalies including malocclusion of teeth, most likely due to abnormal craniofacial form. Malformation of the mandible can precipitate a plethora of complications including disrupting development of the upper jaw and palate, impediment of the airway, and alteration of occlusion necessary for proper mastication. The current paradigm of mandibular development assumes that Meckel’s cartilage (MC) serves as a support or model for mandibular bone formation and as a template for the later forming mandible. If valid, this implies a functional relationship between MC and the forming mandible, so mandibular dysmorphogenesis might be discerned in MC affecting the relationship between MC and mandibular bone. Here we investigate the relationship of MC to mandible development from the early mineralization of the mandible (E13.5) through the initiation of MC degradation at E17.7 using Fgfr2c^C342Y/+ Crouzon syndrome embryos and their unaffected littermates (Fgfr2c^+/+). Differences between genotypes in both MC and mandibular bone are subtle, however MC of Fgfr2c^C342Y/+ embryos is generally longer relative to unaffected littermates at E15.5 with specific aspects remaining relatively large at E17.5. In contrast, mandibular bone is smaller overall in Fgfr2c^C342Y/+ embryos relative to their unaffected littermates at E15.5 with the posterior aspect remaining relatively small at E17.5. At a cellular level, differences are identified between genotypes early (E13.5) followed by reduced proliferation in MC (E15.5) and in the forming mandible (E17.5) in Fgfr2c^C342Y/+ embryos. Activation of the ERK pathways is reduced in the perichondrium of MC in Fgfr2c^C342Y/+ embryos and increased in bone related cells at E15.5. These data reveal that the Fgfr2c C342Y mutation differentially affects cells by type, location, and developmental age indicating a complex set of changes in the cells that make up the lower jaw.

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.

The cranial endo- and dermal skeletons, which comprise the vertebrate skull, evolved independently over 470 million years ago and form separately during embryogenesis. In mammals, much of the cartilaginous chondrocranium is transient, undergoing endochondral ossification or disappearing, so its role in skull morphogenesis is not well studied and it remains an enigmatic structure. We provide complete three-dimensional (3D) reconstructions of the laboratory mouse chondrocranium from embryonic day 13.5 through 17.5 using a novel methodology of uncertainty-guided segmentation of phosphotungstic enhanced 3D microcomputed tomography images with sparse annotation. We evaluate the embryonic mouse chondrocranium and dermatocranium in 3D and delineate the effects of a Fgfr2 variant on embryonic chondrocranial cartilages and on the association with forming dermal bones using the Fgfr2c ^C342Y/+ Crouzon syndrome mouse. We show that the dermatocranium develops outside of and in shapes that conform to the chondrocranium. Results reveal direct effects of the Fgfr2 variant on embryonic cartilage, on chondrocranium morphology, and on the association between chondrocranium and dermatocranium development. Histologically we observe a trend of relatively more chondrocytes, larger chondrocytes, and/or more matrix in the Fgfr2c ^C342Y/+ embryos at all timepoints before the chondrocranium begins to disintegrate at E16.5. The chondrocrania and forming dermatocrania of Fgfr2c ^C342Y/+ embryos are relatively large, but a contrasting trend begins at E16.5 and continues into early postnatal (P0 and P2) timepoints, with the skulls of older Fgfr2c ^C342Y/+ mice reduced in most dimensions compared to Fgfr2c ^+/+ littermates. Our findings have implications for the study and treatment of human craniofacial disease, for understanding the impact of chondrocranial morphology on skull growth, and potentially on the evolution of skull morphology.