Tissue strain amplification at the osteocyte lacuna: a microstructural finite element analysis

Assessment of submicron bone tissue composition in plastic-embedded samples using optical photothermal infrared (O-PTIR) spectral imaging and machine learning

Understanding the composition of bone tissue at the submicron level is crucial to elucidate factors contributing to bone disease and fragility. Here, we introduce a novel approach utilizing optical photothermal infrared (O-PTIR) spectroscopy and imaging coupled with machine learning analysis to assess bone tissue composition at 500 nm spatial resolution. This approach was used to evaluate thick bone samples embedded in typical poly(methyl methacrylate) (PMMA) blocks, eliminating the need for cumbersome thin sectioning. We demonstrate the utility of O-PTIR imaging to assess the distribution of bone tissue mineral and protein, as well as to explore the structure-composition relationship surrounding microporosity at a spatial resolution unattainable by conventional infrared imaging modalities. Using bone samples from wildtype (WT) mice and from a mouse model of osteogenesis imperfecta (OIM), we further showcase the application of O-PTIR spectroscopy to quantify mineral content, crystallinity, and carbonate content in spatially defined regions across the cortical bone. Notably, we show that machine learning analysis using support vector machine (SVM) was successful in identifying bone phenotypes (typical in WT, fragile in OIM) based on input of spectral data, with over 86 % of samples correctly identified when using the collagen spectral range. Our findings highlight the potential of O-PTIR spectroscopy and imaging as valuable tools for exploring bone submicron composition.

Evaluating the influence on osteocyte mechanobiology within the lacunar-canalicular system for varying lacunar equancy and perilacunar elasticity: A multiscale fluid-structure interaction analysis

The lacunar morphology and perilacunar tissue properties of osteocytes in bone can vary under different physiological and pathological conditions. How these alterations collectively change the overall micromechanics of osteocytes in the lacunar-canalicular system (LCS) of an osteon still requires special focus. Therefore, a Haversian canal and LCS-based osteon model was established to evaluate the changes in the hydrodynamic environment around osteocytes under physiological loading using fluid-structure interaction analysis, followed by a sub-modelled finite element analysis to assess the mechanical responses of osteocytes and their components. Osteocytes were modelled with detailed configurations, including cytoplasm, nucleus, and cytoskeleton, and parametric variations in lacunar equancy (L.Eq) and perilacunar elasticity (Pl.E) were considered within the osteon model. The study aimed to conduct a comparative study among osteon models with varying L. Eq and Pl. E to check the resulting differences in osteocyte mechanobiology. The results demonstrated that the average mechanical stimulation of each subcellular component of osteocytes increased with decreases in L. Eq and Pl. E, reflecting conditions typically seen in young, healthy bone as per previous literature. However, hydrodynamic responses, such as fluid flow and fluid shear stress on osteocytes, varied proportionally with the elasticity difference between the bone matrix and the perilacunar region during Pl. E variation. Additionally, the findings revealed that a minimal percentage of energy was used to transmit mechanical responses through microtubules from the cell membrane to the nucleus, and this energy percentage increased with higher L. Eq. The outcomes of the study could help to quantify how the osteocyte microenvironment and its mechanosensitivity within cortical bone changes with L. Eq and Pl. E alterations in different bone conditions, from young to aged and healthy to diseased.

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

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

Application of a combined cancellous lag screw enhances the stability of locking plate fixation of osteoporotic lateral tibial plateau fracture by providing interfragmentary compression force

BACKGROUND

Insufficient interfragmentary compression force (IFCF) frequently leads to unstable fixation of osteoporotic lateral tibial plateau fractures (OLTPFs). A combined cancellous lag screw (CCLS) enhances IFCF; however, its effect on OLTPF fixation stability remains unclear. Therefore, we investigated the effect of CCLS on OLTPF stability using locking plate fixation (LPF).

MATERIALS AND METHODS

Twelve synthetic osteoporotic tibial bones were used to simulate OLTPFs, which were fixed using LPF, LPF-AO cancellous lag screws (LPF-AOCLS), and LPF-CCLS. Subsequently, 10,000 cyclic loadings from 30 to 400 N were performed. The initial axial stiffness (IAS), maximal axial micromotion of the lateral fragment (MAM-LF) measured every 1000 cycles, and failure load after 10,000 cycles were tested. The same three fixations for OLTPF were simulated using finite element analysis (FEA). IFCFs of 0, 225, and 300 N were applied to the LPF, LPF-AOCLS, and LPF-CCLS, respectively, with a 1000-N axial compressive force. The MAM-LF, peak von Mises stress (VMS), peak equivalent elastic strain of the lateral fragment (EES-LF), and nodes of EES-LF > 2% (considered bone destruction) were calculated.

RESULTS

Biomechanical tests revealed the LPF-AOCLS and LPF-CCLS groups to be superior to the LPF group in terms of the IAS, MAM-LF, and failure load (all p < 0.05). FEA revealed that the MAM-LF, peak VMS, peak EES-LF, and nodes with EES-LF > 2% in the LPF were higher than those in the LPF-AOCLS and LPF-CCLS.

CONCLUSION

IFCF was shown to enhance the stability of OLTPFs using LPF. Considering overscrewing, CCLS is preferably recommended, although there were no significant differences between CCLS and AOCLS.

Computational modeling based on confocal imaging predicts changes in osteocyte and dendrite shear stress due to canalicular loss with aging

Exercise and physical activity exert mechanical loading on the bones which induces bone formation. However, the relationship between the osteocyte lacunar-canalicular morphology and mechanical stress experienced locally by osteocytes transducing signals for bone formation is not fully understood. In this study, we used computational modeling to predict the effect of canalicular density, the number of fluid inlets, and load direction on fluid flow shear stress (FFSS) and bone strains and how these might change following the microstructural deterioration of the lacunar-canalicular network that occurs with aging. Four distinct computational models were initially generated of osteocytes with either ten or eighteen dendrites using a fluid-structure interaction method with idealized geometries. Next, a young and a simulated aged osteocyte were developed from confocal images after FITC staining of the femur of a 4-month-old C57BL/6 mouse to estimate FFSS using a computational fluid dynamics approach. The models predicted higher fluid velocities in the canaliculi versus the lacunae. Comparison of idealized models with five versus one fluid inlet indicated that with four more inlets, one-half of the dendrites experienced FFSS greater than 0.8 Pa, which has been associated with osteogenic responses. Confocal image-based models of real osteocytes indicated a six times higher ratio of canalicular to lacunar surface area in the young osteocyte model than the simulated aged model and the average FFSS in the young model (FFSS = 0.46 Pa) was three times greater than the aged model (FFSS = 0.15 Pa). Interestingly, the surface area with FFSS values above 0.8 Pa was 23 times greater in the young versus the simulated aged model. These findings may explain the impaired mechano-responsiveness of osteocytes with aging.

The influence of different shaped osteocyte lacunae on microcrack initiation and propagation

BACKGROUND

The morphology of osteocyte lacunae varies in bones of different ages and bone pathologies. Osteocyte lacunae can cause stress concentration and initiate microcracks. However, the influence of changes in osteocyte lacunar shape on microcrack is unknown. Therefore, the aim of this study was to determine the effects of osteocyte lacunae with different shapes on microcrack initiation and propagation.

METHODS

Osteon models containing osteocyte lacunae with different shapes were established. The progressive damage analysis method, based on computer simulations, was used to study the evolution of microdamage within the osteon, including the processes of intralaminar and interlaminar microdamage.

FINDINGS

Models with larger DoE values can effectively delay or prevent the formation of linear microcracks, which ensures high fracture toughness of cortical bone. It is subjected to stronger mechanical stimulation, making it more sensitive to loads. Models with smaller DoE values increase the load threshold for microdamage generation and reduces its impact on bone mechanical performance, making it less susceptible to microdamage than models with larger DoE values.

INTERPRETATION

These findings enhance the limited knowledge of the influence of the lacunar shape on microdamage and contribute to a better understanding of bone biomechanics.

Osteoporosis and Covid-19: Detected similarities in bone lacunar-level alterations via combined AI and advanced synchrotron testing

