Meckel’s Cartilage in Mandibular Development and Dysmorphogenesis

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.

1 Introduction

Crouzon syndrome, with an estimated birth prevalence of 16.5/1,000,000 (Cohen and Kreiborg, 1992) is possibly the mildest of FGFR2 craniosynostosis syndromes (Azoury et al., 2017). In addition to premature closure of cranial sutures, Crouzon syndrome presents with facial anomalies (retrusion, beaked nose, ocular proptosis), absence of major hand and foot abnormalities (Johnson and Wilkie, 2011; Lu et al., 2020), and mandibular dysgenesis (Costaras-Volarich and Pruzansky, 1984; Boutros et al., 2007; Lu et al., 2019). Change in the spatial relationship of the mandible to the cranial base develops independently of mandibular shape deformity, and may influence functional capabilities including eating and speech (Lu et al., 2019). Additional dysmorphologies of the mandible and teeth remain poorly characterized for this condition, although it has been proposed significant mandibular morphological and growth abnormalities occur in children with Crouzon syndrome (Costaras-Volarich and Pruzansky, 1984; Boutros et al., 2007). Development of the lower jaw begins much earlier than these studies have been able to describe and includes growth of Meckel’s cartilage prior to mineralization of the mandible. We propose embryonic development of the lower jaw, including Meckel’s cartilage and mandibular bone, may be affected in Crouzon syndrome.

The FGFR2c C342Y mutation associated with Crouzon syndrome results in constitutive activation of the receptor and is most commonly associated with up-regulation of osteogenic differentiation (Yu et al., 2003; Perlyn et al., 2006; Heuzé et al., 2014). Mouse models of the FGFR2c C342Y mutation have been used to study the development of the Crouzon phenotype in multiple tissues of the head (Perlyn et al., 2006; Martínez-Abadías et al., 2013; Motch Perrine et al., 2017; Lee et al., 2018; Motch Perrine et al., 2019). However, this mutation can also affect other cell types including mesenchymal stem cells, (Baddoo et al., 2003; Miraoui et al., 2009; Martínez-Abadías et al., 2010; Maruyama et al., 2010; Deckelbaum et al., 2012; Yeh et al., 2012; Martínez-Abadías et al., 2013; Chen et al., 2014; Xu et al., 2017). FGFR activation induces major intracellular signaling pathways including RAS-MAPK (ERK1/2), PI3K-AKT, PLCγ-PKC, and signal transducer and activator of transcription (STAT) (Ornitz and Itoh, 2015). ERK pathway is active in the osteogenic progenitors and peripheral chondrocytes of MC (Parada et al., 2015) and is essential for the formation of the lower jaw. In Wnt1-Cre;Erk2^fl/fl mice, NCC-specific conditional knockout of Erk2 disrupts mandibular osteogenic differentiation resulting in micrognathia (Parada et al., 2015). In addition, conditional overexpression of Fgfr2c causes upregulation of phosphorylates ERK (p-ERK), resulting in craniofacial hypoplasia including reduction of the mandible in murine models (Lee et al., 2018). Precise regulation of ERK pathway appears critical for normal mandibular development, but its involvement in development of the lower jaw of Fgfr2c^C342Y/+ Crouzon mice has not been elucidated.

Meckel’s cartilage develops from the most cranial of the transient pharyngeal arches that comprise a series of repeated outgrowths on the lateral sides of the embryonic head. The most cranial, or first pharyngeal arch (PA1) organizes into two separate compartments in jawed vertebrates: the maxillary process cranially and the mandibular arch caudally. NCC from the mandibular arch condense and then differentiate into chondrocytes to form Meckel’s cartilage (MC), a rod-like structure providing shape to the lower jaw. The same NCC that form MC give rise to the dermal bones of the upper and lower jaws and at least part of two of the three middle ear bones (malleus, incus). The segmental arrangement and structural organization of the pharyngeal arches is highly conserved throughout evolution, though the number of PAs varies (Graham and Richardson, 2012; Shone et al., 2016), suggesting conservation of the molecular mechanisms and signals that govern their development and their contribution to the pharyngeal skeleton (Frisdal and Trainor, 2014).

The current paradigm of mandibular development assumes that Meckel’s cartilage (MC) functions as a model during early mandibular development and as a template for growth of the mandible later in development (Svandova et al., 2020). During development, the distal portion of MC becomes the symphysis, a point of fusion of the left and right hemi-mandibles, and the proximal portion of MC forms cartilages of the malleus and incus, bones that ossify endochondrally (McKenzie, 1958; Amin and Tucker, 2006; Frisdal and Trainor, 2014; Svandova et al., 2020). Mineralization of the mandibular bone begins at E14 laterally along MC and eventually encases it (Ramaesh and Bard, 2003). Part of the Mid and/or Posterior regions of MC (Figure 1) transforms into ligament, but other portions of the Mid region degrade and eventually completely disappear (Ishizeki et al., 2001; Svandova et al., 2020; Nakamura et al., 2021).

FIGURE 1

FIGURE 1. Schematic of MC in mandible identifying regions for analysis. Example histological sections (top) associated with developing mandible and MC (bottom) (after Svandova et al., 2020). In the histological sections, mandibular bone matrix (MB) is stained green and is adjacent to the red-orange MC stained using a standard Safranin-O protocol. Below, the developing mandible (tan), segmented from microCT images, is superimposed with MC (blue) and mandibular teeth (red) segmented from PTA-enhanced microCT images of an E17.5 Fgfr2c^+/+ embryo. Note that the color of MC and MB is not consistent between histological images (top) and 3D models created from PT-enhanced microCT images (bottom). MC is divided into 4 equal regions from anterior to posterior with Anterior and Mid 1 containing the break in MC at E17.5 associated with the root of the incisor (I). The incus is not included in this model.

The scientific premise that MC and mandibular development are tightly coordinated implies that the relationship between MC and the forming mandible is critical and that dysmorphogenesis of either (or both) structure(s) can lead to anatomical abnormalities (Ishizeki et al., 1999; Manocha et al., 2019; Svandova et al., 2020). Typically, cells within each structure behave predictably and in an apparently coordinated fashion to create a fully functional lower jaw. Insults, including variants of the FGF/FGFR signaling pathway, can cause significant changes to cells in MC and the developing mandible that may affect the form and function of the mandible. Here, we use Fgfr2c^C342Y/+ Crouzon syndrome mice (Eswarakumar et al., 2004) that have characteristically shortened faces and commonly exhibit aberrant occlusion of the jaws. These mice are known to have abnormal osteogenesis and chondrocyte proliferation, as well as mandibular dysmorphology at birth (Eswarakumar et al., 2004; Holmes et al., 2018; Motch Perrine et al., 2019) and provide an appropriate model for investigating the direct effects of the FGFR2c C342Y mutation on MC and mandibular development. We hypothesize that this mutation will drive alterations in cell activity and structure within MC and the developing mandible resulting in aberrant mandibular form in embryonic Fgfr2c^C342Y/+ mice.

2 Materials and Methods

2.1 Mouse Model

