Longitudinal effects of Parathyroid Hormone treatment on morphological, densitometric and mechanical properties of mouse tibia

Combining PTH(1-34) and mechanical loading has increased benefit to tibia bone mechanics in ovariectomised mice

Combined treatment with PTH(1-34) and mechanical loading confers increased structural benefits to bone than monotherapies. However, it remains unclear how this longitudinal adaptation affects the bone mechanics. This study quantified the individual and combined longitudinal effects of PTH(1-34) and mechanical loading on the bone stiffness and strength evaluated in vivo with validated micro-finite element (microFE) models. C57BL/6 mice were ovariectomised at 14-week-old and treated either with injections of PTH(1-34), compressive tibia loading or both interventions concurrently. Right tibiae were in vivo microCT-scanned every 2 weeks from 14 until 24-week-old. MicroCT images were rigidly registered to reference tibia and the cortical organ level (whole bone) and tissue level (midshaft) morphometric properties and bone mineral content were quantified. MicroCT images were converted into voxel-based homogeneous, linear elastic microFE models to estimate the bone stiffness and strength. This approach allowed us for the first time to quantify the longitudinal changes in mechanical properties induced by combined treatments in a model of accelerated bone resorption. Both changes of stiffness and strength were higher with co-treatment than with individual therapies, consistent with increased benefits with the tibia bone mineral content and cortical area, properties strongly associated with the tibia mechanics. The longitudinal data shows that the two bone anabolics, both individually and combined, had persistent benefit on estimated mechanical properties, and that benefits (increased stiffness and strength) remained after treatment was withdrawn.

The loading direction dramatically affects the mechanical properties of the mouse tibia

Introduction: The in vivo tibial loading mouse model has been extensively used to evaluate bone adaptation in the tibia after mechanical loading treatment. However, there is a prevailing assumption that the load is applied axially to the tibia. The aim of this in silico study was to evaluate how much the apparent mechanical properties of the mouse tibia are affected by the loading direction, by using a validated micro-finite element (micro-FE) model of mice which have been ovariectomized and exposed to external mechanical loading over a two-week period. Methods: Longitudinal micro-computed tomography (micro-CT) images were taken of the tibiae of eleven ovariectomized mice at ages 18 and 20 weeks. Six of the mice underwent a mechanical loading treatment at age 19 weeks. Micro-FE models were generated, based on the segmented micro-CT images. Three models using unitary loads were linearly combined to simulate a range of loading directions, generated as a function of the angle from the inferior-superior axis (θ, 0°-30° range, 5° steps) and the angle from the anterior-posterior axis (ϕ, 0°: anterior axis, positive anticlockwise, 0°-355° range, 5° steps). The minimum principal strain was calculated and used to estimate the failure load, by linearly scaling the strain until 10% of the nodes reached the critical strain level of -14,420 με. The apparent bone stiffness was calculated as the ratio between the axial applied force and the average displacement along the longitudinal direction, for the loaded nodes. Results: The results demonstrated a high sensitivity of the mouse tibia to the loading direction across all groups and time points. Higher failure loads were found for several loading directions (θ = 10°, ϕ 205°-210°) than for the nominal axial case (θ = 0°, ϕ = 0°), highlighting adaptation of the bone for loading directions far from the nominal axial one. Conclusion: These results suggest that in studies which use mouse tibia, the loading direction can significantly impact the failure load. Thus, the magnitude and direction of the applied load should be well controlled during the experiments.

PTH Treatment Increases Cortical Bone Mass More in Response to Compression than Tension in Mice

Parathyroid hormone (PTH) is an anabolic osteoporosis treatment that increases bone mass and reduces fracture risk. Clinically, the effects of PTH are site-specific, increasing bone mass more at the spine than the hip and not increasing bone mass at the radius. Differences in local loading environment between the spine, hip, and radius may help explain the variation in efficacy, as PTH and mechanical loading have been shown to synergistically increase bone mass. We hypothesized that differences in loading mode might further explain these variations. Owing to the curvature of the mouse tibia, cyclic compression of the hindlimb causes bending at the tibial midshaft, placing the anterior surface under tension and the posterior surface under compression. We investigated the combination of PTH treatment and tibial loading in an osteoblast-specific estrogen receptor-alpha knockout mouse model of low bone mass (pOC-ERαKO) and their littermate controls (LCs) and analyzed bone morphology in the tensile, compressive, and neutral regions of the tibial midshaft. We also hypothesized that pretreating wild-type C57Bl/6J (WT) mice with PTH prior to mechanical loading would enhance the synergistic anabolic effects. Compression was more anabolic than tension, and PTH enhanced the effect of loading, particularly under compression. PTH pretreatment maintained the synergistic anabolic effect for longer durations than concurrent treatment and loading alone. Together these data provide insights into more effective physical therapy and exercise regimens for patients receiving PTH treatment. © 2022 American Society for Bone and Mineral Research (ASBMR).

