Geometric mouse variation: Implications to the axial ulnar loading protocol and animal specific calibration

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

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

Naproxen impairs load-induced bone formation, reduces bone toughness, and diminishes woven bone formation following stress fracture in mice

Debilitating stress fractures are surprisingly common in physically active individuals, including athletes, military recruits, and dancers. These individuals are overrepresented in the 30 million daily users of non-steroidal anti-inflammatory drugs (NSAIDs). We hypothesized that regular use of NSAIDs would predispose habitually loaded bones to stress fracture and delay the repair of these injuries. In this project, we used repetitive axial forelimb compression in mice as a model to test these hypotheses. First, adult mice were subjected to six bouts of forelimb compression over a period of two weeks, with aspirin, naproxen, or vehicle continuously administered through drinking water. Naproxen-treated mice had diminished load-induced bone formation as well as a significant loss in toughness in non-loaded bone, which were not observed in aspirin-treated mice. Furthermore, there were no differences in RANKL/OPG ratio or cortical bone parameters. Picrosirius red staining and second harmonic generation imaging revealed that alterations in bone collagen fibril size and organization were driving the loss of toughness in naproxen-treated mice. Separately, adult mice were subjected to an ulnar stress fracture generated by a single bout of fatigue loading, with NSAIDs provided 24 h before injury. Both aspirin-treated and naproxen-treated mice had normal forelimb use in the week after injury, whereas control mice favored the injured forelimb until day 7. However, woven bone volume was only significantly impaired by naproxen. Both NSAIDs were found to significantly inhibit Ptgs2 and Ngf expression following stress fracture, but only naproxen significantly affected serum PGE2 concentration. Overall, our results suggest that naproxen, but not aspirin, may increase the risk of stress fracture and extend the healing time of these injuries, warranting further clinical evaluation for patients at risk for fatigue injuries.

Comparison of strain measurement in the mouse forearm using subject-specific finite element models, strain gaging, and digital image correlation

Mechanical loading in bone leads to the activation of bone-forming pathways that are most likely associated with a minimum strain threshold being experienced by the osteocyte. To investigate the correlation between cellular response and mechanical stimuli, researchers must develop accurate ways to measure/compute strain both externally on the bone surface and internally at the osteocyte level. This study investigates the use of finite element (FE) models to compute bone surface strains on the mouse forearm. Strains from three FE models were compared to data collected experimentally through strain gaging and digital image correlation (DIC). Each FE model was assigned subject-specific bone properties and consisted of one-dimensional springs representing the interosseous membrane. After three-point bending was performed on the ulnae and radii, moment of inertia was determined from microCT analysis of the bone region between the supports and then used along with standard beam analyses to calculate the Young's modulus. Non-contact strain measurements from DIC were determined to be more suitable for validating numerical results than experimental data obtained through conventional strain gaging. When comparing strain responses in the three ulnae, we observed a 3-14% difference between numerical and DIC strains while the strain gage values were 37-56% lower than numerical values. This study demonstrates a computational approach for capturing bone surface strains in the mouse forearm. Ultimately, strains from these macroscale models can be used as inputs for microscale and nanoscale FE models designed to analyze strains directly in the osteocyte lacunae.

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.