Canalicular Junctions in the Osteocyte Lacuno-Canalicular Network of Cortical Bone

Human cortical bone intrinsic permeability distribution based on 3D canalicular morphology

Bone permeability is a key parameter that drives osteocyte-based mechanobiological modelling and remodelling. While previous experimental and numerical studies have estimated bone permeability based on the morphology of the lacuno-canalicular network, these studies often relied on simplified geometries. In the current study, bone permeability was characterized using more realistic canalicular geometry for the morphological data. Bone samples harvested from 27 human femoral bones were investigated using synchrotron radiation-based nano-computed tomography with a voxel size of 100 nm. After segmenting the canaliculi and lacunae, each canaliculus was investigated individually by applying a distance map and watershed algorithms. Bone permeability based on canalicular morphology was then assessed using the Kozeny relation, which defines the permeability of a porous medium with capillary-like pores. An averaged intrinsic permeability value of 8.8 10-18 m2 was obtained. It should be noted that this study considered an empty canalicular network, however in vivo, both cellular and peri-cellular matrices decrease space for interstitial fluid flow and thus permeability. Furthermore, a voxel size of 100 nm does not allow for the detection of smaller canaliculi, which may also modify average permeability. With the current data set and the analytic process applied, the results showed a heterogeneous permeability distribution within bone tissue, both when comparing osteonal and interstitial tissues and within an individual osteon. A difference was observed between male and female samples, and permeability appeared to significantly decrease with age. Finally, a significant correlation was found between permeability and canalicular length density, defined as canalicular length per unit bone volume. This study proposes a new form of the Kozeny law to express bone canalicular permeability as a proportional relationship with the canalicular length density. Importantly, this parameter can be directly quantified through confocal fluorescence microcopy, which is more convenient than synchrotron radiation-based nano-computed tomography. In conclusion, the current study confirms that confocal microscopy can be serve as a reliable tool to estimate bone permeability. However, the permeability values calculated here are solely based on canalicular morphology and do not consider cellular and peri-cellular intra-canalicular features.

A systematic study of the effect of measurement parameters on determination of osteocyte lacunar properties using laboratory X-ray micro-computed tomography

Accurate 3D characterization of osteocyte lacunae is important when investigating the role of osteocytes under various physiological and pathological conditions but remains a challenge. With the continued development of laboratory X-ray micro-computed tomography, an increasing number of studies employ these techniques beyond traditional bone morphometry to quantify osteocyte lacunae. However, there is a lack of knowledge on the effect of measurement parameters on the image quality and resolution and in turn the osteocyte lacunar quantification. Herein, we have examined the interplay between scan parameters and the resultant lacunar quantification in terms of lacunar size, shape, and density by comparison with a synchrotron benchmark dataset. We summarize our conclusions in a guide for use of μ-CT for osteocyte lacunar quantification: (1) Identification of the measurement requirements to address the research questions. (2) Collection and preparation of suitable sample(s) that fulfills these requirements. (3) Experimental considerations including determination of the required voxel size, in turn dictating the maximum FOV and by extension the maximum size of the sample(s). The experimental parameters chosen should ensure optimal image contrast, sufficient signal to noise, angular sampling etc. Usually, it is advisable to measure as well as possible within the limits of time, budget, data storage and analysis capabilities. (4) Data analysis and reporting of the results, including visual examination of the data at multiple steps in the analysis, to ensure correct feature identification and suitable reporting approaches. (5) Cross study comparisons, which may be unsuitable if the experimental conditions and analysis strategies are not comparable.

Human bone ultrastructure in 3D: Multimodal correlative study combining nanoscale X-ray computed tomography and quantitative polarized Raman spectroscopy

