Network architecture strongly influences the fluid flow pattern through the lacunocanalicular network in human osteons

Influence of intramedullary pressure on Lacuno-Canalicular fluid flow: A systematic review

While most of current models investigating bone remodelling are based on matrix deformation, intramedullary pressure also plays a role. Bone remodelling is orchestrated by the Lacuno-Canalicular Network (LCN) fluid-flow. The aim of this review was hence to assess the influence of intramedullary pressure on the fluid circulation within the LCN. Three databases (Science Direct, Web of Science, and PubMed) were used. The first phase of the search returned 731 articles, of which 9 respected the inclusion/exclusion criteria and were included. These studies confirm the association between intramedullary pressure and fluid dynamics in the LCN. Among the included studies, 7 experimental studies using animal models and 2 numerical models were found. The studies were then ranked according to the nature of the applied loading, either axial compression or direct cyclic intramedullary pressure. The current review revealed that there is an influence of intramedullary pressure on LCN fluid dynamics and that this influence depends on the magnitude and the frequency of the applied pressure. Two studies confirmed that the influence was effective even without bone matrix deformation. While intramedullary pressure is closely associated with LCN fluid, there is a severe lack of studies on this topic. STATEMENT OF SIGNIFICANCE: Since the 1990's, numerical models developed to investigate fluid flow in bone submicrometric porous network are based on the flow induced by matrix deformation. Bone fluid flow is known to be involved in cells stimulation and hence directly influences bone remodeling. Different studies have shown that intramedullary pressure is also associated with bone mechanosensitive adaptation. This pressure is developed in bone due to blood circulation and is increased during loading or muscle stimulation. The current article reviews the studies investigating the influence of this pressure on bone porous fluid flow. They show that fluid flow is involved by this pressure even without bone matrix deformation. The current review article highlights the severe lack of studies about this mechanism.

From breast cancer cell homing to the onset of early bone metastasis: The role of bone (re)modeling in early lesion formation

Breast cancer often metastasizes to bone, causing osteolytic lesions. Structural and biophysical changes are rarely studied yet are hypothesized to influence metastasis. We developed a mouse model of early bone metastasis and multimodal imaging to quantify cancer cell homing, bone (re)modeling, and onset of metastasis. Using tissue clearing and three-dimensional (3D) light sheet fluorescence microscopy, we located enhanced green fluorescent protein-positive cancer cells and small clusters in intact bones and quantified their size and spatial distribution. We detected early bone lesions using in vivo microcomputed tomography (microCT)-based time-lapse morphometry and revealed altered bone (re)modeling in the absence of detectable lesions. With a new microCT image analysis tool, we tracked the growth of early lesions over time. We showed that cancer cells home in all bone compartments, while osteolytic lesions are only detected in the metaphysis, a region of high (re)modeling. Our study suggests that higher rates of (re)modeling act as a driver of lesion formation during early metastasis.

Peripheral canalicular branching is decreased in streptozotocin-induced diabetes and correlates with decreased whole-bone ultimate load and perilacunar elastic work

Osteocytes, the most abundant cell type in bone, play a crucial role in mechanosensation and signaling for bone formation and resorption. These cells reside within a complex lacuno-canalicular network (OLCN). Osteocyte signaling is reduced under diabetic conditions, and both type 1 and type 2 diabetes lead to reduced bone turnover, perturbed bone composition, and increased fracture risk. We hypothesized that this reduced bone turnover, and altered bone composition with diabetes is associated with reduced OLCN architecture and connectivity. This study aimed to elucidate: (1) the sequence of OLCN changes with diabetes related to bone turnover and (2) whether changes to the OLCN are associated with tissue composition and mechanical properties. Twelve- to fourteen-week-old male C57BL/6 mice were administered streptozotocin at 50 mg/kg for 5 consecutive days to induce hyperglycemia, sacrificed at baseline (BL), or after being diabetic for 3 (D3) and 7 (D7) wk with age-matched (C3, C7) controls (n = 10-12 per group). Mineralized femoral sections were infiltrated with rhodamine, imaged with confocal microscopy, then the OLCN morphology and topology were characterized and correlated against bone histomorphometry, as well as local and whole-bone mechanics and composition. D7 mice exhibited a lower number of peripheral branches relative to C7. The total number of canalicular intersections (nodes) was lower in D3 and D7 relative to BL (P < 0.05 for all), and a reduced bone formation rate (BFR) was observed at D7 vs C7. The number of nodes explained only 15% of BFR, but 45% of Ct.BV/TV, and 31% of ultimate load. The number of branches explained 30% and 22% of the elastic work at the perilacunar and intracortical region, respectively. Collectively, the reduction in OLCN architecture and association of OLCN measures with bone turnover, mechanics, and composition highlights the relevance of the osteocyte and the OLCN and a potential therapeutic target for treating diabetic skeletal fragility.

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.

