Bone morphometry strongly predicts cortical bone stiffness and strength, but not toughness, in inbred mouse models of high and low bone mass

A Meta-Analysis Study to Define Variations in Murine Long Bone Biomechanical Testing

A systematic literature search and meta-analysis were performed to evaluate the variability in biomechanical testing of murine long bones, specifically focused on point-bending tests of mice femora. Due to the lack of standardized protocols for these tests, the assessment quantifies the heterogeneity in reported mechanical properties across existing literature. This study followed preferred reporting items for systematic reviews and meta-analyses (PRISMA) and strengthening the reporting of observational studies in epidemiology (STROBE) guidelines to search publicly available databases for relevant studies. After title and abstract screening, full-text reviews identified 73 articles meeting the inclusion criteria. Data was extracted from these studies, including stiffness, maximum load, modulus, and ultimate stress values for both three-point and four-point bending tests. The data were analyzed through ANOVA and metaregression to assess variability caused by age, sex, and genetic strain. The reviewers also assessed the quality of the included studies. The meta-analysis revealed significant heterogeneity in reported mechanical properties, with I2 values ranging from 72% to 100% in the three point-bend tests of pooled genetic strains. This heterogeneity persisted even after accounting for age, sex, and genetic strain differences. The review concludes that nonstandardized testing setups are the likely major source of the observed variability in reported data more than the population characteristics of the mice, highlighting the need for more consistent testing methodologies in future studies.

Do Intramedullary Screws Provide Adequate Fixation for Humeral and Ulnar Components in Total Elbow Arthroplasty? A Cadaveric Analysis

PURPOSE

The objective of this study was to determine the structural properties of the cadaver bone-screw interface for cementless intramedullary screw fixation in the context of total elbow arthroplasty.

METHODS

The intramedullary canals of seven humerus and seven ulna specimens from fresh-frozen cadavers were drilled using custom drill bits until the inner cortex was reached and then hand tapped for the corresponding thread size. Titanium screws were advanced into the tapped holes until securely seated. The bones were potted and then mounted on a uniaxial material testing machine. A tensile load was applied, and end-of-test elongation, failure load, energy absorbed, and stiffness were determined. End-of-test load and elongation were defined as the elongation and load experienced by the structure at 3,000 N or failure. Each specimen was inspected for evidence of pullout, loosening, or visible fractures.

RESULTS

The end-of-test load and elongation for the humerus specimens were 2721 ± 738 N and 3.0 ± 0.9 mm, respectively. The ulna specimens reached 92% of the humerus specimens' end-of-test load at 2,514 ± 678 N and 120% of their end-of-test elongation (3.6 ± 0.6 mm). The stiffness of the humerus specimens was 1,077 ± 336 N/mm, which was 1.3 times greater than the stiffness of the ulna specimens (790 ± 211 N/mm). Lastly, the energy absorbed by the humerus samples was 3.6 ± 1.6 J, which was 92% of the energy absorbed by the ulna samples at 3.9 ± 1.1 J. One humerus and three ulnas failed before the end-of-test load of 3,000 N. Two failures were caused by screw pullout and two by bone fracture.

CONCLUSIONS

Our findings demonstrate that intramedullary screw fixation is successful in withstanding forces that are greater than required for osseointegration.

CLINICAL RELEVANCE

Uncemented fixation may be beneficial in elbow arthroplasty.

Bone marrow dosimetry for mice: exposure from bone-seeking 89,90Sr

Studies of radiobiological effects in murine rodents exposed to internal radiation in the wild or in laboratory experiments require dosimetric support. The main problem of bone marrow (BM) dosimetry for bone-seeking β-emitters is dosimetric modeling, because the bone is a heterogeneous structure with complex microarchitecture. To date, there are several approaches to calculating the absorbed dose in BM, which mostly use rough geometric approximations. Recently, in the framework of studies of people exposed to ⁹⁰Sr in the Urals, a new approach (SPSD) has been developed. The aim of the current study was to test for the first time the possibility of extension of the SPSD approach elaborated for humans to mice. For this, computational phantoms of femur bones of laboratory animals (C57BL/6, C57BL/6 J, BALB/c, BALB/cJ) aged 5-8 weeks (growing) and > 8 weeks (adults) were created. The dose factors DF_Sr-90(BM ← TBV + CBV) to convert the Sr isotope activity concentration in a bone tissue into units of dose rate absorbed in the bone marrow were 1.75 ± 0.42 and 2.57 ± 0.93 μGy day^-1 per Bq g^-1 for growing and adult animals, respectively, while corresponding values for DF_Sr-89(BM ← TBV + CBV) were 1.08 ± 0.27 and 1.66 ± 0.67 μGy day^-1 per Bq g^-1, respectively. These results are about 2.5 times lower than skeleton-average DFs calculated assuming homogenous bone, where source and target coincide. The results of the present study demonstrate the possibility of application of the SPSD approach elaborated for humans to non-human mammals. It is concluded that the study demonstrates the feasibility and appropriateness of application of the SPSD approach elaborated for humans to non-human mammals. This approach opens up new prospects for studying the radiobiological consequences of red bone marrow exposure for both laboratory and wildlife mammals.

Cortical bone relationships are maintained regardless of sex and diet in a large population of LGXSM advanced intercross mice

Introduction

Knowledge of bone structure-function relationships in mice has been based on relatively small sample sets that limit generalizability. We sought to investigate structure-function relationships of long bones from a large population of genetically diverse mice. Therefore, we analyzed previously published data from the femur and radius of male and female mice from the F34 generation of the Large-by-Small advanced intercross line (LGXSM AI), which have over a two-fold continuous spread of bone and body sizes (Silva et al. 2019 JBMR).

Methods

Morphological traits, mechanical properties, and estimated material properties were collected from the femur and radius from 1113 LGXSM AI adult mice (avg. age 25 wks). Males and females fed a low-fat or high-fat diet were evaluated to increase population variation. The data were analyzed using principal component analysis (PCA), Pearson's correlation, and multivariate linear regression.

