Quantitative 3D analysis of the canal network in cortical bone by micro-computed tomography

Vascular study of decellularized porcine long bones: Characterization of a tissue engineering model

INTRODUCTION

Massive bone allografts enable the reconstruction of critical bone defects in numerous conditions (e.g. tumoral, infection or trauma). Unfortunately, their biological integration remains insufficient and the reconstruction may suffer from several postoperative complications. Perfusion-decellularization emerges as a tissue engineering potential solution to enhance osseointegration. Therefore, an intrinsic vascular study of this novel tissue engineering tool becomes essential to understand its efficacy and applicability.

MATERIAL AND METHODS

32 porcine long bones (humeri and femurs) were used to assess the quality of their vascular network prior and after undergoing a perfusion-decellularization protocol. 12 paired bones were used to assess the vascular matrix prior (N = 6) and after our protocol (N = 6) by immunohistochemistry. Collagen IV, Von Willebrand factor and CD31 were targeted then quantified. The medullary macroscopic vascular network was evaluated with 12 bones: 6 were decellularized and the other 6 were, as control, not treated. All 12 underwent a contrast-agent injection through the nutrient artery prior an angio CT-scan acquisition. The images were processed and the length of medullary vessels filled with contrast agent were measured on angiographic cT images obtained in control and decellularized bones by 4 independent observers to evaluate the vascular network preservation. The microscopic cortical vascular network was evaluated on 8 bones: 4 control and 4 decellularized. After injection of gelatinous fluorochrome mixture (calcein green), non-decalcified fluoroscopic microscopy was performed in order to assess the perfusion quality of cortical vascular lacunae.

RESULTS

The continuity of the microscopic vascular network was assessed with Collagen IV immunohistochemistry (p-value = 0.805) while the decellularization quality was observed through CD31 and Von Willebrand factor immunohistochemistry (p-values <0.001). The macroscopic vascular network was severely impaired after perfusion-decellularization; nutrient arteries were still patent but the amount of medullary vascular channels measured was significantly higher in the control group compared to the decellularized group (p-value <0.001). On average, the observers show good agreement on these results, except in the decellularized group where more inter-observer discrepancies were observed. The microscopic vascular network was observed with green fluoroscopic signal in almost every canals and lacunae of the bone cortices, in three different bone locations (proximal metaphysis, diaphysis and distal metaphysis).

CONCLUSION

Despite the aggressiveness of the decellularization protocol on medullary vessels, total porcine long bones decellularized by perfusion retain an acellular cortical microvascular network. By injection through the intact nutrient arteries, this latter vascular network can still be used as a total bone infusion access for bone tissue engineering in order to enhance massive bone allografts prior implantation.

Understanding basic multicellular unit activity in cortical bone through 3D morphological analysis: New methods to define zones of the remodeling space

The activity of basic multicellular units (BMU) in cortical bone is classically described as a sequential order of events- resorption, reversal and formation. This simplified portrayal of the remodeling process is pervasive despite the reported variability in remodeling space morphology. These variations may reflect meaningful nuances in BMU activity but methods to quantify 3D remodeling space morphology within the context of the cellular activity are currently lacking. This study developed new techniques to define zones of BMU activity based on the 3D morphology of remodeling spaces in rabbit cortical bone and integrated morphological data with the BMU longitudinal erosion rate (LER) to elucidate the spatial-temporal coordination of BMUs and estimate mineral apposition rate (MAR). The tibiae of New Zealand white rabbits (n = 5) were imaged in vivo using synchrotron radiation and two weeks later ex vivo with desktop microCT. The in vivo and ex vivo datasets were co-registered, and 27 remodeling spaces were identified at both timepoints. A radial profile representing the 3D morphology was the platform for partitioning the remodeling spaces into resorption, reversal and formation zones. Manual, automated and semi-automated partitioning approaches were compared, and the zone-segmentations were used to calculate the length, change in radius and slope of each zone. The manual approach most accurately defined the zones of idealized remodeling spaces with known dimensions (relative error = 0.9-9.2 %) while the semi-automated method reliably defined the zones in rabbit remodeling spaces (ICC = 0.85-1.00). Combining LER and the manually derived zone dimensions indicated that a BMU passes through a cross-section in approximately 18.8 days with resorption, reversal and formation taking 4.1, 2.2, and 12.5 days, respectively. MAR estimated by the 3D analysis was not significantly different than that determined with classic histomorphometry (p = 0.48). These techniques have the potential to assess dynamic parameters of bone resorption and formation, eliminate the need for fluorochrome labeling and provide a more comprehensive perspective of the remodeling process.

Characterization of rat vertebrae cortical bone microstructures using confocal Raman microscopy combined to tomography and electron microscopy

BACKGROUND

The rat vertebrae is a good model to study bone regeneration after implantation of biomaterials used to treat bone loss, a major problem in oral and dental surgery. However, the precise characterization of bone microstructures in the rat vertebrae has not been reported. Therefore, the aim of this study was to achieve the complete analysis of such bone, at different scales, in order to have a clear model of healthy bone for comparison with regenerated bone.

METHODS

In order to image the cortical bone of rat caudal vertebra, confocal Raman microscopy was combined with high resolution X-ray micro computed tomography (micro-CT), with scanning electron microscopy (SEM) using backscatter electron imaging and with more conventional histology coloration techniques. SEM and Raman microscopy were done in various regions of the cortical bone corresponding to external, middle and internal areas. The spongy bone was imaged in parallel. Micro-CT was performed on the whole vertebra to monitor the network of haversian canals in the cortical bone. Osteonic canals characteristics, and relative chemical composition were analysed in several regions of interest, in cortical and spongy bone. Five rats were included in this study.

RESULTS

On micro-CT images, differences in intensity were observed in the cortical bone, substantiated by SEM. Chemical analysis with Raman spectra confirmed the difference in composition between the different regions of the cortical and spongy bone. PCA and k-mean cluster analysis separated these groups, except for the external and middle cortical bone. Peak intensity ratio confirmed these results with a CO3 to ν2 PO4 ratio significantly different for the internal cortical bone. Grayscale images stack extracted from micro-CT showed that global architecture of cortical bone was characterized by a dense and complex network of haversian osteonic canals, starting from the surface towards the vertebrae center. The mean diameter of the canals was 18.4 µm (SD 8.6 µm) and the mean length was 450 µm (SD 152 µm). Finally, Raman reconstructed images of the lamellar bone showed an enlargement of the lamellar layer width, both in circumferential lamellar bone and around haversian canals.

CONCLUSIONS

Micro-CT and confocal Raman microscopy are good tools to complete classical analysis using optical and electron microscopy. The results and measurements presented in a rat model known for its small inter-individual differences provide the main characteristics of a mature bone. This study will allow the community working on this rat vertebrate model to have a set of characteristics, in particular on the structure of the haversian canals.

Direct Assessment of Rabbit Cortical Bone Basic Multicellular Unit Longitudinal Erosion Rate: A 4D Synchrotron-Based Approach

Cortical bone remodeling is carried out by basic multicellular units (BMUs), which couple resorption to formation. Although fluorochrome labeling has facilitated study of BMU formative parameters since the 1960s, some resorptive parameters, including the longitudinal erosion rate (LER), have remained beyond reach of direct measurement. Indeed, our only insights into this spatiotemporal parameter of BMU behavior come from classical studies that indirectly inferred LER. Here, we demonstrate a 4D in vivo method to directly measure LER through in-line phase contrast synchrotron imaging. The tibias of rabbits (n = 15) dosed daily with parathyroid hormone were first imaged in vivo (synchrotron micro-CT; day 15) and then ex vivo 14 days later (conventional micro-CT; day 29). Mean LER assessed by landmarking the co-registered scans was 23.69 ± 1.73 μm/d. This novel approach holds great promise for the direct study of the spatiotemporal coordination of bone remodeling, its role in diseases such as osteoporosis, as well as related treatments. © 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).