Deep learning to overcome Zernike phase-contrast nanoCT artifacts for automated micro-nano porosity segmentation in bone

Bone material contains a hierarchical network of micro- and nano-cavities and channels, known as the lacuna-canalicular network (LCN), that is thought to play an important role in mechanobiology and turnover. The LCN comprises micrometer-sized lacunae, voids that house osteocytes, and submicrometer-sized canaliculi that connect bone cells. Characterization of this network in three dimensions is crucial for many bone studies. To quantify X-ray Zernike phase-contrast nanotomography data, deep learning is used to isolate and assess porosity in artifact-laden tomographies of zebrafish bones. A technical solution is proposed to overcome the halo and shade-off domains in order to reliably obtain the distribution and morphology of the LCN in the tomographic data. Convolutional neural network (CNN) models are utilized with increasing numbers of images, repeatedly validated by `error loss' and `accuracy' metrics. U-Net and Sensor3D CNN models were trained on data obtained from two different synchrotron Zernike phase-contrast transmission X-ray microscopes, the ANATOMIX beamline at SOLEIL (Paris, France) and the P05 beamline at PETRA III (Hamburg, Germany). The Sensor3D CNN model with a smaller batch size of 32 and a training data size of 70 images showed the best performance (accuracy 0.983 and error loss 0.032). The analysis procedures, validated by comparison with human-identified ground-truth images, correctly identified the voids within the bone matrix. This proposed approach may have further application to classify structures in volumetric images that contain non-linear artifacts that degrade image quality and hinder feature identification.

Osteocyte Mechanotransduction in Orthodontic Tooth Movement

PURPOSE OF REVIEW

Orthodontic tooth movement is characterized by periodontal tissue responses to mechanical loading, leading to clinically relevant functional adaptation of jaw bone. Since osteocytes are significant in mechanotransduction and orchestrate osteoclast and osteoblast activity, they likely play a central role in orthodontic tooth movement. In this review, we attempt to shed light on the impact and role of osteocyte mechanotransduction during orthodontic tooth movement.

RECENT FINDINGS

Mechanically loaded osteocytes produce signaling molecules, e.g., bone morphogenetic proteins, Wnts, prostaglandins, osteopontin, nitric oxide, sclerostin, and RANKL, which modulate the recruitment, differentiation, and activity of osteoblasts and osteoclasts. The major signaling pathways activated by mechanical loading in osteocytes are the wingless-related integration site (Wnt)/β-catenin and RANKL pathways, which are key regulators of bone metabolism. Moreover, osteocytes are capable of orchestrating bone adaptation during orthodontic tooth movement. A better understanding of the role of osteocyte mechanotransduction is crucial to advance orthodontic treatment. The optimal force level on the periodontal tissues for orthodontic tooth movement producing an adequate biological response, is debated. This review emphasizes that both mechanoresponses and inflammation are essential for achieving tooth movement clinically. To fully comprehend the role of osteocyte mechanotransduction in orthodontic tooth movement, more knowledge is needed of the biological pathways involved. This will contribute to optimization of orthodontic treatment and enhance patient outcomes.

Computation of physiological loading induced interstitial fluid motion in muscle standardized femur: Healthy vs. osteoporotic bone

BACKGROUND AND OBJECTIVES

Physiological loading-induced mechanical environments regulate bone modeling and remodeling. Thus, loading-induced normal strain is typically considered a stimulus to osteogenesis. However, several studies noticed new bone formation near the sites of minimal normal strain, e.g., the neutral axis of bending in long bones, which raises a question on how bone mass is maintained near these sites. Secondary mechanical components such as shear strain and interstitial fluid flow also stimulate bone cells and regulate bone mass. However, the osteogenic potential of these components is not well established. Accordingly, the present study estimates the distribution of physiological muscle loading-induced mechanical environments such as normal strain, shear strain, pore pressure, and interstitial fluid flow in long bones.

METHODS

A poroelastic finite element muscle standardized femur (MuscleSF) model is developed to compute the distribution of the mechanical environment as a function of bone porosities associated with osteoporotic and disuse bone loss.

RESULTS

The results indicate the presence of higher shear strain and interstitial fluid motion near the minimal strain sites, i.e., the neutral axis of bending of femoral cross-sections. This suggests that secondary stimuli may maintain the bone mass at these locations. Pore pressure and interstitial fluid motion reduce with the increased porosity associated with bone disorders, possibly resulting in diminished skeletal mechano-sensitivity to exogenous loading.

CONCLUSIONS

These outcomes present a better understanding of mechanical environment-mediated regulation of site-specific bone mass, which can be beneficial in developing prophylactic exercise to prevent bone loss in osteoporosis and muscle disuse.

