Mineralizing surface is the main target of mechanical stimulation independent of age: 3D dynamic in vivo morphometry

Spatial variations in the osteocyte lacuno-canalicular network density and analysis of the connectomic parameters

Osteocyte lacuno-canalicular network (LCN) is comprised of micrometre-sized pores and submicrometric wide channels in bone. Accumulating evidence suggests multiple functions of this network in material transportation, mechanobiological signalling, mineral homeostasis and bone remodelling. Combining rhodamine staining and confocal laser scanning microscopy, the longitudinal cross-sections of six mouse tibiae were imaged, and the connectome of the network was quantified with a focus on the spatial heterogeneities of network density, connectivity and length of canaliculi. In-vivo loading and double calcein labelling on these tibiae allowed differentiating the newly formed bone from the pre-existing regions. The canalicular density of the murine cortical bone varied between 0.174 and 0.243 μm/μm3, and therefore is three times larger than the corresponding value for human femoral midshaft osteons. The spatial heterogeneity of the network was found distinctly more pronounced across the cortex than along the cortex. We found that in regions with a dense network, the LCN conserves its largely tree-like character, but increases the density by including shorter canaliculi. The current study on healthy mice should serve as a motivating starting point to study the connectome of genetically modified mice, including models of bone diseases and of reduced mechanoresponse.

From breast cancer cell homing to the onset of early bone metastasis: The role of bone (re)modeling in early lesion formation

Breast cancer often metastasizes to bone, causing osteolytic lesions. Structural and biophysical changes are rarely studied yet are hypothesized to influence metastasis. We developed a mouse model of early bone metastasis and multimodal imaging to quantify cancer cell homing, bone (re)modeling, and onset of metastasis. Using tissue clearing and three-dimensional (3D) light sheet fluorescence microscopy, we located enhanced green fluorescent protein-positive cancer cells and small clusters in intact bones and quantified their size and spatial distribution. We detected early bone lesions using in vivo microcomputed tomography (microCT)-based time-lapse morphometry and revealed altered bone (re)modeling in the absence of detectable lesions. With a new microCT image analysis tool, we tracked the growth of early lesions over time. We showed that cancer cells home in all bone compartments, while osteolytic lesions are only detected in the metaphysis, a region of high (re)modeling. Our study suggests that higher rates of (re)modeling act as a driver of lesion formation during early metastasis.

Cortical bone adaptation response is region specific, but not peak load dependent: insights from μ CT image analysis and mechanostat simulations of the mouse tibia loading model

The two major aims of the present study were: (i) quantify localised cortical bone adaptation at the surface level using contralateral endpoint imaging data and image analysis techniques, and (ii) investigate whether cortical bone adaptation responses are universal or region specific and dependent on the respective peak load. For this purpose, we re-analyse previously published μ CT data of the mouse tibia loading model that investigated bone adaptation in response to sciatic neurectomy and various peak load magnitudes (F = 0, 2, 4, 6, 8, 10, 12 N). A beam theory-based approach was developed to simulate cortical bone adaptation in different sections of the tibia, using longitudinal strains as the adaptive stimuli. We developed four mechanostat models: universal, surface-based, strain directional-based, and combined surface and strain direction-based. Rates of bone adaptation in these mechanostat models were computed using an optimisation procedure (131,606 total simulations), performed on a single load case (F = 10 N). Subsequently, the models were validated against the remaining six peak loads. Our findings indicate that local bone adaptation responses are quasi-linear and bone region specific. The mechanostat model which accounted for differences in endosteal and periosteal regions and strain directions (i.e. tensile versus compressive) produced the lowest root mean squared error between simulated and experimental data for all loads, with a combined prediction accuracy of 76.6, 55.0 and 80.7% for periosteal, endosteal, and cortical thickness measurements (in the midshaft of the tibia). The largest root mean squared errors were observed in the transitional loads, i.e. F = 2 to 6 N, where inter-animal variability was highest. Finally, while endpoint imaging studies provide great insights into organ level bone adaptation responses, the between animal and loaded versus control limb variability make simulations of local surface-based adaptation responses challenging.

Zoledronate reduces loading-induced microdamage in cortical ulna of ovariectomized rats

As a daily physiological mechanism in bone, microdamage accumulation dissipates energy and helps to prevent fractures. However, excessive damage accumulation might bring adverse effects to bone mechanical properties, which is especially problematic among the osteoporotic and osteopenic patients treated by bisphosphonates. Some pre-clinical studies in the literature applied forelimb loading models to produce well-controlled microdamage in cortical bone. Ovariectomized animals were also extensively studied to assimilate human conditions of estrogen-related bone loss. In the present study, we combined both experimental models to investigate microdamage accumulation in the context of osteopenia and zoledronate treatment. Three-month-old normal and ovariectomized rats treated by saline or zoledronate underwent controlled compressive loading on their right forelimb to create in vivo microdamage, which was then quantified by barium sulfate contrast-enhanced micro-CT imaging. Weekly in vivo micro-CT scans were taken to evaluate bone (re)modeling and to capture microstructural changes over time. After sacrifice, three-point-bending tests were performed to assess bone mechanical properties. Results show that the zoledronate treatment can reduce cortical microdamage accumulation in ovariectomized rats, which might be explained by the enhancement of several bone structural properties such as ultimate force, yield force, cortical bone area and volume. The rats showed increased bone formation volume and surface after the generation of microdamage, especially for the normal and the ovariectomized groups. Woven bone formation was also observed in loaded ulnae, which was most significant in ovariectomized rats. Although all the rats showed strong correlations between periosteal bone formation and microdamage accumulation, the correlation levels were lower for the zoledronate-treated groups, potentially because of their lower levels of microdamage. The present study provides insights to further investigations of pharmaceutical treatments for osteoporosis and osteopenia. The same experimental concept can be applied in future studies on microdamage and drug testing.

Impact loading intensifies cortical bone (re)modeling and alters longitudinal bone growth of pubertal rats

Physical exercise is important for musculoskeletal development during puberty, which builds bone mass foundation for later in life. However, strenuous levels of training might bring adverse effects to bone health, reducing longitudinal bone growth. Animal models with various levels of physical exercise were largely used to provide knowledge to clinical settings. Experiments from our previous studies applied different levels of mechanical loading on rat tibia during puberty accompanied by weekly in vivo micro-CT scans. In the present article, we apply 3D image registration-based methods to retrospectively analyze part of the previously acquired micro-CT data. Longitudinal bone growth, growth plate thickness, and cortical bone (re)modeling were evaluated from rats' age of 28-77 days. Our results show that impact loading inhibited proximal bone growth throughout puberty. We hypothesize that impact loading might bring different growth alterations to the distal and proximal growth plates. High impact loading might lead to pathological consequence of osteochondrosis and catch-up growth due to growth inhibition. Impact loading also increased cortical bone (re)modeling before and after the peak proximal bone growth period of young rats, of which the latter case might be caused by the shift from modeling to remodeling as the dominant activity toward the end of rat puberty. We confirm that the tibial endosteum is more mechano-sensitive than the periosteum in response to mechanical loading. To our knowledge, this is the first study to follow up bone growth and bone (re)modeling of young rats throughout the entire puberty with a weekly time interval.

In vivo microCT-based time-lapse morphometry reveals anatomical site-specific differences in bone (re)modeling serving as baseline parameters to detect early pathological events