Results

Using PCA groupings and hierarchical clustering, we identified a reduced set of traits that span the population variation and are relatively independent of each other. These include three morphometry parameters (cortical area, medullary area, and length), two mechanical properties (ultimate force and post-yield displacement), and one material property (ultimate stress). When comparing traits of the femur to the radius, morphological traits are moderately well correlated (r²: 0.18-0.44) and independent of sex and diet. However, mechanical and material properties are weakly correlated or uncorrelated between the long bones. Ultimate force can be predicted from morphology with moderate accuracy for both long bones independent of variations due to genetics, sex, or diet; however, predictions miss up to 50 % of the variation in the population. Estimated material properties in the femur are moderately to strongly correlated with bone size parameters, while these correlations are very weak in the radius.

Discussion

Our results indicate that variation in cortical bone phenotype in the F34 LGXSM AI mouse population can be adequately described by a reduced set of bone traits. These traits include cortical area, medullary area, bone length, ultimate force, post-yield displacement, and ultimate stress. The weak correlation of mechanical and material properties between the femur and radius indicates that the results from routine three-point bending tests of one long bone (e.g., femur) may not be generalizable to another long bone (e.g., radius). Additionally, these properties could not be fully predicted from bone morphology alone, confirming the importance of mechanical testing. Finally, material properties of the femur estimated based on beam theory equations showed a strong dependence on geometry that was not seen in the radius, suggesting that differences in femur size within a study may confound interpretation of estimated material properties.

Site-specific characteristics of bone and progenitor cells in control and ovariectomized rats

One-third of postmenopausal women experience at least one osteoporotic bone fracture in their lifetime that occurs spontaneously or from low-impact events. However, osteoporosis-associated jaw bone fractures are extremely rare. It was also observed that jaw bone marrow stem cells (BMSCs) have a higher capacity to form mineralized tissues than limb BMSCs. At present, the underlying causes and mechanisms of variations between jaw bone and limb bone during postmenopause are largely unknown. Thus, the objective of the current study was to examine the site-specific effects of estrogen deficiency using comprehensive analysis of bone quantity and quality, and its association with characterization of cellular components of bone. Nine rats (female, 6 months old) for each bilateral sham and ovariectomy (OVX) surgery were obtained and maintained for 2 months after surgery. A hemi-mandible and a femur from each rat were characterized for parameters of volume, mineral density, cortical and trabecular morphology, and static and dynamic mechanical analysis. Another set of 5 rats (female, 9 months old) was obtained for assays of BMSCs. Following cytometry to identify BMSCs, bioassays for proliferation, and osteogenic, adipogenic, chondrogenic differentiation, and cell mitochondrial stress tests were performed. In addition, mRNA expression of BMSCs was analyzed. OVX decreased bone quantity and quality (mineral content, morphology, and energy dissipation) of femur while those of mandible were not influenced. Cellular assays demonstrated that mandible BMSCs showed greater differentiation than femur BMSCs. Gene ontology pathway analysis indicated that the mandibular BMSCs showed most significant differential expression of genes in the regulatory pathways of osteoblast differentiation, SMAD signaling, cartilage development, and glucose transmembrane transporter activity. These findings suggested that active mandibular BMSCs maintain bone formation and mineralization by balancing the rapid bone resorption caused by estrogen deficiency. These characteristics likely help reduce the risk of osteoporotic fracture in postmenopausal jawbone.

Characterizing the Mechanical Behavior of Bone and Bone Surrogates in Compression Using pQCT

Many axial and appendicular skeleton bones are subjected to repetitive loading during daily activities. Until recently, the structural analysis of fractures has been limited to 2D sections, and the dynamic assessment of fracture progression has not been possible. The structural failure was analyzed using step-wise micro-compression combined with time-lapsed micro-computed tomographic imaging. The structural failure was investigated in four different sample materials (two different bone surrogates, lumbar vertebral bodies from bovine and red deer). The samples were loaded in different force steps based on uniaxial compression tests. The micro-tomography images were used to create three-dimensional models from which various parameters were calculated that provide information about the structure and density of the samples. By superimposing two 3D images and calculating the different surfaces, it was possible to precisely analyze which trabeculae failed in which area and under which load. According to the current state of the art, bone mineral density is usually used as a value for bone quality, but the question can be raised as to whether other values such as trabecular structure, damage accumulation, and bone mineralization can predict structural competence better than bone mineral density alone.

Mechanical metric for skeletal biomechanics derived from spectral analysis of stiffness matrix

A new metric for the quantitative and qualitative evaluation of bone stiffness is introduced. It is based on the spectral decomposition of stiffness matrix computed with finite element method. The here proposed metric is defined as an amplitude rescaled eigenvalues of stiffness matrix. The metric contains unique information on the principal stiffness of bone and reflects both bone shape and material properties. The metric was compared with anthropometrical measures and was tested for sex sensitivity on pelvis bone. Further, the smallest stiffness of pelvis was computed under a certain loading condition and analyzed with respect to sex and direction. The metric complements anthropometrical measures and provides a unique information about the smallest bone stiffness independent from the loading configuration and can be easily computed by state-of-the-art subject specified finite element algorithms.

Nonlinear voxel-based finite element model for strength assessment of healthy and metastatic proximal femurs

