Developmental and genetic origins of murine long bone length variation

Perspective: A multi-trait integrative approach to understanding the structural basis of bone fragility for pediatric conditions associated with abnormal bone development

Bone development is a highly orchestrated process that establishes the structural basis of bone strength during growth and functionality across the lifespan. This developmental process is generally robust in establishing mechanical function, being adaptable to many genetic and environmental factors. However, not all factors can be fully accommodated, leading to abnormal bone development and lower bone strength. This can give rise to early-onset bone fragility that negatively impacts bone strength across the lifespan. Current guidelines for assessing bone strength include measuring bone mineral density, but this does not capture the structural details responsible for whole bone strength in abnormally developing bones that would be needed to inform clinicians on how and when to treat to improve bone strength. The clinical consequence of not operationalizing how altered bone development informs decision making includes under-detection and missed opportunities for early intervention, as well as a false positive diagnosis of fragility with possible resultant clinical actions that may actually harm the growing skeleton. In this Perspective, we emphasize the need for a multi-trait, integrative approach to better understand the structural basis of bone growth for pediatric conditions with abnormal bone development. We provide evidence to showcase how this approach might reveal multiple, unique ways in which bone fragility develops across and within an array of pediatric conditions that are associated with abnormal bone development. This Perspective advocates for the development of new translational research aimed at informing better ways to optimize bone growth, prevent fragility fractures, and monitor and treat bone fragility based on the child's skeletal needs.

Early Life Exposure to a Diet With a Supramolecular Lipid Structure Close to That of Mammalian Milk Improves Early Life Growth, Skeletal Development, and Later Life Neurocognitive Function in Individually and Socially Housed Male C57BL/6J Mice

Breastfeeding (duration) can be positively associated with infant growth outcomes as well as improved cognitive functions during childhood and later life stages. (Prolonged) exposure to optimal lipid quantity and quality, i.e., the supramolecular structure of lipids, in mammalian milk, may contribute to these beneficial effects through nutritional early-life programming. In this pre-clinical study, we exposed male C57BL/6J mice from post-natal Days 16 to 42 (i.e., directly following normal lactation), to a diet with large lipid droplets coated with bovine milk fat globule membrane-derived phospholipids, which mimic more closely the supramolecular structure of lipid droplets in mammalian milk. We investigated whether exposure to this diet could affect growth and brain development-related parameters. As these outcomes are also known to be affected by the post-weaning social environment in mice, we included both individually housed and pair-wise housed animals and studied whether effects of diet were modulated by the social environment. After Day 42, all the animals were fed standard semi-synthetic rodent diet. Growth and body composition were assessed, and the mice were subjected to various behavioral tests. Individual housing attenuated adolescent growth, reduced femur length, and increased body fat mass. Adult social interest was increased due to individual housing, while cognitive and behavioral alterations as a result of different housing conditions were modest. The diet increased adolescent growth and femur length, increased lean body mass, reduced adolescent anxiety, and improved adult cognitive performance. These effects of diet exposure were comparable between individually and socially housed mice. Hence, early life exposure to a diet with lipid droplets that mimic the supramolecular structure of those in mammalian milk may improve adolescent growth and alters brain function in both socially and individually housed mice. These findings suggest that lipid structure in infant milk formula may be a relevant target for nutritional solutions, targeting both healthy infants and infants facing growth challenges.

BMP3 Affects Cortical and Trabecular Long Bone Development in Mice

Bone morphogenetic proteins (BMPs) have a major role in tissue development. BMP3 is synthesized in osteocytes and mature osteoblasts and has an antagonistic effect on other BMPs in bone tissue. The main aim of this study was to fully characterize cortical bone and trabecular bone of long bones in both male and female Bmp3^−/− mice. To investigate the effect of Bmp3 from birth to maturity, we compared Bmp3^−/− mice with wild-type littermates at the following stages of postnatal development: 1 day (P0), 2 weeks (P14), 8 weeks and 16 weeks of age. Bmp3 deletion was confirmed using X-gal staining in P0 animals. Cartilage and bone tissue were examined in P14 animals using Alcian Blue/Alizarin Red staining. Detailed long bone analysis was performed in 8-week-old and 16-week-old animals using micro-CT. The Bmp3 reporter signal was localized in bone tissue, hair follicles, and lungs. Bone mineralization at 2 weeks of age was increased in long bones of Bmp3^−/− mice. Bmp3 deletion was shown to affect the skeleton until adulthood, where increased cortical and trabecular bone parameters were found in young and adult mice of both sexes, while delayed mineralization of the epiphyseal growth plate was found in adult Bmp3^−/− mice.

A Na+/K+ ATPase Pump Regulates Chondrocyte Differentiation and Bone Length Variation in Mice

The genetic and developmental mechanisms involved in limb formation are relatively well documented, but how these mechanisms are modulated by changes in chondrocyte physiology to produce differences in limb bone length remains unclear. Here, we used high throughput RNA sequencing (RNAseq) to probe the developmental genetic basis of variation in limb bone length in Longshanks, a mouse model of experimental evolution. We find that increased tibia length in Longshanks is associated with altered expression of a few key endochondral ossification genes such as Npr3, Dlk1, Sox9, and Sfrp1, as well reduced expression of Fxyd2, a facultative subunit of the cell membrane-bound Na⁺/K⁺ ATPase pump (NKA). Next, using murine tibia and cell cultures, we show a dynamic role for NKA in chondrocyte differentiation and in bone length regulation. Specifically, we show that pharmacological inhibition of NKA disrupts chondrocyte differentiation, by upregulating expression of mesenchymal stem cell markers (Prrx1, Serpina3n), downregulation of chondrogenesis marker Sox9, and altered expression of extracellular matrix genes (e.g., collagens) associated with proliferative and hypertrophic chondrocytes. Together, Longshanks and in vitro data suggest a broader developmental and evolutionary role of NKA in regulating limb length diversity.