Therapeutic Targeting Notch2 Protects Bone Micro-Vasculatures from Methotrexate Chemotherapy-Induced Adverse Effects in Rats

Intensive cancer chemotherapy is well known to cause bone vasculature disfunction and damage, but the mechanism is poorly understood and there is a lack of treatment. Using a rat model of methotrexate (MTX) chemotherapy (five once-daily dosses at 0.75 mg/kg), this study investigated the roles of the Notch2 signalling pathway in MTX chemotherapy-induced bone micro-vasculature impairment. Gene expression, histological and micro-computed tomography (micro-CT) analyses revealed that MTX-induced micro-vasculature dilation and regression is associated with the induction of Notch2 activity in endothelial cells and increased production of inflammatory cytokine tumour necrosis factor alpha (TNFα) from osteoblasts (bone forming cells) and bone marrow cells. Blockade of Notch2 by a neutralising antibody ameliorated MTX adverse effects on bone micro-vasculature, both directly by supressing Notch2 signalling in endothelial cells and indirectly via reducing TNFα production. Furthermore, in vitro studies using rat bone marrow-derived endothelial cell revealed that MTX treatment induces Notch2/Hey1 pathway and negatively affects their ability in migration and tube formation, and Notch2 blockade can partially protect endothelial cell functions from MTX damage.

Local anisotropy in mineralized fibrocartilage and subchondral bone beneath the tendon-bone interface

The enthesis allows the insertion of tendon into bone thanks to several remarkable strategies. This complex and clinically relevant location often features a thin layer of fibrocartilage sandwiched between tendon and bone to cope with a highly heterogeneous mechanical environment. The main purpose of this study was to investigate whether mineralized fibrocartilage and bone close to the enthesis show distinctive three-dimensional microstructural features, possibly to enable load transfer from tendon to bone. As a model, the Achilles tendon-calcaneus bone system of adult rats was investigated with histology, backscattered electron imaging and micro-computed tomography. The microstructural porosity of bone and mineralized fibrocartilage in different locations including enthesis fibrocartilage, periosteal fibrocartilage and bone away from the enthesis was characterized. We showed that calcaneus bone presents a dedicated protrusion of low porosity where the tendon inserts. A spatially resolved analysis of the trabecular network suggests that such protrusion may promote force flow from the tendon to the plantar ligament, while partially relieving the trabecular bone from such a task. Focusing on the tuberosity, highly specific microstructural aspects were highlighted. Firstly, the interface between mineralized and unmineralized fibrocartilage showed the highest roughness at the tuberosity, possibly to increase failure resistance of a region carrying large stresses. Secondly, fibrochondrocyte lacunae inside mineralized fibrocartilage, in analogy with osteocyte lacunae in bone, had a predominant alignment at the enthesis and a rather random organization away from it. Finally, the network of subchondral channels inside the tuberosity was highly anisotropic when compared to contiguous regions. This dual anisotropy of subchondral channels and cell lacunae at the insertion may reflect the alignment of the underlying collagen network. Our findings suggest that the microstructure of fibrocartilage may be linked with the loading environment. Future studies should characterize those microstructural aspects in aged and or diseased conditions to elucidate the poorly understood role of bone and fibrocartilage in enthesis-related pathologies.

Quantitative analysis of bone microvasculature in a mouse model using the monogenic signal phase asymmetry and marker-controlled watershed

Thethree-dimensional (3D) imaging and quantitative analysis of bone microvasculature are important to describe angiogenesis involvement in bone metastatic processes. Here, we propose an algorithm based on marker-controlled watershed for the 3D segmentation of vessels and bone in mouse bone imaged with a contrast agent using synchrotron radiation micro-computed tomography (SR-μCT). Markers were generated using hysteresis thresholding and morphological filters, and the control surface was constructed using the monogenic signal phase asymmetry. The accuracy and robustness of the proposed method were evaluated on a series of synthetic volumes generated to mimic the vessel, bone and background structures. Different contrast between different structures, as well as different noise levels were considered. A series of multi-class synthetic volumes were segmented using the proposed method, and the overall segmentation quality was evaluated using the Matthews correlation coefficient (MCC) by comparing to the ground truth. Additionally, we evaluated the segmentation of thin structures under various levels of Gaussian noise. The simulation study indicated that the algorithm was performant in multi-class segmentation with different contrast, noise, and thickness. The algorithm was applied to images of bone from a mouse model of breast cancer bone metastasis acquired using SR-μCT. The segmentation quality was evaluated using the Dice coefficient and the MCC by comparing to manual segmentation. The proposed method performed better than hysteresis thresholding and marker-controlled watershed using the magnitude of the gradient as control surface. Several quantitative parameters on bone and vessels were extracted, including bone volume fraction (BV/TV), vessel volume fraction (VV/TV) and the mean vessel thickness (VTh). The bone volume fraction (BV/TV) was significantly lower in the metastatic group compared to the healthy group. This demonstrated the effectiveness of the algorithm for the study of bone and vessel microstructures in mouse model.

Stressed volume estimated by finite element analysis predicts the fatigue life of human cortical bone: The role of vascular canals as stress concentrators

The fatigue life of cortical bone can vary several orders of magnitude, even in identical loading conditions. A portion of this variability is likely related to intracortical microarchitecture and the role of vascular canals as stress concentrators. The size, spatial distribution, and density of canals determine the peak magnitude and volume of stress concentrations. This study utilized a combination of experimental fatigue testing and image-based finite element (FE) analysis to establish the relationship between the stressed volume (i.e., volume of bone above yield stress) associated with vascular canals and the fatigue life of cortical bone. Thirty-six cortical bone samples were prepared from human femora and tibiae from five donors. Samples were allocated to four loading groups, corresponding to stress ranges of 60, 70, 80, and 90 MPa, then cyclically loaded in zero-compression until fracture. Porosity, canal diameter, canal separation, and canal number for each sample was quantified using X-ray microscopy (XRM) after testing. FE models were created from XRM images and used to calculate the stressed volume. Stressed volume was a good predictor of fatigue life, accounting for 67% of the scatter in fatigue-life measurements. An increase in stressed volume was most strongly associated with higher levels of intracortical porosity and larger canal diameters. The findings from this study suggest that a large portion of the fatigue-life variance of cortical bone in zero-compression is driven by intracortical microarchitecture, and that fatigue failure may be predicted by quantifying the stress concentrations associated with vascular canals.

Intra-skeletal vascular density in a bipedal hopping macropod with implications for analyses of rib histology

Human ribs are thought to be less affected by mechanical strain at the microscopic level than limb bones, implying that rib remodelling better reflects bone physiological homeostasis. Here, we test the hypothesis that rib tissue will be well vascularized and thus enhance susceptibility to metabolic influence. An intra-skeletal comparison of bone vascular canal density was conducted using a macropod animal model adapted to bipedal habitual hopping. The right humerus, ulna, radius, femur, tibia, fibula, a mid-thoracic and upper-thoracic rib of an eastern grey kangaroo (Macropus giganteus) were sectioned at the midshaft, from which histological sections were prepared. Bone vascularity from a maximum of 12 mm² of sub-periosteal parallel-fibred and lamellar bone was recorded, resulting in a total of 2047 counted vessels. Vascular canal density data were corrected by cortical width, maximum length, and midshaft circumference robusticity indices computed for each bone. The fibula consistently had the highest vascular canal density, even when corrected for maximum length, cortical width and midshaft circumference robusticities. This was followed by the mid- and upper-thoracic ribs. Vascularity differences between bones were relatively consistent whether vascular canal density was controlled for by cortical width or midshaft circumference robusticities. Vascular canal density and robusticity indices were also positively and negatively correlated (p < 0.05). Results confirm that the ribs are well vascularized, which facilitates bone metabolic processes such as remodelling, but the fibula also appears to be a well vascularized bone. Future research investigating human bone metabolism will benefit from examining thoracic rib or fibula samples.