Beam theory for rapid strain estimation in the mouse tibia compression model

The mouse tibia compression model is a leading model for studying bone's mechanoadaptive response to load. In studying this mechanoadaptive response, (FE) modelling is often used to determine the stress/strain within the tibia. The development of such models can be challenging and computationally expensive. An alternate approach is to use continuum mechanics based analytical theories, such as beam theory (BT). However, applying BT to the mouse tibia requires the fibula be neglected, introducing error in the stress/strain distribution. While several studies have applied BT to the mouse tibia, no study has explored the accuracy of this approach. To address these questions, this work investigates the use of BT in determining stress/strain within the mouse tibia. By comparing BT against FE modelling, it was found that BT can accurately predict tibial stress/strain if correction factors are applied to account for the effect of the fibula. The 25, 37, 50 and 75% cross sections are studied. Focusing on the 37% cross section, without correction, BT can have errors of approximately 21.6%. With correction, this is reduced to 6.6%. Such correction factors are presented. The developed BT model is applicable in the diaphysis and distal metaphysis, where the assumptions of BT are valid. This work verifies BT for determining localised strains in a mouse tibia compression model. This is anticipated to provide efficiency dividends, allowing for high throughput modelling of the mouse tibia, advancing study of bone's mechanoadaptive response.

Reproducibility of Densitometric and Biomechanical Assessment of the Mouse Tibia From In Vivo Micro-CT Images

Interventions for bone diseases (e.g. osteoporosis) require testing in animal models before clinical translation and the mouse tibia is among the most common tested anatomical sites. In vivo micro-Computed Tomography (microCT) based measurements of the geometrical and densitometric properties are non-invasive and therefore constitute an important tool in preclinical studies. Moreover, validated micro-Finite Element (microFE) models can be used for predicting the bone mechanical properties non-invasively. However, considering that the image processing pipeline requires operator-dependant steps, the reproducibility of these measurements has to be assessed. The aim of this study was to evaluate the intra- and inter-operator reproducibility of several bone parameters measured from microCT images. Ten in vivo microCT images of the right tibia of five mice (at 18 and 22 weeks of age) were processed. One experienced operator (intra-operator analysis) and three different operators (inter-operator) aligned each image to a reference through a rigid registration and selected a volume of interest below the growth plate. From each image the following parameters were measured: total bone mineral content (BMC) and density (BMD), BMC in 40 subregions (ten longitudinal sections, four quadrants), microFE-based stiffness and failure load. Intra-operator reproducibility was acceptable for all parameters (precision error, PE < 3.71%), with lowest reproducibility for stiffness (3.06% at week 18, 3.71% at week 22). The inter-operator reproducibility was slightly lower (PE < 4.25%), although still acceptable for assessing the properties of most interventions. The lowest reproducibility was found for BMC in the lateral sector at the midshaft (PE = 4.25%). Densitometric parameters were more reproducible than most standard morphometric parameters calculated in the proximal trabecular bone. In conclusion, microCT and microFE models provide reproducible measurements for non-invasive assessment of the mouse tibia properties.

Positive interactions of mechanical loading and PTH treatments on spatio-temporal bone remodelling

Osteoporosis is one of the most common skeletal diseases, but current therapies are limited to generalized antiresorptive or anabolic interventions, which do not target regions that would benefit from improvements to skeletal health. To improve the evaluation of treatment plans, we used a spatio-temporal multiscale approach that combines longitudinal in vivo micro-computed tomography (micro-CT) and in silico subject-specific finite element modeling to quantitatively map bone adaptation changes due to disease and treatment at high resolution. Our findings show time and region-dependent modifications in bone remodelling following one and two sets of mechanical loading and/or pharmacological interventions. The multiscale results highlighted that the distal section was unaffected by mechanical loading alone but the proximal tibia had the greatest gain from positive interactions of combined therapies. Mechanical loading abated the catabolic effect of PTH, but the main benefit of combined treatments occurred from the additive interactions of the two therapies in periosteal apposition. These results provide detailed insight into the efficacy of combined treatments, facilitating the optimisation of dosage and treatment duration in preclinical mouse studies, and the development of novel interventions for skeletal diseases. STATEMENT OF SIGNIFICANCE: Combined mechanical loading and pharmacotherapy have the potential to slow osteoporosis-induced bone loss but current therapies do not target the regions in need of strengthening. We show for the first time spatial region-dependant interactions between PTH and mechanical loading treatment in OVX mouse tibiae, highlighting local regions in the tibia that benefitted from separate and combined treatments. Combined experimental-computational analysis also detailed the lasting period of each treatment per location in the tibia, the extent of positive (or negative) interactions of the combined therapies, and the impact of each treatment on the regulation of bone adaptation spatio-temporally. This approach can be used to create hypothesis about the interactions of different treatments to optimise the design of biomaterials and medical interventions.

