Bone quality and bone strength in BXH recombinant inbred mice

Systems genetics in diversity outbred mice inform BMD GWAS and identify determinants of bone strength

Genome-wide association studies (GWASs) for osteoporotic traits have identified over 1000 associations; however, their impact has been limited by the difficulties of causal gene identification and a strict focus on bone mineral density (BMD). Here, we use Diversity Outbred (DO) mice to directly address these limitations by performing a systems genetics analysis of 55 complex skeletal phenotypes. We apply a network approach to cortical bone RNA-seq data to discover 66 genes likely to be causal for human BMD GWAS associations, including the genes SERTAD4 and GLT8D2. We also perform GWAS in the DO for a wide-range of bone traits and identify Qsox1 as a gene influencing cortical bone accrual and bone strength. In this work, we advance our understanding of the genetics of osteoporosis and highlight the ability of the mouse to inform human genetics.

Osteoporosis GWAS faces two challenges, causal gene discovery and a lack of phenotypic diversity. Here, the authors use the Diversity Outbred mouse population to inform human GWAS using networks and map genetic loci for 55 bone traits, identifying new potential bone strength genes.

Overexpression of constitutively active MAP3K7 in ameloblasts causes enamel defects of mouse teeth

OBJECTIVE

Compelling evidence suggests that mitogen-activated protein kinases (Mapks) play an important role in amelogenesis. However, the role of transforming growth factor (TGF)-β-activating kinase 1 (Tak1, Map3k7), which is a known upstream kinase of Mapks, during amelogenesis remains to be determined. The aim of this study was to investigate the possible involvement of Map3k7 in amelogenesis.

DESIGN

We generated transgenic mice that produced constitutively active human MAP3K7 (caMAP3K7) under the control of amelogenin (Amelx) gene promoter. Radiography and micro-computed tomography (μCT) analysis was used to detect the radio-opacity and density of the teeth. The enamel microstructure was observed with a scanning electron microscope. Histological analysis was used to observe the adhesion between ameloblasts and residual organic matrix of the enamel. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to analyze the expression of enamel matrix protein.

RESULTS

The enamel of mandibular molars in caMAP3K7-overexpressing mice displayed pigmentation and a highly irregular structure compared with the wild type littermates. Teeth of transgenic animals underwent rapid attrition due to the brittleness of the enamel layer. The microstructure of enamel, normally a highly ordered arrangement of hydroxyapatite crystals, was completely disorganized. The gross histological appearances of ameloblasts and supporting cellular structures, as well as the expression of the enamel protein amelotin (Amtn) were altered by the overexpression of caMAP3K7.

CONCLUSIONS

Our data demonstrated that protein expression, processing and secretion occurred abnormally in transgenic mice overexpressing caMAP3K7. The overexpression of caMAP3K7 had a profound effect on enamel structure by disrupting the orderly growth of enamel prisms.

Genetic Dissection of a QTL Affecting Bone Geometry

Parameters of bone geometry such as width, length, and cross-sectional area are major determinants of bone strength. Although these traits are highly heritable, few genes influencing bone geometry have been identified. Here, we dissect a major quantitative trait locus (QTL) influencing femur size. This QTL was originally identified in an F2 cross between the C57BL/6J-hg/hg (HG) and CAST/EiJ strains and was referred to as femur length in high growth mice 2 (Feml2). Feml2 was located on chromosome (Chr.) 9 at ∼20 cM. Here, we show that the HG.CAST-(D9Mit249-D9Mit133)/Ucd congenic strain captures Feml2. In an F2 congenic cross, we fine-mapped the location of Feml2 to an ∼6 Mbp region extending from 57.3 to 63.3 Mbp on Chr. 9. We have identified candidates by mining the complete genome sequence of CAST/EiJ and through allele-specific expression (ASE) analysis of growth plates in C57BL/6J × CAST/EiJ F1 hybrids. Interestingly, we also find that the refined location of Feml2 overlaps a cluster of six independent genome-wide associations for human height. This work provides the foundation for the identification of novel genes affecting bone geometry.

Enamel Hypomineralization and Structural Defects in Amelotin-deficient Mice