While advanced imaging strategies have improved the diagnosis of bone-related pathologies, early signs of bone alterations remain difficult to detect. The Covid-19 pandemic has brought attention to the need for a better understanding of bone micro-scale toughening and weakening phenomena. This study used an artificial intelligence-based tool to automatically investigate and validate four clinical hypotheses by examining osteocyte lacunae on a large scale with synchrotron image-guided failure assessment. The findings indicate that trabecular bone features exhibit intrinsic variability related to external loading, micro-scale bone characteristics affect fracture initiation and propagation, osteoporosis signs can be detected at the micro-scale through changes in osteocyte lacunar features, and Covid-19 worsens micro-scale porosities in a statistically significant manner similar to the osteoporotic condition. Incorporating these findings with existing clinical and diagnostic tools could prevent micro-scale damages from progressing into critical fractures.

Evaluating the Role of Canalicular Morphology and Perilacunar Region Properties on Local Mechanical Environment of Lacunar-Canalicular Network Using Finite Element Modeling

Physiological and pathological processes such as aging, diseases, treatments, and lactation can alter lacunar-canalicular network (LCN) morphology and perilacunar region properties. These modifications can impact the mechanical environment of osteocytes which in turn can influence osteocyte mechanosensitivity and the remodeling process. In this study, we aim to evaluate how the modifications in the canalicular morphology, lacunar density, and the perilacunar region properties influence the local mechanical environment of LCN and the apparent bone properties using three-dimensional finite element (FE) modeling. The simulation results showed that a 50% reduction in perilacunar elastic modulus led to about 7% decrease in apparent elastic modulus of the bone. The increase in canalicular density, length, and diameter did not influence the strain amplification in the models but they increased the amount of highly strained bone around LCN. Change in lacunar density did not influence the strain amplification and the amount of highly strained regions on LCN surfaces. Reduction in perilacunar elastic modulus increased both the strain amplification and the volume of highly strained tissue around and on the surface of LCN. The FE models of LCN in this study can be utilized to quantify the influence of modifications in canalicular morphology, lacunar density, and perilacunar region properties on the apparent bone properties and the local mechanical environment of LCN. Although this is a numerical study with idealized models, it provides important information on how mechanical environment of osteocytes is influenced by the modifications in LCN morphology and perilacunar region properties due to physiological and pathological processes.

Oscillating Fluid Flow Activated Osteocyte Lysate-Based Hydrogel for Regulating Osteoblast/Osteoclast Homeostasis to Enhance Bone Repair

As major regulators on bone formation/resorption in response to mechanical stimuli, osteocytes have shown great promise for restoring bone injury. However, due to the unmanageable and unabiding cell functions in unloading or diseased environments, the efficacy of osteogenic induction by osteocytes has been enormously limited. Herein, a facile method of oscillating fluid flow (OFF) loading for cell culture is reported, which enables osteocytes to initiate only osteogenesis and not the osteolysis process. After OFF loading, multiple and sufficient soluble mediators are produced in osteocytes, and the collected osteocyte lysates invariably induce robust osteoblastic differentiation and proliferation while restraining osteoclast generation and activity under unloading or pathological conditions. Mechanistic studies confirm that elevated glycolysis and activation of the ERK1/2 and Wnt/β-catenin pathways are the major contributors to the initiation of osteoinduction functions induced by osteocytes. Moreover, an osteocyte lysate-based hydrogel is designed to establish a stockpile of "active osteocytes" to sustainably deliver bioactive proteins, resulting in accelerated healing through regulation of endogenous osteoblast/osteoclast homeostasis.

Using Finite Element Modeling in Bone Mechanoadaptation

PURPOSE OF THE REVIEW

Bone adapts structure and material properties in response to its mechanical environment, a process called mechanoadpatation. For the past 50 years, finite element modeling has been used to investigate the relationships between bone geometry, material properties, and mechanical loading conditions. This review examines how we use finite element modeling in the context of bone mechanoadpatation.

RECENT FINDINGS

Finite element models estimate complex mechanical stimuli at the tissue and cellular levels, help explain experimental results, and inform the design of loading protocols and prosthetics. FE modeling is a powerful tool to study bone adaptation as it complements experimental approaches. Before using FE models, researchers should determine whether simulation results will provide complementary information to experimental or clinical observations and should establish the level of complexity required. As imaging technics and computational capacity continue increasing, we expect FE models to help in designing treatments of bone pathologies that take advantage of mechanoadaptation of bone.

Both network architecture and micro cracks effects on lacuno-canalicular liquid flow efficiency within the context of multiphysics approach for bone remodeling

When physical forces are applied to bone, its mechanical adaptive behaviors change according to the microarchitecture configuration. This leads to changes in biological and physical thresholds in the remodeling cell population, involving sensor cells (osteocytes) interacting with each other and changes in osteocyte shape due to variation in lacunar shape. The resulting alterations in fluid flow leads to changes in the membrane electrical potential and shear stress. Eventual creation of microcracks, may lead in turn to modify cell activity. In contrast, the redundancy in the lacuno canalicular network (LCN) interconnectivity maintains partial flow. Our goal was to investigate the role of fluid flow in LCN by proposing a model of electro-mechanical energy spread through inhomogeneous microarchitectures. We focused on mechano-sensitivity to changes in load-induced flow impacted by neighboring micro cracks and quantifying its critical role in changing, velocity, shear stress and orientation of liquid mass transportation from one cell to another. To enhance the concept of intricacy LCN micro-structure to fluid flow, we provide a new combined effects factor considered as osteocytes sensor efficiency. We customized an influence function for each osteocyte, coupling: in one hand, the spatial distribution within remodeling influence areas, conducting a significant fluid spread, leading hydro-dynamic behavior and impacted further by presence of micro cracks and; in other hand, the fluid electro kinetic behavior. As an attempt to fill the limitations stated by many of the recent studies, we reveal in numerical simulation, some results which cannot be measured in vitro/in vivo studies. Numerical calculations were performed in order to evaluate, among many others, how liquid flow conditions changes between lacunas, how the orientation and the magnitude of the governing flow in LCN can regulate osteocytes efficiency. In addition to be regulated by osteocytes, a direct effects of fluid flow are also acting on osteoblast activity. In summary, this new approach considers mechano-sensitivity in relation to liquid flow dynamic and suggests additional pathway for Osseo integration via osteoblast regulation. However, this novel modeling approach may help improve the mapping and design bone scaffolds and/or selection of scaffold implantation regions.

Isolating the Role of Bone Lacunar Morphology on Static and Fatigue Fracture Progression through Numerical Simulations

Currently, the onset of bone damage and the interaction of cracks with the surrounding micro-architecture are still black boxes. With the motivation to address this issue, our research targets isolating lacunar morphological and densitometric effects on crack advancement under both static and cyclic loading conditions by implementing static extended finite element models (XFEM) and fatigue analyses. The effect of lacunar pathological alterations on damage initiation and progression is evaluated; the results indicate that high lacunar density considerably reduces the mechanical strength of the specimens, resulting as the most influencing parameter among the studied ones. Lacunar size has a lower effect on mechanical strength, reducing it by 2%. Additionally, specific lacunar alignments play a key role in deviating the crack path, eventually slowing its progression. This could shed some light on evaluating the effects of lacunar alterations on fracture evolution in the presence of pathologies.

Osteocyte Remodeling of the Lacunar-Canalicular System: What's in a Name?

PURPOSE OF REVIEW

Osteocytes directly modify the bone surrounding the expansive lacunar-canalicular system (LCS) through both resorption and deposition. The existence of this phenomenon is now widely accepted, but is referred to as "osteocyte osteolysis," "LCS remodeling," and "perilacunar remodeling," among other names. The uncertainty in naming this physiological process reflects the many persistent questions about why and how osteocytes interact with local bone matrix. The goal of this review is to examine the purpose and nature of LCS remodeling and its impacts on multiscale bone quality.

RECENT FINDINGS