The bone structure is very dynamic and continuously adapts its geometry to external stimuli by modeling and remodeling the mineralized tissue. In vivo microCT-based time-lapse morphometry is a powerful tool to study the temporal and spatial dynamics of bone (re)modeling. Here an advancement in the methodology to detect and quantify site-specific differences in bone (re)modeling of 12-week-old BALB/c nude mice is presented. We describe our method of quantifying new bone surface interface readouts and how these are influenced by bone curvature. This method is then used to compare bone surface (re)modeling in mice across different anatomical regions to demonstrate variations in the rate of change and spatial gradients thereof. Significant differences in bone (re)modeling baseline parameters between the metaphyseal and epiphyseal, as well as cortical and trabecular bone of the distal femur and proximal tibia are shown. These results are validated using conventional static in vivo microCT analysis. Finally, the insights from these new baseline values of physiological bone (re)modeling were used to evaluate pathological bone (re)modeling in a pilot breast cancer bone metastasis model. The method shows the potential to be suitable to detect early pathological events and track their spatio-temporal development in both cortical and trabecular bone. This advancement in (re)modeling surface analysis and defined baseline parameters according to distinct anatomical regions will be valuable to others investigating various disease models with site-distinct local alterations in bone (re)modeling.

Tissue-Level Regeneration and Remodeling Dynamics are Driven by Mechanical Stimuli in the Microenvironment in a Post-Bridging Loaded Femur Defect Healing Model in Mice

Bone healing and remodeling are mechanically driven processes. While the generalized response to mechanical stimulation in bone is well-understood, much less is known about the mechanobiology-regulating tissue-scale bone formation and resorption during the reparative and remodeling phases of fracture healing. In this study, we combined computational approaches in the form of finite element analysis and experimental approaches by using a loaded femoral defect model in mice to investigate the role of mechanical stimulation in the microenvironment of bone. Specifically, we used longitudinal micro-computed tomography to observe temporal changes in bone at different densities and micro-finite element analysis to map the mechanics of the microenvironment to tissue-scale formation, quiescence (no change in bone presence between time points), and resorption dynamics in the late reparative and remodeling phases (post bridging). Increasing levels of effective strain led to increasing conditional probability of bone formation, while decreasing levels of effective strain led to increasing probability of bone resorption. In addition, the analysis of mineralization dynamics showed both a temporal and effective strain level-dependent behavior. A logarithmic-like response was displayed, where the conditional probability of bone formation or resorption increased rapidly and plateaued or fell rapidly and plateaued as mechanical strain increased.

In vivo and in silico monitoring bone regeneration during distraction osteogenesis of the mouse femur

BACKGROUND AND OBJECTIVE

Distraction osteogenesis (DO) is a mechanobiological process of producing new bone by gradual and controlled distraction of the surgically separated bone segments. Mice have been increasingly used to study the role of relevant biological factors in regulating bone regeneration during DO. However, there remains a lack of in silico DO models and related mechano-regulatory tissue differentiation algorithms for mouse bone. This study sought to establish an in silico model based on in vivo experimental data to simulate the bone regeneration process during DO of the mouse femur.

METHODS

In vivo micro-CT, including time-lapse morphometry was performed to monitor the bone regeneration in the distraction gap. A 2D axisymmetric finite element model, with a geometry originating from the experimental data, was created. Bone regeneration was simulated with a fuzzy logic-based two-stage (distraction and consolidation) mechano-regulatory tissue differentiation algorithm, which was adjusted from that used for fracture healing based on our in vivo experimental data. The predictive potential of the model was further tested with varied distraction frequencies and distraction rates.

RESULTS

The computational simulations showed similar bone regeneration patterns to those observed in the micro-CT data from the experiment throughout the DO process. This consisted of rapid bone formation during the first 10 days of the consolidation phase, followed by callus reshaping via bone remodeling. In addition, the computational model predicted a faster and more robust bone healing response as the model's distraction frequency was increased, whereas higher or lower distraction rates were not conducive to healing.

CONCLUSIONS

This in silico model could be used to investigate basic mechanobiological mechanisms involved in bone regeneration or to optimize DO strategies for potential clinical applications.

Validation of an in vivo micro-CT-based method to quantify longitudinal bone growth of pubertal rats

Bone growth is an essential part of skeletal development during childhood and puberty. Accurately characterizing longitudinal bone growth helps to better understand the determining factors of peak bone mass, which has impacts on bone quality later in life. Animal models were largely used to study longitudinal bone growth. However, the commonly used histology-based method is destructive and unable to follow up the growth curve of live animals in longitudinal experiments. In this study, we validated an in vivo micro-CT-based method against the histology-based method to quantify longitudinal bone growth rates of young rats non-destructively. CD (Sprague Dawley) IGS rats aged 35, 49 and 63 days received the same treatments: two series of repeated in vivo micro-CT scans on their proximal hind limb at a five-day interval, and two calcein injections separated by three days. The longitudinal bone growth rate was quantified by registering time-lapse micro-CT images in 3D, calculating the growth distance on registered images, and dividing the distance by the five-day gap. The growth rate was also evaluated by measuring the 2D distance between consecutive calcein fluorescent bands on microscopic images, divided by the three-day gap. The two methods were both validated independently with reproducible repeated measurements, where the micro-CT-based method showed higher precision. They were also validated against each other with low relative errors and a strong Pearson sample correlation coefficient (0.998), showing a significant (p < 0.0001) linear correlation between paired results. We conclude that the micro-CT-based method can serve as an alternative to the histology-based method for the quantification of longitudinal growth. Thanks to its non-invasive nature and true 3D capability, the micro-CT-based method helps to accommodate in vivo longitudinal animal studies with highly reproducible measurements.

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.

Bone adaptation to mechanical loading in mice is affected by circadian rhythms

Physical forces are critical for successful function of many organs including bone. Interestingly, the timing of exercise during the day alters physiology and gene expression in many organs due to circadian rhythms. Circadian clocks in tissues, such as bone, express circadian clock genes that target tissue-specific genes, resulting in tissue-specific rhythmic gene expression (clock-controlled genes). We hypothesized that the adaptive response of bone to mechanical loading is regulated by circadian rhythms. First, mice were sham loaded and sacrificed 8 h later, which amounted to tissues being collected at zeitgeber time (ZT)2, 6, 10, 14, 18, and 22. Cortical bone of the tibiae collected from these mice displayed diurnal expression of core clock genes and key osteocyte and osteoblast-related genes, such as the Wnt-signaling inhibitors Sost and Dkk1, indicating these are clock-controlled genes. Serum bone turnover markers did not display rhythmicity. Second, mice underwent a single bout of in vivo loading at either ZT2 or ZT14 and were sacrificed 1, 8, or 24 h after loading. Loading at ZT2 resulted in Sost upregulation, while loading at ZT14 led to Sost and Dkk1 downregulation. Third, mice underwent daily in vivo tibial loading over 2 weeks administered either in the morning, (ZT2, resting phase) or evening (ZT14, active phase). In vivo microCT was performed at days 0, 5, 10, and 15 and conventional histomorphometry was performed at day 15. All outcome measures indicated a robust response to loading, but only microCT-based time-lapse morphometry showed that loading at ZT14 resulted in a greater endocortical bone formation response compared to mice loaded at ZT2. The decreased Sost and Dkk1 expression coincident with the modest, but significant time-of-day specific increase in adaptive bone formation, suggests that circadian clocks influence bone mechanoresponse.

Mechanical loading and parathyroid hormone effects and synergism in bone vary by site and modeling/remodeling regime