Nonlinear finite element (FE) models can accurately quantify bone strength in healthy and metastatic femurs. However, their use in clinical practice is limited since state-of-the-art implementations using tetrahedral meshes involve a lot of manual work for which specific modelling software and engineering knowledge are required. Voxel-based meshes could enable the transition since they are robust and can be highly automated. Therefore, the aim of this work was to bridge the modelling gap between the tetrahedral and voxel-based approach. Specifically, we validated a nonlinear voxel-based FE method relative to experimental data from 20 femurs with and without artificial metastases that had been mechanically loaded until failure. CT scans of the femurs were segmented and automatically converted into a voxel-based mesh with hexahedral elements. Nonlinear material properties were implemented in an open-source linear voxel-based FE solver by adding an additional loop to the routine such that the material properties could be adapted after each increment. Bone strength, quantified as the maximum force in the force-displacement curve, was evaluated. The results were compared to a previously established nonlinear tetrahedral FE approach as well as to the experimentally measured bone strength. The voxel-based FE model predicted the experimental bone strength very well both for healthy (R² = 0.90, RMSE = 0.88 kN) and metastatic femurs (R² = 0.93, RMSE = 0.64 kN). The model precision and accuracy were very similar to the ones obtained with the tetrahedral model (R² = 0.90/0.93, RMSE = 0.90/0.64 kN for intact/metastatic respectively). The more intuitive voxel-based meshes thus quantified macroscale femoral strength equally well as state-of-the-art tetrahedral models. The robustness, high level of automation and time-efficiency (< 30 min) of the implemented workflow offer great potential for developing FE models to improve fracture risk prediction in clinical practice.

•

A nonlinear voxel-based FE model was evaluated to assess femoral bone strength

•

Both healthy and metastatic femurs were evaluated

•

The FE models predicted bone strength with high accuracy and precision

•

Voxel-based and tetrahedral FE models showed similar accuracy and precision

•

An iterative routine enabled material nonlinearity in a linear FE solver

Isolated Cyclic Loading During Adolescence Improves Tibial Bone Microstructure and Strength at Adulthood

Bone is a unique living tissue, which responds to the mechanical stimuli regularly imposed on it. Adolescence facilitates a favorable condition for the skeleton that enables the exercise to positively influence bone architecture and overall strength. However, it is still dubious for how long the skeletal benefits gained in adolescence is preserved at adulthood. The current study aims to use a rat model to investigate the effects of in vivo low- (LI), medium- (MI), and high- (HI) intensity cyclic loadings applied during puberty on longitudinal bone development, morphometry, and biomechanics during adolescence as well as at adulthood. Forty-two young (4-week-old) male rats were randomized into control, sham, LI, MI, and HI groups. After a 5 day/week for 8 weeks cyclic loading regime applied on the right tibia, loaded rats underwent a subsequent 41-week, normal cage activity period. Right tibias were removed at 52 weeks of age, and a comprehensive assessment was performed using μCT, mechanical testing, and finite element analysis. HI and MI groups exhibited reduced body weight and food intake at the end of the loading period compared with shams, but these effects disappeared afterward. HI cyclic loading increased BMD, bone volume fraction, trabecular thickness, trabecular number, and decreased trabecular spacing after loading. All loading-induced benefits, except BMD, persisted until the end of the normal cage activity period. Moreover, HI loading induced enhanced bone area, periosteal perimeter, and moment of inertia, which remained up to the 52nd week. After the normal cage activity at adulthood, the HI group showed increased ultimate force and stress, stiffness, postyield displacement and energy, and toughness compared with the sham group. Overall, our findings suggest that even though both trabecular and cortical bone drifted through age-related changes during aging, HI cyclic loading performed during adolescence can render lifelong benefits in bone microstructure and biomechanics. © 2020 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

Age- and Strain-Related Differences in Bone Microstructure and Body Composition During Development in Inbred Male Mouse Strains

We explored age- and strain-related differences in bone microstructure and body composition in male C57BL/6J, DBA/2JRj and C3H/J mice. Bone microstructure of the femur, tibia and L4 was assessed by μCT at the age of 8, 16 and 24 weeks. The weight of several muscles and fat depots were measured at the same time points. At all timepoints, C3H/J mice had the thickest cortices followed by DBA/2JRj and C57BL/6J mice. Nevertheless, C57BL/6J mice had higher Tb.BV/TV and Tb.N, and lower Tb.Sp than DBA/2JRj and C3H/J mice at least at 24 weeks of age. Skeletal development patterns differed among strains. C57BL/6J and DBA/2JRj mice, but not C3H/J mice, experienced significant increases in the sum of the masses of 6 individual muscles by 24 weeks of age. In C57BL/6J and DBA/2JRj mice, the mass of selected fat depots reached highest values at 24 weeks, whist, in C3H/J mice, the highest values of fat depots masses were achieved at 16 weeks. Early strain differences in muscle and fat masses were largely diminished by 24 weeks of age. C3H/J and C57BL/6J mice displayed the most favorable cortical and trabecular bone parameters, respectively. Strain differences in body composition were less overt than strain specificity in bone microstructure, however, they possibly influenced aspects of skeletal development.

High Impact Exercise Improves Bone Microstructure and Strength in Growing Rats

Physical activity is beneficial for skeletal development. However, impact sports during adolescence, leading to bone growth retardation and/or bone quality improvement, remains unexplained. This study investigated the effects of in vivo low (LI), medium (MI), and high (HI) impact loadings applied during puberty on bone growth, morphometry and biomechanics using a rat model. 4-week old rats (n = 30) were divided into control, sham, LI, MI, and HI groups. The impact was applied on the right tibiae, 5 days/week for 8 weeks mimicking walking (450 µε), uphill running (850 µε) and jumping (1250 µε) conditions. Trabecular and cortical parameters were determined by micro-CT, bone growth rate by calcein labeling and toluidine blue staining followed by histomorphometry. Bio-mechanical properties were evaluated from bending tests. HI group reduced rat body weight and food consumption compared to shams. Bone growth rate also decreased in MI and HI groups despite developing thicker hypertrophic and proliferative zone heights. HI group showed significant increment in bone mineral density, trabecular thickness, cortical and total surface area. Ultimate load and stiffness were also increased in MI and HI groups. We conclude that impact loading during adolescence reduces bone growth moderately but improves bone quality and biomechanics at the end of the growing period.