Towards a 2D cortical osseous tissue representation and generation at micro scale. A computational model for bone simulations

BACKGROUND AND OBJECTIVE

the acquisition of microscopic images of human bones is a complex and expensive process. Moreover, the objective of obtaining a large data bank with microscopic images in order to carry out massive studies or to train automatic generation algorithms is not an option. Consequently, most of the current work focuses on the analysis of small regions captured by a microscope. The aim is the development of a tool to represent bone tissue at microscopic levels which is suitable for performing physical simulations, as well as for the diagnosis of various diseases. This work includes the whole process from the digitization of a human bone to the generation of bone tissue in a determined area of the bone selected through a cutting plane.

METHODS

based on the anatomy of the bone structure, the parameters that allow the representation of the bone tissue at mesoscale level have been analyzed. Although the models are randomly generated, they are based on statistical parameters. The model generator is based on the analysis of images of bone tissue and its parameters, performing a representation of each of its relevant structures in a way that fulfils these parameters.

RESULTS

the tool is useful for the virtual generation of bone tissue that satisfies the main characteristics of the cortical bone. The models obtained have been favorably evaluated in two stages. In the first stage, a scientific group has examined a set of images, in which images of the models generated were mixed with images obtained through traditional methods. Then, the physical characteristics of the generated tissue have been compared with the morphology of the bone tissue.

CONCLUSIONS

the model generator allows us to perform precise simulations in order to obtain realistic images with physical characteristics in accordance with reality. It is necessary to emphasize that even though the most relevant structures are included, the proposed model generator can be expanded to include new parameters or elements, so that it can be adapted to new needs. It could even break down randomness and parameterize it completely in order to allow the recreation of the tissue conditions of other studies.

Development of limb bone laminarity in the homing pigeon (Columba livia)

Background

Birds show adaptations in limb bone shape that are associated with resisting locomotor loads. Whether comparable adaptations occur in the microstructure of avian cortical bone is less clear. One proposed microstructural adaptation is laminar bone in which the proportion of circumferentially-oriented vascular canals (i.e., laminarity) is large. Previous work on adult birds shows elevated laminarity in specific limb elements of some taxa, presumably to resist torsion-induced shear strain during locomotion. However, more recent analyses using improved measurements in adult birds and bats reveal lower laminarity than expected in bones associated with torsional loading. Even so, there may still be support for the resistance hypothesis if laminarity increases with growth and locomotor maturation.

Methods

Here, we tested that hypothesis using a growth series of 17 homing pigeons (15–563 g). Torsional rigidity and laminarity of limb bones were measured from histological sections sampled from midshaft. Ontogenetic trends in laminarity were assessed using principal component analysis to reduce dimensionality followed by beta regression with a logit link function.

Results

We found that torsional rigidity of limb bones increases disproportionately with growth, consistent with rapid structural compensation associated with locomotor maturation. However, laminarity decreases with maturity, weakening the hypothesis that high laminarity is a flight adaptation at least in the pigeon. Instead, the histological results suggest that low laminarity, specifically the relative proportion of longitudinal canals aligned with peak principal strains, may better reflect the loading history of a bone.

Association between intracortical microarchitecture and the compressive fatigue life of human bone: A pilot study

Many mechanical properties of cortical bone are largely governed by the underlying microarchitecture; however, the influence of microarchitecture on the fatigue life of bone is poorly understood. Furthermore, imaging-based studies investigating intracortical microarchitecture may expose bone samples to large doses of radiation that may compromise fatigue resistance. The purpose of this pilot study was to 1) investigate the relationship between intracortical microarchitecture and the fatigue life of human bone in compression and 2) examine the effects of synchrotron irradiation on fatigue life measurements. Cortical samples were prepared from the femoral and tibial shafts of three cadaveric donors. A subset of samples was imaged using synchrotron X-ray microCT to quantify microarchitecture, including porosity, canal diameter, lacunar density, lacunar volume, and lacunar orientation. A second group of control samples was not imaged and used only for mechanical testing. Fatigue life was quantified by cyclically loading both groups in zero-compression until failure. Increased porosity and larger canal diameter were both logarithmically related to a shorter fatigue life, whereas lacunar density demonstrated a positive linear relationship with fatigue life (r² = 45–73%, depending on measure). Irradiation from microCT scanning reduced fatigue life measurements by 91%, but relationships with microarchitecture measurements remained. Additional research is needed to support the findings of this pilot study and fully establish the relationship between intracortical microarchitecture and the compressive fatigue life of bone.

•

Increased porosity and larger canal diameters were associated with a shorter compressive fatigue life.

•

A higher lacunar density was related to a longer compressive fatigue life.

•

Irradiation from synchrotron X-ray microCT scanning reduced fatigue life by 91%.

•

The influence of microarchitecture on fatigue life exhibited similar trends for both irradiated and non-irradiated bone.

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

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

Osteoderm microstructure of Riostegotherium yanei, the oldest Xenarthra

Abstract: Riostegotherium yanei from the Itaboraí Basin, Brazil, is the oldest known Xenarthra. This paper aims to describe the internal morphology of the osteoderms of Riostegotherium yanei from the perspective of histology and micro-CT approaches, expanding the available data on cingulate osteoderm microstructure. Seven osteoderms of R. yanei were used for the internal microstructure description and eight of Dasypus novemcinctus for comparison. The osteoderms of Riostegotherium yanei lacks the diploë-like structure typical of glyptodonts but has a three-layered structure composed of two layers of non-Haversian compact bone enclosing a central layer of primary and secondary osteons. This internal organization is distinct from other Astegotheriini of comparable age, but similar to Dasypus. The 3D reconstruction of Riostegotherium yanei revealed two patterns of internal organization. Pattern 1 of movable osteoderm is composed of large remodeled areas at the base and a more compact bone at the tongue; in Pattern 2 (both movable and buckler), the internal cavities are much smaller, more numerous, and more interconnected to each other. In one buckler osteoderm, the cavities are organized somewhat radially with a compact central region (Pattern 1). Pattern 1 of both movable and buckler osteoderms resemble that of Dasypus.

Cortical bone vessel identification and quantification on contrast-enhanced MR images

Background

Cortical bone porosity is a major determinant of bone strength. Despite the biomechanical importance of cortical bone porosity, the biological drivers of cortical porosity are unknown. The content of cortical pore space can indicate pore expansion mechanisms; both of the primary components of pore space, vessels and adipocytes, have been implicated in pore expansion. Dynamic contrast-enhanced MRI (DCE-MRI) is widely used in vessel detection in cardiovascular studies, but has not been applied to visualize vessels within cortical bone. In this study, we have developed a multimodal DCE-MRI and high resolution peripheral QCT (HR-pQCT) acquisition and image processing pipeline to detect vessel-filled cortical bone pores.

Methods

For this in vivo human study, 19 volunteers (10 males and 9 females; mean age =63±5) were recruited. Both distal and ultra-distal regions of the non-dominant tibia were imaged by HR-pQCT (82 µm nominal resolution) for bone structure segmentation and by 3T DCE-MRI (Gadavist; 9 min scan time; temporal resolution =30 sec; voxel size 230×230×500 µm³) for vessel visualization. The DCE-MRI was registered to the HR-pQCT volume and the voxels within the MRI cortical bone region were extracted. Features of the DCE data were calculated and voxels were categorized by a 2-stage hierarchical kmeans clustering algorithm to determine which voxels represent vessels. Vessel volume fraction (volume ratio of vessels to cortical bone), vessel density (average vessel count per cortical bone volume), and average vessel volume (mean volume of vessels) were calculated to quantify the status of vessel-filled pores in cortical bone. To examine spatial resolution and perform validation, a virtual phantom with 5 channel sizes and an applied pseudo enhancement curve was processed through the proposed image processing pipeline. Overlap volume ratio and Dice coefficient was calculated to measure the similarity between the detected vessel map and ground truth.

