Phenotypic integration of skeletal traits during growth buffers genetic variants affecting the slenderness of femora in inbred mouse strains

Mitofusin 2 plays a critical role in maintaining the functional integrity of the neuromuscular-skeletal axis

PURPOSE

Mitofusin 2 (Mfn2) is one of two mitofusins involved in regulating mitochondrial size, shape and function, including mitophagy, an important cellular mechanism to limit oxidative stress. Reduced expression of Mfn2 has been associated with impaired osteoblast differentiation and function and a reduction in the number of viable osteocytes in bone. We hypothesized that the genetic absence of Mfn2 in these cells would increase their susceptibility to aging-associated metabolic stress, leading to a progressive impairment in skeletal homeostasis over time.

METHODS

Mfn2 was selectively deleted in vivo at three different stages of osteoblast lineage commitment by crossing mice in which the Mfn2 gene was floxed with transgenic mice expressing Cre under the control of the promoter for Osterix (OSX), collagen1a1, or DMP1 (Dentin Matrix Acidic Phosphoprotein 1).

RESULTS

Mice in which Mfn2 was deleted using DMP1-cre demonstrated a progressive and dramatic decline in bone mineral density (BMD) beginning at 10 weeks of age (n = 5 for each sex and each genotype from age 10 to 20 weeks). By 15 weeks, there was evidence for a functional decline in muscle performance as assessed using a rotarod apparatus (n = 3; 2 males/ 1 female for each genotype), accompanied by a decline in lean body mass. A marked reduction in trabecular bone mass was evident on bone histomorphometry, and biomechanical testing at 25 weeks (k/o: 2 male/1 female, control 2 male/2 female) revealed severely impaired femur strength. Extensive regional myofiber atrophy and degeneration was observed on skeletal muscle histology. Electron microscopy showed progressive disruption of cellular architecture, with disorganized sarcomeres and a bloated mitochondrial reticulum. There was also evidence of neurodegeneration within the ventral horn and roots of the lumbar spinal cord, which was accompanied by myelin loss and myofiber atrophy. Deletion of Mfn2 using OSX-cre or Col1a1-cre did not result in a musculoskeletal phenotype. Where possible, male and female animals were analyzed separately, but small numbers of animals in each group limited statistical power. For other outcomes, where sex was not considered, small sample sizes might still limit the strength of the observation.

CONCLUSION

Despite known functional overlap of Mfn1 and Mfn2 in some tissues, and their co-expression in bone, muscle and spinal cord, deletion of Mfn2 using the 8 kB DMP1 promoter uncovered an important non-redundant role for Mfn2 in maintaining the neuromuscular/bone axis.

Effects of long-term voluntary wheel running and selective breeding for wheel running on femoral nutrient canals

The nutrient artery provides ~50%-70% of the total blood volume to long bones in mammals. Studying the functional characteristics of this artery in vivo can be difficult and expensive, so most researchers have measured the nutrient foramen, an opening on the outer surface of the bone that served as the entry point for the nutrient artery during development and bone ossification. Others have measured the nutrient canal (i.e., the passage which the nutrient artery once occupied), given that the external dimensions of the foramen do not necessarily remain uniform from the periosteal surface to the medullary cavity. The nutrient canal, as an indicator of blood flow to long bones, has been proposed to provide a link to studying organismal activity (e.g., locomotor behavior) from skeletal morphology. However, although external loading from movement and activity causes skeletal remodeling, it is unclear whether it affects the size or configuration of nutrient canals. To investigate whether nutrient canals can exhibit phenotypic plasticity in response to physical activity, we studied a mouse model in which four replicate high runner (HR) lines have been selectively bred for high voluntary wheel-running behavior. The selection criterion is the average number of wheel revolutions on days 5 and 6 of a 6-day period of wheel access as young adults (~6-8 weeks old). An additional four lines are bred without selection to serve as controls (C). For this study, 100 female mice (half HR, half C) from generation 57 were split into an active group housed with wheels and a sedentary group housed without wheels for 12 weeks starting at ~24 days of age. Femurs were collected, soft tissues were removed, and femora were micro-computed tomography scanned at a resolution of 12 μm. We then imported these scans into AMIRA and created 3D models of femoral nutrient canals. We tested for evolved differences in various nutrient canal traits between HR and C mice, plastic changes resulting from chronic exercise, and the selection history-by-exercise interaction. We found few differences between the nutrient canals of HR versus C mice, or between the active and sedentary groups. We did find an interaction between selection history and voluntary exercise for the total number of nutrient canals per femur, in which wheel access increased the number of canals in C mice but decreased it in HR mice. Our results do not match those from an earlier study, conducted at generation 11, which was prior to the HR lines reaching selection limits for wheel running. The previous study found that mice from the HR lines had significantly larger total canal cross-sectional areas compared to those from C lines. However, this discrepancy is consistent with studies of other skeletal traits, which have found differences between HR and C mice to be somewhat inconsistent across generations, including the loss of some apparent adaptations with continued selective breeding after reaching a selection limit for wheel-running behavior.

Lifelong Excess in GH Elicits Sexually Dimorphic Effects on Skeletal Morphology and Bone Mechanical Properties

Excess in growth hormone (GH) levels, seen in patients with acromegaly, is associated with increases in fractures. This happens despite wider bones and independent of bone mineral density. We used the bovine GH (bGH) transgenic mice, which show constitutive excess in GH and insulin-like growth factor 1 (IGF-1) in serum and tissues, to study how lifelong increases in GH and IGF-1 affect skeletal integrity. Additionally, we crossed the acid labile subunit (ALS) null (ALSKO) to the bGH mice to reduce serum IGF-1 levels. Our findings indicate sexually dimorphic effects of GH on cortical and trabecular bone. Male bGH mice showed enlarged cortical diameters, but with marrow cavity expansion and thin cortices as well as increased vascular porosity that were associated with reductions in diaphyseal strength and stiffness. In contrast, female bGH mice presented with significantly smaller-diameter diaphysis, with greater cortical bone thickness and with a slightly reduced tissue elastic modulus (by microindentation), ultimately resulting in overall stronger, stiffer bones. We found increases in C-terminal telopeptide of type 1 collagen and procollagen type 1 N propeptide in serum, independent of circulating IGF-1 levels, indicating increased bone remodeling with excess GH. Sexual dimorphism in response to excess GH was also observed in the trabecular bone compartment, particularly at the femur distal metaphysis. Female bGH mice preserved their trabecular architecture during aging, whereas trabecular bone volume in male bGH mice significantly reduced and was associated with thinning of the trabeculae. We conclude that pathological excess in GH results in sexually dimorphic changes in bone architecture and gains in bone mass that affect whole-bone mechanical properties, as well as sex-specific differences in bone material properties. © 2022 American Society for Bone and Mineral Research (ASBMR).

Components of the Gut Microbiome That Influence Bone Tissue-Level Strength