Methods for quantitative characterization of bone injury from computed-tomography images

Computed tomography (CT) images can potentially provide insights into bone structure for diagnosis of disorders and diseases. However, evaluation of trabecular bone structure and whole bone shape is often qualitative or semi-quantitative. This limits inter-study comparisons and the ability to detect subtle bone quality variations during early disease onset or in response to new treatments. In this work, we enable quantitative characterization of bone diseases through bone morphometry, texture analysis, and shape analysis methods. The potential of our analysis methods to identify the impact of hemophilia is validated in a mouse femur wound model. In our results, shape localizes and characterizes the formation of spurious bone, and our texture and bone morphometry analysis results provide extra information about the composition of that bone. Some of our one-dimensional (1D) textural features were able to significantly differentiate our injured femurs from our healthy femurs, even with this small sample size demonstrating the potential of the proposed analysis framework. While trabecular bone morphometrics have been a pillar in 3D microCT bone research for decades, the proposed analysis framework augments how we define and understand phenotypical presentation of bone disease. The contributed open source software is exposed to the medical image analysis community through 3D Slicer extensions to ensure both robustness and reproducibility.

Mechanical Characterization of Bone: State of the Art in Experimental Approaches-What Types of Experiments Do People Do and How Does One Interpret the Results?

PURPOSE OF REVIEW

The mechanical integrity of bone is determined by the direct measurement of bone mechanical properties. This article presents an overview of the current, most common, and new and upcoming experimental approaches for the mechanical characterization of bone. The key outcome variables of mechanical testing, as well as interpretations of the results in the context of bone structure and biology are also discussed.

RECENT FINDINGS

Quasi-static tests are the most commonly used for determining the resistance to structural failure by a single load at the organ (whole bone) level. The resistance to crack initiation or growth by fracture toughness testing and fatigue loading offers additional and more direct characterization of tissue material properties. Non-traditional indentation techniques and in situ testing are being increasingly used to probe the material properties of bone ultrastructure. Destructive ex vivo testing or clinical surrogate measures are considered to be the gold standard for estimating fracture risk. The type of mechanical test used for a particular investigation depends on the length scale of interest, where the outcome variables are influenced by the interrelationship between bone structure and composition. Advancement in the sensitivity of mechanical characterization techniques to detect changes in bone at the levels subjected to modifications by aging, disease, and/or pharmaceutical treatment is required. As such, a number of techniques are now available to aid our understanding of the factors that contribute to fracture risk.

Improving results in rat fracture models: enhancing the efficacy of biomechanical testing by a modification of the experimental setup

Background

Animal fracture models, primarily performed in rats, are crucial to investigate normal and pathological bone healing. However, results of biomechanical testing representing a major outcome measure show high standard deviations often precluding statistical significance. Therefore, the aim of our study was a systematical examination of biomechanical characteristics of rat femurs during three-point bending. Furthermore, we tried to reduce variation of results by individually adapting the span of bearing and loading areas to the bone’s length.

Methods

We examined 40 paired femurs of male Wistar-rats by DXA (BMD and BMC of the whole femur) and pQCT-scans at the levels of bearing and loading areas of the subsequent biomechanical three-point bending test. Individual adjustment of bearing and loading bars was done respecting the length of each specimen. Subgroups of light (< 400 g, n = 22) and heavy (> 400 g, n = 18) animals were formed and analysed separately. We furthermore compared the results of the individualised bending-setting to 20 femurs tested with a fix span of 15 mm.

Results

Femurs showed a length range of 34 to 46 mm. The failure loads ranged from 116 to 251 N (mean 175.4 ± 45.2 N; heavy animals mean 221 ± 18.9 N; light animals mean 138.1 ± 16.4 N) and stiffness ranged from 185 N/mm to 426 N/mm (mean 315.6 ± 63 N/mm; heavy animals mean 358.1 ± 34.64 N/mm; light animals mean 280.8 ± 59.85 N/mm). The correlation of densitometric techniques and failure loads was high (DXA R² = 0.89 and pQCT R² = 0.88). In comparison to femurs tested with a fix span, individual adaptation of biomechanical testing homogenized our data significantly. Most notably, the standard deviation of failure loads (221 ± 18.95 N individualized setting vs. 205.5 ± 30.36 N fixed) and stiffness (358.1 ± 34.64 N/mm individualized setting vs. 498.5 ± 104.8 N/mm fixed) was reduced by at least one third.

Conclusions

Total variation observed in any trait reflects biological and methodological variation. Precision of the method hence affects the statistical power of the study. By simply adapting the setting of the biomechanical testing, interindividual variation could be reduced, which improves the precision of the method significantly.

Electronic supplementary material

The online version of this article (10.1186/s12891-018-2155-y) contains supplementary material, which is available to authorized users.

Relationships of bone characteristics in MYO9B deficient femurs

The objective of this study was to examine relationships among a variety of bone characteristics, including volumetric, mineral density, geometric, dynamic mechanical analysis, and static fracture mechanical properties. As MYO9B is an unconventional myosin in bone cells responsible for normal skeletal growth, bone characteristics of wild-type (WT), heterozygous (HET), and MYO9B knockout (KO) mice groups were compared as an animal model to express different bone quantity and quality. Forty-five sex-matched 12-week-old mice were used in this study. After euthanization, femurs were isolated and scanned using microcomputed tomography (micro-CT) to assess bone volumetric, tissue mineral density (TMD), and geometric parameters. Then, a non-destructive dynamic mechanical analysis (DMA) was performed by applying oscillatory bending displacement on the femur. Finally, the same femur was subject to static fracture testing. KO group had significantly lower length, bone mineral density (BMD), bone mass and volume, dynamic and static stiffness, and strength than WT and HET groups (p < 0.019). On the other hand, TMD parameters of KO group were comparable with those of WT group. HET group showed volumetric, geometric, and mechanical properties similar to WT group, but had lower TMD (p < 0.014). Non-destructive micro-CT and DMA parameters had significant positive correlations with strength (p < 0.015) without combined effect of groups and sex on the correlations (p > 0.077). This comprehensive characterization provides a better understanding of interactive behavior between the tissue- and organ-level of the same femur. The current findings elucidate that MYO9B is responsible for controlling bone volume to determine the growth rate and fracture risk of bone.