The Role of the Loading Condition in Predictions of Bone Adaptation in a Mouse Tibial Loading Model

The in vivo mouse tibial loading model is used to evaluate the effectiveness of mechanical loading treatment against skeletal diseases. Although studies have correlated bone adaptation with the induced mechanical stimulus, predictions of bone remodeling remained poor, and the interaction between external and physiological loading in engendering bone changes have not been determined. The aim of this study was to determine the effect of passive mechanical loading on the strain distribution in the mouse tibia and its predictions of bone adaptation. Longitudinal micro-computed tomography (micro-CT) imaging was performed over 2 weeks of cyclic loading from weeks 18 to 22 of age, to quantify the shape change, remodeling, and changes in densitometric properties. Micro-CT based finite element analysis coupled with an optimization algorithm for bone remodeling was used to predict bone adaptation under physiological loads, nominal 12N axial load and combined nominal 12N axial load superimposed to the physiological load. The results showed that despite large differences in the strain energy density magnitudes and distributions across the tibial length, the overall accuracy of the model and the spatial match were similar for all evaluated loading conditions. Predictions of densitometric properties were most similar to the experimental data for combined loading, followed closely by physiological loading conditions, despite no significant difference between these two predicted groups. However, all predicted densitometric properties were significantly different for the 12N and the combined loading conditions. The results suggest that computational modeling of bone's adaptive response to passive mechanical loading should include the contribution of daily physiological load.

There is direct relationship between bone bridge length and coracoclavicular fixation resistance to failure: Biomechanical study in a porcine model

BACKGROUND

This study aims to evaluate the relation between coracoclavicular resistance to failure and the distance between clavicular tunnels. The hypothesis is that a greater clavicular bone bridge between tunnels achieves a stronger coracoclavicular fixation.

METHODS

Descriptive Laboratory Study. Thirty-six (36) coracoclavicular models were constructed utilizing porcine metatarsals. Coracoclavicular stabilizations were performed using a subcoracoid loop fixation configuration through two clavicular tunnels, tied at the clavicle's superior cortex using a locking knot. Models were randomly assigned to 1 of 3 experimental groups of variable bone bridge length between clavicular tunnels: 5 mm, 10 mm, and 15 mm. Each group had 12 models. Fixation resistance was assessed through the ultimate failure point under an axial load to failure trial. Failure patterns were documented. A one-way ANOVA test was used, and a Tukey post hoc as needed (P < 0.05).

FINDINGS

Mean strength per bone bridge length: 5 mm = 312 N (Range: 182-442 N); 10 mm = 430 N (Range: 368-595 N); 15 mm = 595 N (Range: 441-978 N). The 15 mm group had a significantly higher ultimate failure point than the other two groups: 5 mm (P < 0.001) and 10 mm (P < 0.001). All fixations systematically failed by a superior cortex clavicle fracture at the midpoint between tunnels.

INTERPRETATION

A direct relationship between bone bridge length and coracoclavicular resistance to failure was demonstrated, being the 15 mm length a significantly higher strength construct in a tied loop model.

High resolution imaging in bone tissue research-review

This review article focuses on imaging of bone tissue to understand skeletal health with regards to bone quality. Skeletal fragility fractures are due to bone diseases such as osteoporosis which result in low bone mass and bone mineral density (BMD) leading to high risk of fragility fractures. Recent advances in imaging and analysis technologies have highly benefitted the field of biological sciences. In particular, their application in skeletal health has been of significant importance in understanding bone mechanical behavior (structure and properties) at the tissue level. While synchrotron based microCT technique has remained the gold standard for non-destructive evaluation of structure in material and biological sciences, several lab based microCT systems have been developed to provide high resolution imaging of specimens with greater access, and ease of use in laboratory settings. Lab based microCT scanners are widely used in the bone field as a standard tool to evaluate three-dimensional (3D) morphologies of bone structure at image resolutions appropriate for bone samples from small animals to bone biopsy specimens from humans. Both synchrotron and standard lab based microCT systems provide high resolution imaging ex vivo for a small sized specimen. A few X-ray based systems are also commercially available for in vivo scanning at relatively low image resolutions. Synchrotron-based CT microscopy is being used for various ultra-high-resolution image analyses using complex 3D software. However, the synchrotron-based CT technology is in high demand, allows only limited numbers of specimens, expensive, requires complex additional instrumentation, and is not easily available to researchers as it requires access to a synchrotron source which is always limited. Therefore, desktop laboratory scanners (microXCT, Zeiss/Xradia, Scanco, SkyScan. etc.), mimicking the synchrotron based CT technology or image resolution, have been developed to solve the accessibility issues. These lab based scanners have helped both material science, and the bone field to investigate bone tissue morphologies at submicron mage resolutions. Considerable progress has been made in both in vivo and ex vivo imaging towards providing high resolution images of bone tissue. Both clinical and research imaging technologies will continue to improve and help understand osteoporosis and other related skeletal issues in order to develop targeted treatments for bone fragility. This review summarizes the high resolution imaging work in bone research.