While LCS remodeling is clearly important for systemic calcium mobilization, this process may have additional potential drivers and may impact the ability of bone to resist fracture. There is abundant evidence that the osteocyte can resorb and replace bone mineral and does so outside of extreme challenges to mineral homeostasis. The impacts of the osteocyte on organic matrix are less certain, especially regarding whether osteocytes produce osteoid. Though multiple lines of evidence point towards osteocyte production of organic matrix, definitive work is needed. Recent high-resolution imaging studies demonstrate that LCS remodeling influences local material properties. The role of LCS remodeling in the maintenance and deterioration of bone matrix quality in aging and disease are active areas of research. In this review, we highlight current progress in understanding why and how the osteocyte removes and replaces bone tissue and the consequences of these activities to bone quality. We posit that answering these questions is essential for evaluating whether, how, when, and why LCS remodeling may be manipulated for therapeutic benefit in managing bone fragility.

On the fracture behavior of cortical bone microstructure: The effects of morphology and material characteristics of bone structural components

Bone encompasses a complex arrangement of materials at different length scales, which endows it with a range of mechanical, chemical, and biological capabilities. Changes in the microstructure and characteristics of the material, as well as the accumulation of microcracks, affect the bone fracture properties. In this study, two-dimensional finite element models of the microstructure of cortical bone were considered. The eXtended Finite Element Method (XFEM) developed by Abaqus software was used for the analysis of the microcrack propagation in the model as well as for local sensitivity analysis. The stress-strain behavior obtained for the different introduced models was substantially different, confirming the importance of bone tissue microstructure for its failure behavior. Considering the role of interfaces, the results highlighted the effect of cement lines on the crack deflection path and global fracture behavior of the bone microstructure. Furthermore, bone micromorphology and areal fraction of cortical bone tissue components such as osteons, cement lines, and pores affected the bone fracture behavior; specifically, pores altered the crack propagation path since increasing porosity reduced the maximum stress needed to start crack propagation. Therefore, cement line structure, mineralization, and areal fraction are important parameters in bone fracture. The parameter-wise sensitivity analysis demonstrated that areal fraction and strain energy release rate had the greatest and the lowest effect on ultimate strength, respectively. Furthermore, the component-wise sensitivity analysis revealed that for the areal fraction parameter, pores had the greatest effect on ultimate strength, whereas for the other parameters such as elastic modulus and strain energy release rate, cement lines had the most important effect on the ultimate strength. In conclusion, the finding of the current study can help to predict the fracture mechanisms in bone by taking the morphological and material properties of its microstructure into account.

Large-scale osteocyte lacunar morphological analysis of transiliac bone in normal and osteoporotic premenopausal women

Bone's ability to adapt is governed by the network of embedded osteocytes, which inhabit individual pores called lacunae. The morphology of these lacunae and their resident osteocytes are known to change with age and diseases such as postmenopausal osteoporosis. However, it is unclear whether alterations in lacunar morphology are present in younger populations with osteoporosis. To investigate this, we implemented a previously validated methodology to image and quantify the three-dimensional morphometries of lacunae on a large scale with ultra-high-resolution micro-computed tomography (microCT) in transiliac bone biopsies from three groups of premenopausal women: control n = 39; idiopathic osteoporosis (IOP) n = 45; idiopathic low BMD (ILBMD) n = 19. Lacunar morphometric parameters were measured in both trabecular and cortical bone such as lacunar density (Lc.N/BV), lacunar volume (Lc.V), and lacunar sphericity (Lc.Sr). These were then compared against each other and also with previously measured tissue morphometries such as bone volume density (BV/TV), trabecular separation (Tb.Sp), trabecular number (Tb.N), and others. We detected no differences in lacunar morphology between the IOP, ILBMD and healthy premenopausal women. In contrast, we did find significant differences between lacunar morphologies including Lc.N/BV, Lc. V, and Lc. Sr in cortical and trabecular regions within all three groups (p < 0.001), which was consistent with our previous findings on a subgroup of the healthy group. Furthermore, we discovered strong correlations between Lc. Sr from trabecular regions with the measured BV/TV (R = -0.90, p < 0.05). The findings and comprehensive lacunar dataset we present here will be a crucial foundation for future investigations of the relationship between osteocyte lacunar morphology and disease.

A Comparative Study on the Multiscale Mechanical Responses of Human Femoral Neck Between the Young and the Elderly Using Finite Element Method

Background: Femoral neck fracture (FNF) is the most serious bone disease in the elderly population. The multiscale mechanical response is a key to predicting the strength of the femoral neck, assessing the risk of FNF, and exploring the role of mechanosensation and mechanotransmission in bone remodeling, especially in the context of aging bone.

Methods: Multiscale finite element (FE) models of the proximal femur for both young and elderly people were developed. The models included organ scale (proximal femur), tissue scale (cortical bone), tissue element scale (osteon), and cell scale [osteocyte lacuna-canalicular network (LCN) and extracellular matrix (ECM), OLCEM]. The mechanical responses of cortical bone and osteocytes in the mid-femoral neck and the differences in mechanical responses between these two scales were investigated.

Results: The mechanical responses of cortical bone and osteocyte showed significant differences between the elderly and the young. The minimum principal strains and mean SEDs of cortical bone in the elderly were 2.067–4.708 times and 3.093–14.385 times of the values in the young, respectively; the minimum principal strains and mean SEDs of osteocyte in the elderly were 1.497–3.246 times and 3.044–12 times of the values in the young, respectively; the amplification factors of minimum principal strain in the inferior (Inf), anterior (Ant), and posterior (Post) quadrants in the young were 1.241–1.804 times of the values in the elderly, but the amplification factor of minimum principal strain in the superior (Sup) quadrant was 87.4% of the value in the elderly; the amplification factors of mean SED in the young were 1.124–9.637 times of the values in the elderly.

Conclusion: The mass and bone mineral density (BMD) of cortical bone in the femoral neck is closely related to the mechanical response of osteocytes, which provides a new idea for improving cortical bone quality. Perhaps cortical bone quality could be improved by stimulating osteocytes. Quadrantal differences of bone quality in the mid-femoral neck should be considered to improve fracture risk prediction in the future.

Finite Element Models of Osteocytes and Their Load-Induced Activation

PURPOSE OF REVIEW

Osteocytes are the conductors of bone adaptation and remodelling. Buried inside the calcified matrix, they sense mechanical cues and signal osteoclasts in case of low activity, and osteoblasts when stresses are high. How do osteocytes detect mechanical stress? What physical signal do they perceive? Finite element analysis is a useful tool to address these questions as it allows calculating stresses, strains and fluid flow where they cannot be measured. The purpose of this review is to evaluate the capabilities and challenges of finite element models of bone, in particular the osteocytes and load-induced activation mechanisms.

RECENT FINDINGS

High-resolution imaging and increased computational power allow ever more detailed modelling of osteocytes, either in isolation or embedded within the mineralised matrix. Over the years, homogeneous models of bone and osteocytes got replaced by heterogeneous and microstructural models, including, e.g. the lacuno-canalicular network and the cytoskeleton. The lacuno-canalicular network induces strain amplifications and the osteocyte protrusions seem to be stimulated much more than the cell body, both by strain and fluid flow. More realistic cell geometries, like minute constrictions of the canaliculi, increase this effect. Microstructural osteocyte models describe the transduction of external stimuli to the nucleus. Supracellular multiscale models (e.g. of a tunnelling osteon) allow to study differential loading of osteocytes and to distinguish between strain and fluid flow as the pivotal stimulatory cue. In the future, the finite element models may be enhanced by including chemical transport and intercellular communication between osteocytes, osteoclasts and osteoblasts.

Perilacunar bone tissue exhibits sub-micrometer modulus gradation which depends on the recency of osteocyte bone formation in both young adult and early-old-age female C57Bl/6 mice