Sex-specific effects of neonatal oral sucrose treatment on growth and liver choline and glucocorticoid metabolism in adulthood

Hospitalized preterm infants experience painful medical procedures. Oral sucrose is the nonpharmacological standard of care for minor procedural pain relief. Infants are treated with numerous doses of sucrose, raising concerns about potential long-term effects. The objective of this study was to determine the long-term effects of neonatal oral sucrose treatment on growth and liver metabolism in a mouse model. Neonatal female and male mice were randomly assigned to one of two oral treatments (n = 7-10 mice/group/sex): sterile water or sucrose. Pups were treated 10 times/day for the first 6 days of life with 0.2 mg/g body wt of respective treatments (24% solution; 1-4 μL/dose) to mimic what is given to preterm infants. Mice were weaned at age 3 wk onto a control diet and fed until age 16 wk. Sucrose-treated female and male mice gained less weight during the treatment period and were smaller at weaning than water-treated mice (P ≤ 0.05); no effect of sucrose treatment on body weight was observed at adulthood. However, adult sucrose-treated female mice had smaller tibias and lower serum insulin-like growth factor-1 than adult water-treated female mice (P ≤ 0.05); these effects were not observed in males. Lower liver S-adenosylmethionine, phosphocholine, and glycerophosphocholine were observed in adult sucrose-treated compared with water-treated female and male mice (P ≤ 0.05). Sucrose-treated female, but not male, mice had lower liver free choline and higher liver betaine compared with water-treated female mice (P < 0.01). Our findings suggest that repeated neonatal sucrose treatment has long-term sex-specific effects on growth and liver methionine and choline metabolism.

A New Pipeline to Automatically Segment and Semi-Automatically Measure Bone Length on 3D Models Obtained by Computed Tomography

The characterization of developmental phenotypes often relies on the accurate linear measurement of structures that are small and require laborious preparation. This is tedious and prone to errors, especially when repeated for the multiple replicates that are required for statistical analysis, or when multiple distinct structures have to be analyzed. To address this issue, we have developed a pipeline for characterization of long-bone length using X-ray microtomography (XMT) scans. The pipeline involves semi-automated algorithms for automatic thresholding and fast interactive isolation and 3D-model generation of the main limb bones, using either the open-source ImageJ plugin BoneJ or the commercial Mimics Innovation Suite package. The tests showed the appropriate combination of scanning conditions and analysis parameters yields fast and comparable length results, highly correlated with the measurements obtained via ex vivo skeletal preparations. Moreover, since XMT is not destructive, the samples can be used afterward for histology or other applications. Our new pipelines will help developmental biologists and evolutionary researchers to achieve fast, reproducible and non-destructive length measurement of bone samples from multiple animal species.

Long bone histomorphogenesis of the naked mole-rat: Histodiversity and intraspecific variation

Lacking fur, living in eusocial colonies and having the longest lifespan of any rodent, makes naked mole-rats (NMRs) rather peculiar mammals. Although they exhibit a high degree of polymorphism, skeletal plasticity and are considered a novel model to assess the effects of delayed puberty on the skeletal system, scarce information on their morphogenesis exists. Here, we examined a large ontogenetic sample (n = 76) of subordinate individuals to assess the pattern of bone growth and bone microstructure of fore- and hindlimb bones by using histomorphological techniques. Over 290 undecalcified thin cross-sections from the midshaft of the humerus, ulna, femur, and tibia from pups, juveniles and adults were analyzed with polarized light microscopy. Similar to other fossorial mammals, NMRs exhibited a systematic cortical thickening of their long bones, which clearly indicates a conserved functional adaptation to withstand the mechanical strains imposed during digging, regardless of their chisel-tooth predominance. We describe a high histodiversity of bone matrices and the formation of secondary osteons in NMRs. The bones of pups are extremely thin-walled and grow by periosteal bone formation coupled with considerable expansion of the medullary cavity, a process probably tightly regulated and adapted to optimize the amount of minerals destined for skeletal development, to thus allow the female breeder to produce a higher number of pups, as well as several litters. Subsequent cortical thickening in juveniles involves high amounts of endosteal bone apposition, which contrasts with the bone modeling of other mammals where a periosteal predominance exists. Adults have bone matrices predominantly consisting of parallel-fibered bone and lamellar bone, which indicate intermediate to slow rates of osteogenesis, as well as the development of poorly vascularized lamellar-zonal tissues separated by lines of arrested growth (LAGs) and annuli. These features reflect the low metabolism, low body temperature and slow growth rates reported for this species, as well as indicate a cyclical pattern of osteogenesis. The presence of LAGs in captive individuals was striking and indicates that postnatal osteogenesis and its consequent cortical stratification most likely represents a plesiomorphic thermometabolic strategy among endotherms which has been suggested to be regulated by endogenous rhythms. However, the generalized presence of LAGs in this and other subterranean taxa in the wild, as well as recent investigations on variability of environmental conditions in burrow systems, supports the hypothesis that underground environments experience seasonal fluctuations that may influence the postnatal osteogenesis of animals by limiting the extension of burrow systems during the unfavorable dry seasons and therefore the finding of food resources. Additionally, the intraspecific variation found in the formation of bone tissue matrices and vascularization suggested a high degree of developmental plasticity in NMRs, which may help explaining the polymorphism reported for this species. The results obtained here represent a valuable contribution to understanding the relationship of several aspects involved in the morphogenesis of the skeletal system of a mammal with extraordinary adaptations.

Individual housing of male C57BL/6J mice after weaning impairs growth and predisposes for obesity