Osteocyte lacunae in transiliac bone biopsy samples across life span

Osteocytes act as bone mechanosensors, regulators of osteoblast/osteoclast activity and mineral homeostasis, however, knowledge about their functional/morphological changes throughout life is limited. We used quantitative backscattered electron imaging (qBEI) to investigate osteocyte lacunae sections (OLS) as a 2D-surrogate characterizing the osteocytes. OLS characteristics, the density of mineralized osteocyte lacunae (i.e., micropetrotic osteocytes, md.OLS-Density in nb/mm²) and the average degree of mineralization (Ca_Mean in weight% calcium) of cortex and spongiosa were analyzed in transiliac biopsy samples from healthy individuals under 30 (n=59) and over 30 years (n=50) (i.e., before and after the age of peak bone mass, respectively). We found several differences in OLS-characteristics: 1). Inter-individually between the age groups: OLS-Density and OLS-Porosity were reduced by about 20% in older individuals in spongiosa and in cortex versus younger probands (both, p < 0.001). 2). Intra-individually between bone compartments: OLS-Density was higher in the cortex, +18.4%, p < 0.001 for younger and +7.6%, p < 0.05 for older individuals. Strikingly, the most frequent OLS nearest-neighbor distance was about 30 µm in both age groups and at both bone sites revealing a preferential organization of osteocytes in clusters. OLS-Density was negatively correlated with Ca_Mean in both spongiosa and cortex (both, p < 0.001). Few mineralized OLS were found in young individuals along with an increase of md.OLS-Density with age. In summary, this transiliac bone sample analysis of 200000 OLS from 109 healthy individuals throughout lifespan reveals several age-related differences in OLS characteristics. Moreover, our study provides reference data from healthy individuals for different ages to be used for diagnosis of bone abnormalities in diseases. STATEMENT OF SIGNIFICANCE: Osteocytes are bone cells embedded in lacunae within the mineralized bone matrix and have a key role in the bone metabolism and the mineral homeostasis. Not easily accessible, we used quantitative backscattered electron imaging to determine precisely number and shape descriptors of the osteocyte lacunae in 2D. We analyzed transiliac biopsy samples from 109 individuals with age distributed from 2 to 95 years. Compact cortical bone showed constantly higher lacunar density than cancellous bone but the lacunar density in both bone tissue decreased with age before the peak bone mass age at 30 years and stabilized or even increased after this age. This extensive study provides osteocyte lacunae reference data from healthy individuals usable for bone pathology diagnosis.

Subcanalicular Nanochannel Volume Is Inversely Correlated With Calcium Content in Human Cortical Bone

The spatial distribution of mineralization density is an important signature of bone growth and remodeling processes, and its alterations are often related to disease. The extracellular matrix of some vertebrate mineralized tissues is known to be perfused by a lacunocanalicular network (LCN), a fluid-filled unmineralized structure that harbors osteocytes and their fine processes and transports extracellular fluid and its constituents. The current report provides evidence for structural and compositional heterogeneity at an even smaller, subcanalicular scale. The work reveals an extensive unmineralized three-dimensional (3D) network of nanochannels (~30 nm in diameter) penetrating the mineralized extracellular matrix of human femoral cortical bone and encompassing a greater volume fraction and surface area than these same parameters of the canaliculi comprising the LCN. The present study combines high-resolution focused ion beam-scanning electron microscopy (FIB-SEM) to investigate bone ultrastructure in 3D with quantitative backscattered electron imaging (qBEI) to estimate local bone mineral content. The presence of nanochannels has been found to impact qBEI measurements fundamentally, such that volume percentage (vol%) of nanochannels correlates inversely with weight percentage (wt%) of calcium. This mathematical relationship between nanochannel vol% and calcium wt% suggests that the nanochannels could potentially provide space for ion and small molecule transport throughout the bone matrix. Collectively, these data propose a reinterpretation of qBEI measurements, accounting for nanochannel presence in human bone tissue in addition to collagen and mineral. Further, the results yield insight into bone mineralization processes at the nanometer scale and present the possibility for a potential role of the nanochannel system in permitting ion and small molecule diffusion throughout the extracellular matrix. Such a possible function could thereby lead to the sequestration or occlusion of the ions and small molecules within the extracellular matrix. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

On computational predictions of fluid flow and its effects on bone healing in dental implant treatments: an investigation of spatiotemporal fluid flow in cyclic loading

