Tomography-Based Quantification of Regional Differences in Cortical Bone Surface Remodeling and Mechano-Response

Aging decreases osteocyte peri-lacunar-canalicular system turnover in female C57BL/6JN mice

Osteocytes engage in bone resorption and mineralization surrounding their expansive lacunar-canalicular system (LCS) through peri-LCS turnover. However, fundamental questions persist about where, when, and how often osteocytes engage in peri-LCS turnover and how these processes change with aging. Furthermore, whether peri-LCS turnover is associated with natural variation in cortical tissue strain remains unexplored. To address these questions, we utilized confocal scanning microscopy, immunohistochemistry, and scanning electron microscopy to characterize osteocyte peri-LCS turnover in the cortical (mid-diaphysis) and cancellous (metaphysis) regions of femurs from young adult (5 mo) and early-old-age (22 mo) female C57BL/6JN mice. LCS bone mineralization was measured by the presence of perilacunar fluorochrome labels. LCS bone resorption was measured by immunohistochemical marker of bone resorption. The dynamics of peri-LCS turnover were estimated from serial fluorochrome labeling, where each mouse was administered two labels between 2 and 16 days before euthanasia. Osteocyte participation in mineralizing their surroundings is highly abundant in both cortical and cancellous bone of young adult mice but significantly decreases with aging. LCS bone resorption also decreases with aging. Aging has a greater impact on peri-LCS turnover dynamics in cancellous bone than in cortical bone. Lacunae with recent peri-LCS turnover are larger in both age groups. While peri-LCS turnover is associated with variation in tissue strain between cortical quadrants and intracortical location for 22 mo mice, these associations were not seen for 5 mo mice. The impact of aging on decreasing peri-LCS turnover may have significant implications for bone quality and mechanosensation.

Developing long bones respond to surrounding tissues by trans-pairing of periosteal osteoclasts and endocortical osteoblasts

Developing long bones alter their shape while maintaining uniform cortical thickness via coordinated activity of bone-forming osteoblasts and bone-resorbing osteoclasts at periosteal and endosteal surfaces, a process we designate trans-pairing. Two types of trans-pairing shift cortical bone in opposite orientations: peri-forming trans-pairing (peri-t-p) increases bone marrow space and endo-forming trans-pairing (endo-t-p) decreases it, via paired activity of bone resorption and formation across the cortex. Here, we focused on endo-t-p in growing bones. Analysis of endo-t-p activity in the cortex of mouse fibulae revealed osteoclasts under the periosteum compressed by muscles, and expression of RANKL in periosteal cells of the cambium layer. Furthermore, mature osteoblasts were localized on the endosteum, while preosteoblasts were at the periosteum and within cortical canals. X-ray tomographic microscopy revealed the presence of cortical canals more closely associated with endo- than with peri-t-p. Sciatic nerve transection followed by muscle atrophy and unloading induced circumferential endo-t-p with concomitant spread of cortical canals. Such canals likely supply the endosteum with preosteoblasts from the periosteum under endo-t-p, allowing bone shape to change in response to mechanical stress or nerve injury.

Mechanoregulation analysis of bone formation in tissue engineered constructs requires a volumetric method using time-lapsed micro-computed tomography