Results

In the human study, mean vessel volume fraction was 2.2%±1.0%, mean vessel density was 0.68±0.27 vessel/mm³, and mean average vessel volume was 0.032±0.012 mm³/vessel. Signal intensity for detected vessel voxels increased during the scan, while signal for non-vessel voxels within pores did not enhance. In the validation phantom, channels with diameter 250 µm or greater were detected successfully, with volume ratio equal to 1 and Dice coefficient above 0.6. Both statistics decreased dramatically for channel sizes less than 250 µm.

Conclusions

We have a developed a multi-modal image acquisition and processing pipeline that successfully detects vessels within cortical bone pores. The performance of this technique degrades for vessel diameters below the in-plane spatial resolution of the DCE-MRI acquisition. This approach can be applied to investigate the biological systems associated with cortical pore expansion.

Towards a Connectomic Description of the Osteocyte Lacunocanalicular Network in Bone

PURPOSE OF REVIEW

Osteocytes are the most abundant bone cells. They are completely encased in mineralized tissue, sitting inside lacunae that are connected by a multitude of canaliculi. In recent years, the osteocyte network has been shown to fulfill endocrine functions and to communicate with a number of other organs. This review addresses emerging knowledge on the connectome of the lacunocanalicular network in different types of bone tissue.

RECENT FINDINGS

Recent advances in three-dimensional imaging technology started to reveal parameters that are well known from general theory to characterize the function of networks, such as network density, degree of nodes, or shortest path length through the network. The connectome of the lacunocanalicular network differs in some aspects between lamellar and woven bone and seems to change with age. More research is needed to relate network structure to function, such as intercellular transport or communication and its role in mechanosensation, as well as to understand the effect of diseases.

Intracortical Bone Mechanics Are Related to Pore Morphology and Remodeling in Human Bone

During aging and in osteoporosis, cortical bone becomes more porous, making it more fragile and susceptible to fractures. The aim of this study was to investigate the intracortical compression- induced strain energy distribution, and determine whether intracortical pores associated with high strain energy density (SED) in the surrounding bone matrix have a different morphology and distribution, as well as different remodeling characteristics than matrix with normal SED. Fibular diaphyseal specimens from 20 patients undergoing a jaw reconstruction (age range 41 to 75 years; 14 men and 6 women) were studied. Bone specimens were µCT-scanned, plastic embedded, and sectioned for histology. Three-dimensional microfinite element models of each specimen were tested in compression, and the SED of the bone immediately surrounding the intracortical pores was calculated within a plane of interest corresponding to the histological sections. The SED of a pore, relative to the distribution of the SED of all pores in each specimen, was used to classify pores as either a high or normal SED pore. Pores with high SED were larger, less circular, and were located closer to the endosteal surface of the cortex than normal SED pores (p < 0.001). Histological analysis of the remodeling events generating the pores revealed that the high SED pores compared with normal SED pores had 13.3-fold higher odds of being an erosive (70%) or formative (7%) pore versus a quiescent pore (p < 0.001), 5.9-fold higher odds of resulting from remodeling upon existing pores (type 2 pore) versus remodeling generating new pores (type 1 pore) (p < 0.001), and 3.2-fold higher odds of being a coalescing type 2 pore versus a noncoalescing type 2 pore (p < 0.001). Overall, the study demonstrates a strong relationship between cortical bone mechanics and pore morphology, distribution, and remodeling characteristics in human fibular bone. © 2018 American Society for Bone and Mineral Research.

Identifying animal taxa used to manufacture bone tools during the Middle Stone Age at Sibudu, South Africa: Results of a CT-rendered histological analysis

This paper presents the histological characterisation of a selection of worked bone artefacts from Middle Stone Age layers at Sibudu cave, South Africa. Histographic rendering is achieved using high-resolution Computed Tomography, which is non-destructive and facilitates three-dimensional histologic analysis. Excellent congruency in image quality was achieved with previous studies using this method. The results show that most of the artefact fragments contain mostly primary lamellar tissue, which is the bone tissue best adapted to withstand impact stresses. This indicates that bone with greater elastic properties was chosen. Histological characterisation allows the identification of animal taxa. Based on the sample analysed in this paper, Perissodactyla bone was used predominantly in the older layers at the site. Artiodactyla are represented throughout but appear far more frequently in the later (post-Howiesons Poort onwards) layers. Some of the Artiodactyla specimens have high proportions of Haversian tissue, reducing elasticity. The higher percentages of Haversian tissue in the post-Howiesons Poort artefacts relative to Holocene examples from southern Africa suggests that people may have started experimenting with bone from different animal taxa at this time and had not yet learned to eliminate the mechanically weaker secondary tissue. Apart from mechanical considerations, possible cultural constraints governing raw material selection is also explored.

A three-dimensional geometric quantification of human cortical canals using an innovative method with micro-computed tomographic data

The complex architecture of bone has been investigated for several decades. Some pioneer works proved an existing link between microstructure and external mechanical loading applied on bone. Due to sinuous network of canals and limitations of experimental acquisition technique, there has been little quantitative analysis of three-dimensional description of cortical network. The aim of this study is to provide an algorithmic process, using Python 3.5, in order to identify 3D geometrical characteristics of voids considered as canals. This script is based on micro-computed tomographic slices of two bone samples harvested from the humerus and femur of male cadaveric subject. Slice images are obtained from 2.94 μm isotropic resolution. This study provides a generic method of image processing which considers beam hardening artefact so as to avoid heuristic choice of global threshold value. The novelty of this work is the quantification of numerous three-dimensional canals features, such as orientation or canal length, but also connectivity features, such as opening angle, and the accurate definition of canals as voids which ranges from connectivity to possibly another intersection. The script was applied to one humeral and one femoral samples in order to analyse the difference in architecture between bearing and non-bearing cortical bones. This preliminary study reveals that the femoral specimen is more porous than the humeral one whereas the canal network is denser and more connected.

The effect of growth rate on the three-dimensional orientation of vascular canals in the cortical bone of broiler chickens

Vascular canals in cortical bone during growth and development typically show an anisotropic pattern with canals falling into three main categories: circumferential, radial, and longitudinal. Two major hypotheses attempt to explain the preferred orientations in bone: that vascular canal orientation is optimized to resist a predominant strain direction from functional loading, or that it reflects growth requirements and velocity. We use a controlled growth experiment in broiler chickens to investigate the effect of growth rate on vascular canal orientation. Using feed restriction we set up a fast growing control group and a slow growing restricted group. We compared the microstructure in the humerus and the femur at 42 days of age using synchrotron micro‐computed tomography (micro‐CT), a three‐dimensional (3D) method that visualizes the full canal network. We measured the 3D orientation of each canal in the whole cross‐section of the bone cortex using a set of custom imagej scripts. Using these orientations we compute laminar, radial, and longitudinal indices that measure the proportion of circumferential, radial, and longitudinal canals, by unit of length, in the cortex. Following previous studies we hypothesized that vascular canal orientation is related to growth, with radial canals linked to a faster growth rate and related to functional loading through a high laminar index in flight bones which reflects torsional loading resulting from active flight. The control group had final body weights that were nearly twice the final weights of the restricted group and higher absolute growth rates. We found consistent patterns in the comparison between the humerus and the femur in both groups, with the humerus having higher laminar and longitudinal indices, and a lower radial index than the femur. The control group had higher radial indices and lower laminar and longitudinal indices in both the humerus and the femur than the restricted group. The higher radial indices in our control group point to a link between radial canals and faster growth, and between laminar canals and slower growth, while the higher laminar indices in the humerus point to a link between circumferential canals and torsional loading. Overall, our results indicate that the orientation of the cortical canal network in a bone is the consequence of a complex interaction between the growth rate of that bone and functional loading environment.

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

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

