Search

Asymmetries in craniofacial anomalies are commonly observed. In the facial skeleton, the left side is more commonly and/or severely affected than the right. Such asymmetries complicate treatment options. Mechanisms underlying variation in disease severity between individuals as well as within individuals (asymmetries) are still relatively unknown.Developmental reductions in fibroblast growth factor 8 (Fgf8) have a dosage dependent effect on jaw size, shape, and symmetry. Further, Fgf8 mutants have directionally asymmetric jaws with the left side being more affected than the right. Defects in lower jaw development begin with disruption to Meckel's cartilage, which is discontinuous. All skeletal elements associated with the proximal condensation are dysmorphic, exemplified by a malformed and misoriented malleus. At later stages, Fgf8 mutants exhibit syngnathia, which falls into two broad categories: bony fusion of the maxillary and mandibular alveolar ridges and zygomatico-mandibular fusion. All of these morphological defects exhibit both inter- and intra-specimen variation.We hypothesize that these asymmetries are linked to heart development resulting in higher levels of Fgf8 on the right side of the face, which may buffer the right side to developmental perturbations. This mouse model may facilitate future investigations of mechanisms underlying human syngnathia and facial asymmetry.

Morphogenesis requires highly coordinated, complex interactions between cellular processes: proliferation, migration, and apoptosis, along with physical tissue interactions. How these cellular and tissue dynamics drive morphogenesis remains elusive. Three dimensional (3D) microscopic imaging poses great promise, and generates elegant images. However, generating even moderate through-put quantified images is challenging for many reasons. As a result, the association between morphogenesis and cellular processes in 3D developing tissues has not been fully explored. To address this critical gap, we have developed an imaging and image analysis pipeline to enable 3D quantification of cellular dynamics along with 3D morphology for the same individual embryo. Specifically, we focus on how 3D distribution of proliferation relates to morphogenesis during mouse facial development. Our method involves imaging with light-sheet microscopy, automated segmentation of cells and tissues using machine learning-based tools, and quantification of external morphology via geometric morphometrics. Applying this framework, we show that changes in proliferation are tightly correlated to changes in morphology over the course of facial morphogenesis. These analyses illustrate the potential of this pipeline to investigate mechanistic relationships between cellular dynamics and morphogenesis during embryonic development.