Bone can adapt its microstructure to mechanical loads through mechanoregulation of the (re)modeling process. This process has been investigated in vivo using time-lapsed micro-computed tomography (micro-CT) and micro-finite element (FE) analysis using surface-based methods, which are highly influenced by surface curvature. Consequently, when trying to investigate mechanoregulation in tissue engineered bone constructs, their concave surfaces make the detection of mechanoregulation impossible when using surface-based methods. In this study, we aimed at developing and applying a volumetric method to non-invasively quantify mechanoregulation of bone formation in tissue engineered bone constructs using micro-CT images and FE analysis. We first investigated hydroxyapatite scaffolds seeded with human mesenchymal stem cells that were incubated over 8 weeks with one mechanically loaded and one control group. Higher mechanoregulation of bone formation was measured in loaded samples with an area under the curve for the receiver operating curve (AUCformation) of 0.633-0.637 compared to non-loaded controls (AUCformation: 0.592-0.604) during culture in osteogenic medium (p < 0.05). Furthermore, we applied the method to an in vivo mouse study investigating the effect of loading frequencies on bone adaptation. The volumetric method detected differences in mechanoregulation of bone formation between loading conditions (p < 0.05). Mechanoregulation in bone formation was more pronounced (AUCformation: 0.609-0.642) compared to the surface-based method (AUCformation: 0.565-0.569, p < 0.05). Our results show that mechanoregulation of formation in bone tissue engineered constructs takes place and its extent can be quantified with a volumetric mechanoregulation method using time-lapsed micro-CT and FE analysis. STATEMENT OF SIGNIFICANCE: Many efforts have been directed towards optimizing bone scaffolds for tissue growth. However, the impact of the scaffolds mechanical environment on bone growth is still poorly understood, requiring accurate assessment of its mechanoregulation. Existing surface-based methods were unable to detect mechanoregulation in tissue engineered constructs, due to predominantly concave surfaces in scaffolds. We present a volumetric approach to enable the precise and non-invasive quantification and analysis of mechanoregulation in bone tissue engineered constructs by leveraging time-lapsed micro-CT imaging, image registration, and finite element analysis. The implications of this research extend to diverse experimental setups, encompassing culture conditions, and material optimization, and investigations into bone diseases, enabling a significant stride towards comprehensive advancements in bone tissue engineering and regenerative medicine.

Large scale analysis of osteocyte lacunae in klotho hypomorphic mice using high-resolution micro-computed tomography

BACKGROUND

Osteocytes are the most abundant cell type in adult bone, and the morphological characteristics of osteocytes and their lacunae appear to influence bone mass and fragility. Although conventional computed tomography (CT) has contributed greatly to advances in bone morphometry, capturing details of the entire hierarchical assembly, e.g., osteocyte lacuna parameters, has been limited by the analytical performance of CT (> 1 µm resolution).

METHODS

We used high-resolution (700 nm) micro-CT to evaluate and compare the osteocyte lacuna parameters over a large scale, i.e., in a maximum of about 45,700 lacunae (average), in tibial metaphyseal cortical bones of wild-type (WT) and αKlotho-hypomorphic (kl/kl) mice, the latter a model that exhibits osteopenia and aberrant osteocytes.

RESULTS

Of osteocyte lacuna parameters, lacunar surface per lacunar volume were significantly lower and lacuna diameter were significantly larger in kl/kl mice compared to WT mice. By analysis of individual osteocyte lacunae, we found that lacunar sphericity in kl/kl mice was higher than that in WT mice, and the diameters in the major and the minor axes were respectively lower and higher in kl/kl mice, especially at the proximal site of the region of interest.

CONCLUSION

We successfully assessed osteocyte lacuna parameters on the largest scale in mice reported to date and found that the shape of osteocyte lacunae of kl/kl mice are significantly different from those of WT mice. Although the mechanisms underlying the lacunar shape differences observed are not yet clear, changes in lacunar geometry are known to affect the transitions of strains to the osteocyte microenvironment and likely local osteocyte response(s). Thus, the fact that the differences are limited to the mesial region near the primary spongiosa suggests the likelihood of site-specific anomalies in mechanosensitive effects in kl/kl osteocytes with consequent site-specific effects bone metabolism and function.

Mechanostat parameters estimated from time-lapsed in vivo micro-computed tomography data of mechanically driven bone adaptation are logarithmically dependent on loading frequency