Osteocytes resorb and replace bone local to the lacunar-canalicular system (LCS). However, whether osteocyte remodeling impacts bone quality adjacent to the LCS is not understood. Further, while aging is well-established to decrease osteocyte viability and truncate LCS geometry, it is unclear if aging also decreases perilacunar bone quality. To address these questions, we employed atomic force microscopy (AFM) to generate nanoscale-resolution modulus maps for cortical femur osteocyte lacunae from young (5-month) and early-old-age (22-month) female C57Bl/6 mice. AFM-mapped lacunae were also imaged with confocal laser scanning microscopy to determine which osteocytes recently deposited bone as determined by the presence of fluorochrome labels administered 2d and 8d before euthanasia. Modulus gradation with distance from the lacunar wall was compared for labeled (i.e., bone forming) and non-labeled lacunae in both young and aged mice. All mapped lacunae showed sub-microscale modulus gradation, with peak modulus values 200-400 nm from the lacunar wall. Perilacunar modulus gradations depended on the recency of osteocyte bone formation (i.e., the presence of labels). For both ages, 2d-labeled perilacunar bone had lower peak and bulk modulus compared to non-labeled perilacunar bone. Lacunar length reduced with age, but lacunar shape and size were not strong predictors of modulus gradation. Our findings demonstrate for the first time that osteocyte perilacunar remodeling impacts bone tissue modulus, one contributor to bone quality. Given the immense scale of the LCS, differences in perilacunar modulus resulting from osteocyte remodeling activity may affect the quality of a substantial amount of bone tissue.

Impaired bone quality in the superolateral femoral neck occurs independent of hip geometry and bone mineral density

Skeletal adaptation is substantially influenced by mechanical loads. Osteocytes and their lacuno-canalicular network have been identified as a key player in load sensation and bone quality regulation. In the femoral neck, one of the most common fracture sites, a complex loading pattern with lower habitual loading in the superolateral neck and higher compressive stresses in the inferomedial neck is present. Variations in the femoral neck-shaft angle (NSA), i.e., coxa vara or coxa valga, provide the opportunity to examine the influence of loading patterns on bone quality. We obtained femoral neck specimens of 28 osteoarthritic human subjects with coxa vara, coxa norma and coxa valga during total hip arthroplasty. Bone mineral density (BMD) was assessed preoperatively by dual energy X-ray absorptiometry (DXA). Cortical and trabecular microstructure and three-dimensional osteocyte lacunar characteristics were assessed in the superolateral and inferomedial neck using ex vivo high resolution micro-computed tomography. Additionally, BMD distribution and osteocyte lacunar characteristics were analyzed by quantitative backscattered electron imaging (qBEI). All groups presented thicker inferomedial than superolateral cortices. Furthermore, the superolateral site exhibited a lower osteocyte lacunar density along with lower lacunar sphericity than the inferomedial site, independent of NSA. Importantly, BMD and corresponding T-scores correlated with microstructural parameters at the inferomedial but not superolateral neck. In conclusion, we provide micromorphological evidence for fracture vulnerability of the superolateral neck, which is independent of NSA and BMD. The presented bone qualitative data provide an explanation why DXA may be insufficient to predict a substantial proportion of femoral neck fractures. STATEMENT OF SIGNIFICANCE: The femoral neck, one of the most common fracture sites, is subject to a complex loading pattern. Site-specific differences (i.e., superolateral vs. inferomedial) in bone quality influence fracture risk, but it is unclear how this relates to hip geometry and bone mineral density (BMD) measurements in vivo. Here, we examine femoral neck specimens using a variety of high-resolution imaging techniques and demonstrate impaired bone quality in the superolateral compared to the inferomedial neck. Specifically, we found impaired cortical and trabecular microarchitecture, mineralization, and osteocyte properties, regardless of neck-shaft angle. Since BMD correlated with bone quality of the inferomedial but not the superolateral neck, our results illustrate why bone densitometry may not predict a substantial proportion of femoral neck fractures.

The Mechanosensory Role of Osteocytes and Implications for Bone Health and Disease States

Bone homeostasis is a dynamic equilibrium between bone-forming osteoblasts and bone-resorbing osteoclasts. This process is primarily controlled by the most abundant and mechanosensitive bone cells, osteocytes, that reside individually, within chambers of porous hydroxyapatite bone matrix. Recent studies have unveiled additional functional roles for osteocytes in directly contributing to local matrix regulation as well as systemic roles through endocrine functions by communicating with distant organs such as the kidney. Osteocyte function is governed largely by both biochemical signaling and the mechanical stimuli exerted on bone. Mechanical stimulation is required to maintain bone health whilst aging and reduced level of loading are known to result in bone loss. To date, both in vivo and in vitro approaches have been established to answer important questions such as the effect of mechanical stimuli, the mechanosensors involved, and the mechanosensitive signaling pathways in osteocytes. However, our understanding of osteocyte mechanotransduction has been limited due to the technical challenges of working with these cells since they are individually embedded within the hard hydroxyapatite bone matrix. This review highlights the current knowledge of the osteocyte functional role in maintaining bone health and the key regulatory pathways of these mechanosensitive cells. Finally, we elaborate on the current therapeutic opportunities offered by existing treatments and the potential for targeting osteocyte-directed signaling.

Structural role of osteocyte lacunae on mechanical properties of bone matrix: A cohesive finite element study

Despite the extensive studies on biological function of osteocytes, there are limited studies that evaluated the structural role of osteocyte lacunae on local mechanical properties of the bone matrix. As a result, the goal of this study was to elucidate the independent contribution of osteocyte lacunae structure on mechanical properties and fracture behavior of the bone matrix uncoupled from its biological effects and bone tissue composition variation. This study combined cohesive finite element modeling with experimental data from a lactation rat model to evaluate the influence of osteocyte lacunar area porosity, density, size, axis ratio, and orientation on the elastic modulus, ultimate strength, and ultimate strain of the bone matrix as well as on local crack formation and propagation. It also performed a parametric study to isolate the influence of a single osteocyte lacunae structural property on the mechanical properties of the bone matrix. The experimental measurements demonstrated statistically significant differences in lacunar size between ovariectomized rats with lactation history and virgin groups (both ovariectomized and intact) and in axis ratio between rats with lactation history and virgins. There were no differences in mechanical properties between virgin and lactation groups as determined by the finite element simulations. However, there were statistically significant linear relationships between the physiological range of osteocyte lacunar area porosity, density, size, and orientation and the elastic modulus and ultimate strength of the bone matrix in virgin and lactation rats. The parametric study also revealed similar but stronger relationships between elastic modulus and ultimate strength and lacunar density, size, and orientation. The simulations also demonstrated that the osteocyte lacunae guided the crack propagation through local stress concentrations. In summary, this study enhanced the limited knowledge on the structural role of osteocyte lacunae on local mechanical properties of the bone matrix. These data are important in gaining a better understanding of the mechanical implications of the local modifications due to osteocytes in the bone matrix.

Computational study of the mechanical influence of lacunae and perilacunar zones in cortical bone microcracking

The mechanical behavior of cortical bone is influenced by microstructural components such as osteons, Haversian canals, and osteocyte lacunae that arise from biological remodeling processes. This study takes a computational approach to investigate the role of the perilacunar zones formed by the local remodeling processes of lacunar-dwelling osteocytes by utilizing phase-field finite element models based on histological imaging of human bone. The models simulated the microdamage accumulation that occurs in cortical bone under transverse compression in bone without lacunae, with lacunae, and with a perilacunar zone surrounding lacunae in order to investigate the role of these features. The results of the simulations found that while lacunae create stress concentration which initiate further damage, perilacunar regions can delay or prevent the emergence and growth of microcracks.

Induction of somatopause in adult mice compromises bone morphology and exacerbates bone loss during aging

Somatopause refers to the gradual declines in growth hormone (GH) and insulin‐like growth factor‐1 throughout aging. To define how induced somatopause affects skeletal integrity, we used an inducible GH receptor knockout (iGHRKO) mouse model. Somatopause, induced globally at 6 months of age, resulted in significantly more slender bones in both male and female iGHRKO mice. In males, induced somatopause was associated with progressive expansion of the marrow cavity leading to significant thinning of the cortices, which compromised bone strength. We report progressive declines in osteocyte lacunar number, and increases in lacunar volume, in iGHRKO males, and reductions in lacunar number accompanied by ~20% loss of overall canalicular connectivity in iGHRKO females by 30 months of age. Induced somatopause did not affect mineral/matrix ratio assessed by Raman microspectroscopy. We found significant increases in bone marrow adiposity and high levels of sclerostin, a negative regulator of bone formation in iGHRKO mice. Surprisingly, however, despite compromised bone morphology, osteocyte senescence was reduced in the iGHRKO mice. In this study, we avoided the confounded effects of constitutive deficiency in the GH/IGF‐1 axis on the skeleton during growth, and specifically dissected its effects on the aging skeleton. We show here, for the first time, that induced somatopause compromises bone morphology and the bone marrow environment.