For (metabolic) research models using mice, singly housing is widely used for practical purposes to study e.g. energy balance regulation and derangements herein. Mouse (social) housing practices could however influence study results by modulating (metabolic) health outcomes. To study the effects of the social housing condition, we assessed parameters for energy balance regulation and proneness to (diet induced) obesity in male C57Bl/6J mice that were housed individually or socially (in pairs) directly after weaning, both at standard ambient temperature of 21°C. During adolescence, individually housed mice had reduced growth rate, while energy intake and energy expenditure were increased compared to socially housed counterparts. At 6 weeks of age, these mice had reduced lean body mass, but significantly higher white adipose tissue mass compared to socially housed mice, and higher UCP-1 mRNA expression in brown adipose tissue. During adulthood, body weight gain of individually housed animals exceeded that of socially housed mice, with elevations in both energy intake and expenditure. At 18 weeks of age, individually housed mice showed higher adiposity and higher mRNA expression of UCP-1 in inguinal white but not in brown adipose tissue. Exposure to an obesogenic diet starting at 6 weeks of age further amplified body weight gain and adipose tissue deposition and caused strong suppression of inguinal white adipose tissue mRNA UCP-1 expression. This study shows that post-weaning individual housing of male mice impairs adolescent growth and results in higher susceptibility to obesity in adulthood with putative roles for thermoregulation and/or affectiveness.

RNA-seq analysis of chondrocyte transcriptome reveals genetic heterogeneity in LG/J and SM/J murine strains

OBJECTIVE

To investigate the transcriptomic differences in chondrocytes obtained from LG/J (large, healer) and SM/J (small, non-healer) murine strains in an attempt to discern the molecular pathways implicated in cartilage regeneration and susceptibility to osteoarthritis (OA).

DESIGN

We performed RNA-sequencing on chondrocytes derived from LG/J (n = 16) and SM/J (n = 16) mice. We validated the expression of candidate genes and compared single nucleotide polymorphisms (SNPs) between the two mouse strains. We also examined gene expression of positional candidates for ear pinna regeneration and long bone length quantitative trait loci (QTLs) that display differences in cartilaginous expression.

RESULTS

We observed a distinct genetic heterogeneity between cells derived from LG/J and SM/J mouse strains. We found that gene ontologies representing cell development, cartilage condensation, and regulation of cell differentiation were enriched in LG/J chondrocytes. In contrast, gene ontologies enriched in the SM/J chondrocytes were mainly related to inflammation and degeneration. Moreover, SNP analysis revealed that multiple validated genes vary in sequence between LG/J and SM/J in coding and highly conserved noncoding regions. Finally, we showed that most QTLs have 20-30% of their positional candidates displaying differential expression between the two mouse strains.

CONCLUSIONS

While the enrichment of pathways related to cell differentiation, cartilage development and cartilage condensation infers superior healing potential of LG/J strain, the enrichment of pathways related to cytokine production, immune cell activation and inflammation entails greater susceptibility of SM/J strain to OA. These data provide novel insights into chondrocyte transcriptome and aid in identification of the quantitative trait genes and molecular differences underlying the phenotypic differences associated with individual QTLs.

Variation in mouse pelvic morphology maps to locations enriched in Sox9 Class II and Pitx1 regulatory features

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high-resolution quantitative trait locus (QTL) mapping in the F₃₄ generation of an advanced intercross experiment (LG,SM-G₃₄ ) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome- or chromosome-wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation-seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

Endochondral ossification and the evolution of limb proportions

Mammals have remarkably diverse limb proportions hypothesized to have evolved adaptively in the context of locomotion and other behaviors. Mechanistically, evolutionary diversity in limb proportions is the result of differential limb bone growth. Longitudinal limb bone growth is driven by the process of endochondral ossification, under the control of the growth plates. In growth plates, chondrocytes undergo a tightly orchestrated life cycle of proliferation, matrix production, hypertrophy, and cell death/transdifferentiation. This life cycle is highly conserved, both among the long bones of an individual, and among homologous bones of distantly related taxa, leading to a finite number of complementary cell mechanisms that can generate heritable phenotype variation in limb bone size and shape. The most important of these mechanisms are chondrocyte population size in chondrogenesis and in individual growth plates, proliferation rates, and hypertrophic chondrocyte size. Comparative evidence in mammals and birds suggests the existence of developmental biases that favor evolutionary changes in some of these cellular mechanisms over others in driving limb allometry. Specifically, chondrocyte population size may evolve more readily in response to selection than hypertrophic chondrocyte size, and extreme hypertrophy may be a rarer evolutionary phenomenon associated with highly specialized modes of locomotion in mammals (e.g., powered flight, ricochetal bipedal hopping). Physical and physiological constraints at multiple levels of biological organization may also have influenced the cell developmental mechanisms that have evolved to produce the highly diverse limb proportions in extant mammals. This article is categorized under: Establishment of Spatial and Temporal Patterns > Regulation of Size, Proportion, and Timing Comparative Development and Evolution > Regulation of Organ Diversity Comparative Development and Evolution > Organ System Comparisons Between Species.

How a growing organismal perspective is adding new depth to integrative studies of morphological evolution

Over the past half century, the field of Evolutionary Developmental Biology, or Evo-devo, has integrated diverse fields of biology into a more synthetic understanding of morphological diversity. This has resulted in numerous insights into how development can evolve and reciprocally influence morphological evolution, as well as generated several novel theoretical areas. Although comparative by default, there remains a great gap in our understanding of adaptive morphological diversification and how developmental mechanisms influence the shape and pattern of phenotypic variation. Herein we highlight areas of research that are in the process of filling this void, and areas, if investigated more fully, that will add new insights into the diversification of morphology. At the centre of our discussion is an explicit awareness of organismal biology. Here we discuss an organismal framework that is supported by three distinct pillars. First, there is a need for Evo-devo to adopt a high-resolution phylogenetic approach in the study of morphological variation and its developmental underpinnings. Secondly, we propose that to understand the dynamic nature of morphological evolution, investigators need to give more explicit attention to the processes that generate evolutionarily relevant variation at the population level. Finally, we emphasize the need to address more thoroughly the processes that structure variation at micro- and macroevolutionary scales including modularity, morphological integration, constraint, and plasticity. We illustrate the power of these three pillars using numerous examples from both invertebrates and vertebrates to emphasize that many of these approaches are already present within the field, but have yet to be formally integrated into many research programs. We feel that the most exciting new insights will come where the traditional experimental approaches to Evo-devo are integrated more thoroughly with the principles of this organismal framework.