Modifications to the constituents of the gut microbiome influence bone density and tissue-level strength, but the specific microbial components that influence tissue-level strength in bone are not known. Here, we selectively modify constituents of the gut microbiota using narrow-spectrum antibiotics to identify components of the microbiome associated with changes in bone mechanical and material properties. Male C57BL/6J mice (4 weeks) were divided into seven groups (n = 7-10/group) and had taxa within the gut microbiome removed through dosing with: (i) ampicillin; (ii) neomycin; (iii) vancomycin; (iv) metronidazole; (v) a cocktail of all four antibiotics together (with zero-calorie sweetener to ensure intake); (vi) zero-calorie sweetener only; or (vii) no additive (untreated) for 12 weeks. Individual antibiotics remove only some taxa from the gut, while the cocktail of all four removes almost all microbes. After accounting for differences in geometry, whole bone strength was reduced in animals with gut microbiome modified by neomycin (-28%, p = 0.002) and was increased in the group in which the gut microbiome was altered by sweetener alone (+39%, p < 0.001). Analysis of the fecal microbiota detected seven lower-ranked taxa differentially abundant in animals with impaired tissue-level strength and 14 differentially abundant taxa associated with increased tissue-level strength. Histological and serum markers of bone turnover and trabecular bone volume per tissue volume (BV/TV) did not differ among groups. These findings demonstrate that modifications to the taxonomic components of the gut microbiome have the potential to decrease or increase tissue-level strength of bone independent of bone quantity and without noticeable changes in bone turnover. © 2021 American Society for Bone and Mineral Research (ASBMR).

Articular Cartilage Friction, Strain, and Viability Under Physiological to Pathological Benchtop Sliding Conditions

In vivo, articular cartilage is exceptionally resistant to wear, damage, and dysfunction. However, replicating cartilage's phenomenal in vivo tribomechanics (i.e., high fluid load support, low frictions and strains) and mechanobiology on the benchtop has been difficult, because classical testing approaches tend to minimize hydrodynamic contributors to tissue function. Our convergent stationary contact area (cSCA) configuration retains the ability for hydrodynamically-mediated processes to contribute to interstitial hydration recovery and tribomechanical function via 'tribological rehydration'. Using the cSCA, we investigated how in situ chondrocyte survival is impacted by the presence of tribological rehydration during the reciprocal sliding of a glass counterface against a compressively loaded equine cSCA cartilage explant. When tribological rehydration was compromised during testing, by slow-speed sliding, 'pathophysiological' tribomechanical environments and high surface cell death were observed. When tribological rehydration was preserved, by high-speed sliding, 'semi-physiological' sliding environments and suppressed cell death were realized. Inclusion of synovial fluid during testing fostered 'truly physiological' sliding outcomes consistent with the in vivo environment but had limited influence on cell death compared to high-speed sliding in PBS. Subsequently, path analysis identified friction as a primary driver of cell death, with strain an indirect driver, supporting the contention that articulation mediated rehydration can benefit both the biomechanical properties and biological homeostasis of cartilage.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12195-021-00671-2.

External Bone Size Is a Key Determinant of Strength-Decline Trajectories of Aging Male Radii

Given prior work showing associations between remodeling and external bone size, we tested the hypothesis that wide bones would show a greater negative correlation between whole-bone strength and age compared with narrow bones. Cadaveric male radii (n = 37 pairs, 18 to 89 years old) were evaluated biomechanically, and samples were sorted into narrow and wide subgroups using height-adjusted robustness (total area/bone length). Strength was 54% greater (p < 0.0001) in wide compared with narrow radii for young adults (<40 years old). However, the greater strength of young-adult wide radii was not observed for older wide radii, as the wide (R² = 0.565, p = 0.001), but not narrow (R² = 0.0004, p = 0.944) subgroup showed a significant negative correlation between strength and age. Significant positive correlations between age and robustness (R² = 0.269, p = 0.048), cortical area (Ct.Ar; R² = 0.356, p = 0.019), and the mineral/matrix ratio (MMR; R² = 0.293, p = 0.037) were observed for narrow, but not wide radii (robustness: R² = 0.015, p = 0.217; Ct.Ar: R² = 0.095, p = 0.245; MMR: R² = 0.086, p = 0.271). Porosity increased with age for the narrow (R² = 0.556, p = 0.001) and wide (R² = 0.321, p = 0.022) subgroups. The wide subgroup (p < 0.0001) showed a significantly greater elevation of a new measure called the Cortical Pore Score, which quantifies the cumulative effect of pore size and location, indicating that porosity had a more deleterious effect on strength for wide compared with narrow radii. Thus, the divergent strength-age regressions implied that narrow radii maintained a low strength with aging by increasing external size and mineral content to mechanically offset increases in porosity. In contrast, the significant negative strength-age correlation for wide radii implied that the deleterious effect of greater porosity further from the centroid was not offset by changes in outer bone size or mineral content. Thus, the low strength of elderly male radii arose through different biomechanical mechanisms. Consideration of different strength-age regressions (trajectories) may inform clinical decisions on how best to treat individuals to reduce fracture risk. © 2019 American Society for Bone and Mineral Research.

Dissecting the Genetics of Osteoporosis using Systems Approaches

Osteoporosis is a condition characterized by low bone mineral density (BMD) and an increased risk of fracture. Traits contributing to osteoporotic fracture are highly heritable, indicating that a comprehensive understanding of bone requires a thorough understanding of the genetic basis of bone traits. Towards this goal, genome-wide association studies (GWASs) have identified over 500 loci associated with bone traits. However, few of the responsible genes have been identified, and little is known of how these genes work together to influence systems-level bone function. In this review, we describe how systems genetics approaches can be used to fill these knowledge gaps.

Ontogeny of hallucal metatarsal rigidity and shape in the rhesus monkey (Macaca mulatta) and chimpanzee (Pan troglodytes)

Life history variables including the timing of locomotor independence, along with changes in preferred locomotor behaviors and substrate use during development, influence how primates use their feet throughout ontogeny. Changes in foot function during development, in particular the nature of how the hallux is used in grasping, can lead to different structural changes in foot bones. To test this hypothesis, metatarsal midshaft rigidity [estimated from the polar second moment of area (J) scaled to bone length] and cross-sectional shape (calculated from the ratio of maximum and minimum second moments of area, Imax /Imin ) were examined in a cross-sectional ontogenetic sample of rhesus macaques (Macaca mulatta; n = 73) and common chimpanzees (Pan troglodytes; n = 79). Results show the hallucal metatarsal (Mt1) is relatively more rigid (with higher scaled J-values) in younger chimpanzees and macaques, with significant decreases in relative rigidity in both taxa until the age of achieving locomotor independence. Within each age group, Mt1 rigidity is always significantly higher in chimpanzees than macaques. When compared with the lateral metatarsals (Mt2-5), the Mt1 is relatively more rigid in both taxa and across all ages; however, this difference is significantly greater in chimpanzees. Length and J scale with negative allometry in all metatarsals and in both species (except the Mt2 of chimpanzees, which scales with positive allometry). Only in macaques does Mt1 midshaft shape significantly change across ontogeny, with older individuals having more elliptical cross-sections. Different patterns of development in metatarsal diaphyseal rigidity and shape likely reflect the different ways in which the foot, and in particular the hallux, functions across ontogeny in apes and monkeys.

Femoral Neck External Size but not aBMD Predicts Structural and Mass Changes for Women Transitioning Through Menopause