Multiscale mechanical responses of young and elderly human femurs: A finite element investigation

BACKGROUND

Bone remodeling in the elderly is no longer balanced. As a result, the morphologies and mechanical properties of bone at different scales will change. These changes would affect the mechanical responses of bone, which might exacerbate the imbalance of bone remodeling and even cause age-related bone diseases.

METHODS

Considering those changes, multiscale finite element (FE) models of bone in the young and the elderly were developed that included macroscale (proximal femur), mesoscale (cortical bone), microscale (Haversian system) and sub-microscale (osteocyte-lacuna-canaliculus-extracellular matrix system, OLCES). The stress and strain distributions at different scales and transmissions among different scales were investigated.

RESULTS

The stresses of the elderly at macroscale, mesoscale and microscale were higher than those in the young by 23.7%, 62.5% and 8.0%, respectively, and the stresses of the elderly and the young at sub-microscale were almost the same. The strain of the elderly at macroscale, mesoscale, microscale and sub-microscale were higher than those in the young by 48.6%, 56.8%, 11.9% and 25.1%, respectively. The stress and strain transmission rates (ησand ηε) from mesoscale to microscale were decreased by 1.8%, and 2.5% than those from macroscale to mesoscale in the elderly, respectively; but increased by 13.8%, and 4.7% in the young, respectively. ηε from microscale to sub-microscale in the elderly was higher than that in the young by 21.3%.

CONCLUSIONS

Degeneration of cortical bone mechanical property in the elderly causes increases in stress and strain at macroscale and mesoscale. The reduction of lacunar number in the elderly is not conducive to the mechanical transmission from mesoscale to microscale. The differences in stress and strain at microscale between the young and the elderly are smaller than those at macroscale or mesoscale. The strain stimulus sensed by osteocyte in the elderly is not weakened compared with that in the young.

Tunnels in the rock: Dynamics of osteocyte morphogenesis

Osteocytes are dynamic, bone matrix-remodeling cells that form an intricate network of interconnected projections through the bone matrix, called the lacunar-canalicular system. Osteocytes are the dominant mechanosensory cells in bone and their mechanosensory and mechanotransductive functions follow their morphological form. During osteocytogenesis and development of the osteocyte lacunar-canalicular network, osteocytes must dramatically remodel both their cytoskeleton and their extracellular matrix. In this review, we summarize our current understanding of the mechanisms that govern osteocyte differentiation, cytoskeletal morphogenesis, mechanotransduction, and matrix remodeling. We postulate that the physiologic activation of matrix remodeling in adult osteocytes, known as perilacunar/canalicular remodeling (PLR) represents a re-activation of the developmental program by which the osteocyte network is first established. While much of osteocyte biology remains unclear, new tools and approaches make the present moment a particularly fruitful and exciting time to study the development of these remarkable cells.

Alterations in osteocyte lacunar morphology affect local bone tissue strains

Osteocytes are capable of remodeling their perilacunar bone matrix, which causes considerable variations in the shape and size of their lacunae. If these variations in lacunar morphology cause changes in the mechanical environment of the osteocytes, in particular local strains, they would subsequently affect bone mechanotransduction, since osteocytes are likely able to directly sense these strains. The purpose of this study is to quantify the effect of alterations in osteocyte lacunar morphology on peri-lacunar bone tissue strains. To this end, we related the actual lacunar shape in fibulae of six young-adult (5-month) and six old (23-month) mice, quantified by high-resolution micro-computed tomography, to microscopic strains, analyzed by micro-finite element modeling. We showed that peak effective strain increased by 12.6% in osteocyte cell bodies (OCYs), 9.6% in pericellular matrix (PCM), and 5.3% in extra cellular matrix (ECM) as the lacunae volume increased from 100-200 μm³ to 500-600 μm³. Lacunae with a larger deviation (>8°) in orientation from the longitudinal axis of the bone are exposed to 8% higher strains in OCYs, 6.5% in PCM, 4.2% in ECM than lacunae with a deviation in orientation below 8°. Moreover, increased lacuna sphericity from 0 to 0.5 to 0.7-1 led to 25%, 23%, and 13% decrease in maximum effective strains in OCYs, PCM, and ECM, respectively. We further showed that due to the presence of smaller and more round lacunae in old mice, local bone tissue strains are on average 5% lower in the vicinity of lacunae and their osteocytes of old mice compared to young. Understanding how changes in lacunar morphology affect the micromechanical environment of osteocytes presents a first step in unraveling their potential role in impaired bone mechanoresponsiveness with e.g. aging.

Cementocyte alterations associated with experimentally induced cellular cementum apposition in hyp mice

BACKGROUND

Cellular cementum, a mineralized tissue covering apical tooth roots, grows by apposition to maintain the tooth in its occlusal position. We hypothesized that resident cementocytes would show morphological changes in response to cementum apposition, possibly implicating a role in cementum biology.

METHODS

Mandibular first molars were induced to super-erupt (EIA) by extraction of maxillary molars, promoting rapid new cementum formation. Tissue and cell responses were analyzed at 6 and/or 21 days post-procedure (dpp).

RESULTS

High-resolution micro-computed tomography (micro-CT) and confocal laser scanning microscopy showed increased cellular cementum by 21 dpp. Transmission electron microscopy (TEM) revealed that cementocytes under EIA were 50% larger than control cells, supported by larger pore sizes detected by micro-CT. Cementocytes under EIA displayed ultrastructural changes consistent with increased activity, including increased cytoplasm and nuclear size. We applied EIA to Hyp mutant mice, where cementocytes have perilacunar hypomineralization defects, to test cell and tissue responses in an altered mechanoresponsive milieu. Hyp and WT molars displayed similar super-eruption, with Hyp molars exhibiting 28% increased cellular cementum area versus 22% in WT mice at 21 dpp. Compared to control, Hyp cementocytes featured well-defined, disperse euchromatin and a thick layer of peripherally condensed heterochromatin in nuclei, indicating cellular activity. Immunohistochemistry (IHC) for cementum markers revealed intense dentin matrix protein-1 expression and abnormal osteopontin deposition in Hyp mice. Both WT and Hyp cementocytes expressed gap junction protein, connexin 43.

CONCLUSION

This study provides new insights into the EIA model and cementocyte activity in association with new cementum formation.

Bone physiological microenvironment and healing mechanism: Basis for future bone-tissue engineering scaffolds

Bone-tissue defects affect millions of people worldwide. Despite being common treatment approaches, autologous and allogeneic bone grafting have not achieved the ideal therapeutic effect. This has prompted researchers to explore novel bone-regeneration methods. In recent decades, the development of bone tissue engineering (BTE) scaffolds has been leading the forefront of this field. As researchers have provided deep insights into bone physiology and the bone-healing mechanism, various biomimicking and bioinspired BTE scaffolds have been reported. Now it is necessary to review the progress of natural bone physiology and bone healing mechanism, which will provide more valuable enlightenments for researchers in this field. This work details the physiological microenvironment of the natural bone tissue, bone-healing process, and various biomolecules involved therein. Next, according to the bone physiological microenvironment and the delivery of bioactive factors based on the bone-healing mechanism, it elaborates the biomimetic design of a scaffold, highlighting the designing of BTE scaffolds according to bone biology and providing the rationale for designing next-generation BTE scaffolds that conform to natural bone healing and regeneration.

Osteon: Structure, Turnover, and Regeneration

Bone is composed of dense and solid cortical bone and honeycomb-like trabecular bone. Although cortical bone provides the majority of mechanical strength for a bone, there are few studies focusing on cortical bone repair or regeneration. Osteons (the Haversian system) form structural and functional units of cortical bone. In recent years, emerging evidences have shown that the osteon structure (including osteocytes, lamellae, lacunocanalicular network, and Haversian canals) plays critical roles in bone mechanics and turnover. Therefore, reconstruction of the osteon structure is crucial for cortical bone regeneration. This article provides a systematic summary of recent advances in osteons, including the structure, function, turnover, and regenerative strategies. First, the hierarchical structure of osteons is illustrated and the critical functions of osteons in bone dynamics are introduced. Next, the modeling and remodeling processes of osteons at a cellular level and the turnover of osteons in response to mechanical loading and aging are emphasized. Furthermore, several bioengineering approaches that were recently developed to recapitulate the osteon structure are highlighted. Impact statement This review provides a comprehensive summary of recent advances in osteons, especially the roles in bone formation, remodeling, and regeneration. Besides introducing the hierarchical structure and critical functions of osteons, we elucidate the modeling and remodeling of osteons at a cellular level. Specifically, we highlight the bioengineering approaches that were recently developed to mimic the hierarchical structure of osteons. We expect that this review will provide informative insights and attract increasing attentions in orthopedic community, shedding light on cortical bone regeneration in the future.