The Fgfr2c^C342Y/+ mice were a generous gift from Dr. Eswarakumar, who previously described the generation of this model (Eswarakumar et al., 2004). The model was maintained on a CD1 background for viability and breeding. Embryos were collected at embryonic days E13.5, 15.5, and 17.5. PCR for sex and genotype was performed using a tail snipped at collection. We do not anticipate any differences due to sexual dimorphism based on our previous work on multiple mouse models (Motch Perrine et al., 2017; Motch Perrine et al., 2019; Lesciotto et al., 2022), still we have indicated the balance of males and females used in this study in Table 1. Mouse litters were produced in compliance with animal welfare guidelines approved by Icahn School of Medicine at Mount Sinai and Pennsylvania State University Animal Care and Use Committees (ISMMS #07-0757 and PSU #46558). To more precisely define developmental age, E13.5 and E15.5 embryos were staged using the Embryonic Mouse Ontogenetic Staging System, EMOSS (Musy et al., 2018). E17.5 embryos have developed beyond the capabilities of this system so are considered to be an average of 17.5 days post conception. Average developmental age at E13.5 was 341.14 h for Fgfr2c^C342Y/+ embryos and 338.4 h for Fgfr2c^+/+ embryos. At E15.5, average developmental age was 371.93 h for Fgfr2c^C342Y/+ embryos and 372.11 h for Fgfr2c^+/+ embryos.

TABLE 1

TABLE 1. Specimen demographics.

2.2 Imaging Protocols

Phosphotungstic acid (PTA)-enhanced Microcomputed Tomography (microCT) images for MC and microCT images for bone analyses were acquired by the Center for Quantitative Imaging at the Pennsylvania State University (https://iee.psu.edu/labs/center-quantitative-imaging) using the General Electric v|tom|x L300 nano/microCT system. Five Fgfr2c^C342Y/+ and 5 Fgfr2c^+/+ embryos were scanned at E13.5, 5 Fgfr2c^C342Y/+ and 6 Fgfr2c^+/+ embryos were scanned at E15.5, and 4 Fgfr2c^C342Y/+ and 5 Fgfr2c^+/+ embryos were scanned at E17.5. Image data were reconstructed on a 2,024 × 2,024 pixel grid as a 32-bit volume, and were reduced to 16-bit volume for image analysis using Avizo 2019.3 and 2020.2 (Thermo Fisher Scientific). Scanning parameters varied from 60–100 kV and 75–170 μA, to accommodate age group and type of scan performed. Voxel sizes ranged from 6.9 to 15 microns (µm) for bone scans and 4.5–8 µm for PTA-enhanced scans.

2.3 Segmentation and Quantification of Cartilage

3D models of MC were created from PTA-enhanced microCT (Lesciotto et al., 2020) images of specimens at E13.5, E15.5 and E17.5 (Table 1) using a previously reported automatic segmentation approach (Zheng et al., 2020). Due to the complex 3D cartilage structures and variations in our images as well as large volumes of image data, an effective and efficient automatic method is needed for our cartilage segmentation task. Recently, deep learning methods have achieved great successes in biomedical image segmentation (Ronneberger et al., 2015). Compared with common CT images, each of our 3D microCT volumes is very large (about 1,500 × 2,000 × 2,000 voxels per image), and consequently, it will be extremely difficult to rely on manual annotation to produce a sufficiently large amount of labeled data for training a deep learning segmentation model. Hence, our segmentation approach must be able to accommodate a large image size and make the most out of sparsely labeled training data. First, we chose a very sparse subset of 2D image slices (e.g., 2%–10%) in our 3D microCT training volumes for manual annotation that represents and covers the rest (unannotated) 2D slices of the training volumes well. Then, sparse manual annotations were conducted by experts (MKP and SMMP) using Avizo 2020.2 (Thermo Fisher Scientific). Next, the annotated 2D slices were used to train a judiciously designed fully convolutional network (FCN) model to identify target objects (called pseudo-labels or PLs) in all the unannotated 2D slices of the training volumes. In this way, we sought to bridge the large gap between our sparse annotation and full annotation of the entire 3D training volumes. Since our FCN model is not trained with a standard full annotation protocol (due to the sparse annotation used), its identified PLs on the unannotated 2D slices may be unreliable and noisy. Thus, we estimated the reliability of the PLs by computing the associated uncertainty maps of the PLs, which quantify the pixel-wise prediction confidence. With the expanded training set (i.e., 3D images with both manual labels and pseudo labels) and guided by the uncertainty maps, the FCN model was trained iteratively to distill more stable knowledge about the training data, thus becoming more robust and generalizing better to unseen data. Finally, the segmentation results were integrated along the three orthogonal planes (orthogonal to each of the x, y, and z directions) to further boost the overall segmentation accuracy (by average ensemble of segmentation results along the three directions). Results of automated segmentation were reviewed by an expert anatomist and segmentation was corrected manually if necessary. Volume and area were estimated using the Material Statistics of Avizo 2020.2 (ThermoFisher Scientific). Volume and area were compared between the genotypes within age categories with independent-samples Mann Whitney U test performed in IBM SPSS Statistics v 27 (IBM Corp). Once created, the 3D models of MC were quantified by collecting the 3D coordinates of landmarks designed for age-specific analyses (Figure 2; Table 2). Landmark collection as well as volume and area calculations were performed in Avizo 2019.3 and 2020.2 (Thermo Fisher Scientific). Each sample was digitized twice by the same observer, checked and corrected for gross error, and remaining measurement error was minimized by averaging the coordinates of the two trials. A maximum of 5% intraobserver error in landmark placement was accepted. Morphological differences were assessed as described in Section 2.5. PTA-enhanced microCT images were used to hand segment tooth germs visualized in Figure 1.

FIGURE 2

FIGURE 2. 3D landmarks on developing MC and mandible. (A) Two medial landmarks and 8, 10, or 14 bilateral landmarks shown on MC at E13.5 (top), E15.5 (middle) and E17.5 (bottom) embryos respectively. Lateral views (buccal) are displayed on the left and superior views are on the right. (B) Ten bilateral landmarks are shown on the left hemi-mandible of an E17.5 embryo. Lateral views are displayed on the left (buccal view (top) and lingual view (bottom)) and superior view of the left and right hemi-mandibles are shown on the right. Scale = 1.0 mm.

TABLE 2

TABLE 2. Definitions of MC landmarks.

2.4 Segmentation and Quantification of Bone

Isosurfaces of mandibular bone were created by automated segmentation of bone from 16-bit microCT images using a minimum threshold of 70–100 mg/cm³ partial density hydroxyapatite based on phantoms scanned with the samples using Avizo 2020.2 (ThermoFisher Scientific). Data were processed using a median filter in Avizo 2020.2 (ThermoFisher Scientific) to remove noise and reviewed by an expert anatomist. No manual corrections to form were necessary for MB. Volume and area were estimated for each specimen using the Material Statistics module of Avizo 2020.2 (Thermo Fisher Scientific). Volume and area were compared between the genotypes within age categories with independent-samples Mann Whitney U test performed in IBM SPSS Statistics v 27 (IBM Corp). These isosurfaces were used to assess mandibular form by recording the 3D coordinates of 14 (7 bilateral) landmarks at E15.5 and 20 (10 bilateral) landmarks at E17.5 (Figure 2; Table 3) in Avizo 2019.3 (Thermo Fisher Scientific). Each sample was digitized twice by the same observer, checked and corrected for gross error, and measurement error was minimized by averaging the coordinates of the two trials. A maximum of 5% intraobserver error in landmark placement was accepted. Morphological differences were assessed as described in Section 2.5. Sample sizes for each genotype and embryonic age are listed in Table 1.

TABLE 3

TABLE 3. Definitions of mandible bone landmarks.

2.5 3D Morphometric Analysis of Meckel’s Cartilage and Mandible

To statistically determine differences in MC and mandibular form between genotypes, we used Euclidian Distance Matrix Analysis (EDMA) (Lele and Richtsmeier, 1995; Lele and Richtsmeier, 2001). EDMA converts 3D landmark data into a matrix of all possible linear distances between unique landmark pairs and tests for statistical significance of differences between shapes using a boot-strapped hypothesis testing procedure. Differences of specific linear distances are statistically evaluated by a 90% confidence interval produced through a non-parametric bootstrapping procedure (Lele and Richtsmeier, 1995). Rejection of the null hypothesis of similarity for linear distances enables localization of differences to specific dimensions. EDMA analyses were performed using WinEDMA (Cole, 2002).

2.6 Superimposition of PTA-Enhanced and Bone Targeting MicroCT Images

Scans targeting soft tissues such as MC (PTA-enhanced microCTs) and scans targeting hard tissues such as MB were performed on each specimen as described in Section 2.2. MB is very difficult to fully segment from PTA-enhanced microCTs, so therefore both scan types were used. For figures with superimposition of MB and MC, isosurfaces of MC and MB were created as described in Sections 2.3,2.4, respectively. Landmarks 1, 2, and 15 (Table 3) were taken on MB in both scans in Avizo 2020.2 (ThermoFisher Scientific) for each individual for which a superimposition was required. A rigid transformation, which moves the points of the first set as close as possible onto the points of the second set (the sum of the squared distances between corresponding points is minimized) was computed from these landmarks for aligning the model (bone microCT scan) to the reference (PTA-enhanced microCT scan) in Avizo 2020.2 (ThermoFisher Scientific). The superimposition was adjusted slightly by hand if necessary following landmark transformation of images. Superimpositions were used only for figure creation. Analysis of individual tissues was performed on individual scans targeting the tissue of interest.

2.7 Histology

Specimen were collected (n = 3 per age per genotype) and fixed in 4% paraformaldehyde for 24 h. Care was taken to balance between sexes, and E17.5 specimens were decalcified in 0.25 M EDTA at pH 7.4 for 12 days. Isolated heads were processed for paraffin-based histology per established protocol (Behringer et al., 2013; Motch Perrine et al., 2021). Using a rotary microtome, 7 µm serial, coronal sections were cut and mounted on Superfrost-Plus slides (Thermo Fisher Scientific), and the range of slides containing Meckel’s cartilage was determined for each specimen. Slide ranges were then divided into four equal regions (Anterior, Mid 1, Mid 2, Posterior, Figure 1), and one slide was selected per region per specimen for each stain. For immunohistochemistry, slides were subjected to epitope retrieval using an autoclave (121°C, 10 min) and sodium citrate buffer pH 6.0 before blocking with 3% hydrogen peroxide and 1% goat serum/bovine serum albumin. Sections were incubated with primary antibodies at room temperature for 2 h, Proliferating Cell Nuclear Antigen (PCNA, AbCam, ab18197, 1:2,000) or for 1 h, alkaline phosphatase (ALP, ab108377, 1:500). After washing in PBS, slides were incubated with HRP conjugated secondary antibody for one hour at room temperature (ab6721, 1:500) and Diaminobenzidine (DAB, Vector Laboratories) chromagen was used to identify immunoreactive structures (Durham et al., 2019). Safranin-O and tartrate-resistant acid phosphatase (TRAP) staining were performed per established protocol (Howie et al., 2018). Slides were imaged using a Leica BX50 microscope, DFC450 camera, and LAS-X imaging software (Leica Biosystems). Regions of interest (MC, mandible) were identified and analyzed using ImageJ color deconvolution and masks to count stained areas by color (Supplementary Figure S1) (Ruifrok and Johnston, 2001; Haub and Meckel, 2015; Durham et al., 2018; Landini et al., 2021). Careful gating for cell size (3–30 µm) was used to reduce noise in the form of stained areas smaller or larger than a cell. The same gating parameters were applied to every image analyzed. At least two images were captured per slide for each stain for each specimen investigated. Osteoclasts were identified by TRAP positivity and multiple nuclei and were counted by a single investigator (ELD) two times with significant correlation between counts. Counts were averaged for investigation. Non-parametric Mann-Whitney U tests were used where necessary to compare genotypes at each age, and then by each region using SPSS 25 software (IBM).

2.8 Immunofluorescence

Specimens (n = 3 per genotype) collected at E15.5 were fixed in 4% paraformaldehyde overnight, embedded in OCT compound, and frozen. Samples were sectioned at 10 µm of thickness on the cryostat. For immunofluorescence, sections were permeabilized in PBST (0.2% Triton X 100 in PBS) for 5 min. Antigen retrieval was performed at 100°C for 10 min in citrate buffer (Sigma-Aldrich, C9999). Sections were then blocked with 3% bovine serum albumin in PBS. Anti-p-ERK1/2 (Cell Signaling Technology, 4370) was diluted at 1:100 with PBS and incubated on the sections for 1 h. After washing with PBS, secondary antibody (1:500, Thermofisher, A32732) was incubated for 1 h. The sections were rinsed briefly in PBS and stained with 0.5 μg/ml 4′,6-Diamidino-2-Phenylindole (DAPI; Thermo Fisher Scientific) in PBS before mounting in anti-fade mounting medium (Vector Labs, H-1700). Slides were imaged using a Nikon Eclipse E 600 fluorescence microscope. Images were analyzed using ImageJ.

3 Results

3.1 Orientation of Meckel’s Cartilage and Mandible During Development

Visual inspection of superimpositions of bone and cartilage segmentations indicated only subtle morphological differences between Fgfr2c^C342Y/+ and Fgfr2c^+/+ embryos from E13.5 through E17.5 (Figure 3). Very thin areas of bone with little mineralization may not be detected at the indicated bone threshold in microCT images, as compared to traditional histological staining techniques.

FIGURE 3

FIGURE 3. Relative orientation of MC and developing mandible in Fgfr2c^+/+ and Fgfr2c^C342Y/+ embryos. Segmentations of multimodal images of MC (blue) and mandibular bone (tan) of developing mandible at E13.5 (A), E15.5 (B), and E17.5 (C) with representative Fgfr2c^+/+ embryos displayed in the top row and Fgfr2c^C342Y/+ Crouzon syndrome embryos displayed in the bottom row of each age group. Lateral views (buccal) are displayed on the left and superior views on the right. Not enough bone is mineralized at E13.5 to be detected by microCT. Scale = 1.0 mm.

3.2 Impact of Fgfr2c C342Y Mutation on Meckel’s Cartilage

At E13.5, Fgfr2c^C342Y/+ and Fgfr2c^+/+ embryos do not exhibit significant differences in MC form. However, at E15.5 MC volume is significantly larger in these embryos as compared to Fgfr2c^+/+ littermates (Table 4) and MC is significantly longer in the anteroposterior dimension (Figure 4A). MC of Fgfr2c^C342Y/+ embryos is not larger overall at E17.5 but linear distances characterizing the most anterior and most posterior aspects of MC are larger relative to unaffected littermates (Figure 4B).

TABLE 4

TABLE 4. Results of Mann-Whitney tests and descriptive statistics for cartilage and bone volume by age and genotype.

FIGURE 4

FIGURE 4. Morphological differences of MC in Fgfr2c^C342Y/+ embryos and Fgfr2c^+/+ littermates at E15.5 and E17.5. Linear distances that are statistically significantly different between genotypes by confidence interval testing (α = 0.10) are shown for E15.5 (A) and E17.5 (B). Thin yellow lines indicate linear distances that are between 5% and 10% larger in Fgfr2c^C342Y/+ embryos as compared to Fgfr2c^+/+ embryos and thicker orange lines indicate linear distances that are between 10% and 20% larger in Fgfr2c^C342Y/+ embryos as compared to Fgfr2c^+/+ embryos. For orientation, superimpositions of MC and mandibular bone are inset. Lateral views (buccal) are displayed on the left and superior views are on the right. Scale = 1.0 mm.

Effects of the Fgfr2c C342Y mutation on chondrocytes within MC are not qualitatively apparent, especially at the earliest time points investigated here (Figure 5A). Since asymmetry was not noted in our assessment of MC form, we did not separate left and right in our histological assessment. The area of MC measured in coronal sections of E13.5 embryos revealed increased area of MC in Fgfr2c^C342Y/+ Crouzon syndrome embryos (p = 0.001) as compared to Fgfr2c^+/+ littermates. This was driven by significantly increased area in Mid 1 (p = 0.003) and Posterior (p ≤ 0.001) regions of MC in Fgfr2c^C342Y/+ embryos at E13.5. No statistically significant differences in cross-sectional area were identified at the later timepoints (E15.5, E17.5) (Figure 5B). Since cartilage growth can result from increased matrix deposition (Figure 5C), net increase in cells (Figure 5D), or an increase in chondrocyte hypertrophy (cell diameter) (Figure 5E), we investigated these characteristics in MC and identified more matrix area in Fgfr2c^C342Y/+ embryos at E13.5 (p ≤ 0.001) relative to Fgfr2c^+/+littermates. This was the result of greater matrix in all regions in Fgfr2c^C342Y/+ embryos, but Mid 1 (p = 0.003), Mid 2 (p = 0.001), and Posterior (p = 0.001) regions had statistically significantly more cartilage matrix as compared to Fgfr2c^+/+ littermates (Figure 5C). Interestingly, there were no differences in chondrocyte number between genotypes at E13.5 or E15.5. The number of chondrocytes within MC was statistically significantly decreased in Posterior of Fgfr2c^C342Y/+ embryos at E17.5 (Figure 5D). Our assessment of the diameter of chondrocytes within MC indicated no differences by genotype at E13.5 and E15.5, however there were significantly smaller cells in Anterior of MC in Fgfr2c^C342Y/+ embryos at E17.5 (Figure 5E). This result correlates with the trend toward more cells in Fgfr2c^C342Y/+ embryos in Anterior at E17.5 (Figure 5D).

FIGURE 5

FIGURE 5. Histomorphometric analysis of MC in Fgfr2c^C342Y/+ embryos and Fgfr2c^+/+ littermates, E13.5, E15.5, E17.5. (A) Representative Safranin O-stained sections of MC and surrounding mandible bone (MB) in the coronal plane from Fgfr2c^+/+ and Fgfr2c^C342Y/+ embryos at E13.5 (top two rows), E15.5 (middle two rows) and E17.5 (bottom two rows). Chondrocytes and matrix are identified by the red-orange Safranin O stain with green bone matrix surrounding (MB). The incisor (I) is adjacent to MC and MB in Anterior and Mid 1. Scale = 0.10 mm. (B) Assessment of cross-sectional area of MC indicates some differences by genotype early in development (E13.5, top) with Fgfr2c^C342Y/+ embryos having a relatively greater area overall (box plot, left) and at specific regions in the antero-posterior dimension (line plots, right). (C) Similarly, more cartilage matrix was identified at E13.5 (top) in Fgfr2c^C342Y/+ embryos with no differences by genotype being noted later in development [E15.5 (middle), E17.5 (bottom)]. (D) The number of chondrocytes within MC did not vary by genotype at either E13.5 (top) or E15.5 (middle), but fewer chondrocytes were counted in Posterior of Fgfr2c^C342Y/+ embryos at E17.5 (bottom). (E) The average diameter of chondrocytes varied little from anterior to posterior aspects of MC and not at all by genotype at E13.5 (top) and E15.5 (middle). Chondrocytes in Anterior of E17.5 Fgfr2c^C342Y/+ embryos were smaller than those in Fgfr2c^+/+ embryos (bottom two rows). White bar in box plot identifies median. X in box plot identifies mean. n = 3 individuals per genotype per age with at least 4 images per region per embryo being analyzed. **p ≤ 0.01, ***p ≤ 0.001.

3.3 Fgfr2c Mutation Affects Proliferation in and Around Meckel’s Cartilage

Because we identified some significant differences in cell number between genotypes, we investigated proliferation to see if the mutation influenced chondrocyte and/or bone lineage cells during initial ossification. E15.5 and E17.5 embryos were used because the bone was robust enough at these ages for analysis. Differences in proliferation in MC are subtle (Figure 6A; Supplementary Figure S1), however an overall decrease in proliferation in chondrocytes from Fgfr2c^C342Y/+ embryos was identified at E15.5 (Figure 6B, top) relative to Fgfr2c^+/+ littermates. This relationship was no longer significant by E17.5 (Figure 6B, row 2). Within the forming mandible adjacent to MC, proliferation was increased in Anterior, and decreased in Posterior of E15.5 Fgfr2c^C342Y/+ embryos relative to Fgfr2c^+/+ littermates (Figure 6B, row 3). At E17.5 an overall decrease in proliferation in the mandibles of Fgfr2c^C342Y/+ embryos was driven by significantly less proliferation in Posterior (Figure 6B, bottom).

FIGURE 6

FIGURE 6. Effect of Fgfr2c C342Y mutation on proliferation in and around MC in E15.5 and E17.5 embryos. (A) Representative coronal sections of MC and surrounding mandible bone (MB) in Fgfr2c^+/+ and Fgfr2c^C342Y/+ embryos stained with Proliferating Cell Nuclear Antigen (PCNA, brown stain) counterstained with hematoxylin (purple) at E15.5 (top rows) and E17.5 (bottom rows). Scale = 0.1 mm. Note that due to variation in the location of the natural break in MC, only a small remnant of MC remains in Mid 1 at E17.5 in the displayed image of an Fgfr2c^+/+ embryo. (B) Quantification of proliferating cells in MC (top 2 rows) and mandible (bottom 2 rows) standardized by the total number of cells within each region. Proliferation was reduced in MC of Fgfr2c^C342Y/+ embryos at E15.5 and in mandible of Fgfr2c^C342Y/+ embryos at E17.5. White bar in box plot identifies median. X in box plot identifies mean. n = 3 individuals per age per genotype with at least 4 images per region per embryo being analyzed. *p ≤ 0.05. I = incisor; MO = molar.

3.4 Impact of Fgfr2c Mutation on Mandible

As demonstrated by EDMA of 14 anatomical landmarks (Figure 2B; Table 3), the mandible is generally smaller in the antero-posterior dimension of Fgfr2c^C342Y/+ embryos relative to Fgfr2c^+/+ littermates at E15.5 (Figure 7A). EDMA of 20 anatomical landmarks (Figure 2B; Table 3) indicated differences between genotypes are concentrated in the posterior of the developing mandible corresponding to an area that stretches from the posterior portion of Mid 1 through the rostral portion of Posterior of MC that includes the coronoid, condylar, and angular processes (Figure 7B). There was no overall difference in bone volume at E15.5 or E17.5 (Table 4).

FIGURE 7

FIGURE 7. Morphological differences in Fgfr2c^C342Y/+ and Fgfr2c^+/+ mandibular bone at E15.5 and E17.5. Results of EDMA using 3D landmark coordinates of 14 anatomical landmarks at E15.5 (A) and 20 anatomical landmarks at E17.5 (B) show linear distances within each model that are significantly smaller by at least 5% between genotypes. Lateral views (buccal) at left and superior views at right. Thin pink lines indicate linear distances that are decreased by 5%–10% in Fgfr2c^C342Y/+ embryos, and thick purple lines indicate linear distances that are reduced by 10%–23% in Fgfr2c^C342Y/+ embryos. For orientation, superimpositions of MC and mandibular bone are inset. Scalebar = 1.0 mm.

To assess differentiation of osteoblasts, which are responsible for morphological changes in the mandible of Fgfr2c^C342Y/+ embryos, we investigated the levels of ALP-positive osteoblasts, which are responsible for bone mineralization (Hall, 2014), but identified no significant effects of this mutation on the ratio of alkaline phosphatase positive osteoblasts to total cells in the mandible (Figures 8A,B). In contrast to the lack of effect on the number of bone building osteoblasts, bone resorbing osteoclasts (measured by TRAP staining and multi-nucleation) were reduced in Fgfr2c^C342Y/+ embryos at E17.5 relative to Fgfr2c^+/+ littermates. This overall reduction in osteoclasts was driven by a significant reduction in osteoclasts in Mid 1 of Fgfr2c^C342Y/+ embryos (Figures 8 C,D).

FIGURE 8

FIGURE 8. Immunohiostochemical analysis of osteoblasts and osteoclasts in developing mandible of Fgfr2c^C342Y/+ embryos and Fgfr2c^+/+ littermates. (A,C) Representative coronal sections of MC and surrounding mandible in Fgfr2c^+/+ (top rows) and Fgfr2c^C342Y/+ (bottom rows) embryos stained with alkaline phosphatase (ALP, (A) and Tartrate-Resistant Acid Phosphatase (TRAP, (C) at E15.5 (top) and E17.5 (bottom). Reduced TRAP staining can be noted in Fgfr2c^C342Y/+ images especially at Mid 1. (B) Quantification of ALP positive osteoblasts in the developing mandible standardized by the total number of cells within bone matrix. (D) Count of TRAP positive multi-nucleate osteoclasts in the mandible bone matrix. White bar in box plot identifies median. X in box plot identifies mean. n = 3 individuals per age per genotype with at least 4 images per region per embryo being analyzed. * = p ≤ 0.05.

3.5 Differential Effects of Fgfr2c C342Y Mutation on ERK Pathway Activation in Meckel’s Cartilage and Mandible

To assess how Fgfr2c C342Y alters the ERK pathway in the lower jaw, we analyzed the expression of p-ERK1/2 in MC and mandible at E15.5, the earliest timepoint when the phenotypic effects of Fgfr2c C342Y mutation on MC and mandible were observed, prior to significant degradation in MC (Sakakura et al., 2007). p-ERK1/2 was detected in the mandible, tongue, tooth, peripheral chondrocytes in MC and perichondrium of MC by immunostaining (Figure 9). Little phosphorylated ERK was identified within MC, however there was significant activation in the perichondrium particularly in Anterior, Mid 2, and Posterior of MC, though this activation was reduced in E15.5 Fgfr2c^C342Y/+ embryos as compared to Fgfr2c^+/+ littermates. For mandibular bone, the expression level and pattern of p-ERK1/2 in Anterior and Mid 1 of Fgfr2c^C342Y/+ embryos were similar to Fgfr2c^+/+ littermates, however, ERK activation is increased in Mid 2 and Posterior. Our assessment of the downstream effects of the Fgfr2c mutation indicates that the ERK pathway is activated particularly in the posterior portion of the mandible of Fgfr2c^C342Y/+ embryos at E15.5 (Figure 9), matching the area of greatest morphological differences at that age (Figure 7).

FIGURE 9

FIGURE 9. ERK pathway activation in the lower jaw of Fgfr2c^+/+ and Fgfr2c^C342Y/+ embryos at E15.5. Immunofluorescence for p-ERK1/2 in coronal sections of Fgfr2c^+/+ (top) and Fgfr2c^C342Y/+ (bottom) embryos at E15.5. Anterior, Mid 1, Mid 2, and Posterior indicate the regions of the lower jaw identified in Figure 1. Images in MB ROI column show higher magnification of the boxed areas in Mid 2. Arrows mark the locations of the perichondrium of MC. I, incisor; MB, mandible; MC, Meckel’s cartilage; T, tongue. n = 3 per genotype with at least 4 images per individual being analyzed. Scale = 100 µm.

4 Discussion

Though the morphological differences between MC and the developing mandible of Fgfr2c^C342Y/+ and Fgfr2c^+/+ embryos are subtle, our data demonstrate that this mutation has differential significant effects by tissue type and time at the cellular and morphological levels. We identified a larger, aberrantly formed MC in Fgfr2c^C342Y/+ embryos after E14 when mineralization of the mandible has already begun (Ramaesh and Bard, 2003). However, our data reveal a relatively small mandible in Fgfr2c^C342Y/+ embryos at E15.5 and E17.5. Similarly, investigations of newborn and postnatal Fgfr2c^C342Y/+ mandibles found them to be relatively small (Khominsky et al., 2018; Motch Perrine et al., 2019). These observations suggest a larger MC gives way to a smaller mandible during the latter days of prenatal development.

Overall, the significant differences in MC between genotypes identified by histomorphometry were intermittent in time and patchy in anatomic distribution. Our findings of increased cross-sectional area and matrix area in Fgfr2c^C342Y/+ E13.5 embryos were not associated with any significant changes in MC morphology at E13.5, but Fgfr2c^C342Y/+ MC becomes relatively larger later in development. These data suggest that even subtle changes early in development can produce significant morphological effects later, or that there are additional/alternate cellular processes at work.

Observations of growth plate cartilages in long bones identify chondrocyte proliferation, hypertrophy, and matrix deposition as the cellular processes that contribute to cartilage growth (Breur et al., 1991; Wilsman et al., 2008; Cooper et al., 2013). No differences between genotypes in MC area, contribution of matrix, chondrocyte number, or chondrocyte size were identified at E15.5, and proliferation was reduced in Fgfr2c^C342Y/+ embryos at E15.5. Few morphological and cellular differences between genotypes were identified at E17.5 as well. Though it is possible that statistically insignificant changes in cell size, cell number, and contribution of matrix can have biologically significant effects on development through additive and/or complementary mechanisms, it is also possible that MC increases in size using processes other than those identified in growth plate cartilages (Kaucka et al., 2017). Further, compensatory or redundant growth mechanisms for cartilage and bone might also explain why we did not identify dramatic effects of this mutation. The lack of consistent differences in MC morphology and cell behaviors between genotypes over time calls for alternate approaches focused on additional, potentially influential processes that produce affects that differ by cell type and developmental stage.

Investigation of the downstream targets of FGFR2 activation confirmed our assessment that the effects of the Fgfr2c C342Y mutation are cell type specific. ERK1/2 is the main effector of FGF signaling (Brewer et al., 2016; Ray et al., 2020). ERK pathway promotes osteoblast differentiation and bone formation in vitro and in vivo (Xiao et al., 2000; Xiao et al., 2002; Kim et al., 2019). In a previous study, Fgfr2c^C342Y/+ mice exhibited upregulated p-ERK in osteogenic fronts of the coronal suture compared with controls (Shibata et al., 2014; Lee et al., 2018). Our results show that ERK phosphorylation is increased in mandibular bone related cells and decreased in cartilage related perichondrium in the Fgfr2c^C342Y/+ embryos. From these assessments, we can conclude that the effects of constitutive activation of the FGF pathway are cell type specific.

The assertion that MC functions as a model for developing bony tissues of the mandible early in development and as a template for growth later in development is not corroborated by our analysis. Not only do the embryonic Fgfr2c^C342Y/+ mandible and MC differ in opposite ways relative to their unaffected littermates, but we see few consistent changes in cellular activity in the four lower jaw regions identified across the embryonic times assayed. Likewise, in mice in which the Sox9 gene, a transcription factor essential for endochondral bone formation, was conditionally inactivated in cranial neural crest cells, cartilage derived from these cells and endochondral bones were totally absent. However, the mandible still formed in these mice though it was small and misshapen (Mori-Akiyama et al., 2003). Though we know MC lengthens and increases in diameter with growth, the cell processes responsible for these and other morphological changes in MC have been studied by only a few investigators (e.g., Sakakura et al., 2007; Shibata et al., 2014; Shibata et al., 2014; Kaucka et al., 2017) and are not well defined.

Maintaining the balance between proliferation and differentiation is crucial for osteogenesis. Alkaline phosphatase, which identifies osteoprogenitors, pre-osteoblasts and differentiated osteoblasts (Huang et al., 2007) can be detected in the mandible as early as E11 (Shibata and Yokohama-Tamaki, 2008). Though no significant differences between genotypes were identified by alkaline phosphatase staining, reduced proliferation and increased ERK activity in the mandible indicate a shift from proliferation to osteoblast differentiation in the mandible, particularly posteriorly at E15.5 (Figure 9). Interestingly, this result correlated with the location of the most significant morphological differences consistent across the age groups revealed by EDMA (Figure 7). Our results suggest that subtle tissue-specific imbalance between proliferation and differentiation triggered by the Fgfr2c C342Y mutation may contribute to the small morphological changes of MC and mandible. In addition, the tissue, time, and region-specific effects of this mutation can impact the relationship between developing MC and mandible. Our results prompt the hypothesis that seemingly small effects combine additively or by interaction to subtly change MC and mandible morphology. Further, Modifying the relationship between these important structures may produce aberrant scaffolding that contributes to dysmorphogenesis associated with this mutation.

Data Availability Statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation. Bone microCT and PTA-enhanced microCT data used for this study are also available on Penn State University Libraries ScholarSphere; doi:10.26207/qgke-r185.

Ethics Statement

The animal study was reviewed and approved by the Icahn School of Medicine at Mount Sinai and Pennsylvania State University Animal Care and Use Committees (PSU #46558 and ISMMS #07-0757).

Author Contributions

MP: Conceptualization, methodology, validation, formal analysis, investigation, data curation, writing-original draft, review and editing, visualization. ED: Conceptualization, methodology, validation, formal analysis, investigation, data curation, writing-original draft, review and editing, visualization. NR: Formal analysis, investigation, review and editing. JS: Investigation, review and editing. AC: Investigation, review and editing. HZ: Software, data curation, writing—review and editing. DC: Software, data curation, writing—review and editing. KK: Conceptualization, visualization, writing—review and editing. EJ: Resources, supervision, funding acquisition, writing—review and editing. JR: Resources, supervision, project administration, funding acquisition, writing—review and editing. MW: Conceptualization, methodology, validation, formal analysis, investigation, data curation, writing-original draft, review and editing, visualization, project administration SP: Conceptualization, methodology, validation, formal analysis, investigation, data curation, writing-original draft, review and editing, visualization, project administration.

Funding

This work was supported in part by NIH/NICDR grants R01 DE027677 to JR; P01 HD078233 to JR and EJ; R01 DE029832 and R01 DE031439 to JR and EJ with support for SP and MW; R01 DE029322 to EJ and MW; NSF grant CCF-1617735 to DC.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.871927/full#supplementary-material

References

Amin, S., and Tucker, A. S. (2006). Joint Formation in the Middle Ear: Lessons from the Mouse and guinea Pig. Dev. Dyn. 235, 1326–1333. doi:10.1002/dvdy.20666

PubMed Abstract | CrossRef Full Text | Google Scholar

Azoury, S. C., Reddy, S., Shukla, V., and Deng, C.-X. (2017). Fibroblast Growth Factor Receptor 2 (FGFR2) Mutation Related Syndromic Craniosynostosis. Int. J. Biol. Sci. 13, 1479–1488. doi:10.7150/ijbs.22373

PubMed Abstract | CrossRef Full Text | Google Scholar

Baddoo, M., Hill, K., Wilkinson, R., Gaupp, D., Hughes, C., Kopen, G. C., et al. (2003). Characterization of Mesenchymal Stem Cells Isolated from Murine Bone Marrow by Negative Selection. J. Cel. Biochem. 89, 1235–1249. doi:10.1002/jcb.10594

CrossRef Full Text | Google Scholar

Behringer, R., Gertsenstein, M., Nagy, K., and Nagy, A. (2013). Manipulating the Mouse Embryo: A Laboratory Manual. Fourth edition (Cold Spring Harbor, New York: Cold Spring Harbor Laboratory Press). Revised edition.

Google Scholar

Boutros, S., Shetye, P. R., Ghali, S., Carter, C. R., McCarthy, J. G., and Grayson, B. H. (2007). Morphology and Growth of the Mandible in Crouzon, Apert, and Pfeiffer Syndromes. J. Craniofac. Surg. 18, 146–150. doi:10.1097/01.scs.0000248655.53405.a7

PubMed Abstract | CrossRef Full Text | Google Scholar

Breur, G. J., VanEnkevort, B. A., Farnum, C. E., and Wilsman, N. J. (1991). Linear Relationship between the Volume of Hypertrophic Chondrocytes and the Rate of Longitudinal Bone Growth in Growth Plates. J. Orthop. Res. 9, 348–359. doi:10.1002/jor.1100090306

PubMed Abstract | CrossRef Full Text | Google Scholar

Brewer, J. R., Mazot, P., and Soriano, P. (2016). Genetic Insights into the Mechanisms of Fgf Signaling. Genes Dev. 30, 751–771. doi:10.1101/gad.277137.115

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, P., Zhang, L., Weng, T., Zhang, S., Sun, S., Chang, M., et al. (2014). A Ser252Trp Mutation in Fibroblast Growth Factor Receptor 2 (FGFR2) Mimicking Human Apert Syndrome Reveals an Essential Role for FGF Signaling in the Regulation of Endochondral Bone Formation. PLoS One 9, e87311. doi:10.1371/journal.pone.0087311

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, M. M., and Krelborg, S. (1992). Birth Prevalence Studies of the Crouzon Syndrome: Comparison of Direct and Indirect Methods. Clin. Genet. 41, 12–15. doi:10.1111/j.1399-0004.1992.tb03620.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cole, T. (2002). WinEDMA: Software for Euclidean Distance Matrix Analysis. Kansas City: University of Missouri, School of Medicine.

Google Scholar

Cooper, K. L., Oh, S., Sung, Y., Dasari, R. R., Kirschner, M. W., and Tabin, C. J. (2013). Multiple Phases of Chondrocyte Enlargement Underlie Differences in Skeletal Proportions. Nature 495, 375–378. doi:10.1038/nature11940

PubMed Abstract | CrossRef Full Text | Google Scholar

Costaras-Volarich, M., and Pruzansky, S. (1984). Is the Mandible Intrinsically Different in Apert and Crouzon Syndromes? Am. J. Orthod. 85, 475–487. doi:10.1016/0002-9416(84)90087-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Deckelbaum, R. A., Holmes, G., Zhao, Z., Tong, C., Basilico, C., and Loomis, C. A. (2012). Regulation of Cranial Morphogenesis and Cell Fate at the Neural Crest-Mesoderm Boundary by Engrailed 1. Development 139, 1346–1358. doi:10.1242/dev.076729

PubMed Abstract | CrossRef Full Text | Google Scholar

Durham, E. L., Howie, R. N., Houck, R., Oakes, B., Grey, Z., Hall, S., et al. (2018). Involvement of Calvarial Stem Cells in Healing: A Regional Analysis of Large Cranial Defects. Wound Rep. Reg. 26, 359–365. doi:10.1111/wrr.12658

PubMed Abstract | CrossRef Full Text | Google Scholar

Durham, E., Howie, R. N., Larson, N., LaRue, A., and Cray, J. (2019). Pharmacological Exposures May Precipitate Craniosynostosis through Targeted Stem Cell Depletion. Stem Cel Res. 40, 101528. doi:10.1016/j.scr.2019.101528

PubMed Abstract | CrossRef Full Text | Google Scholar

Eswarakumar, V. P., Horowitz, M. C., Locklin, R., Morriss-Kay, G. M., and Lonai, P. (2004). A Gain-Of-Function Mutation of Fgfr2c Demonstrates the Roles of This Receptor Variant in Osteogenesis. Proc. Natl. Acad. Sci. U.S.A. 101, 12555–12560. doi:10.1073/pnas.0405031101

PubMed Abstract | CrossRef Full Text | Google Scholar

Frisdal, A., and Trainor, P. A. (2014). Development and Evolution of the Pharyngeal Apparatus. Wires Dev. Biol. 3, 403–418. doi:10.1002/wdev.147

PubMed Abstract | CrossRef Full Text | Google Scholar

Graham, A., and Richardson, J. (2012). Developmental and Evolutionary Origins of the Pharyngeal Apparatus. EvoDevo 3, 24. doi:10.1186/2041-9139-3-24

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, B. K. (2014). “Bones and Cartilage,” in Developmental and Evolutionary Skeletal Biology. Second Edition (Amsterdam: Elsevier/AP).

Google Scholar

Haub, P., and Meckel, T. (2015). A Model Based Survey of Colour Deconvolution in Diagnostic Brightfield Microscopy: Error Estimation and Spectral Consideration. Sci. Rep. 5, 12096. doi:10.1038/srep12096

PubMed Abstract | CrossRef Full Text | Google Scholar

Heuzé, Y., Holmes, G., Peter, I., Richtsmeier, J. T., and Jabs, E. W. (2014). Closing the gap: Genetic and Genomic Continuum from Syndromic to Nonsyndromic Craniosynostoses. Curr. Genet. Med. Rep. 2, 135–145. doi:10.1007/s40142-014-0042-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Holmes, G., O'Rourke, C., Perrine, S. M. M., Lu, N., van Bakel, H., Richtsmeier, J. T., et al. (2018). Midface and Upper Airway Dysgenesis in FGFR2-Craniosynostosis Involves Multiple Tissue-specific and Cell Cycle Effects. Development 145, dev166488. doi:10.1242/dev.166488

PubMed Abstract | CrossRef Full Text | Google Scholar

Howie, R. N., Herberg, S., Durham, E., Grey, Z., Bennfors, G., Elsalanty, M., et al. (2018). Selective Serotonin Re-uptake Inhibitor Sertraline Inhibits Bone Healing in a Calvarial Defect Model. Int. J. Oral Sci. 10, 25. doi:10.1038/s41368-018-0026-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, W., Yang, S., Shao, J., and Li, Y.-P. (2007). Signaling and Transcriptional Regulation in Osteoblast Commitment and Differentiation. Front. Biosci. 12, 3068–3092. doi:10.2741/2296

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishizeki, K., Saito, H., Shinagawa, T., Fujiwara, N., and Nawa, T. (1999). Histochemical and Immunohistochemical Analysis of the Mechanism of Calcification of Meckel's Cartilage during Mandible Development in Rodents. J. Anat. 194 (Pt 2), 265–277. doi:10.1046/j.1469-7580.1999.19420265.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishizeki, K., Takahashi, N., and Nawa, T. (2001). Formation of the Sphenomandibular Ligament by Meckel's Cartilage in the Mouse: Possible Involvement of Epidermal Growth Factor as Revealed by Studies In Vivo and In Vitro. Cel Tissue Res. 304, 67–80. doi:10.1007/s004410100354

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, D., and Wilkie, A. O. M. (2011). Craniosynostosis. Eur. J. Hum. Genet. 19, 369–376. doi:10.1038/ejhg.2010.235

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaucka, M., Zikmund, T., Tesarova, M., Gyllborg, D., Hellander, A., Jaros, J., et al. (2017). Oriented Clonal Cell Dynamics Enables Accurate Growth and Shaping of Vertebrate Cartilage. eLife 6, e25902. doi:10.7554/eLife.25902

PubMed Abstract | CrossRef Full Text | Google Scholar

Khominsky, A., Yong, R., Ranjitkar, S., Townsend, G., and Anderson, P. J. (2018). Extensive Phenotyping of the Orofacial and Dental Complex in Crouzon Syndrome. Arch. Oral Biol. 86, 123–130. doi:10.1016/j.archoralbio.2017.10.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J.-M., Yang, Y.-S., Park, K. H., Oh, H., Greenblatt, M. B., and Shim, J.-H. (2019). The ERK MAPK Pathway Is Essential for Skeletal Development and Homeostasis. Ijms 20, 1803. doi:10.3390/ijms20081803

PubMed Abstract | CrossRef Full Text | Google Scholar

Landini, G., Martinelli, G., and Piccinini, F. (2021). Colour Deconvolution: Stain Unmixing in Histological Imaging. Bioinformatics 37, 1485–1487. doi:10.1093/bioinformatics/btaa847

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, K. K. L., Peskett, E., Quinn, C. M., Aiello, R., Adeeva, L., Moulding, D. A., et al. (2018). Overexpression of Fgfr2c Causes Craniofacial Bone Hypoplasia and Ameliorates Craniosynostosis in the Crouzon Mouse. Dis. Model. Mech. 11, dmm035311. doi:10.1242/dmm.035311

PubMed Abstract | CrossRef Full Text | Google Scholar

Lele, S., and Richtsmeier, J. T. (1995). Euclidean Distance Matrix Analysis: Confidence Intervals for Form and Growth Differences. Am. J. Phys. Anthropol. 98, 73–86. doi:10.1002/AJPA.1330980107

PubMed Abstract | CrossRef Full Text | Google Scholar

Lele, S. R., and Richtsmeier, J. T. (2001). An Invariant Approach to Statistical Analysis of Shapes. New York: Chapman and Hall/CRC. doi:10.1201/9781420036176

CrossRef Full Text | Google Scholar

Lesciotto, K. M., Motch Perrine, S. M., Kawasaki, M., Stecko, T., Ryan, T. M., Kawasaki, K., et al. (2020). Phosphotungstic Acid‐enhanced microCT: Optimized Protocols for Embryonic and Early Postnatal Mice. Dev. Dyn. 249, 573–585. doi:10.1002/dvdy.136

PubMed Abstract | CrossRef Full Text | Google Scholar

Lesciotto, K. M., Tomlinson, L., Leonard, S., and Richtsmeier, J. T. (2022). Embryonic and Early Postnatal Cranial Bone Volume and Tissue mineral Density Values for C57BL/6J Laboratory Mice. Dev. Dyn., 1–13. doi:10.1002/dvdy.458

CrossRef Full Text | Google Scholar

Lu, X., Sawh-Martinez, R., Forte, A. J., Wu, R., Cabrejo, R., Wilson, A., et al. (2019). Mandibular Spatial Reorientation and Morphological Alteration of Crouzon and Apert Syndrome. Ann. Plast. Surg. 83, 568–582. doi:10.1097/SAP.0000000000001811

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, X., Sawh-Martinez, R., Forte, A. J., Wu, R., Cabrejo, R., Wilson, A., et al. (2020). Classification of Subtypes of Crouzon Syndrome Based on the Type of Vault Suture Synostosis. J. Craniofac. Surg. 31, 678–684. doi:10.1097/SCS.0000000000006173

PubMed Abstract | CrossRef Full Text | Google Scholar

Manocha, S., Farokhnia, N., Khosropanah, S., Bertol, J. W., Santiago, J., and Fakhouri, W. D. (2019). Systematic Review of Hormonal and Genetic Factors Involved in the Nonsyndromic Disorders of the Lower Jaw. Dev. Dyn. 248, 162–172. doi:10.1002/dvdy.8

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez-Abadías, N., Percival, C., Aldridge, K., Hill, C. A., Ryan, T., Sirivunnabood, S., et al. (2010). Beyond the Closed Suture in Apert Syndrome Mouse Models: Evidence of Primary Effects of FGFR2 Signaling on Facial Shape at Birth. Dev. Dyn. 239, 3058–3071. doi:10.1002/dvdy.22414

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez-Abadías, N., Motch, S. M., Pankratz, T. L., Wang, Y., Aldridge, K., Jabs, E. W., et al. (2013). Tissue-specific Responses to Aberrant FGF Signaling in Complex Head Phenotypes. Dev. Dyn. 242, 80–94. doi:10.1002/dvdy.23903

PubMed Abstract | CrossRef Full Text | Google Scholar

Maruyama, T., Mirando, A. J., Deng, C.-X., and Hsu, W. (2010). The Balance of WNT and FGF Signaling Influences Mesenchymal Stem Cell Fate during Skeletal Development. Sci. Signal. 3, ra40. doi:10.1126/scisignal.2000727

PubMed Abstract | CrossRef Full Text | Google Scholar

McKenzie, J. (1958). The First Arch Syndrome. Arch. Dis. Child. 33, 477–486. doi:10.1136/adc.33.171.477

PubMed Abstract | CrossRef Full Text | Google Scholar

Miraoui, H., Oudina, K., Petite, H., Tanimoto, Y., Moriyama, K., and Marie, P. J. (2009). Fibroblast Growth Factor Receptor 2 Promotes Osteogenic Differentiation in Mesenchymal Cells via ERK1/2 and Protein Kinase C Signaling. J. Biol. Chem. 284, 4897–4904. doi:10.1074/jbc.M805432200

PubMed Abstract | CrossRef Full Text | Google Scholar

Mori-Akiyama, Y., Akiyama, H., Rowitch, D. H., and de Crombrugghe, B. (2003). Sox9 Is Required for Determination of the Chondrogenic Cell Lineage in the Cranial Neural Crest. Proc. Natl. Acad. Sci. U.S.A. 100, 9360–9365. doi:10.1073/pnas.1631288100

PubMed Abstract | CrossRef Full Text | Google Scholar

Motch Perrine, S. M., Stecko, T., Neuberger, T., Jabs, E. W., Ryan, T. M., and Richtsmeier, J. T. (2017). Integration of Brain and Skull in Prenatal Mouse Models of Apert and Crouzon Syndromes. Front. Hum. Neurosci. 11, 369. doi:10.3389/fnhum.2017.00369

PubMed Abstract | CrossRef Full Text | Google Scholar

Motch Perrine, S. M., Wu, M., Stephens, N. B., Kriti, D., van Bakel, H., Jabs, E. W., et al. (2019). Mandibular Dysmorphology Due to Abnormal Embryonic Osteogenesis in FGFR2-Related Craniosynostosis Mice. Dis. Models Mech. 12, dmm038513. doi:10.1242/dmm.038513

CrossRef Full Text | Google Scholar

Motch Perrine, S. M., Pitirri, M. K., Durham, E. L., Kawasaki, M., Zheng, H., Chen, D. Z., et al. (2021). Over 400 Million Years of Cooperation: Untangling the Chondrocranium-Dermatocranium Connection. bioRxiv 2021.11.24.469914. doi:10.1101/2021.11.24.469914

CrossRef Full Text | Google Scholar

Musy, M., Flaherty, K., Raspopovic, J., Robert-Moreno, A., Richtsmeier, J., and Sharpe, J. (2018). A Quantitative Method for Staging Mouse Embryos Based on Limb Morphometry. Development 145, dev154856. doi:10.1242/dev.154856

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakamura, M., Yang, M.-C., Ashida, K., Mayanagi, M., and Sasano, Y. (2021). Calcification and resorption of mouse Meckel's cartilage analyzed by von Kossa and tartrate-resistant acid phosphatase histochemistry and scanning electron microscopy/energy-dispersive X-ray spectrometry. Anat. Sci. Int. 97, 213–220. doi:10.1007/s12565-021-00643-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ornitz, D. M., and Itoh, N. (2015). The Fibroblast Growth Factor Signaling Pathway. WIREs Dev. Biol. 4, 215–266. doi:10.1002/wdev.176

PubMed Abstract | CrossRef Full Text | Google Scholar

Parada, C., Han, D., Grimaldi, A., Sarrión, P., Park, S. S., Pelikan, R., et al. (2015). Disruption of the ERK/MAPK Pathway in Neural Crest Cells as a Potential Cause of Pierre Robin Sequence. Development 142, 3734–3745. doi:10.1242/dev.125328

PubMed Abstract | CrossRef Full Text | Google Scholar

Perlyn, C. A., Morriss-Kay, G., Darvann, T., Tenenbaum, M., and Ornitz, D. M. (2006). Model for the Pharmacologic Treatment of Crouzon Syndrome. Neurosurgery 59, 210–215. doi:10.1227/01.NEU.0000224323.53866.1E

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramaesh, T., and Bard, J. B. L. (2003). The Growth and Morphogenesis of the Early Mouse Mandible: A Quantitative Analysis. J. Anat. 203, 213–222. doi:10.1046/j.1469-7580.2003.00210.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ray, A. T., Mazot, P., Brewer, J. R., Catela, C., Dinsmore, C. J., and Soriano, P. (2020). FGF Signaling Regulates Development by Processes beyond Canonical Pathways. Genes Dev. 34, 1735–1752. doi:10.1101/gad.342956.120

PubMed Abstract | CrossRef Full Text | Google Scholar

Ronneberger, O., Fischer, P., and Brox, T. (2015). “U-net: Convolutional Networks for Biomedical Image Segmentation,” in Medical Image Computing and Computer-Assisted Intervention – MICCAI 2015. Editors N. Navab, J. Hornegger, W. M. Wells, and A. F. Frangi (Cham: Springer International Publishing), 234–241. doi:10.1007/978-3-319-24574-4_28

CrossRef Full Text | Google Scholar

Ruifrok, A. C., and Johnston, D. A. (2001). Quantification of Histochemical Staining by Color Deconvolution. Anal. Quant. Cytol. Histol. 23 (4), 291–299.

PubMed Abstract | Google Scholar

Sakakura, Y., Hosokawa, Y., Tsuruga, E., Irie, K., and Yajima, T. (2007). In Situ localization of Gelatinolytic Activity during Development and Resorption of Meckel's Cartilage in Mice. Eur. J. Oral Sci. 115, 212–223. doi:10.1111/j.1600-0722.2007.00447.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shibata, S., and Yokohama-Tamaki, T. (2008). Anin Situhybridization Study ofRunx2,Osterix, andSox9in the Anlagen of Mouse Mandibular Condylar Cartilage in the Early Stages of Embryogenesis. J. Anat. 213, 274–283. doi:10.1111/j.1469-7580.2008.00934.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shibata, S., Sakamoto, Y., Yokohama-Tamaki, T., Murakami, G., and Cho, B. H. (2014). Distribution of Matrix Proteins in Perichondrium and Periosteum during the Incorporation of Meckel's Cartilage into Ossifying Mandible in Midterm Human Fetuses: An Immunohistochemical Study. Anat. Rec. 297, 1208–1217. doi:10.1002/ar.22911

PubMed Abstract | CrossRef Full Text | Google Scholar

Shone, V., Oulion, S., Casane, D., Laurenti, P., and Graham, A. (2016). Mode of Reduction in the Number of Pharyngeal Segments within the Sarcopterygians. Zoolog. Lett 2, 6. doi:10.1186/s40851-016-0043-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Svandova, E., Anthwal, N., Tucker, A. S., and Matalova, E. (2020). Diverse Fate of an Enigmatic Structure: 200 Years of Meckel's Cartilage. Front. Cel Dev. Biol. 8, 821. doi:10.3389/fcell.2020.00821

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilsman, N. J., Bernardini, E. S., Leiferman, E., Noonan, K., and Farnum, C. E. (2008). Age and Pattern of the Onset of Differential Growth Among Growth Plates in Rats. J. Orthop. Res. 26, 1457–1465. doi:10.1002/jor.20547

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, G., Jiang, D., Thomas, P., Benson, M. D., Guan, K., Karsenty, G., et al. (2000). MAPK Pathways Activate and Phosphorylate the Osteoblast-specific Transcription Factor, Cbfa1. J. Biol. Chem. 275, 4453–4459. doi:10.1074/jbc.275.6.4453

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, G., Jiang, D., Gopalakrishnan, R., and Franceschi, R. T. (2002). Fibroblast Growth Factor 2 Induction of the Osteocalcin Gene Requires MAPK Activity and Phosphorylation of the Osteoblast Transcription Factor, Cbfa1/Runx2. J. Biol. Chem. 277, 36181–36187. doi:10.1074/jbc.M206057200

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, W., Luo, F., Wang, Q., Tan, Q., Huang, J., Zhou, S., et al. (2017). Inducible Activation of FGFR2 in Adult Mice Promotes Bone Formation after Bone Marrow Ablation. J. Bone Miner Res. 32, 2194–2206. doi:10.1002/jbmr.3204

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeh, E., Atique, R., Ishiy, F. A. A., Fanganiello, R. D., Alonso, N., Matushita, H., et al. (2012). FGFR2 Mutation Confers a Less Drastic Gain of Function in Mesenchymal Stem Cells Than in Fibroblasts. Stem Cel Rev Rep 8, 685–695. doi:10.1007/s12015-011-9327-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, K., Xu, J., Liu, Z., Sosic, D., Shao, J., Olson, E. N., et al. (2003). Conditional Inactivation of FGF Receptor 2 Reveals an Essential Role for FGF Signaling in the Regulation of Osteoblast Function and Bone Growth. Development 130, 3063–3074. doi:10.1242/dev.00491

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, H., Motch Perrine, S. M., Pitirri, M. K., Kawasaki, K., Wang, C., Richtsmeier, J. T., et al. (2020). Cartilage Segmentation in High-Resolution 3D Micro-CT Images via Uncertainty-Guided Self-Training with Very Sparse Annotation. Med. Image Comput. Comput. Assist. Interv. 12261, 802–812. doi:10.1007/978-3-030-59710-8_78

PubMed Abstract | CrossRef Full Text | Google Scholar