Mechanical loading is a key factor governing bone adaptation. Both preclinical and clinical studies have demonstrated its effects on bone tissue, which were also notably predicted in the mechanostat theory. Indeed, existing methods to quantify bone mechanoregulation have successfully associated the frequency of (re)modeling events with local mechanical signals, combining time-lapsed in vivo micro-computed tomography (micro-CT) imaging and micro-finite element (micro-FE) analysis. However, a correlation between the local surface velocity of (re)modeling events and mechanical signals has not been shown. As many degenerative bone diseases have also been linked to impaired bone (re)modeling, this relationship could provide an advantage in detecting the effects of such conditions and advance our understanding of the underlying mechanisms. Therefore, in this study, we introduce a novel method to estimate (re)modeling velocity curves from time-lapsed in vivo mouse caudal vertebrae data under static and cyclic mechanical loading. These curves can be fitted with piecewise linear functions as proposed in the mechanostat theory. Accordingly, new (re)modeling parameters can be derived from such data, including formation saturation levels, resorption velocity moduli, and (re)modeling thresholds. Our results revealed that the norm of the gradient of strain energy density yielded the highest accuracy in quantifying mechanoregulation data using micro-finite element analysis with homogeneous material properties, while effective strain was the best predictor for micro-finite element analysis with heterogeneous material properties. Furthermore, (re)modeling velocity curves could be accurately described with piecewise linear and hyperbola functions (root mean square error below 0.2 µm/day for weekly analysis), and several (re)modeling parameters determined from these curves followed a logarithmic relationship with loading frequency. Crucially, (re)modeling velocity curves and derived parameters could detect differences in mechanically driven bone adaptation, which complemented previous results showing a logarithmic relationship between loading frequency and net change in bone volume fraction over 4 weeks. Together, we expect this data to support the calibration of in silico models of bone adaptation and the characterization of the effects of mechanical loading and pharmaceutical treatment interventions in vivo.

Molecular Identification of Spatially Distinct Anabolic Responses to Mechanical Loading in Murine Cortical Bone

Osteoporosis affects over 200 million women worldwide, one-third of whom are predicted to suffer from an osteoporotic fracture in their lifetime. The most promising anabolic drugs involve administration of expensive antibodies. Because mechanical loading stimulates bone formation, our current data, using a mouse model, replicates the anabolic effects of loading in humans and may identify novel pathways amenable to oral treatment. Murine tibial compression produces axially varying deformations along the cortical bone, inducing highest strains at the mid-diaphysis and lowest at the metaphyseal shell. To test the hypothesis that load-induced transcriptomic responses at different axial locations of cortical bone would vary as a function of strain magnitude, we loaded the left tibias of 10-week-old female C57Bl/6 mice in vivo in compression, with contralateral limbs as controls. Animals were euthanized at 1, 3, or 24 hours post-loading or loaded for 1 week (n = 4-5/group). Bone marrow and cancellous bone were removed, cortical bone was segmented into the metaphyseal shell, proximal diaphysis, and mid-diaphysis, and load-induced differential gene expression and enriched biological processes were examined for the three segments. At each time point, the mid-diaphysis (highest strain) had the greatest transcriptomic response. Similarly, biological processes regulating bone formation and turnover increased earlier and to the greatest extent at the mid-diaphysis. Higher strain induced greater levels of osteoblast and osteocyte genes, whereas expression was lower in osteoclasts. Among the top differentially expressed genes at 24-hours post-loading, 17 had known functions in bone biology, of which 12 were present only in osteoblasts, 3 exclusively in osteoclasts, and 2 were present in both cell types. Based on these results, we conclude that murine tibial loading induces spatially unique transcriptomic responses correlating with strain magnitude in cortical bone. © 2022 American Society for Bone and Mineral Research (ASBMR).

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.

Mechanoregulated trabecular bone adaptation: Progress report on in silico approaches

