Experimental and finite element analysis of the rat ulnar loading model-correlations between strain and bone formation following fatigue loading

Computer model of non-Newtonian canalicular fluid flow in lacunar-canalicular system of bone tissue

Brittle bone diseases are a global healthcare problem for orthopaedic clinicians, that reduces bone strength and promotes bone fracture risk. In vivo studies reported that loading-induced fluid flow through the lacunar-canalicular channel (LCS) of bone tissue inhibit such bone loss and encourages osteogenesis i.e. new bone formation. Canalicular fluid flow converts mechanical signals into biological signals and regulates bone reconstruction by releasing signalling molecules responsible for mechanotransduction. In-silico model mostly considers canalicular fluid is Newtonian, however, physiological canalicular fluid may be non-Newtonian in nature as it contains nutrients and supplements. Accordingly, this study attempts to develop a two-dimensional in-silico model to compute loading-induced non-Newtonian canalicular fluid flow in a complex LCS of bone tissue. Moreover, canalicular fluid is considered as a Jeffery fluid, that can easily be reduced to Newtonian fluid as a special case. The results show that physiological loading modulates the canalicular fluid flow, wall shear stress (WSS) and streamline in bone LCS. Fluid velocity and WSS increases with increase in non-dimensional frequency and non-Newtonian parameter (Jeffery fluid parameters) and reduce with change in permeability. The outcomes of this study may provide new insights in the role of mechanical loading-induced non-Newtonian canalicular fluid flow dynamics in bone LCS. The key findings of this study can be used to improve the understanding of osteocyte mechanobiology involved inside the bone tissue.

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.

Optogenetic-induced muscle loading leads to mechanical adaptation of the Achilles tendon enthesis in mice

Skeletal shape depends on the transmission of contractile muscle forces from tendon to bone across the enthesis. Loss of muscle loading impairs enthesis development, yet little is known if and how the postnatal enthesis adapts to increased loading. Here, we studied adaptations in enthesis structure and function in response to increased loading, using optogenetically induced muscle contraction in young (i.e., growth) and adult (i.e., mature) mice. Daily bouts of unilateral optogenetic loading in young mice led to radial calcaneal expansion and warping. This also led to a weaker enthesis with increased collagen damage in young tendon and enthisis, with little change in adult mice. We then used RNA sequencing to identify the pathways associated with increased mechanical loading during growth. In tendon, we found enrichment of glycolysis, focal adhesion, and cell-matrix interactions. In bone, we found enrichment of inflammation and cell cycle. Together, we demonstrate the utility of optogenetic-induced muscle contraction to elicit in vivo adaptation of the enthesis.

Optogenetic-Induced Muscle Loading Leads to Mechanical Adaptation of the Achilles Tendon Enthesis in Mice

The growth of the skeleton depends on the transmission of contractile muscle forces from tendon to bone across the extracellular matrix-rich enthesis. Loss of muscle loading leads to significant impairments in enthesis development. However, little is known about how the enthesis responds to increased loading during postnatal growth. To study the cellular and matrix adaptations of the enthesis in response to increased muscle loading, we used optogenetics to induce skeletal muscle contraction and unilaterally load the Achilles tendon and enthesis in young (i.e., during growth) and adult (i.e., mature) mice. In young mice, daily bouts of unilateral optogenetic loading led to expansion of the calcaneal apophysis and growth plate, as well as increased vascularization of the normally avascular enthesis. Daily loading bouts, delivered for 3 weeks, also led to a mechanically weaker enthesis with increased molecular-level accumulation of collagen damage in young mice. However, adult mice did not exhibit impaired mechanical properties or noticeable structural adaptations to the enthesis. We then focused on the transcriptional response of the young tendon and bone following optogenetic-induced loading. After 1 or 2 weeks of loading, we identified, in tendon, transcriptional activation of canonical pathways related to glucose metabolism (glycolysis) and inhibited pathways associated with cytoskeletal remodeling (e.g., RHOA and CREB signaling). In bone, we identified activation of inflammatory signaling (e.g., NFkB and STAT3 signaling) and inhibition of ERK/MAPK and PTEN signaling. Thus, we have demonstrated the utility of optogenetic-induced skeletal muscle contraction to elicit structural, functional, and molecular adaptation of the enthesis in vivo especially during growth.

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

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

Animal models of impaired long bone healing and tissue engineering- and cell-based in vivo interventions

Bone healing after injury typically follows a systematic process and occurs spontaneously under appropriate physiological conditions. However, impaired long bone healing is still quite common and may require surgical intervention. Various complications can result in different forms of impaired bone healing including nonunion, critical-size defects, or stress fractures. While a nonunion may occur due to impaired biological signaling and/or mechanical instability, a critical-size defect exhibits extensive bone loss that will not spontaneously heal. Comparatively, a stress fracture occurs from repetitive forces and results in a non-healing crack or break in the bone. Clinical standards of treatment vary between these bone defects due to their pathological differences. The use of appropriate animal models for modeling healing defects is critical to improve current treatment methods and develop novel rescue therapies. This review provides an overview of these clinical bone healing impairments and current animal models available to study the defects in vivo. The techniques used to create these models are compared, along with the outcomes, to clarify limitations and future objectives. Finally, rescue techniques focused on tissue engineering and cell-based therapies currently applied in animal models are specifically discussed to analyze their ability to initiate healing at the defect site, providing information regarding potential future therapies. In summary, this review focuses on the current animal models of nonunion, critical-size defects, and stress fractures, as well as interventions that have been tested in vivo to provide an overview of the clinical potential and future directions for improving bone healing.

Flow inside a bone scaffold: Visualization using 3D phase contrast MRI and comparison with numerical simulations

We report on results of experimental flow measurements inside a bone scaffold model, subjected to a uniform incoming flow (applied perfusion). Understanding the flow behavior inside a tissue engineered scaffold is essential for mechanistic studies of mechanobiology, particularly flow-sensitive bone cells. Nearly all existing studies that quantify interstitial flow inside engineered bone scaffolds have been based on numerical results, in part due to the difficulties associated with quantitative measurements and visualization of flow inside large, opaque bone or bone mimics. Thus, an experimental platform to complement and validate in silico studies is needed. Therefore, we developed a flow visualization method using Phase-Contrast Magnetic Resonance Imaging (PC-MRI) to measure flow velocities within a 3D-printed microCT-based rendering of a bone scaffold. We designed and built a non-magnetic recirculating water tunnel to apply uniform perfusion to the 3D-printed model and we measured flow distribution within the scaffold and compared these experimental results with CFD results. Both magnitude and distribution of flow velocities observed at different slices of the scaffold were in quantitative agreement numerically and experimentally. This experimental approach can be used to both validate numerical studies and provide insight into the flow behavior inside tissue-engineered scaffolds for a range of applications, including fundamental mechanobiology of healthy cells, and in the context of diseases, such as cancer.

Implant resonance and the mechanostat theory: Applications of therapeutic ultrasound for porous metallic scaffolds