Unique mechanical properties of cortical bone are defined by the arrangement and ratio of its organic and inorganic constituents. This arrangement can be influenced by ageing and disease, urging the understanding of normal and deviant morphological patterns down to the nanoscale level, as much as the exploration of techniques able to grant that knowledge. Here, the ultrastructure and composition of seven samples taken from the femoral neck cortical bone of a single donor (52 y.o. female, no metabolic bone disease) is assessed with emerging characterization techniques. Laboratory-based nanoscale X-ray computed tomography providing ∼50 nm spatial resolution at (16 nm)3 voxel size resolves not only the lacuno-canalicular network but also the mineral ellipsoids associated with mineralized collagen fibrils (MCF). Site-matching 3D data with quantitative polarized Raman spectroscopy provides, in turn, complementary information on relative mineral and organic composition, while both techniques allow to quantify the MCF orientation. Bone matrix composition and lacuna-canalicular network organization are shown to vary between the osteonal and interstitial zones. Both plywood and gradual oscillating motifs of bone lamellation are observed, in line with existing theories. By combining these two methods, future studies can concentrate on other bone ultrastructural units of interest like interlamellar and cement interfaces, the structure of MCF around lacunae and near Haversian channels, as well as the influence of metabolic diseases on bone ultrastructure. STATEMENT OF SIGNIFICANCE: This study provides new insights into bone hierarchical organization, revealing local composition and lacuno-canalicular network organization within osteonal and interstitial bone zones, as well as their mineralized collagen fiber (MCF) orientation within the lamella. Synchrotron-like resolution was achieved on a laboratory-based nano-CT by exposing the volumes of interest from the bulk sample and applying machine learning segmentation algorithms. Site-matched analysis with quantitative Polarized Raman spectroscopy (qPRS) provided indirect access to relative mineral and organic composition variations and local MCF out-of-plane angle, with good agreement between the two methods. The proposed correlative experiment workflow greatly facilitates the characterization of bone ultrastructure and can be applied to other fields dealing with ordered hierarchical materials of similar feature sizes.

The Lacunocanalicular Network is Denser in C57BL/6 Compared to BALB/c Mice

The lacunocanalicular network (LCN) is an intricate arrangement of cavities (lacunae) and channels (canaliculi), which permeates the mineralized bone matrix. In its porosity, the LCN accommodates the cell network of osteocytes. These two nested networks are attributed a variety of essential functions including transport, signaling, and mechanosensitivity due to load-induced fluid flow through the LCN. For a more quantitative assessment of the networks' function, the three-dimensional architecture has to be known. For this reason, we aimed (i) to quantitatively characterize spatial heterogeneities of the LCN in whole mouse tibial cross-sections of BALB/c mice and (ii) to analyze differences in LCN architecture by comparison with another commonly used inbred mouse strain, the C57BL/6 mouse. Both tibiae of five BALB/c mice (female, 26-week-old) were stained using rhodamine 6G and whole tibiae cross-sections were imaged using confocal laser scanning microscopy. Using image analysis, the LCN was quantified in terms of density and connectivity and lacunar parameters, such as lacunar degree, volume, and shape. In the same tibial cross-sections, the calcium content was measured using quantitative backscattered electron imaging (qBEI). A structural analysis of the LCN properties showed that spatially denser parts of the LCN are mainly due to a higher density of branching points in the network. While a high intra-individual variability of network density was detected within the cortex, the inter-individual variability between different mice was low. In comparison to C57BL/6J mice, BALB/c mice showed a distinct lower canalicular density. This reduced network was already detectable on a local network level with fewer canaliculi emanating from lacunae. Spatial correlation with qBEI images demonstrated that bone modeling resulted in disruptions in the network architecture. The spatial heterogeneity and differences in density of the LCN likely affects the fluid flow within the network and therefore bone's mechanoresponse to loading.

A comparative study on the effects of biodegradable high-purity magnesium screw and polymer screw for fixation in epiphyseal trabecular bone

