Measurement and estimation of osteocyte mechanical strain

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.

Biodegradation-affected fatigue behavior of extrusion-based additively manufactured porous iron-manganese scaffolds

Additively manufactured (AM) biodegradable porous iron-manganese (FeMn) alloys have recently been developed as promising bone-substituting biomaterials. However, their corrosion fatigue behavior has not yet been studied. Here, we present the first study on the corrosion fatigue behavior of an extrusion-based AM porous Fe35Mn alloy under cyclic loading in air and in the revised simulated body fluid (r-SBF), including the fatigue crack morphology and distribution in the porous structure. We hypothesized that the fatigue behavior of the architected AM Fe35Mn alloy would be strongly affected by the simultaneous biodegradation process. We defined the endurance limit as the maximum stress at which the scaffolds could undergo 3 million loading cycles without failure. The endurance limit of the scaffolds was determined to be 90 % of their yield strength in air, but only 60 % in r-SBF. No notable crack formation in the specimens tested in air was observed even after loading up to 90 % of their yield strength. As for the specimens tested in r-SBF, however, cracks formed in the specimens subjected to loads exceeding 60 % of their yield strength appeared to initiate on the periphery and propagate toward the internal struts. Altogether, the results show that the extrusion-based AM porous Fe35Mn alloy is capable of tolerating up to 60 % of its yield strength for up to 3 million cycles, which corresponds to 1.5 years of use of load-bearing implants subjected to repetitive gait cycles. The fatigue performance of the alloy thus further enhances its potential for trabecular bone substitution subjected to cyclic compressive loading. STATEMENT OF SIGNIFICANCE: Fatigue behavior of extrusion-based AM porous Fe35Mn alloy scaffolds in air and revised simulated body fluid was studied. The Fe35Mn alloy scaffolds endured 90 % of their yield strength for up to 3 × 106 loading cycles in air. Moreover, the scaffolds tolerated 3 × 106 loading cycles at 60 % of their yield strength in revised simulated body fluid. The Fe35Mn alloy scaffolds exhibited a capacity of withstanding 1.5-year physiological loading when used as bone implants.

Design and manufacturing of biomimetic scaffolds for bone repair inspired by bone trabeculae

Porous scaffold (PorS) implants, particularly those that mimic the structural features of natural cancellous bone (NCanB), are increasingly essential for the treatment of large-area bone defects. However, the mechanical properties of NCanB-based bionic bone scaffold (BioS) and its performance as a bone repair material have not been fully explored. This study investigates the effect of bionic structure parameters on the mechanical properties and bone reconstruction performance of BioS. Using laser powder bed fusion (L-PBF) technology, different BioS with various structural parameters were created and evaluated using Micro-CT, compression testing, Finite Element (FE) Simulation, and computational fluid dynamics (CFD), and compared to commonly used clinical PorS. Assess the capacity of the BioS scaffold to support and enhance bone reconstruction following implantation through the evaluation of its mechanical properties, permeability, and fluid shear stress (FSS). BioS-85-90 and BioS-80-50 showed suitable mechanical properties, performed well in FE simulation of implantation, demonstrated outstanding abilities for osteoinductive ingrowth and bone tissue differentiation, and proved to be reliable materials for the reconstruction of bone defects. Therefore, BioS shows significant potential for clinical application as a bone reconstruction material, providing a solid foundation for the integration of tissue engineering and bionic design.

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.

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.

Histochemical characteristics on minimodeling-based bone formation induced by anabolic drugs for osteoporotic treatment

Modeling, the changes of bone size and shape, often takes place at the developmental stages, whereas bone remodeling-replacing old bone with new bone-predominantly occurs in adults. Unlike bone remodeling, bone formation induced by modeling i.e., minimodeling (microscopic modeling in cancellous bone) is independent of osteoclastic bone resorption. Although recently-developed drugs for osteoporotic treatment could induce minimodeling-based bone formation in addition to remodeling-based bone formation, few reports have demonstrated the histological aspects of minimodeling-based bone formation. After administration of eldecalcitol or romosozumab, unlike teriparatide treatment, mature osteoblasts formed new bone by minimodeling, without developing thick preosteoblastic layers. The histological characteristics of minimodeling-based bone formation is quite different from remodeling, as it is not related to osteoclastic bone resorption, resulting in convex-shaped new bone and smooth cement lines called arrest lines. In this review, we will show histological properties of minimodeling-based bone formation by osteoporotic drugs.

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

Histomorphometric analysis of osteocyte lacunae in human and pig: exploring its potential for species discrimination