Adaptation is the process by which bone responds to changes in loading environment and modulates its properties and spatial organization to meet the mechanical demands. Adaptation in trabecular bone is achieved through increase in bone mass and alignment of trabecular-bone morphology along the loading direction. This transformation of internal microstructure is governed by mechanical stimuli sensed by mechanosensory cells in the bone matrix. Realisation of adaptation in the form of local bone-resorption and -formation activities as a function of mechanical stimuli is still debated. In silico modelling is a useful tool for simulation of various scenarios that cannot be investigated in vivo and particularly well suited for prediction of trabecular bone adaptation. This progress report presents the recent advances in in silico modelling of mechanoregulated adaptation at the scale of trabecular bone tissue. Four well-established bone-adaptation models are reviewed in terms of their recent improvements and validation. They consider various mechanical factors: (i) strain energy density, (ii) strain and damage, (iii) stress nonuniformity and (iv) daily stress. Contradictions of these models are discussed and their ability to describe adequately a real-life mechanoregulation process in bone is compared.

Osteocytes regulate bone anabolic response to mechanical loading in male mice via activation of integrin α5

Physical mechanical stimulation can maintain and even increase bone mass. Here, we report an important role of osteocytic integrin α5 in regulating the anabolic response of bone to mechanical loading using an Itga5 conditional gene knockout (cKO) mouse model. Integrin α5 gene deletion increased apoptotic osteocytes and reduced cortical anabolic responses to tibial compression including decreased endosteal osteoblasts and bone formation, and increased endosteal osteoclasts and bone resorption, contributing to the decreased bone area fraction and biomechanical properties, leading to an enlarged bone marrow area in cKO mice. Similar disruption of anabolic responses to mechanical loading was also detected in cKO trabecular bone. Moreover, integrin α5 deficiency impeded load-induced Cx43 hemichannel opening, and production and release of PGE2, an anabolic factor, resulting in attenuated effects of the loading on catabolic sclerostin (SOST) reduction and anabolic β-catenin increase. Together, this study shows an indispensable role of integrin α5 in osteocytes in the anabolic action of mechanical loading on skeletal tissue through activation of hemichannels and PGE2-evoked gene expression. Integrin α5 could act as a potential new therapeutic target for bone loss, especially in the elderly population with impeded mechanical sensitivity.

An in silico model for woven bone adaptation to heavy loading conditions in murine tibia

Existing in silico models for lamellar bone adaptation to mechanical loading are unsuitable for predicting woven bone growth. This anomaly is due to the difference in mechanobiology of the woven bone with respect to that of the lamellar bone. The present study is aimed at developing an in silico bone-adaptation model for woven bone at cellular and tissue levels. The diffusion of Ca²⁺ ions reaching lining cells from the osteocytic network and the bone cortex in response to a mechanical loading on the cortical bone has been considered as a stimulus. The diffusion of ions within osteocytic network has been computed with a lacunar-canalicular network (LCN) in which bone cells are uniformly arranged. Strain energy density is assumed to regulate ion flow within the network when the induced normal strain is above a threshold level. If the induced strain exceeds another higher threshold level, then the strain with a power constant is additionally assumed to regulate the stimulus. The intracellular flow of Ca²⁺ ions within the LCN has been simulated using Fick's laws of diffusion, using a finite element method. The ion diffusion from bone cortex to vesicles has been formulated as a normal strain with a power constant. The stimuli reaching the surface cells are assumed to form the new bone. The mathematical model closely predicts woven bone growth in mouse and rat tibia for various in vivo loading conditions. This model is the first to predict woven bone growth at tissue and cellular levels in response to heavy mechanical loading.

Connexin hemichannels with prostaglandin release in anabolic function of bone to mechanical loading