Intermittent injections of parathyroid hormone (PTH) and mechanical loading are both known to effect a net increase in bone mass. Fundamentally, bone metabolism can be divided into modeling (uncoupled formation or resorption) and remodeling (subsequent formation biologically coupled to resorption in space and time). Methods to delineate the bone response between these regimes are scant but have garnered recent attention and acceptance, and will be critical tools to properly assess short- and long-term efficacy of osteoporosis treatments. To this end, we employ a time-lapse micro-computed tomography strategy to quantify and localize modeling and remodeling volumes over 4 weeks of concurrent PTH treatment and mechanical loading. Modeled and remodeled volumes are probed for differences with respect to treatment, loading, and interactions thereof in trabecular and cortical bone compartments, which were further separated by plate/rod microarchitecture and periosteal/endosteal surfaces, respectively. Loading effects are further considered independently with regard to localized strain environments. Our findings indicate that in trabecular bone, PTH and loading stimulate anabolic modeling additively, and remodeling synergistically. PTH tends to lead to bone accumulation indiscriminate of trabecular microarchitecture, whereas loading tends to more strongly affect plates than rods. The cortical surfaces responded uniquely to PTH and loading, with synergistic effects on the periosteal surface for anabolic modeling, and on the endosteal surface for catabolic modeling. The increase in catabolic modeling due to loading, which is enhanced by PTH, is concentrated to areas of the endosteal surface under low strain and to our knowledge has not previously been reported. Taken together, the effects of PTH, loading, and their interactions, are shown to be dependent on the specific bone compartment and metabolic regime; this may explain some discrepancies in previously-reported findings.

Real-time finite element analysis allows homogenization of tissue scale strains and reduces variance in a mouse defect healing model

Mechanical loading allows both investigation into the mechano-regulation of fracture healing as well as interventions to improve fracture-healing outcomes such as delayed healing or non-unions. However, loading is seldom individualised or even targeted to an effective mechanical stimulus level within the bone tissue. In this study, we use micro-finite element analysis to demonstrate the result of using a constant loading assumption for all mouse femurs in a given group. We then contrast this with the application of an adaptive loading approach, denoted real time Finite Element adaptation, in which micro-computed tomography images provide the basis for micro-FE based simulations and the resulting strains are manipulated and targeted to a reference distribution. Using this approach, we demonstrate that individualised femoral loading leads to a better-specified strain distribution and lower variance in tissue mechanical stimulus across all mice, both longitudinally and cross-sectionally, while making sure that no overloading is occurring leading to refracture of the femur bones.

Non-invasive prediction of the mouse tibia mechanical properties from microCT images: comparison between different finite element models

New treatments for bone diseases require testing in animal models before clinical translation, and the mouse tibia is among the most common models. 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 experimental data. Different modelling techniques can be used to estimate the bone properties, and the accuracy associated with each is unclear. The aim of this study was to evaluate the ability of different microCT-based microFE models to predict the mechanical properties of the mouse tibia under compressive load. Twenty tibiae were microCT scanned at 10.4 µm voxel size and subsequently compressed at 0.03 mm/s until failure. Stiffness and failure load were measured from the load-displacement curves. Different microFE models were generated from each microCT image, with hexahedral or tetrahedral mesh, and homogeneous or heterogeneous material properties. Prediction accuracy was comparable among models. The best correlations between experimental and predicted mechanical properties, as well as lower errors, were obtained for hexahedral models with homogeneous material properties. Experimental stiffness and predicted stiffness were reasonably well correlated (R2 = 0.53-0.65, average error of 13-17%). A lower correlation was found for failure load (R2 = 0.21-0.48, average error of 9-15%). Experimental and predicted mechanical properties normalized by the total bone mass were strongly correlated (R2 = 0.75-0.80 for stiffness, R2 = 0.55-0.81 for failure load). In conclusion, hexahedral models with homogeneous material properties based on in vivo microCT images were shown to best predict the mechanical properties of the mouse tibia.

Age and Sex Differences in Load-Induced Tibial Cortical Bone Surface Strain Maps

Bone adapts its architecture to the applied load; however, it is still unclear how bone mechano‐adaptation is coordinated and why potential for adaptation adjusts during the life course. Previous animal models have suggested strain as the mechanical stimulus for bone adaptation, but yet it is unknown how mouse cortical bone load‐related strains vary with age and sex. In this study, full‐field strain maps (at 1 N increments up to 12 N) on the bone surface were measured in young, adult, and old (aged 10, 22 weeks, and 20 months, respectively), male and female C57BL/6J mice with load applied using a noninvasive murine tibial model. Strain maps indicate a nonuniform strain field across the tibial surface, with axial compressive loads resulting in tension on the medial side of the tibia because of its curved shape. The load‐induced surface strain patterns and magnitudes show sexually dimorphic changes with aging. A comparison of the average and peak tensile strains indicates that the magnitude of strain at a given load generally increases during maturation, with tibias in female mice having higher strains than in males. The data further reveal that postmaturation aging is linked to sexually dimorphic changes in average and maximum strains. The strain maps reported here allow for loading male and female C57BL/6J mouse legs in vivo at the observed ages to create similar increases in bone surface average or peak strain to more accurately explore bone mechano‐adaptation differences with age and sex. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC. on behalf of American Society for Bone and Mineral Research.

Mechanical loading prevents bone destruction and exerts anti-tumor effects in the MOPC315.BM.Luc model of myeloma bone disease

Bone continually adapts to changing external loading conditions via (re)modeling (modeling and remodeling) processes. While physical activity is known to beneficially enhance bone mass in healthy individuals, little is known in how physical stimuli affect osteolytic bone destruction in patients suffering from multiple myeloma bone disease. Multiple myeloma (MM) is caused by malignant plasma cells in the bone marrow, shifting the balance in bone remodeling towards massive resorption. We hypothesized that in vivo tibial mechanical loading has anabolic effects in mice with locally injected MOPC315.BM.Luc cells. Conventional microCT analysis revealed enhanced cortical bone mass and microstructure in loaded compared to nonloaded mice. State-of-the-art time-lapse microCT based image analysis demonstrated bone (re)modeling processes at the endosteal and periosteal surfaces as the underlying causes of increased bone mass. Loading prevented the progression and development of osteolytic destruction. Physical stimuli also diminished local MM cell growth and dissemination evidenced by quantification of MM cell-specific immunoglobulin A levels in the serum of mice and by bioluminescence analysis. These data indicate that mechanical loading not only rescues the bone phenotype, but also exerts cell-extrinsic anti-myeloma effects in the MOPC315.BM.Luc model. In conclusion, the use of physical stimuli should be further investigated as an anabolic treatment for osteolytic bone destruction in patients with MM.

Cortical bone adaptation to a moderate level of mechanical loading in male Sost deficient mice

Loss-of-function mutations in the Sost gene lead to high bone mass phenotypes. Pharmacological inhibition of Sost/sclerostin provides a new drug strategy for treating osteoporosis. Questions remain as to how physical activity may affect bone mass under sclerostin inhibition and if that effect differs between males and females. We previously observed in female Sost knockout (KO) mice an enhanced cortical bone formation response to a moderate level of applied loading (900 με at the tibial midshaft). The purpose of the present study was to examine cortical bone adaptation to the same strain level applied to male Sost KO mice. Strain-matched in vivo compressive loading was applied to the tibiae of 10-, 26- and 52-week-old male Sost KO and littermate control (LC) mice. The effect of tibial loading on bone (re)modeling was measured by microCT, 3D time-lapse in vivo morphometry, 2D histomorphometry and gene expression analyses. As expected, Sost deficiency led to high cortical bone mass in 10- and 26-week-old male mice as a result of increased bone formation. However, the enhanced bone formation associated with Sost deficiency did not appear to diminish with skeletal maturation. An increase in bone resorption was observed with skeletal maturation in male LC and Sost KO mice. Two weeks of in vivo loading (900 με at the tibial midshaft) induced only a mild anabolic response in 10- and 26-week-old male mice, independent of Sost deficiency. A decrease in the Wnt inhibitor Dkk1 expression was observed 3 h after loading in 52-week-old Sost KO and LC mice, and an increase in Lef1 expression was observed 8 h after loading in 10-week-old Sost KO mice. The current results suggest that long-term inhibition of sclerostin in male mice does not influence the adaptive response of cortical bone to moderate levels of loading. In contrast with our previous strain-matched study in females showing enhanced bone responses with Sost ablation, these results in males indicate that the influence of Sost deficiency on the cortical bone formation response to a moderate level of loading differs between males and females. Clinical studies examining antibodies to inhibit sclerostin may need to consider that the efficacy of additional physical activity regimens may be sex dependent.