Evidence for Genetic Contribution to Variation in Posttraumatic Osteoarthritis in Mice

OBJECTIVE

Recombinant inbred mouse strains generated from an LG/J and SM/J intercross offer a unique resource to study complex genetic traits such as osteoarthritis (OA). We undertook this study to determine the susceptibility of 14 strains to various phenotypes characteristic of posttraumatic OA. We hypothesized that phenotypic variability is associated with genetic variability.

METHODS

Ten-week-old male mice underwent surgical destabilization of the medial meniscus (DMM) to induce posttraumatic OA. Mice were killed 8 weeks after surgery, and knee joints were processed for histology to score cartilage degeneration and synovitis. Micro-computed tomography was used to analyze trabecular bone parameters including subchondral bone plate thickness and synovial ectopic calcifications. Gene expression in the knees was assessed using a QuantiGene Plex assay.

RESULTS

Broad-sense heritability ranged from 0.18 to 0.58, which suggested that the responses to surgery were moderately heritable. The LGXSM-33, LGXSM-5, LGXSM-46, and SM/J strains were highly susceptible to OA, while the LGXSM-131b, LGXSM-163, LGXSM-35, LGXSM-128a, LGXSM-6, and LG/J strains were relatively OA resistant. This study was the first to accomplish measurement of genetic correlations of phenotypes that are characteristic of posttraumatic OA. Cartilage degeneration was significantly positively associated with synovitis (r = 0.83-0.92), and subchondral bone plate thickness was negatively correlated with ectopic calcifications (r = -0.59). Moreover, we showed that 40 of the 78 genes tested were significantly correlated with various OA phenotypes. However, unlike the OA phenotypes, there was no evidence for genetic variation in differences in gene expression levels between DMM-operated and sham-operated knees.

CONCLUSION

For these mouse strains, various characteristics of posttraumatic OA varied with genetic composition, which demonstrated a genetic basis for susceptibility to posttraumatic OA. The heritability of posttraumatic OA was established. Phenotypes exhibited various degrees of correlations; cartilage degeneration was positively correlated with synovitis, but not with the formation of ectopic calcifications. Further investigation of the genome regions that contain genes implicated in OA, as well as further investigation of gene expression data, will be useful for studying mechanisms of OA and identifying therapeutic targets.

Developmental Bias and Evolution: A Regulatory Network Perspective

Phenotypic variation is generated by the processes of development, with some variants arising more readily than others-a phenomenon known as "developmental bias." Developmental bias and natural selection have often been portrayed as alternative explanations, but this is a false dichotomy: developmental bias can evolve through natural selection, and bias and selection jointly influence phenotypic evolution. Here, we briefly review the evidence for developmental bias and illustrate how it is studied empirically. We describe recent theory on regulatory networks that explains why the influence of genetic and environmental perturbation on phenotypes is typically not uniform, and may even be biased toward adaptive phenotypic variation. We show how bias produced by developmental processes constitutes an evolving property able to impose direction on adaptive evolution and influence patterns of taxonomic and phenotypic diversity. Taking these considerations together, we argue that it is not sufficient to accommodate developmental bias into evolutionary theory merely as a constraint on evolutionary adaptation. The influence of natural selection in shaping developmental bias, and conversely, the influence of developmental bias in shaping subsequent opportunities for adaptation, requires mechanistic models of development to be expanded and incorporated into evolutionary theory. A regulatory network perspective on phenotypic evolution thus helps to integrate the generation of phenotypic variation with natural selection, leaving evolutionary biology better placed to explain how organisms adapt and diversify.

Genetic Dissection of a Supergene Implicates Tfap2a in Craniofacial Evolution of Threespine Sticklebacks

In nature, multiple adaptive phenotypes often coevolve and can be controlled by tightly linked genetic loci known as supergenes. Dissecting the genetic basis of these linked phenotypes is a major challenge in evolutionary genetics. Multiple freshwater populations of threespine stickleback fish (Gasterosteus aculeatus) have convergently evolved two constructive craniofacial traits, longer branchial bones and increased pharyngeal tooth number, likely as adaptations to dietary differences between marine and freshwater environments. Prior QTL mapping showed that both traits are partially controlled by overlapping genomic regions on chromosome 21 and that a regulatory change in Bmp6 likely underlies the tooth number QTL. Here, we mapped the branchial bone length QTL to a 155 kb, eight-gene interval tightly linked to, but excluding the coding regions of Bmp6 and containing the candidate gene Tfap2a Further recombinant mapping revealed this bone length QTL is separable into at least two loci. During embryonic and larval development, Tfap2a was expressed in the branchial bone primordia, where allele specific expression assays revealed the freshwater allele of Tfap2a was expressed at lower levels relative to the marine allele in hybrid fish. Induced loss-of-function mutations in Tfap2a revealed an essential role in stickleback craniofacial development and show that bone length is sensitive to Tfap2a dosage in heterozygotes. Combined, these results suggest that closely linked but genetically separable changes in Bmp6 and Tfap2a contribute to a supergene underlying evolved skeletal gain in multiple freshwater stickleback populations.

Artificial selection sheds light on developmental mechanisms of limb elongation