Fluid flow in (porous) bone plays an important role in its maintenance, adaptation, and healing after an injury. Experimental and computational studies apply mechanical loading on bone to predict fluid flow development and/or to find its material properties. In most cases, mechanical loading is applied as a linear function in time. Multiple loading functions-with identical peak load and loading frequency-were used to investigate load-induced fluid flow and predict bone healing surrounding a dental implant. Implementing an instantaneous healing stimulus led to major differences in healing predictions for slightly different loading functions. Load-induced fluid flow was found to be displacement-rate dependent with complex spatial-temporal variations and not necessarily symmetrical during loading and unloading phases. Haversine loading resulted in more numerical stability compared to ramped/triangular loading, providing the opportunity for further investigation of the effects of various physiological masticatory loadings. It was concluded that using the average healing stimulus during cyclic loading gives the most robust bone healing predictions.

Mechanoregulated trabecular bone adaptation: Progress report on in silico approaches

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

A morphometric analysis of the osteocyte canaliculus using applied automatic semantic segmentation by machine learning

INTRODUCTION

Osteocytes play a role as mechanosensory cells by sensing flow-induced mechanical stimuli applied on their cell processes. High-resolution imaging of osteocyte processes and the canalicular wall are necessary for the analysis of this mechanosensing mechanism. Focused ion beam-scanning electron microscopy (FIB-SEM) enabled the visualization of the structure at the nanometer scale with thousands of serial-section SEM images. We applied machine learning for the automatic semantic segmentation of osteocyte processes and canalicular wall and performed a morphometric analysis using three-dimensionally reconstructed images.

MATERIALS AND METHODS

Six-week-old-mice femur were used. Osteocyte processes and canaliculi were observed at a resolution of 2 nm/voxel in a 4 × 4 μm region with 2000 serial-section SEM images. Machine learning was used for automatic semantic segmentation of the osteocyte processes and canaliculi from serial-section SEM images. The results of semantic segmentation were evaluated using the dice similarity coefficient (DSC). The segmented data were reconstructed to create three-dimensional images and a morphological analysis was performed.

RESULTS

The DSC was > 83%. Using the segmented data, a three-dimensional image of approximately 3.5 μm in length was reconstructed. The morphometric analysis revealed that the median osteocyte process diameter was 73.8 ± 18.0 nm, and the median pericellular fluid space around the osteocyte process was 40.0 ± 17.5 nm.

CONCLUSION

We used machine learning for the semantic segmentation of osteocyte processes and canalicular wall for the first time, and performed a morphological analysis using three-dimensionally reconstructed images.

Effects of Osteocyte Shape on Fluid Flow and Fluid Shear Stress of the Loaded Bone

This study was conducted to better understand the specific behavior of the intraosseous fluid flow. We calculated the number and distribution of bone canaliculi around the osteocytes based on the varying shapes of osteocytes. We then used these calculated parameters and other bone microstructure data to estimate the anisotropy permeability of the lacunar-canalicular network. Poroelastic finite element models of the osteon were established, and the influence of the osteocyte shape on the fluid flow properties of osteons under an axial displacement load was analyzed. Two types of boundary conditions (BC) that might occur in physiological environments were considered on the cement line of the osteon. BC1 allows free fluid passage from the outer elastic restraint boundary, and BC2 is impermeable and allows no free fluid passage from outer displacement constrained boundary. They both have the same inner boundary conditions that allow fluid to pass through. Changes in the osteocyte shape altered the maximum value of pressure gradient (PG), pore pressure (PP), fluid velocity (FV), and fluid shear stress (FSS) relative to the reference model (spherical osteocytes). The maximum PG, PP, FV, and FSS in BC2 were nearly 100% larger than those in BC1, respectively. It is found that the BC1 was closer to the real physiological environment. The fluid flow along different directions in the elongated osteocyte model was more evident than that in other models, which may have been due to the large difference in permeability along different directions. Changes in osteocyte shape significantly affect the degrees of anisotropy of fluid flow and porous media of the osteon. The model presented in this study can accurately quantify fluid flow in the lacunar-canalicular network.

Bone permeability and mechanotransduction: Some current insights into the function of the lacunar-canalicular network

Lacunar-canalicular (LC) permeability involves the passage of fluids, nutrients, oxygen, ions, and signalling molecules through bone tissue, facilitating the maintenance of bone vitality and function and responses to various physiological conditions and diseases. LC permeability and fluid flow-shear stress/drag force play important roles in mechanotransduction in bone tissue by inducing mechanical stimuli in osteocytes, modulating cellular functions, and determining bone adaptation. Alterations in LC structure may therefore influence the fluid flow pattern through the LC network, thereby affecting the ability of osteocytes to sense and translate mechanical signals and possibly contributing to bone remodelling. Several bone-health conditions are associated with changes in LC structure and function and may affect mechanotransduction and responses, although the mechanisms underlying these associations are still not fully understood. In this review, recent studies of LC networks, their formation and transfer mechanical stimuli, and changes in structure, functional permeability, and mechanotransduction that result from age, pathology, and mechanical loading are discussed. Additionally, applications of vibration and low-intensity pulsed ultrasound in bone healthcare and regeneration fields are also presented.

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.