Amelotin (AMTN) is a relatively recently discovered enamel protein that is predominantly expressed by ameloblasts during the maturation stage of amelogenesis and is present at lower levels in the junctional epithelium of erupted teeth. Previous studies have suggested a function of this protein in enamel mineralization and cell attachment. Genetic mouse models have been instrumental in defining the role of many enamel-related proteins, but a genetic mouse model lacking the Amtn gene has not been reported. Here, we describe the generation of amelotin-deficient mice and the analysis of their enamel phenotype in comparison with that of wild-type animals. Ablation of AMTN expression resulted in mechanically inferior enamel of mandibular incisors that showed chipping and fractures at the incisal edge. Enamel mineralization was delayed, resulting in hypomineralized inner enamel and structural defects in the outer enamel. Erupted enamel close to the gingival margin showed increased surface roughness. The expression levels of the enamel matrix proteins AMEL, AMBN, ENAM, and ODAM and the enamel proteases MMP-20 and KLK-4 were not significantly altered, although the expression of KLK-4 was delayed. The morphology of ameloblasts showing prominent Tomes' processes during the secretory stage was not altered, and there was no indication of disruption of cell structures or activities, but a residual layer, presumably consisting of organic material, remained at the enamel surface close to the gingival margin. The integrity of the dentogingival attachment at the junctional epithelium appeared unaffected by AMTN deficiency. These observations indicate that AMTN plays a subtle yet critical role in enamel biomineralization, particularly during the establishment of the outer and surface enamel layers. This role appears to be largely independent of other enamel proteins.

High-throughput quantification of the mechanical competence of murine femora--a highly automated approach for large-scale genetic studies

Animal models are widely used to gain insight into the role of genetics on bone structure and function. One of the main strategies to map the genes regulating specific traits is called quantitative trait loci (QTL) analysis, which generally requires a very large number of animals (often more than 1000) to reach statistical significance. QTL analysis for mechanical traits has been mainly based on experimental mechanical testing, which, in view of the large number of animals, is time consuming. Hence, the goal of the present work was to introduce an automated method for large-scale high-throughput quantification of the mechanical properties of murine femora. Specifically, our aims were, first, to develop and validate an automated method to quantify murine femoral bone stiffness. Second, to test its high-throughput capabilities on murine femora from a large genetic study, more specifically, femora from two growth hormone (GH) deficient inbred strains of mice (B6-lit/lit and C3.B6-lit/lit) and their first (F1) and second (F2) filial offsprings. Automated routines were developed to convert micro-computed tomography (micro-CT) images of femora into micro-finite element (micro-FE) models. The method was experimentally validated on femora from C57BL/6J and C3H/HeJ mice: for both inbred strains the micro-FE models closely matched the experimentally measured bone stiffness when using a single tissue modulus of 13.06 GPa. The mechanical analysis of the entire dataset (n=1990) took approximately 44 CPU hours on a supercomputer. In conclusion, our approach, in combination with QTL analysis could help to locate genes directly involved in controlling bone mechanical competence.

The importance of the intracortical canal network for murine bone mechanics

As shown by recent data bone strength estimation can greatly be improved by including microarchitectural parameters in the analysis. Our previous results showed that intracortical canals (the living space of the vasculature and/or remodeling units) are a major contributor to cortical tissue porosity, and therefore, can be linked to mechanical bone properties. Consequently, the goal of this study was to investigate the importance of the intracortical canal network for murine bone mechanics. To study intracortical canals within murine femoral bone, we used a mouse model, including two mouse strains, C57BL/6J-Ghrhr(lit)/J (B6-lit/+) and C3.B6-Ghrhr(lit)/J (C3.B6-lit/+) representing low and high bone mass, respectively. The intracortical canal network was assessed by synchrotron radiation-based micro-computed tomography and the mechanical bone properties were derived from three-point bending experiments. Multiple linear regression models were built to explain the variation in ultimate force, work to fracture, and stiffness in terms of the morphometric parameters. The power to explain the variation in bone mechanics was increased significantly for most mechanical measures when including morphometric parameters of intracortical canals in addition to macroscopic morphometric measures. Specifically, we could derive generalized (mouse strain-independent) models for ultimate force, where the incorporation of intracortical canals in addition to macroscopic bone measures improved the explained variation in ultimate force considerably, which was confirmed by an increase in adjusted R(2) of 73% and 8% for B6-lit/+ and C3.B6-lit/+, respectively. Further, we observed that the heterogeneity of the morphometric measures for the individual canal branches play an important role for explaining the variation in ultimate force. Finally, the current study provides strong evidence that work to fracture of murine bone, which is triggered critically by microcracks, is affected by intracortical canals. In summary, the study suggests that the intracortical canal network is important for bone mechanics.

Identification of quantitative trait loci influencing skeletal architecture in mice: emergence of Cdh11 as a primary candidate gene regulating femoral morphology