Species diversity in limb lengths and proportions is thought to have evolved adaptively in the context of locomotor and habitat specialization, but the heritable cellular processes that drove this evolution within species are poorly understood. In this study, we take a novel "micro-evo-devo" approach, using artificial selection on relative limb length to amplify phenotypic variation in a population of mice, known as Longshanks, to examine the cellular mechanisms of postnatal limb development that contribute to intraspecific limb length variation. Cross-sectional growth data indicate that differences in bone length between Longshanks and random-bred controls are not due to prolonged growth, but to accelerated growth rates. Histomorphometric and cell proliferation assays on proximal tibial growth plates show that Longshanks' increased limb bone length is associated with an increased number of proliferative chondrocytes. In contrast, we find no differences in other growth plate cellular features known to underlie interspecific differences in limb bone size and shape, such as the rates of chondrocyte proliferation or the size and number of hypertrophic cells in the growth plate. These data suggest that small differences among individuals in the number of proliferating chondrocytes are a potentially important determinant of selectable intraspecific variation in individual limb bone lengths, independent of body size.

Alterations to the Gut Microbiome Impair Bone Strength and Tissue Material Properties

Alterations in the gut microbiome have been associated with changes in bone mass and microstructure, but the effects of the microbiome on bone biomechanical properties are not known. Here we examined bone strength under two conditions of altered microbiota: (1) an inbred mouse strain known to develop an altered gut microbiome due to deficits in the immune system (the Toll-like receptor 5-deficient mouse [TLR5KO]); and (2) disruption of the gut microbiota (ΔMicrobiota) through chronic treatment with selected antibiotics (ampicillin and neomycin). The bone phenotypes of TLR5KO and WT (C57Bl/6) mice were examined after disruption of the microbiota from 4 weeks to 16 weeks of age as well as without treatment (n = 7 to 16/group, 39 animals total). Femur bending strength was less in ΔMicrobiota mice than in untreated animals and the reduction in strength was not fully explained by differences in bone cross-sectional geometry, implicating impaired bone tissue material properties. Small differences in whole-bone bending strength were observed between WT and TLR5KO mice after accounting for differences in bone morphology. No differences in trabecular bone volume fraction were associated with genotype or disruption of gut microbiota. Treatment altered the gut microbiota by depleting organisms from the phyla Bacteroidetes and enriching for Proteobacteria, as determined from sequencing of fecal 16S rRNA genes. Differences in splenic immune cell populations were also observed; B and T cell populations were depleted in TLR5KO mice and in ΔMicrobiota mice (p < 0.001), suggesting an association between alterations in bone tissue material properties and immune cell populations. We conclude that alterations in the gut microbiota for extended periods during growth may lead to impaired whole-bone mechanical properties in ways that are not explained by bone geometry. © 2017 American Society for Bone and Mineral Research.

Genetics of Skeletal Evolution in Unusually Large Mice from Gough Island

Organisms on islands often undergo rapid morphological evolution, providing a platform for understanding mechanisms of phenotypic change. Many examples of evolution on islands involve the vertebrate skeleton. Although the genetic basis of skeletal variation has been studied in laboratory strains, especially in the house mouse Mus musculus domesticus, the genetic determinants of skeletal evolution in natural populations remain poorly understood. We used house mice living on the remote Gough Island-the largest wild house mice on record-to understand the genetics of rapid skeletal evolution in nature. Compared to a mainland reference strain from the same subspecies (WSB/EiJ), the skeleton of Gough Island mice is considerably larger, with notable expansions of the pelvis and limbs. The Gough Island mouse skeleton also displays changes in shape, including elongations of the skull and the proximal vs. distal elements in the limbs. Quantitative trait locus (QTL) mapping in a large F₂ intercross between Gough Island mice and WSB/EiJ reveals hundreds of QTL that control skeletal dimensions measured at 5, 10, and/or 16 weeks of age. QTL exhibit modest, mostly additive effects, and Gough Island alleles are associated with larger skeletal size at most QTL. The QTL with the largest effects are found on a few chromosomes and affect suites of skeletal traits. Many of these loci also colocalize with QTL for body weight. The high degree of QTL colocalization is consistent with an important contribution of pleiotropy to skeletal evolution. Our results provide a rare portrait of the genetic basis of skeletal evolution in an island population and position the Gough Island mouse as a model system for understanding mechanisms of rapid evolution in nature.

Characterizing phenotypic divergence using three-dimensional geometric morphometrics in four populations of threespine stickleback (Gasterosteus aculeatus; Pisces: Gasterosteidae) in Katmai National Park and Preserve, Alaska

The threespine stickleback (Gasterosteus aculeatus L., 1758) is a vertebrate model for the study of the relationship between phenotype and environment in facilitating rapid evolutionary change. Using four populations from a system of lakes in Katmai National Park and Preserve, Alaska, and microcomputed tomography and three-dimensional geometric morphometrics, we test the hypothesis that stickleback populations inhabiting freshwater environments display cranial phenotypes that are intermediate between the putative ancestral form and the low-plated freshwater populations that demonstrate substantial divergence toward new phenotypic optima. We further test the hypothesis that phenotypic covariance structure is disrupted in the context of such putatively recent adaptive events. We report significant phenotypic differences among all four populations that includes a component of sexual dimorphism. Furthermore, we show evidence of disrupted phenotypic covariance structure among these populations. Taken together, these findings indicate the importance of phenotypic quantification as a key step in elucidating both the ecological processes responsible for rapid adaptive radiations and the role of developmental mechanisms in biasing evolutionary change.

Isometric Scaling in Developing Long Bones Is Achieved by an Optimal Epiphyseal Growth Balance