Simulation study on the effect of resistance exercise on the hydrodynamic microenvironment of osteocytes in microgravity

Osteoporosis occurs in astronauts after long-term space flight owing to the lack of gravity. The mechanical microenvironment of osteocytes in load-bearing bone are changed during resistance exercise, which prevents massive bone loss in the human body. A cylindrical fluid-structure coupling finite element model for osteons with a two-stage pore structure (i.e., Haversian canal, lacunar-canalicular system) was established with the software COMSOL. In the Earth's gravity field and in microgravity, considering the effects of pulsating pressure of arterioles, a comparative study was performed on the changes in hydrodynamic microenvironment of osteocytes during human body high-intensity exercise at different frequencies (defined as causing bone to produce 3000 με) and the body is at rest. Positive and negative liquid pressure (with respect to one atmosphere pressure) alternately acted on osteocytes during human exercising, but only positive pressure acted on osteocytes during human resting. The variation range of liquid pressure acted on osteocytes during human exercising was significantly higher than that during resting. The liquid flow velocity around osteocytes during body exercise was about four orders of magnitude higher than that during resting. In microgravity, moderate physical exercise can obviously improve the hydrodynamic microenvironment of osteocytes in load-bearing bone, which could compensate for the lack of mechanical stimulation to osteocytes caused by the lack of gravity, thereby promoting the normal physiological function of osteocytes. To a certain extent, these results revealed the biomechanical mechanism by which exercise has an effect in fighting osteoporosis in astronauts.

Anisotropic solid-state PLA foaming templated by crystal phase pre-oriented with 3D printing: Cell supporting structures with directional capillary transfer function

Bones represent a superb biomaterial that combines high mechanical stiffness with nutrition delivery to its osteocyte cells through the microscopical Haversian canals and bone canaliculi. Such structure is hard to...

Multiscale modeling of bone tissue mechanobiology

Mechanical environment has a crucial role in our organism at the different levels, ranging from cells to tissues and our own organs. This regulatory role is especially relevant for bones, given their importance as load-transmitting elements that allow the movement of our body as well as the protection of vital organs from load impacts. Therefore bone, as living tissue, is continuously adapting its properties, shape and repairing itself, being the mechanical loads one of the main regulatory stimuli that modulate this adaptive behavior. Here we review some key results of bone mechanobiology from computational models, describing the effect that changes associated to the mechanical environment induce in bone response, implant design and scaffold-driven bone regeneration.

High-resolution image-based simulation reveals membrane strain concentration on osteocyte processes caused by tethering elements

Osteocytes are vital for regulating bone remodeling by sensing the flow-induced mechanical stimuli applied to their cell processes. In this mechanosensing mechanism, tethering elements (TEs) connecting the osteocyte process with the canalicular wall potentially amplify the strain on the osteocyte processes. The ultrastructure of the osteocyte processes and canaliculi can be visualized at a nanometer scale using high-resolution imaging via ultra-high voltage electron microscopy (UHVEM). Moreover, the irregular shapes of the osteocyte processes and the canaliculi, including the TEs in the canalicular space, should considerably influence the mechanical stimuli applied to the osteocytes. This study aims to characterize the roles of the ultrastructure of osteocyte processes and canaliculi in the mechanism of osteocyte mechanosensing. Thus, we constructed a high-resolution image-based model of an osteocyte process and a canaliculus using UHVEM tomography and investigated the distribution and magnitude of flow-induced local strain on the osteocyte process by performing fluid–structure interaction simulation. The analysis results reveal that local strain concentration in the osteocyte process was induced by a small number of TEs with high tension, which were inclined depending on the irregular shapes of osteocyte processes and canaliculi. Therefore, this study could provide meaningful insights into the effect of ultrastructure of osteocyte processes and canaliculi on the osteocyte mechanosensing mechanism.

Closing the osteon: Do osteocytes sense strain rate rather than fluid flow?

Osteons are cylindrical structures of bone created by matrix resorbing osteoclasts, followed by osteoblasts that deposit new bone. Osteons align with the principal loading direction and it is thought that the osteoclasts are directed by osteocytes, the mechanosensitive cells that reside inside the bone matrix. These osteocytes are presumably controlled by interstitial fluid flow, induced by the physiological loading of bones. Here I consider the stimulation of osteocytes while the osteon is closed by osteoblasts. In a conceptual finite element model, bone is considered a poro-elastic material and subjected to locomotion-induced loading conditions. It appears that the magnitude of flow is constant along the closing cone, while shear strain rate in the bone matrix diminishes linearly with the deposition of bone. This suggests that shear strain rate, rather than fluid flow, is the physical cue that controls osteocytes and bone deposition in newly formed osteons.