Whole bone testing in small animals: systematic characterization of the mechanical properties of different rodent bones available for rat fracture models

Objectives

Rat fracture models are extensively used to characterize normal and pathological bone healing. Despite, systematic research on inter- and intra-individual differences of common rat bones examined is surprisingly not available. Thus, we studied the biomechanical behaviour and radiological characteristics of the humerus, the tibia and the femur of the male Wistar rat—all of which are potentially available in the experimental situation—to identify useful or detrimental biomechanical properties of each bone and to facilitate sample size calculations.

Methods

40 paired femura, tibiae and humeri of male Wistar rats (10–38 weeks, weight between 240 and 720 g) were analysed by DXA, pQCT scan and three-point-bending. Bearing and loading bars of the biomechanical setup were adapted percentually to the bone’s length. Subgroups of light (skeletal immature) rats under 400 g (N = 11, 22 specimens of each bone) and heavy (mature) rats over 400 g (N = 9, 18 specimens of each bone) were formed and evaluated separately.

Results

Radiologically, neither significant differences between left and right bones, nor a specific side preference was evident. Mean side differences of the BMC were relatively small (1–3% measured by DXA and 2.5–5% by pQCT). Over all, bone mineral content (BMC) assessed by DXA and pQCT (TOT CNT, CORT CNT) showed high correlations between each other (BMC vs. TOT and CORT CNT: R² = 0.94–0.99). The load–displacement diagram showed a typical, reproducible curve for each type of bone. Tibiae were the longest bones (mean 41.8 ± 4.12 mm) followed by femurs (mean 38.9 ± 4.12 mm) and humeri (mean 29.88 ± 3.33 mm). Failure loads and stiffness ranged from 175.4 ± 45.23 N / 315.6 ± 63.00 N/mm for the femurs, 124.6 ± 41.13 N / 260.5 ± 59.97 N/mm for the humeri to 117.1 ± 33.94 N / 143.8 ± 36.99 N/mm for the tibiae. Smallest interindividual differences were observed in failure loads of the femurs (CV% 8.6) and tibiae (CV% 10.7) of heavy animals, light animals showed good consistency in failure loads of the humeri (CV% 7.7). Most consistent results of both sides (left vs. right) in failure loads were provided by the femurs of light animals (mean difference 4.0 ± 2.8%); concerning stiffness, humeri of heavy animals were most consistent (mean difference of 6.2 ± 5%). In general, the failure loads showed strong correlations to the BMC (R² = 0.85–0.88) whereas stiffness correlated only moderate, except for the humerus (BMC vs. stiffness: R² = 0.79).

Discussion

Altogether, the rat’s femur of mature specimens showed the most accurate and consistent radiological and biomechanical results. In synopsis with the common experimental use enabling comparison among different studies, this bone offers ideal biomechanical conditions for three point bending experiments. This can be explained by the combination of a superior aspect ratio and a round and long, straight morphology, which satisfies the beam criteria more than other bones tested.

Electronic supplementary material

The online version of this article (10.1186/s40001-018-0307-z) contains supplementary material, which is available to authorized users.

Detection of bone loss via subchondral bone analysis

To date, there is no single sign, symptom, or test that can clearly diagnose early stages of Temporomandibular Joint Osteoarthritis (TMJ OA). However, it has been observed that changes in the bone occur in early stages of this disease, involving structural changes both in the texture and morphometry of the bone marrow and the subchondral cortical plate. In this paper we present a tool to detect and highlight subtle variations in subchondral bone structure obtained from high resolution Cone Beam Computed Tomography (hr-CBCT) in order to help with detecting early TMJ OA. The proposed tool was developed in ITK and 3DSlicer and it has been disseminated as open-source software tools. We have validated both our texture analysis and morphometry analysis biomarkers for detection of TMJ OA comparing hr-CBCT to μCT. Our initial statistical results using the multidimensional features computed with our tool indicate that it is possible to classify areas of demonstrated loss of trabecular bone in both μCT and hr-CBCT. This paper describes the first steps to alleviate the current inability of radiological changes to diagnose TMJ OA before morphological changes are too advanced by quantifying subchondral bone biomarkers. This paper indicates that texture based and morphometry based biomarkers have the potential to identify OA patients at risk for further bone destruction.

The RhoGAP Myo9b Promotes Bone Growth by Mediating Osteoblastic Responsiveness to IGF-1

The Ras homolog A (RhoA) subfamily of Rho guanosine triphosphatases (GTPases) regulates actin-based cellular functions in bone such as differentiation, migration, and mechanotransduction. Polymorphisms or genetic ablation of RHOA and some of its regulatory guanine exchange factors (GEFs) have been linked to poor bone health in humans and mice, but the effects of RhoA-specific GTPase-activating proteins (GAPs) on bone quality have not yet been identified. Therefore, we examined the consequences of RhoGAP Myo9b gene knockout on bone growth, phenotype, and cellular activity. Male and female mice lacking both alleles demonstrated growth retardation and decreased bone formation rates during early puberty. These mice had smaller, weaker bones by 4 weeks of age, but only female KOs had altered cellular numbers, with fewer osteoblasts and more osteoclasts. By 12 weeks of age, bone quality in KOs worsened. In contrast, 4-week-old heterozygotes demonstrated bone defects that resolved by 12 weeks of age. Throughout, Myo9b ablation affected females more than males. Osteoclast activity appeared unaffected. In primary osteogenic cells, Myo9b was distributed in stress fibers and focal adhesions, and its absence resulted in poor spreading and eventual detachment from culture dishes. Similarly, MC3T3-E1 preosteoblasts with transiently suppressed Myo9b levels spread poorly and contained decreased numbers of focal adhesions. These cells also demonstrated reduced ability to undergo IGF-1-induced spreading or chemotaxis toward IGF-1, though responses to PDGF and BMP-2 were unaffected. IGF-1 receptor (IGF1R) activation was normal in cells with diminished Myo9b levels, but the activated receptor was redistributed from stress fibers and focal adhesions into nuclei, potentially affecting receptor accessibility and gene expression. These results demonstrate that Myo9b regulates a subset of RhoA-activated processes necessary for IGF-1 responsiveness in osteogenic cells, and is critical for normal bone formation in growing mice. © 2017 American Society for Bone and Mineral Research.