Optimization of the failure criterion in micro-Finite Element models of the mouse tibia for the non-invasive prediction of its failure load in preclinical applications

New treatments against osteoporosis require testing in animal models and the mouse tibia is among the most common studied anatomical sites. In vivo micro-Computed Tomography (microCT) based micro-Finite Element (microFE) models can be used for predicting the bone strength non-invasively, after proper validation against experiments. The aim of this study was to evaluate the ability of different microCT-based bone parameters and microFE models to predict tibial structural mechanical properties in compression. Twenty tibiae were scanned at 10.4 μm voxel size and subsequently tested in uniaxial compression at 0.03 mm/s until failure. Stiffness and failure load were measured from the load-displacement curves. Standard morphometric parameters were measured from the microCT images. The spatial distribution of bone mineral content (BMC) was evaluated by dividing the tibia into 40 regions. MicroFE models were generated by converting each microCT image into a voxel-based mesh with homogeneous isotropic material properties. Failure load was estimated by using different failure criteria, and the optimized parameters were selected by minimising the errors with respect to experimental measurements. Experimental and predicted stiffness were moderately correlated (R² = 0.65, error = 14% ± 8%). Normalized failure load was best predicted by microFE models (R² = 0.81, error = 9% ± 6%). Failure load was not correlated to the morphometric parameters and weakly correlated with some geometrical parameters (R² < 0.37). In conclusion, microFE models can improve the current estimation of the mouse tibia structural properties and in this study an optimal failure criterion has been defined. Since it is a non-invasive method, this approach can be applied longitudinally for evaluating temporal changes in the bone strength.

Bone remodelling in the mouse tibia is spatio-temporally modulated by oestrogen deficiency and external mechanical loading: A combined in vivo/in silico study

Osteoporosis disrupts the healthy remodelling process in bone and affects its mechanical properties. Mechanical loading has been shown to be effective in stimulating bone formation to mitigate initial bone loss. However, no study has investigated the effects of repeated mechanical loading, with a pause of one week in between, in the mouse tibia with oestrogen deficiency. This study uses a combined experimental and computational approach, through longitudinal monitoring with micro-computed tomography, to evaluate the effects of loading on bone adaptation in the tibiae of ovariectomised (OVX) C57BL/6 mice from 14 to 22 weeks of age. Micro-FE models coupled with bone adaptation algorithms were used to estimate changes in local tissue strains due to OVX and mechanical loading, and to quantify the relationship between local strain and remodelling. The first in vivo mechanical loading increased apposition, by 50-150%, while resorption decreased by 50-60%. Both endosteal and periosteal resorption increased despite the second mechanical loading, and periosteal resorption was up to 70% higher than that after the first loading. This was found to correlate with an initial decrease in average strain energy density after the first loading, which was lower and more localised after the second loading. Predictions of bone adaptation showed that between 50 and 90% of the load-induced bone apposition is linearly strain driven at the organ-level, but resorption is more biologically driven at the local level. The results imply that a systematic increase in peak load or loading rate may be required to achieve a similar bone adaptation rate in specific regions of interests.

In silico clinical trials: concepts and early adoptions

Innovations in information and communication technology infuse all branches of science, including life sciences. Nevertheless, healthcare is historically slow in adopting technological innovation, compared with other industrial sectors. In recent years, new approaches in modelling and simulation have started to provide important insights in biomedicine, opening the way for their potential use in the reduction, refinement and partial substitution of both animal and human experimentation. In light of this evidence, the European Parliament and the United States Congress made similar recommendations to their respective regulators to allow wider use of modelling and simulation within the regulatory process. In the context of in silico medicine, the term 'in silico clinical trials' refers to the development of patient-specific models to form virtual cohorts for testing the safety and/or efficacy of new drugs and of new medical devices. Moreover, it could be envisaged that a virtual set of patients could complement a clinical trial (reducing the number of enrolled patients and improving statistical significance), and/or advise clinical decisions. This article will review the current state of in silico clinical trials and outline directions for a full-scale adoption of patient-specific modelling and simulation in the regulatory evaluation of biomedical products. In particular, we will focus on the development of vaccine therapies, which represents, in our opinion, an ideal target for this innovative approach.

PTH(1-34) treatment and/or mechanical loading have different osteogenic effects on the trabecular and cortical bone in the ovariectomized C57BL/6 mouse