The impact of adult bone traits on changes in bone structure and mass during aging is not well understood. Having shown that intracortical remodeling correlates with external size of adult long bones led us to hypothesize that age-related changes in bone traits also depend on external bone size. We analyzed hip dual-energy X-ray absorptiometry images acquired longitudinally over 14 years for 198 midlife women transitioning through menopause. The 14-year change in bone mineral content (BMC, R²  = 0.03, p = 0.015) and bone area (R²  = 0.13, p = 0.001), but not areal bone mineral density (aBMD, R²  = 0.00, p = 0.931) correlated negatively with baseline femoral neck external size, adjusted for body size using the residuals from a linear regression between baseline bone area and height. The dependence of the 14-year changes in BMC and bone area on baseline bone area remained significant after adjusting for race/ethnicity, postmenopausal hormone use, the 14-year change in weight, and baseline aBMD, weight, height, and age. Women were sorted into tertiles using the baseline bone area-height residuals. The 14-year change in BMC (p = 0.009) and bone area (p = 0.001) but not aBMD (p = 0.788) differed across the tertiles. This suggested that women showed similar changes in aBMD for different structural and biological reasons: women with narrow femoral necks showed smaller changes in BMC but greater increases in bone area compared to women with wide femoral necks who showed greater losses in BMC but without large compensatory increases in bone area. This finding is opposite to expectations that periosteal expansion acts to mechanically offset bone loss. Thus, changes in femoral neck structure and mass during menopause vary widely among women and are predicted by baseline external bone size but not aBMD. How these different structural and mass changes affect individual strength-decline trajectories remains to be determined. © 2017 American Society for Bone and Mineral Research.

Canalization Leads to Similar Whole Bone Mechanical Function at Maturity in Two Inbred Strains of Mice

Previously, we showed that cortical mineralization is coordinately adjusted to mechanically offset external bone size differences between A/J (narrow) and C57BL/6J (wide) mouse femora to achieve whole bone strength equivalence at adulthood. The identity of the genes and their interactions that are responsible for establishing this homeostatic state (ie, canalization) remain unknown. We hypothesize that these inbred strains, whose interindividual differences in bone structure and material properties mimic that observed among humans, achieve functional homeostasis by differentially adjusting key molecular pathways regulating external bone size and mineralization throughout growth. The cortices of A/J and C57BL/6J male mouse femora were phenotyped and gene expression levels were assessed across growth (ie, ages 2, 4, 6, 8, 12, 16 weeks). A difference in total cross-sectional area (p < 0.01) and cortical tissue mineral density were apparent between mouse strains by age 2 weeks and maintained at adulthood (p < 0.01). These phenotypic dissimilarities corresponded to gene expression level differences among key regulatory pathways throughout growth. A/J mice had a 1.55- to 7.65-fold greater expression among genes inhibitory to Wnt pathway induction, whereas genes involved in cortical mineralization were largely upregulated 1.50- to 3.77-fold to compensate for their narrow diaphysis. Additionally, both mouse strains showed an upregulation among Wnt pathway antagonists corresponding to the onset of adult ambulation (ie, increased physiological loads). This contrasts with other studies showing an increase in Wnt pathway activation after functionally isolated, experimental in vivo loading regimens. A/J and C57BL/6J long bones provide a model to develop a systems-based approach to identify individual genes and the gene-gene interactions that contribute to trait differences between the strains while being involved in the process by which these traits are coordinately adjusted to establish similar levels of mechanical function, thus providing insight into the process of canalization. © 2017 American Society for Bone and Mineral Research.

The development of inter-strain variation in cortical and trabecular traits during growth of the mouse lumbar vertebral body

How cortical and trabecular bone co-develop to establish a mechanically functional structure is not well understood. Comparing early postnatal differences in morphology of lumbar vertebral bodies for three inbred mouse strains identified coordinated changes within and between cortical and trabecular traits. These early coordinate changes defined the phenotypic differences among the inbred mouse strains.

INTRODUCTION

Age-related changes in cortical and trabecular traits have been well studied; however, very little is known about how these bone tissues co-develop from day 1 of postnatal growth to establish functional structures by adulthood. In this study, we aimed to establish how cortical and trabecular tissues within the lumbar vertebral body change during growth for three inbred mouse strains that express wide variation in adult bone structure and function.

METHODS

Bone traits were quantified for lumbar vertebral bodies of female A/J, C57BL/6J (B6), and C3H/HeJ (C3H) inbred mouse strains from 1 to 105 days of age (n = 6-10 mice/age/strain).

RESULTS

Inter-strain differences in external bone size were observed as early as 1 day of age. Reciprocal and rapid changes in the trabecular bone volume fraction and alignment in the direction of axial compression were observed by 7 days of age. Importantly, the inter-strain difference in adult trabecular bone volume fraction was established by 7 days of age. Early variation in external bone size and trabecular architecture was followed by progressive increases in cortical area between 28 and 105 days of age, with the greatest increases in cortical area seen in the mouse strain with the lowest trabecular mass.

CONCLUSION

Establishing the temporal changes in bone morphology for three inbred mouse strains revealed that genetic variation in adult trabecular traits were established early in postnatal development. Early variation in trabecular architecture preceded strain-specific increases in cortical area and changes in cortical thickness. This study established the sequence of how cortical and trabecular traits co-develop during growth, which is important for identifying critical early ages to further focus on intervention studies that optimize adult bone strength.

HBM Mice Have Altered Bone Matrix Composition and Improved Material Toughness

The G171V mutation in the low-density lipoprotein receptor-related protein 5 (LRP5) leads to a high bone mass (HBM) phenotype. Studies using HBM transgenic mouse models have consistently found increased bone mass and whole-bone strength, but little attention has been paid to the composition of the bone matrix. The current study sought to determine if the cortical bone matrix composition differs in HBM and wild-type mice and to determine how much of the variance in bone material properties is explained by variance in matrix composition. Consistent with previous studies, HBM mice had greater cortical area, moment of inertia, ultimate force, bending stiffness, and energy to failure than wild-type animals. The increased energy to failure was primarily caused by a large increase in post-yield behavior, with no difference in pre-yield behavior. The HBM mice had increased mineral-to-matrix and collagen cross-link ratios, and decreased crystallinity, carbonate, and acid phosphate substitution as measured by Fourier transform infrared microspectroscopy, but no differences in crystal length, intra-fibular strains, and mineral spacing compared to wild-type controls, as measured by X-ray scattering. The largest between genotype difference in material properties was a twofold increase in the modulus of toughness in HBM mice. Step-wise regression analyses showed that the specific matrix compositional parameters most closely associated with material properties varied between the wild-type and HBM genotypes. Although the mechanisms controlling the paradoxical combination of more mineralized yet tougher bone in HBM mice remain to be fully explained, the findings suggest that LRP5 represents a target to not only build bone mass but also to improve bone quality.

Genetic regulation of bone strength: a review of animal model studies