The mechanoresponse of bone is closely related to the osteocyte lacunocanalicular network architecture

The explanation of how bone senses and adapts to mechanical stimulation still relies on hypotheses. The fluid flow hypothesis claims that a load-induced fluid flow through the lacunocanalicular network can be sensed by osteocytes, which reside within the network structure. We show that considering the network architecture results in a better prediction of bone remodeling than mechanical strain alone. This was done by calculating the fluid flow through the lacunocanalicular network in bone volumes covering the complete cross-sections of mouse tibiae, which underwent controlled in vivo loading. The established relationship between mechanosensitivity and network architecture in individual animals implies possibilities for patient-specific therapies. A new connectomics approach to analyze lacunocanalicular network properties is necessary to understand skeletal mechanobiology.

Organisms rely on mechanosensing mechanisms to adapt to changes in their mechanical environment. Fluid-filled network structures not only ensure efficient transport but can also be employed for mechanosensation. The lacunocanalicular network (LCN) is a fluid-filled network structure, which pervades our bones and accommodates a cell network of osteocytes. For the mechanism of mechanosensation, it was hypothesized that load-induced fluid flow results in forces that can be sensed by the cells. We use a controlled in vivo loading experiment on murine tibiae to test this hypothesis, whereby the mechanoresponse was quantified experimentally by in vivo micro-computed tomography (µCT) in terms of formed and resorbed bone volume. By imaging the LCN using confocal microscopy in bone volumes covering the entire cross-section of mouse tibiae and by calculating the fluid flow in the three-dimensional (3D) network, we could perform a direct comparison between predictions based on fluid flow velocity and the experimentally measured mechanoresponse. While local strain distributions estimated by finite-element analysis incorrectly predicts preferred bone formation on the periosteal surface, we demonstrate that additional consideration of the LCN architecture not only corrects this erroneous bias in the prediction but also explains observed differences in the mechanosensitivity between the three investigated mice. We also identified the presence of vascular channels as an important mechanism to locally reduce fluid flow. Flow velocities increased for a convergent network structure where all of the flow is channeled into fewer canaliculi. We conclude that, besides mechanical loading, LCN architecture should be considered as a key determinant of bone adaptation.

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.

Human and mouse bones physiologically integrate in a humanized mouse model while maintaining species-specific ultrastructure

Humanized mouse models are increasingly studied to recapitulate human-like bone physiology. While human and mouse bone architectures differ in multiple scales, the extent to which chimeric human-mouse bone physiologically interacts and structurally integrates remains unknown. Here, we identify that humanized bone is formed by a mosaic of human and mouse collagen, structurally integrated within the same bone organ, as shown by immunohistochemistry. Combining this with materials science techniques, we investigate the extracellular matrix of specific human and mouse collagen regions. We show that human-like osteocyte lacunar-canalicular network is retained within human collagen regions and is distinct to that of mouse tissue. This multiscale analysis shows that human and mouse tissues physiologically integrate into a single, functional bone tissue while maintaining their species-specific ultrastructural differences. These results offer an original method to validate and advance tissue-engineered human-like bone in chimeric animal models, which grow to be eloquent tools in biomedical research.

Effects of long-term in vivo micro-CT imaging on hallmarks of osteopenia and frailty in aging mice

In vivo micro-CT has already been used to monitor microstructural changes of bone in mice of different ages and in models of age-related diseases such as osteoporosis. However, as aging is accompanied by frailty and subsequent increased sensitivity to external stimuli such as handling and anesthesia, the extent to which longitudinal imaging can be applied in aging studies remains unclear. Consequently, the potential of monitoring individual mice during the entire aging process-from healthy to frail status-has not yet been exploited. In this study, we assessed the effects of long-term in vivo micro-CT imaging-consisting of 11 imaging sessions over 20 weeks-on hallmarks of aging both on a local (i.e., static and dynamic bone morphometry) and systemic (i.e., frailty index (FI) and body weight) level at various stages of the aging process. Furthermore, using a premature aging model (PolgA(D257A/D257A)), we assessed whether these effects differ between genotypes. The 6th caudal vertebrae of 4 groups of mice (PolgA(D257A/D257A) and PolgA(+/+)) were monitored by in vivo micro-CT every 2 weeks. One group was subjected to 11 scans between weeks 20 and 40 of age, whereas the other groups were subjected to 5 scans between weeks 26-34, 32-40 and 40-46, respectively. The long-term monitoring approach showed small but significant changes in the static bone morphometric parameters compared to the other groups. However, no interaction effect between groups and genotype was found, suggesting that PolgA mutation does not render bone more or less susceptible to long-term micro-CT imaging. The differences between groups observed in the static morphometric parameters were less pronounced in the dynamic morphometric parameters. Moreover, the body weight and FI were not affected by more frequent imaging sessions. Finally, we observed that longitudinal designs including baseline measurements at young adult age are more powerful at detecting effects of in vivo micro-CT imaging on hallmarks of aging than cross-sectional comparisons between multiple groups of aged mice subjected to fewer imaging sessions.

Longitudinal in vivo monitoring of callus remodeling in BMP-7- and Zoledronate-treated fractures

Pharmacological interventions that combine pro-anabolic and anti-catabolic drugs to treat recalcitrant fractures have shown remarkable efficacy in augmenting the regenerative response. Specifically, in rodent models of fracture repair, treatment with BMP-7 and Zoledronate (ZA) has almost uniformally resulted in complete union. However, delayed remodeling may be problematic for ZA-treated fractures. The increase in newly formed bone is substantial but if translated in humans, delayed remodeling may delay functional recovery. Our objective was to determine if, and to what extent, bone morphogenetic protein (BMP) (in synergistically administered BMP-7 + ZA) can modulate the delayed hard callus remodeling caused by ZA. Callus remodeling in BMP-7-only and BMP-7 + ZA-treated osteotomies were monitored using in vivo µCT to follow the progression of healing at 6-week intervals over 24 weeks in an open femoral fracture rat model. None of the groups recovered baseline cortical bone volumes within 24 weeks post-osteotomy. Treatment prolonged the remodeling phase but the kinetics of remodeling appeared to differ between BMP and BMP + ZA groups. However, the mechanical characteristics were largely restored. Callus/bone volumes in BMP-only treated fractures peaked as early as week 3 suggesting that remodeling is stimulated prematurely. However, this rate of remodeling was not maintained as BMP-7 was found to exhibit negligible changes in callus/bone volumes between weeks 6 and 18, whereas declines in callus/bone volumes were present at these time points in the BMP-7 + ZA group. Our findings suggest that inclusion of ZA as an anti-catabolic agent may not be detrimental to the regenerative process despite a prolonged remodeling phase.

Hallmarks of frailty and osteosarcopenia in prematurely aged PolgA(D257A/D257A) mice

BACKGROUND