In preclinical mouse models, a synergistic anabolic response to PTH(1–34) and tibia loading was shown. Whether combined treatment improves bone properties with oestrogen deficiency, a cardinal feature of osteoporosis, remains unknown. This study quantified the individual and combined longitudinal effects of PTH(1–34) and loading on the bone morphometric and densitometric properties in ovariectomised mice. C57BL/6 mice were ovariectomised at 14-weeks-old and treated either with injections of PTH(1–34); compressive loading of the right tibia; both interventions concurrently; or both interventions on alternating weeks. Right tibiae were microCT-scanned from 14 until 24-weeks-old. Trabecular metaphyseal and cortical midshaft morphometric properties, and bone mineral content (BMC) in 40 different regions of the tibia were measured. Mice treated only with loading showed the highest trabecular bone volume fraction at week 22. Cortical thickness was higher with co-treatment than in the mice treated with PTH alone. In the mid-diaphysis, increases in BMC were significantly higher with loading than PTH. In ovariectomised mice, the osteogenic benefits of co-treatment on the trabecular bone were lower than loading alone. However, combined interventions had increased, albeit regionally-dependent, benefits to cortical bone. Increased benefits were largest in the mid-diaphysis and postero-laterally, regions subjected to higher strains under compressive loads.

Mechanical regulation of bone formation and resorption around implants in a mouse model of osteopenic bone

Although mechanical stimulation is considered a promising approach to accelerate implant integration, our understanding of load-driven bone formation and resorption around implants is still limited. This lack of knowledge may delay the development of effective loading protocols to prevent implant loosening, especially in osteoporosis. In healthy bone, formation and resorption are mechanoregulated processes. In the intricate context of peri-implant bone regeneration, it is not clear whether bone (re)modelling can still be load-driven. Here, we investigated the mechanical control of peri-implant bone (re)modelling with a well-controlled mechanobiological experiment. We applied cyclic mechanical loading after implant insertion in tail vertebrae of oestrogen depleted mice and we monitored peri-implant bone response by in vivo micro-CT. Experimental data were combined with micro-finite element simulations to estimate local tissue strains in (re)modelling locations. We demonstrated that a substantial increase in bone mass around the implant could be obtained by loading the entire bone. This augmentation could be attributed to a large reduction in bone resorption rather than to an increase in bone formation. We also showed that following implantation, mechanical regulation of bone (re)modelling was transiently lost. Our findings should help to clarify the role of mechanical stimulation on the maintenance of peri-implant bone mass.

Regional Nanoindentation Properties in Different Locations on the Mouse Tibia From C57BL/6 and Balb/C Female Mice

The local spatial heterogeneity of the material properties of the cortical and trabecular bone extracted from the mouse tibia is not well-known. Nevertheless, its characterization is fundamental to be able to study comprehensively the effect of interventions and to generate computational models to predict the bone strength preclinically. The goal of this study was to evaluate the nanoindentation properties of bone tissue extracted from two different mouse strains across the tibia length and in different sectors. Left tibiae were collected from four female mice, two C57BL/6, and two Balb/C mice. Nanoindentations with maximum 6 mN load were performed on different microstructures, regions along the axis of the tibiae, and sectors (379 in total). Reduced modulus (Er) and hardness (H) were computed for each indentation. Trabecular bone of Balb/C mice was 21% stiffer than that of C57BL/6 mice (20.8 ± 4.1 GPa vs. 16.5 ± 7.1 GPa). Moreover, the proximal regions of the bones were 13-36% less stiff than the mid-shaft and distal regions of the same bones. No significant differences were found for the different sectors for E _r and H for Balb/C mice. The bone in the medial sector was found to be 8-14% harder and stiffer than the bone in the anterior or posterior sectors for C57BL/6 mice. In conclusion, this study showed that the nanoindentation properties of the mouse tibia are heterogeneous across the tibia length and the trabecular bone properties are different between Balb/C and C57BL/6 mice. These results will help the research community to identify regions where to characterize the mechanical properties of the bone during preclinical optimisation of treatments for skeletal diseases.

ChronoMID-Cross-modal neural networks for 3-D temporal medical imaging data

ChronoMID-neural networks for temporally-varying, hence Chrono, Medical Imaging Data-makes the novel application of cross-modal convolutional neural networks (X-CNNs) to the medical domain. In this paper, we present multiple approaches for incorporating temporal information into X-CNNs and compare their performance in a case study on the classification of abnormal bone remodelling in mice. Previous work developing medical models has predominantly focused on either spatial or temporal aspects, but rarely both. Our models seek to unify these complementary sources of information and derive insights in a bottom-up, data-driven approach. As with many medical datasets, the case study herein exhibits deep rather than wide data; we apply various techniques, including extensive regularisation, to account for this. After training on a balanced set of approximately 70000 images, two of the models-those using difference maps from known reference points-outperformed a state-of-the-art convolutional neural network baseline by over 30pp (> 99% vs. 68.26%) on an unseen, balanced validation set comprising around 20000 images. These models are expected to perform well with sparse data sets based on both previous findings with X-CNNs and the representations of time used, which permit arbitrarily large and irregular gaps between data points. Our results highlight the importance of identifying a suitable description of time for a problem domain, as unsuitable descriptors may not only fail to improve a model, they may in fact confound it.

