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.

[1]  Matthew J. Silva,et al.  Age-Related Changes in Bone Structure and Strength in Female and Male BALB/c Mice , 2010, Calcified Tissue International.

[2]  Felix Eckstein,et al.  Non-invasive axial loading of mouse tibiae increases cortical bone formation and modifies trabecular organization: a new model to study cortical and cancellous compartments in a single loaded element. , 2005, Bone.

[3]  R. Zernicke,et al.  Strain Gradients Correlate with Sites of Exercise‐Induced Bone‐Forming Surfaces in the Adult Skeleton , 1997, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[4]  David L. Wilson,et al.  Voxel size and measures of individual resorption cavities in three-dimensional images of cancellous bone. , 2009, Bone.

[5]  S. Cowin,et al.  A model for the excitation of osteocytes by mechanical loading-induced bone fluid shear stresses. , 1994, Journal of biomechanics.

[6]  Craig R. Slyfield,et al.  Mechanical failure begins preferentially near resorption cavities in human vertebral cancellous bone under compression. , 2012, Bone.

[7]  Ling Qin,et al.  3D image registration is critical to ensure accurate detection of longitudinal changes in trabecular bone density, microstructure, and stiffness measurements in rat tibiae by in vivo microcomputed tomography (μCT). , 2013, Bone.

[8]  Alexander G Robling,et al.  Improved Bone Structure and Strength After Long‐Term Mechanical Loading Is Greatest if Loading Is Separated Into Short Bouts , 2002, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[9]  P. Prendergast,et al.  A mechano-regulation model for tissue differentiation during fracture healing: analysis of gap size and loading. , 2002, Journal of biomechanics.

[10]  Georg N Duda,et al.  Monitoring in vivo (re)modeling: a computational approach using 4D microCT data to quantify bone surface movements. , 2015, Bone.

[11]  Sandra J. Shefelbine,et al.  Ex vivo determination of bone tissue strains for an in vivo mouse tibial loading model , 2014, Journal of biomechanics.

[12]  R. Müller,et al.  Experimental and finite element analysis of the mouse caudal vertebrae loading model: prediction of cortical and trabecular bone adaptation , 2012, Biomechanics and modeling in mechanobiology.

[13]  Kyle E. Niemeyer,et al.  Three‐dimensional surface texture visualization of bone tissue through epifluorescence‐based serial block face imaging , 2009, Journal of microscopy.

[14]  Ralph Müller,et al.  Mineralization kinetics in murine trabecular bone quantified by time-lapsed in vivo micro-computed tomography. , 2013, Bone.

[15]  Matthew J Silva,et al.  Aging diminishes lamellar and woven bone formation induced by tibial compression in adult C57BL/6. , 2014, Bone.

[16]  A. Pitsillides,et al.  Using digital image correlation to determine bone surface strains during loading and after adaptation of the mouse tibia. , 2010, Journal of biomechanics.

[17]  Kent D. Butz,et al.  Characterization of cancellous and cortical bone strain in the in vivo mouse tibial loading model using microCT-based finite element analysis. , 2014, Bone.

[18]  Melissa L Knothe Tate,et al.  Periosteum, bone's “smart” bounding membrane, exhibits direction‐dependent permeability , 2013, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[19]  Georg N Duda,et al.  Skeletal maturity leads to a reduction in the strain magnitudes induced within the bone: a murine tibia study. , 2015, Acta biomaterialia.

[20]  S. Cowin,et al.  A model for strain amplification in the actin cytoskeleton of osteocytes due to fluid drag on pericellular matrix. , 2001, Journal of biomechanics.

[21]  Craig R. Slyfield,et al.  Three-dimensional characterization of resorption cavity size and location in human vertebral trabecular bone. , 2012, Bone.

[22]  C. Hernandez,et al.  The effects of misalignment during in vivo loading of bone: techniques to detect the proximity of objects in three-dimensional models. , 2014, Journal of biomechanics.

[23]  Kevin W Eliceiri,et al.  NIH Image to ImageJ: 25 years of image analysis , 2012, Nature Methods.

[24]  Stuart J Warden,et al.  Cortical and trabecular bone adaptation to incremental load magnitudes using the mouse tibial axial compression loading model. , 2013, Bone.

[25]  Sara Checa,et al.  Tomography-Based Quantification of Regional Differences in Cortical Bone Surface Remodeling and Mechano-Response , 2017, Calcified Tissue International.

[26]  Matthew J. Silva,et al.  Tibial Loading Increases Osteogenic Gene Expression and Cortical Bone Volume in Mature and Middle-Aged Mice , 2012, PloS one.

[27]  Matthew J Silva,et al.  Experimental and finite element analysis of strains induced by axial tibial compression in young-adult and old female C57Bl/6 mice. , 2014, Journal of biomechanics.

[28]  Ralph Müller,et al.  In vivo micro-computed tomography allows direct three-dimensional quantification of both bone formation and bone resorption parameters using time-lapsed imaging. , 2011, Bone.

[29]  Georg N Duda,et al.  Aging Leads to a Dysregulation in Mechanically Driven Bone Formation and Resorption , 2015, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[30]  S. J. Shefelbine,et al.  Predicting cortical bone adaptation to axial loading in the mouse tibia , 2015, Journal of The Royal Society Interface.

[31]  G J Kazakia,et al.  Automated high‐resolution three‐dimensional fluorescence imaging of large biological specimens , 2007, Journal of microscopy.

[32]  Craig R. Slyfield,et al.  Three‐dimensional dynamic bone histomorphometry , 2012, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[33]  E. V. Tkachenko,et al.  Anti-resorptive agents reduce the size of resorption cavities: a three-dimensional dynamic bone histomorphometry study. , 2013, Bone.

[34]  Mary Beth Tribble,et al.  μCT-based, in vivo dynamic bone histomorphometry allows 3D evaluation of the early responses of bone resorption and formation to PTH and alendronate combination therapy. , 2015, Bone.

[35]  R. Müller,et al.  Longitudinal Assessment of In Vivo Bone Dynamics in a Mouse Tail Model of Postmenopausal Osteoporosis , 2011, Calcified Tissue International.

[36]  L. E. Lanyon,et al.  Mechanical loading-related changes in osteocyte sclerostin expression in mice are more closely associated with the subsequent osteogenic response than the peak strains engendered , 2011, Osteoporosis International.

[37]  A. Pitsillides,et al.  Characterizing a novel and adjustable noninvasive murine joint loading model. , 2011, Arthritis and rheumatism.

[38]  C. Rubin,et al.  Strain Gradients Correlate with Sites of Periosteal Bone Formation , 1997, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[39]  Georg N Duda,et al.  The Periosteal Bone Surface is Less Mechano-Responsive than the Endocortical , 2016, Scientific Reports.

[40]  André F. Pereira Cortical Bone Adaptation: A Finite-Element Study of the Mouse Tibia , 2014 .

[41]  Andrew H. Gee,et al.  Cortical Thickness Mapping to Identify Focal Osteoporosis in Patients with Hip Fracture , 2012, PloS one.

[42]  R. Main,et al.  Load-induced changes in bone stiffness and cancellous and cortical bone mass following tibial compression diminish with age in female mice , 2014, Journal of Experimental Biology.

[43]  Georg N Duda,et al.  Mineralizing surface is the main target of mechanical stimulation independent of age: 3D dynamic in vivo morphometry. , 2014, Bone.