Population- and family-based studies have established that fragility fracture risk is heritable; yet, the genome-wide association studies published to date have only accounted for a small fraction of the known variation for fracture risk of either the femur or the lumbar spine. Much work has been carried out using animal models toward finding genetic loci that are associated with bone strength. Studies using animal models overcome some of the issues associated with using patient data, but caution is needed when interpreting the results. In this review, we examine the types of tests that have been used for forward genetics mapping in animal models to identify loci and/or genes that regulate bone strength and discuss the limitations of these test methods. In addition, we present a summary of the quantitative trait loci that have been mapped for bone strength in mice, rats and chickens. The majority of these loci co-map with loci for bone size and/or geometry and thus likely dictate strength via modulating bone size. Differences in bone matrix composition have been demonstrated when comparing inbred strains of mice, and these matrix differences may be associated with differences in bone strength. However, additional work is needed to identify loci that act on bone strength at the materials level.

Establishing biomechanical mechanisms in mouse models: practical guidelines for systematically evaluating phenotypic changes in the diaphyses of long bones

Mice are widely used in studies of skeletal biology, and assessment of their bones by mechanical testing is a critical step when evaluating the functional effects of an experimental perturbation. For example, a gene knockout may target a pathway important in bone formation and result in a "low bone mass" phenotype. But how well does the skeleton bear functional loads; eg, how much do bones deform during loading and how resistant are bones to fracture? By systematic evaluation of bone morphological, densitometric, and mechanical properties, investigators can establish the "biomechanical mechanisms" whereby an experimental perturbation alters whole-bone mechanical function. The goal of this review is to clarify these biomechanical mechanisms and to make recommendations for systematically evaluating phenotypic changes in mouse bones, with a focus on long-bone diaphyses and cortical bone. Further, minimum reportable standards for testing conditions and outcome variables are suggested that will improve the comparison of data across studies. Basic biomechanical principles are reviewed, followed by a description of the cross-sectional morphological properties that best inform the net cellular effects of a given experimental perturbation and are most relevant to biomechanical function. Although morphology is critical, whole-bone mechanical properties can only be determined accurately by a mechanical test. The functional importance of stiffness, maximum load, postyield displacement, and work-to-fracture are reviewed. Because bone and body size are often strongly related, strategies to adjust whole-bone properties for body mass are detailed. Finally, a comprehensive framework is presented using real data, and several examples from the literature are reviewed to illustrate how to synthesize morphological, tissue-level, and whole-bone mechanical properties of mouse long bones.

Characterization of complex, co-adapted skeletal biomechanics phenotypes: a needed paradigm shift in the genetics of bone structure and function

The genetic architecture of skeletal biomechanical performance has tremendous potential to advance our knowledge of the biological mechanisms that drive variation in skeletal fragility and osteoporosis risk. Research using traditional approaches that focus on specific gene pathways is increasing our understanding of how and to what degree those pathways may affect population-level variation in fracture susceptibility, and shows that known pathways may affect bone fragility through unsuspected mechanisms. Non-traditional approaches that incorporate a new appreciation for the degree to which bone traits co-adapt to functional loading environments, using a wide variety of redundant compensatory mechanisms to meet both physiological and mechanical demands, represent a radical departure from the dominant reductionist paradigm and have the potential to rapidly advance our understanding of bone fragility and identification of new targets for therapeutic intervention.

Altered lacunar and vascular porosity in osteogenesis imperfecta mouse bone as revealed by synchrotron tomography contributes to bone fragility

Osteogenesis imperfecta (brittle bone disease) is caused by mutations in the collagen genes and results in skeletal fragility. Changes in bone porosity at the tissue level indicate changes in bone metabolism and alter bone mechanical integrity. We investigated the cortical bone tissue porosity of a mouse model of the disease, oim, in comparison to a wild type (WT-C57BL/6), and examined the influence of canal architecture on bone mechanical performance. High-resolution 3D representations of the posterior tibial and the lateral humeral mid-diaphysis of the bones were acquired for both mouse groups using synchrotron radiation-based computed tomography at a nominal resolution of 700nm. Volumetric morphometric indices were determined for cortical bone, canal network and osteocyte lacunae. The influence of canal porosity architecture on bone mechanics was investigated using microarchitectural finite element (μFE) models of the cortical bone. Bright-field microscopy of stained sections was used to determine if canals were vascular. Although total cortical porosity was comparable between oim and WT bone, oim bone had more numerous and more branched canals (p<0.001), and more osteocyte lacunae per unit volume compared to WT (p<0.001). Lacunae in oim were more spherical in shape compared to the ellipsoidal WT lacunae (p<0.001). Histology revealed blood vessels in all WT and oim canals. μFE models of cortical bone revealed that small and branched canals, typical of oim bone, increase the risk of bone failure. These results portray a state of compromised bone quality in oim bone at the tissue level, which contributes to its deficient mechanical properties.

ENU mutation mapped to a distal region of chromosome 11 is a major determinant of bone size

Using a phenotype driven n-ethyl-nitrosourea (ENU) screen in growth hormone-deficient mice, we have identified a mutant (named 14104) that exhibited a smaller bone size. Phenotype measurements by microcomputed tomography revealed that mutant mice exhibited a 43 and 34% reduction in tissue area and bone area, respectively at the femur middiaphysis. Dynamic histomorphometry revealed a 30 and 15% lower bone formation rate at the periosteal and endosteal surface, respectively. Breaking strength of the femur was reduced by 30% in the mutant mice. To determine if the 14104 locus is involved in a mechanical loading signaling pathway, the skeletal anabolic response to tibia axial loading was evaluated. The increase in trabecular response in the loaded region was severely compromised by the 14014 mutation. To identify the location of mutation, we performed linkage analysis using 62 polymorphic markers in the B6-DBA/2J F2 mice. The genome-wide linkage analysis identified the location of the mutation to a 72 to 83 cM region on chromosome 11 with peak logarithm of the odds scores of 15 for periosteal circumference at marker D11mit338. Sequence analysis revealed no mutation in the coding region of 11 potential candidate genes. Based on these data and published data on the skeletal phenotype of genes in this region, we concluded that the 109-119 Mb region of chromosome 11 harbors a bone size gene that regulates periosteal bone formation. The mutant strain developed in this study provides an important tool to identify a novel mechanosensitive gene that determines bone size during postnatal development.

Intracortical bone remodeling variation shows strong genetic effects

Intracortical microstructure influences crack propagation and arrest within bone cortex. Genetic variation in intracortical remodeling may contribute to mechanical integrity and, therefore, fracture risk. Our aim was to determine the degree to which normal population-level variation in intracortical microstructure is due to genetic variation. We examined right femurs from 101 baboons (74 females, 27 males; aged 7-33 years) from a single, extended pedigree to determine osteon number, osteon area (On.Ar), haversian canal area, osteon population density, percent osteonal bone (%On.B), wall thickness (W.Th), and cortical porosity (Ct.Po). Through evaluation of the covariance in intracortical properties between pairs of relatives, we quantified the contribution of additive genetic effects (heritability [h (2)]) to variation in these traits using a variance decomposition approach. Significant age and sex effects account for 9 % (Ct.Po) to 21 % (W.Th) of intracortical microstructural variation. After accounting for age and sex, significant genetic effects are evident for On.Ar (h (2) = 0.79, p = 0.002), %On.B (h (2) = 0.82, p = 0.003), and W.Th (h (2) = 0.61, p = 0.013), indicating that 61-82 % of the residual variation (after accounting for age and sex effects) is due to additive genetic effects. This corresponds to 48-75 % of the total phenotypic variance. Our results demonstrate that normal, population-level variation in cortical microstructure is significantly influenced by genes. As a critical mediator of crack behavior in bone cortex, intracortical microstructural variation provides another mechanism through which genetic variation may affect fracture risk.