The development of treatment strategies for improving secondary stability at the bone-implant interface is a challenge. Porous implants are one solution for improving long-term implant stability, but the osteoconduction process of implants into the bone can be slow. Strain-driven osteogenesis from the mechanostat theory offers insight into pathways for post-operative treatment but mechanisms to deliver strain to the bone-implant interface need refinement. In this work, the use of therapeutic ultrasound is simulated to induce resonance into a porous implant structure. Local strains through the scaffold are measured by varying systemic variables such as damping ratio, applied vibrational force, primary bone-implant stability, and input frequency. At the natural frequency of the system with applied forces of 0.5 N and a damping ratio of 0.5%, roughly half of the nodes in the simulated environment exceed the microstrain threshold of 1000 με required for new bone formation. A high degree of sensitivity was noted upon changing input frequency, with minor sensitivities arising from damping ratio and applied vibrational force. These findings suggest that the application of therapeutic resonance to improve osseointegration of the bone-implant interface may be viable for applications including dental implants or segmental bone defects.

Relating Bone Strain to Local Changes in Radius Microstructure Following 12 Months of Axial Forearm Loading in Women

Work in animal models suggests that bone structure adapts to local bone strain, but this relationship has not been comprehensively studied in humans. Here, we quantified the influence of strain magnitude and gradient on bone adaptation in the forearm of premenopausal women performing compressive forearm loading (n = 11) and nonloading controls (n = 10). High resolution peripheral quantitative computed tomography (HRpQCT) scans of the distal radius acquired at baseline and 12 months of a randomized controlled experiment were used to identify local sites of bone formation and resorption. Bone strain was estimated using validated finite element (FE) models. Trabecular strain magnitude and gradient were higher near (within 200 μm) formation versus resorption (p < 0.05). Trabecular formation and resorption occurred preferentially near very high (>95th percentile) versus low (<5th percentile) strain magnitude and gradient elements, and very low strain elements were more likely to be near resorption than formation (p < 0.05). In the cortical compartment, strain gradient was higher near formation versus resorption (p < 0.05), and both formation and resorption occurred preferentially near very high versus low strain gradient elements (p < 0.05). At most, 54% of very high and low strain elements were near formation or resorption only, and similar trends were observed in the control and load groups. These findings suggest that strain, likely in combination with other physiological factors, influences adaptation under normal loads and in response to a novel loading intervention, and represents an important step toward defining exercise interventions to maximize bone strength.

Phenotype characterization of crossed beaks in Beijing-You chickens based on morphological observation

The prevalence of crossed beaks ranging from 0.2 to 7.4% was documented in at least 12 chicken strains. Previous studies focused largely on candidate molecules, whereas the morphological observation was missing. This study reported a detailed phenotype and prevalence of crossed beaks based on morphological observation in nine thousand nine hundred 1-day-old female Beijing-You chicks. Affected chicks were classified into 2 categories based on the direction of the mandibular deformation: left and right. Each category was selected to sacrifice for the measurement of length, width, and thickness of the bilateral mandibular ramus (MR). The normal chicks were used as controls. Paraffin section was made for the bilateral MR of a crossed beak and a normal control for histology analysis. A total of 97 out of 9,900 chickens showed beak deformity including 71 crossed beaks (0.72%) and 26 side beaks (0.26%) for which the upper and lower beak were both bent in the same direction. There was no difference in the direction of the bend of the lower beak in crossed beaks (P > 0.05). The incidence of crossed beaks increased quickly from 0 to 56 d and no new incidence after 56 d. The angle of the crossed beaks was below 5° in the first week and had grown more severe with age until 56 d. The mandible structure showed that condyle served as a growth center for the MR extension. The short-side MR of crossed beaks was thicker than normal ones (P < 0.05) and caused the mandible deviated to the same direction. Meanwhile, the short-side MR prevented the occlusion, leading the jugal arch deformity, which in turn resulted in a bent maxillary horizontally. Similarly, chicks with side beaks also had asymmetry in MR length and the deformities of the jugal arch after dissection. In summary, asymmetric growth of bilateral MR induced crossed beaks and side beaks; the mandibular condyle could be an ideal sample for the related molecular mechanism studies underlying this trait.

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.

Strain Distribution Evaluation of Rat Tibia under Axial Compressive Load by Combining Strain Gauge Measurement and Finite Element Analysis