One of the major challenges that developing organs face is scaling, that is, the adjustment of physical proportions during the massive increase in size. Although organ scaling is fundamental for development and function, little is known about the mechanisms that regulate it. Bone superstructures are projections that typically serve for tendon and ligament insertion or articulation and, therefore, their position along the bone is crucial for musculoskeletal functionality. As bones are rigid structures that elongate only from their ends, it is unclear how superstructure positions are regulated during growth to end up in the right locations. Here, we document the process of longitudinal scaling in developing mouse long bones and uncover the mechanism that regulates it. To that end, we performed a computational analysis of hundreds of three-dimensional micro-CT images, using a newly developed method for recovering the morphogenetic sequence of developing bones. Strikingly, analysis revealed that the relative position of all superstructures along the bone is highly preserved during more than a 5-fold increase in length, indicating isometric scaling. It has been suggested that during development, bone superstructures are continuously reconstructed and relocated along the shaft, a process known as drift. Surprisingly, our results showed that most superstructures did not drift at all. Instead, we identified a novel mechanism for bone scaling, whereby each bone exhibits a specific and unique balance between proximal and distal growth rates, which accurately maintains the relative position of its superstructures. Moreover, we show mathematically that this mechanism minimizes the cumulative drift of all superstructures, thereby optimizing the scaling process. Our study reveals a general mechanism for the scaling of developing bones. More broadly, these findings suggest an evolutionary mechanism that facilitates variability in bone morphology by controlling the activity of individual epiphyseal plates.

A novel computational approach for studying bone morphogenesis reveals that the longitudinal proportions of developing long bones are accurately maintained throughout elongation by the balance between proximal and distal growth rates.

One of the major challenges that developing organs face is scaling, that is, the adjustment of physical proportions during the massive increase in size. Bone superstructures are projections that typically serve for tendon and ligament insertion or articulation. Therefore, superstructure position along the bone is crucial for musculoskeletal functionality. As bones are rigid structures that elongate only from their ends, it is unclear how superstructure positions are regulated during growth to end up in the right locations.

Here, by analyzing a massive database of micro-CT images of developing mouse long bones, we show that all superstructures maintain their relative positions throughout development. It has been suggested that during development, superstructures are continuously reconstructed and relocated along the shaft, a process known as drift. However, our analysis reveals that most superstructures did not drift at all, implying the involvement of another mechanism. Indeed, we identify a novel mechanism for bone scaling, whereby each bone exhibits a specific and unique balance between the growth rates from its two ends, which accurately maintains the relative position of its superstructures. Moreover, we show mathematically that this mechanism minimizes the cumulative drift of all superstructures, thereby optimizing the scaling process.

Mind the gap: genetic manipulation of basicranial growth within synchondroses modulates calvarial and facial shape in mice through epigenetic interactions

Phenotypic integration patterns in the mammalian skull have long been a focus of intense interest as a result of their suspected influence on the trajectory of hominid evolution. Here we test the hypothesis that perturbation of cartilage growth, which directly affects only the chondrocranium during development, will produce coordinated shape changes in the adult calvarium and face regardless of mechanism. Using two murine models of cartilage undergrowth that target two very different mechanisms, we show that strong reduction in cartilage growth produces a short, wide, and more flexed cranial base. This in turn produces a short, wide face in both models. Cranial base and face are already correlated early in ontogeny, and the relationship between these modules gains structure through postnatal growth and development. These results provide further evidence that there exist physical interactions between developing parts of the phenotype that produce variation at a distance from the actual locus upon which a particular selective pressure is acting. Phenotypic changes observed over the course of evolution may not all require adaptationist explanations; rather, it is likely that a substantial portion of observed phenotypic variation over the history of a clade is not directly adaptive but rather a secondary consequence of some local response to selection.

Two developmentally temporal quantitative trait loci underlie convergent evolution of increased branchial bone length in sticklebacks

In convergent evolution, similar phenotypes evolve repeatedly in independent populations, often reflecting adaptation to similar environments. Understanding whether convergent evolution proceeds via similar or different genetic and developmental mechanisms offers insight towards the repeatability and predictability of evolution. Oceanic populations of threespine stickleback fish, Gasterosteus aculeatus, have repeatedly colonized countless freshwater lakes and streams, where new diets lead to morphological adaptations related to feeding. Here, we show that heritable increases in branchial bone length have convergently evolved in two independently derived freshwater stickleback populations. In both populations, an increased bone growth rate in juveniles underlies the convergent adult phenotype, and one population also has a longer cartilage template. Using F2 crosses from these two freshwater populations, we show that two quantitative trait loci (QTL) control branchial bone length at distinct points in development. In both populations, a QTL on chromosome 21 controls bone length throughout juvenile development, and a QTL on chromosome 4 controls bone length only in adults. In addition to these similar developmental profiles, these QTL show similar chromosomal locations in both populations. Our results suggest that sticklebacks have convergently evolved longer branchial bones using similar genetic and developmental programmes in two independently derived populations.

Environmental temperature impact on bone and cartilage growth

Environmental temperature can have a surprising impact on extremity growth in homeotherms, but the underlying mechanisms have remained elusive for over a century. Limbs of animals raised at warm ambient temperature are significantly and permanently longer than those of littermates housed at cooler temperature. These remarkably consistent lab results closely resemble the ecogeographical tenet described by Allen's "extremity size rule," that appendage length correlates with temperature and latitude. This phenotypic growth plasticity could have adaptive significance for thermal physiology. Shortened extremities help retain body heat in cold environments by decreasing surface area for potential heat loss. Homeotherms have evolved complex mechanisms to maintain tightly regulated internal temperatures in challenging environments, including "facultative extremity heterothermy" in which limb temperatures can parallel ambient. Environmental modulation of tissue temperature can have direct and immediate consequences on cell proliferation, metabolism, matrix production, and mineralization in cartilage. Temperature can also indirectly influence cartilage growth by modulating circulating levels and delivery routes of essential hormones and paracrine regulators. Using an integrated approach, this article synthesizes classic studies with new data that shed light on the basis and significance of this enigmatic growth phenomenon and its relevance for treating human bone elongation disorders. Discussion centers on the vasculature as a gateway to understanding the complex interconnection between direct (local) and indirect (systemic) mechanisms of temperature-enhanced bone lengthening. Recent advances in imaging modalities that enable the dynamic study of cartilage growth plates in vivo will be key to elucidating fundamental physiological mechanisms of long bone growth regulation.