Fully automated segmentation of callus by micro-CT compared to biomechanics

Background

A high percentage of closed femur fractures have slight comminution. Using micro-CT (μCT), multiple fragment segmentation is much more difficult than segmentation of unfractured or osteotomied bone. Manual or semi-automated segmentation has been performed to date. However, such segmentation is extremely laborious, time-consuming and error-prone. Our aim was to therefore apply a fully automated segmentation algorithm to determine μCT parameters and examine their association with biomechanics.

Methods

The femura of 64 rats taken after randomised inhibitory or neutral medication, in terms of the effect on fracture healing, and controls were closed fractured after a Kirschner wire was inserted. After 21 days, μCT and biomechanical parameters were determined by a fully automated method and correlated (Pearson’s correlation).

Results

The fully automated segmentation algorithm automatically detected bone and simultaneously separated cortical bone from callus without requiring ROI selection for each single bony structure. We found an association of structural callus parameters obtained by μCT to the biomechanical properties. However, results were only explicable by additionally considering the callus location.

Conclusions

A large number of slightly comminuted fractures in combination with therapies that influence the callus qualitatively and/or quantitatively considerably affects the association between μCT and biomechanics. In the future, contrast-enhanced μCT imaging of the callus cartilage might provide more information to improve the non-destructive and non-invasive prediction of callus mechanical properties. As studies evaluating such important drugs increase, fully automated segmentation appears to be clinically important.

Development of an efficient screening system to identify novel bone metabolism-related genes using the exchangeable gene trap mutagenesis mouse models

Despite numerous genetic studies on bone metabolism, understanding of the specific mechanisms is lacking. We developed an efficient screening system to identify novel genes involved in bone metabolism using mutant mouse strains registered with the Exchangeable Gene Trap Clones (EGTC) database. From 1278 trap clones in the EGTC database, 52 candidate lines were selected in the first screening, determined based on "EST profile", "X-gal", "Related article", and "Novel gene". For the second screening, bone morphometric analysis, biomechanical strength analysis, bone X-gal staining, etc. were performed on candidate lines. Forty-two male trap lines (80.8%) showed abnormalities with either bone morphometric analysis or biomechanical strength analysis. In the screening process, X-gal staining was significantly efficient (P = 0.0057). As examples, Lbr and Nedd4 trap lines selected using the screening system showed significant bone decrease and fragility, suggesting a relationship with osteoblast differentiation. This screening system using EGTC mouse lines is extremely efficient for identifying novel genes involved in bone metabolism. The gene trap lines identified as abnormal using this screening approach are highly likely to trap important genes for bone metabolism. These selected trap mice will be valuable for use as novel bio-resources in bone research.

Follistatin-like 3 is a mediator of exercise-driven bone formation and strengthening

Exercise is vital for maintaining bone strength and architecture. Follistatin-like 3 (FSTL3), a member of follistatin family, is a mechanosensitive protein upregulated in response to exercise and is involved in regulating musculoskeletal health. Here, we investigated the potential role of FSTL3 in exercise-driven bone remodeling. Exercise-dependent regulation of bone structure and functions was compared in mice with global Fstl3 gene deletion (Fstl3-/-) and their age-matched Fstl3+/+ littermates. Mice were exercised by low-intensity treadmill walking. The mechanical properties and mineralization were determined by μCT, three-point bending test and sequential incorporation of calcein and alizarin complexone. ELISA, Western-blot analysis and qRT-PCR were used to analyze the regulation of FSTL3 and associated molecules in the serum specimens and tissues. Daily exercise significantly increased circulating FSTL3 levels in mice, rats and humans. Compared to age-matched littermates, Fstl3-/- mice exhibited significantly lower fracture tolerance, having greater stiffness, but lower strain at fracture and yield energy. Furthermore, increased levels of circulating FSTL3 in young mice paralleled greater strain at fracture compared to the lower levels of FSTL3 in older mice. More significantly, Fstl3-/- mice exhibited loss of mechanosensitivity and irresponsiveness to exercise-dependent bone formation as compared to their Fstl3+/+ littermates. In addition, FSTL3 gene deletion resulted in loss of exercise-dependent sclerostin regulation in osteocytes and osteoblasts, as compared to Fstl3+/+ osteocytes and osteoblasts, in vivo and in vitro. The data identify FSTL3 as a critical mediator of exercise-dependent bone formation and strengthening and point to its potential role in bone health and in musculoskeletal diseases.

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.

Bone morphology in 46 BXD recombinant inbred strains and femur-tibia correlation

We examined the bone properties of BXD recombinant inbred (RI) mice by analyzing femur and tibia and compared their phenotypes of different compartments. 46 BXD RI mouse strains were analyzed including progenitor C57BL/6J (n = 16) and DBA/2J (n = 15) and two first filial generations (D2B6F1 and B6D2F1). Strain differences were observed in bone quality and structural properties (P < 0.05) in each bone profile (whole bone, cortical bone, or trabecular bone). It is well known that skeletal phenotypes are largely affected by genetic determinants and genders, such as bone mineral density (BMD). While genetics and gender appear expectedly as the major determinants of bone mass and structure, significant correlations were also observed between femur and tibia. More importantly, positive and negative femur-tibia associations indicated that genetic makeup had an influence on skeletal integrity. We conclude that (a) femur-tibia association in bone morphological properties significantly varies from strain to strain, which may be caused by genetic differences among strains, and (b) strainwise variations were seen in bone mass, bone morphology, and bone microarchitecture along with bone structural property.