Functional integration of skeletal traits: an intraskeletal assessment of bone size, mineralization, and volume covariance

Understanding the functional integration of skeletal traits and how they naturally vary within and across populations will benefit assessments of functional adaptation directed towards interpreting bone stiffness in contemporary and past humans. Moreover, investigating how these traits intraskeletally vary will guide us closer towards predicting fragility from a single skeletal site. Using an osteological collection of 115 young adult male and female African-Americans, we assessed the functional relationship between bone robustness (i.e. total area/length), cortical tissue mineral density (Ct.TMD), and cortical area (Ct.Ar) for the upper and lower limbs. All long bones demonstrated significant trait covariance (p < 0.005) independent of body size, with slender bones having 25-50% less Ct.Ar and 5-8% higher Ct.TMD compared to robust bones. Robustness statistically explained 10.2-28% of Ct.TMD and 26.6-64.6% of Ct.Ar within male and female skeletal elements. This covariance is systemic throughout the skeleton, with either the slender or robust phenotype consistently represented within all long bones for each individual. These findings suggest that each person attains a unique trait set by adulthood that is both predictable by robustness and partially independent of environmental influences. The variation in these functionally integrated traits allows for the maximization of tissue stiffness and minimization of mass so that regardless of which phenotype is present, a given bone is reasonably stiff and strong, and sufficiently adapted to perform routine, habitual loading activities. Covariation intrinsic to functional adaptation suggests that whole bone stiffness depends upon particular sets of traits acquired during growth, presumably through differing levels of cellular activity, resulting in differing tissue morphology and composition. The outcomes of this intraskeletal examination of robustness and its correlates may have significant value in our progression towards improved clinical assessments of bone strength and fragility.

The genetic profile of bone repair

Bone repair following a fracture is a complex, well orchestrated, physiological process in response to injury. Even though the exact number of the genes and expressed proteins involved in fracture healing remains unknown, the molecular complexity of the repair process has been demonstrated, and it involves numerous genes and molecules, such as extracellular matrix genes, growth and differentiation factors, matrix metalloproteinases, angiogenic factors and others. Discrepancies in fracture healing responses and final outcome seen in the clinical practice may be attributed among other factors to biological variations between patients and different genetic "profiles", resulting in "altered" signalling pathways that regulate the bone repair process. Preliminary human studies support a "genetic" component in the pathophysiology of impaired bone repair seen in atrophic non-unions by correlating genetic variations of specific molecules regulating fracture healing with non-union. However, the role of the genetic "profile" of each individual in fracture healing and final outcome, and its possible interaction with other exogenous factors remains a topic of extensive research.

Systematic evaluation of skeletal mechanical function

Many genetic and environmental perturbations lead to measurable changes in bone morphology, matrix composition, and matrix organization. The skeletal system is highly adaptive, such that genetic and environmental perturbations affecting one trait are often accompanied by compensatory changes in other traits. This adaptation should be considered when attempting to draw meaningful conclusions about the role of a gene, for example. The challenge is to figure out which alterations can be attributed to the perturbation and which result from adaptive changes associated with establishing mechanical function. The goal of this chapter is to describe straight-forward biomechanical methods that can be used to determine whether a genetic or environmental perturbation affected bone strength. We introduce a systematic method for evaluating how bone strength was altered in the context of morphology and tissue-level mechanical properties, which are determined in large part from matrix composition, matrix organization, and porosity. We present this work as a first step toward screening mice for a phenotypic effect and for establishing the associated biomechanical mechanism by which function was altered. The outcome of these analyses generally provides insight into the next set of experiments required to further connect the cellular perturbation with the functional changes. The protocols were written to enable researchers without a background in engineering to conduct the assays or to enable researchers to better understand the outcomes of similar assays conducted by colleagues knowledgeable in engineering.

Genetic predisposition to non-union: evidence today

Atrophic non-union represents a complex clinical condition and research is ongoing in an effort to elucidate its pathophysiology and to offer new and more efficient treatment modalities. Differences seen in fracture healing responses and final outcome may be attributed among other factors to biological variations between patients resulting in a "disturbed" signalling pathway and an "inert or deficient local biology with reduced potentials for bone regeneration". The genetic contribution with or without the interaction of other exogenous factors in cases of impaired fracture healing, is yet to be elucidated. However, preliminary animal and human studies demonstrate the molecular basis of fracture non-unions and correlate genetic variants of the molecules regulating fracture healing and their expression patterns with impaired bone healing and fracture non-union. Further research is needed to clarify the genetic component and its role and interaction with other risk factors that may result in increased susceptibility of a patient to develop this complication.

The influence of age at menarche on cross-sectional geometry of bone in young adulthood

Elucidating the somatic and maturational influences on the biomechanical properties of bone in children is crucial for a proper understanding of bone strength and quality in childhood and later life, and has significant potential for predicting adult fracture and osteoporosis risks. The ability of a long bone to resist bending and torsion is primarily a function of its cross-sectional geometric properties, and is negatively impacted by smaller external bone diameter. In pubescent girls, elevated levels of estrogen impede subperiosteal bone growth and increase endosteal bone deposition, resulting in bones averaging a smaller external and internal diameter relative to boys. In addition, given a well-documented secular trend for an earlier menarche, the age at which the rate of subperiosteal bone deposition decreases may also be younger in more recent cohorts of girls. In this study we examined the relationship between pubertal timing and subsequent bone strength in girls. Specifically, we investigated the effects of age at menarche on bone strength indicators (polar moment of inertia and section modulus) determined from cross-sectional geometry of the second metacarpal (MC2) using data derived from serial hand-wrist radiographs of female participants (N=223) in the Fels Longitudinal Study, with repeated measures of MC2 between the ages of 7 and 35 years. Using multivariate regression models, we evaluated the effects of age at menarche on associations between measures of bone strength in early adulthood and the same measures at a prepubertal age. Results indicate that later age at menarche is associated with stronger adult bone (in torsion and bending) when controlling for prepubertal bone strength (R(2) ranged between 0.54 and 0.70, p<0.001). Since cross-sectional properties of bone in childhood may have long lasting implications, they should be considered along with pubertal timing in assessing risk for future fracture and in clinical recommendations.

Dissection of complex adult traits in a mouse synthetic population

Finding the causative genetic variations that underlie complex adult traits is a significant experimental challenge. The unbiased search strategy of genome-wide association (GWAS) has been used extensively in recent human population studies. These efforts, however, typically find only a minor fraction of the genetic loci that are predicted to affect variation. As an experimental model for the analysis of adult polygenic traits, we measured a mouse population for multiple phenotypes and conducted a genome-wide search for effector loci. Complex adult phenotypes, related to body size and bone structure, were measured as component phenotypes, and each subphenotype was associated with a genomic spectrum of candidate effector loci. The strategy successfully detected several loci for the phenotypes, at genome-wide significance, using a single, modest-sized population (N = 505). The effector loci each explain 2%–10% of the measured trait variation and, taken together, the loci can account for over 25% of a trait's total population variation. A replicate population (N = 378) was used to confirm initially observed loci for one trait (femur length), and, when the two groups were merged, the combined population demonstrated increased power to detect loci. In contrast to human population studies, our mouse genome-wide searches find loci that individually explain a larger fraction of the observed variation. Also, the additive effects of our detected mouse loci more closely match the predicted genetic component of variation. The genetic loci discovered are logical candidates for components of the genetic networks having evolutionary conservation with human biology.