Mechanical stimulation, such as physical exercise, is essential for bone formation and health. Here, we demonstrate the critical role of osteocytic Cx43 hemichannels in anabolic function of bone in response to mechanical loading. Two transgenic mouse models, R76W and Δ130–136, expressing dominant-negative Cx43 mutants in osteocytes were adopted. Mechanical loading of tibial bone increased cortical bone mass and mechanical properties in wild-type and gap junction-impaired R76W mice through increased PGE₂, endosteal osteoblast activity, and decreased sclerostin. These anabolic responses were impeded in gap junction/hemichannel-impaired Δ130–136 mice and accompanied by increased endosteal osteoclast activity. Specific inhibition of Cx43 hemichannels by Cx43(M1) antibody suppressed PGE₂ secretion and impeded loading-induced endosteal osteoblast activity, bone formation and anabolic gene expression. PGE₂ administration rescued the osteogenic response to mechanical loading impeded by impaired hemichannels. Together, osteocytic Cx43 hemichannels could be a potential new therapeutic target for treating bone loss and osteoporosis.

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.

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.

Bursa-Derived Cells Show a Distinct Mechano-Response to Physiological and Pathological Loading in vitro

The mechano-response of highly loaded tissues such as bones or tendons is well investigated, but knowledge regarding the mechano-responsiveness of adjacent tissues such as the subacromial bursa is missing. For a better understanding of the physiological role of the bursa as a friction-reducing structure in the joint, the study aimed to analyze whether and how bursa-derived cells respond to physiological and pathological mechanical loading. This might help to overcome some of the controversies in the field regarding the role of the bursa in the development and healing of shoulder pathologies. Cells of six donors seeded on collagen-coated silicon dishes were stimulated over 3 days for 1 or 4 h with 1, 5, or 10% strain. Orientation of the actin cytoskeleton, YAP nuclear translocation, and activation of non-muscle myosin II (NMM-II) were evaluated for 4 h stimulations to get a deeper insight into mechano-transduction processes. To investigate the potential of bursa-derived cells to adapt their matrix formation and remodeling according to mechanical loading, outcome measures included cell viability, gene expression of extracellular matrix and remodeling markers, and protein secretions. The orientation angle of the actin cytoskeleton increased toward a more perpendicular direction with increased loading and lowest variations for the 5% loading group. With 10% tension load, cells were visibly stressed, indicated by loss in actin density and slightly reduced cell viability. A significantly increased YAP nuclear translocation occurred for the 1% loading group with a similar trend for the 5% group. NMM-II activation was weak for all stimulation conditions. On the gene expression level, only the expression of TIMP2 was down-regulated in the 1 h group compared to control. On the protein level, collagen type I and MMP2 increased with higher/longer straining, respectively, whereas TIMP1 secretion was reduced, resulting in an MMP/TIMP imbalance. In conclusion, this study documents for the first time a clear mechano-responsiveness in bursa-derived cells with activation of mechano-transduction pathways and thus hint to a physiological function of mechanical loading in bursa-derived cells. This study represents the basis for further investigations, which might lead to improved treatment options of subacromial bursa-related pathologies in the future.

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.

Quantifying Asymmetry in Gait: The Weighted Universal Symmetry Index to Evaluate 3D Ground Reaction Forces

Though gait asymmetry is used as a metric of functional recovery in clinical rehabilitation, there is no consensus on an ideal method for its evaluation. Various methods have been proposed to analyze single bilateral signals but are limited in scope, as they can often use only positive signals or discrete values extracted from time-scale data as input. By defining five symmetry axioms, a framework for benchmarking existing methods was established and a new method was described here for the first time: the weighted universal symmetry index (wUSI), which overcomes limitations of other methods. Both existing methods and the wUSI were mathematically compared to each other and in respect to their ability to fulfill the proposed symmetry axioms. Eligible methods that fulfilled these axioms were then applied using both discrete and continuous approaches to ground reaction force (GRF) data collected from healthy gait, both with and without artificially induced asymmetry using a single instrumented elbow crutch. The wUSI with a continuous approach was the only symmetry method capable of identifying GRF asymmetry differences in different walking conditions in all three planes of motion. When used with a continuous approach, the wUSI method was able to detect asymmetries while avoiding artificial inflation, a common problem reported in other methods. In conclusion, the wUSI is proposed as a universal method to quantify three-dimensional GRF asymmetries, which may also be expanded to other biomechanical signals.

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.