A novel algorithm to predict bone changes in the mouse tibia properties under physiological conditions

Understanding how bone adapts to mechanical stimuli is fundamental for optimising treatments against musculoskeletal diseases in preclinical studies, but the contribution of physiological loading to bone adaptation in mouse tibia has not been quantified so far. In this study, a novel mechanistic model to predict bone adaptation based on physiological loading was developed and its outputs were compared with longitudinal scans of the mouse tibia. Bone remodelling was driven by the mechanical stimuli estimated from micro-FEA models constructed from micro-CT scans of C57BL/6 female mice (N = 5) from weeks 14 and 20 of age, to predict bone changes in week 16 or 22. Parametric analysis was conducted to evaluate the sensitivity of the models to subject-specific or averaged parameters, parameters from week 14 or week 20, and to strain energy density (SED) or maximum principal strain (ε_maxprinc). The results at week 20 showed no significant difference in bone densitometric properties between experimental and predicted images across the tibia for both stimuli, and 59% and 47% of the predicted voxels matched with the experimental sites in apposition and resorption, respectively. The model was able to reproduce regions of bone apposition in both periosteal and endosteal surfaces (70% and 40% for SED and ε_maxprinc, respectively), but it under-predicted the experimental sites of resorption by over 85%. This study shows for the first time the potential of a subject-specific mechanoregulation algorithm to predict bone changes in a mouse model under physiological loading. Nevertheless, the weak predictions of resorption suggest that a combined stimulus or biological stimuli should be accounted for in the model.

Effect of repeated in vivo microCT imaging on the properties of the mouse tibia

In longitudinal studies, in vivo micro-Computed Tomography (microCT) imaging is used to investigate bone changes over time due to interventions in mice. However, ionising radiation can provoke significant variations in bone morphometric parameters. In a previous study, we evaluated the effect of reducing the integration time on the properties of the mouse tibia measured from microCT images. A scanning procedure (100 ms integration time, 256 mGy nominal radiation dose) was selected as the best compromise between image quality and radiation dose induced on the animal. In this work, the effect of repeated in vivo scans has been evaluated using the selected procedure. The right tibia of twelve female C57BL/6 (six wild type, WT, six ovariectomised, OVX) and twelve BALB/c (six WT, six OVX) mice was scanned every two weeks, starting at week 14 of age. At week 24, mice were sacrificed and both tibiae were scanned. Standard trabecular and cortical morphometric parameters were calculated. The spatial distribution of densitometric parameters (e.g. bone mineral content) was obtained by dividing each tibia in 40 partitions. Stiffness and strength in compression were estimated using homogeneous linear elastic microCT-based micro-Finite Element models. Differences between right (irradiated) and left (non-irradiated control) tibiae were evaluated for each parameter. The irradiated tibiae had higher Tb.Th (+3.3%) and Tb.Sp (+11.6%), and lower Tb.N (-14.2%) compared to non-irradiated tibiae, consistently across both strains and intervention groups. A reduction in Tb.BV/TV (-14.9%) was also observed in the C57BL/6 strain. In the OVX group, a small reduction was also observed in Tt.Ar (-5.0%). In conclusion, repeated microCT scans (at 256 mGy, 5 scans, every two weeks) had limited effects on the mouse tibia, compared to the expected changes induced by bone treatments. Therefore, the selected scanning protocol is acceptable for measuring the effect of bone interventions in vivo.

The longitudinal effects of ovariectomy on the morphometric, densitometric and mechanical properties in the murine tibia: A comparison between two mouse strains