Lacunar-canalicular bone remodeling: Impacts on bone quality and tools for assessment

Osteocytes can resorb as well as replace bone adjacent to the expansive lacunar-canalicular system (LCS). Suppressed LCS remodeling decreases bone fracture toughness, but it is unclear how altered LCS remodeling impacts bone quality. The first goal of this review is to assess how LCS remodeling impacts LCS morphology as well as the composition and mechanical properties of surrounding bone tissue. The second goal is to compare tools available for the assessment of bone quality at length-scales that are physiologically-relevant to LCS remodeling. We find that changes to LCS morphology occur in response to a variety of physiological conditions and diseases and can be classified in two general phenotypes. In the 'aging phenotype', seen in aging and in some disuse models, the LCS is truncated and osteocytes apoptosis is increased. In the 'osteocytic osteolysis' phenotype, which is adaptive in some physiological settings and possibly maladaptive in others, the LCS enlarges and osteocytes generally maintain viability. Bone composition and mechanical properties vary near the osteocyte and change with at least some conditions that alter LCS morphology. However, few studies have evaluated bone composition and mechanical properties close to the LCS and so the impacts of LCS remodeling phenotypes on bone tissue quality are still undetermined. We summarize the current understanding of how LCS remodeling impacts LCS morphology, tissue-scale bone composition and mechanical properties, and whole-bone material properties. Tools are compared for assessing tissue-scale bone properties, as well as the resolution, advantages, and limitations of these techniques.

The Osteocyte: New Insights

Osteocytes are an ancient cell, appearing in fossilized skeletal remains of early fish and dinosaurs. Despite its relative high abundance, even in the context of nonskeletal cells, the osteocyte is perhaps among the least studied cells in all of vertebrate biology. Osteocytes are cells embedded in bone, able to modify their surrounding extracellular matrix via specialized molecular remodeling mechanisms that are independent of the bone forming osteoblasts and bone-resorbing osteoclasts. Osteocytes communicate with osteoclasts and osteoblasts via distinct signaling molecules that include the RankL/OPG axis and the Sost/Dkk1/Wnt axis, among others. Osteocytes also extend their influence beyond the local bone environment by functioning as an endocrine cell that controls phosphate reabsorption in the kidney, insulin secretion in the pancreas, and skeletal muscle function. These cells are also finely tuned sensors of mechanical stimulation to coordinate with effector cells to adjust bone mass, size, and shape to conform to mechanical demands.

A 3D, Dynamically Loaded Hydrogel Model of the Osteochondral Unit to Study Osteocyte Mechanobiology

Osteocytes are mechanosensitive cells that orchestrate signaling in bone and cartilage across the osteochondral unit. The mechanisms by which osteocytes regulate osteochondral homeostasis and degeneration in response to mechanical cues remain unclear. This study introduces a novel 3D hydrogel bilayer composite designed to support osteocyte differentiation and bone matrix deposition in a bone-like layer and to recapitulate key aspects of the osteochondral unit's complex loading environment. The bilayer hydrogel is fabricated with a soft cartilage-like layer overlaying a stiff bone-like layer. The bone-like layer contains a stiff 3D-printed hydrogel structure infilled with a soft, degradable, cellular hydrogel. The IDG-SW3 cells embedded within the soft hydrogel mature into osteocytes and produce a mineralized collagen matrix. Under dynamic compressive strains, near-physiological levels of strain are achieved in the bone layer (≤ 0.08%), while the cartilage layer bears the majority of the strains (>99%). Under loading, the model induces an osteocyte response, measured by prostaglandin E2, that is frequency, but not strain, dependent: a finding attributed to altered fluid flow within the composite. Overall, this new hydrogel platform provides a novel approach to study osteocyte mechanobiology in vitro in an osteochondral tissue-mimetic environment.

Large osteocyte lacunae in iliac crest infantile bone are not associated with impaired mineral distribution or signs of osteocytic osteolysis

The enlargement of osteocyte lacunae via osteocytic osteolysis was previously detected in situations of increased calcium demand (e.g., lactation, vitamin D deficiency). However, it is unclear whether similar processes occur also in the growing infantile skeleton and how this is linked to the mineral distribution within the bone matrix. Human iliac crest biopsies of 30 subjects (0-6 months, n = 14; 2-8 years, n = 6 and 18-25 years, n = 10) were acquired. Bone microarchitecture was assessed by micro-CT, while cellular bone histomorphometry was performed on undecalcified histological sections. Quantitative backscattered electron imaging (qBEI) was conducted to determine the bone mineral density distribution (BMDD) as well as osteocyte lacunar size and density. We additionally evaluated cathepsin K positive osteocytes using immunohistochemistry. Infantile bone was characterized by various signs of ongoing bone development such as higher bone (re)modeling, lower cortical and trabecular thickness compared to young adults. Importantly, a significantly higher osteocyte lacunar density and increased lacunar area were detected. Large osteocyte lacunae were associated with a more heterogeneous bone mineral density distribution of the trabecular bone matrix due to the presence of hypermineralized cartilage remnants, whereas the mean mineralization (i.e., CaMean) was not different in infantile bone. Absence of cathepsin K expression in osteocyte lacunae indicated nonexistent osteocytic osteolysis. Taken together, we demonstrated that the overall mineralization distribution in infantile bone is not altered compared to young adults besides high trabecular mineralization heterogeneity. Our study also provides important reference values for bone microstructure, BMDD and osteocyte characteristics in infants, children and young adults. Infantile bone displays large osteocyte lacunae indicating a developmental phenomenon rather than osteocytic osteolysis. Larger osteocytes may have superior mechanosensory abilities to enable bone adaption during growth.

Osteocyte lacunar strain determination using multiscale finite element analysis

Osteocytes are thought to be the primary mechanosensory cells within bone, regulating both osteoclasts and osteoblasts to control load induced changes in bone resorption and formation. Osteocytes initiate intracellular responses including activating the Wnt/β-catenin signaling pathway after experiencing mechanical forces. In response to changing mechanical loads (strain) the osteocytes signal to cells on the bone surface. However, this process of osteocyte activation appears heterogeneous since it occurs in sub-populations of osteocytes, even within regions predicted to be experiencing similar global strain magnitudes determined based on traditional finite element modeling approaches. Several studies have investigated the strain responses of osteocyte lacunae using finite element (FE) models, but many were limited by the use of idealized geometries (e.g., ellipsoids) and analysis of a single osteocyte. Finite element models by other groups included more details, such as canaliculi, but all were done on models consisting of a single osteocyte. We hypothesized that variation in size and orientation of the osteocyte lacunae within bone would give rise to micro heterogeneity in the strain fields that could better explain the observed patterns of osteocyte activation following load. The osteocytes in our microscale and nanoscale models have an idealized oval shape and some are based on confocal scans. However, all the FE models in this preliminary study consist of multiple osteocytes. The number of osteocytes in the 3D confocal scan models ranged from five to seventeen. In this study, a multi-scale computational approach was used to first create an osteocyte FE model at the microscale level to examine both the theoretical lacunar and perilacunar strain responses based on two parameters: 1) lacunar orientation and 2) lacunar size. A parametric analysis was performed by steadily increasing the perilacunar modulus (5, 10, 15, and 20 GPa). Secondly, a nanoscale FE model was built using known osteocyte dimensions to determine the predicted strains in the perilacunar matrix, fluid space, and cell body regions. Finally, 3-D lacunar models were created using confocal image stacks from mouse femurs to determine the theoretical strain in the lacunae represented by realistic geometries. Overall, lacunar strains decreased by 14% in the cell body, 15% in the fluid space region and 25% in the perilacunar space as the perilacunar modulus increased, indicating a stress shielding effect. Lacunar strains were lower for the osteocytes aligned along the loading axis compared to those aligned perpendicular to axis. Increases in lacuna size also led to increased lacunar strains. These finite element model findings suggest that orientation and lacunar size may contribute to the heterogeneous initial pattern of osteocyte strain response observed in bone following in vivo applied mechanical loads. A better understanding of how mechanical stimuli directly affect the lacunae and perilacunar tissue strains may ultimately lead to a better understanding of the process of osteocyte activation in response to mechanical loading.