This study is aimed at providing an effective method for determining strain-load relationship and at quantifying the strain distribution within the whole tibia under axial compressive load on rats. Rat tibial models with axial compressive load were designed. Strains in three directions (0°, 45°, and 90°) at the proximal shaft of the tibia were measured by using a strain gauge rosette, which was used to calculate the maximum and minimum principal strains. Moreover, the strain at the midshaft of the tibia was measured by a single-element strain gauge. The slopes of the strain-load curves with different peak loads were calculated to assess the stability of the strain gauge measurement. Mechanical environment in the whole tibia by the axial compressive load was quantified using finite element analysis (FEA) based on microcomputed tomography images. The von Mises elastic strain distributions of the whole tibiae were evaluated. Slopes of the strain-load curves showed no significant differences among different peak loads (ANOVA; P > 0.05), indicating that the strain-load relationship obtained from the strain gauge measurement was reasonable and stable. The FEA results corresponded to the experimental results with an error smaller than 15% (paired Student's t-test, P > 0.05), signifying that the FEA can simulate the experiment reasonably. FEA results showed that the von Mises elastic strain was the lowest in the middle and gradually increased to both sides along the lateral direction, with the maximal von Mises elastic strain being observed on the posterior side under the distal tibiofibular synostosis. The method of strain gauge measurements and FEA used in this study can provide a feasible way to obtain the mechanical environment of the tibiae under axial compressive load on the rats and serve as a reference for further exploring the mechanical response of the bone by axial compressive load.

Prolonged high force high repetition pulling induces osteocyte apoptosis and trabecular bone loss in distal radius, while low force high repetition pulling induces bone anabolism

We have an operant rat model of upper extremity reaching and grasping in which we examined the impact of performing a high force high repetition (High-ForceHR) versus a low force low repetition (Low-ForceHR) task for 18weeks on the radius and ulna, compared to age-matched controls. High-ForceHR rats performed at 4 reaches/min and 50% of their maximum voluntary pulling force for 2h/day, 3days/week. Low-ForceHR rats performed at 6% maximum voluntary pulling force. High-ForceHR rats showed decreased trabecular bone volume in the distal metaphyseal radius, decreased anabolic indices in this same bone region (e.g., decreased osteoblasts and bone formation rate), and increased catabolic indices (e.g., microcracks, increased osteocyte apoptosis, secreted sclerostin, RANKL, and osteoclast numbers), compared to controls. Distal metaphyseal trabeculae in the ulna of High-ForceHR rats showed a non-significant decrease in bone volume, some catabolic indices (e.g., decreased trabecular numbers) yet also some anabolic indices (e.g., increased osteoblasts and trabecular thickness). In contrast, the mid-diaphyseal region of High-ForceHR rats' radial and ulnar bones showed few to no microarchitecture differences and no changes in apoptosis, sclerostin or RANKL levels, compared to controls. In further contrast, Low-ForceHR rats showed increased trabecular bone volume in the radius in the distal metaphysis and increased cortical bone area its mid-diaphysis. These changes were accompanied by increased anabolic indices, no microcracks or osteocyte apoptosis, and decreased RANKL in each region, compared to controls. Ulnar bones of Low-ForceHR rats also showed increased anabolic indices, although fewer than in the adjacent radius. Thus, prolonged performance of an upper extremity reaching and grasping task is loading-, region-, and bone-dependent, with high force loads at high repetition rates inducing region-specific increases in bone degradative changes that were most prominent in distal radial trabeculae, while low force task loads at high repetition rates induced adaptive bone responses.

Spatiotemporal characterization of microdamage accumulation in rat ulnae in response to uniaxial compressive fatigue loading

Repetitive fatigue loading can induce microdamage accumulation in bone matrix, which results in impaired mechanical properties and increased fracture susceptibility. However, the spatial distribution and time-variant process of microdamage accumulation in fatigue-loaded skeleton, especially for linear microcracks which are known to initiate bone remodeling, remain not fully understood. In this study, the time-varying process of the morphology and distribution of microcracks in rat ulnae subjected to uniaxial compressive fatigue loading was investigated. Right forelimbs of thirty four-month-old male Sprague-Dawley rats were subjected to one bout of cyclic ramp loading with 0.67 Hz at a normalized peak force of 0.055 N/g body weight for 6000 cycles, and the contralateral left ulnae were not loaded as the control samples. Ten rats were randomly euthanized on Days 3, 5, and 7 post fatigue loading. Our findings via two-dimensional histomorphometric measurements based on basic fuchsin staining and three-dimensional quantifications using contrast-enhanced micro-computed tomography (MicroCT) with precipitated BaSO₄ staining demonstrated that the accumulation of linear microcracks (increase in the amount of linear microcracks) on Day 5 was significantly higher than that on Day 3 and Day 7 post fatigue loading. Our histological and histomorphometric results revealed that linear microcrack density (Cr.Dn) in the tensile cortex at Days 3, 5 and 7 post fatigue loading was significantly higher than that in the compressive side, whereas linear microcrack length (Cr.Le) in the tensile cortex at Day 3 was significantly lower than that in the compressive cortex. Our findings revealed that microcrack accumulation exhibited a non-linear time-varying process at 3, 5 and 7 days post axial compressive fatigue loading (with observable peak Cr.Dn at Day 5). Our findings also revealed distinct distribution of microcrack density and morphology in rat ulnae with tensile and compressive strains, as characterized by more microcracks accumulated in tensile cortices, and longer cracks shown in compressive cortices.

Temporal mechanically-induced signaling events in bone and dorsal root ganglion neurons after in vivo bone loading

Mechanical signals play an integral role in the regulation of bone mass and functional adaptation to bone loading. The osteocyte has long been considered the principle mechanosensory cell type in bone, although recent evidence suggests the sensory nervous system may play a role in mechanosensing. The specific signaling pathways responsible for functional adaptation of the skeleton through modeling and remodeling are not clearly defined. In vitro studies suggest involvement of intracellular signaling through mitogen-activated protein kinase (MAPK), phosphatidylinositol 3-kinase (PI3K)/protein kinase B (Akt), and mammalian target of rapamycin (mTOR). However, anabolic signaling responses to bone loading using a whole animal in vivo model have not been studied in detail. Therefore, we examined mechanically-induced signaling events at five time points from 0 to 24 hours after loading using the rat in vivo ulna end-loading model. Western blot analysis of bone for MAPK's, PI3K/Akt, and mTOR signaling, and quantitative reverse transcription polymerase chain reaction (qRT-PCR) to estimate gene expression of calcitonin gene-related protein alpha (CGRP-α), brain-derived neurotrophic factor (BDNF), nerve growth factor (NGF), c-jun, and c-fos in dorsal root ganglion (DRG) of the brachial intumescence were performed. There was a significant increase in signaling through MAPK's including extracellular signal-related kinase (ERK) and c-Jun N-terminal kinase (JNK) in loaded limbs at 15 minutes after mechanical loading. Ulna loading did not significantly influence expression of the genes of interest in DRG neurons. Bone signaling and DRG gene expression from the loaded and contralateral limbs was correlated (SR>0.40, P<0.05). However, bone signaling did not correlate with expression of the genes of interest in DRG neurons. These results suggest that signaling through the MAPK pathway may be involved in load-induced bone formation in vivo. Further characterization of the molecular events involved in regulation of bone adaptation is needed to understand the timing and impact of loading events, and the contribution of the neuronal signaling to functional adaptation of bone.

The fracture toughness of small animal cortical bone measured using arc-shaped tension specimens: Effects of bisphosphonate and deproteinization treatments

Small animal models, and especially transgenic models, have become widespread in the study of bone mechanobiology and metabolic bone disease, but test methods for measuring fracture toughness on multiple replicates or at multiple locations within a single small animal bone are lacking. Therefore, the objective of this study was to develop a method to measure cortical bone fracture toughness in multiple specimens and locations along the diaphysis of small animal bones. Arc-shaped tension specimens were prepared from the mid-diaphysis of rabbit ulnae and loaded to failure to measure the radial fracture toughness in multiple replicates per bone. The test specimen dimensions, crack length, and maximum load met requirements for measuring the plane strain fracture toughness. Experimental groups included a control group, bisphosphonate treatment group, and an ex vivo deproteinization treatment following bisphosphonate treatment (5 rabbits/group and 15 specimens/group). The fracture toughness of ulnar cortical bone from rabbits treated with zoledronic acid for six months exhibited no difference compared with the control group. Partially deproteinized specimens exhibited significantly lower fracture toughness compared with both the control and bisphosphonate treatment groups. The deproteinization treatment increased tissue mineral density (TMD) and resulted in a negative linear correlation between the measured fracture toughness and TMD. Fracture toughness measurements were repeatable with a coefficient of variation of 12-16% within experimental groups. Retrospective power analysis of the control and deproteinization treatment groups indicated a minimum detectable difference of 0.1MPa·m^1/2. Therefore, the overall results of this study suggest that arc-shaped tension specimens offer an advantageous new method for measuring the fracture toughness in small animal bones.

Effect of athletic fatigue damage and the associated bone targeted remodeling in the rat ulna

Background

Fatigue damage of the long bones is prevalent in running athletes and military recruits due to vigorous mid- and long-term physical activity. The current study attempted to know the features of bony athletic fatigue damage and to explore the mechanism of fatigue damage repair through bone targeted remodeling process.

Methods

Right ulnae of the Wistar rats were fatigue loaded on an INSTRON 5865 to construct the athletic fatigue damage model, and several time points (i.e. experimental days: 0, 7, 13 and 19) were selected to simulate physiological status, preliminary, mid-term and perennial stage during continuous high-intensive training, respectively. The multi-level responses of rat ulnae under the athletic fatigue loading, including cellular protein expression, micro damage or micro-crack and macro mechanical properties, were tested and statistically analyzed.

Results

Wistar rats, subjected to the athletic fatigue loading protocol, experienced a decrease of ulna fatigue mechanical properties and an active bone resorption of the loaded ulnae in the early stage, whereafter, a hyperactive bone formation and significant improvements of ulnae fatigue mechanical properties were detected. However, a deterioration of quasi-static mechanical properties in the subsequent period implied limitations of bone remodeling to maintain the bearing capacity of bone during long-term strenuous exercise.

Conclusions

In summary, after athletic fatigue loading, bone targeted remodeling is activated and proceeds to repair fatigue damage, but only to a certain extent.

Bone mechanobiology in mice: toward single-cell in vivo mechanomics

Mechanically driven bone (re)modeling is a multiscale process mediated through complex interactions between multiple cell types and their microenvironments. However, the underlying mechanisms of how cells respond to mechanical signals are still unclear and are at the focus of the field of bone mechanobiology. Traditionally, this complex process has been addressed by reducing the system to single scales and cell types. It is only recently that more integrative approaches have been established to study bone mechanobiology across multiple scales in which mechanical load at the organ level is related to molecular responses at the cellular level. The availability of mouse loading models and imaging techniques with improved spatial and temporal resolution has made it possible to track dynamic bone (re)modeling at the tissue and cellular level in vivo. Coupled with advanced computational models, the (re)modeling activities at the tissue scale can be associated with the mechanical microenvironment. However, methods are lacking to link the molecular responses of different cell types to their local mechanical microenvironment and bone (re)modeling activities occurring at the tissue scale. With recent improvements in "omics" technologies and single-cell molecular biology, it is now possible to sequence the complete genome and transcriptome of single cells. These technologies offer unique opportunities to comprehensively investigate the cellular transcriptional profiles within their specific microenvironment. By combining single-cell "omics" technologies with well-established tissue-scale models of bone mechanobiology, we propose a mechanomics approach to locally analyze the transcriptome of single cells with respect to their local 3D mechanical in vivo environment.

Development of an in vivo bone fatigue damage model using axial compression of the rabbit forelimb

Many nontraumatic fractures seen clinically in patients with metabolic bone disorders or on antiresorptive treatment show an increased incidence of microdamage accumulation and impaired intracortical remodeling. However, the lack of basal remodeling and Haversian bone in rodents limits their translatability in studying bone damage repair mechanisms. The work presented here demonstrates the development of the forelimb loading model in rabbits, the smallest mammal with intracortical Haversian remodeling. The forelimbs of post-mortem female New Zealand white rabbits were loaded in axial end compression to determine their basic monotonic and fatigue properties. Following time zero characterization, stress fractures were created in vivo and animals were allowed to recover for a period of two to five weeks. The rabbit forelimb when loaded in axial compression demonstrates a consistent mid-diaphyseal fracture location characterized by a local mixed compression-bending loading environment. Forelimb apparent stiffness, when fatigue loaded, demonstrates a progressive increase until macrocrack formation, at which time apparent stiffness rapidly declines until failure. Stress fractures in the rabbit ulna display robust periosteal expansion and woven bone formation two weeks following fracture. Subsequent healing at five weeks post-fracture is marked by woven bone densification, resorption and intracortical remodeling along the stress fracture line. The rabbit forelimb fatigue model is a promising new platform by which bone׳s response to damage may be studied.

Development of an in vivo rabbit ulnar loading model

Ulnar and tibial cyclic compression in rats and mice have become the preferred animal models for investigating the effects of mechanical loading on bone modeling/remodeling. Unlike rodents, rabbits provide a larger bone volume and normally exhibit intracortical Haversian remodeling, which may be advantageous for investigating mechanobiology and pharmaceutical interventions in cortical bone. Therefore, the objective of this study was to develop and validate an in vivo rabbit ulnar loading model. Ulnar tissue strains during loading of intact forelimbs were characterized and calibrated to applied loads using strain gauge measurements and specimen-specific finite element models. Periosteal bone formation in response to varying strain levels was measured by dynamic histomorphometry at the location of maximum strain in the ulnar diaphysis. Ulnae loaded at 3000 microstrain did not exhibit periosteal bone formation greater than the contralateral controls. Ulnae loaded at 3500, 4000, and 4500 microstrain exhibited a dose-dependent increase in periosteal mineralizing surface (MS/BS) compared with contralateral controls during the second week of loading. Ulnae loaded at 4500 microstrain exhibited the most robust response with significantly increased MS/BS at multiple time points extending at least 2weeks after loading was ceased. Ulnae loaded at 5250 microstrain exhibited significant woven bone formation. Rabbits required greater strain levels to produce lamellar and woven bone on periosteal surfaces compared with rats and mice, perhaps due to lower basal levels of MS/BS. In summary, bone adaptation during rabbit ulnar loading was tightly controlled and may provide a translatable model for human bone biology in preclinical investigations of metabolic bone disease and pharmacological treatments.

Strain energy density gradients in bone marrow predict osteoblast and osteoclast activity: a finite element study

Huiskes et al. hypothesized that mechanical strains sensed by osteocytes residing in trabecular bone dictate the magnitude of load-induced bone formation. More recently, the mechanical environment in bone marrow has also been implicated in bone׳s response to mechanical stimulation. In this study, we hypothesize that trabecular load-induced bone formation can be predicted by mechanical signals derived from an integrative µFE model, incorporating a description of both the bone and marrow phase. Using the mouse tail loading model in combination with in vivo micro-computed tomography (µCT) we tracked load induced changes in the sixth caudal vertebrae of C57BL/6 mice to quantify the amount of newly mineralized and eroded bone volumes. To identify the mechanical signals responsible for adaptation, local morphometric changes were compared to micro-finite element (µFE) models of vertebrae prior to loading. The mechanical parameters calculated were strain energy density (SED) on trabeculae at bone forming and resorbing surfaces, SED in the marrow at the boundary between bone forming and resorbing surfaces, along with SED in the trabecular bone and marrow volumes. The gradients of each parameter were also calculated. Simple regression analysis showed mean SED gradients in the trabecular bone matrix to significantly correlate with newly mineralized and eroded bone volumes R(2)=0.57 and 0.41, respectively, p<0.001). Nevertheless, SED gradients in the marrow were shown to be the best predictor of osteoblastic and osteoclastic activity (R(2)=0.83 and 0.60, respectively, p<0.001). These data suggest that the mechanical environment of the bone marrow plays a significant role in determining osteoblast and osteoclast activity.