Mechano-Regulation of Trabecular Bone Adaptation Is Controlled by the Local in vivo Environment and Logarithmically Dependent on Loading Frequency

It is well-established that cyclic, but not static, mechanical loading has anabolic effects on bone. However, the function describing the relationship between the loading frequency and the amount of bone adaptation remains unclear. Using a combined experimental and computational approach, this study aimed to investigate whether trabecular bone mechano-regulation is controlled by mechanical signals in the local in vivo environment and dependent on loading frequency. Specifically, by combining in vivo micro-computed tomography (micro-CT) imaging with micro-finite element (micro-FE) analysis, we monitored the changes in microstructural as well as the mechanical in vivo environment [strain energy density (SED) and SED gradient] of mouse caudal vertebrae over 4 weeks of either cyclic loading at varying frequencies of 2, 5, or 10 Hz, respectively, or static loading. Higher values of SED and SED gradient on the local tissue level led to an increased probability of trabecular bone formation and a decreased probability of trabecular bone resorption. In all loading groups, the SED gradient was superior in the determination of local bone formation and resorption events as compared to SED. Cyclic loading induced positive net (re)modeling rates when compared to sham and static loading, mainly due to an increase in mineralizing surface and a decrease in eroded surface. Consequently, bone volume fraction increased over time in 2, 5, and 10 Hz (+15%, +21% and +24%, p ≤ 0.0001), while static loading led to a decrease in bone volume fraction (-9%, p ≤ 0.001). Furthermore, regression analysis revealed a logarithmic relationship between loading frequency and the net change in bone volume fraction over the 4 week observation period (R 2 = 0.74). In conclusion, these results suggest that trabecular bone adaptation is regulated by mechanical signals in the local in vivo environment and furthermore, that mechano-regulation is logarithmically dependent on loading frequency with frequencies below a certain threshold having catabolic effects, and those above anabolic effects. This study thereby provides valuable insights toward a better understanding of the mechanical signals influencing trabecular bone formation and resorption in the local in vivo environment.

Anatomical variations in cortical bone surface permeability: Tibia versus femur

Cortical bone surfaces (periosteal and endosteal) exhibit differential (re)modelling response to mechanical loading. This poses a serious challenge in establishing an in silico model to predict site-specific new bone formation as a function of mechanical stimulus. In this regard, mechanical loading-induced fluid motion in lacunar-canalicular system (LCS) is assumed osteogenic. Micro-architectural properties, especially permeability regulate canalicular fluid motion within the bone. The knowledge of these properties is required to compute flow distribution. Along the same line, it is possible that cortical surfaces may experience differential fluid distribution due to anatomical variations in microarchitectural properties which may induce distinct new bone response at cortical surfaces. Nevertheless, these properties are not well reported for cortical surfaces in the literature. Accordingly, the present study aims to measure microarchitectural properties especially permeability at different anatomical locations (medial, lateral, anterior, and posterior) of periosteal and endosteal surfaces using nanoindentation. A standard poroelastic optimization technique was used to estimate permeability, shear modulus, and Poisson's ratio. The properties are also compared for two weight-bearing bones i.e. tibia and femur. Endosteal surface was found more permeable as compared to the periosteal surface. Tibial endosteal surface had shown greater permeability values at most of the anatomical locations as compared to femoral endosteal surface. The outcomes may be used to precisely predict site-specific osteogenesis in cortical bone as a function of canalicular flow distribution. This work may ultimately be beneficial in designing the loading parameters to stimulate desired new bone response for the prevention and the cure of bone loss.

Mechano-adaptive Responses of Alveolar Bone to Implant Hyper-loading in a pre-clinical in vivo model

OBJECTIVES

Oral implants transmit biting forces to peri-implant bone. In turn, those forces subject peri-implant bone to mechanical stresses and strains. Here, our objective was to understand how peri-implant bone responded to conditions of normal versus hyper-loading in a mouse model.