Simultaneous visualisation of calcified bone microstructure and intracortical vasculature using synchrotron X-ray phase contrast-enhanced tomography

3D imaging of the bone vasculature is of key importance in the understanding of skeletal disease. As blood vessels in bone are deeply encased in the calcified matrix, imaging techniques that are applicable to soft tissues are generally difficult or impossible to apply to the skeleton. While canals in cortical bone can readily be identified and characterised in X-ray computed tomographic data in 3D, the soft tissue comprising blood vessels that are putatively contained within the canal structures does not provide sufficient image contrast necessary for image segmentation. Here, we report an approach that allows for rapid, simultaneous visualisation of calcified bone tissue and the vasculature within the calcified bone matrix. Using synchrotron X-ray phase contrast-enhanced tomography we show exemplar data with intracortical capillaries uncovered at sub-micrometre level without the need for any staining or contrast agent. Using the tibiofibular junction of 15 week-old C57BL/6 mice post mortem, we show the bone cortical porosity simultaneously along with the soft tissue comprising the vasculature. Validation with histology confirms that we can resolve individual capillaries. This imaging approach could be easily applied to other skeletal sites and transgenic models, and could improve our understanding of the role the vasculature plays in bone disease.

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.

Accuracy and reproducibility of mouse cortical bone microporosity as quantified by desktop microcomputed tomography

Bone's microporosity plays important roles in bone biology and bone mechanical quality. In this study, we explored the accuracy and reproducibility of nondestructive desktop μCT for 3D visualization and subsequent morphometric analysis of mouse cortical bone microporosity including the vascular canal network and osteocyte lacunae. The accuracy of measurements was evaluated in five murine fibula using confocal laser scanning microscopy (CLSM) in conjunction with Fluorescein isothiocyanate (FITC) staining as the reference method. The reproducibility of μCT-derived cortical bone microstructural indices was examined in 10 fibulae of C57Bl/6J male mice at a nominal resolution of 700 nanometer. Three repeated measurements were made on different days. An excellent correlation between μCT and CLSM was observed for both mean lacuna volume (r = 0.98, p = 0.002) and for mean lacuna orientation (r = 0.93, p = 0.02). Whereas the two techniques showed no significant differences for these parameters, the mean lacuna sphericity acquired from μCT was significantly higher than CLSM (p = 0.01). Reproducibility was high, with precision errors (PE) of 1.57-4.69% for lacuna parameters, and of 1.01-9.45% for vascular canal parameters. Intraclass correlation coefficient (ICC) showed a high reliability of the measurements, ranging from 0.998-1.000 for cortical parameters, 0.973-0.999 for vascular canal parameters and 0.755-0.991 for lacuna parameters. In conclusion, desktop μCT is a valuable tool to quantify the 3D characteristics of bone vascular canals as well as lacunae which can be applied to intact murine bones with high accuracy and reproducibility. Thus, μCT might be an important tool to improve our understanding of the physiological and biomechanical significance of these cannular and lacunar structure in cortical bone.

3D characterization of morphological changes in the intervertebral disc and endplate during aging: A propagation phase contrast synchrotron micro-tomography study

A better understanding of functional changes in the intervertebral disc (IVD) and interaction with endplate is essential to elucidate the pathogenesis of IVD degeneration disease (IDDD). To date, the simultaneous depiction of 3D micro-architectural changes of endplate with aging and interaction with IVD remains a technical challenge. We aim to characterize the 3D morphology changes of endplate and IVD during aging using PPCST. The lumbar vertebral level 4/5 IVDs harvested from 15-day-, 4- and 24-month-old mice were initially evaluated by PPCST with histological sections subsequently analyzed to confirm the imaging efficiency. Quantitative assessments of age-related trends after aging, including mean diameter, volume fraction and connectivity of the canals, and endplate porosity and thickness, reached a peak at 4 months and significantly decreased at 24 months. The IVD volume consistently exhibited same trend of variation with the endplate after aging. In this study, PPCST simultaneously provided comprehensive details of 3D morphological changes of the IVD and canal network in the endplate and the interaction after aging. The results suggest that PPCST has the potential to provide a new platform for attaining a deeper insight into the pathogenesis of IDDD, providing potential therapeutic targets.

Evaluating differential nuclear DNA yield rates and osteocyte numbers among human bone tissue types: A synchrotron radiation micro-CT approach

Molecular human identification has conventionally focused on DNA sampling from dense, weight-bearing cortical bone tissue, typically from femora or tibiae. A comparison of skeletal elements from three contemporary individuals demonstrated that elements with high quantities of cancellous bone yielded nuclear DNA at the highest rates, suggesting that preferentially sampling cortical bone may be suboptimal (Mundorff & Davoren, 2014). Despite these findings, the reason for the differential DNA yields between cortical and cancellous bone tissues remains unknown. The primary goal of this work is to ascertain whether differences in bone microstructure can be used to explain differential nuclear DNA yield among bone tissue types observed by Mundorff and Davoren (2014), with a focus on osteocytes and the three-dimensional (3D) quantification of their associated lacunae. Osteocytes and other bone cells are recognized to house DNA in bone tissue, thus examining the density of their lacunae may explain why nuclear DNA yield rates differ among bone tissue types. Lacunae were visualized and quantified using synchrotron radiation-based micro-Computed Tomographic imaging (SR micro-CT). Volumes of interest (VOIs) from cortical and cancellous bone tissues (n=129) were comparatively analyzed from the three skeletons sampled for Mundorff and Davoren's (2014) study. Analyses tested the primary hypothesis that the abundance and density of osteocytes (inferred from their lacunar spaces) vary between cortical and cancellous bone tissue types. Results demonstrated that osteocyte lacunar abundance and density vary between cortical and cancellous bone tissue types, with cortical bone VOIs containing a higher lacunar abundance and density. We found that the osteocyte lacunar density values are independent of nuclear DNA yield, suggesting an alternative explanation for the higher nuclear DNA yields from bones with greater quantities of cancellous bone tissue. The use of SR micro-CT allowed for a scale of analysis that revealed a high range of variation in lacunar abundance in both tissue types. Moreover, high-resolution SR micro-CT imaging revealed potential soft tissue remnants within marrow spaces not visible macroscopically. It is hypothesized that soft tissue remnants observed among the trabeculae of skeletal elements with high quantities of cancellous bone tissue are responsible for the high nuclear DNA yields. These findings have significant implications for bone-sample selection for nuclear DNA analysis in a forensic context when skeletal remains are recovered from the ground surface.

Asymmetric bone remodeling in mandibular and maxillary tori

OBJECTIVES

Tori are frequent paucisymptomatic bony outgrowths of the oral cavity in three locations: torus palatinus (TP), mandibularis (TM), and maxillaris (TMax). Their usually described histological characteristics are unspecific: normal cortical bone with more or less trabecular bone. The aim of this study was to describe tori's specific morphological and histomorphometric characteristics.

MATERIALS AND METHODS

Histological characteristics in a series of 18 tori collected after surgical removal were analyzed. Microcomputed tomography provided a 3D analysis. Mineral apposition rate (MAR) was measured after double tetracycline labeling. Osteoid tissue was identified by Goldner's trichrome and osteoclasts by the tartrate resistant acid phosphatase identification in undecalcified sections. Iron and aluminum were detected by histochemical staining methods. Osteoid thickness and MAR were determined at the outer surface of the torus and in the Haversian canals.

RESULTS

Tori appeared made of lamellar Haversian bone in 16/18 cases. Osteoid thickness did not differ between the outer surface and within the canals. An asymmetric bone remodeling was observed in the Haversian canals of 15 tori: osteoid seams were deposited on the side close to the free torus surface, and Howship's lacunae with numerous osteoclasts were observed on the opposite side. A high MAR was found at the surface of the tori and within the canals. There was no iron or aluminum deposit.