Biological constraints that limit compensation of a common skeletal trait variant lead to inequivalence of tibial function among healthy young adults

Having a better understanding of how complex systems like bone compensate for the natural variation in bone width to establish mechanical function will benefit efforts to identify traits contributing to fracture risk. Using a collection of pQCT images of the tibial diaphysis from 696 young adult women and men, we tested the hypothesis that bone cells cannot surmount the nonlinear relationship between bone width and whole bone stiffness to establish functional equivalence across a healthy population. Intrinsic cellular constraints limited the degree of compensation, leading to functional inequivalence relative to robustness, with slender tibias being as much as two to three times less stiff relative to body size compared with robust tibias. Using Path Analysis, we identified a network of compensatory trait interactions that explained 79% of the variation in whole-bone bending stiffness. Although slender tibias had significantly less cortical area relative to body size compared with robust tibias, it was the limited range in tissue modulus that was largely responsible for the functional inequivalence. Bone cells coordinately modulated mineralization as well as the cortical porosity associated with internal bone multicellular units (BMU)-based remodeling to adjust tissue modulus to compensate for robustness. Although anecdotal evidence suggests that functional inequivalence is tolerated under normal loading conditions, our concern is that the functional deficit of slender tibias may contribute to fracture susceptibility under extreme loading conditions, such as intense exercise during military training or falls in the elderly. Thus, we show the natural variation in bone robustness was associated with predictable functional deficits that were attributable to cellular constraints limiting the amount of compensation permissible in human long bone. Whether these cellular constraints can be circumvented prophylactically to better equilibrate function among individuals remains to be determined.

Variation in childhood skeletal robustness is an important determinant of cortical area in young adults

A better understanding of bone growth will benefit efforts to reduce fracture incidence, because variation in elderly bone traits is determined primarily by adulthood. The natural variation in robustness was used as a model to understand how variable growth patterns define adult bone morphology. Longitudinally acquired hand radiographs of 29 boys and 30 girls were obtained from the Bolton-Brush study for 6 time points spanning 8 to 18 years of age. Segregating individuals into tertiles based on robustness revealed that the biological activity underlying bone growth varied significantly with the natural variation in robustness. For boys, slender metacarpals used an osteoblast-dependent growth pattern to establish function, whereas robust metacarpals used an osteoclast-dependent growth pattern. In contrast, differences in biological activity between girls with slender and robust metacarpals were largely based on the age at which the marrow surface changed from expansion to infilling. Importantly, cortical area for slender metacarpals was as much as 19.7% and 32.2% lower than robust metacarpals for boys and girls, respectively, indicating that robustness was a major determinant of adult cortical area. Finally, after accounting for robustness and body weight effects, we found that the inter-individual variation in cortical area was established as early as 8 years of age. While variation in the amount of bone acquired during growth has primarily been attributed to factors like nutrition, exercise, and genetic background, we showed that the natural variation in robustness was also a major determinant of cortical area, which is an important determinant of bone mass. This predictable relationship between robustness and cortical area should be incorporated into clinical diagnostic measures and experimental studies.

Functional interactions among morphologic and tissue quality traits define bone quality

BACKGROUND

Advances in diagnostic and treatment regimens that aim to reduce fracture incidence will benefit from a better understanding of how bone morphology and tissue quality define whole-bone mechanical properties.

QUESTIONS/PURPOSES

The goal of this article was to review what is known about the interactions among morphologic and tissue quality traits and how these interactions contribute to bone quality (ie, whole-bone mechanical function). Several questions were addressed. First, how do interactions among morphology and tissue quality traits relate to functional adaptation? Second, what are the emergent patterns of functionally adapted trait sets in long bones? Third, how effective is phenotypic integration at establishing function across a population? Fourth, what are the emergent patterns of functionally adapted trait sets in corticocancellous structures? Fifth, how do functional interactions change with aging?

METHODS

A literature review was conducted with papers identified primarily through citations listed in reference sections as well as general searches using Google Scholar and PubMed.

RESULTS

The interactions among adult traits or phenotypic integration are an emergent property of the compensatory mechanisms complex systems used to establish function or homeostasis. Traits are not regulated independently but vary simultaneously (ie, covary) in specific ways to establish function. This covariation results in individuals acquiring unique sets of traits to establish bone quality.

CONCLUSIONS AND CLINICAL RELEVANCE

Biologic constraints imposed on the skeletal system result in a population showing a pattern of trait sets that is predictable based on external bone size and that can be used to identify individuals with reduced bone quality relative to their bone size and body size.

Comprehensive skeletal phenotyping and linkage mapping in an intercross of recombinant congenic mouse strains HcB-8 and HcB-23

Bone biomechanical performance is a complex trait or, more properly, an ensemble of complex traits. Biomechanical performance incorporates flexibility under loading, yield and failure load, and energy to failure; all are important measures of bone function. To date, the vast majority of work has focused on yield and failure load and its surrogate, bone mineral density. We performed a reciprocal intercross of the mouse strains HcB-8 and HcB-23 to map and ultimately identify genes that contribute to differences in biomechanical performance. Mechanical testing was performed by 3-point bending of the femora. We measured femoral diaphysis cross-sectional anatomy from photographs of the fracture surfaces. We used beam equations to calculate material level mechanical properties. We performed a principal component (PC) analysis of normalized whole bone phenotypes (17 input traits). We measured distances separating mandibular landmarks from calibrated digital photographs and performed linkage analysis. Experiment-wide α = 0.05 significance thresholds were established by permutation testing. Three quantitative trait loci (QTLs) identified in these studies illustrate the advantages of the comprehensive phenotyping approach. A pleiotropic QTL on chromosome 4 affected multiple whole bone phenotypes with LOD scores as large as 17.5, encompassing size, cross-sectional ellipticity, stiffness, yield and failure load, and bone mineral density. This locus was linked to 3 of the PCs but unlinked to any of the tissue level phenotypes. From this pattern, we infer that the QTL operates by modulating the proliferative response to mechanical loading. On this basis, we successfully predicted that this locus also affects the length of a specific region of the mandible. A pleiotropic locus on chromosome 10 with LOD scores displays opposite effects on failure load and toughness with LOD scores of 4.5 and 5.5, respectively, so that the allele that increases failure load decreases toughness. A chromosome 19 QTL for PC2 with an LOD score of 4.8 was not detected with either the whole bone or tissue level phenotypes. We conclude that first, comprehensive, system-oriented phenotyping provides much information that could not be obtained by focusing on bone mineral density alone. Second, mechanical performance includes inherent trade-offs between strength and brittleness. Third, considering the aggregate phenotypic data allows prediction of novel QTLs.

Material and mechanical properties of bones deficient for fibrillin-1 or fibrillin-2 microfibrils