Bone adaptation to cyclic loading in murine caudal vertebrae is maintained with age and directly correlated to the local micromechanical environment

The ability of the skeleton to adapt to mechanical stimuli (mechanosensitivity) has most often been investigated at the whole-bone level, but less is known about the local mechanoregulation of bone remodeling at the bone surface, especially in context of the aging skeleton. The aim of this study was to determine the local and global mechanosensitivity of the sixth caudal vertebra during cyclic loading (8 N, three times per week, for six weeks) in mice aged 15, 52, and 82 weeks at the start of loading. Bone adaptation was monitored with in vivo micro-computed tomography. Strain energy density (SED), assumed as the mechanical stimulus for bone adaptation, was determined with micro-finite element models. Mechanical loading had a beneficial effect on the bone microstructure and bone stiffness in all age groups. Mineralizing surface was on average 13% greater (p<0.05) in loaded than control groups in 15- and 82-week-old mice, but not for 52-week-old mice. SED at the start of loading correlated to the change in bone volume fraction in the following 6 weeks for loaded groups (r(2)=0.69-0.85) but not control groups. At the local level, SED was 14-20% greater (p<0.01) at sites of bone formation, and 15-20% lower (p<0.01) at sites of bone resorption compared to quiescent bone surfaces for all age groups, indicating SED was a stimulus for bone adaptation. Taken together, these results support that mechanosensitivity is maintained with age in caudal vertebrae of mice at a local and global level. Since age-related bone loss was not observed in caudal vertebrae, results from the current study might not be translatable to aged humans.