In recent years, several studies have focused on species discrimination of bone fragments by histological analysis. According to literature, the most consistent distinguishing features are Haversian canal and Haversian system areas. Nonetheless, there is a consistent overlap between human and non-human secondary osteon dimensions. One of the features that have never been analyzed for the purpose of species discrimination is the osteocyte lacuna, a small oblong cavity in which the osteocyte is locked in. The aim of this study is to verify whether there are significant quantitative differences between human and pig lacunae within secondary osteons with similar areas. Study sample comprises the midshaft of long bones (humerus, radius, ulna, femur, tibia, and fibula) of a medieval human adult and a juvenile pig. Sixty-eight secondary osteons with similar areas have been selected for each species and a total of 1224 osteocyte lacunae have been measured. For each osteon, the total number of lacunae was counted, and the following measurements were taken: minimum and maximum diameter, area, perimeter, and circularity of nine lacunae divided between inner, intermediate, and outer lacunae. Statistical analysis showed minimal differences between human and pig in the number of lacunae per osteons and in the minimum diameter (P > 0.05). On the contrary, a significant difference (P < 0.001) has been observed in the maximum diameter, perimeter, area, and circularity. Although there is the need for further research on different species and larger sample, these results highlighted the potential for the use of osteocyte lacunae as an additional parameter for species discrimination. Concerning the difference between the dimensions of osteocyte lacunae based on their position within the osteon (inner, intermediate, and outer lacunae), results showed that their size decreases from the cement line towards the Haversian canal both in human and pig.

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.

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.

Diet and Exercise: a Match Made in Bone

PURPOSE OF REVIEW

Multiple dietary components have the potential to positively affect bone mineral density in early life and reduce loss of bone mass with aging. In addition, regular weight-bearing physical activity has a strong positive effect on bone through activation of osteocyte signaling. We will explore possible synergistic effects of dietary components and mechanical stimuli for bone health by identifying dietary components that have the potential to alter the response of osteocytes to mechanical loading.

RECENT FINDINGS

Several (sub)cellular aspects of osteocytes determine their signaling towards osteoblasts and osteoclasts in response to mechanical stimuli, such as the osteocyte cytoskeleton, estrogen receptor α, the vitamin D receptor, and the architecture of the lacunocanalicular system. Potential modulators of these features include 1,25-dihydroxy vitamin D3, several forms of vitamin K, and the phytoestrogen genistein. Multiple dietary components potentially affect osteocyte function and therefore may have a synergistic effect on bone health when combined with a regime of physical activity.

Aging, Osteocytes, and Mechanotransduction

PURPOSE OF REVIEW

The bone is able to adapt its structure to mechanical signals via the bone remodeling process governed by mechanosensitive osteocytes. With aging, an imbalance in bone remodeling results in osteoporosis. In this review, we hypothesized that changes in lacunar morphology underlie the decreased bone mechanoresponsiveness to mechanical loading with aging.

RECENT FINDINGS

Several studies have reported considerable variations in the shape of osteocytes and their lacunae with aging. Since osteocytes can sense matrix strain directly via their cell bodies, the variations in osteocyte morphology may cause changes in osteocyte mechanosensitivity. As a consequence, the load-adaptive response of osteocytes may change with aging, even when mechanical loading would remain unchanged. Though extensive quantitative data is lacking, evidence exists that the osteocyte lacunae are becoming smaller and more spherical with aging. Future dedicated studies might reveal whether these changes would affect osteocyte mechanosensation and the subsequent biological response, and whether this is (one of) the pathways involved in age-related bone loss.

The response of anosteocytic bone to controlled loading

The bones of the skeleton of most advanced teleost fish do not contain osteocytes. Considering the pivotal role assigned to osteocytes in the process of modeling and remodeling (the adaptation of external and internal bone structure and morphology to external loads and the repair of areas with micro-damage accumulation, respectively) it is unclear how, and even whether, their skeleton can undergo modeling and remodeling. Here, we report on the results of a study of controlled loading of the anosteocytic opercula of tilapia (Oreochromis aureus). Using a variety of microscopy techniques we show that the bone of the anosteocytic tilapia actively adapts to applied loads, despite the complete absence of osteocytes. We show that in the directly loaded area, the response involves a combination of bone resorption and bone deposition; we interpret these results and the structure of the resultant bone tissue to mean that both modeling and remodeling are taking place in response to load. We further show that adjacent to the loaded area, new bone is deposited in an organized, layered manner, typical of a modeling process. The material stiffness of the newly deposited bone is higher than that of the bone which was present prior to loading. The absence of osteocytes requires another candidate cell for mechanosensing and coordinating the modeling process, with osteoblasts seeming the most likely candidates.

TECHNOLOGIES FOR STRAIN ASSESSMENT FROM WHOLE BONE TO MINERALIZED OSTEOID LEVEL: A CRITICAL REVIEW