CONCLUSIONS

Tori may be characterized by a specific asymmetric bone remodeling which seems to determine their shape.

CLINICAL RELEVANCE

This finding could constitute a specific histological feature allowing to differentiate tori from exostoses. Graphical abstract.

Multiscale and multimodality computed tomography for cortical bone analysis

In clinical studies, high resolution peripheral quantitative computed tomography (HR-pQCT) is used to separately evaluate cortical bone and trabecular bone with an isotropic voxel of 82 µm³, and typical cortical parameters are cortical density (D.comp), thickness (Ct.Th), and porosity (Ct.Po). In vitro, micro-computed tomography (micro-CT) is used to explore the internal cortical bone micro-structure with isotropic voxels and high resolution synchrotron radiation (SR); micro-CT is considered the 'gold standard'. In 16 tibias and 8 femurs, HR-pQCT measurements were compared to conventional micro-CT measurements. To test modality effects, conventional micro-CT measurements were compared to SR micro-CT measurements at 7.5 µm³; SR micro-CT measurements were also tested at different voxel sizes for the femurs, specifically, 7.5 µm³ versus 2.8 µm³. D.comp (r  =  -0.88, p  <  10^-3) was the parameter best correlated with porosity (Po.V/TV). The correlation was not affected by the removal of pores under 130 µm. Ct.Th was also significantly highly correlated (r  =  -0.89 p  <  10^-3), while Ct.Po was correlated with its counterpart Po.V/TV (r  =  0.74, p  <  10^-3). From SR micro-CT and conventional micro-CT at 7.5 µm³ in matching areas, Po.V/TV and pore diameter were underestimated in conventional micro-CT with mean  ±  standard deviation (SD) biases of  -2.5  ±  1.9% and  -0.08  ±  0.08 mm, respectively. In contrast, pore number (Po.N) and pore separation (Po.Sp) were overestimated with mean  ±  SD biases of  +0.03  ±  0.04 mm^-1 and  +0.02  ±  0.04 mm, respectively. The results from the tibia and femur were similar when the results of SR micro-CT at 7.5 µm³ and 2.8 µm³ were compared. Po.V/TV, specific surface of pores (Po.S/Po.V), and Po.N were underestimated with mean biases of  -1.7  ±  0.9%, -4.6  ±  4.4 mm^-1, and  -0.26  ±  0.15 mm^-1, respectively. In contrast, pore spacing was overestimated at 7.5 µm³ compared to 2.8 µm³ with mean biases of 0.05  ±  0.03 mm. Cortical bone measurements from HR-pQCT images provided consistent results compared to those obtained using conventional micro-CT at the distal tibia. D.comp was highly correlated to Po.V/TV because it considers both the micro-porosity (Haversian systems) and macro-porosity (resorption lacunae) of cortical bone. The complexity of canal organization, (including shape, connectivity, and surface) are not fully considered in conventional micro-CT in relation to beam hardening and cone beam reconstruction artifacts. With the exception of Po.V/TV measurements, morphological and topological measurements depend on the characteristics of the x-ray beam, and to a lesser extent, on image resolution.

Cortical Bone Porosity: What Is It, Why Is It Important, and How Can We Detect It?

There is growing recognition of the role of micro-architecture in osteoporotic bone loss and fragility. This trend has been driven by advances in imaging technology, which have enabled a transition from measures of mass to micro-architecture. Imaging trabecular bone has been a key research focus, but advances in resolution have also enabled the detection of cortical bone micro-architecture, particularly the network of vascular canals, commonly referred to as 'cortical porosity.' This review aims to provide an overview of what this level of porosity is, why it is important, and how it can be characterized by imaging. Moving beyond a 'trabeculocentric' view of bone loss holds the potential to improve diagnosis and monitoring of interventions. Furthermore, cortical porosity is intimately linked to the remodeling process, which underpins bone loss, and thus a larger potential exists to improve our fundamental understanding of bone health through imaging of both humans and animal models.

Three-dimensional reconstruction of Haversian systems in human cortical bone using synchrotron radiation-based micro-CT: morphology and quantification of branching and transverse connections across age

This study uses synchrotron radiation-based micro-computed tomography (CT) scans to reconstruct three-dimensional networks of Haversian systems in human cortical bone in order to observe and analyse interconnectivity of Haversian systems and the development of total Haversian networks across different ages. A better knowledge of how Haversian systems interact with each other is essential to improve understanding of remodeling mechanisms and bone maintenance; however, previous methodological approaches (e.g. serial sections) did not reveal enough detail to follow the specific morphology of Haversian branching, for example. Accordingly, the aim of the present study was to identify the morphological diversity of branching patterns and transverse connections, and to understand how they change with age. Two types of branching morphologies were identified: lateral branching, resulting in small osteon branches bifurcating off of larger Haversian canals; and dichotomous branching, the formation of two new osteonal branches from one. The reconstructions in this study also suggest that Haversian systems frequently target previously existing systems as a path for their course, resulting in a cross-sectional morphology frequently referred to as 'type II osteons'. Transverse connections were diverse in their course from linear to oblique to curvy. Quantitative assessment of age-related trends indicates that while in younger human individuals transverse connections were most common, in older individuals more evidence of connections resulting from Haversian systems growing inside previously existing systems was found. Despite these changes in morphological characteristics, a relatively constant degree of overall interconnectivity is maintained throughout life. Altogether, the present study reveals important details about Haversian systems and their relation to each other that can be used towards a better understanding of cortical bone remodeling as well as a more accurate interpretation of morphological variants of osteons in cross-sectional microscopy. Permitting visibility of reversal lines, synchrotron radiation-based micro-CT is a valuable tool for the reconstruction of Haversian systems, and future analyses have the potential to further improve understanding of various important aspects of bone growth, maintenance and health.

Number and distribution of nutrient foramina within the femoral neck and their relationship to the retinacula of Weitbrecht: an anatomical study

Precise knowledge of the vascular supply of the femoral head is critical when contemplating surgery around the femoral head and neck junction. To determine the blood supply to the femoral neck, 2417 nutrient foramina from 76 cadaveric specimens were analyzed based on size, number, and distribution. Within the subcapital, transcervical, and basicervical regions of the femoral neck, the largest numbers of nutrient foramina were found on the superior (lateral) surface, followed by the anterior and posterior surfaces, and then the inferior (medial) surface (all p < 0.001). The diameters of most of the nutrient foramina were less than 1 mm. For the posterior and superior surfaces, the nutrient foramina in the basicervical region were significantly larger than those within the transcervical or subcapital regions (nutrient foramina >2 mm posteriorly: 23.6, 12.7, and 9.0 % in the basicervical, transcervical, and subcapital regions, respectively; superiorly: 23.7 vs. 15.4 vs. 16.8 %, respectively). In conclusion, neither the anterior nor the inferior surfaces in the basicervical, transcervical, and subcapital regions showed any significant differences in nutrient foraminal size. The areas containing densely distributed nutrient foramina were consistent with the regions covered by the retinacula of Weitbrecht.

Inter-individual changes in cortical bone three-dimensional microstructure and elastic coefficient have opposite effects on radial sound speed

Knowledge about simultaneous contributions of tissue microstructure and elastic properties on ultrasound speed in cortical bone is limited. In a previous study, porosities and elastic coefficients of cortical bone in human femurs were shown to change with age. In the present study, influences of inter-individual and site-dependent variation in cortical bone microstructure and elastic properties on radial speed of sound (SOS; at 4, 6, and 8 MHz) were investigated using three-dimensional (3D) finite difference time domain modeling. Models with fixed (nominal model) and sample-specific (sample-specific model) values of radial elastic coefficients were compared. Elastic coefficients and microstructure for samples (n = 24) of human femoral shafts (n = 6) were derived using scanning acoustic microscopy and micro-computed tomography images, respectively. Porosity-related SOS varied more extensively in nominal models than in sample-specific models. Linear correlation between pore separation and SOS was similar (R = 0.8, p < 0.01, for 4 MHz) for both models. The determination coefficient (R(2)= 0.75, p < 0.05) between porosity and radial SOS, especially at 4 MHz, was highest in the posterior quadrant. The determination coefficient was lower for models with sample-specific values of radial elastic coefficient implemented (R(2) < 0.33, p < 0.05), than for nominal models (0.48 < R(2)< 0.63, p < 0.05). This information could be useful in in vivo pulse-echo cortical thickness measurements applying constant SOS.