MATERIAL AND METHODS

Sixty-six mice were randomly assigned to 2 groups; both groups underwent bilateral maxillary first molar extraction followed by complete healing. Titanium alloy implants were placed in healed sites and positioned below the occlusal plane. After osseointegration, a composite crown was affixed to the implant so masticatory loading would ensue. In controls, the remaining dentition was left intact but in the hyper-loaded (test) group, the remaining molars were extracted. 3D finite element analysis (FEA) calculated peri-implant strains resulting from normal and hyper-loading. Peri-implant tissues were analyzed at multiple time points using micro-computed tomography (µCT) imaging, histology, enzymatic assays of bone remodeling, and vital dye labeling to evaluate bone accrual.

RESULTS

Compared to controls, hyper-loaded implants experienced a 3.6-fold increase in occlusal force, producing higher peri-implant strains. Bone formation and resorption were both significantly elevated around hyper-loaded implants, eventually culminating in a significant increase in peri-implant bone volume/total volume (BV/TV). In our mouse model, masticatory hyper-loading of an osseointegrated implant was associated with increased peri-implant strain, increased peri-implant bone remodeling, and a net gain in bone deposition.

CONCLUSION

Hyper-loading results in bone strain with catabolic and anabolic bone responses, leading to a net gain in bone deposition.

Interspecies Comparison of Alveolar Bone Biology, Part I: Morphology and Physiology of Pristine Bone

INTRODUCTION

Few interspecies comparisons of alveolar bone have been documented, and this knowledge gap raises questions about which animal models most accurately represent human dental conditions or responses to surgical interventions.

OBJECTIVES

The objective of this study was to employ state-of-the-art quantitative metrics to directly assess and compare the structural and functional characteristics of alveolar bone among humans, mini pigs, rats, and mice.

METHODS

The same anatomic location (i.e., the posterior maxillae) was analyzed in all species via micro-computed tomographic imaging, followed by quantitative analyses, coupled with histology and immunohistochemistry. Bone remodeling was evaluated with alkaline phosphatase activity and tartrate-resistant acid phosphatase staining to identify osteoblast and osteoclast activities. In vivo fluorochrome labeling was used as a means to assess mineral apposition rates.

RESULTS

Collectively, these analyses demonstrated that bone volume differed among the species, while bone mineral density was equal. All species showed a similar density of alveolar osteocytes, with a highly conserved pattern of collagen organization. Collagen maturation was equal among mouse, rat, and mini pig. Bone remodeling was a shared feature among the species, with morphologically indistinguishable hemiosteonal appearances, osteocytic perilacunar remodeling, and similar mineral apposition rates in alveolar bone.

CONCLUSIONS

Our analyses demonstrated equivalencies among the 4 species in a plurality of the biological features of alveolar bone. Despite contradictory results from older studies, we found no evidence for the superiority of pig models over rodent models in representing human bone biology.

KNOWLEDGE TRANSFER STATEMENT

Animal models are extensively used to evaluate bone tissue engineering strategies, yet there are few state-of-the-art studies that rigorously compare and quantify the factors influencing selection of a given animal model. Consequently, there is an urgent need to assess preclinical animal models for their predictive value to dental research. Our article addresses this knowledge gap and, in doing so, provides a foundation for more effective standardization among animal models commonly used in dentistry.

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.

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.

In Vivo and In Vitro Mechanical Loading of Mouse Achilles Tendons and Tenocytes-A Pilot Study