•

A multi-scale computational approach used to first create multiple osteocyte FE model at the microscale level

•

3-D Lacuna model created using confocal image stacks from a mouse femur to determine the theoretical strain in the lacunae.

•

Lacunar strains decreased as the perilacunar modulus increased.

•

Lacunar strains were lower for the osteocytes aligned along the loading axis compared to those aligned perpendicular to axis.

•

Increases in lacuna size also led to increased lacunar strains

Sequential Treatment of Estrogen Deficient, Osteopenic Rats with Alendronate, Parathyroid Hormone (1-34), or Raloxifene Alters Cortical Bone Mineral and Matrix Composition

Anti-resorptive and anabolic treatments can be used sequentially to treat osteoporosis, but their effects on bone composition are incompletely understood. Osteocytes may influence bone tissue composition with sequential therapies because bisphosphonates diffuse into the canalicular network and anabolic treatments increase osteocyte lacunar size. Cortical bone composition of osteopenic, ovariectomized (OVX) rats was compared to that of Sham-operated rats and OVX rats given monotherapy or sequential regimens of single approved anti-osteoporosis medications. Adult female Sprague-Dawley rats were OVX (N = 37) or Sham-OVXd (N = 6). After 2 months, seven groups of OVX rats were given three consecutive 3-month periods of treatment with vehicle (V), h-PTH (1-34) (P), alendronate (A), or raloxifene (R), using the following orders: VVV, PVV, RRR, RPR, AAA, AVA, and APA. Compositional properties around osteocyte lacunae of the left tibial cortex were assessed from Raman spectra in perilacunar and non-perilacunar bone matrix regions. Sequential treatments involving parathyroid hormone (PTH) caused lower mean collagen maturity relative to monotherapies. Mean mineral:matrix ratio was 2.2% greater, mean collagen maturity was 1.4% greater, and mean carbonate:phosphate ratio was 2.2% lower in the perilacunar than in the non-perilacunar bone matrix region (all P < 0.05). These data demonstrate cortical bone tissue composition differences around osteocytes caused by sequential treatment with anti-osteoporosis medications. We speculate that the region-specific differences demonstrate the ability of osteocytes to alter bone tissue composition adjacent to lacunae.

What is the influence of two strain rates on the relationship between human cortical bone toughness and micro-structure?

Cortical bone fracture mechanisms are well studied under quasi-static loading. The influence of strain rate on crack propagation mechanisms needs to be better understood, however. We have previously shown that several aspects of the bone micro-structure are involved in crack propagation, such as the complete porosity network, including the Haversian system and the lacunar network, as well as biochemical aspects, such as the maturity of collagen cross-links. The aim of this study is to investigate the influence of strain rate on the toughness of human cortical bone with respect to its microstructure and organic non-collagenous composition. Two strain rates will be considered: quasi-static loading (10−4 s−1), a standard condition, and a higher loading rate (10−1 s−1), representative of a fall. Cortical bone samples were extracted from eight female donors (age 50–91 years). Three-point bending tests were performed until failure. Synchrotron radiation micro-computed tomography imaging was performed to assess bone microstructure including the Haversian system and the lacunar system. Collagen enzymatic cross-link maturation was measured using a high performance liquid chromatography column. Results showed that that under quasi-static loading, the elastic contribution of the fracture process is correlated to both the collagen cross-links maturation and the microstructure, while the plastic contribution is correlated only to the porosity network. Under fall-like loading, bone organization appears to be less linked to crack propagation.

Development of protocols for the first serial block-face scanning electron microscopy (SBF SEM) studies of bone tissue

There is an unmet need for a high-resolution three-dimensional (3D) technique to simultaneously image osteocytes and the matrix in which these cells reside. In serial block-face scanning electron microscopy (SBF SEM), an ultramicrotome mounted within the vacuum chamber of a microscope repeatedly sections a resin-embedded block of tissue. Backscattered electron scans of the block face provide a stack of high-resolution two-dimensional images, which can be used to visualise and quantify cells and organelles in 3D. High-resolution 3D images of biological tissues from SBF SEM have been exploited considerably to date in the neuroscience field. However, non-brain samples, in particular hard biological tissues, have appeared more challenging to image by SBF SEM due to the difficulties of sectioning and rendering the samples conductive. We have developed and propose protocols for bone tissue preparation using SBF SEM, for imaging simultaneously soft and hard bone tissue components in 3D. We review the state of the art in high-resolution imaging of osteocytes, provide a historical perspective of SBF SEM, and we present first SBF SEM proof-of-concept studies for murine and human tissue. The application of SBF SEM to hard tissues will facilitate qualitative and quantitative 3D studies of tissue microstructure and ultrastructure in bone development, ageing and pathologies such as osteoporosis and osteoarthritis.

Effects of occlusal forces on the peri-implant-bone interface stability

The occlusal forces and their influence on the initiation of peri-implant bone loss or their relationship with peri-implantitis have created discussion during the past 30 years given the discrepancies observed in clinical, animal, and finite element analysis studies. Beyond these contradictions, in the case of an osseointegrated implant, the occlusal forces can influence the implant-bone interface and the cells responsible for the bone remodeling in different ways that may result in the maintenance or loss of the osseointegration. This comprehensive review focuses on the information available about the forces transmitted through the implant-crown system to the implant-bone interface and the mechano-transduction phenomena responsible for the bone cells' behavior and their interactions. Knowledge of the basic molecular biology of the peri-implant bone would help clinicians to understand the complex phenomenon of occlusal forces and their effects on the implant-bone interface, and would allow better control of the negative effects of mechanical stresses, leading to therapy with fewer risks and complications.

The effects of mechanically loaded osteocytes and inflammation on bone remodeling in a bisphosphonate-induced environment

Bisphosphonate-related osteonecrosis of the jaw is a disease appearing after tooth removal in patients undergoing bisphosphonate treatment for metastasizing cancers and osteoporosis. The complexity of the condition requires a multicellular model to address the net effects of two key risk factors: mechanical trauma (pathologic overload) and inflammation. In this work, a system comprised of a polydimethylsiloxane chip and mechanical loading device is used to expose bisphosphonate-treated osteocytes to mechanical trauma. Specifically, osteocytes are treated with the potent nitrogen-containing bisphosphonate, zoledronic acid, and exposed to short-term pathologic overload via substrate stretch. During bone remodeling, osteocyte apoptosis plays a role in attracting pre-osteoclasts to sites of damage; as such, lactate dehydrogenase activity, cell death and protein expression are evaluated as functions of load. Additionally, the effects of osteocyte soluble factors on osteoclast and osteoblast functional activity are quantified. Osteoclast activity and bone resorption are quantified in the presence and absence of inflammatory components, lipopolysaccharide and interferon gamma. Results suggest that inflammation associated with bacterial infection may hinder bone resorption by osteoclasts. In addition, osteocytes may respond to overload by altering expression of soluble signals that act on osteoblasts to attenuate bone formation. These findings give insight into the multicellular interactions implicated in bisphosphonate-related osteonecrosis of the jaw.

Examining the influence of PTH(1-34) on tissue strength and composition