Mapping the natural variation in whole bone stiffness and strength across skeletal sites

Traits of the skeletal system are coordinately adjusted to establish mechanical homeostasis in response to genetic and environmental factors. Prior work demonstrated that this 'complex adaptive' process is not perfect, revealing a two-fold difference in whole bone stiffness of the tibia across a population. Robustness (specifically, total cross-sectional area relative to length) varies widely across skeletal sites and between sexes. However, it is unknown whether the natural variation in whole bone stiffness and strength also varies across skeletal sites and between men and women. We tested the hypotheses that: 1) all major long bones of the appendicular skeleton demonstrate inherent, systemic constraints in the degree to which morphological and compositional traits can be adjusted for a given robustness; and 2) these traits covary in a predictable manner independent of body size and robustness. We assessed the functional relationships among robustness, cortical area (Ct.Ar), cortical tissue mineral density (Ct.TMD), and bone strength index (BSI) across the long bones of the upper and lower limbs of 115 adult men and women. All bones showed a significant (p<0.001) positive regression between BSI and robustness after adjusting for body size, with slender bones being 1.7-2.3 times less stiff and strong in men and 1.3-2.8 times less stiff and strong in women compared to robust bones. Our findings are the first to document the natural inter-individual variation in whole bone stiffness and strength that exist within populations and that is predictable based on skeletal robustness for all major long bones. Documenting and further understanding this natural variation in strength may be critical for differentially diagnosing and treating skeletal fragility.

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009; Meyer 2011; Nosil 2012). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009), or pollinator-mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. (2014) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early-stage divergence in Midas cichlids.

Genomic correlates of relationship QTL involved in fore- versus hind limb divergence in mice

Divergence of serially homologous elements of organisms is a common evolutionary pattern contributing to increased phenotypic complexity. Here, we study the genomic intervals affecting the variational independence of fore- and hind limb traits within an experimental mouse population. We use an advanced intercross of inbred mouse strains to map the loci associated with the degree of autonomy between fore- and hind limb long bone lengths (loci affecting the relationship between traits, relationship quantitative trait loci [rQTL]). These loci have been proposed to interact locally with the products of pleiotropic genes, thereby freeing the local trait from the variational constraint due to pleiotropic mutations. Using the known polymorphisms (single nucleotide polymorphisms [SNPs]) between the parental strains, we characterized and compared the genomic regions in which the rQTL, as well as their interaction partners (intQTL), reside. We find that these two classes of QTL intervals harbor different kinds of molecular variation. SNPs in rQTL intervals more frequently reside in limb-specific cis-regulatory regions than SNPs in intQTL intervals. The intQTL loci modified by the rQTL, in contrast, show the signature of protein-coding variation. This result is consistent with the widely accepted view that protein-coding mutations have broader pleiotropic effects than cis-regulatory polymorphisms. For both types of QTL intervals, the underlying candidate genes are enriched for genes involved in protein binding. This finding suggests that rQTL effects are caused by local interactions among the products of the causal genes harbored in rQTL and intQTL intervals. This is the first study to systematically document the population-level molecular variation underlying the evolution of character individuation.

Deficiency of zebrafish fgf20a results in aberrant skull remodeling that mimics both human cranial disease and evolutionarily important fish skull morphologies

The processes that direct skull remodeling are of interest to both human-oriented studies of cranial dysplasia and evolutionary studies of skull divergence. There is increasing awareness that these two fields can be mutually informative when natural variation mimics pathology. Here we describe a zebrafish mutant line, devoid of blastema (dob), which does not have a functional fgf20a protein, and which also presents cranial defects similar to both adaptive and clinical variation. We used geometric morphometric methods to provide quantitative descriptions of the effects of the dob mutation on skull morphogenesis. In combination with "whole-mount in situ hybridization" labeling of normal fgf20a expression and assays for osteoblast and osteoclast activity, the results of these analyses indicate that cranial dysmorphologies in dob zebrafish are generated by aberrations in post-embryonic skull remodeling via decreased osteoblasotgenesis and increased osteoclastogenesis. Mutational effects include altered skull vault geometries and midfacial hypoplasia that are consistent with key diagnostic signs for multiple human craniofacial syndromes. These phenotypic shifts also mimic changes in the functional morphology of fish skulls that have arisen repeatedly in several highly successful radiations (e.g., damselfishes and East-African rift-lake cichlids). Our results offer the dob/fgf20a mutant as an experimentally tractable model with which to examine post-embryonic skull development as it relates to human disease and vertebrate evolution.

Macroevolutionary diversity of amniote limb proportions predicted by developmental interactions

Mammals, birds, and reptiles exhibit a remarkable diversity of limb proportions. These evolved differences are thought to reflect selection for biomechanical, postural, and locomotor requirements primarily acting on independent variation in later fetal and postnatal segmental growth. However, earlier conserved developmental events also have the potential to impact the evolvability of limb proportions by limiting or biasing initial variation among segments. Notably, proximo-distal patterning of the amniote limb through activation-inhibition dynamics predicts that initial proportions of segments should exhibit both tradeoffs between stylopod and autopod and a diagnostic reduction in variance of the zeugopod. Here it is demonstrated that this developmental "design rule" predicts patterns of macroevolutionary diversity despite the effects of variation in segmental growth over ontogeny, lineage-specific differences in phylogenetic history, or functional adaptation. These results provide critical comparative evidence of a conserved Turing-like mechanism in proximo-distal limb segmentation, and suggest that development has played a previously unrecognized role in the evolvability of limb proportions in a wide range of amniote taxa.