Frailty is a geriatric syndrome characterized by increased susceptibility to adverse health outcomes. One major determinant thereof is the gradual weakening of the musculoskeletal system and the associated osteosarcopenia. To improve our understanding of the underlying pathophysiology and, more importantly, to test potential interventions aimed at counteracting frailty, suitable animal models are needed.

METHODS

To evaluate the relevance of prematurely aged PolgA(D257A/D257A) mice as a model for frailty and osteosarcopenia, we quantified the clinical mouse frailty index in PolgA(D257A/D257A) and wild-type littermates (PolgA(+/+) , WT) with age and concertedly assessed the quantity and quality of bone and muscle tissue. Lastly, the anabolic responsiveness of skeletal muscle, muscle progenitors, and bone was assessed.

RESULTS

PolgA(D257A/D257A) accumulated health deficits at a higher rate compared with WT, resulting in a higher frailty index at 40 and 46 weeks of age (+166%, +278%, P < 0.0001), respectively, with no differences between genotypes at 34 weeks. Concomitantly, PolgA(D257A/D257A) displayed progressive musculoskeletal deterioration such as reduced bone and muscle mass as well as impaired functionality thereof. In addition to lower muscle weights (-14%, P < 0.05, -23%, P < 0.0001) and fibre area (-20%, P < 0.05, -22%, P < 0.0001) at 40 and 46 weeks, respectively, PolgA(D257A/D257A) showed impairments in grip strength and concentric muscle forces (P < 0.05). PolgA(D257A/D257A) mutation altered the acute response to various anabolic stimuli in skeletal muscle and muscle progenitors. While PolgA(D257A/D257A) muscles were hypersensitive to eccentric contractions as well as leucine administration, shown by larger downstream signalling response of the mechanistic target of rapamycin complex 1, myogenic progenitors cultured in vitro showed severe anabolic resistance to leucine and robust impairments in cell proliferation. Longitudinal micro-computed tomography analysis of the sixth caudal vertebrae showed that PolgA(D257A/D257A) had lower bone morphometric parameters (e.g. bone volume fraction, trabecular, and cortical thickness, P < 0.05) as well as reduced remodelling activities (e.g. bone formation and resorption rate, P < 0.05) compared with WT. When subjected to 4 weeks of cyclic loading, young but not aged PolgA(D257A/D257A) caudal vertebrae showed load-induced bone adaptation, suggesting reduced mechanosensitivity with age.

CONCLUSIONS

PolgA(D257A/D257A) mutation leads to hallmarks of age-related frailty and osteosarcopenia and provides a powerful model to better understand the relationship between frailty and the aging musculoskeletal system.

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.

Effects of Long-Term Sclerostin Deficiency on Trabecular Bone Mass and Adaption to Limb Loading Differ in Male and Female Mice

A new therapeutic option to treat osteoporosis is focused on Wnt signaling and its inhibitor sclerostin, a product of the Sost gene. In this work, we study the effect of sclerostin deficiency on trabecular bone formation and resorption in male and female mice and whether it affects mechano-responsiveness. Male and female 10- and 26-week-old Sost knockout (KO) and littermate controls (LCs) were subjected to in vivo mechanical loading of the left tibia for 2 weeks. The right tibia served as internal control. The mice were imaged using in vivo micro-computed tomography at days 0, 5, 10, and 15 and tibiae were collected for histomorphometric analyses after euthanasia. Histomorphometry and micro-CT-based 3D time-lapse morphometry revealed an anabolic and anti-catabolic effect of Sost deficiency although increased trabecular bone resorption accompanied by diminished trabecular bone formation occurred with age. Loading led to diminished resorption in adult female but not in male mice. A net gain in bone volume could be achieved with mechanical loading in Sost KO adult female mice, which occurred through a further reduction in resorbed bone volume. Our data show that sclerostin deficiency has a particularly positive effect in adult female mice. Sclerostin antibodies are approved to treat postmenopausal women with high risk of osteoporotic fractures. Further studies are required to clarify whether both sexes benefit equally from sclerostin inhibition.

Decoding rejuvenating effects of mechanical loading on skeletal aging using in vivo μCT imaging and deep learning

Throughout the process of aging, dynamic changes of bone material, micro- and macro-architecture result in a loss of strength and therefore in an increased likelihood of fragility fractures. To date, precise contributions of age-related changes in bone (re)modeling and (de)mineralization dynamics to this fragility increase are not completely understood. Here, we present an image-based deep learning approach to quantitatively describe the effects of short-term aging and adaptive response to cyclic loading applied to proximal mouse tibiae and fibulae. Our approach allowed us to perform an end-to-end age prediction based on μCT imaging to determine the dynamic biological process of aging during a two week period, therefore permitting short-term bone aging analysis with 95% accuracy in predicting time points. In a second application, our deep learning analysis reveals that two weeks of in vivo mechanical loading are associated with an underlying rejuvenating effect of 5 days. Additionally, by quantitatively analyzing the learning process, we could, for the first time, identify the localization of the age-relevant encoded information and demonstrate 89% load-induced similarity of these locations in the loaded tibia with younger control bones. These data therefore suggest that our method enables identifying a general prognostic phenotype of a certain skeletal age as well as a temporal and localized loading-treatment effect on this apparent skeletal age for the studied mouse tibia and fibula. Future translational applications of this method may provide an improved decision-support method for osteoporosis treatment at relatively low cost. STATEMENT OF SIGNIFICANCE: Bone is a highly complex and dynamic structure that undergoes changes during the course of aging as well as in response to external stimuli, such as loading. Automatic assessment of "age" and "state" of the bone may lead to early prognosis of deceases such as osteoporosis and enables evaluating the effects of certain treatments. Here, we present an artificial intelligence-based method capable of automatically predicting the skeletal age from μCT images with 95% accuracy. Additionally, we utilize it to demonstrate the rejuvenation effects of in-vivo loading treatment on bones. We further, for the first time, break down aging-related local changes in bone by quantitatively analyzing "what the age assessment model has learned" and use this information to investigate the structural details of rejuvenation process.

Lasting organ-level bone mechanoadaptation is unrelated to local strain

Bones adapt to mechanical forces according to strict principles predicting straight shape. Most bones are, however, paradoxically curved. To solve this paradox, we used computed tomography-based, four-dimensional imaging methods and computational analysis to monitor acute and chronic whole-bone shape adaptation and remodeling in vivo. We first confirmed that some acute load-induced structural changes are reversible, adhere to the linear strain magnitude regulation of remodeling activities, and are restricted to bone regions in which marked antiresorptive actions are evident. We make the novel observation that loading exerts significant lasting modifications in tibial shape and mass across extensive bone regions, underpinned by (re)modeling independent of local strain magnitude, occurring at sites where the initial response to load is principally osteogenic. This is the first report to demonstrate that bone loading stimulates nonlinear remodeling responses to strain that culminate in greater curvature adjusted for load predictability without sacrificing strength.

Bone adaptation: Safety factors and load predictability in shaping skeletal form