Bone has distinctive structures and mechanical properties at the whole bone, perilacunar and mineralized osteoid levels. A systematic understanding of bone strain magnitudes at different anatomical levels and their internal interactions is the prerequisite to advances in bone mechanobiology. However, due to the intrinsic shortcomings of the strain-measuring technologies, the systematic assessment of bone strain at different anatomical levels under physiological conditions and a deep understanding of their internal interactions are still restricted. To promote technological advances and provide systematic and valuable information for mechanical engineers and bone biomechanical researchers, the most useful methods for measuring bone strain at different anatomical levels are demonstrated in this review, and suggestions for the future development of the technologies and their potential integrated applications are proposed.

The Cementocyte-An Osteocyte Relative?

Cementum is a mineralized tissue covering the tooth root that functions in tooth attachment and posteruptive adjustment of tooth position. During formation of the apically located cellular cementum, some cementoblasts become embedded in the cementoid matrix and become cementocytes. As apparently terminally differentiated cells embedded in a mineralized extracellular matrix, cementocytes are part of a select group of specialized cells, also including osteocytes, hypertrophic chondrocytes, and odontoblasts. The differentiation and potential function(s) of cementocytes are virtually unknown, and the question may be posed whether the cementocyte is a dynamic actor in cementum in comparable fashion with the osteocyte in the skeleton, responding to changing tooth functions and endocrine signals and actively directing local cementum metabolism. This review summarizes the literature with regard to cementocytes, comparing them to their closest "cousins," the osteocytes, where insights gained from osteocyte studies serve to inform the critical examination of cementocytes. The review identifies important unanswered questions about these cells regarding their origins, differentiation, morphology and lacuno-canalicular system, selective markers, and potential functions.

Bone strain magnitude is correlated with bone strain rate in tetrapods: implications for models of mechanotransduction

Hypotheses suggest that structural integrity of vertebrate bones is maintained by controlling bone strain magnitude via adaptive modelling in response to mechanical stimuli. Increased tissue-level strain magnitude and rate have both been identified as potent stimuli leading to increased bone formation. Mechanotransduction models hypothesize that osteocytes sense bone deformation by detecting fluid flow-induced drag in the bone's lacunar-canalicular porosity. This model suggests that the osteocyte's intracellular response depends on fluid-flow rate, a product of bone strain rate and gradient, but does not provide a mechanism for detection of strain magnitude. Such a mechanism is necessary for bone modelling to adapt to loads, because strain magnitude is an important determinant of skeletal fracture. Using strain gauge data from the limb bones of amphibians, reptiles, birds and mammals, we identified strong correlations between strain rate and magnitude across clades employing diverse locomotor styles and degrees of rhythmicity. The breadth of our sample suggests that this pattern is likely to be a common feature of tetrapod bone loading. Moreover, finding that bone strain magnitude is encoded in strain rate at the tissue level is consistent with the hypothesis that it might be encoded in fluid-flow rate at the cellular level, facilitating bone adaptation via mechanotransduction.

Childhood cortical porosity is related to microstructural properties of the bone-muscle junction

Childhood cortical porosity is attributable to giant asymmetrical drifting osteonal canals that arise predominantly along the primary-secondary bone interface (PSBI). Bone from the external iliac crest cortex of 92 subjects aged 0 to 25 years was examined histomorphometrically for differences in microstructural properties between primary and secondary bone that might account for features of drifting osteonal canals. Primary compared with secondary bone showed greater numbers of osteocyte lacunae, thinner collagen lamellae, and a scaffold of elastic perforating fibers (PFs). The greater number of osteocyte lacunae compounded by known perilacunar strain amplification and the presence of elastic PFs are expected to be associated with greater bone tissue strain in primary than in secondary bone and thus with strain gradients at the PSBI. Strain gradients may lead local osteocytes to originate resorption canals and to promote transverse drift of the resorption front into lower-strain secondary bone, thus creating giant asymmetrical drifting osteonal canals that remodel primary to secondary bone. PFs extended from muscle fibers through periosteum and primary bone to the PSBI, where they were resorbed by origination of drifting canals. Growth modeling by periosteal osteoblasts proceeds in the gaps between PFs. Through the direct connection between muscle and the PSBI via PFs, muscle forces may influence not only modeling by raising strain but also remodeling of primary to secondary bone by increasing strain gradients at the PSBI. With reduction in primary bone width after the mid-teens, numbers of drifting canals and porosity declined. Differences in microstructural properties between primary and secondary bone are expected to generate strain gradients at the PSBI that contribute to site, transverse drift, asymmetry and large size of drifting canals, and, hence, to cortical porosity. Cortical porosity in children is a physiological feature of bone growth in width. Advisability of therapeutic intervention remains to be defined.