The lacunar-canaliculi system is a network of channels that is created and maintained by osteocytes as they are embedded throughout cortical bone. As osteocytes modify their lacuna space, the local tissue composition and tissue strength are subject to change. Although continual exposure to parathyroid hormone (PTH) can induce adaptation at the lacunar wall, the impact of intermittent PTH treatment on perilacunar adaptation remains unclear. Therefore, the primary objective of this study was to establish how intermittent PTH(1-34) treatment influences perilacunar adaptation with respect to changes in tissue composition. We hypothesized that local changes in tissue composition following PTH(1-34) are associated with corresponding gains in tissue strength and resistance to microdamage at the whole bone level. Adult male C57BL/6J mice were treated daily with PTH(1-34) or vehicle for 3 weeks. In response to PTH(1-34), Raman spectroscopy revealed a significant decrease in the carbonate-to-phosphate ratio and crystallinity across the entire tissue, while the mineral-to-matrix ratio demonstrated a significant decrease in just the perilacunar region. The shift in perilacunar composition largely explained the corresponding increase in tissue strength, while the degree of new tissue added at the endosteum and periosteum did not produce any significant changes in cortical area or moment of inertia that would explain the increase in tissue strength. Furthermore, fatigue testing revealed a greater resistance to crack formation within the existing tissue following PTH(1-34) treatment. As a result, the shift in perilacunar composition presents a unique mechanism by which PTH(1-34) produces localized differences in tissue quality that allow more energy to be dissipated under loading, thereby increasing tissue strength and resistance to microdamage. In addition, our findings demonstrate the potential for PTH(1-34) to amplify osteocytes' mechanotransduction by producing a more compliant tissue. Overall, the present study demonstrates that changes in tissue composition localized at the lacuna wall contribute to the strength and fatigue resistance of cortical bone gained in response to intermittent PTH(1-34) treatment.

Mechanical Stimuli in the Local In Vivo Environment in Bone: Computational Approaches Linking Organ-Scale Loads to Cellular Signals

PURPOSE OF REVIEW

Connecting organ-scale loads to cellular signals in their local in vivo environment is a current challenge in the field of bone (re)modelling. Understanding this critical missing link would greatly improve our ability to anticipate mechanotransduction during different modes of stimuli and the resultant cellular responses. This review characterises computational approaches that could enable coupling links across the multiple scales of bone.

RECENT FINDINGS

Current approaches using strain and fluid shear stress concepts have begun to link organ-scale loads to cellular signals; however, these approaches fail to capture localised micro-structural heterogeneities. Furthermore, models that incorporate downstream communication from osteocytes to osteoclasts, bone-lining cells and osteoblasts, will help improve the understanding of (re)modelling activities. Incorporating this potentially key information in the local in vivo environment will aid in developing multiscale models of mechanotransduction that can predict or help describe resultant biological events related to bone (re)modelling. Progress towards multiscale determination of the cell mechanical environment from organ-scale loads remains elusive. Construction of organ-, tissue- and cell-scale computational models that include localised environmental variation, strain amplification and intercellular communication mechanisms will ultimately help couple the hierarchal levels of bone.

Age-related changes in female mouse cortical bone microporosity

Osteocyte lacunae are small cavities within the bone matrix. Their dimensions and spatial arrangement affect bone mechanical properties. Furthermore, their size and shape affect the strain in bone tissue close to the lacunae; hence, they are supposed to affect the mechanosensory function of the osteocytes residing in the lacunae. It was the purpose of this study to quantify the morphological features of osteocyte lacunae, whether these are affected by aging and whether these vary among different anatomical location. In addition, we aimed at quantifying the vascular canals as these affect bone's microporosity too. We quantified the microporosity in the fibular midshaft of young-adult and old female C57BL/6 mice. Using micro-computed tomography (μCT), we found that advanced age was associated with a significantly decreased vascular canal number and volume density. In aged mice, the mean volume of the lacuna was significantly smaller than in young animals and they were more round. Lacuna number density close to the neutral axis of the fibula was higher in older mice than in young ones. The characterization of bone microporosity presents a first step in further unraveling their potential role in age-related reductions in bone strength.

Effect of osteoporosis treatment agents on the cortical bone osteocyte microenvironment in adult estrogen-deficient, osteopenic rats

Though osteoporosis is a significant cause of disability worldwide, treatment with pharmacologic agents decreases risk of fragility fracture. Though these treatments act through the bone remodeling system to improve bone mass, it is unclear if they alter the response of bone to mechanical loading at the level of the osteocyte. This pre-clinical study determined the relationship between microstructural bone tissue properties and osteocyte lacunar size and density to strain around osteocytes with standard osteoporosis treatment or sequential therapies. Six-month-old female ovariectomized (OVX) Sprague-Dawley rats were cycled through various sequences of pharmacological treatments [alendronate (Aln), raloxifene (Ral) and human parathyroid hormone-1,34 (PTH)] for three month intervals, over nine months. Linear nanoindentation mapping was used to determine Young's modulus in perilacunar and bone matrix regions around cortical bone osteocyte lacunae. Measurements of lacunar diameter and density were completed. Treatment-related differences in Young's modulus in the perilacunar and bone matrix regions were not observed. We confirmed previous data that showed that the bone matrix region was stiffer than the perilacunar matrix region. Whole bone material properties were correlated to perilacunar matrix stiffness. Finite element models predicted a range of mechanical strain amplification factors estimated at the osteocyte across treatment groups. In summary, though the perilacunar matrix near cortical osteocyte lacuna is not as stiff as bone matrix further away, osteoporosis treatment agents do not affect the stiffness of bone tissue near osteocyte lacunae.

•

Monotherapy with osteoporosis treatment agents does not affect the stiffness of bone tissue around osteocyte lacunae.

•

Sequential use of osteoporosis treatment agents does not affect bone tissue stiffness around osteocyte lacunae.

•

Perilacunar cortical bone tissue is not as stiff as bone matrix further from osteocyte lacunae.

•

Whole bone material properties are negatively correlated to the stiffness of perilacunar bone tissue.

A mathematical multiscale model of bone remodeling, accounting for pore space-specific mechanosensation

While bone tissue is a hierarchically organized material, mathematical formulations of bone remodeling are often defined on the level of a millimeter-sized representative volume element (RVE), "smeared" over all types of bone microstructures seen at lower observation scales. Thus, there is no explicit consideration of the fact that the biological cells and biochemical factors driving bone remodeling are actually located in differently sized pore spaces: active osteoblasts and osteoclasts can be found in the vascular pores, whereas the lacunar pores host osteocytes - bone cells originating from former osteoblasts which were then "buried" in newly deposited extracellular bone matrix. We here propose a mathematical description which considers size and shape of the pore spaces where the biological and biochemical events take place. In particular, a previously published systems biology formulation, accounting for biochemical regulatory mechanisms such as the rank-rankl-opg pathway, is cast into a multiscale framework coupled to a poromicromechanical model. The latter gives access to the vascular and lacunar pore pressures arising from macroscopic loading. Extensive experimental data on the biological consequences of this loading strongly suggest that the aforementioned pore pressures, together with the loading frequency, are essential drivers of bone remodeling. The novel approach presented here allows for satisfactory simulation of the evolution of bone tissue under various loading conditions, and for different species; including scenarios such as mechanical dis- and overuse of murine and human bone, or in osteocyte-free bone.

Computational Investigation on the Biomechanical Responses of the Osteocytes to the Compressive Stimulus: A Poroelastic Model

Osteocytes, the major type of bone cells embedded in the bone matrix and surrounded by the lacunar and canalicular system, can serve as biomechanosensors and biomechanotranducers of the bone. Theoretical analytical methods have been employed to investigate the biomechanical responses of osteocytes in vivo; the poroelastic properties have not been taken into consideration in the three-dimensional (3D) finite element model. In this study, a 3D poroelastic idealized finite element model was developed and was used to predict biomechanical behaviours (maximal principal strain, pore pressure, and fluid velocity) of the osteocyte-lacunar-canalicular system under 150-, 1000-, 3000-, and 5000-microstrain compressive loads, respectively, representing disuse, physiological, overuse, and pathological overload loading stimuli. The highest local strain, pore pressure, and fluid velocity were found to be highest at the proximal region of cell processes. These data suggest that the strain, pore pressure, and fluid velocity of the osteocyte-lacunar-canalicular system increase with the global loading and that the poroelastic material property affects the biomechanical responses to the compressive stimulus. This new model can be used to predict the mechanobiological behaviours of osteocytes under the four different compressive loadings and may provide an insight into the mechanisms of mechanosensation and mechanotransduction of the bone.