Much is known about skeletal adaptation in relation to the mechanical functions that bones serve. This includes how bone adapts to mechanical loading during an individual's lifetime as well as over evolutionary time. Although controlled loading in animal models allows us to observe short-term bone adaptation (epigenetic mechanobiology), examining an assemblage of extant vertebrate bones or a group of fossils' bony structures can reveal the combined effects of long-term trends in loading history and the effects of natural selection. In this survey we examine adaptations that take place over both time scales and highlight a few of the extraordinary insights first published by John Currey. First, we provide a historical perspective on bone adaptation control mechanisms, followed by a discussion of safety factors in bone. We then summarize examples of structural- and material-level adaptations and mechanotransduction, and analyze the relationship between these structural- and material-level adaptations observed in situations where loading modes are either predictable or unpredictable. We argue that load predictability is a major consideration for bone adaptation broadly across an evolutionary timescale, but that its importance can also be seen during ontogenetic growth trajectories, which are subject to natural selection as well. Furthermore, we suggest that bones with highly predictable load patterns demonstrate more precise design with lower safety factors, while bones that experience less predictable loads or those that are less capable of repair and adaptation are designed with a higher safety factor. Finally, exposure to rare loading events with high potential costs of failure leads to design of structures with very high safety factor compared to everyday loading experience. Understanding bone adaptations at the structural and material levels, which take place over an individual's lifetime or over evolutionary time has numerous applications in translational and clinical research to understand and treat musculoskeletal diseases, as well as to permit the furthering of human extraterrestrial exploration in environments with altered gravity.

Murine Axial Compression Tibial Loading Model to Study Bone Mechanobiology: Implementing the Model and Reporting Results

In vivo, tibial loading in mice is increasingly used to study bone adaptation and mechanotransduction. To achieve standardized and defined experimental conditions, loading parameters and animal-related factors must be considered when performing in vivo loading studies. In this review, we discuss these loading and animal-related experimental conditions, present methods to assess bone adaptation, and suggest reporting guidelines. This review originated from presentations by each of the authors at the workshop "Developing Best Practices for Mouse Models of In Vivo Loading" during the Preclinical Models Section at the Orthopaedic Research Society Annual Meeting, San Diego, CA, March 2017. Following the meeting, the authors engaged in detailed discussions with consideration of relevant literature. The guidelines and recommendations in this review are provided to help researchers perform in vivo loading experiments in mice, and thus further our knowledge of bone adaptation and the mechanisms involved in mechanotransduction. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 38:233-252, 2020.

Dipyridamole-loaded 3D-printed bioceramic scaffolds stimulate pediatric bone regeneration in vivo without disruption of craniofacial growth through facial maturity

This study investigates a comprehensive model of bone regeneration capacity of dypiridamole-loaded 3D-printed bioceramic (DIPY-3DPBC) scaffolds composed of 100% beta-tricalcium phosphate (β -TCP) in an immature rabbit model through the time of facial maturity. The efficacy of this construct was compared to autologous bone graft, the clinical standard of care in pediatric craniofacial reconstruction, with attention paid to volume of regenerated bone by 3D reconstruction, histologic and mechanical properties of regenerated bone, and long-term safety regarding potential craniofacial growth restriction. Additionally, long-term degradation of scaffold constructs was evaluated. At 24 weeks in vivo, DIPY-3DPBC scaffolds demonstrated volumetrically significant osteogenic regeneration of calvarial and alveolar defects comparable to autogenous bone graft with favorable biodegradation of the bioactive ceramic component in vivo. Characterization of regenerated bone reveals osteogenesis of organized, vascularized bone with histologic and mechanical characteristics comparable to native bone. Radiographic and histologic analyses were consistent with patent craniofacial sutures. Lastly, through application of 3D morphometric facial surface analysis, our results support that DIPY-3DPBC scaffolds do not cause premature closure of sutures and preserve normal craniofacial growth. Based on this novel evaluation model, this DIPY-3DPBC scaffold strategy is a promising candidate as a safe, efficacious pediatric bone tissue engineering strategy.

Remodelling of human bone on the chorioallantoic membrane of the chicken egg: De novo bone formation and resorption

Traditionally used as an angiogenic assay, the chorioallantoic membrane (CAM) assay of the chick embryo offers significant potential as an in vivo model for xenograft organ culture. Viable human bone can be cultivated on the CAM and increases in bone volume are evident; however, it remains unclear by what mechanism this change occurs and whether this reflects the physiological process of bone remodelling. In this study we tested the hypothesis that CAM-induced bone remodelling is a consequence of host and graft mediated processes. Bone cylinders harvested from femoral heads post surgery were placed on the CAM of green fluorescent protein (GFP)-chick embryos for 9 days, followed by micro computed tomography (μCT) and histological analysis. Three-dimensional registration of consecutive μCT-scans showed newly mineralised tissue in CAM-implanted bone cylinders, as well as new osteoid deposition histologically. Immunohistochemistry demonstrated the presence of bone resorption and formation markers (Cathepsin K, SOX9 and RUNX2) co-localising with GFP staining, expressed by avian cells only. To investigate the role of the human cells in the process of bone formation, decellularised bone cylinders were implanted on the CAM and comparable increases in bone volume were observed, indicating that avian cells were responsible for the bone mineralisation process. Finally, CAM-implantation of acellular collagen sponges, containing bone morphogenetic protein 2, resulted in the deposition of extracellular matrix and tissue mineralisation. These studies indicate that the CAM can respond to osteogenic stimuli and support formation or resorption of implanted human bone, providing a humanised CAM model for regenerative medicine research and a novel short-term in vivo model for tissue engineering and biomaterial testing.

Reproducibility and Radiation Effect of High-Resolution In Vivo Micro Computed Tomography Imaging of the Mouse Lumbar Vertebra and Long Bone

A moderate radiation dose, in vivo µCT scanning protocol was developed and validated for long-term monitoring of multiple skeletal sites (femur, tibia, vertebra) in mice. A customized, 3D printed mouse holder was designed and utilized to minimize error associated with animal repositioning, resulting in good to excellent reproducibility in most cortical and trabecular bone microarchitecture and density parameters except for connectivity density. Repeated in vivo µCT scans of mice were performed at the right distal femur and the 4th lumbar vertebra every 3 weeks until euthanized at 9 weeks after the baseline scan. Comparing to the non-radiated counterparts, no radiation effect was found on trabecular bone volume fraction, osteoblast and osteoblast number/surface, or bone formation rate at any skeletal site. However, trabecular number, thickness, and separation, and structure model index were sensitive to ionizing radiation associated with the µCT scans, resulting in subtle but significant changes over multiple scans. Although the extent of radiation damage on most trabecular bone microarchitecture measures are comparable or far less than the age-related changes during the monitoring period, additional considerations need to be taken to minimize the confounding radiation factors when designing experiments using in vivo µCT imaging for long-term monitoring of mouse bone.

Site-Specific Load-Induced Expansion of Sca-1+Prrx1+ and Sca-1-Prrx1+ Cells in Adult Mouse Long Bone Is Attenuated With Age

Aging is associated with significant bone loss and increased fracture risk, which has been attributed to a diminished response to anabolic mechanical loading. In adults, skeletal progenitors proliferate and differentiate into bone‐forming osteoblasts in response to increasing mechanical stimuli, though the effects of aging on this response are not well‐understood. Here we show that both adult and aged mice exhibit load‐induced periosteal bone formation, though the response is significantly attenuated with age. We also show that the acute response of adult bone to loading involves expansion of Sca‐1⁺Prrx1⁺ and Sca‐1⁻Prrx1⁺ cells in the periosteum. On the endosteal surface, loading enhances proliferation of both these cell populations, though the response is delayed by 2 days relative to the periosteal surface. In contrast to the periosteum and endosteum, the marrow does not exhibit increased proliferation of Sca‐1⁺Prrx1⁺ cells, but only of Sca‐1⁻Prrx1⁺ cells, underscoring fundamental differences in how the stem cell niche in distinct bone envelopes respond to mechanical stimuli. Notably, the proliferative response to loading is absent in aged bone even though there are similar baseline numbers of Prrx1 + cells in the periosteum and endosteum, suggesting that the proliferative capacity of progenitors is attenuated with age, and proliferation of the Sca‐1⁺Prrx1⁺ population is critical for load‐induced periosteal bone formation. These findings provide a basis for the development of novel therapeutics targeting these cell populations to enhance osteogenesis for overcoming age‐related bone loss. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

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.