Modeling microdamage behavior of cortical bone

Bone is a complex material which exhibits several hierarchical levels of structural organization. At the submicron-scale, the local tissue porosity gives rise to discontinuities in the bone matrix which have been shown to influence damage behavior. Computational tools to model the damage behavior of bone at different length scales are mostly based on finite element (FE) analysis, with a range of algorithms developed for this purpose. Although the local mechanical behavior of bone tissue is influenced by microstructural features such as bone canals and osteocyte lacunae, they are often not considered in FE damage models due to the high computational cost required to simulate across several length scales, i.e., from the loads applied at the organ level down to the stresses and strains around bone canals and osteocyte lacunae. Hence, the aim of the current study was twofold: First, a multilevel FE framework was developed to compute, starting from the loads applied at the whole bone scale, the local mechanical forces acting at the micrometer and submicrometer level. Second, three simple microdamage simulation procedures based on element removal were developed and applied to bone samples at the submicrometer-scale, where cortical microporosity is included. The present microdamage algorithm produced a qualitatively analogous behavior to previous experimental tests based on stepwise mechanical compression combined with in situ synchrotron radiation computed tomography. Our results demonstrate the feasibility of simulating microdamage at a physiologically relevant scale using an image-based meshing technique and multilevel FE analysis; this allows relating microdamage behavior to intracortical bone microstructure.

Differential effects of collagen prolyl 3-hydroxylation on skeletal tissues

Mutations in the genes encoding cartilage associated protein (CRTAP) and prolyl 3-hydroxylase 1 (P3H1 encoded by LEPRE1) were the first identified causes of recessive Osteogenesis Imperfecta (OI). These proteins, together with cyclophilin B (encoded by PPIB), form a complex that 3-hydroxylates a single proline residue on the α1(I) chain (Pro986) and has cis/trans isomerase (PPIase) activity essential for proper collagen folding. Recent data suggest that prolyl 3-hydroxylation of Pro986 is not required for the structural stability of collagen; however, the absence of this post-translational modification may disrupt protein-protein interactions integral for proper collagen folding and lead to collagen over-modification. P3H1 and CRTAP stabilize each other and absence of one results in degradation of the other. Hence, hypomorphic or loss of function mutations of either gene cause loss of the whole complex and its associated functions. The relative contribution of losing this complex's 3-hydroxylation versus PPIase and collagen chaperone activities to the phenotype of recessive OI is unknown. To distinguish between these functions, we generated knock-in mice carrying a single amino acid substitution in the catalytic site of P3h1 (Lepre1^H662A). This substitution abolished P3h1 activity but retained ability to form a complex with Crtap and thus the collagen chaperone function. Knock-in mice showed absence of prolyl 3-hydroxylation at Pro986 of the α1(I) and α1(II) collagen chains but no significant over-modification at other collagen residues. They were normal in appearance, had no growth defects and normal cartilage growth plate histology but showed decreased trabecular bone mass. This new mouse model recapitulates elements of the bone phenotype of OI but not the cartilage and growth phenotypes caused by loss of the prolyl 3-hydroxylation complex. Our observations suggest differential tissue consequences due to selective inactivation of P3H1 hydroxylase activity versus complete ablation of the prolyl 3-hydroxylation complex.

The prolyl 3-hydroxylase complex serves to hydroxylate a single residue in type I collagen and also serves as a collagen chaperone. The complex is comprised of prolyl 3-hydroxylase 1, cartilage associated protein, and cyclophilin B. Mutations have been identified in the genes encoding the complex members in patients with recessive Osteogenesis Imperfecta. Recent data suggest that prolyl 3-hydroxylation of collagen does not alter the stability of collagen but may rather mediate protein-protein interactions. Additionally, the collagen chaperoning function of the complex is an important rate limiting step in the modification of type I collagen. Irrespective of whether patients with mutations in the genes encoding the members of the prolyl 3-hydroxylase complex have hypomorphic or complete loss of function alleles, either circumstance leads to the loss of both functions of the prolyl 3-hydroxylase complex. Thus, it is unknown how collagen chaperoning and/or hydroxylation affect bone and cartilage homeostasis. In this study, we generated a mouse model lacking the prolyl 3-hydroxylation activity of the complex while maintaining the chaperoning function. We found that the hydroxylase mutant mice have normal cartilage and normal cortical bone but decreased trabecular bone, suggesting that there is a differential requirement for hydroxylation in different tissues.

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.

Effects of extensive circumferential periosteal stripping on the microstructure and mechanical properties of the murine femoral cortex

Extensive periosteal stripping (PS) is a risk factor for post-radiation pathologic fracture following surgery for extremity soft tissue tumors. The purpose of this study was to determine the effects of PS on bone structure and mechanical properties. Thirty-one skeletally mature mice underwent PS, with circumferential removal of periosteum from an 8-mm segment of the mid-diaphysis of the left femur. Thirty-one control mice underwent sham surgery in which the femur was isolated without manipulation of the periosteum. At 2, 6, 12, or 26 weeks following surgery, the left femora were examined by micro-CT to quantify cortical thickness (CtTh), cross-sectional area (CSA), bone volume (BV), and polar moment of inertia (PMI). Three-point mechanical bend testing was performed and peak load, stiffness, and energy to failure were determined. PS resulted in significantly decreased CtTh, CSA, BV, and PMI at all time points. Peak load, stiffness, and energy to failure were significantly reduced at 2, 6, and 12 weeks. There were no significant differences in mechanical properties at 26 weeks. In this mouse model, extensive circumferential PS resulted in sustained changes in bone structure that were still evident after 6 months, accompanied by reductions in bone strength that persisted for at least 3 months.