Impaired 1,25 dihydroxyvitamin D3 action and hypophosphatemia underlie the altered lacuno-canalicular remodeling observed in the Hyp mouse model of XLH

Osteocytes remodel the perilacunar matrix and canaliculi. X-linked hypophosphatemia (XLH) is characterized by elevated serum levels of fibroblast growth factor 23 (FGF23), leading to decreased 1,25 dihydroxyvitamin D₃ (1,25D) production and hypophosphatemia. Bones from mice with XLH (Hyp) have enlarged osteocyte lacunae, enhanced osteocyte expression of genes of bone remodeling, and impaired canalicular structure. The altered lacuno-canalicular (LCN) phenotype is improved with 1,25D or anti-FGF23 antibody treatment, pointing to roles for 1,25D and/or phosphate in regulating this process. To address whether impaired 1,25D action results in LCN alterations, the LCN phenotype was characterized in mice lacking the vitamin D receptor (VDR) in osteocytes (VDR^f/f;DMP1Cre+). Mice lacking the sodium phosphate transporter NPT2a (NPT2aKO) have hypophosphatemia and high serum 1,25D levels, therefore the LCN phenotype was characterized in these mice to determine if increased 1,25D compensates for hypophosphatemia in regulating LCN remodeling. Unlike Hyp mice, neither VDR^f/f;DMP1Cre+ nor NPT2aKO mice have dramatic alterations in cortical microarchitecture, allowing for dissecting 1,25D and phosphate specific effects on LCN remodeling in tibial cortices. Histomorphometric analyses demonstrate that, like Hyp mice, tibiae and calvariae in VDR^f/f;DMP1Cre+ and NPT2aKO mice have enlarged osteocyte lacunae (tibiae: 0.15±0.02μm²(VDR^f/f;DMP1Cre-) vs 0.19±0.02μm²(VDR^f/f;DMP1Cre+), 0.12±0.02μm²(WT) vs 0.18±0.0μm²(NPT2aKO), calvariae: 0.09±0.02μm²(VDR^f/f;DMP1Cre-) vs 0.11±0.02μm²(VDR^f/f;DMP1Cre+), 0.08±0.02μm²(WT) vs 0.13±0.02μm²(NPT2aKO), p<0.05 all comparisons) and increased immunoreactivity of bone resorption marker Cathepsin K (Ctsk). The osteocyte enriched RNA isolated from tibiae in VDR^f/f;DMP1Cre+ and NPT2aKO mice have enhanced expression of matrix resorption genes that are classically expressed by osteoclasts (Ctsk, Acp5, Atp6v0d2, Nhedc2). Treatment of Ocy454 osteocytes with 1,25D or phosphate inhibits the expression of these genes. Like Hyp mice, VDR^f/f;DMP1Cre+ and NPT2aKO mice have impaired canalicular organization in tibia and calvaria. These studies demonstrate that hypophosphatemia and osteocyte-specific 1,25D actions regulate LCN remodeling. Impaired 1,25D action and low phosphate levels contribute to the abnormal LCN phenotype observed in XLH.

Physiological Loading-Induced Interstitial Fluid Dynamics in Osteon of Osteogenesis Imperfecta Bone

Osteogenesis imperfecta (OI), also known as "brittle bone disease," is a genetic bone disorder. OI bones experience frequent fractures. Surgical procedures are usually followed by clinicians in the management of OI. It has been observed physical activity is equally beneficial in reducing OI bone fractures in both children and adults as mechanical stimulation improves bone mass and strength. Loading-induced mechanical strain and interstitial fluid flow stimulate bone remodeling activities. Several studies have characterized strain environment in OI bones, whereas very few studies attempted to characterize the interstitial fluid flow. OI significantly affects bone micro-architecture. Thus, this study anticipates that canalicular fluid flow reduces in OI bone in comparison to the healthy bone in response to physiological loading due to altered poromechanical properties. This work attempts to understand the canalicular fluid distribution in single osteon models of OI and healthy bone. A poromechanical model of osteon is developed to compute pore-pressure and interstitial fluid flow as a function of gait loading pattern reported for OI and healthy subjects. Fluid distribution patterns are compared at different time-points of the stance phase of the gait cycle. It is observed that fluid flow significantly reduces in OI bone. Additionally, flow is more static than dynamic in OI osteon in comparison to healthy subjects. This work attempts to identify the plausible explanation behind the diminished mechanotransduction capability of OI bone. This work may further be extended for designing better biomechanical therapies to enhance the fluid flow in order to improve osteogenic activities in OI bone.