Experience in the Adaptive Immunity Impacts Bone Homeostasis, Remodeling, and Healing

Bone formation as well as bone healing capacity is known to be impaired in the elderly. Although bone formation is outpaced by bone resorption in aged individuals, we hereby present a novel path that considerably impacts bone formation and architecture: Bone formation is substantially reduced in aged individual owing to the experience of the adaptive immunity. Thus, immune-aging in addition to chronological aging is a potential risk factor, with an experienced immune system being recognized as more pro-inflammatory. The role of the aging immune system on bone homeostasis and on the bone healing cascade has so far not been considered. Within this study mice at different age and immunological experience were analyzed toward bone properties. Healing was assessed by introducing an osteotomy, immune cells were adoptively transferred to disclose the difference in biological vs. chronological aging. In vitro studies were employed to test the interaction of immune cell products (cytokines) on cells of the musculoskeletal system. In metaphyseal bone, immune-aging affects bone homeostasis by impacting bone formation capacity and thereby influencing mass and microstructure of bone trabeculae leading to an overall reduced mechanical competence as found in bone torsional testing. Furthermore, bone formation is also impacted during bone regeneration in terms of a diminished healing capacity observed in young animals who have an experienced human immune system. We show the impact of an experienced immune system compared to a naïve immune system, demonstrating the substantial differences in the healing capacity and bone homeostasis due to the immune composition. We further showed that in vivo mechanical stimulation changed the immune system phenotype in young mice toward a more naïve composition. While this rescue was found to be significant in young individuals, aged mice only showed a trend toward the reconstitution of a more naïve immune phenotype. Considering the immune system's experience level in an individual, will likely allow one to differentiate (stratify) and treat (immune-modulate) patients more effectively. This work illustrates the relevance of including immune diagnostics when discussing immunomodulatory therapeutic strategies for the progressively aging population of the industrial countries.

Mechanically stimulated ATP release from murine bone cells is regulated by a balance of injury and repair

Bone cells sense and actively adapt to physical perturbations to prevent critical damage. ATP release is among the earliest cellular responses to mechanical stimulation. Mechanical stimulation of a single murine osteoblast led to the release of 70 ± 24 amole ATP, which stimulated calcium responses in neighboring cells. Osteoblasts contained ATP-rich vesicles that were released upon mechanical stimulation. Surprisingly, interventions that promoted vesicular release reduced ATP release, while inhibitors of vesicular release potentiated ATP release. Searching for an alternative ATP release route, we found that mechanical stresses induced reversible cell membrane injury in vitro and in vivo. Ca²⁺/PLC/PKC-dependent vesicular exocytosis facilitated membrane repair, thereby minimizing cell injury and reducing ATP release. Priming cellular repair machinery prior to mechanical stimulation reduced subsequent membrane injury and ATP release, linking cellular mechanosensitivity to prior mechanical exposure. Thus, our findings position ATP release as an integrated readout of membrane injury and repair.

Athletes' skeletons get stronger with training, while bones weaken in people who cannot move or in astronauts experiencing weightlessness. This is because bone cells thrive when exposed to forces. When a bone cell is exposed to a physical force, the first thing that happens is the release of the energy-rich molecule called ATP into the space outside the cell. This molecule then binds to the neighboring cell to unleash a cascade of responses. ATP can exit the cell either through special canals in the cell membrane or released in tiny pod-like structures called vesicles. It is known that strong forces can injure the cell membrane and cause ATP to spill out. However, it is less clear how ATP is released when cells are subjected to regular forces.

Mikolajewicz et al. investigated whether ATP exits through injured membranes of cells experiencing regular forces. Bone cells grown in the laboratory were gently poked with a glass needle or placed in a turbulent fluid to simulate forces experienced in the body. Dyes and fluorescent imaging techniques were used to observe the movement of vesicles and calculate the concentration of ATP in these cells.

The experiments showed that regular forces in the body do indeed injure the cell membranes and cause ATP to spill out. But importantly, the cells repaired the injuries quickly by releasing vesicles that patch the wound. As soon as the membrane is sealed, ATP stops coming out. From the first injury, cells adapted and quickly strengthened their membrane and repair system to be more resilient against future forces. This process was also seen in the shin bones of mice.

These results are important because knowing how bone cells sense, respond and convert physical forces can help us develop treatments for astronauts, the injured and aged.

Variability in strain distribution in the mice tibia loading model: A preliminary study using digital volume correlation

It is well known that bone has an enormous adaptive capacity to mechanical loadings, and to this extent, several in vivo studies on mouse tibia use established cyclic compressive loading protocols to investigate the effects of mechanical stimuli. In these experiments, the applied axial load is well controlled but the positioning of the hind-limb between the loading endcaps may dramatically affect the strain distribution induced on the tibia. In this study, the full field strain distribution induced by a typical in vivo setup on mouse tibiae was investigated through a combination of in situ compressive testing, µCT scanning and a global digital volume correlation (DVC) approach. The precision of the DVC method and the effect of repositioning on the strain distributions were evaluated. Acceptable uncertainties of the DVC approach for the analysis of loaded tibiae (411 ± 58µɛ) were found for nodal spacing of approximately 50 voxels (520 µm). When pairs of in situ preloaded and loaded images were registered, low variability of the strain distributions within the tibia were seen (range of mean differences in principal strains: 585-1800µɛ). On contrary, larger differences were seen after repositioning (range of mean differences in principal strains: 2500-5500µɛ). To conclude, these preliminary results on thee specimens showed that the DVC approach applied to the mouse tibia can be precise enough to evaluate local strain distributions under loads, and that repositioning of the hind-limb within the testing machine can induce large differences in the strain distributions that should be accounted for when modelling this system.

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.

Bone healing response in cyclically loaded implants: Comparing zero, one, and two loading sessions per day

When bone implants are loaded, they are inevitably subjected to displacement relative to bone. Such micro-motion generates stress/strain states at the interface that can cause beneficial or detrimental sequels. The objective of this study is to better understand the mechanobiology of bone healing at the tissue-implant interface during repeated loading. Machined screw shaped Ti implants were placed in rat tibiae in a hole slightly bigger than the implant diameter. Implants were held stable by a specially-designed bone plate that permits controlled loading. Three loading regimens were applied, (a) zero loading, (b) one daily loading session of 60 cycles with an axial force of 1.5 N/cycle for 7 days, and (c) two such daily sessions with the same axial force also for 7 days. Finite element analysis was used to characterize the mechanobiological conditions produced by the loading sessions. After 7 days, the implants with surrounding interfacial tissue were harvested and processed for histological, histomorphometric and DNA microarray analyses. Histomorphometric analyses revealed that the group subjected to repeated loading sessions exhibited a significant decrease in bone-implant contact and increase in bone-implant distance, as compared to unloaded implants and those subjected to only one loading session. Gene expression profiles differed during osseointegration between all groups mainly with respect to inflammatory and unidentified gene categories. The results indicate that increasing the daily cyclic loading of implants induces deleterious changes in the bone healing response, most likely due to the accumulation of tissue damage and associated inflammatory reaction at the bone-implant interface.

Transient peak-strain matching partially recovers the age-impaired mechanoadaptive cortical bone response