Oestrogen deficiency-related bone loss in the ovariectomized (OVX) mouse is a common model for osteoporosis. However, a comprehensive in vivo assessment of intervention-related changes in multiple bone properties, and in multiple mouse strains, is required in order to identify an appropriate model for future evaluation of novel anti-osteoporotic therapies. The aim of this study was to evaluate the effect of OVX on the morphometric and densitometric properties measured in the microCT images and the mechanical properties estimated with finite element models of the tibia in two mouse strains, C57BL/6 and BALB/c. 14-weeks-old female C57BL/6 and BALB/c mice were divided into two groups per strain: (1) ovariectomized, (2) non-operated control. The right tibia was scanned at baseline (14 weeks) and then every two weeks thereafter, until 24-weeks-old, using in vivo microCT. Changes in trabecular and cortical bone morphometry, spatiotemporal changes in densitometric properties and in mechanical properties (from micro-finite element (μFE) analysis) were computed. Differences between OVX and non-operated controls were evaluated by ANCOVA, adjusted for 14-weeks baseline. In morphometry, trabecular bone mass was significantly reduced in both C57BL/6 and BALB/c from four weeks following surgery. Though the OVX-effect was transient in BALB/c as bone mass reached skeletal homeostasis. OVX inhibited the age-related thickening of cortical bone only in C57BL/6. In both strains, increments in bone mineral content were significantly lower with OVX only in the proximal tibia, with intervention-related differences increasing with time. OVX had no effect on μFE estimates of stiffness nor failure load in either strain. The results of this study show strain-, time- and region-(trabecular or cortical) dependent changes in morphometric and densitometric properties. These findings highlight the importance of choosing an appropriate mouse model and time points for research of treatments against accelerated bone resorption.

A new method to monitor bone geometry changes at different spatial scales in the longitudinal in vivo μCT studies of mice bones

Longitudinal studies of bone adaptation in mice using in vivo micro-computed tomography (μCT) have been commonly used for pre-clinical evaluation of physical and pharmacological interventions. The main advantage of this approach is to use each mouse as its own control, reducing considerably the sample size required by statistical power analysis. To date, multi-scale estimation of bone adaptations become essential since the bone activity that takes place at different scales may be associated with different bone mechanisms. Measures of bone adaptations at different time scales have been attempted in a previous study. This paper extends quantification of bone activity at different spatial scales with a proposition of a novel framework. The method involves applying level-set method (LSM) to track the geometric changes from the longitudinal in vivo μCT scans of mice tibia. Bone low- and high-spatial frequency patterns are then estimated using multi-resolution analysis. The accuracy of the framework is quantified by applying it to two times separated scanned images with synthetically manipulated global and/or local activity. The Root Mean Square Deviation (RMSD) was approximately 1.5 voxels or 0.7 voxels for the global low-spatial frequency or local high-spatial frequency changes, respectively. The framework is further applied to the study of bone changes in longitudinal datasets of wild-type mice tibiae over time and space. The results demonstrate the ability for the spatio-temporal quantification and visualisation of bone activity at different spatial scales in longitudinal studies thus providing further insight into bone adaptation mechanisms.

Effect of parathyroid hormone on the structural, densitometric and failure behaviors of mouse tibia in the spatiotemporal space

Parathyroid hormone (PTH) is an anabolic bone drug approved by the US Food and Drug Administration (FDA) to treat osteoporosis. However, previous studies using cross-sectional designs have reported variable and sometimes contradictory results. The aim of the present study was to quantify the localized effect of PTH on the structural and densitometric behaviors of mouse tibia and their links with the global mechanical behavior of bone using a novel spatiotemporal image analysis approach and a finite element analysis technique. Twelve female C57BL/6J mice were divided into two groups: the control and PTH treated groups. The entire right tibiae were imaged using an in vivo micro-computed tomography (μCT) system eight consecutive times. Next, the in vivo longitudinal tibial μCT images were rigidly registered and divided into 10 compartments across the entire tibial space. The bone volume (BV), bone mineral content (BMC), bone tissue mineral density (TMD), and tibial endosteal and periosteal areas (TEA and TPA) were quantified in each compartment. Additionally, finite element models of all the tibiae were generated to analyze the failure behavior of the tibia. It was found that both the BMC and BV started to increase in the proximal tibial region, and then the increases extended to the entire tibial region after two weeks of treatment (p < 0.05). PTH intervention significantly reduced the TEA in most tibial compartments after two weeks of treatment, and the TPA increased in most tibial regions after four weeks of treatment (p < 0.05). Tibial failure loads significantly increased after three weeks of PTH treatment (p < 0.01). The present study provided the first evidence of the localized effect of PTH on bone structural and densitometric properties, as well as their links with the global mechanical behaviors of bone, which are important pieces of information for unveiling the mechanism of PTH intervention.

Stochastic analysis of a heterogeneous micro-finite element model of a mouse tibia

Finite element (FE) analysis can be used to predict bone mechanical environments that can be used for many important applications, such as the understanding of bone mechano-regulation mechanisms. However, when defining the FE models, uncertainty in bone material properties may lead to marked variations in the predicted mechanical environment. The aim of this study is to investigate the influence of uncertainty in bone material property on the mechanical environment of bone. A heterogeneous FE model of a mouse tibia was created from micro computed tomography images. Axial compression loading was applied, and all possible bone density-modulus relationships were considered through stochastic analysis. The 1st and 3rd principal strains (ε₁ and ε₃) and the strain energy density (SED) were quantified in the tibial volume of interest (VOI). The bounds of ε₁, ε₃, and SED were determined by the bounds of the density-modulus relationship; the bone mechanical environment (ε₁, ε₃, and SED) and the bone density-modulus relationship exhibit the same trend of change; the relative percentage differences caused by bone material uncertainty are up to 28%, 28%, and 21% for ε₁, ε₃, and SED, respectively. These data provide guidelines on the adoption of bone density-modulus relationship in heterogeneous FE models.