Structural and mechanical repair of diffuse damage in cortical bone in vivo

Physiological wear and tear causes bone microdamage at several hierarchical levels, and these have different biological consequences. Bone remodeling is widely held to be the mechanism by which bone microdamage is repaired. However, recent studies showed that unlike typical linear microcracks, small crack damage, the clusters of submicron-sized matrix cracks also known as diffuse damage (Dif.Dx), does not activate remodeling. Thus, the fate of diffuse damage in vivo is not known. To examine this, we induced selectively Dif.Dx in rat ulnae in vivo by using end-load ulnar bending creep model. Changes in damage content were assessed by histomorphometry and mechanical testing immediately after loading (ie, acute loaded) or at 14 days after damage induction (ie, survival ulnae). Dif.Dx area was markedly reduced over the 14-day survival period after loading (p < 0.02). We did not observe any intracortical resorption, and there was no increase in cortical bone area in survival ulnae. The reduction in whole bone stiffness in acute loaded ulnae was restored to baseline levels in survival ulnae (p > 0.6). Microindentation studies showed that Dif.Dx caused a highly localized reduction in elastic modulus in diffuse damage regions of the ulnar cortex. Moduli in these previously damaged bone areas were restored to control values by 14 days after loading. Our current findings indicate that small crack damage in bone can be repaired without bone remodeling, and they suggest that alternative repair mechanisms exist in bone to deal with submicron-sized matrix cracks. Those mechanisms are currently unknown and further investigations are needed to elucidate the mechanisms by which this direct repair occurs.

Experimental and finite element analysis of dynamic loading of the mouse forearm

Bone formation is reported to initiate in osteocytes by mechanotransduction due to dynamic loading of bone. The first step towards this is to characterize the dynamic strain fields in the overall bone. Here, the previously developed mouse forearm ulna-radius model, subjected to static loading, has been further enhanced by incorporating a loading cap and applying a cyclic dynamic load to more closely approximate experimental biological conditions. This study also incorporates data obtained from strain gauging both the ulna and radius simultaneously. Based on separate experiments, the elastic modulus of the ulna and radius were determined to be 13.8 and 9.9 GPa, respectively. Another novel aspect of the numerical model is the inclusion of the interosseous membrane in the FE model with membrane stiffness ranging from 5-15 N/mm that have been found to give strain values closer to that from the experiments. Interestingly, the inclusion of the interosseous membrane helped to equalize the peak strain magnitudes in the ulna and radius (∼1800 at 2 N load and ∼3200 at 3.5 N), which was also observed experimentally. This model represents a significant advance towards being able to simulate through FE analysis the strain fields generated in vivo upon mechanical loading of the mouse forearm.

Comparison of three methods of calculating strain in the mouse ulna in exogenous loading studies

Axial compression of mouse limbs is commonly used to induce bone formation in a controlled, non-invasive manner. Determination of peak strains caused by loading is central to interpreting results. Load-strain calibration is typically performed using uniaxial strain gauges attached to the diaphyseal, periosteal surface of a small number of sacrificed animals. Strain is measured as the limb is loaded to a range of physiological loads known to be anabolic to bone. The load-strain relationship determined by this subgroup is then extrapolated to a larger group of experimental mice. This method of strain calculation requires the challenging process of strain gauging very small bones which is subject to variability in placement of the strain gauge. We previously developed a method to estimate animal-specific periosteal strain during axial ulnar loading using an image-based computational approach that does not require strain gauges. The purpose of this study was to compare the relationship between load-induced bone formation rates and periosteal strain at ulnar midshaft using three different methods to estimate strain: (A) Nominal strain values based solely on load-strain calibration; (B) Strains calculated from load-strain calibration, but scaled for differences in mid-shaft cross-sectional geometry among animals; and (C) An alternative image-based computational method for calculating strains based on beam theory and animal-specific bone geometry. Our results show that the alternative method (C) provides comparable correlation between strain and bone formation rates in the mouse ulna relative to the strain gauge-dependent methods (A and B), while avoiding the need to use strain gauges.

Torsion and antero-posterior bending in the in vivo human tibia loading regimes during walking and running

Bending, in addition to compression, is recognized to be a common loading pattern in long bones in animals. However, due to the technical difficulty of measuring bone deformation in humans, our current understanding of bone loading patterns in humans is very limited. In the present study, we hypothesized that bending and torsion are important loading regimes in the human tibia. In vivo tibia segment deformation in humans was assessed during walking and running utilizing a novel optical approach. Results suggest that the proximal tibia primarily bends to the posterior (bending angle: 0.15°-1.30°) and medial aspect (bending angle: 0.38°-0.90°) and that it twists externally (torsion angle: 0.67°-1.66°) in relation to the distal tibia during the stance phase of overground walking at a speed between 2.5 and 6.1 km/h. Peak posterior bending and peak torsion occurred during the first and second half of stance phase, respectively. The peak-to-peak antero-posterior (AP) bending angles increased linearly with vertical ground reaction force and speed. Similarly, peak-to-peak torsion angles increased with the vertical free moment in four of the five test subjects and with the speed in three of the test subjects. There was no correlation between peak-to-peak medio-lateral (ML) bending angles and ground reaction force or speed. On the treadmill, peak-to-peak AP bending angles increased with walking and running speed, but peak-to-peak torsion angles and peak-to-peak ML bending angles remained constant during walking. Peak-to-peak AP bending angle during treadmill running was speed-dependent and larger than that observed during walking. In contrast, peak-to-peak tibia torsion angle was smaller during treadmill running than during walking. To conclude, bending and torsion of substantial magnitude were observed in the human tibia during walking and running. A systematic distribution of peak amplitude was found during the first and second parts of the stance phase.

The Contribution of Experimental in vivo Models to Understanding the Mechanisms of Adaptation to Mechanical Loading in Bone

Changing loading regimens by natural means such as exercise, with or without interference such as osteotomy, has provided useful information on the structure:function relationship in bone tissue. However, the greatest precision in defining those aspects of the overall strain environment that influence modeling and remodeling behavior has been achieved by relating quantified changes in bone architecture to quantified changes in bones' strain environment produced by direct, controlled artificial bone loading. Jiri Hert introduced the technique of artificial loading of bones in vivo with external devices in the 1960s using an electromechanical device to load rabbit tibiae through transfixing stainless steel pins. Quantifying natural bone strains during locomotion by attaching electrical resistance strain gages to bone surfaces was introduced by Lanyon, also in the 1960s. These studies in a variety of bones in a number of species demonstrated remarkable uniformity in the peak strains and maximum strain rates experienced. Experiments combining strain gage instrumentation with artificial loading in sheep, pigs, roosters, turkeys, rats, and mice has yielded significant insight into the control of strain-related adaptive (re)modeling. This diversity of approach has been largely superseded by non-invasive transcutaneous loading in rats and mice, which is now the model of choice for many studies. Together such studies have demonstrated that over the physiological strain range, bone's mechanically adaptive processes are responsive to dynamic but not static strains; the size and nature of the adaptive response controlling bone mass is linearly related to the peak loads encountered; the strain-related response is preferentially sensitive to high strain rates and unresponsive to static ones; is most responsive to unusual strain distributions; is maximized by remarkably few strain cycles, and that these are most effective when interrupted by short periods of rest between them.