With mechanical strength close to cortical bone, biodegradable and osteopromotive properties, magnesium (Mg)-based implants are promising biomaterials for orthopedic applications. However, during the degradation of such implants, there are still concerns on the potential adverse effects such as formation of cavities, osteolytic phenomena and chronic inflammation. Therefore, to transform Mg-based implants into clinical practice, the present study evaluated the local effects of high-purity Mg screws (HP-Mg, 99.99 wt%) by comparing with clinically approved polylactic acid (PLA) screws in epiphyseal trabecular bone of rabbits. After implantation of screws at the rabbit distal femur, bone microstructural, histomorphometric and biomechanical properties were measured at various time points (weeks 4, 8 and 16) using micro-CT, histology and histomorphometry, micro-indentation and scanning electron microscope. HP-Mg screws promoted peri-implant bone ingrowth with higher bone mass (BV/TV at week 4: 0.189 ± 0.022 in PLA group versus 0.313 ± 0.053 in Mg group), higher biomechanical properties (hardness at week 4: 35.045 ± 1.000 HV in PLA group versus 51.975 ± 2.565 HV in Mg group), more mature osteocyte LCN architecture, accelerated bone remodeling process and alleviated immunoreactive score (IRS of Ram11 at week 4: 5.8 ± 0.712 in PLA group versus 3.75 ± 0.866 in Mg group) as compared to PLA screws. Furthermore, we conducted finite element analysis to validate the superiority of HP-Mg screws as orthopedic implants by demonstrating reduced stress concentration and uniform stress distribution around the bone tunnel, which led to lower risks of trabecular microfractures. In conclusion, HP-Mg screws demonstrated greater osteogenic bioactivity and limited inflammatory response compared to PLA screws in the epiphyseal trabecular bone of rabbits. Our findings have paved a promising way for the clinical application of Mg-based implants.

3D quantification of the lacunocanalicular network on human femoral diaphysis through synchrotron radiation-based nanoCT

Osteocytes are the major actors in bone mechanobiology. Within bone matrix, they are trapped close together in a submicrometric interconnected network: the lacunocanalicular network (LCN). The interstitial fluid circulating within the LCN transmits the mechanical information to the osteocytes that convert it into a biochemical signal. Understanding the interstitial fluid dynamics is necessary to better understand the bone mechanobiology. Due to the submicrometric dimensions of the LCN, making it difficult to experimentally investigate fluid dynamics, numerical models appear as a relevant tool for such investigation. To develop such models, there is a need for geometrical and morphological data on the human LCN. This study aims at providing morphological data on the human LCN from measurement of 27 human femoral diaphysis bone samples using synchrotron radiation nano-computed tomography with an isotropic voxel size of 100 nm. Except from the canalicular diameter, the canalicular morphological parameters presented a high variability within one sample. Some differences in terms of both lacunar and canalicular morphology were observed between the male and female populations. But it has to be highlighted that all the canaliculi cannot be detected with a voxel size of 100 nm. Hence, in the current study, only a specific population of large canaliculi that could be characterize. Still, to the authors knowledge, this is the first time such a data set was introduced to the community. Further processing will be achieved in order to provide new insight on the LCN permeability.

Crossing length scales: X-ray approaches to studying the structure of biological materials