Bone strength is influenced by many properties intrinsic to bone, including its mass, geometry, and mineralization. To further advance our understanding of the genetic basis of bone-strength-related traits, we used a large (n = 815), moderately (G(4) ) advanced intercross line (AIL) of mice derived from a high-runner selection line (HR) and the C57BL/6J inbred strain. In total, 16 quantitative trait loci (QTLs) were identified that affected areal bone mineral density (aBMD) and femoral length and width. Four significant (p < .05) and one suggestive (p < .10) QTLs were identified for three aBMD measurements: total body, vertebral, and femoral. A QTL on chromosome (Chr.) 3 influenced all three aBMD measures, whereas the other four QTLs were unique to a single measure. A total of 10 significant and one suggestive QTLs were identified for femoral length (FL) and two measures of femoral width, anteroposterior (AP) and mediolateral (ML). FL QTLs were distinct from loci affecting AP and ML width, and of the 7 AP QTLs, only three affected ML. A QTL on Chr. 8 that explained 7.1% and 4.0% of the variance in AP and ML, respectively, was mapped to a 6-Mb region harboring 12 protein-coding genes. The pattern of haplotype diversity across the QTL region and expression profiles of QTL genes suggested that of the 12, cadherin 11 (Cdh11) was most likely the causal gene. These findings, when combined with existing data from gene knockouts, identify Cdh11 as a strong candidate gene within which genetic variation may affect bone morphology.

Phenotypic integration among trabecular and cortical bone traits establishes mechanical functionality of inbred mouse vertebrae

Conventional approaches to identifying quantitative trait loci (QTLs) regulating bone mass and fragility are limited because they examine cortical and trabecular traits independently. Prior work examining long bones from young adult mice and humans indicated that skeletal traits are functionally related and that compensatory interactions among morphological and compositional traits are critical for establishing mechanical function. However, it is not known whether trait covariation (i.e., phenotypic integration) also is important for establishing mechanical function in more complex, corticocancellous structures. Covariation among trabecular, cortical, and compositional bone traits was examined in the context of mechanical functionality for L(4) vertebral bodies across a panel of 16-wk-old female AXB/BXA recombinant inbred (RI) mouse strains. The unique pattern of randomization of the A/J and C57BL/6J (B6) genome among the RI panel provides a powerful tool that can be used to measure the tendency for different traits to covary and to study the biology of complex traits. We tested the hypothesis that genetic variants affecting vertebral size and mass are buffered by changes in the relative amounts of cortical and trabecular bone and overall mineralization. Despite inheriting random sets of A/J and B6 genomes, the RI strains inherited nonrandom sets of cortical and trabecular bone traits. Path analysis, which is a multivariate analysis that shows how multiple traits covary simultaneously when confounding variables like body size are taken into consideration, showed that RI strains that tended to have smaller vertebrae relative to body size achieved mechanical functionality by increasing mineralization and the relative amounts of cortical and trabecular bone. The interdependence among corticocancellous traits in the vertebral body indicated that variation in trabecular bone traits among inbred mouse strains, which is often thought to arise from genetic factors, is also determined in part by the adaptive response to variation in traits describing the cortical shell. The covariation among corticocancellous traits has important implications for genetic analyses and for interpreting the response of bone to genetic and environmental perturbations.

Strain-specific effects of rosiglitazone on bone mass, body composition, and serum insulin-like growth factor-I

Activation of peroxisome proliferator activated receptor-gamma (PPARG) is required for the differentiation of marrow mesenchymal stem cell into adipocytes and is associated with the development of age-related marrow adiposity in mice. Thiazolidinediones are agonists for PPARG and have a heterogeneous effect on bone mineral density (BMD). We postulated that genetic determinants influence the skeletal response to thiazolidinediones. We examined the effects of rosiglitazone (3 mg/kg . d for 8 wk) on BMD, body composition, and serum IGF-I in adult female mice from four inbred strains. C3H/HeJ mice showed the most significant response to treatment, exhibiting decreased femoral and vertebral BMD, reduced distal femoral bone volume fraction and a decrease in serum IGF-I. In DBA/2J, there were no changes in femoral BMD or bone volume fraction, but there was a decrease in vertebral BMD. C57BL/6J mice showed increases in marrow adiposity, without associated changes in trabecular bone volume; the skeletal effects from rosiglitazone in A/J mice were minimal. No association between trabecular bone volume and marrow adiposity was found. The effect of rosiglitazone on gene expression in the femur was then examined in the C3H/HeJ and C57BL/6J strains by microarray. Increased gene expression was observed in the PPARG signaling pathway and fatty acid metabolism in both C3H/HeJ and C57BL/6J, but a significant down-regulation of genes associated with cell cycle was noted only in the C3H/HeJ strain. The divergent skeletal responses to rosiglitazone in this study suggest the existence of a strong genetic background effect.

Effects of tissue preservation on murine bone mechanical properties