Magnitude of loads influences the site of failure of highly curved bones

The structure and material properties of bones along with applied boundary conditions determine the region of peak stresses, where fracture is expected to occur. As the site of peak stresses is not influenced by the magnitude of applied load, the fracture site is not expected to change during fatigue loading of whole bone at different loads. However, in a highly curved bone such as the rat ulna, the magnitude of applied loads was found to influence the fracture site. Fatigue loading was conducted under load control on intact rat forearms and on excised ulnae. The distance to the site of failure from the proximal olecranon process of ulnae was determined. In intact forearms, the site of failure demonstrated a linear progression distally, towards sites with lower moment of inertia (or sites exhibiting lower section modulus). Intact rat forearms and excised ulnae loaded to failure at low loads fractured 2-3mm distal to where they failed when applying high loads. This indicates a shift in the site of failure by approximately 10% of whole bone length just by varying the applied load magnitude. The site of failure in excised ulnae was similar when loading at 2Hz or at 4Hz, suggesting that this was frequency independent in this range and indicating that strain rate was not an important contributing factor. Creep loading of excised ulnae also demonstrated similar changes in the site of failure, indicating that magnitude of loads, and not type of loading were important in determining the site of failure. This has important implications with regards to the volume of bone that undergoes damage under physiological loading, before it fails.

Experimental and finite element analysis of strains induced by axial tibial compression in young-adult and old female C57Bl/6 mice

Axial compression of the mouse tibia is used to study strain-adaptive bone (re)modeling. In some studies, comparisons between mice of different ages are of interest. We characterized the tibial deformation and force-strain relationships in female C57Bl/6 mice at 5-, 12- and 22-months age. A three-gauge experimental method was used to determine the strain distribution at the mid-diaphysis, while specimen-specific finite element analysis was used to examine strain distribution along the tibial length. The peak strains in the tibial mid-diaphyseal cross-section are compressive and occur at the postero-lateral apex. The magnitudes of these peak compressive strains are 1.5 to 2 times those on the opposite, antero-medial face (a site often used for strain gauge placement). For example, -10 N force applied to a 5-months old mouse engenders a peak compressive strain of -2800 µε and a tensile strain on the antero-medial face of +1450 µε. The orientation of the neutral axis at the mid-diaphysis did not differ with age (p=0.46), indicating a similar deformation mode in young and old tibiae. On the other hand, from 5- to 22-months there is a 25% reduction in cortical thickness and moment of inertia (p<0.05), resulting in significantly greater tibial strain magnitudes in older mice for equivalent applied force (p<0.05). We conclude that comparisons of tibial loading responses in young-adult and old C57Bl/6 tibiae are facilitated by similar deformation pattern across ages, but that modest adjustment of force levels is required to engender matching peak strains.

A three-dimensional tissue culture model of bone formation utilizing rotational co-culture of human adult osteoblasts and osteoclasts

Living bone is a complex, three-dimensional composite material consisting of numerous cell types spatially organized within a mineralized extracellular matrix. To date, mechanistic investigation of the complex cellular level cross-talk between the major bone-forming cells involved in the response of bone to mechanical and biochemical stimuli has been hindered by the lack of a suitable in vitro model that captures the "coupled" nature of this response. Using a novel rotational co-culture approach, we have generated large (>4mm diameter), three-dimensional mineralized tissue constructs from a mixture of normal human primary osteoblast and osteoclast precursor cells without the need for any exogenous osteoconductive scaffolding material that might interfere with such cell-cell interactions. Mature, differentiated bone constructs consist of an outer region inhabited by osteoclasts and osteoblasts and a central region containing osteocytes encased in a self-assembled, porous mineralized extracellular matrix. Bone constructs exhibit morphological, mineral and biochemical features similar to remodeling human trabecular bone, including the expression of mRNA for SOST, BGLAP, ACP5, BMP-2, BMP-4 and BMP-7 within the construct and the secretion of BMP-2 protein into the medium. This "coupled" model of bone formation will allow the future investigation of various stimuli on the process of normal bone formation/remodeling as it relates to the cellular function of osteoblasts, osteoclasts and osteocytes in the generation of human mineralized tissue.

Adaptive and Injury Response of Bone to Mechanical Loading

Bone responds to supraphysiological mechanical loads by increasing bone formation. Depending on the applied strain magnitude (and other loading parameters) the response can be either adaptive (mostly lamellar bone) or injury (mostly woven bone). Seminal studies of Hert, Lanyon, and Rubin originally established the basic "rules" of bone mechanosensitivity. These were reinforced by subsequent studies using non-invasive rodent loading models, most notably by Turner et al. More recent work with these models have been able to explore the structural, transcriptional, and molecular mechanisms which distinguish the two responses (lamellar vs. woven). Wnt/Lrp signaling has emerged as a key mechanoresponsive pathway for lamellar bone. However, there is still much to study with regard to effects of ageing, osteocytes, other signaling pathways, and the molecular regulation that modulates lamellar vs. woven bone formation. This review summarizes not only the historical findings but also the current data for these topics.

Functional adaptation in female rats: the role of estrogen signaling

Background

Sex steroids have direct effects on the skeleton. Estrogen acts on the skeleton via the classical genomic estrogen receptors alpha and beta (ERα and ERβ), a membrane ER, and the non-genomic G-protein coupled estrogen receptor (GPER). GPER is distributed throughout the nervous system, but little is known about its effects on bone. In male rats, adaptation to loading is neuronally regulated, but this has not been studied in females.

Methodology/Principal Findings

We used the rat ulna end-loading model to induce an adaptive modeling response in ovariectomized (OVX) female Sprague-Dawley rats. Rats were treated with a placebo, estrogen (17β-estradiol), or G-1, a GPER-specific agonist. Fourteen days after OVX, rats underwent unilateral cyclic loading of the right ulna; half of the rats in each group had brachial plexus anesthesia (BPA) of the loaded limb before loading. Ten days after loading, serum estrogen concentrations, dorsal root ganglion (DRG) gene expression of ERα, ERβ, GPER, CGRPα, TRPV1, TRPV4 and TRPA1, and load-induced skeletal responses were quantified. We hypothesized that estrogen and G-1 treatment would influence skeletal responses to cyclic loading through a neuronal mechanism. We found that estrogen suppresses periosteal bone formation in female rats. This physiological effect is not GPER-mediated. We also found that absolute mechanosensitivity in female rats was decreased, when compared with male rats. Blocking of adaptive bone formation by BPA in Placebo OVX females was reduced.

Conclusions

Estrogen acts to decrease periosteal bone formation in female rats in vivo. This effect is not GPER-mediated. Gender differences in absolute bone mechanosensitivity exist in young Sprague-Dawley rats with reduced mechanosensitivity in females, although underlying bone formation rate associated with growth likely influences this observation. In contrast to female and male rats, central neuronal signals had a diminished effect on adaptive bone formation in estrogen-deficient female rats.

Probabilistic finite element analysis of a craniofacial finite element model