Does 3D orientation account for variation in osteon morphology assessed by 2D histology?

The primary microstructural unit of cortical bone, the secondary osteon or Haversian system, is widely assumed to have a cylindrical shape. It is generally accepted that osteons are roughly circular in cross-section and deviations from circularity have been attributed to deviations from longitudinal orientation. To our knowledge this idealized geometric relationship, which assumes osteons are perfect cylinders, has not been rigorously explored. As such, we sought to explore two research questions: (i) Does the orientation of osteons in 3D explain variation in shapes visualized in 2D? (ii) Can differences in osteon 3D orientation explain previously reported age-related differences observed in their 2D cross-sectional shape (e.g. more circular shape and decreased area with age)? To address these questions we utilized a combination of 2D histology to identify osteon shape and superimposed micro-computed tomography data to assess osteon orientation in 3D based upon the osteonal canal. Shape was assessed by the inverse of Aspect Ratio (On.AspR(-1), based on a fitted ellipse) - which ranged from 0 (infinitely elongated shape) to 1 (perfectly circular). A sample (n = 27) of human female anterior femoral cortical bone samples from across the human lifespan (20-87 years) were included in the analysis, which involved 1418 osteons. The overall mean measure of On.AspR(-1) was 0.703 (1.42 Aspect Ratio). Mean osteon orientation was 79.1° (90° being longitudinal). While we anticipated a positive relation between orientation and On.AspR(-1), we found the opposite - a weak negative correlation (with more oblique 3D osteon alignment, the 2D shape became more circular as reflected by increased On.AspR(-1)). When analysis of covariance was performed with age and orientation as covariates, the negative relation with orientation was replaced by a significant relation with age alone. This relation with age accounted for 41% of the variation of On.AspR(-1). The results revealed that osteons, on average, are not circular in cross-section and that 3D orientation cannot account for deviation from circular shape. Osteons thus are strictly speaking not cylinders, as they tend to have elliptical cross-sections. We observed that osteons did become less elliptical in cross-section with age independent of orientation - suggesting this is a real change in morphology.

Poromicromechanics reveals that physiological bone strains induce osteocyte-stimulating lacunar pressure

Mechanical loads which are macroscopically acting onto bony organs, are known to influence the activities of biological cells located in the pore spaces of bone, in particular so the signaling and production processes mediated by osteocytes. The exact mechanisms by which osteocytes are actually able to “feel” the mechanical loading and changes thereof, has been the subject of numerous studies, and, while several hypotheses have been brought forth over time, this topic has remained a matter of debate. Relaxation times reported in a recent experimental study of Gardinier et al. (Bone 46(4):1075–1081, 2010) strongly suggest that the lacunar pores are likely to experience, during typical physiological load cycles, not only fluid transport, but also undrained conditions. The latter entail the buildup of lacunar pore pressures, which we here quantify by means of a thorough multiscale modeling approach. In particular, the proposed model is based on classical poroelasticity theory, and able to account for multiple pore spaces. First, the model reveals distinct nonlinear dependencies of the resulting lacunar (and vascular) pore pressures on the underlying bone composition, highlighting the importance of a rigorous multiscale approach for appropriate computation of the aforementioned pore pressures. Then, the derived equations are evaluated for macroscopic (uniaxial as well as hydrostatic) mechanical loading of physiological magnitude. The resulting model-predicted pore pressures agree very well with the pressures that have been revealed, by means of in vitro studies, to be of adequate magnitude for modulating the responses of biological cells, including osteocytes. This underlines that osteocytes may respond to many types of loading stimuli at the same time, in particular so to fluid flow and hydrostatic pressure.

Spatial distribution of intracortical porosity varies across age and sex

Cortical bone porosity is a major determinant of strength, stiffness, and fracture toughness of cortical tissue. The goal of this work was to investigate changes in spatial distribution and microstructure of cortical porosity associated with aging in men and women. The specific aims were to: 1) develop an automated technique for spatial analysis of cortical microstructure based on HR-pQCT data, and; 2) apply this technique to explore sex- and age-specific spatial distribution and microstructure of porosity within the cortex. We evaluated HR-pQCT images of the distal tibia from a cross-sectional cohort of 145 individuals, characterizing detectable pores as being in the endosteal, midcortical, or periosteal layers of the cortex. Metrics describing porosity, pore number, and pore size were quantified within each layer and compared across sexes, age groups, and cortical layers. The elderly cohort (65-78 years, n=22) displayed higher values than the young cohort (20-29 years, n=29) for all parameters both globally and within each layer. While all three layers displayed significant age-related porosity increases, the greatest difference in porosity between the young and elderly cohort was in the midcortical layer (+344%, p<0.001). Similarly, the midcortical layer reflected the greatest differences between young and elderly cohorts in both pore number (+243%, p<0.001) and size (+28%, p<0.001). Females displayed greater age-related changes in porosity and pore number than males. Females and males displayed comparable small to non-significant changes with age in pore size. In summary, considerable variability exists in the spatial distribution of detectable cortical porosity at the distal tibia, and this variability is dependent on age and sex. Intracortical pore distribution analysis may ultimately provide insight into both mechanisms of pore network expansion and biomechanical consequences of pore distribution.

Assessment of bone vascularization and its role in bone remodeling

Bone is a composite organ that fulfils several interconnected functions, which may conflict with each other in pathological conditions. Bone vascularization is at the interface between these functions. The roles of bone vascularization are better documented in bone development, growth and modeling than in bone remodeling. However, every bone remodeling unit is associated with a capillary in both cortical and trabecular envelopes. Here we summarize the most recent data on vessel involvement in bone remodeling, and we present the characteristics of bone vascularization. Finally, we describe the various techniques used for bone vessel imaging and quantitative assessment, including histology, immunohistochemistry, microtomography and intravital microscopy. Studying the role of vascularization in adult bone should provide benefits for the understanding and treatment of metabolic bone diseases.

A comparison of methods for in vivo assessment of cortical porosity in the human appendicular skeleton

The recent advent of high-resolution peripheral quantitative computed tomography (HR-pQCT) provides new opportunities to measure in vivo human bone microarchitecture. Increasingly, cortical porosity (CtPo) is of particular interest due to its relationship with bone quality and turnover. The two approaches that have emerged to measure CtPo from HR-pQCT are threshold-based and density-based methods, and the purpose of this work was to compare the performance of each against a gold-standard synchrotron radiation micro-computed tomography (SRμCT) measurement. Human cadaveric cortical bone specimens (N=23) were measured by SRμCT and HR-pQCT, and high correlations were found for both methods. The density-based approach had an r2=0.939 (95% confidence interval (CI) of +6.17% to +20.99%) and consistently overestimated porosity as measured by SRμCT, while the threshold-based approach had an r2=0.977 and consistently underestimated porosity (95% CI of -2.60% to -10.76%). The density-based approach is prone to beam hardening artifacts and susceptible to natural variations of tissue mineral density (TMD), but is less affected by motion artifacts that may occur in in vivo scans. The threshold-based method has the advantage that it provides structural information that complements the cortical porosity measure, such as number of pores and connectivity, and can accurately detect the larger pores which are the most relevant to bone biomechanical strength. With the first generation HR-pQCT systems the accuracy of detecting pores larger than 140 μm diameter is excellent (r2=0.983; 95% CI of -4.88% to +2.45%). The accuracy of the threshold-based method will improve as new HR-pQCT systems emerge and provide a robust quantitative approach to measure cortical porosity.