Murine bone specimens are used extensively in skeletal research to assess the effects of environmental, physiologic and pathologic factors on their mechanical properties. Given the destructive nature of mechanical testing, it is normally performed as a terminal procedure, where specimens must be preserved without affecting their mechanical properties. To this end, we aimed to study the effects of tissue preservation (freezing and formalin fixation) on the elastic and viscoelastic mechanical properties of murine femur and vertebrae. A total of 120 femurs and 180 vertebral bodies (L3-L5) underwent non-destructive cyclic loading to assess their viscoelastic properties followed by mono-cyclic loading to failure to assess their elastic properties. All specimens underwent re-hydration in 0.9% saline for 30min prior to mechanical testing. Analysis indicated that stiffness, modulus of elasticity, yield load, yield strength, ultimate load and ultimate strength of frozen and formalin-fixed femurs and vertebrae were not different from fresh specimens. Cyclic loading of both femurs and vertebrae indicated that loss, storage and dynamic moduli were not affected by freezing. However, formalin fixation altered their viscoelastic properties. Our findings suggest that freezing and formalin fixation over a 2-week period do not alter the elastic mechanical properties of murine femurs and vertebrae, provided that specimens are re-hydrated for at least half an hour prior to testing. However, formalin fixation weakened the viscoelastic properties of murine bone by reducing its ability to dissipate viscous energy. Future studies should address the long-term effects of both formalin fixation and freezing on the mechanical properties of murine bone.

Fourier transform infrared imaging microspectroscopy and tissue-level mechanical testing reveal intraspecies variation in mouse bone mineral and matrix composition

Fracture susceptibility is heritable and dependent upon bone morphology and quality. However, studies of bone quality are typically overshadowed by emphasis on bone geometry and bone mineral density. Given that differences in mineral and matrix composition exist in a variety of species, we hypothesized that genetic variation in bone quality and tissue-level mechanical properties would also exist within species. Sixteen-week-old female A/J, C57BL/6J (B6), and C3H/HeJ (C3H) inbred mouse femora were analyzed using Fourier transform infrared imaging and tissue-level mechanical testing for variation in mineral composition, mineral maturity, collagen cross-link ratio, and tissue-level mechanical properties. A/J femora had an increased mineral-to-matrix ratio compared to B6. The C3H mineral-to-matrix ratio was intermediate of A/J and B6. C3H femora had reduced acid phosphate and carbonate levels and an increased collagen cross-link ratio compared to A/J and B6. Modulus values paralleled mineral-to-matrix values, with A/J femora being the most stiff, B6 being the least stiff, and C3H having intermediate stiffness. In addition, work-to-failure varied among the strains, with the highly mineralized and brittle A/J femora performing the least amount of work-to-failure. Inbred mice are therefore able to differentially modulate the composition of their bone mineral and the maturity of their bone matrix in conjunction with tissue-level mechanical properties. These results suggest that specific combinations of bone quality and morphological traits are genetically regulated such that mechanically functional bones can be constructed in different ways.

Small animal bone biomechanics

Animal models, in particular mice, offer the possibility of naturally achieving or genetically engineering a skeletal phenotype associated with disease and conducting destructive fracture tests on bone to determine the resulting change in bone's mechanical properties. Several recent developments, including nano- and micro-indentation testing, microtensile and microcompressive testing, and bending tests on notched whole bone specimens, offer the possibility to mechanically probe small animal bone and investigate the effects of aging, therapeutic treatments, disease, and genetic variation. In contrast to traditional strength tests on small animal bones, fracture mechanics tests display smaller variation and therefore offer the possibility of reducing sample sizes. This article provides an analysis of what such tests measure and proposes methods to reduce errors associated with testing smaller than ideal specimens.

Differential effects of bone structural and material properties on bone competence in C57BL/6 and C3H/He inbred strains of mice

The femoral neck is a relevant and sensitive site for studying the degree of osteopenia. Engineering principles predict that bone structural parameters, like cross-sectional geometry, are important determinants of bone mechanical parameters. Mechanical parameters are also directly affected by the material properties of the bone tissue. However, the relative importance of structural and material properties is still unknown. The aim of this study was to compare bone competence and structural parameters between a murine strain showing a low bone mass phenotype, C57BL/6 (B6), and another one showing a high bone mass phenotype, C3H/He (C3H), in order to better determine the role of bone structure and geometry in bone failure behavior. Murine femora of 12- and 16-week-old B6 and 12- and 16-week-old C3H inbred strains were mechanically tested under axial loading of the femoral head. In order to assess the structural properties, we performed three-dimensional morphometric analyses in five different compartments of the mouse femur using micro-computed tomography. The mechanical tests revealed that B6 femora became stiffer, stronger, and tougher at 12-16 weeks, while bone brittleness stayed constant. C3H bone stiffness increased, but strength remained constant, work to failure decreased, and bone became more brittle. These age effects indicated that B6 did not reach peak bone properties at 16 weeks of age and C3H did reach maximal skeletal biomechanical properties before 16 weeks of age. Our investigations showed that 83% of the strength of the femoral neck in the B6 strain was explained by cortical thickness at this location; in contrast, in C3H none of the mechanical properties of the femoral neck was explained by bone structural parameters. The relative contributions of bone structural and material properties on bone strength are different in B6 and C3H. We hypothesize that these different contributions are related to differences at the ultrastructural level of bone that affect bone failure.