Secondary osteon variants and remodeling in human bone

Histomorphometric analysis of human cortical bone has documented the occurrence of secondary osteon variants. These include drifting osteons which form tails as they move erratically through the cortex and Type II osteons which show partial resorption and redeposition within the cement line of the osteon. Little is known about the biological significance of these variants. Prior studies suggested correlations with age, biomechanics, diet, and mineral homeostasis. No study has yet tested for osteon variant associations with static measures of bone remodeling. In this study, thin sections (n = 112) of the posterior femur representing a late English Medieval adult human osteological collection, subdivided by age, sex, and socio-economic status, were examined to determine whether remodeling indicators reconstructed from osteon parameters (area, diameter, area ratios) and densities differed between categories of presence or absence of Type II and drifting osteon variants. Of the 112 sections, 33 presented with Type II osteons, and 38 had drifting osteons. Sporadic statistically significant results were identified. Haversian canal:osteon area ratio differed (p =  0.017) with Type II osteon presence, Type II osteons were more prevalent in males than females (p =  0.048), and drifting osteons were associated with smaller osteon (p =  0.049) and Haversian canal area (p =  0.05). These results may be explained through some biological (sex) and social (status) processes such as a period of physiological recovery (e.g., following lactation, malnutrition). However, the general lack of consistent relationships between osteon variants and remodeling indicators suggests they occur as a result of natural variation.

Modeling cortical bone adaptation using strain gradients

Cyclic and low-magnitude loading promotes osteogenesis (i.e. new bone formation). Normal strain, strain energy density and fatigue damage accumulation are typically considered as osteogenic stimuli in computer models to predict site-specific new bone formation. These models however had limited success in explaining osteogenesis near the sites of minimal normal strain, for example, neutral axis of bending. Other stimuli such as fluid motion or strain gradient also stimulate bone formation. In silico studies modeled the new bone formation as a function of fluid motion, however, computation of fluid motion involves complex mathematical calculations. Strain gradients drive fluid flow and thus can also be established as the stimulus. Osteogenic potential of strain gradients is however not well established. The present study establishes strain gradients as osteogenic stimuli. Bending-induced strain gradients are computed at cortical bone cross-sections reported in animal loading in vivo studies. Correlation analysis between strain gradients and site of osteogenesis is analyzed. In silico model is also developed to test the osteogenic potential of strain gradients. The model closely predicts in vivo new bone distribution as a function of strain gradients. The outcome establishes strain gradient as computationally easy and robust stimuli to predict site-specific osteogenesis. The present study may be useful in the development of biomechanical approaches to mitigate bone loss.

Simulation of the mechanical behavior of osteons using artificial gravity devices in microgravity

Aviation medical research shows that disuse osteoporosis will occur after long-term space flight. Even with countermeasures such as exercise and drug treatments, this outcome cannot be avoided in flight. In recent years, the application of artificial gravity devices that change the mechanical microenvironment of bone in microgravity have shown promise in mitigating the risk of disuse osteoporosis. Considering the existence of osteocytes, a fluid-solid coupling finite element model for osteons with two-stage pore structure (Haversian canal, lacunar-canalicular system) was established. In order to study the changes in the mechanical behavior of osteocytes under the action of various artificial gravity (AG) devices, including long-arm centrifuge (LAC), short-arm centrifuge (SAC), and a lower body negative pressure (LBNP) chamber. In addition, the difference in pulsating pressure and static pressure stress caused by the gravity gradient under the AG devices was examined. The simulation results showed that the AG devices could effectively improve the stress level of osteocytes in microgravity. The mechanical microenvironment of osteocytes that was provided by the LAC was closest to that of the Earth's gravitational field. The mechanical stimulation on osteocytes was not significantly improved by the SAC, but from a practical viewpoint, it occupied less space than the LAC. The LBNP chamber created a higher level of stress for osteocytes. Therefore, the LAC was an ideal device for replacing Earth's gravitational field, except for the practical limitations of its physical size. In contrast, the LBNP device had the greatest application potential in training for its expansibility and convenience.

The mechanoresponse of bone is closely related to the osteocyte lacunocanalicular network architecture

The explanation of how bone senses and adapts to mechanical stimulation still relies on hypotheses. The fluid flow hypothesis claims that a load-induced fluid flow through the lacunocanalicular network can be sensed by osteocytes, which reside within the network structure. We show that considering the network architecture results in a better prediction of bone remodeling than mechanical strain alone. This was done by calculating the fluid flow through the lacunocanalicular network in bone volumes covering the complete cross-sections of mouse tibiae, which underwent controlled in vivo loading. The established relationship between mechanosensitivity and network architecture in individual animals implies possibilities for patient-specific therapies. A new connectomics approach to analyze lacunocanalicular network properties is necessary to understand skeletal mechanobiology.