The contribution of non-collagenous components of the extracellular matrix to bone strength is largely undefined. Here we report that deficiency of fibrillin-1 or fibrillin-2 microfibrils causes distinct changes in bone material and mechanical properties. Morphometric examination of mice with hypomorphic or null mutations in fibrillin-1 or fibrillin-2, respectively, revealed appreciable differences in the postnatal shaping and growth of long bones. Fourier transform infrared imaging spectroscopy indicated that fibrillin-1 plays a predominantly greater role than fibrillin-2 in determining the material properties of bones. Biomechanical tests demonstrated that fibrillin-2 exerts a greater positive influence on the mechanical properties of bone than fibrillin-1 assemblies. Published evidence indirectly supports the notion that the above findings are mostly, if not exclusively, related to the differential control of TGFβ family signaling by fibrillin proteins. Our study therefore advances our understanding of the role that extracellular microfibrils play in bone physiology and implicitly, in the pathogenesis of bone loss in human diseases caused by mutations in fibrillin-1 or -2.

Serum IGF-1 affects skeletal acquisition in a temporal and compartment-specific manner

Insulin-like growth factor-1 (IGF-1) plays a critical role in the development of the growing skeleton by establishing both longitudinal and transverse bone accrual. IGF-1 has also been implicated in the maintenance of bone mass during late adulthood and aging, as decreases in serum IGF-1 levels appear to correlate with decreases in bone mineral density (BMD). Although informative, mouse models to date have been unable to separate the temporal effects of IGF-1 depletion on skeletal development. To address this problem, we performed a skeletal characterization of the inducible LID mouse (iLID), in which serum IGF-1 levels are depleted at selected ages. We found that depletion of serum IGF-1 in male iLID mice prior to adulthood (4 weeks) decreased trabecular bone architecture and significantly reduced transverse cortical bone properties (Ct.Ar, Ct.Th) by 16 weeks (adulthood). Likewise, depletion of serum IGF-1 in iLID males at 8 weeks of age, resulted in significantly reduced transverse cortical bone properties (Ct.Ar, Ct.Th) by 32 weeks (late adulthood), but had no effect on trabecular bone architecture. In contrast, depletion of serum IGF-1 after peak bone acquisition (at 16 weeks) resulted in enhancement of trabecular bone architecture, but no significant changes in cortical bone properties by 32 weeks as compared to controls. These results indicate that while serum IGF-1 is essential for bone accrual during the postnatal growth phase, depletion of IGF-1 after peak bone acquisition (16 weeks) is compartment-specific and does not have a detrimental effect on cortical bone mass in the older adult mouse.

Linkage mapping of principal components for femoral biomechanical performance in a reciprocal HCB-8 × HCB-23 intercross

Studies of bone genetics have addressed an array of related phenotypes, including various measures of biomechanical performance, bone size, bone, shape, and bone mineral density. These phenotypes are not independent, resulting in redundancy of the information they provide. Principal component (PC) analysis transforms multiple phenotype data to a new set of orthogonal "synthetic" phenotypes. We performed PC analysis on 17 femoral biomechanical, anatomic, and body size phenotypes in a reciprocal intercross of HcB-8 and HcB-23, accounting for 80% of the variance in 4 PCs. Three of the 4 PCs were mapped in the cross. The linkage analysis revealed a quantitative trait locus (QTL) with LOD = 4.7 for PC2 at 16 cM on chromosome 19 that was not detected using the directly measured phenotypes. The chromosome 19 QTL falls within a ~10 megabase interval, with Osf1 as a positional candidate gene. PC QTLs were also found on chromosomes 1, 2, 4, 6, and 10 that coincided with those identified for directly measured or calculated material property phenotypes. The novel chromosome 19 QTL illustrates the power advantage that attends use of PC phenotypes for linkage mapping. Constraint of the chromosome 19 candidate interval illustrates an important advantage of experimental crosses between recombinant congenic mouse strains.

Functional significance of genetic variation underlying limb bone diaphyseal structure

Limb bone diaphyseal structure is frequently used to infer hominin activity levels from skeletal remains, an approach based on the well-documented ability of bone to adjust to its loading environment during life. However, diaphyseal structure is also determined in part by genetic factors. This study investigates the possibility that genetic variation underlying diaphyseal structure is influenced by the activity levels of ancestral populations and might also have functional significance in an evolutionary context. We adopted an experimental evolution approach and tested for differences in femoral diaphyseal structure in 1-week-old mice from a line that had been artificially selected (45 generations) for high voluntary wheel running and non-selected controls. As adults, selected mice are significantly more active on wheels and in home cages, and have thicker diaphyses. Structural differences at 1 week can be assumed to primarily reflect the effects of selective breeding rather than direct mechanical stimuli, given that the onset of locomotion in mice is shortly after Day 7. We hypothesized that if genetically determined diaphyseal structure reflects the activity patterns of members of a lineage, then selected animals will have relatively larger diaphyseal dimensions at 1 week compared to controls. The results provide strong support for this hypothesis and suggest that limb bone cross sections may not always only reflect the activity levels of particular fossil individuals, but also convey an evolutionary signal providing information about hominin activity in the past.

Genetically determined phenotype covariation networks control bone strength

To identify genes affecting bone strength, we studied how genetic variants regulate components of a phenotypic covariation network that was previously shown to accurately characterize the compensatory trait interactions involved in functional adaptation during growth. Quantitative trait loci (QTLs) regulating femoral robustness, morphologic compensation, and mineralization (tissue quality) were mapped at three ages during growth using AXB/BXA Recombinant Inbred (RI) mouse strains and adult B6-i(A) Chromosome Substitution Strains (CSS). QTLs for robustness were identified on chromosomes 8, 12, 18, and 19 and confirmed at all three ages, indicating that genetic variants established robustness postnatally without further modification. A QTL for morphologic compensation, which was measured as the relationship between cortical area and body weight, was identified on chromosome 8. This QTL limited the amount of bone formed during growth and thus acted as a setpoint for diaphyseal bone mass. Additional QTLs were identified from the CSS analysis. QTLs for robustness and morphologic compensation regulated bone structure independently (ie, in a nonpleiotropic manner), indicating that each trait may be targeted separately to individualize treatments aiming to improve strength. Multiple regression analyses showed that variation in morphologic compensation and tissue quality, not bone size, determined femoral strength relative to body weight. Thus an individual inheriting slender bones will not necessarily inherit weak bones unless the individual also inherits a gene that impairs compensation. This systems genetic analysis showed that genetically determined phenotype covariation networks control bone strength, suggesting that incorporating functional adaptation into genetic analyses will advance our understanding of the genetic basis of bone strength.

Genetic variation in the structural pattern of osteoclast activity during post-natal growth of mouse femora