We employed a probabilistic finite element analysis (FEA) method to determine how variability in material property values affects stress and strain values in a finite model of a Macaca fascicularis cranium. The material behavior of cortical bone varied in three ways: isotropic homogeneous, isotropic non-homogeneous, and orthotropic non-homogeneous. The material behavior of the trabecular bone and teeth was always treated as isotropic and homogeneous. All material property values for the cranium were randomized with a Gaussian distribution with either coefficients of variation (CVs) of 0.2 or with CVs calculated from empirical data. Latin hypercube sampling was used to determine the values of the material properties used in the finite element models. In total, four hundred and twenty six separate deterministic FE simulations were executed. We tested four hypotheses in this study: (1) uncertainty in material property values will have an insignificant effect on high stresses and a significant effect on high strains for homogeneous isotropic models; (2) the effect of variability in material property values on the stress state will increase as non-homogeneity and anisotropy increase; (3) variation in the in vivo shear strain values reported by Strait et al. (2005) and Ross et al. (2011) is not only due to variations in muscle forces and cranial morphology, but also due to variation in material property values; (4) the assumption of a uniform coefficient of variation for the material property values will result in the same trend in how moderate-to-high stresses and moderate-to-high strains vary with respect to the degree of non-homogeneity and anisotropy as the trend found when the coefficients of variation for material property values are calculated from empirical data. Our results supported the first three hypotheses and falsified the fourth. When material properties were varied with a constant CV, as non-homogeneity and anisotropy increased the level of variability in the moderate-to-high strains decreased while the level of variability in the moderate-to-high stresses increased. However, this is not the pattern observed when CVs calculated from empirical data were applied to the material properties where the lowest level of variability in both stresses and strains occurred when the cranium was modeled with a low level of non-homogeneity and anisotropy. Therefore, when constant material property variability is assumed, inaccurate trends in the level of variability present in modest-to-high magnitude stresses and strains are produced. When the cranium is modeled with the highest level of accuracy (high non-homogeneity and anisotropy) and when randomness in the material properties is calculated from empirical data, there is a large level of variability in the significant strains (CV=0.369) and a low level of variability in the modest-to-high magnitude stresses (CV=0.150). This result may have important implications with regard to the mechanical signals driving bone remodeling and adaptation through natural selection.

Load/strain distribution between ulna and radius in the mouse forearm compression loading model

Finite element analysis (FEA) of the mouse forearm compression loading model is used to relate strain distributions with downstream changes in bone formation and responses of bone cells. The objective of this study was to develop two FEA models - the first one with the traditional ulna only and the second one in which both the ulna and radius are included, in order to examine the effect of the inclusion of the radius on the strain distributions in the ulna. The entire mouse forearm was scanned using microCT and images were converted into FEA tetrahedral meshes using a suite of software programs. The performance of both linear and quadratic tetrahedral elements and coarse and fine meshes were studied. A load of 2N was applied to the ulna/radius model and a 1.3N load (based on previous investigations of load sharing between the ulna and radius in rats) was applied to the ulna only model for subsequent simulations. The results showed differences in the cross sectional strain distributions and magnitude within the ulna for the combined ulna/radius model versus the ulna only model. The maximal strain in the combined model occurred about 4mm toward the distal end from the ulna mid-shaft in both models. Results from the FEA model simulations were also compared to experimentally determined strain values. We conclude that inclusion of the radius in FE models to predict strains during in vivo forearm loading increases the magnitude of the estimated ulna strains compared to those predicted from a model of the ulna alone but the distribution was similar. This has important ramifications for future studies to understand strain thresholds needed to activate bone cell responses to mechanical loading.

Role of calcitonin gene-related peptide in bone repair after cyclic fatigue loading

Background

Calcitonin gene related peptide (CGRP) is a neuropeptide that is abundant in the sensory neurons which innervate bone. The effects of CGRP on isolated bone cells have been widely studied, and CGRP is currently considered to be an osteoanabolic peptide that has effects on both osteoclasts and osteoblasts. However, relatively little is known about the physiological role of CGRP in-vivo in the skeletal responses to bone loading, particularly fatigue loading.

Methodology/Principal Findings

We used the rat ulna end-loading model to induce fatigue damage in the ulna unilaterally during cyclic loading. We postulated that CGRP would influence skeletal responses to cyclic fatigue loading. Rats were fatigue loaded and groups of rats were infused systemically with 0.9% saline, CGRP, or the receptor antagonist, CGRP_8–37, for a 10 day study period. Ten days after fatigue loading, bone and serum CGRP concentrations, serum tartrate-resistant acid phosphatase 5b (TRAP5b) concentrations, and fatigue-induced skeletal responses were quantified. We found that cyclic fatigue loading led to increased CGRP concentrations in both loaded and contralateral ulnae. Administration of CGRP_8–37 was associated with increased targeted remodeling in the fatigue-loaded ulna. Administration of CGRP or CGRP_8–37 both increased reparative bone formation over the study period. Plasma concentration of TRAP5b was not significantly influenced by either CGRP or CGRP_8–37 administration.

Conclusions

CGRP signaling modulates targeted remodeling of microdamage and reparative new bone formation after bone fatigue, and may be part of a neuronal signaling pathway which has regulatory effects on load-induced repair responses within the skeleton.

Experimental and finite element analysis of the mouse caudal vertebrae loading model: prediction of cortical and trabecular bone adaptation

In this study, we attempt to predict cortical and trabecular bone adaptation in the mouse caudal vertebrae loading model using knowledge of bone's local mechanical environment at the onset of loading. In a previous study, we demonstrated appreciable 25.9 and 11% increases in both trabecular and cortical bone volume density, respectively, when subjecting the fifth caudal vertebrae (C5) of C57BL/6 (B6) mice to an acute loading regime (amplitude of 8N, 3000 cycles, 10 Hz, 3 times a week for 4 weeks). We have also established a validated finite element (FE) model of the C5 vertebra using micro-computed tomography (micro-CT), which characterizes, in 3D, the micro-mechanical strains present in both cortical and trabecular compartments due to the applied loads. To investigate the relationship between load-induced bone adaptation and mechanical strains in-vivo and in-silico data sets were compared. Using data from the previous cross-sectional study, we divided cortical and trabecular compartments into 15 subregions and determined, for each region, a bone formation parameter ΔBV/BS (a cross-sectional measure of the bone volume added to cortical and trabecular surfaces following the described loading regime). Linear regression was then used to correlate mean regional values of ΔBV/BS with mean values of mechanical strains derived from the FE models which were similarly regionalized. The mechanical parameters investigated were strain energy density (SED), the orthogonal strains (e (x), e (y), e (z)) and the three shear strains (e (xy), e (yz), e (zx)). For cortical regions, regression analysis showed SED to correlate extremely well with ΔBV/BS (R (2) = 0.82) and e (z) (R (2) = 0.89). Furthermore, SED was found to predict expansion of the cortical shell correlating significantly with the regional percentage increases in cortical tissue volume (R (2) = 0.92), cortical marrow volume (R (2) = 0.91) and cortical thickness (R (2) = 0.56). For trabecular regions, FE parameters were found not to correlate with load-induced trabecular bone morphology. These results indicate that load-induced cortical morphology can be predicted from population data, whereas the prediction of trabecular morphology requires subject-specific micro- architecture.

Experimentally tested computer modeling of stress fractures in rats