Organisms rely on mechanosensing mechanisms to adapt to changes in their mechanical environment. Fluid-filled network structures not only ensure efficient transport but can also be employed for mechanosensation. The lacunocanalicular network (LCN) is a fluid-filled network structure, which pervades our bones and accommodates a cell network of osteocytes. For the mechanism of mechanosensation, it was hypothesized that load-induced fluid flow results in forces that can be sensed by the cells. We use a controlled in vivo loading experiment on murine tibiae to test this hypothesis, whereby the mechanoresponse was quantified experimentally by in vivo micro-computed tomography (µCT) in terms of formed and resorbed bone volume. By imaging the LCN using confocal microscopy in bone volumes covering the entire cross-section of mouse tibiae and by calculating the fluid flow in the three-dimensional (3D) network, we could perform a direct comparison between predictions based on fluid flow velocity and the experimentally measured mechanoresponse. While local strain distributions estimated by finite-element analysis incorrectly predicts preferred bone formation on the periosteal surface, we demonstrate that additional consideration of the LCN architecture not only corrects this erroneous bias in the prediction but also explains observed differences in the mechanosensitivity between the three investigated mice. We also identified the presence of vascular channels as an important mechanism to locally reduce fluid flow. Flow velocities increased for a convergent network structure where all of the flow is channeled into fewer canaliculi. We conclude that, besides mechanical loading, LCN architecture should be considered as a key determinant of bone adaptation.

Human and mouse bones physiologically integrate in a humanized mouse model while maintaining species-specific ultrastructure

Humanized mouse models are increasingly studied to recapitulate human-like bone physiology. While human and mouse bone architectures differ in multiple scales, the extent to which chimeric human-mouse bone physiologically interacts and structurally integrates remains unknown. Here, we identify that humanized bone is formed by a mosaic of human and mouse collagen, structurally integrated within the same bone organ, as shown by immunohistochemistry. Combining this with materials science techniques, we investigate the extracellular matrix of specific human and mouse collagen regions. We show that human-like osteocyte lacunar-canalicular network is retained within human collagen regions and is distinct to that of mouse tissue. This multiscale analysis shows that human and mouse tissues physiologically integrate into a single, functional bone tissue while maintaining their species-specific ultrastructural differences. These results offer an original method to validate and advance tissue-engineered human-like bone in chimeric animal models, which grow to be eloquent tools in biomedical research.

Using confocal imaging approaches to understand the structure and function of osteocytes and the lacunocanalicular network

Although overlooked in the past, osteocytes have come to the forefront of skeletal biology and are now recognized as a key cell type that integrates hormonal, mechanical and other signals to control bone mass through regulation of both osteoblast and osteoclast activity. With the surge of recent interest in osteocytes as bone regulatory cells and the discovery that they also function as endocrine regulators of phosphate homeostasis, there has been renewed interest in understanding the structure and function of these unique and relatively inaccessible cells. Osteocytes are embedded within the mineralized bone matrix and are housed within a complex lacunocanalicular system which connects them with the circulation and with other organ systems. This has presented unique challenges for imaging these cells. This review summarizes recent advances in confocal imaging approaches for visualizing osteocytes and their lacunocanalicular networks in both living and fixed bone specimens and discusses how computational approaches can be combined with live and fixed cell imaging techniques to generate quantitative outputs and predictive models. The integration of advanced imaging with computational approaches promises to lead to a more in depth understanding of the structure and function of osteocyte networks and the lacunocanalicular system in the healthy and aging state as well as in pathological conditions in bone.

Signalling molecule transport analysis in lacunar-canalicular system

Mechanical loading-induced fluid flow in lacunar-canalicular space (LCS) of bone excites osteocyte cells to release signalling molecules which initiate osteo-activities. Theoretical models considered canaliculi as a uniform and symmetrical space/channel in bone. However, experimental studies reported that canalicular walls are irregular and curvy resulting in inhomogeneous fluid motion which may influence the molecular transport. Therefore, a new mathematical model of LCS with curvy canalicular walls is developed to characterize cantilever bending-induced canalicular flow behaviour in terms of pore-pressure, fluid velocity, and streamlines. The model also analyses the mobility of signalling molecules involved in bone mechanotransduction as a function of loading frequency and permeability of LCS. Inhomogeneous flow is observed at higher loading frequency which amplifies mechanotransduction; nevertheless, it also promotes trapping of signalling molecules. The effects of shape and size of signalling molecules on transport behaviour are also studied. Trivially, signalling molecules larger in size and weight move slower as compared to molecules small in size and weight which validates the findings of the present study. The outcomes will ultimately be useful in designing better biomechanical exercise in combination with pharmaceutical agents to improve the bone health.