Mechanical force is a key factor for the maintenance, adaptation, and function of tendons. Investigating the impact of mechanical loading in tenocytes and tendons might provide important information on in vivo tendon mechanobiology. Therefore, the study aimed at understanding if an in vitro loading set up of tenocytes leads to similar regulations of cell shape and gene expression, as loading of the Achilles tendon in an in vivo mouse model. In vivo: The left tibiae of mice (n = 12) were subject to axial cyclic compressive loading for 3 weeks, and the Achilles tendons were harvested. The right tibiae served as the internal non-loaded control. In vitro: tenocytes were isolated from mice Achilles tendons and were loaded for 4 h or 5 days (n = 6 per group) based on the in vivo protocol. Histology showed significant differences in the cell shape between in vivo and in vitro loading. On the molecular level, quantitative real-time PCR revealed significant differences in the gene expression of collagen type I and III and of the matrix metalloproteinases (MMP). Tendon-associated markers showed a similar expression profile. This study showed that the gene expression of tendon markers was similar, whereas significant changes in the expression of extracellular matrix (ECM) related genes were detected between in vivo and in vitro loading. This first pilot study is important for understanding to which extent in vitro stimulation set-ups of tenocytes can mimic in vivo characteristics.

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.

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.

Sost deficiency leads to reduced mechanical strains at the tibia midshaft in strain-matched in vivo loading experiments in mice

Sclerostin, a product of the Sost gene, is a Wnt-inhibitor and thus negatively regulates bone accrual. Canonical Wnt/β-catenin signalling is also known to be activated in mechanotransduction. Sclerostin neutralizing antibodies are being tested in ongoing clinical trials to target osteoporosis and osteogenesis imperfecta but their interaction with mechanical stimuli on bone formation remains unclear. Sost knockout (KO) mice were examined to gain insight into how long-term Sost deficiency alters the local mechanical environment within the bone. This knowledge is crucial as the strain environment regulates bone adaptation. We characterized the bone geometry at the tibial midshaft of young and adult Sost KO and age-matched littermate control (LC) mice using microcomputed tomography imaging. The cortical area and the minimal and maximal moment of inertia were higher in Sost KO than in LC mice, whereas no difference was detected in either the anterior-posterior or medio-lateral bone curvature. Differences observed between age-matched genotypes were greater in adult mice. We analysed the local mechanical environment in the bone using finite-element models (FEMs), which showed that strains in the tibiae of Sost KO mice are lower than in age-matched LC mice at the diaphyseal midshaft, a region commonly used to assess cortical bone formation and resorption. Our FEMs also suggested that tissue mineral density is only a minor contributor to the strain distribution in tibial cortical bone from Sost KO mice compared to bone geometry. Furthermore, they indicated that although strain gauging experiments matched strains at the gauge site, strains along the tibial length were not comparable between age-matched Sost KO and LC mice or between young and adult animals within the same genotype.

Adaptive changes in micromechanical environments of cancellous and cortical bone in response to in vivo loading and disuse

The skeleton accommodates changes in mechanical environments by increasing bone mass under increased loads and decreasing bone mass under disuse. However, little is known about the adaptive changes in micromechanical behavior of cancellous and cortical tissues resulting from loading or disuse. To address this issue, in vivo tibial loading and hindlimb unloading experiments were conducted on 16-week-old female C57BL/6J mice. Changes in bone mass and tissue-level strains in the metaphyseal cancellous and midshaft cortical bone of the tibiae, resulting from loading or unloading, were determined using microCT and finite element (FE) analysis, respectively. We found that loading- and unloading-induced changes in bone mass were more pronounced in the cancellous than cortical bone. Simulated FE-loading showed that a greater proportion of elements experienced relatively lower longitudinal strains following load-induced bone adaptation, while the opposite was true in the disuse model. While the magnitudes of maximum or minimum principal strains in the metaphyseal cancellous and midshaft cortical bone were not affected by loading, strains oriented with the long axis were reduced in the load-adapted tibia suggesting that loading-induced micromechanical benefits were aligned primarily in the loading direction. Regression analyses demonstrated that bone mass was a good predictor of bone tissue strains for the cortical bone but not for the cancellous bone, which has complex microarchitecture and spatially-variant strain environments. In summary, loading-induced micromechanical benefits for cancellous and cortical tissues are received primarily in the direction of force application and cancellous bone mass may not be related to the micromechanics of cancellous bone.

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

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.