Mechanoadaptation maintains bone mass and architecture; its failure underlies age-related decline in bone strength. It is unclear whether this is due to failure of osteocytes to sense strain, osteoblasts to form bone or insufficient mechanical stimulus. Mechanoadaptation can be restored to aged bone by surgical neurectomy, suggesting that changes in loading history can rescue mechanoadaptation. We use non-biased, whole-bone tibial analyses, along with characterisation of surface strains and ensuing mechanoadaptive responses in mice at a range of ages, to explore whether sufficient load magnitude can activate mechanoadaptation in aged bone. We find that younger mice adapt when imposed strains are lower than in mature and aged bone. Intriguingly, imposition of short-term, high magnitude loading effectively primes cortical but not trabecular bone of aged mice to respond. This response was regionally-matched to highest strains measured by digital image correlation and to osteocytic mechanoactivation. These data indicate that aged bone’s loading response can be partially recovered, non-invasively by transient, focal high strain regions. Our results indicate that old murine bone does respond to load when the loading is of sufficient magnitude, and bones’ age-related adaptation failure may be due to insufficient mechanical stimulus to trigger mechanoadaptation.

Spatial relationship between bone formation and mechanical stimulus within cortical bone: Combining 3D fluorochrome mapping and poroelastic finite element modelling

Bone is a dynamic tissue and adapts its architecture in response to biological and mechanical factors. Here we investigate how cortical bone formation is spatially controlled by the local mechanical environment in the murine tibia axial loading model (C57BL/6). We obtained 3D locations of new bone formation by performing ‘slice and view’ 3D fluorochrome mapping of the entire bone and compared these sites with the regions of high fluid velocity or strain energy density estimated using a finite element model, validated with ex-vivo bone surface strain map acquired ex-vivo using digital image correlation. For the comparison, 2D maps of the average bone formation and peak mechanical stimulus on the tibial endosteal and periosteal surface across the entire cortical surface were created. Results showed that bone formed on the periosteal and endosteal surface in regions of high fluid flow. Peak strain energy density predicted only the formation of bone periosteally. Understanding how the mechanical stimuli spatially relates with regions of cortical bone formation in response to loading will eventually guide loading regime therapies to maintain or restore bone mass in specific sites in skeletal pathologies.

•

3D spatial representation of new bone formation after loading is shown by fluorochrome mapping of the entire mouse tibia

•

Regions of new bone formation spatially associate with regions of high strain and fluid mechanical stimulus in a FE model

•

The FE model was validated with the strains on the bone surface determined ex-vivo using digital image correlation

•

Regions of new bone formation co-localize in sites of peak fluid flow, both endosteally and periosteally

•

Peak strain energy density was able to predict only periosteal bone formation

Bone remodeling and mechanobiology around implants: Insights from small animal imaging

Anchorage of orthopedic implants depends on the interfacial bonding between the implant and the host bone as well as on the mass and microstructure of peri-implant bone, with all these factors being continuously regulated by the biological process of bone (re)modeling. In osteoporotic bone, implant integration may be jeopardized not only by lower peri-implant bone quality but also by reduced intrinsic regeneration ability. The first aim of this review is to provide a critical overview of the influence of osteoporosis on bone regeneration post-implantation. Mechanical stimulation can trigger bone formation and inhibit bone resorption; thus, judicious administration of mechanical loading can be used as an effective non-pharmacological treatment to enhance implant anchorage. Our second aim is to report recent achievements on the application of external mechanical stimulation to improve the quantity of peri-implant bone. The review focuses on peri-implant bone changes in osteoporotic conditions and following mechanical loading, prevalently using small animals and in vivo monitoring approaches. We intend to demonstrate the necessity to reveal new biological information on peri-implant bone mechanobiology to better target implant anchorage and fracture fixation in osteoporotic conditions. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:584-593, 2018.

Osteotomy models - the current status on pain scoring and management in small rodents

Fracture healing is a complex regeneration process which produces new bone tissue without scar formation. However, fracture healing disorders occur in approximately 10% of human patients and cause severe pain and reduced quality of life. Recently, the development of more standardized, sophisticated and commercially available osteosynthesis techniques reflecting clinical approaches has increased the use of small rodents such as rats and mice in bone healing research dramatically. Nevertheless, there is no standard for pain assessment, especially in these species, and consequently limited information regarding the welfare aspects of osteotomy models. Moreover, the selection of analgesics is restricted for osteotomy models since non-steroidal anti-inflammatory drugs (NSAIDs) are known to affect the initial, inflammatory phase of bone healing. Therefore, opioids such as buprenorphine and tramadol are often used. However, dosage data in the literature are varied. Within this review, we clarify the background of osteotomy models, explain the current status and challenges of animal welfare assessment, and provide an example score sheet including model specific parameters. Furthermore, we summarize current refinement options and present a brief outlook on further 3R research.

Impaired bone formation in ovariectomized mice reduces implant integration as indicated by longitudinal in vivo micro-computed tomography

Although osteoporotic bone, with low bone mass and deteriorated bone architecture, provides a less favorable mechanical environment than healthy bone for implant fixation, there is no general agreement on the impact of osteoporosis on peri-implant bone (re)modeling, which is ultimately responsible for the long term stability of the bone-implant system. Here, we inserted an implant in a mouse model mimicking estrogen deficiency-induced bone loss and we monitored with longitudinal in vivo micro-computed tomography the spatio-temporal changes in bone (re)modeling and architecture, considering the separate contributions of trabecular, endocortical and periosteal surfaces. Specifically, 12 week-old C57BL/6J mice underwent OVX/SHM surgery; 9 weeks after we inserted special metal-ceramics implants into the 6^th caudal vertebra and we measured bone response with in vivo micro-CT weekly for the following 6 weeks. Our results indicated that ovariectomized mice showed a reduced ability to increase the thickness of the cortical shell close to the implant because of impaired peri-implant bone formation, especially at the periosteal surface. Moreover, we observed that healthy mice had a significantly higher loss of trabecular bone far from the implant than estrogen depleted animals. Such behavior suggests that, in healthy mice, the substantial increase in peri-implant bone formation which rapidly thickened the cortex to secure the implant may raise bone resorption elsewhere and, specifically, in the trabecular network of the same bone but far from the implant. Considering the already deteriorated bone structure of estrogen depleted mice, further bone loss seemed to be hindered. The obtained knowledge on the dynamic response of diseased bone following implant insertion should provide useful guidelines to develop advanced treatments for osteoporotic fracture fixation based on local and selective manipulation of bone turnover in the peri-implant region.

Sost deficiency led to a greater cortical bone formation response to mechanical loading and altered gene expression

Bone adaptation optimizes mass and structure, but the mechano-response is already reduced at maturation. Downregulation of sclerostin was believed to be a mandatory step in mechano-adaptation, but in young mice it was shown that load-induced formation can occur independent of sclerostin, a product of the Sost gene. We hypothesized that the bone formation and resorption response to loading is not affected by Sost deficiency, but is age-specific. Our findings indicate that the anabolic response to in vivo tibial loading was reduced at maturation in Sost Knockout (KO) and littermate control (LC) mice. Age affected all anabolic and catabolic parameters and altered Sost and Wnt target gene expression. While load-induced cortical resorption was similar between genotypes, loading-induced gains in mineralizing surface was enhanced in Sost KO compared to LC mice. Loading led to a downregulation in expression of the Wnt inhibitor Dkk1. Expression of Dkk1 was greater in both control and loaded limbs of Sost KO compared to LC mice suggesting a compensatory role in the absence of Sost. These data suggest physical activity could enhance bone mass concurrently with sclerostin-neutralizing antibodies, but treatment strategies should consider the influence of age on ultimate load-induced bone mass gains.