Background Postoperative cognitive decline (PCD) requires microglial activation. Voltage-gated Kv1.3 potassium channels are involved in microglial activation. We determined the role of Kv1.3 in PCD and the efficacy and safety of inhibiting Kv1.3 with phenoxyalkoxypsoralen-1 (PAP-1) in preventing PCD in a mouse model. Methods After institutional approval, we assessed whether Kv1.3-deficient mice (Kv1.3–/–) exhibited PCD, evidenced by tibial-fracture surgery-induced decline in aversive freezing behaviour, and whether PAP-1 could prevent PCD and postoperative neuroinflammation in PCD-vulnerable diet-induced obese (DIO) mice. We also evaluated whether PAP-1 altered either postoperative peripheral inflammation or tibial-fracture healing. Results Freezing behaviour was unaltered in postoperative Kv1.3–/– mice. In DIO mice, PAP-1 prevented postoperative (i) attenuation of freezing behaviour (54 [17.3]% vs 33.4 [12.7]%; P=0.03), (ii) hippocampal microglial activation by size (130 [31] pixels vs 249 [49]; P Conclusions Microglial-mediated PCD requires Kv1.3 activity, determined by genetic and pharmacological targeting approaches. Phenoxyalkoxypsoralen-1 blockade of Kv1.3 prevented surgery-induced hippocampal microglial activation and neuroinflammation in mice known to be vulnerable to PCD. Regarding perioperative safety, these beneficial effects of PAP-1 treatment occurred without impacting fracture healing. Kv1.3 blockers, currently undergoing clinical trials for other conditions, may represent an effective and safe intervention to prevent PCD.

Background The frontonasal ectodermal zone (FEZ) is a signaling center that regulates patterned development of the upper jaw, and Sonic hedgehog (SHH) mediates FEZ activity. Induction of SHH expression in the FEZ results from SHH-dependent signals from the brain and neural crest cells. Given the role of miRNAs in modulating gene expression, we investigated the extent to which miRNAs regulate SHH expression and FEZ signaling. Results In the FEZ, the miR-199 family appears to be regulated by SHH-dependent signals from the brain; expression of this family increased from HH18 to HH22, and upon activation of SHH signaling in the brain. However, the miR-199 family is more broadly expressed in the mesenchyme of the frontonasal process and adjacent neuroepithelium. Downregulating the miR-199 genes expanded SHH expression in the FEZ, resulting in wider faces, while upregulating miR-199 genes resulted in decreased SHH expression and narrow faces. Hypoxia inducible factor 1 alpha (HIF1A) and mitogen-activated protein kinase kinase kinase 4 (MAP3K4) appear to be potential targets of miR-199b. Reduction of MAP3K4 altered beak development but increased apoptosis, while reducing HIF1A reduced expression of SHH in the FEZ and produced malformations independent of apoptosis. Conclusions Our results demonstrate that this miRNA family appears to participate in regulating SHH expression in the FEZ; however, specific molecular mechanisms remain unknown.

The majority of fractures, especially in elderly and osteoporotic patients, occurs in metaphyseal bone. However, only a few experimental models exist to study metaphyseal bone healing in mice. Currently used mouse models of metaphyseal fracture healing are either based on drill hole defects, lacking adequate biomechanical stimulation at the site of fracture and therefore endochondral ossification in the fracture callus, or are introduced into the distal part of the mouse femur stabilized by a locking plate, which is challenging due to the small specimen size. Therefore, the aim of the current study was to develop a new mouse model to study metaphyseal fracture healing of the proximal femur. We chose a combination between an open osteotomy and a closed intramedullary stabilization. A 24 G needle was inserted into the femur in a closed manner, then an osteotomy was made with a 0.4-mm Gigli wire saw between the third and the lesser trochanter of the femur using an open approach. Fractured femurs were analyzed using microcomputed tomography and histology at days 14 and 21 after surgery. No animals were lost due to surgery or anesthesia. All animals displayed normal limb loading and a physiological gait pattern within the first three days after fracture. We found robust endochondral ossification during the fracture healing process with high expression of late chondrocyte and early osteogenic markers at day 14 (d14). By day 21 (d21), all fractures had a bony bridging score of 3 or more, indicating successful healing. Callus volume significantly decreased from d14 to d21, whereas high numbers of osteoclasts appeared at the fracture callus until d21, indicating that callus remodeling had already started at d21. In conclusion, we successfully developed a novel mouse model to study endochondral fracture healing of the proximal femur. This model might be useful for future studies using transgenic animals to unravel molecular mechanisms of osteoporotic metaphyseal fracture healing.

Complex morphological traits are the product of many genes with transient or lasting developmental effects that interact in anatomical context. Mouse models are a key resource for disentangling such effects, because they offer myriad tools for manipulating the genome in a controlled environment. Unfortunately, phenotypic data are often obtained using laboratory-specific protocols, resulting in self-contained datasets that are difficult to relate to one another for larger scale analyses. To enable meta-analyses of morphological variation, particularly in the craniofacial complex and brain, we created MusMorph, a database of standardized mouse morphology data spanning numerous genotypes and developmental stages, including E10.5, E11.5, E14.5, E15.5, E18.5, and adulthood. To standardize data collection, we implemented an atlas-based phenotyping pipeline that combines techniques from image registration, deep learning, and morphometrics. Alongside stage-specific atlases, we provide aligned micro-computed tomography images, dense anatomical landmarks, and segmentations (if available) for each specimen (N=10,056). Our workflow is open-source to encourage transparency and reproducible data collection. The MusMorph data and scripts are available on FaceBase (www.facebase.org, doi.org/10.25550/3-HXMC) and GitHub (https://github.com/jaydevine/MusMorph).

Fractures in the elderly represent a significant and rising socioeconomic problem. Although aging has been associated with delays in healing, there is little direct clinical data isolating the effects of aging on bone healing from the associated comorbidities that are frequently present in elderly populations. Basic research has demonstrated that all of the components of fracture repair -cells, extracellular matrix, blood supply, and molecules and their receptors— are negatively impacted by the aging process, which likely explains poorer clinical outcomes. Improved understanding of age-related fracture healing should aid in the development of novel treatment strategies, technologies, and therapies to improve bone repair in elderly patients.

Background and ObjectiveA variety of genetic mutations are known to affect cell proliferation and apoptosis during organism development, leading to structural birth defects such as facial clefting. Yet, the mechanisms how these alterations influence the development of the face remain unclear. Cell proliferation and its relation to shape variation can be studied in high detail using Light-Sheet Microscopy (LSM) imaging across a range of developmental time points. However, the large number of LSM images captured at cellular resolution precludes manual analysis. Thus, the aim of this work was to develop and evaluate automatic methods to segment tissues and proliferating cells in these images in an accurate and efficient way.MethodsWe developed, trained, and evaluated convolutional neural networks (CNNs) for segmenting tissues, cells, and specifically proliferating cells in LSM datasets. We compared the automatically extracted tissue and cell annotations to corresponding manual segmentations for three specific applications: (i) tissue segmentation (neural ectoderm and mesenchyme) in nuclear-stained LSM images, (ii) cell segmentation in nuclear-stained LSM images, and (iii) segmentation of proliferating cells in Phospho-Histone H3 (PHH3)-stained LSM images.ResultsThe automatic CNN-based tissue segmentation method achieved a macro-average F-score of 0.84 compared to a macro-average F-score of 0.89 comparing corresponding manual segmentations from two observers. The automatic cell segmentation method in nuclear-stained LSM images achieved an F-score of 0.57, while comparing the manual segmentations resulted in an F-score of 0.39. Finally, the automatic segmentation method of proliferating cells in the PHH3-stained LSM datasets achieved an F-score of 0.56 for the automated method, while comparing the manual segmentations resulted in an F-score of 0.45.ConclusionsThe proposed automatic CNN-based framework for tissue and cell segmentation leads to results comparable to the inter-observer agreement, accelerating the LSM image analysis. The trained CNN models can also be applied for shape or morphological analysis of embryos, and more generally in other areas of cell biology.

The avian embryo has been used as a model system for more than a century and has led to fundamental understanding of vertebrate development. One of the strengths of this model system is that the effect of, and interaction among, tissues can be directly assessed in chimeric embryos. We have previously shown that signals from the forebrain contribute to facial morphogenesis by regulating the shape of the expression domain of Sonic hedgehog (SHH) in the Frontonasal Ectodermal Zone (FEZ). In this article, the method of generating the forebrain chimeras and provide illustrations of the outcomes of these experiments is described.

Geometric morphometrics is the statistical analysis of landmark-based shape variation and its covariation with other variables. Over the past two decades, the gold standard of landmark data acquisition has been manual detection by a single observer. This approach has proven accurate and reliable in small-scale investigations. However, big data initiatives are increasingly common in biology and morphometrics. This requires fast, automated, and standardized data collection. We combine techniques from image registration, geometric morphometrics, and deep learning to automate and optimize anatomical landmark detection. We test our method on high-resolution, micro-computed tomography images of adult mouse skulls. To ensure generalizability, we use a morphologically diverse sample and implement fundamentally different deformable registration algorithms. Compared to landmarks derived from conventional image registration workflows, our optimized landmark data show up to a 39.1% reduction in average coordinate error and a 36.7% reduction in total distribution error. In addition, our landmark optimization produces estimates of the sample mean shape and variance-covariance structure that are statistically indistinguishable from expert manual estimates. For biological imaging datasets and morphometric research questions, our approach can eliminate the time and subjectivity of manual landmark detection whilst retaining the biological integrity of these expert annotations.

There are currently no pharmacological approaches in fracture healing designed to therapeutically stimulate endochondral ossification. In this study, we test nerve growth factor (NGF) as an understudied therapeutic for fracture repair. We first characterized endogenous expression of Ngf and its receptor tropomyosin receptor kinase A (TrkA) during tibial fracture repair, finding that they peak during the cartilaginous phase. We then tested two injection regimens and found that local β-NGF injections during the endochondral/cartilaginous phase promoted osteogenic marker expression. Gene expression data from β-NGF stimulated cartilage callus explants show a promotion in markers associated with endochondral ossification such as Ihh, Alpl, and Sdf-1. Gene ontology enrichment analysis revealed the promotion of genes associated with Wnt activation, PDGF- and integrin-binding. Subsequent histological analysis confirmed Wnt activation following local β-NGF injections. Finally, we demonstrate functional improvements to bone healing following local β-NGF injections which resulted in a decrease in cartilage and increase of bone volume. Moreover, the newly formed bone contained higher trabecular number, connective density, and bone mineral density. Collectively, we demonstrate β-NGF’s ability to promote endochondral repair in a murine model and uncover mechanisms that will serve to further understand the molecular switches that occur during cartilage to bone transformation.

Realistic mappings of genes to morphology are inherently multivariate on both sides of the equation. The importance of coordinated gene effects on morphological phenotypes is clear from the intertwining of gene actions in signaling pathways, gene regulatory networks, and developmental processes underlying the development of shape and size. Yet, current approaches tend to focus on identifying and localizing the effects of individual genes and rarely leverage the information content of high dimensional phenotypes. Here, we explicitly model the joint effects of biologically coherent collections of genes on a multivariate trait—craniofacial shape — in a sample of n = 1,145 mice from the Diversity Outbred (DO) experimental line. We use biological process gene ontology (GO) annotations to select skeletal and facial development gene sets and solve for the axis of shape variation that maximally covaries with gene set marker variation. We use our process-centered, multivariate genotype-phenotype (MGP) approach to determine the overall contributions to craniofacial variation of genes involved in relevant processes and how variation in different processes corresponds to multivariate axes of shape variation. Further, we compare the directions of effect in phenotype space of mutations to the primary axis of shape variation associated with broader pathways within which they are thought to function. Finally, we leverage the relationship between mutational and pathway-level effects to predict phenotypic effects beyond craniofacial shape in specific mutants. We also introduce an online application which provides users the means to customize their own process-centered craniofacial shape analyses in the DO. The process-centered approach is generally applicable to any continuously varying phenotype and thus has wide-reaching implications for complex-trait genetics.

Background Trauma is the leading cause of death and disability worldwide, especially in the young population. Cardiac injuries are an independent predictor for a poor overall outcome after trauma. The aim of the present study was to analyze systemic inflammation as well as local cardiac inflammation after experimental limb-, neuro- and combined trauma in mice. Methods Male C57BL/6 mice received either a closed tibia fracture (Fx), isolated traumatic brain injury (TBI) or a combination of both (Fx ​+ ​TBI). Control animals underwent sham procedure. After 6 and 24 ​h, systemic levels of inflammatory mediators were analyzed, respectively. Locally, cardiac inflammation and cardiac structural alterations were investigated in left ventricular tissue of mice 6 and 24 ​h after trauma. Results Mice showed enhanced systemic inflammation after combined trauma, which was manifested by increased levels of KC, MCP-1 and G-CSF. Locally, mice exhibited increased expression of inflammatory cytokines (IL-1β, TNF) in heart tissue, which was probably mediated via toll-like receptor (TLR) signaling. Furthermore, mice demonstrated a redistribution of connexin 43 in cardiac tissue, which appeared predominantly after combined trauma. Besides inflammation and structural cardiac alterations, expression of glucose transporter 4 (GLUT4) mRNA was increased in the heart early after TBI and after combination of TBI and limb fracture, indicating a modification of energy metabolism. Early after combination of TBI and tibia fracture, nitrosative stress was increased, manifested by elevation of nitrotyrosine in cardiac tissue. Finally, mice showed a trend of increased systemic levels of cardiac troponin I and heart-fatty acid binding protein (HFABP) after combined trauma, which was associated with a significant decrease of troponin I and HFABP mRNA expression in cardiac tissue after TBI and combination of TBI and limb fracture. Conclusion Mice exhibited early cardiac alterations as well as alterations in cardiac glucose transporter expression, indicating a modification of energy metabolism, which might be linked to increased systemic- and local cardiac inflammation after limb-, neuro- and combined trauma. These cardiac alterations might predispose individuals for secondary cardiac damage after trauma that might compromise cardiac function after TBI and long bone fracture. Translational potential statement Injuries to the head and extremities frequently occur after severe trauma. In our study, we analyzed the effects of closed tibia fracture, isolated TBI, and the combination of both injuries with regard to the development of post-traumatic secondary cardiac injuries.

Chondrocytes within the fracture callus transform into osteoblasts during bone regeneration, but the molecular mechanisms regulating this process are unknown. Wnt ligands are expressed within the fracture callus, and hypertrophic chondrocytes undergoing transformation to osteoblasts exhibit nuclear localization of β-catenin, indicating active Wnt signaling in these cells. Here, we show that conditional knock out (cKO) of β-catenin in chondrocytes inhibits the transformation of chondrocytes to osteoblasts, while stabilization of β-catenin in chondrocytes accelerates this process. After cKO, chondrocyte-derived cells were located in the bone marrow cavity and upon re-fracture formed cartilage. Lineage tracing in wild type mice revealed that in addition to osteoblasts, chondrocytes give rise to stem cells that contribute to repair of subsequent fractures. These data indicate that Wnt signaling directs cell fate choices of chondrocytes during fracture healing by stimulating transformation of chondrocytes to osteoblasts, and provide a framework for developing Wnt-therapies to stimulate repair.

The fields of developmental biology and tissue engineering have been revolutionized in recent years by technological advancements, expanded understanding, and biomaterials design, leading to the emerging paradigm of "developmental" or "biomimetic" tissue engineering. While developmental biology and tissue engineering have long overlapping histories, the fields have largely diverged in recent years at the same time that crosstalk opportunities for mutual benefit are more salient than ever. In this perspective article, we will use musculoskeletal development and tissue engineering as a platform on which to discuss these emerging crosstalk opportunities and will present our opinions on the bright future of these overlapping spheres of influence. The multicellular programs that control musculoskeletal development are rapidly becoming clarified, represented by shifting paradigms in our understanding of cellular function, identity, and lineage specification during development. Simultaneously, advancements in bioartificial matrices that replicate the biochemical, microstructural, and mechanical properties of developing tissues present new tools and approaches for recapitulating development in tissue engineering. Here, we introduce concepts and experimental approaches in musculoskeletal developmental biology and biomaterials design and discuss applications in tissue engineering as well as opportunities for tissue engineering approaches to inform our understanding of fundamental biology. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:2356-2368, 2017.

Quantitative analysis of morphogenesis aids our understanding of developmental processes by providing a method to link changes in shape with cellular and molecular processes. Over the last decade, many methods have been developed for 3D imaging of embryos using microCT scanning to quantify the shape of embryos during development. These methods generally involve a powerful, cross-linking fixative such as paraformaldehyde to limit shrinkage during the CT scan. However, the extended time frames that these embryos are incubated in such fixatives prevent use of the tissues for molecular analysis after microCT scanning. This is a significant problem because it limits the ability to correlate variation in molecular data with morphology at the level of individual embryos. Here we outline a novel method that allows RNA, DNA, or protein isolation following CT scan while also allowing imaging of different tissue layers within the developing embryo. We show shape differences early in craniofacial development (E11.5) between common mouse genetic backgrounds, and demonstrate that we are able to generate RNA from these embryos after CT scanning that is suitable for downstream real time PCR (RT-PCR) and RNAseq analyses. Developmental Dynamics 246:431-436, 2017. © 2017 Wiley Periodicals, Inc.

The elderly population suffers from higher rates of complications during fracture healing that result in increased morbidity and mortality. Inflammatory dysregulation is associated with increased age and is a contributing factor to the myriad of age‐related diseases. Therefore, we investigated age‐related changes to an important cellular regulator of inflammation, the macrophage, and the impact on fracture healing outcomes. We demonstrated that old mice (24 months) have delayed fracture healing with significantly less bone and more cartilage compared to young mice (3 months). The quantity of infiltrating macrophages into the fracture callus was similar in old and young mice. However, RNA‐seq analysis demonstrated distinct differences in the transcriptomes of macrophages derived from the fracture callus of old and young mice, with an up‐regulation of M1/pro‐inflammatory genes in macrophages from old mice as well as dysregulation of other immune‐related genes. Preventing infiltration of the fracture site by macrophages in old mice improved healing outcomes, with significantly more bone in the calluses of treated mice compared to age‐matched controls. After preventing infiltration by macrophages, the macrophages remaining within the fracture callus were collected and examined via RNA‐seq analysis, and their transcriptome resembled macrophages from young calluses. Taken together, infiltrating macrophages from old mice demonstrate detrimental age‐related changes, and depleting infiltrating macrophages can improve fracture healing in old mice., Delayed fracture healing in old mice is associated with transcriptomic differences between old and young macrophages infiltrating the fracture callus. An aged macrophage phenotype was characterized by increased pro‐inflammatory gene expression compared to young. Inhibition of macrophage infiltration in old mice improved fracture healing outcomes and demonstrated a resident macrophages population with a transcriptome more similar to young mice.

Background Cleft lip and palate is one of the most common human birth defects, but the underlying etiology is poorly understood. The A/WySn mouse is a spontaneously occurring model of multigenic clefting in which 20% to 30% of individuals develop an orofacial cleft. Recent work has shown altered methylation at a specific retrotransposon insertion downstream of the Wnt9b locus in clefting animals, which results in decreased Wnt9b expression. Results Using a newly developed protocol that allows us to measure morphology, gene expression, and DNA methylation in the same embryo, we relate gene expression in an individual embryo directly to its three-dimensional morphology for the first time. We find that methylation at the retrotransposon relates to Wnt9b expression and morphology. IAP methylation relates to shape of the nasal process in a manner consistent with clefting. Embryos with low IAP methylation exhibit increased among-individual variance in facial shape. Conclusions Methylation and gene expression relate nonlinearly to nasal process morphology. Individuals at one end of a continuum of phenotypic states display a clinical phenotype and increased phenotypic variation. Variable penetrance and expressivity in this model is likely determined both by among-individual variation in methylation and changes in phenotypic robustness along the underlying liability distribution for orofacial clefting.

BackgroundPolarity is necessary for epithelial cells to perform distinct functions at their apical and basal surfaces. Oral epithelial cell-derived ameloblasts at secretory stage (SABs) synthesize large amounts of enamel matrix proteins (EMPs), largely amelogenins. EMPs are unidirectionally secreted into the enamel space through their apical cytoplasmic protrusions, or Tomes’ processes (TPs), to guide the enamel formation. Little is known about the transcriptional regulation underlying the establishment of cell polarity and unidirectional secretion of SABs.ResultsThe higher-order chromatin architecture of eukaryotic genome plays important roles in cell- and stage-specific transcriptional programming. A genome organizer, special AT-rich sequence-binding protein 1 (SATB1), was discovered to be significantly upregulated in ameloblasts compared to oral epithelial cells using a whole-transcript microarray analysis. TheSatb1−/−mice possessed deformed ameloblasts and a thin layer of hypomineralized and non-prismatic enamel. Remarkably,Satb1−/−ameloblasts at the secretory stage lost many morphological characteristics found at the apical surface of wild-type (wt)SABs, including the loss of Tomes’ processes, defective inter-ameloblastic adhesion, and filamentous actin architecture. As expected, the secretory function ofSatb1−/−SABs was compromised as amelogenins were largely retained in cells. We found the expression of epidermal growth factor receptor pathway substrate 8 (Eps8), a known regulator for actin filament assembly and small intestinal epithelial cytoplasmic protrusion formation, to be SATB1 dependent. In contrast towtSABs, EPS8 could not be detected at the apical surface ofSatb1−/−SABs.Eps8expression was greatly reduced in small intestinal epithelial cells inSatb1−/−mice as well, displaying defective intestinal microvilli.ConclusionsOur data show that SATB1 is essential for establishing secretory ameloblast cell polarity and for EMP secretion. In line with the deformed apical architecture, amelogenin transport to the apical secretory front and secretion into enamel space were impeded inSatb1−/−SABs resulting in a massive cytoplasmic accumulation of amelogenins and a thin layer of hypomineralized enamel. Our studies strongly suggest that SATB1-dependentEps8expression plays a critical role in cytoplasmic protrusion formation in both SABs and in small intestines. This study demonstrates the role of SATB1 in the regulation of amelogenesis and the potential application of SATB1 in ameloblast/enamel regeneration.

The biology of how faces are built and come to differ from one another is complex. Discovering the genes that contribute to differences in facial morphology is one key to untangling this complexity, with important implications for medicine and evolutionary biology. This study maps quantitative trait loci (QTL) for skeletal facial shape using Diversity Outbred (DO) mice. The DO is a randomly outcrossed population with high heterozygosity that captures the allelic diversity of eight inbred mouse lines from three subspecies. The study uses a sample of 1147 DO animals (the largest sample yet employed for a shape QTL study in mouse), each characterized by 22 three-dimensional landmarks, 56,885 autosomal and X-chromosome markers, and sex and age classifiers. We identified 37 facial shape QTL across 20 shape principal components (PCs) using a mixed effects regression that accounts for kinship among observations. The QTL include some previously identified intervals as well as new regions that expand the list of potential targets for future experimental study. Three QTL characterized shape associations with size (allometry). Median support interval size was 3.5 Mb. Narrowing additional analysis to QTL for the five largest magnitude shape PCs, we found significant overrepresentation of genes with known roles in growth, skeletal development, and sensory organ development. For most intervals, one or more of these genes lies within 0.25 Mb of the QTL’s peak. QTL effect sizes were small, with none explaining more than 0.5% of facial shape variation. Thus, our results are consistent with a model of facial diversity that is influenced by key genes in skeletal and facial development and, simultaneously, is highly polygenic.Author SummaryThe mammalian face is a complex structure serving many functions. We studied the genetic basis for facial skeletal diversity in a large sample of mice from an experimental population designed for the study of complex traits. We quantified the contribution of genetic variation to variation in three-dimensional facial shape across more than 55,000 genetic markers spread throughout the mouse genome. We found 37 genetic regions which are very likely to contribute to differences in facial shape. We then conducted a more detailed analysis of the genetic regions associated with the most variable aspects of facial shape. For these regions, a disproportionately large number of genes are known to be important to growth and to skeletal and facial development. The magnitude of these genetic contributions to differences in facial shape are consistently small. Our results therefore support the notion that facial skeletal diversity is influenced by many genes of small effect, but that some of these small effects may be related to genes that are fundamental to skeletal and facial development.

The elderly population suffers from higher rates of complications during fracture healing that result in increased morbidity and mortality. Inflammatory dysregulation is associated with increased age and is a contributing factor to the myriad of age-related diseases. Therefore, we investigated age-related changes to an important cellular regulator of inflammation, the macrophage, and the impact on fracture healing outcomes. We demonstrated that old mice (24 months) have delayed fracture healing with significantly less bone and more cartilage compared to young mice (3 months). The quantity of infiltrating macrophages into the fracture callus was similar in old and young mice. However, RNA-seq analysis demonstrated distinct differences in the transcriptomes of macrophages derived from the fracture callus of old and young mice, with an upregulation of M1/pro-inflammatory genes in macrophages from old mice as well as dysregulation of other immune-related genes. Preventing infiltration of the fracture site by macrophages in old mice improved healing outcomes, with significantly more bone in the calluses of treated mice compared to age-matched controls. After preventing infiltration by macrophages, the macrophages within the fracture callus were collected and examined via RNA-seq analysis, and their transcriptome resembled macrophages from young calluses. Taken together, infiltrating macrophages from old mice demonstrate detrimental age-related changes, and depleting infiltrating macrophages can improve fracture healing in old mice.

Allometry refers to the ways in which organismal shape is associated with size. It is a special case of integration, or the tendency for traits to covary, in that variation in size is ubiquitous and evolutionarily important. Allometric variation is so commonly observed that it is routinely removed from morphometric analyses or invoked as an explanation for evolutionary change. In this case, familiarity is mistaken for understanding because rarely do we know the mechanisms by which shape correlates with size or understand their significance. As with other forms of integration, allometric variation is generated by variation in developmental processes that affect multiple traits, resulting in patterns of covariation. Given this perspective, we can dissect the genetic and developmental determinants of allometric variation. Our work on the developmental and genetic basis for allometric variation in craniofacial shape in mice and humans has revealed that allometric variation is highly polygenic. Different measures of size are associated with distinct but overlapping patterns of allometric variation. These patterns converge in part on a common genetic basis. Finally, environmental modulation of size often generates variation along allometric trajectories, but the timing of genetic and environmental perturbations can produce deviations from allometric patterns when traits are differentially sensitive over developmental time. These results question the validity of viewing allometry as a singular phenomenon distinct from morphological integration more generally.

Chronic psychosocial stress/trauma represents an increasing burden in our modern society and a risk factor for the development of mental disorders, including posttraumatic stress disorder (PTSD). PTSD, in turn, is highly comorbid with a plethora of inflammatory disorders and has been associated with increased bone fracture risk. Since a balanced inflammatory response after fracture is crucial for successful bone healing, we hypothesize that stress/trauma alters the inflammatory response after fracture and, consequently, compromises fracture healing. Here we show, employing the chronic subordinate colony housing (CSC) paradigm as a clinically relevant mouse model for PTSD, that mice subjected to CSC displayed increased numbers of neutrophils in the early fracture hematoma, whereas T lymphocytes and markers for cartilage-to-bone transition and angiogenesis were reduced. At late stages of fracture healing, CSC mice were characterized by decreased bending stiffness and bony bridging of the fracture callus. Strikingly, a single systemic administration of the β-adrenoreceptor (AR) blocker propranolol before femur osteotomy prevented bone marrow mobilization of neutrophils and invasion of neutrophils into the fracture hematoma, both seen in the early phase after fracture, as well as a compromised fracture healing in CSC mice. We conclude that chronic psychosocial stress leads to an imbalanced immune response after fracture via β-AR signaling, accompanied by disturbed fracture healing. These findings offer possibilities for clinical translation in patients suffering from PTSD and fracture.