Modalities for Visualization of Cortical Bone Remodeling: The Past, Present, and Future

Bone's ability to respond to load-related phenomena and repair microdamage is achieved through the remodeling process, which renews bone by activating groups of cells known as basic multicellular units (BMUs). The products of BMUs, secondary osteons, have been extensively studied via classic two-dimensional techniques, which have provided a wealth of information on how histomorphology relates to skeletal structure and function. Remodeling is critical in maintaining healthy bone tissue; however, in osteoporotic bone, imbalanced resorption results in increased bone fragility and fracture. With increasing life expectancy, such degenerative bone diseases are a growing concern. The three-dimensional (3D) morphology of BMUs and their correlation to function, however, are not well-characterized and little is known about the specific mechanisms that initiate and regulate their activity within cortical bone. We believe a key limitation has been the lack of 3D information about BMU morphology and activity. Thus, this paper reviews methodologies for 3D investigation of cortical bone remodeling and, specifically, structures associated with BMU activity (resorption spaces) and the structures they create (secondary osteons), spanning from histology to modern ex vivo imaging modalities, culminating with the growing potential of in vivo imaging. This collection of papers focuses on the theme of "putting the 'why' back into bone architecture." Remodeling is one of two mechanisms "how" bone structure is dynamically modified and thus an improved 3D understanding of this fundamental process is crucial to ultimately understanding the "why."

Double incretin receptor knock-out (DIRKO) mice present with alterations of trabecular and cortical micromorphology and bone strength

A role for gut hormone in bone physiology has been suspected. We evidenced alterations of microstructural morphology (trabecular and cortical) and bone strength (both at the whole-bone--and tissue-level) in double incretin receptor knock-out (DIRKO) mice as compared to wild-type littermates. These results support a role for gut hormones in bone physiology.

INTRODUCTION

The two incretins, glucose-dependent insulinotropic polypeptide (GIP) and glucagon-like peptide-1 (GLP-1), have been shown to control bone remodeling and strength. However, lessons from single incretin receptor knock-out mice highlighted a compensatory mechanism induced by elevated sensitivity to the other gut hormone. As such, it is unclear whether the bone alterations observed in GIP or GLP-1 receptor deficient animals resulted from the lack of a functional gut hormone receptor, or by higher sensitivity for the other gut hormone. The aims of the present study were to investigate the bone microstructural morphology, as well as bone tissue properties, in double incretin receptor knock-out (DIRKO) mice.

METHODS

Twenty-six-week-old DIRKO mice were age- and sex-matched with wild-type (WT) littermates. Bone microstructural morphology was assessed at the femur by microCT and quantitative X-ray imaging, while tissue properties were investigated by quantitative backscattered electron imaging and Fourier-transformed infrared microscopy. Bone mechanical response was assessed at the whole-bone- and tissue-level by 3-point bending and nanoindentation, respectively.

RESULTS

As compared to WT animals, DIRKO mice presented significant augmentations in trabecular bone mass and trabecular number whereas bone outer diameter, cortical thickness, and cortical area were reduced. At the whole-bone-level, yield stress, ultimate stress, and post-yield work to fracture were significantly reduced in DIRKO animals. At the tissue-level, only collagen maturity was reduced by 9 % in DIRKO mice leading to reductions in maximum load, hardness, and dissipated energy.

CONCLUSIONS

This study demonstrated the critical role of gut hormones in controlling bone microstructural morphology and tissue properties.

Bone hierarchical structure in three dimensions

Bone is a complex hierarchically structured family of materials that includes a network of cells and their interconnected cell processes. New insights into the 3-D structure of various bone materials (mainly rat and human lamellar bone and minipig fibrolamellar bone) were obtained using a focused ion beam electron microscope and the serial surface view method. These studies revealed the presence of two different materials, the major material being the well-known ordered arrays of mineralized collagen fibrils and associated macromolecules, and the minor component being a relatively disordered material composed of individual collagen fibrils with no preferred orientation, with crystals inside and possibly between fibrils, and extensive ground mass. Significantly, the canaliculi and their cell processes are confined within the disordered material. Here we present a new hierarchical scheme for several bone tissue types that incorporates these two materials. The new scheme updates the hierarchical scheme presented by Weiner and Wagner (1998). We discuss the structures at different hierarchical levels with the aim of obtaining further insights into structure-function-related questions, as well as defining some remaining unanswered questions.

Structural analysis of cortical porosity applied to HR-pQCT data

PURPOSE

The investigation of cortical porosity is an important aspect of understanding biological, pathoetiological, and biomechanical processes occurring within the skeleton. With the emergence of HR-pQCT as a noninvasive tool suitable for clinical use, cortical porosity at appendicular sites can be directly visualized in vivo. The aim of this study was to introduce a novel topological analysis of the cortical pore network for HR-pQCT data and determine the influence of resolution on measures of cortical pore network microstructure and topology.

METHODS

Cadaveric radii were scanned using HR-pQCT at two different voxel sizes (41 and 82 μm) and also using μCT at a voxel size of 18 μm. HR-pQCT and μCT image sets were spatially coregistered. Segmentation and quantification of cortical porosity (Ct.Po) and mean pore diameter (Ct.Po.Dm) were achieved using an established extended cortical analysis technique. Topological classification of individual pores was performed using topology-preserving skeletonization and multicolor dilation algorithms. Based on the pore skeleton topological classification, the following parameters were quantified: total number of planar surface-skeleton canals (N.Slabs), tubular curve-skeleton canals (N.Tubes), and junction elements (N.Junctions), mean slab volume (Slab.Vol), mean tube volume (Tube.Vol), mean slab orientation (Slab.θ), mean tube orientation (Tube.θ), N.Slabs/N.Tubes, and integral (total) slab volume/integral tube volume (iSlab.Vol/iTube.Vol). An in vivo reproducibility study was also conducted to assess short-term precision of the topology parameters. Precision error was characterized using root mean square coefficient of variation (RMSCV%).

RESULTS

Correlations to μCT values for Ct.Po were significant for both the 41 and 82 μm HR-pQCT data (41: r(2) = 0.82, p < 0.001, 82: r(2) = 0.75, p < 0.001). For Ct.Po.Dm, only the 41 μm data were significantly predictive of μCT values (r(2) = 0.72, p < 0.01) Data at both HR-pQCT voxel sizes were strongly predictive of the μCT values for N.Slabs (41: r(2) = 0.93, p < 0.001; 82: r(2) = 0.84, p < 0.001), N.Tubes (41: r(2) = 0.94, p < 0.001; 82: r(2) = 0.84, p < 0.001), and N.Junctions (41: r(2) = 0.93, p < 0.001; 82: r(2) = 0.78, p < 0.001), though proportional bias was evident in these correlations. Weak correlations were seen for iSlab.Vol/iTube.Vol at both voxel sizes (41: r(2) = 0.52, p < 0.01; 82: r(2) = 0.39, p < 0.05). Slab.Vol was significantly correlated to μCT data at 41 μm (r(2) = 0.60, p < 0.01) but not at 82 μm, while Tube.Vol was significantly correlated at both voxel sizes (41: r(2) = 0.79, p < 0.001; 82: r(2) = 0.68, p < 0.01). In vivo precision error for these parameters ranged from 2.31 to 9.68 RMSCV%.

CONCLUSIONS

Strong correlations between μCT- and HR-pQCT-derived measurements were found, particularly in HR-pQCT images obtained at 41 μm. These data are in agreement with our previous study investigating the effect of voxel size on standard HR-pQCT metrics of trabecular and cortical microstructure, and extend our previous findings to include topological descriptors of the cortical pore network.