Mouse tail vertebrae adapt to cyclic mechanical loading by increasing bone formation rate and decreasing bone resorption rate as shown by time-lapsed in vivo imaging of dynamic bone morphometry

It is known that mechanical loading leads to an increase in bone mass through a positive shift in the balance between bone formation and bone resorption. How the remodeling sites change over time as an effect of loading remains, however, to be clarified. The purpose of this paper was to investigate how bone formation and resorption sites are modulated by mechanical loading over time by using a new imaging technique that extracts three dimensional formation and resorption parameters from time-lapsed in vivo micro-computed tomography images. To induce load adaptation, the sixth caudal vertebra of C57BL/6 mice was cyclically loaded through pins in the adjacent vertebrae at either 8 N or 0 N (control) three times a week for 5 min (3000 cycles) over a total of 4 weeks. The results showed that mechanical loading significantly increased trabecular bone volume fraction by 20% (p<0.001) and cortical area fraction by 6% (p<0.001). The bone formation rate was on average 23% greater (p<0.001) and the bone resorption rate was on average 25% smaller (p<0.001) for the 8 N group than for the 0 N group. The increase in bone formation rate for the 8 N group was mostly an effect of a significantly increased surface of bone formation sites (on average 16%, p<0.001), while the thickness of bone formation packages was less affected (on average 5% greater, p<0.05). At the same time the surface of bone resorption sites was significantly reduced (on average 15%, p<0.001), while the depth of resorption pits remained the same. For the 8 N group, the strength of the whole bone increased significantly by 24% (p<0.001) over the loading period, while the strain energy density in the trabecular bone decreased significantly by 24% (p<0.001). In conclusion, mouse tail vertebrae adapt to mechanical loading by increasing the surface of formation sites and decreasing the surface of resorption sites, leading to an overall increase in bone strength. This new imaging technique will provide opportunities to investigate in vivo bone remodeling in the context of disease and treatment options, with the added value that both bone formation and bone resorption parameters can be nondestructively calculated over time.

The role of GH/IGF-I-mediated mechanisms in sex differences in cortical bone size in mice

Cortical bone dimensions are important determinants of bone strength. Gender differences in cortical bone size caused by greater periosteal expansion in males than in females during the pubertal growth spurt are well established both in humans and in experimental animal models. However, the mechanism by which gender influences cortical bone size is still a matter of investigation. The role of androgens and estrogen in pubertal bone growth has been examined in human disorders as well as animal models, such as gonadectomized or sex steroid receptor knockout mice. Based on the findings that growth hormone (GH) and insulin-like growth factor I (IGF-I) are major regulators of postnatal skeletal growth, we and others have predicted that sex hormones interact with the GH/IGF-I axis to regulate cortical bone size. However, studies conflict as to whether estrogen and androgens impact cortical bone size through the canonical pathway, through GH without IGF-I mediation, through IGF-I without GH stimulation, or independent of GH/IGF-I. We review recent data on the impact of sex steroids and components of the GH/IGF axis on sexual dimorphism in bone size. While the GH/IGF-I axis is a major player in regulating peak bone size, the relative contribution of GH/IGF-dependent mechanisms to sex differences in cortical bone size remains to be established.

Swimming training increases the post-yield energy of bone in young male rats

The purpose of this study is to investigate the effects of non-weight-bearing exercise on growing bone. Male Wistar rats (7 week-old) were assigned to one baseline control group, one control group and two swimming training groups, which were trained with 2 and 4% body-weight mass added, respectively. After an 8-week training period, three groups showed significant development compared to the baseline control group. Among the three 15-week-old groups, swimming-trained rats were lower in body weight (BW), densitometry and size-related measurements. In femoral biomechanical testing, swimming training groups were significantly lower in yield moment and ultimate moment, which may be due to a significantly lower long bone cross-sectional moment of inertia. However, the two swimming groups were higher in post-yield energy absorption and displacement. Further, in estimated tissue-level biomaterial properties, no differences were shown in yield stress, strain or toughness among the three groups. Using BW as a covariate, results of ANCOVA showed no differences in size-related parameters among the three groups, and some parameters were even higher in the two swimming groups. Regarding Pearson's correlation, size-related parameters correlated well to BW and whole bone strength but not to tissue post-yield behaviors. In conclusion, when compared to age-matched control group, swimming rats showed lower bone strength and lower yield energy absolutely at the structural level, but similar yield stress and yield toughness at the tissue level. Moreover, swimming training benefited growing bone in post-yield behaviors. Further studies should investigate the parameters that contribute to this exercise-induced post-yield behavior.

Hierarchical microimaging of bone structure and function

With recent advances in molecular medicine and disease treatment in osteoporosis, quantitative image processing of three-dimensional bone structures is critical in the context of bone quality assessment. Biomedical imaging technology such as MRI or CT is readily available, but few attempts have been made to expand the capabilities of these systems by integrating quantitative analysis tools and by exploring structure-function relationships in a hierarchical fashion. Nevertheless, such quantitative end points are an important factor for success in basic research and in the development of novel therapeutic strategies. CT is key to these developments, as it images and quantifies bone in three dimensions and provides multiscale biological imaging capabilities with isotropic resolutions of a few millimeters (clinical CT), a few tens of micrometers (microCT) and even as high as 100 nanometers (nanoCT). The technology enables the assessment of the relationship between microstructural and ultrastructural measures of bone quality and certain diseases or therapies. This Review focuses on presenting strategies for three-dimensional approaches to hierarchical biomechanical imaging in the study of microstructural and ultrastructural bone failure. From this Review, it can be concluded that biomechanical imaging is extremely valuable for the study of bone failure mechanisms at different hierarchical levels.