Biological materials have outstanding properties. With ease, challenging mechanical, optical or electrical properties are realised from comparatively `humble' building blocks. The key strategy to realise these properties is through extensive hierarchical structuring of the material from the millimetre to the nanometre scale in 3D. Though hierarchical structuring in biological materials has long been recognized, the 3D characterization of such structures remains a challenge. To understand the behaviour of materials, multimodal and multi-scale characterization approaches are needed. In this review, we outline current X-ray analysis approaches using the structures of bone and shells as examples. We show how recent advances have aided our understanding of hierarchical structures and their functions, and how these could be exploited for future research directions. We also discuss current roadblocks including radiation damage, data quantity and sample preparation, as well as strategies to address them.

2D vs. 3D Evaluation of Osteocyte Lacunae - Methodological Approaches, Recommended Parameters, and Challenges: A Narrative Review by the European Calcified Tissue Society (ECTS)

PURPOSE OF REVIEW

Quantification of the morphology of osteocyte lacunae has become a powerful tool to investigate bone metabolism, pathologies and aging. This review will provide a brief overview of 2D and 3D imaging methods for the determination of lacunar shape, orientation, density, and volume. Deviations between 2D-based and 3D-based lacunar volume estimations are often not sufficiently addressed and may give rise to contradictory findings. Thus, the systematic error arising from 2D-based estimations of lacunar volume will be discussed, and an alternative calculation proposed. Further, standardized morphological parameters and best practices for sampling and segmentation are suggested.

RECENT FINDINGS

We quantified the errors in reported estimation methods of lacunar volume based on 2D cross-sections, which increase with variations in lacunar orientation and histological cutting plane. The estimations of lacunar volume based on common practice in 2D imaging methods resulted in an underestimation of lacunar volume of up to 85% compared to actual lacunar volume in an artificial dataset. For a representative estimation of lacunar size and morphology based on 2D images, at least 400 lacunae should be assessed per sample.

Spatial variations in the osteocyte lacuno-canalicular network density and analysis of the connectomic parameters

Osteocyte lacuno-canalicular network (LCN) is comprised of micrometre-sized pores and submicrometric wide channels in bone. Accumulating evidence suggests multiple functions of this network in material transportation, mechanobiological signalling, mineral homeostasis and bone remodelling. Combining rhodamine staining and confocal laser scanning microscopy, the longitudinal cross-sections of six mouse tibiae were imaged, and the connectome of the network was quantified with a focus on the spatial heterogeneities of network density, connectivity and length of canaliculi. In-vivo loading and double calcein labelling on these tibiae allowed differentiating the newly formed bone from the pre-existing regions. The canalicular density of the murine cortical bone varied between 0.174 and 0.243 μm/μm3, and therefore is three times larger than the corresponding value for human femoral midshaft osteons. The spatial heterogeneity of the network was found distinctly more pronounced across the cortex than along the cortex. We found that in regions with a dense network, the LCN conserves its largely tree-like character, but increases the density by including shorter canaliculi. The current study on healthy mice should serve as a motivating starting point to study the connectome of genetically modified mice, including models of bone diseases and of reduced mechanoresponse.

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.

Multiscale morphological analysis of bone microarchitecture around Mg-10Gd implants

The utilization of biodegradable magnesium (Mg)-based implants for restoration of bone function following trauma represents a transformative approach in orthopaedic application. One such alloy, magnesium-10 weight percent gadolinium (Mg-10Gd), has been specifically developed to address the rapid degradation of Mg while enhancing its mechanical properties to promote bone healing. Previous studies have demonstrated that Mg-10Gd exhibits favorable osseointegration; however, it exhibits distinct ultrastructural adaptation in comparison to conventional implants like titanium (Ti). A crucial aspect that remains unexplored is the impact of Mg-10Gd degradation on the bone microarchitecture. To address this, we employed hierarchical three-dimensional imaging using synchrotron radiation in conjunction with image-based finite element modelling. By using the methods outlined, the vascular porosity, lacunar porosity and the lacunar-canaliculi network (LCN) morphology of bone around Mg-10Gd in comparison to Ti in a rat model from 4 weeks to 20 weeks post-implantation was investigated. Our investigation revealed that within our observation period, the degradation of Mg-10Gd implants was associated with significantly lower (p < 0.05) lacunar density in the surrounding bone, compared to Ti. Remarkably, the LCN morphology and the fluid flow analysis did not significantly differ for both implant types. In summary, a more pronounced lower lacunae distribution rather than their morphological changes was detected in the surrounding bone upon the degradation of Mg-10Gd implants. This implies potential disparities in bone remodelling rates when compared to Ti implants. Our findings shed light on the intricate relationship between Mg-10Gd degradation and bone microarchitecture, contributing to a deeper understanding of the implications for successful osseointegration.

Correlative study of liquid in human bone by 3D neutron microscopy and lab-based X-ray μCT

Liquid plays an important role in bone that has a complex 3D hierarchical pore structure. However, liquid (water) is difficult to discern from e.g. an organic matrix by X-ray imaging. Therefore, we use a correlative approach using both high resolution X-ray and neutron imaging. Human femoral bone with liquid adsorbed into some of the pores was imaged with both the Neutron Microscope at the ICON beamline, SINQ at PSI, and by lab-based μCT using 2.7 μm voxel size. Segmentation of the two datasets showed that, even though the liquid was clearly distinguishable in the neutron data and not in the X-ray data, it remained challenging to segment it from bone due to overlaps of peaks in the gray level histograms. In consequence, segmentations from X-ray and neutron data varied significantly. To address this issue, the segmented X-ray porosities was overlaid on the neutron data, making it possible to localize the liquid in the vascular porosities of the bone sample and use the neutron attenuation to identify it as H2O. The contrast in the neutron images was lowered slightly between the bone and the liquid compared to the bone and the air. This correlative study shows that the complementary use of X-rays and neutrons is very favorable, since H2O is very distinct in the neutron data, while D2O, H2O, and organic matter can barely be distinguished from air in the X-ray data.

Micron-resolution Imaging of Cortical Bone under 14 T Ultrahigh Magnetic Field

Compact, mineralized cortical bone tissues are often concealed on magnetic resonance (MR) images. Recent development of MR instruments and pulse techniques has yielded significant advances in acquiring anatomical and physiological information from cortical bone despite its poor 1 H signals. This work demonstrates the first MR research on cortical bones under an ultrahigh magnetic field of 14 T. The 1 H signals of different mammalian species exhibit multi-exponential decays of three characteristic T2 or T2 * values: 0.1-0.5 ms, 1-4 ms, and 4-8 ms. Systematic sample comparisons attribute these T2 /T2 * value ranges to collagen-bound water, pore water, and lipids, respectively. Ultrashort echo time (UTE) imaging under 14 T yielded spatial resolutions of 20-80 microns, which resolves the 3D anatomy of the Haversian canals. The T2 * relaxation characteristics further allow spatial classifications of collagen, pore water and lipids in human specimens. The study achieves a record of the spatial resolution for MR imaging in bone and shows that ultrahigh-field MR has the unique ability to differentiate the soft and organic compartments in bone tissues.

Detailing the influence of PEO-coated biodegradable Mg-based implants on the lacuno-canalicular network in sheep bone: A pilot study

An increasing prevalence of bone-related injuries and aging geriatric populations continue to drive the orthopaedic implant market. A hierarchical analysis of bone remodelling after material implantation is necessary to better understand the relationship between implant and bone. Osteocytes, which are housed and communicate through the lacuno-canalicular network (LCN), are integral to bone health and remodelling processes. Therefore, it is essential to examine the framework of the LCN in response to implant materials or surface treatments. Biodegradable materials offer an alternative solution to permanent implants, which may require revision or removal surgeries. Magnesium alloys have resurfaced as promising materials due to their bone-like properties and safe degradation in vivo. To further tailor their degradation capabilities, surface treatments such as plasma electrolytic oxidation (PEO) have demonstrated to slow degradation. For the first time, the influence of a biodegradable material on the LCN is investigated by means of non-destructive 3D imaging. In this pilot study, we hypothesize noticeable variations in the LCN caused by altered chemical stimuli introduced by the PEO-coating. Utilising synchrotron-based transmission X-ray microscopy, we have characterised morphological LCN differences around uncoated and PEO-coated WE43 screws implanted into sheep bone. Bone specimens were explanted after 4, 8, and 12 weeks and regions near the implant surface were prepared for imaging. Findings from this investigation indicate that the slower degradation of PEO-coated WE43 induces healthier lacunar shapes within the LCN. However, the stimuli perceived by the uncoated material with higher degradation rates induces a greater connected LCN better prepared for bone disturbance.

Synchrotron X-ray Studies of the Structural and Functional Hierarchies in Mineralised Human Dental Enamel: A State-of-the-Art Review

Hard dental tissues possess a complex hierarchical structure that is particularly evident in enamel, the most mineralised substance in the human body. Its complex and interlinked organisation at the Ångstrom (crystal lattice), nano-, micro-, and macro-scales is the result of evolutionary optimisation for mechanical and functional performance: hardness and stiffness, fracture toughness, thermal, and chemical resistance. Understanding the physical-chemical-structural relationships at each scale requires the application of appropriately sensitive and resolving probes. Synchrotron X-ray techniques offer the possibility to progress significantly beyond the capabilities of conventional laboratory instruments, i.e., X-ray diffractometers, and electron and atomic force microscopes. The last few decades have witnessed the accumulation of results obtained from X-ray scattering (diffraction), spectroscopy (including polarisation analysis), and imaging (including ptychography and tomography). The current article presents a multi-disciplinary review of nearly 40 years of discoveries and advancements, primarily pertaining to the study of enamel and its demineralisation (caries), but also linked to the investigations of other mineralised tissues such as dentine, bone, etc. The modelling approaches informed by these observations are also overviewed. The strategic aim of the present review was to identify and evaluate prospective avenues for analysing dental tissues and developing treatments and prophylaxis for improved dental health.

Insights into the Molecular and Hormonal Regulation of Complications of X-Linked Hypophosphatemia

X-linked hypophosphatemia (XLH) is characterized by mutations in the PHEX gene, leading to elevated serum levels of FGF23, decreased production of 1,25 dihydroxyvitamin D3 (1,25D), and hypophosphatemia. Those affected with XLH manifest impaired growth and skeletal and dentoalveolar mineralization as well as increased mineralization of the tendon–bone attachment site (enthesopathy), all of which lead to decreased quality of life. Many molecular and murine studies have detailed the role of mineral ions and hormones in regulating complications of XLH, including how they modulate growth and growth plate maturation, bone mineralization and structure, osteocyte-mediated mineral matrix resorption and canalicular organization, and enthesopathy development. While these studies have provided insight into the molecular underpinnings of these skeletal processes, current therapies available for XLH do not fully prevent or treat these complications. Therefore, further investigations are needed to determine the molecular pathophysiology underlying the complications of XLH.

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).

Bone hierarchical structure: spatial variation across length scales

Bone is a complex hierarchical biomineralized material, which is special amongst biominerals because it is replete with cells, namely, osteocytes. While bone has been scrutinized for centuries, many questions remain open and new research hints that the ultrastructure of bone, encompassing both the bone matrix itself and the embedded cell network, is much more heterogeneous than hitherto realized. A number of these new findings have been made thanks to the enormous developments in X-ray imaging that have occurred in recent decades, and there is promise that they will also allow many of the remaining open questions to be addressed. X-ray absorption or phase imaging affords high three-dimensional (3D) resolution and allows traversing the length scales of bone all the way down to the fine details of the lacuno-canalicular network housing the osteocytes. Multimodal X-ray imaging provides combined information covering both the length scales defined by the size of the measured volume and tomographic resolution, as well as those probed by the signal that is measured. In X-ray diffraction computed tomography (XRD-CT), for example, diffraction signals can be reconstructed tomographically, which offers detailed information about the spatial variations in the crystallographic properties of the bone biomineral. Orientational information can be obtained by tensor tomography. The combination of both small-angle X-ray scattering (SAXS) and wide-angle X-ray scattering (WAXS) tensor tomography gives information on the orientation of bone nanostructure and crystals, respectively. These new technical developments promise that great strides towards understanding bone structure can be expected in the near future. In this review, recent findings that have resulted from X-ray imaging are highlighted and speculation is given on what can be expected to follow.

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.

Evaluation of imaging setups for quantitative phase contrast nanoCT of mineralized biomaterials

X-ray nano-tomography with phase contrast (nanoCT) using synchrotron radiation is a powerful tool to non-destructively investigate 3D material properties at the nanoscale. In large bone lesions, such as severe bone fractures, bone cancer or other diseases, bone grafts substituting the lost bone might be necessary. Such grafts can be of biological origin or be composed of a synthetic bone substitute. The long-term functioning of artificial bone substitutes depends on many factors. Synchrotron nanoCT imaging has great potential to contribute to further the understanding of integration of implants into bone tissue by imaging the spatial interaction between bone tissue and implant, and by accessing the interface between implant material and bone tissue. With this aim, a methodology for evaluating the image quality is presented for in-line phase contrast nanoCT images of bone scaffold samples. A PMMA-embedded tricalcium phosphate scaffold was used with both a closed and an open porosity structure and bone ingrowths as a representative system of three known materials. Parameters such as spatial resolution and signal-to-noise ratio were extracted and used to explore and quantitatively compare the effects of implementation choices in the imaging setup, such as camera technology and imaging energy, on the resulting image quality. Increasing the X-ray energy from 17.5 keV to 29.6 keV leads to a notable improvement in image quality regardless of the camera technology used, with the two tested camera setups performing at a comparable level when the recorded intensity was kept constant.

Opportunities for biomineralization research using multiscale computed X-ray tomography as exemplified by bone imaging

Biominerals typically have complex hierarchical structures traversing many length scales. This makes their structural characterization complicated, since it requires 3D techniques that can probe full specimens at down to nanometer-resolution, a combination that is difficult - if not impossible - to achieve simultaneously. One challenging example is bone, a mineralized tissue with a highly complex architecture that is replete with a network of cells. X-ray computed tomography techniques enable multiscale structural characterization through the combination of various equipment and emerge as promising tools for characterizing biominerals. Using bone as an example, we discuss how combining different X-ray imaging instruments allow characterizing bone structures from the nano- to the organ-scale. In particular, we compare and contrast human and rodent bone, emphasize the importance of the osteocyte lacuno-canalicular network in bone, and finally illustrate how combining synchrotron X-ray imaging with laboratory instrumentation for computed tomography is especially helpful for multiscale characterization of biominerals.

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.

Disrupted osteocyte connectivity and pericellular fluid flow in bone with aging and defective TGF-β signaling

Skeletal fragility in the elderly does not simply result from a loss of bone mass. However, the mechanisms underlying the concurrent decline in bone mass, quality, and mechanosensitivity with age remain unclear. The important role of osteocytes in these processes and the age-related degeneration of the intricate lacunocanalicular network (LCN) in which osteocytes reside point to a primary role for osteocytes in bone aging. Since LCN complexity severely limits experimental dissection of these mechanisms in vivo, we used two in silico approaches to test the hypothesis that LCN degeneration, due to aging or an osteocyte-intrinsic defect in transforming growth factor beta (TGF-β) signaling (TβRIIocy-/-), is sufficient to compromise essential osteocyte responsibilities of mass transport and exposure to mechanical stimuli. Using reconstructed confocal images of bone with fluorescently labeled osteocytes, we found that osteocytes from aged and TβRIIocy-/- mice had 33 to 45% fewer, and more tortuous, canaliculi. Connectomic network analysis revealed that diminished canalicular density is sufficient to impair diffusion even with intact osteocyte numbers and overall LCN architecture. Computational fluid dynamics predicts that the corresponding drop in shear stress experienced by aged or TβRIIocy-/- osteocytes is highly sensitive to canalicular surface area but not tortuosity. Simulated expansion of the osteocyte pericellular space to mimic osteocyte perilacunar/canalicular remodeling restored predicted shear stress for aged osteocytes to young levels. Overall, these models show how loss of LCN volume through LCN pruning may lead to impaired fluid dynamics and osteocyte exposure to mechanostimulation. Furthermore, osteocytes emerge as targets of age-related therapeutic efforts to restore bone health and function.

Hierarchical Nature of Nanoscale Porosity in Bone Revealed by Positron Annihilation Lifetime Spectroscopy

Bone is a hierarchical material primarily composed of collagen, water, and mineral that is organized into discrete molecular, nano-, micro-, and macroscale structural components. In contrast to the structural knowledge of the collagen and mineral domains, the nanoscale porosity of bone is poorly understood. In this study, we introduce a well-established pore characterization technique, positron annihilation lifetime spectroscopy (PALS), to probe the nanoscale size and distribution of each component domain by analyzing pore sizes inherent to hydrated bone together with pores generated by successive removal of water and then organic matrix (including collagen and noncollagenous proteins) from samples of cortical bovine femur. Combining the PALS results with simulated pore size distribution (PSD) results from collagen molecule and microfibril structure, we identify pores with diameter of 0.6 nm that suggest porosity within the collagen molecule regardless of the presence of mineral and water. We find that water occupies three larger domain size regions with nominal mean diameters of 1.1, 1.9, and 4.0 nm-spaces that are hypothesized to associate with intercollagen molecular spaces, terminal segments (d-spacing) within collagen microfibrils, and interface spacing between collagen and mineral structure, respectively. Subsequent removal of the organic matrix determines a structural pore size of 5-6 nm for deproteinized bone-suggesting the average spacing between mineral lamella. An independent method to deduce the average mineral spacing from specific surface area (SSA) measurements of the deproteinized sample is presented and compared with the PALS results. Together, the combined PALS and SSA results set a range on the mean mineral lamella thickness of 4-8 nm.

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.

Bone metabolism and evolutionary origin of osteocytes: Novel application of FIB-SEM tomography

Lacunae and canaliculi spaces of osteocytes are remarkably well preserved in fossilized bone and serve as an established proxy for bone cells. The earliest bone in the fossil record is acellular (anosteocytic), followed by cellular (osteocytic) bone in the jawless relatives of jawed vertebrates, the osteostracans, about 400 million years ago. Virtually nothing is known about the physiological pressures that would have initially favored osteocytic over anosteocytic bone. We apply focused ion beam-scanning electron microscopy tomography combined with machine learning for cell detection and segmentation to image fossil cell spaces. Novel three-dimensional high-resolution images reveal areas of low density around osteocyte lacunae and their canaliculi in osteostracan bone. This provides evidence for demineralization that would have occurred in vivo as part of osteocytic osteolysis, a mechanism of mineral homeostasis, supporting the hypothesis that a physiological demand for phosphorus was the principal driver in the initial evolution of osteocytic bone.

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.

Heterogeneity of the osteocyte lacuno-canalicular network architecture and material characteristics across different tissue types in healing bone

Various tissue types, including fibrous connective tissue, bone marrow, cartilage, woven and lamellar bone, coexist in healing bone. Similar to most bone tissue type, healing bone contains a lacuno-canalicular network (LCN) housing osteocytes. These cells are known to orchestrate bone remodeling in healthy bone by sensing mechanical strains and translating them into biochemical signals. The structure of the LCN is hypothesized to influence mineralization processes. Hence, the aim of the present study was to visualize and match spatial variations in the LCN topology with mineral characteristics, within and at the interfaces of the different tissue types that comprise healing bone. We applied a correlative multi-method approach to visualize the LCN architecture and quantify mineral particle size and orientation within healing femoral bone in a mouse osteotomy model (26 weeks old C57BL/6 mice). This approach revealed structural differences across several length scales during endochondral ossification within the following regions: calcified cartilage, bony callus, cortical bone and a transition zone between the cortical and callus region analyzed 21 days after the osteotomy. In this transition zone, we observed a continuous convergence of mineral characteristics and osteocyte lacunae shape as well as discontinuities in the lacunae volume and LCN connectivity. The bony callus exhibits a 34% higher lacunae number density and 40% larger lacunar volume compared to cortical bone. The presented correlations between LCN architecture and mineral characteristics improves our understanding of how bone develops during healing and may indicate a contribution of osteocytes to bone (re)modeling.

Ellipsoidal mesoscale mineralization pattern in human cortical bone revealed in 3D by plasma focused ion beam serial sectioning

Visualizing bone mineralization and collagen fibril organization at intermediate scales between the nanometer and the hundreds of microns range, is still an important challenge. Similarly, visualizing cellular components which locally affect the tissue structure requires a precision of a few tens of nanometers at maximum while spanning several tens of micrometers. In the last decade, gallium focused ion beam (FIB) equipped with a scanning electron microscope (SEM) proved to be an extremely valuable structural tool to meet those ends. In this study, we assess the capability of a recent plasma FIB-SEM technology which provides a potential increase in measurement speed over gallium FIB-SEM, thus paving the way to larger volume analysis. Nanometer-scale layers of demineralized and mineralized unstained human femoral lamellar bone were sequentially sectioned over volumes of 6-16,000 μm³. Analysis of mineralized tissue revealed prolate ellipsoidal mineral clusters measuring approximately 1.1 µm in length by 700 nm at their maximum diameter. Those features, suggested by others in high resolution studies, appear here as a ubiquitous motif in mineralized lamellar bone over thousands of microns cubed, suggesting a heterogeneous and yet regular pattern of mineral deposition past the single collagen fibril level. This large scale view retained sufficient resolution to visualize the collagen fibrils while also partly visualizing the lacuno-canalicular network in three-dimensions. These findings are strong evidence for suitability of PFIB as a bone analysis tool and the need to revisit bone mineralization over multi-length scales with mineralized tissue.

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.

Bone Biomineral Properties Vary across Human Osteonal Bone

The biomineralization of bone remains a puzzle. During Haversian remodeling in the dense human cortical bone, osteoclasts excavate a tunnel that is then filled in by osteoblasts with layers of bone of varying fibril orientations, resulting in a lamellar motif. Such bone represents an excellent possibility to increase our understanding of bone as a material as well as bone biomineralization by studying spatio/temporal variations in the biomineral across an osteon. To this end, fluorescence computed tomography and diffraction scattering computed tomography with sub-micrometer resolution is applied to obtain position resolved fluorescence spectra and diffraction patterns in a 3D volume. The microstructural properties of the apatite biomineral are not homogeneous but depend critically on the time point at which it was laid down. This indicates that the nature of bone biomineral is highly dependent on the microenvironment during bone formation and remodeling.