From bed to bench: How in silico medicine can help ageing research

Driven by the raising ethical concerns surrounding animal experimentation, there is a growing interest for non-animal methods, in vitro or in silico technologies that can be used to reduce, refine, and replace animal experimentation. In addition, animal experimentation is being critically revised in regard to its ability to predict clinical outcomes. In this manuscript we describe an initial exploration where a set of in vivo imaging based subject-specific technologies originally developed to predict the risk of femoral strength and hip fracture in osteoporotic patients, were adapted to assess the efficacy of bone drugs pre-clinically on mice. The CT2S technology we developed generates subject-specific models based on Computed Tomography that can separate fractured and non-fractured patients with an accuracy of 82%. When used in mouse experiments the use of in vivo imaging and modelling was found to improve the reproducibility of Bone Mineral Content measurements to a point where up to 63% less mice would be required to achieve the same statistical power of a conventional cross-sectional study. We also speculate about a possible approach where animal-specific and patient-specific models could be used to better translate the observation made on animal models into predictions of response in humans.

Validation of finite element models of the mouse tibia using digital volume correlation

The mouse tibia is a common site to investigate bone adaptation. Micro-Finite Element (microFE) models based on micro-Computed Tomography (microCT) images can estimate bone mechanical properties non-invasively but their outputs need to be validated with experiments. Digital Volume Correlation (DVC) can provide experimental measurements of displacements over the whole bone volume. In this study we applied DVC to validate the local predictions of microFE models of the mouse tibia in compression. Six mouse tibiae were stepwise compressed within a microCT system. MicroCT images were acquired in four configurations with applied compression of 0.5 N (preload), 6.5 N, 13.0 N and 19.5 N. Failure load was measured after the last scan. A global DVC algorithm was applied to the microCT images in order to obtain the displacement field over the bone volume. Homogeneous, isotropic linear hexahedral microFE models were generated from the images collected in the preload configuration with boundary conditions interpolated from the DVC displacements at the extremities of the tibia. Experimental displacements from DVC and numerical predictions were compared at corresponding locations in the middle of the bone. Stiffness and strength were also estimated from each model and compared with the experimental measurements. The magnitude of the displacement vectors predicted by microFE models was highly correlated with experimental measurements (R² >0.82). Higher but still reasonable errors were found for the Cartesian components. The models tended to overestimate local displacements in the longitudinal direction (R² = 0.69-0.90, slope of the regression line=0.50-0.97). Errors in the prediction of structural mechanical properties were 14% ± 11% for stiffness and 9% ± 9% for strength. In conclusion, the DVC approach has been applied to the validation of microFE models of the mouse tibia. The predictions of the models for both structural and local properties have been found reasonable for most preclinical applications.

EVALUATING THE THEORY OF BONE MECHANOREGULATION IN THE PHYSIOLOGICAL LOADING SCENARIO

In this paper, the theory of bone mechanoregulation under physiological loading was evaluated. The entire right tibiae of wild type (WT, [Formula: see text]) and parathyroid hormone (PTH, [Formula: see text]) treated C57BL/6J female mice were scanned using an in vivo [Formula: see text]CT imaging system at 14, 16, 17, 18, 19, 20, 21, and 22 weeks. The PTH intervention started from week 18 until week 22. Subject-specific finite element (FE) models were created from the [Formula: see text]CT images and physiological loading condition was defined in the FE models. The rates of changes in bone mineral content (BMC), bone mineral density (BMD), and bone tissue density (TMD) were quantified over 40 anatomical compartments across the entire mouse tibia. The resulting values were then correlated to the average 1st principal tensile strain ([Formula: see text]) and the strain energy density (SED) for every compartment at weeks 18, 20, and 22. It was found that: in both groups, [Formula: see text] had a minimal effect on the variability of [Formula: see text]BMC ([Formula: see text]); SED had a significant effect on the variability of [Formula: see text]BMC only in the WT group ([Formula: see text]); [Formula: see text] had a significant effect on the variability of [Formula: see text]BMD only in the PTH group ([Formula: see text]); SED had a significant effect on the variability of [Formula: see text]BMD in both groups ([Formula: see text]); neither SED nor [Formula: see text] had a significant effect on the variability of [Formula: see text]TMD ([Formula: see text]). These results are the first to reveal the mechanism of bone mechanoregulation in the physiological loading scenario.