The objective of this study was to develop a finite-element (FE) modeling methodology for studying the etiology of a stress fracture (SF). Several variants of three-dimensional FE models of a rat hindlimb, which differed in length or stiffness of tissues, enabling the analyses of mechanical strains and stress in the tibia, were created. We compared the occurrence of SFs in an animal model to validate locations of peak strains/stresses in the FE models. Four Sprague-Dawley male rats, age ∼7 wk, were subjected to mechanical cyclic loads of 1.2 Hz and ∼6 N, which were delivered to their hindlimb for 30 min, 3 times/wk, up to 12 wk, by using a specially designed apparatus. The results showed that 1) FE modeling predicted the maximal strains/stresses (∼220,0 με and ∼29 MPa, respectively) between the mid- and proximal thirds of the tibia; 2) in a longer shin, greater and more inhomogeneous tensile strains/stresses were evident, at the same location; 3) anatomical variants in shin length influenced the strain/stress distributions to a greater extent with respect to changes in mechanical properties of tissues; and 4) bone stiffness was more dominant than muscle stiffness in affecting the strain/stress distributions. In the animal study, 35,000 loading cycles were associated with the formation of a SF. The location of the identified SF in the rat limb verified the FE model. We find the suggested model a valuable tool in studying various aspects of SFs.

Hierarchies of damage induced loss of mechanical properties in calcified bone after in vivo fatigue loading of rat ulnae

During fatigue loading of whole bone, damage to bone tissue accumulates, coalesces and leads to fractures. Whether damage affects tissue material properties similarly at the nanoscale (less than 1 μm), microscale (less than 1 mm), and whole bone scale has not been fully evaluated. Therefore, in this study, we examine scale-dependent loss of calcified tissue material properties in rat ulnae, after fatigue loading of rat forearms using the forearm compression model. In vivo fatigue loading was conducted on the right forearms until a displacement end-point was reached. The non-fatigued left forearms served as contralateral controls. Subsequently, three-point bending tests to failure on excised ulnae demonstrated a 41% and 49% reduction in the stiffness and ultimate strength as compared to contralateral control ulnae, respectively. Depth-sensing microindentation demonstrated an average decrease in material properties, such as elastic modulus and hardness, of 28% and 29% respectively. Nanoindentation measured elastic modulus and hardness were reduced by 26% and 29% in damaged bone relative to contralateral controls, respectively. The increased loss of whole bone material properties compared to tissue material properties measured using indentation is mainly attributed to the presence of a macrocrack located in the medial compressive region at the site of peak strains. The similar magnitude of changes in material properties by microindentation and nanoindentation is attributed to damage that may originate at an even smaller scale, as inferred from 10% differences in connectivity of osteocyte canaliculi in damaged bone.

Activation of bone remodeling after fatigue: differential response to linear microcracks and diffuse damage

Recent experiments point to two predominant forms of fatigue microdamage in bone: linear microcracks (tens to a few hundred microns in length) and "diffuse damage" (patches of diffuse stain uptake in fatigued bone comprised of clusters of sublamellar-sized cracks). The physiological relevance of diffuse damage in activating bone remodeling is not known. In this study microdamage amount and type were varied to assess whether linear or diffuse microdamage has similar effects on the activation of intracortical resorption. Activation of resorption was correlated to the number of linear microcracks (Cr.Dn) in the bone (R(2)=0.60, p<0.01). In contrast, there was no activation of resorption in response to diffuse microdamage alone. Furthermore, there was no significant change in osteocyte viability in response to diffuse microdamage, suggesting that osteocyte apoptosis, which is known to activate remodeling at typical linear microcracks in bone, does not result from sublamellar damage. These findings indicate that inability of diffuse microdamage to activate resorption may be due to lack of a focal injury response. Finally, we found that duration of loading does not affect the remodeling response. In conclusion, our data indicate that osteocytes activate resorption in response to linear microcracks but not diffuse microdamage, perhaps due to lack of a focal injury-induced apoptotic response.

Systemic effects of ulna loading in male rats during functional adaptation

Functional skeletal adaptation is thought to be a local phenomenon controlled by osteoctyes. However, the nervous system also may have regulatory effects on adaptation. The aim of this study was to determine the effects of loading of a single bone on adaptation of other appendicular long bones and whether these responses were neuronally regulated. Young male Sprague-Dawley rats were used. The right ulna was loaded to induce a modeling response. In other rats, a second regimen was used to induce bone fatigue with a mixed modeling/remodeling response; a proportion of rats from each group received brachial plexus anesthesia to induce temporary neuronal blocking during bone loading. Sham groups were included. Left and right long bones (ulna, humerus, tibia, and femur) from each rat were examined histologically 10 days after loading. In fatigue- and sham-loaded animals, blood plasma concentrations of TNF-α, RANKL, OPG, and TRAP5b were determined. We found that loading the right ulna induced an increase in bone formation in distant long bones that were not loaded and that this effect was neuronally regulated. Distant effects were most evident in the rats that received loading without bone fatigue. In the fatigue-loaded animals, neuronal blocking induced a significant decrease in plasma TRAP5b at 10 days. Histologically, bone resorption was increased in both loaded and contralateral ulnas in fatigue-loaded rats and was not significantly blocked by brachial plexus anesthesia. In young, growing male rats we conclude that ulna loading induced increased bone formation in multiple bones. Systemic adaptation effects were, at least in part, neuronally regulated.

An approach to quantifying bone overloading and hypertrophy with applications to multiple experimental studies

Many studies have investigated mechanically induced bone formation in mice and rats by applying loads to the long bones, and measuring changes in periosteal cortical bone apposition rates. However, the results are difficult to compare among each other because the loading schemes are generally different. The purpose of the present study was to develop a theoretical framework for evaluating the mechanical stimulus based on the bone daily strain stimulus, which is a function of loading cycles and bone strains. The daily strain stimulus would act as a single unifying parameter for directly comparing data from existing in vivo experiments, and is applied here to twenty previous rat and mouse studies. To calculate the daily strain stimulus, we determined the periosteal daily strain stimulus necessary for bone maintenance (xi(peri,0)) and the strain-cycle weighting exponent (m). In the first approach, we applied data from Rubin and Lanyon's bone maintenance studies. We calculated xi(peri,0) to be 2793 microstrain/day, and m to be 4.5. In the second approach, we used Fritton et al. 's strain gage recordings to calculate xi(peri,0) to be 1496 microstrain/day, and used an m value of 11.88, equal to human bone compressive fatigue properties. Fatigue data provided physiological relevance, and was useful for applying an established remodeling theory to in vivo studies. For both approaches, xi(peri,0) was below the fracture level. We then analyzed the applied strains, cycles, and periosteal bone apposition rates from the previous studies. The range of daily strain stimuli calculated using the first approach was much larger than the range using the second approach (2793-17312 microstrain/day compared to 1496-7681 microstrain/day). None of the studies applied a daily strain stimulus above the complete fatigue failure level, but some studies applied loading that could result in major fatigue microdamage. Bone apposition rates generally increased with increasing daily strain stimulus, which was consistent with previous theoretical models. The results suggest that the daily strain stimulus may be a reasonable first approximation for predicting bone apposition rates in a consistent manner. The use of the daily strain stimulus may be helpful for improving the design of future bone loading studies.