While the spatial activity of osteoblasts has been associated with modeling of bones during development, few studies have examined if variation in the spatial activity of osteoclasts also contributes to the morphogenesis of skeletal tissues. We examined this question by histomorphometric analysis and reconstructing the three-dimensional spatial distribution of osteoclasts in the femora of three inbred strains of male mice (A/J, C57BL/6J [B6], and C3H/HeJ [C3H]) that have differing skeletal, structural, and material properties. Our data show that total osteoclast surface area and osteoclast numbers are related to the overall bone density, but not related to the development of bone diameter or overall cortical area. The analysis of the spatial distribution of the osteoclasts showed that the asymmetrical mid-diaphyseal distribution of osteoclasts in A/J and B6 compared to the more uniform distribution of these cells around the circumference in the C3H mice was consistent with the more ellipsoid shape of A/J and B6 femora compared to the more circular mid-diaphyseal shape of the femora in the C3H mice. The statistically 2- to 3-fold fewer cells on the periosteal surface in the C3H compared to either the B6 or A/J mice is also consistent with the greater cortical thickness that is seen for the C3H mice compared to either B6 or A/J strains. In vitro studies of osteoclastogenesis and the expression of numerous phenotypic properties of osteoclasts prepared from the three strains of mice showed that A/J and B6 mice developed statistically greater numbers of tartrate resistant acid phosphatase (TRAP) positive cells and expressed statistically higher levels of multiple mRNAs that are unique to differentiated osteoclasts than those isolated from the C3H strain. In summary, the 3D reconstructions and histomorphometric analysis suggest that genetic differences lead to spatial variation in the distribution of osteoclasts. These variations in spatial distribution of osteoclasts in turn contribute in part to the development of the structural variations of the femora that are seen in the three strains of mice. In vitro studies suggest that intrinsic genetic variation in osteoclastogenesis and their phenotypic expression may contribute to the differences in their functional activities that give rise to the unique spatial distributions of these cells in bones.

Linkage mapping of femoral material properties in a reciprocal intercross of HcB-8 and HcB-23 recombinant mouse strains

Skeletal fragility is an important health problem with a large genetic component. We performed a 603 animal F2 reciprocal intercross of the recombinant congenic strains HcB-8 and HcB-23 to genetically map quantitative trait loci (QTLs) for tissue-level femoral biomechanical performance. These included elastic and post-yield strain, Young's modulus, elastic and maximum stress, and toughness and were calculated from 3-point bend testing of femora by the application of standard beam equations. We mapped these with R/qtl and QTL Cartographer and established significance levels empirically by permutation testing. Significant QTLs for at least one trait are present on chromosomes 1, 6, and 10 in the full F2 population, with additional QTLs evident in subpopulations defined by sex and cross direction. On chromosome 10, we find a QTL for post-yield strain and toughness, phenotypes that have not been mapped previously. Notably, the HcB-8 allele at this QTL increases post-yield strain and toughness, but decreases bone mineral density (BMD), while the material property QTLs on chromosomes 1, 6, and at a second chromosome 10 QTL are independent of BMD. We find significant sex x QTL and cross direction x QTL interactions. A robust, pleiotropic chromosome 4 QTL that we previously reported at the whole-bone level showed no evidence of linkage at the tissue-level, supporting our interpretation that modeling capacity is its primary phenotype. Our data demonstrate an inverse relationship between femoral perimeter and Young's modulus, with R(2)=0.27, supporting the view that geometric and material bone properties are subject to an integrated set of regulatory mechanisms. Mapping QTLs for tissue-level biomechanical performance advances understanding of the genetic basis of bone quality.

Interindividual variation in functionally adapted trait sets is established during postnatal growth and predictable based on bone robustness

Adults acquire unique sets of morphological and tissue-quality bone traits that are predictable based on robustness and deterministic of strength and fragility. How and when individual trait sets arise during growth has not been established. Longitudinal structural changes of the metacarpal diaphysis were measured for boys and girls from 3 mo to 8 yr of age using hand radiographs obtained from the Bolton-Brush collection. Robustness varied approximately 2-fold among boys and girls, and individual values were established by 2 yr of age, indicating that genetic and environmental factors controlling the relationship between growth in width and growth in length were established early during postnatal growth. Significant negative correlations between robustness and relative cortical area and a significant positive correlation between robustness and a novel measure capturing the efficiency of growth indicated that coordination of the subperiosteal and endocortical surfaces was responsible for this population acquiring a narrow range of trait sets that was predictable based on robustness. Boys and girls with robust diaphyses had proportionally thinner cortices to minimize mass, whereas children with slender diaphyses had proportionally thicker cortices to maximize stiffness. Girls had more slender metacarpals with proportionally thicker cortices compared with boys at all prepubertal ages. Although postnatal growth patterns varied in fundamentally different ways with sex and robustness, the dependence of trait sets on robustness indicated that children sustained variants affecting subperiosteal growth because they shared a common biological factor regulating functional adaptation. Considering the natural variation in acquired trait sets may help identify determinants of fracture risk, because age-related bone loss and gain will affect slender and robust structures differently.

Serum IGF-1 determines skeletal strength by regulating subperiosteal expansion and trait interactions

Strong correlations between serum IGF-1 levels and fracture risk indicate that IGF-1 plays a critical role in regulating bone strength. However, the mechanism by which serum IGF-1 regulates bone structure and fracture resistance remains obscure and cannot be determined using conventional approaches. Previous analysis of adult liver-specific IGF-1-deficient (LID) mice, which exhibit 75% reductions in serum IGF-1 levels, showed reductions in periosteal circumference, femoral cross-sectional area, cortical thickness, and total volumetric BMD. Understanding the developmental sequences and the resultant anatomical changes that led to this adult phenotype is the key for understanding the complex relationship between serum IGF-1 levels and fracture risk. Here, we identified a unique developmental pattern of morphological and compositional traits that contribute to bone strength. We show that reduced bone strength associated with low levels of IGF-1 in serum (LID mice) result in impaired subperiosteal expansion combined with impaired endosteal apposition and lack of compensatory changes in mineralization throughout growth and aging. We show that serum IGF-1 affects cellular activity differently depending on the cortical surface. Last, we show that chronic reductions in serum IGF-1 indirectly affect bone strength through its effect on the marrow myeloid progenitor cell population. We conclude that serum IGF-1 not only regulates bone size, shape, and composition during ontogeny, but it plays a more fundamental role-that of regulating an individual's ability to adapt its bone structure to mechanical loads during growth and development.

A delay in pubertal onset affects the covariation of body weight, estradiol, and bone size

The skeletal system functions as a locomotive organ and a mineral reservoir and combinations of genetic and environmental factors affect the skeletal system. Although delayed puberty is associated with compromised bone mass, suppression of estrogen should be beneficial to cortical strength. The purpose was to employ path analysis to study bone strength and delayed puberty. Forty-five female rats were randomly assigned to a control group (n = 15) and an experimental group (n = 30) that received injections of gonadotropin releasing hormone antagonist (GnRH-a). Causal models were constructed by specifying directed paths between bone traits. The first model tested the hypothesis that the functional relationships between bone traits and body weight were altered by a delay in pubertal onset. GnRH-a injections during puberty altered the covariation between body weight and bone size. The second model was constructed to test the hypothesis that variability in stiffness was causally related to variability in body weight. The model also tested the relationship between the periosteal and endocortical surfaces and their relationship to stiffness. There was no change in the relationship between the surfaces in the GnRH-a group. The third model determined the effect of estradiol on both total area and relative cortical area in both groups. The relationship between periosteal surface and serum estradiol levels was only significant during estrogen suppression. These data suggest that increases in body weight during or prior to puberty may not be protective of bone strength.