Convergent evolution of sexual dimorphism in skull shape using distinct developmental strategies

Studies integrating evolutionary and developmental analyses of morphological variation are of growing interest to biologists as they promise to shed fresh light on the mechanisms of morphological diversification. Sexually dimorphic traits tend to be incredibly divergent across taxa. Such diversification must arise through evolutionary modifications to sex differences during development. Nevertheless, few studies of dimorphism have attempted to synthesize evolutionary and developmental perspectives. Using geometric morphometric analysis of head shape for 50 Anolis species, we show that two clades have converged on extreme levels of sexual dimorphism through similar, male-specific changes in facial morphology. In both clades, males have evolved highly elongate faces whereas females retain faces of more moderate proportion. This convergence is accomplished using distinct developmental mechanisms; one clade evolved extreme dimorphism through the exaggeration of a widely shared, potentially ancestral, developmental strategy whereas the other clade evolved a novel developmental strategy not observed elsewhere in the genus. Together, our analyses indicate that both shared and derived features of development contribute to macroevolutionary patterns of morphological diversity among Anolis lizards.

Heritability of articular cartilage regeneration and its association with ear wound healing in mice

OBJECTIVE

Emerging evidence suggests that genetic components contribute significantly to cartilage degeneration in osteoarthritis pathophysiology, but little information is available on the genetics of cartilage regeneration. Therefore, this study was undertaken to investigate cartilage regeneration in genetic murine models using common inbred strains and a set of recombinant inbred (RI) lines generated from LG/J (healer of ear wounds) and SM/J (nonhealer) inbred mouse strains.

METHODS

An acute full-thickness cartilage injury was introduced in the trochlear groove of 8-week-old mice (n=265) through microsurgery. Mouse knee joints were sagittally sectioned and stained with toluidine blue to evaluate regeneration. For the ear wound phenotype, a bilateral 2-mm through-and-through puncture was created in 6-week-old mice (n=229), and healing outcomes were measured after 30 days. Broad-sense heritability and genetic correlations were calculated for both phenotypes.

RESULTS

Time-course analysis of the RI mouse lines showed no significant regeneration until 16 weeks after surgery; at that time, the strains could be segregated into 3 categories: good, intermediate, and poor healers. Analysis of heritability (H2) showed that both cartilage regeneration (H2=26%; P=0.006) and ear wound closure (H2=53%; P<0.00001) were significantly heritable. The genetic correlations between the two healing phenotypes for common inbred mouse strains (r=0.92) and RI mouse lines (r=0.86) were found to be extremely high.

CONCLUSION

Our findings indicate that articular cartilage regeneration in mice is heritable, the differences between the mouse lines are due to genetic differences, and a strong genetic correlation between the two phenotypes exists, indicating that they plausibly share a common genetic basis. We therefore surmise that LG/J by SM/J intercross mice can be used to dissect the genetic basis of variation in cartilage regeneration.

The developmental renaissance in adaptationism

From an adaptation perspective, unoccupied patches of morphological space are inferred to be empty because they are of low fitness and selected against. These inferences hinge on venturesome assumptions, because emptiness is explained by low fitness and low fitness is inferred from emptiness. Moreover, non-adaptive factors, such as developmental constraint, could also plausibly account for empty morphospace. In response, biologists increasingly study ontogeny to test the assumption that unobserved phenotypes could be produced if selection were to favor them; finding that empty space morphologies can be readily produced in development helps reject constraint and lends support to adaptive hypotheses. This developmental approach to adaptation calls on manifold techniques, including embryology, artificial selection and comparative methods. Belying their diversity, all of these methods examine the causes of empty morphospace and mark a return of development, long excluded from traditional evolutionary biology, to adaptationist practice.

Repeated modification of early limb morphogenesis programmes underlies the convergence of relative limb length in Anolis lizards

The independent evolution of similar morphologies has long been a subject of considerable interest to biologists. Does phenotypic convergence reflect the primacy of natural selection, or does development set the course of evolution by channelling variation in certain directions? Here, we examine the ontogenetic origins of relative limb length variation among Anolis lizard habitat specialists to address whether convergent phenotypes have arisen through convergent developmental trajectories. Despite the numerous developmental processes that could potentially contribute to variation in adult limb length, our analyses reveal that, in Anolis lizards, such variation is repeatedly the result of changes occurring very early in development, prior to formation of the cartilaginous long bone anlagen.

Genetic factors and diet affect long-bone length in the F34 LG,SM advanced intercross

Previous studies on the LG,SM advanced intercross line have identified approximately 40 quantitative trait loci (QTL) for long -bone (humerus, ulna, femur, and tibia) lengths. In this study, long-bone-length QTL were fine-mapped in the F(34) generation (n = 1424) of the LG,SM advanced intercross. Environmental effects were assessed by dividing the population by sex between high-fat and low-fat diets, producing eight sex/diet cohorts. We identified 145 individual bone-length QTL comprising 45 pleiotropic QTL; 69 replicated QTL from previous studies, 35 were new traits significant at previously identified loci, and 41 were novel QTL. Many QTL affected only a subset of the population based on sex and/or diet. Eight of ten known skeletal growth genes were upregulated in 3-week-old LG/J male proximal tibial growth plates relative to SM/J. The sequences of parental strains LG/J and SM/J indicated the presence of over half a million polymorphisms in the confidence intervals of these 45 QTL. We examined 526 polymorphisms and found that 97 represented radical changes to amino acid composition while 40 were predicted to be deleterious to protein function. Additional experimentation is required to understand how changes in gene regulation or protein function can alter the genetic architecture and interact with the environment to produce phenotypic variation.