A quantitative framework for the 3D characterization of the osteocyte lacunar system

Deep learning models to map osteocyte networks can successfully distinguish between young and aged bone

Osteocytes, the most abundant and mechanosensitive cells in bone tissue, play a pivotal role in bone homeostasis and mechano-responsiveness, orchestrating the intricate balance between bone formation and resorption under daily activity. Studying osteocyte connectivity and understanding their intricate arrangement within the lacunar canalicular network (LCN) is essential for unraveling bone physiology. This is particularly true as our bones age, which is associated with decreased integrity of the osteocyte network, disrupted mass transport, and lower sensitivity to the mechanical stimuli that allow the skeleton to adapt to changing demands. Much work has been carried out to investigate this relationship, often involving high resolution microscopy of discrete fragments of this network, alongside advanced computational modelling of individual cells. However, traditional methods of segmenting and measuring osteocyte connectomics are time-consuming and labour-intensive, often hindered by human subjectivity and limited throughput. In this study, we explore the application of deep learning and computer vision techniques to automate the segmentation and measurement of osteocyte connectomics, enabling more efficient and accurate analysis. We compare several state-of-the-art computer vision models (U-Nets and Vision Transformers) to successfully segment the LCN, finding that an Attention U-Net model can accurately segment and measure 81.8% of osteocytes and 42.1% of dendritic processes, when compared to manual labelling. While further development is required, we demonstrate that this degree of accuracy is already sufficient to distinguish between bones of young (2 month old) and aged (36 month old) mice, as well as capturing the degeneration induced by genetic modification of osteocytes. By harnessing the power of these advanced technologies, further developments can unravel the complexities of osteocyte networks in unprecedented detail, revolutionising our understanding of bone health and disease.

An in silico micro-multiphysics agent-based approach for simulating bone regeneration in a mouse femur defect model

Bone defects represent a challenging clinical problem as they can lead to non-union. In silico models are well suited to study bone regeneration under varying conditions by linking both cellular and systems scales. This paper presents an in silico micro-multiphysics agent-based (micro-MPA) model for bone regeneration following an osteotomy. The model includes vasculature, bone, and immune cells, as well as their interaction with the local environment. The model was calibrated by time-lapsed micro-computed tomography data of femoral osteotomies in C57Bl/6J mice, and the differences between predicted bone volume fractions and the longitudinal in vivo measurements were quantitatively evaluated using root mean square error (RMSE). The model performed well in simulating bone regeneration across the osteotomy gap, with no difference (5.5% RMSE, p = 0.68) between the in silico and in vivo groups for the 5-week healing period - from the inflammatory phase to the remodelling phase - in the volume spanning the osteotomy gap. Overall, the proposed micro-MPA model was able to simulate the influence of the local mechanical environment on bone regeneration, and both this environment and cytokine concentrations were found to be key factors in promoting bone regeneration. Further, the validated model matched clinical observations that larger gap sizes correlate with worse healing outcomes and ultimately simulated non-union. This model could help design and guide future experimental studies in bone repair, by identifying which are the most critical in vivo experiments to perform.

Utility and Limitations of TALLYHO/JngJ as a Model for Type 2 Diabetes-Induced Bone Disease

Type 2 diabetes mellitus (T2DM) increases risk of fractures due to bone microstructural and material deficits, though the mechanisms remain unclear. Preclinical models mimicking diabetic bone disease are required to further understand its pathogenesis. The TALLYHO/JngJ (TH) mouse is a polygenic model recapitulating adolescent-onset T2DM in humans. Due to incomplete penetrance of the phenotype ~25% of male TH mice never develop hyperglycemia, providing a strain-matched nondiabetic control. We performed a comprehensive characterization of the metabolic and skeletal phenotype of diabetic TH mice and compared them to either their nondiabetic TH controls or the recommended SWR/J controls to evaluate their suitability to study diabetic bone disease in humans. Compared to both controls, male TH mice with T2DM exhibited higher blood glucose levels, weight along with impaired glucose tolerance and insulin sensitivity. TH mice with/without T2DM displayed higher cortical bone parameters and lower trabecular bone parameters in the femurs and vertebrae compared to SWR/J. The mechanical properties remained unchanged for all three groups except for a low-energy failure in TH mice with T2DM only compared to SWR/J. Histomorphometry analyses only revealed higher number of osteoclasts and osteocytes for SWR/J compared to both groups of TH. Bone turnover markers procollagen type 1 N-terminal propeptide (P1NP) and tartrate-resistant acid phosphatase (TRAP) were low for both groups of TH mice compared to SWR/J. Silver nitrate staining of the femurs revealed low number of osteocyte lacunar and dendrites in TH mice with T2DM. Three-dimensional assessment showed reduced lacunar parameters in trabecular and cortical bone. Notably, osteocyte morphology changed in TH mice with T2DM compared to SWR/J. In summary, our study highlights the utility of the TH mouse to study T2DM, but not necessarily T2DM-induced bone disease, as there were no differences in bone strength and bone cell parameters between diabetic and non-diabetic TH mice. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Microvascular Disease Associates with Larger Osteocyte Lacunae in Cortical Bone in Type 2 Diabetes Mellitus

Clinical studies indicate that microvascular disease (MVD) affects bone microstructure and decreases bone strength in type 2 diabetes mellitus (T2D). Osteocytes are housed in small voids within the bone matrix and lacunae and act as sensors of mechanical forces in bone. These cells regulate osteoclastic bone resorption and osteoblastic bone formation as well as osteocytic perilacunar remodeling. We hypothesized that MVD changes morphometric osteocyte lacunar parameters in individuals with T2D. We collected iliac crest bone biopsies from 35 individuals (10 female, 25 male) with T2D with MVD (15%) or without MVD (21%) with a median age of 67 years (interquartile range [IQR] 62-72 years). The participants were included based on c-peptide levels >700 pmol L-1, absence of anti-GAD65 antibodies, and glycated hemoglobin (HbA1c) levels between 40 and 82 mmol mol-1 or 5.8% and 9.7%, respectively. We assessed osteocyte lacunar morphometric parameters in trabecular and cortical bone regions using micro-computed tomography (micro-CT) at a nominal resolution of 1.2 μm voxel size. The cortical osteocyte lacunar volume (Lc.V) was 7.7% larger (p = 0.05) and more spherical (Lc.Sr, p < 0.01) in the T2D + MVD group. Using linear regression, we found that lacunar density (Lc.N/BV) in trabecular but not cortical bone was associated with HbA1c (p < 0.05, R 2 = 0.067) independently of MVD. Furthermore, Lc.V was larger and Lc.Sr higher in the center than in the periphery of the trabecular and cortical bone regions (p < 0.05). In conclusion, these data imply that MVD may impair skeletal integrity, possibly contributing to increased skeletal fragility in T2D complicated by MVD. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

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

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

Osteocyte lacunae in transiliac bone biopsy samples across life span

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

Abnormal morphological features of osteocyte lacunae in adolescent idiopathic scoliosis: A large-scale assessment by ultra-high-resolution micro-computed tomography

AIM

Abnormal osteocyte lacunar morphology in adolescent idiopathic scoliosis (AIS) has been reported while the results were limited by the number of osteocyte lacunae being quantified. The present study aimed to validate previous findings through (a) comparing morphological features of osteocyte lacunae between AIS patients and controls in spine and ilium using a large-scale assessment, and (b) investigating whether there is an association between the acquired morphological features of osteocyte lacunae and disease severity in AIS.

METHOD

Trabecular bone tissue of the facet joint of human vertebrae on both concave and convex sides at the apex of the scoliotic curve were collected from 4 AIS and 5 congenital scoliosis (CS) patients, and also at the same anatomic site from 3 non-scoliosis (NS) subjects intraoperatively. Trabecular bone tissue from ilium was obtained from 12 AIS vs 9 NS subjects during surgery. Osteocyte lacunae were assessed using ultra-high-resolution micro-computed tomography. Clinical information such as age, body mass index (BMI) and radiological Cobb angle of the major curve were collected.

RESULTS

There was no significant difference between density of osteocyte lacuna and bone volume fraction (BV/TV) between groups. A total of 230,076 and 78,758 osteocyte lacunae from facet joints of apical vertebra of scoliotic curve and iliac bone were included in the analysis, respectively. In facet joint bone biopsies, lacunar stretch (Lc.St) was higher, and lacunar equancy (Lc.Eq), lacunar oblateness (Lc.Ob), and lacunar sphericity (Lc.Sr) were lower in AIS and CS groups when compared with NS group. CA side was associated with higher Lc.St when compared with CX side. In iliac bone biopsies, Lc.Ob was higher and lacunar surface area (Lc.S) was lower in AIS group than NS group. Median values of Lc.St, Lc.Eq and Lc.Sr were significantly associated with radiological Cobb angle with adjustment for age and BMI (R-squared: 0.576, 0.558 and 0.543, respectively).

CONCLUSIONS

This large-scale assessment of osteocyte lacunae confirms that AIS osteocyte lacunae are more oblate in iliac bone that is less influenced by asymmetric loading of the deformed spine than the vertebrae. Shape of osteocyte lacunae in iliac bone is associated with radiological Cobb angle of the major curve in AIS patients, suggesting the likelihood of systemic abnormal osteocyte morphology in AIS. Osteocyte lacunae from concave side of scoliotic curves were more stretched in both AIS and CS groups, which is likely secondary to asymmetric mechanical loading.

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

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

Large-scale quantification of human osteocyte lacunar morphological biomarkers as assessed by ultra-high-resolution desktop micro-computed tomography

Ultra-high-resolution imaging of the osteocyte lacuno-canalicular network (LCN) three-dimensionally (3D) in a high-throughput fashion has greatly improved the morphological knowledge about the constituent structures - positioning them as potential biomarkers. Technologies such as serial focused ion beam/scanning electron microscopy (FIB/SEM) and confocal scanning laser microscopy (CLSM) can image in extremely high resolution, yet only capture a small number of lacunae. Synchrotron radiation computed tomography (SR-CT) can image with both high resolution and high throughput but has a limited availability. Desktop micro-computed tomography (micro-CT) provides an attractive balance: high-throughput imaging on the micron level without the restrictions of SR-CT availability. In this study, accuracy, reproducibility, and sensitivity of large-scale quantification of human osteocyte lacunar morphometries were assessed by ultra-high-resolution desktop micro-computed tomography. For this purpose, thirty-one transiliac human bone biopsies containing trabecular and cortical regions were imaged using ultra-high-resolution desktop micro-CT at a nominal isotropic voxel resolution of 1.2 µm. The resulting 3D images were segmented, component labeled, and the following morphometric parameters of 7.71 million lacunae were measured: Lacunar number (Lc.N), density (Lc.N/BV), porosity (Lc.TV/BV), volume (Lc.V), surface area (Lc.S), surface area to volume ratio (Lc.S/Lc.V), stretch (Lc.St), oblateness (Lc.Ob), sphericity (Lc.Sr), equancy (Lc.Eq), and angle (Lc.θ). Accuracy was quantified by comparing automated lacunar identification to manual identification. Mean true positive rate (TPR), false positive rate (FPR), and false negative rate (FNR) were 89.0%, 3.4%, and 11.0%, respectively. Regarding the reproducibility of lacunar morphometry from repeated measurements, precision errors were low (0.2-3.0%) and intraclass correlation coefficients were high (0.960-0.999). Significant differences between cortical and trabecular regions (p<0.001) existed for Lc.N/BV, Lc.TV/BV, local lacunar surface area (<Lc.S>), and local lacunar volume (<Lc.V>), all of which demonstrate the sensitivity of the method and are possible biomarker candidates. This study provides the foundation required for future large-scale morphometric studies using ultra-high-resolution desktop micro-CT and high-throughput analysis of millions of osteocyte lacunae in human bone samples.

Down to the Bone: A Novel Bio-Inspired Design Concept

The solutions provided through natural evolution of living creatures serve as an ingenious source of inspiration for many technological and applicative fields. Along these lines, bone-inspired concepts lead to fascinating advances in product design, architecture and garments, thanks to the bone’s exceptional combination of strength, toughness and lightness. Structural applications are inspired by the bone’s ability to resist fracture under a large spectrum of forces, while the high surface area and pore connectivity of bone architecture present exciting opportunities from an aesthetic point of view. Behind these inspirations, a disruptive common belief emerges: “down to the bone”, a journey in search of equality, universality and substantiality. Herein, we explore the current state of the art in bone-inspired applications in these fields, considering the two major categories of structural and aesthetic inspirations and discussing further technological developments.

Microstructural heterogeneity of the collagenous network in the loaded and unloaded periodontal ligament and its biomechanical implications

The periodontal ligament (PDL) is a highly heterogeneous fibrous connective tissue and plays a critical role in distributing occlusal forces and regulating tissue remodeling. Its mechanical properties are largely determined by the extracellular matrix, comprising a collagenous fiber network interacting with the capillary system as well as interstitial fluid containing proteoglycans. While the phase-contrast micro-CT technique has portrayed the 3D microscopic heterogeneity of PDL, the topological parameters of its network, which is crucial to understanding the multiscale constitutive behavior of this tissue, has not been characterized quantitatively. This study aimed to provide new understanding of such microscopic heterogeneity of the PDL with quantifications at both tissue and collagen network levels in a spatial manner, by combining phase-contrast micro-CT imaging and a purpose-built image processing algorithm for fiber analysis. Both variations within a PDL and among the PDL with different shapes, i.e. round-shaped and kidney-shaped PDLs, are described in terms of tissue thickness, fiber distribution, local fiber densities, and fiber orientation (namely azimuthal and elevation angles). Furthermore, the tissue and collagen fiber network responses to mechanical loading were evaluated in a similar manner. A 3D helical alignment pattern was observed in the fiber network, which appears to regulate and adapt a screw-like tooth motion under occlusion. The microstructural heterogeneity quantified here allows development of sample-specific constitutive models to characterize the PDL's functional and pathological loading responses, thereby providing a new multiscale framework for advancing our knowledge of this complex limited mobility soft-hard tissue interface.

Ten-Year Simulation of the Effects of Denosumab on Bone Remodeling in Human Biopsies

Postmenopausal osteoporosis is a disease manifesting in degradation of bone mass and microarchitecture, leading to weakening and increased risk of fracture. Clinical trials are an essential tool for evaluating new treatments and may provide further mechanistic understanding of their effects in vivo. However, the histomorphometry from clinical trials is limited to 2D images and reflects single time points. Biochemical markers of bone turnover give global insight into a drug's action, but not the local dynamics of the bone remodeling process and the cells involved. Additionally, comparative trials necessitate separate treatment groups, meaning only aggregated measures can be compared. In this study, in silico modeling based on histomorphometry and pharmacokinetic data was used to assess the effects of treatment versus control on μCT scans of the same biopsy samples over time, matching the changes in bone volume fraction observed in biopsies from denosumab and placebo groups through year 10 of the FREEDOM Extension trial. In the simulation, treatment decreased osteoclast number, which led to a modest increase in trabecular thickness and osteocyte stress shielding. Long-term bone turnover suppression led to increased RANKL production, followed by a small increase in osteoclast number at the end of the 6-month-dosing interval, especially at the end of the Extension study. Lack of treatment led to a significant loss of bone mass and structure. The study's results show how in silico models can generate predictions of denosumab cellular action over a 10-year period, matching static and dynamic morphometric measures assessed in clinical biopsies. The use of in silico models with clinical trial data can be a method to gain further insight into fundamental bone biology and how treatments can perturb this. With rigorous validation, such models could be used for informing the design of clinical trials, such that the number of participants could be reduced to a minimum to show efficacy. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

A Review on Multiscale Bone Damage: From the Clinical to the Research Perspective

The investigation of bone damage processes is a crucial point to understand the mechanisms of age-related bone fractures. In order to reduce their impact, early diagnosis is key. The intricate architecture of bone and the complexity of multiscale damage processes make fracture prediction an ambitious goal. This review, supported by a detailed analysis of bone damage physical principles, aims at presenting a critical overview of how multiscale imaging techniques could be used to implement reliable and validated numerical tools for the study and prediction of bone fractures. While macro- and meso-scale imaging find applications in clinical practice, micro- and nano-scale imaging are commonly used only for research purposes, with the objective to extract fragility indexes. Those images are used as a source for multiscale computational damage models. As an example, micro-computed tomography (micro-CT) images in combination with micro-finite element models could shed some light on the comprehension of the interaction between micro-cracks and micro-scale bone features. As future insights, the actual state of technology suggests that these models could be a potential substitute for invasive clinical practice for the prediction of age-related bone fractures. However, the translation to clinical practice requires experimental validation, which is still in progress.

Multithreaded two-pass connected components labelling and particle analysis in ImageJ

Sequential region labelling, also known as connected components labelling, is a standard image segmentation problem that joins contiguous foreground pixels into blobs. Despite its long development history and widespread use across diverse domains such as bone biology, materials science and geology, connected components labelling can still form a bottleneck in image processing pipelines. Here, I describe a multithreaded implementation of classical two-pass sequential region labelling and introduce an efficient collision resolution step, 'bucket fountain'. Code was validated on test images and against commercial software (Avizo). It was performance tested on images from 2 MB (161 particles) to 6.5 GB (437 508 particles) to determine whether theoretical linear scaling (O(n)) had been achieved, and on 1-40 CPU threads to measure speed improvements due to multithreading. The new implementation achieves linear scaling (b = 0.905-1.052, time ∝ pixelsb ; R 2 = 0.985-0.996), which improves with increasing thread number up to 8-16 threads, suggesting that it is memory bandwidth limited. This new implementation of sequential region labelling reduces the time required from hours to a few tens of seconds for images of several GB, and is limited only by hardware scale. It is available open source and free of charge in BoneJ.

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

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

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.

Quantification of the bone lacunocanalicular network from 3D X-ray phase nanotomography images

There is a growing interest in developing 3D microscopy for the exploration of thick biological tissues. Recently, 3D X-ray nanocomputerised tomography has proven to be a suitable technique for imaging the bone lacunocanalicular network. This interconnected structure is hosting the osteocytes which play a major role in maintaining bone quality through remodelling processes. 3D images have the potential to reveal the architecture of cellular networks, but their quantitative analysis remains a challenge due to the density and complexity of nanometre sized structures and the need to handle and process large datasets, for example, 2048³ voxels corresponding to 32 GB per individual image in our case. In this work, we propose an efficient image processing approach for the segmentation of the network and the extraction of characteristic parameters describing the 3D structure. These parameters include the density of lacunae, the porosity of lacunae and canaliculi, and morphological features of lacunae (volume, surface area, lengths, anisotropy etc.). We also introduce additional parameters describing the local environment of each lacuna and its canaliculi. The method is applied to analyse eight human femoral cortical bone samples imaged by magnified X-ray phase nanotomography with a voxel size of 120 nm, which was found to be a good compromise to resolve canaliculi while keeping a sufficiently large field of view of 246 μm in 3D. The analysis was performed on a total of 2077 lacunae showing an average length, width and depth of 17.1 μm × 9.2 μm × 4.4 μm, with an average number of 58.2 canaliculi per lacuna and a total lacuno-canalicular porosity of 1.12%. The reported descriptive parameters provide information on the 3D organisation of the lacuno-canalicular network in human bones.

Quantifying intracortical bone microstructure: A critical appraisal of 2D and 3D approaches for assessing vascular canals and osteocyte lacunae

Describing and quantifying vascular canal orientation and volume of osteocyte lacunae in bone is important in studies of bone growth, mechanics, health and disease. It is also an important element in analysing fossil bone in palaeohistology, key to understanding the growth, life and death of extinct animals. Often, bone microstructure is studied using two-dimensional (2D) sections, and three-dimensional (3D) shape and orientation of structures are estimated by modelling the structures using idealised geometries based on information from their cross sections. However, these methods rely on structures meeting strict geometric assumptions. Recently, 3D methods have been proposed which could provide a more accurate and robust approach to bone histology, but these have not been tested in direct comparison with their 2D counterparts in terms of accuracy and sensitivity to deviations from model assumptions. We compared 2D and 3D methodologies for estimating key microstructural traits using a combination of experimental and idealised test data sets. We generated populations of cylinders (canals) and ellipsoids (osteocyte lacunae), varying the cross-sectional aspect ratios of cylinders and orientation of ellipsoids to test sensitivity to deviations from cylindricality and longitudinal orientation, respectively. Using published methods, based on 2D sections and 3D data sets, we estimated cylinder orientation and ellipsoid volume. We applied the same methods to six CT data sets of duck cortical bone, using the full volumes for 3D measurements and single CT slices to represent 2D sections. Using in silico test data sets that did deviate from ideal cylinders and ellipsoids resulted in inaccurate estimates of cylinder or canal orientation, and reduced accuracy in estimates of ellipsoid and lacunar volume. These results highlight the importance of using appropriate 3D imaging and quantitative methods for quantifying volume and orientation of 3D structures and offer approaches to significantly enhance our understanding of bone physiology based on accurate measures for bone microstructures.

Genetic Ablation of Osteopontin in Osteomalacic Hyp Mice Partially Rescues the Deficient Mineralization Without Correcting Hypophosphatemia

PHEX is predominantly expressed by bone and tooth-forming cells, and its inactivating mutations in X-linked hypophosphatemia (XLH) lead to renal phosphate wasting and severe hypomineralization of bones and teeth. Also present in XLH are hallmark hypomineralized periosteocytic lesions (POLs, halos) that persist despite stable correction of serum phosphate (P_i ) that improves bulk bone mineralization. In XLH, mineralization-inhibiting osteopontin (OPN, a substrate for PHEX) accumulates in the extracellular matrix of bone. To investigate how OPN functions in Hyp mice (a model for XLH), double-null (Hyp;Opn^-/- ) mice were generated. Undecalcified histomorphometry performed on lumbar vertebrae revealed that Hyp;Opn^-/- mice had significantly reduced osteoid area/bone area (OV/BV) and osteoid thickness of trabecular bone as compared to Hyp mice, despite being as hypophosphatemic as Hyp littermate controls. However, tibias examined by synchrotron radiation micro-CT showed that mineral lacunar volumes remained abnormally enlarged in these double-null mice. When Hyp;Opn^-/- mice were fed a high-P_i diet, serum P_i concentration increased, and OV/BV and osteoid thickness normalized, yet mineral lacunar area remained abnormally enlarged. Enpp1 and Ankh gene expression were increased in double-null mice fed a high-P_i diet, potentially indicating a role for elevated inhibitory pyrophosphate (PP_i ) in the absence of OPN. To further investigate the persistence of POLs in Hyp mice despite stable correction of serum P_i , immunohistochemistry for OPN on Hyp mice fed a high-P_i diet showed elevated OPN in the osteocyte pericellular lacunar matrix as compared to Hyp mice fed a control diet. This suggests that POLs persisting in Hyp mice despite correction of serum P_i may be attributable to the well-known upregulation of mineralization-inhibiting OPN by P_i , and its accumulation in the osteocyte pericellular matrix. This study shows that OPN contributes to osteomalacia in Hyp mice, and that genetic ablation of OPN in Hyp mice improves the mineralization phenotype independent of systemic P_i -regulating factors. © 2020 American Society for Bone and Mineral Research.

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

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

Osteocyte lacunar strain determination using multiscale finite element analysis

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

•

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

•

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

•

Lacunar strains decreased as the perilacunar modulus increased.

•

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

•

Increases in lacuna size also led to increased lacunar strains

Decoding rejuvenating effects of mechanical loading on skeletal aging using in vivo μCT imaging and deep learning

Throughout the process of aging, dynamic changes of bone material, micro- and macro-architecture result in a loss of strength and therefore in an increased likelihood of fragility fractures. To date, precise contributions of age-related changes in bone (re)modeling and (de)mineralization dynamics to this fragility increase are not completely understood. Here, we present an image-based deep learning approach to quantitatively describe the effects of short-term aging and adaptive response to cyclic loading applied to proximal mouse tibiae and fibulae. Our approach allowed us to perform an end-to-end age prediction based on μCT imaging to determine the dynamic biological process of aging during a two week period, therefore permitting short-term bone aging analysis with 95% accuracy in predicting time points. In a second application, our deep learning analysis reveals that two weeks of in vivo mechanical loading are associated with an underlying rejuvenating effect of 5 days. Additionally, by quantitatively analyzing the learning process, we could, for the first time, identify the localization of the age-relevant encoded information and demonstrate 89% load-induced similarity of these locations in the loaded tibia with younger control bones. These data therefore suggest that our method enables identifying a general prognostic phenotype of a certain skeletal age as well as a temporal and localized loading-treatment effect on this apparent skeletal age for the studied mouse tibia and fibula. Future translational applications of this method may provide an improved decision-support method for osteoporosis treatment at relatively low cost. STATEMENT OF SIGNIFICANCE: Bone is a highly complex and dynamic structure that undergoes changes during the course of aging as well as in response to external stimuli, such as loading. Automatic assessment of "age" and "state" of the bone may lead to early prognosis of deceases such as osteoporosis and enables evaluating the effects of certain treatments. Here, we present an artificial intelligence-based method capable of automatically predicting the skeletal age from μCT images with 95% accuracy. Additionally, we utilize it to demonstrate the rejuvenation effects of in-vivo loading treatment on bones. We further, for the first time, break down aging-related local changes in bone by quantitatively analyzing "what the age assessment model has learned" and use this information to investigate the structural details of rejuvenation process.

Three-dimensional characterization of osteocyte volumes at multiple scales, and its relationship with bone biology and genome evolution in ray-finned fishes

Osteocytes, cells embedded within the bone mineral matrix, inform on key aspects of vertebrate biology. In particular, a relationship between volumes of the osteocytes and bone growth and/or genome size has been proposed for several tetrapod lineages. However, the variation in osteocyte volume across different scales is poorly characterized and mostly relies on incomplete, two-dimensional information. In this study, we characterize the variation of osteocyte volumes in ray-finned fishes (Actinopterygii), a clade including more than half of modern vertebrate species in which osteocyte biology is poorly known. We use X-ray synchrotron micro-computed tomography (SRµCT) to achieve a three-dimensional visualization of osteocyte lacunae and direct measurement of their size (volumes). Our specimen sample is designed to characterize variation in osteocyte lacuna morphology at three scales: within a bone, among the bones of one individual and among species. At the intra-bone scale, we find that osteocyte lacunae vary noticeably in size between zones of organized and woven bone (being up to six times larger in woven bone), and across cyclical bone deposition. This is probably explained by differences in bone deposition rate, with larger osteocyte lacunae contained in bone that deposits faster. Osteocyte lacuna volumes vary 3.5-fold among the bones of an individual, and this cannot readily be explained by variation in bone growth rate or other currently observable factors. Finally, we find that genome size provides the best explanation of variation in osteocyte lacuna volume among species: actinopterygian taxa with larger genomes (polyploid taxa in particular) have larger osteocyte lacunae (with a ninefold variation in median osteocyte volume being measured). Our findings corroborate previous two-dimensional studies in tetrapods that also observed similar patterns of intra-individual variation and found a correlation with genome size. This opens new perspectives for further studies on bone evolution, physiology and palaeogenomics in actinopterygians, and vertebrates as a whole.

Assessment of the human bone lacuno-canalicular network at the nanoscale and impact of spatial resolution

Recently, increasing attention has been given to the study of osteocytes, the cells that are thought to play an important role in bone remodeling and in the mechanisms of bone fragility. The interconnected osteocyte system is deeply embedded inside the mineralized bone matrix and lies within a closely fitted porosity known as the lacuno-canalicular network. However, quantitative data on human samples remain scarce, mostly measured in 2D, and there are gaps to be filled in terms of spatial resolution. In this work, we present data on femoral samples from female donors imaged with isotropic 3D spatial resolution by magnified X-ray phase nano computerized-tomography. We report quantitative results on the 3D structure of canaliculi in human femoral bone imaged with a voxel size of 30 nm. We found that the lacuno-canalicular porosity occupies on average 1.45% of the total tissue volume, the ratio of the canalicular versus lacunar porosity is about 37.7%, and the primary number of canaliculi stemming from each lacuna is 79 on average. The examination of this number at different distances from the surface of the lacunae demonstrates branching in the canaliculi network. We analyzed the impact of spatial resolution on quantification by comparing parameters extracted from the same samples imaged with 120 nm and 30 nm voxel sizes. To avoid any bias related to the analysis region, the volumes at 120 nm and 30 nm were registered and cropped to the same field of view. Our results show that the measurements at 120 and 30 nm are strongly correlated in our data set but that the highest spatial resolution provides more accurate information on the canaliculi network and its branching properties.

Spaceflight-Induced Bone Tissue Changes that Affect Bone Quality and Increase Fracture Risk

PURPOSE OF REVIEW

Bone mineral density and systemic factors are used to assess skeletal health in astronauts. Yet, even in a general population, these measures fail to accurately predict when any individual will fracture. This review considers how long-duration human spaceflight requires evaluation of additional bone structural and material quality measures that contribute to microgravity-induced skeletal fragility.

RECENT FINDINGS

In both humans and small animal models following spaceflight, bone mass is compromised via reduced bone formation and elevated resorption levels. Concurrently, bone structural quality (e.g., trabecular microarchitecture) is diminished and the quality of bone material is reduced via impaired tissue mineralization, maturation, and maintenance (e.g., mediated by osteocytes). Bone structural and material quality are both affected by microgravity and may, together, jeopardize astronaut operational readiness and lead to increased fracture risk upon return to gravitational loading. Future studies need to directly evaluate how bone quality combines with diminished bone mass to influence bone strength and toughness (e.g., resistance to fracture). Bone quality assessment promises to identify novel biomarkers and therapeutic targets.

YAP and TAZ Mediate Osteocyte Perilacunar/Canalicular Remodeling

Bone fragility fractures are caused by low bone mass or impaired bone quality. Osteoblast/osteoclast coordination determines bone mass, but the factors that control bone quality are poorly understood. Osteocytes regulate osteoblast and osteoclast activity on bone surfaces but can also directly reorganize the bone matrix to improve bone quality through perilacunar/canalicular remodeling; however, the molecular mechanisms remain unclear. We previously found that deleting the transcriptional regulators Yes-associated protein (YAP) and transcriptional co-activator with PDZ-motif (TAZ) from osteoblast-lineage cells caused lethality in mice due to skeletal fragility. Here, we tested the hypothesis that YAP and TAZ regulate osteocyte-mediated bone remodeling by conditional ablation of both YAP and TAZ from mouse osteocytes using 8 kb-DMP1-Cre. Osteocyte-conditional YAP/TAZ deletion reduced bone mass and dysregulated matrix collagen content and organization, which together decreased bone mechanical properties. Further, YAP/TAZ deletion impaired osteocyte perilacunar/canalicular remodeling by reducing canalicular network density, length, and branching, as well as perilacunar flourochrome-labeled mineral deposition. Consistent with recent studies identifying TGF-β as a key inducer of osteocyte expression of matrix-remodeling enzymes, YAP/TAZ deletion in vivo decreased osteocyte expression of matrix proteases MMP13, MMP14, and CTSK. In vitro, pharmacologic inhibition of YAP/TAZ transcriptional activity in osteocyte-like cells abrogated TGF-β-induced matrix protease gene expression. Together, these data show that YAP and TAZ control bone matrix accrual, organization, and mechanical properties by regulating osteocyte-mediated bone remodeling. Elucidating the signaling pathways that control perilacunar/canalicular remodeling may enable future therapeutic targeting of bone quality to reverse skeletal fragility. © 2019 American Society for Bone and Mineral Research.

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.

The contribution of the pericanalicular matrix to mineral content in human osteonal bone

The osteocyte lacunar-canalicular network (LCN) penetrates bone and houses the osteocytes and their processes. Despite its rather low volume fraction, the LCN represents an outstanding large surface that is possibly used by the osteocytes to interact with the surrounding mineralized bone matrix thereby contributing to mineral homeostasis. The aim of this study was to quantitatively describe such contributions by spatially correlating the local density of the LCN with the mineral content at the same location in micrometer-sized volume elements in human osteons. For this purpose, 65 osteons from the femur midshaft from healthy adults (n = 4) and children (n = 2) were structurally characterized with two different techniques. The 3D structure of the LCN in the osteons was imaged with confocal laser scanning microscopy after staining the bone samples with rhodamine. Subsequent image analysis provided the canalicular length density, i.e. the total length of the canaliculi per unit volume (μm/μm³). Quantitative information on the mineral content (wt%Ca) from the identical regions was obtained using quantitative backscattered electron imaging. As the LCN-porosity lowers the mineral content, a negative correlation between Ca content and network density was expected. Calculations predict a reduction of around -0.97 fmol Ca per μm of network. However, the experiment revealed for 62 out of 65 osteons a positive correlation resulting in an average additional Ca loading of +1.15 fmol per μm of canalicular network, i.e. an accumulation of mineral has occurred at dense network regions. We hypothesize that this accumulation happens in the close vicinity of canaliculi forming mineral reservoirs that can be utilized by osteocytes. Significant differences found between individuals indicate that the extent of mineral loading of the reservoir zone reflects an important parameter for mineral homeostasis.

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

PURPOSE OF REVIEW

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

RECENT FINDINGS

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

Age-related changes in female mouse cortical bone microporosity

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

A new open-source tool for measuring 3D osteocyte lacunar geometries from confocal laser scanning microscopy reveals age-related changes to lacunar size and shape in cortical mouse bone

Osteocytes can participate in systemic mineral homeostasis through perilacunar maintenance and remodeling, where changes to osteocyte lacunar morphology may affect bone structural integrity, tissue strains, and osteocyte mechanosensitivity. Though aging is associated with both decreased bone quality and altered mineral metabolism, it is not known if osteocyte lacunae undergo age-related changes in geometry. In order to survey lacunar changes with age, we developed an open-source program whereby 3D osteocyte lacunae are automatically segmented and then subsequently reconstructed from confocal laser scanning microscopy (CLSM) depth stacks for quantitative analysis of geometry and orientation. This approach takes advantage of the availability and speed of CLSM while avoiding time-consuming and bias-prone manual segmentation. Unlike conventional approaches used to quantify osteocyte lacunar morphology, CLSM enables facile analysis in three-dimensions with clear identification of osteocyte lacunae. We report that 3D osteocyte lacunae measured by CLSM become smaller, more spherical, more oblate, more spatially disorganized, and more sparsely populated with increased age in C57Bl/6 mouse cortical bone in groups spanning 6-24 months old. Critically, these age-related changes are in large part not observed in 2D analyses from the same samples. These results (1) demonstrate proof-of-concept of an efficient method to quantitatively assess osteocyte lacunae in 3D for application to a wide range of studies and (2) motivate further inquiry into how changes to osteocyte lacunar geometries and perilacunar material contribute to diminished bone quality in aging.

Computational 3D imaging to quantify structural components and assembly of protein networks

Traditionally, protein structures have been described by the secondary structure architecture and fold arrangement. However, the relatively novel method of 3D confocal microscopy of fluorescent-protein-tagged networks in living cells allows resolving the detailed spatial organization of these networks. This provides new possibilities to predict network functionality, as structure and function seem to be linked at various scales. Here, we propose a quantitative approach using 3D confocal microscopy image data to describe protein networks based on their nano-structural characteristics. This analysis is constructed in four steps: (i) Segmentation of the microscopic raw data into a volume model and extraction of a spatial graph representing the protein network. (ii) Quantifying protein network gross morphology using the volume model. (iii) Quantifying protein network components using the spatial graph. (iv) Linking these two scales to obtain insights into network assembly. Here, we quantitatively describe the filamentous temperature sensitive Z protein network of the moss Physcomitrella patens and elucidate relations between network size and assembly details. Future applications will link network structure and functionality by tracking dynamic structural changes over time and comparing different states or types of networks, possibly allowing more precise identification of (mal) functions or the design of protein-engineered biomaterials for applications in regenerative medicine.

STATEMENT OF SIGNIFICANCE

Protein networks are highly complex and dynamic structures that play various roles in biological environments. Analyzing the detailed spatial structure of these networks may lead to new insight into biological functions and malfunctions. Here, we propose a tool set that extracts structural information at two scales of the protein network and allows therefore to address questions such as "how is the network built?" or "how networks grow?".

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

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

•

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

•

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

•

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

•

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

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

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

Sost Deficiency does not Alter Bone's Lacunar or Vascular Porosity in Mice

SCLEROSTIN (Sost) is expressed predominantly in osteocytes acting as a negative regulator of bone formation. In humans, mutations in the SOST gene lead to skeletal overgrowth and increased bone mineral density, suggesting that SCLEROSTIN is a key regulator of bone mass. The function of SCLEROSTIN as an inhibitor of bone formation is further supported by Sost knockout (KO) mice which display a high bone mass with elevated bone formation. Previous studies have indicated that Sost exerts its effect on bone formation through Wnt-mediated regulation of osteoblast differentiation, proliferation, and activity. Recent in vitro studies have also suggested that SCLEROSTIN regulates angiogenesis and osteoblast-to-osteocyte transition. Despite this wealth of knowledge of the mechanisms responsible for SCLEROSTIN action, no previous studies have examined whether SCLEROSTIN regulates osteocyte and vascular configuration in cortices of mouse tibia. Herein, we image tibiae from Sost KO mice and their wild-type (WT) counterparts with high-resolution CT to examine whether lack of SCLEROSTIN influences the morphometric properties of lacunae and vascular canal porosity relating to osteocytes and vessels within cortical bone. Male Sost KO and WT mice (n = 6/group) were sacrificed at 12 weeks of age. Fixed tibiae were analyzed using microCT to examine cortical bone mass and architecture. Then, samples were imaged by using benchtop and synchrotron nano-computed tomography at the tibiofibular junction. Our data, consistent with previous studies show that, Sost deficiency leads to significant enhancement of bone mass by cortical thickening and bigger cross-sectional area and we find that this occurs without modifications of tibial ellipticity, a measure of bone shape. In addition, our data show that there are no significant differences in any lacunar or vascular morphometric or geometric parameters between Sost KO mouse tibia and WT counterparts. We, therefore, conclude that the significant increases in bone mass induced by Sost deficiency are not accompanied by any significant modification in the density, organization, or shape of osteocyte lacunae or vascular content within the cortical bone. These data may imply that SCLEROSTIN does not modify the frequency of osteocytogenic recruitment of osteoblasts to initiate terminal osteocytic differentiation in mice.

Coalignment of osteocyte canaliculi and collagen fibers in human osteonal bone

During bone formation osteocytes get connected with each other via a dense network of canaliculi within the mineralized bone matrix. Important functions attributed to the osteocyte network include the control of bone remodeling and a contribution to mineral homeostasis. To detect structural clues of the formation and functionality of the network, this study analyzes the structure and orientation of the osteocyte lacuno-canalicular network (OLCN), specifically in relation to the concentric bone lamellae within human osteons. The network structure within 49 osteons from four samples of cortical bone from the femoral midshaft of middle-aged healthy women was determined by a combination of rhodamine staining and confocal laser scanning microscopy followed by computational image analysis. A quantitative evaluation showed that 64±1% of the canalicular length has an angle smaller than 30° to the direction towards the osteon center, while the lateral network - defined by an orientation angle larger than 60° - comprises 16±1%. With the same spatial periodicity as the bone lamellae, both radial and lateral network show variations in the network density and order. However, only the preferred orientation of the lateral network twists when crossing a lamella. This twist agrees with the preferred orientation of the fibrous collagen matrix. The chirality of the twist was found to be individual-specific. The coalignment between network and matrix extends to the orientation of the elongated osteocyte lacunae. The intimate link between OLCN and collagen matrix implies an interplay between osteocyte processes and the arrangement of the surrounding collagen fibers during osteoid formation.

Effect of Macroanatomic Bone Type and Estrogen Loss on Osteocyte Lacunar Properties in Healthy Adult Women

This is the first study to examine clinical human bone specimens by three-dimensional imaging to characterize osteocyte lacunar properties as a function of macroanatomic bone type and estrogen loss. We applied laboratory-based instrumentation [3D X-ray microscope (3DXRM), MicroXCT-200; Carl Zeiss/Xradia, Inc.] that reaches the same resolution as synchrotron microscopy. We used serial transiliac bone biopsy specimens to examine the effect of macroanatomic bone type and estrogen status on osteocyte lacunar properties. These properties include lacunar size (volume, axes lengths of the ellipsoidal lacunar voids), distribution (density, average near-neighbor lacunar distance), and shape factors (sphericity ratio, average eigenvalues, degree of equancy, elongation, and flatness) in both cortical and trabecular bone tissue. The lacunar properties (volume, surface area, density, near-neighbor distance, etc.) and the shape factors (E1, L1, L2, degree of equancy, degree of elongation) were different between cortical and trabecular bone regardless of estrogen status. In cortical bone and trabecular nodes, the lacunar void volume and surface area were either smaller or tended to be smaller in postmenopausal as compared to premenopausal women. The void volume-to-bone volume ratio of cortical bone showed declining trends with estrogen loss. While there were differences between trabecular and cortical bone tissue, the lacunar void sphericity ratio for trabecular struts shows decreasing trends in postmenopausal women. These data suggest that using 3DXRM can provide new insight into osteocyte lacunar properties in transiliac bone biopsies from patients with various skeletal disease/conditions and pharmaceutical treatments.

Porotic paradox: distribution of cortical bone pore sizes at nano- and micro-levels in healthy vs. fragile human bone

Bone is a remarkable biological nanocomposite material showing peculiar hierarchical organization from smaller (nano, micro) to larger (macro) length scales. Increased material porosity is considered as the main feature of fragile bone at larger length-scales. However, there is a shortage of quantitative information on bone porosity at smaller length-scales, as well as on the distribution of pore sizes in healthy vs. fragile bone. Therefore, here we investigated how healthy and fragile bones differ in pore volume and pore size distribution patterns, considering a wide range of mostly neglected pore sizes from nano to micron-length scales (7.5 to 15000 nm). Cortical bone specimens from four young healthy women (age: 35 ± 6 years) and five women with bone fracture (age: 82 ± 5 years) were analyzed by mercury porosimetry. Our findings showed that, surprisingly, fragile bone demonstrated lower pore volume at the measured scales. Furtnermore, pore size distribution showed differential patterns between healthy and fragile bones, where healthy bone showed especially high proportion of pores between 200 and 15000 nm. Therefore, although fragile bones are known for increased porosity at macroscopic level and level of tens or hundreds of microns as firmly established in the literature, our study with a unique assessment range of nano-to micron-sized pores reveal that osteoporosis does not imply increased porosity at all length scales. Our thorough assessment of bone porosity reveals a specific distribution of porosities at smaller length-scales and contributes to proper understanding of bone structure which is important for designing new biomimetic bone substitute materials.

Spatial heterogeneity in the canalicular density of the osteocyte network in human osteons

Osteocytes interconnect with each other forming an intricate cell network within the mineralized bone matrix. One important function of the osteocyte network is the mechano-regulation of bone remodeling, where a possible mechanism includes the fluid flow through the porosity housing the cell network - the osteocyte lacuno-canalicular network (OLCN). In our study the OLCN in human osteons was three-dimensionally imaged with the aim to obtain a quantitative description of the canalicular density and spatial variations of this quantity within osteons. The topology of the OLCN was determined by first staining the bone samples with rhodamine, then imaging the OLCN with confocal laser scanning microscopy and finally using image analysis to obtain a skeletonized version of the network for further analysis. In total 49 osteons were studied from the femoral cortical bone of four different middle-aged healthy women. The mean canalicular density given as length of the canaliculi in a unit volume was 0.074 ± 0.015 μm/μm³ (corresponding to 74 km/cm³). No correlation was found between the canalicular density and neither the size of the osteon nor the volume fraction occupied by osteocyte lacunae. Within osteons the canalicular density varied substantially with larger regions without any network. On average the canalicular density decreases when moving from the Haversian canal outwards towards the cement line. We hypothesize that a decrease in accessible canaliculi with tissue age as a result of micropetrosis can reduce the local mechanosensitivity of the bone. Systematic future studies on age- and disease-related changes on the topology of the OLCN have to demonstrate the diagnostic potential of the presented characterization method.

•

accurate three-dimensional representation of the canalicular network topology using confocal microscopy and image analysis

•

quantitative network description with an average canalicular density of 74 km per cubic centimeter in healthy human osteons

•

substantial variability of the canalicular density within osteons with 10% of the bone having a density twice the average

•

in osteons large regions lacking an accessible network were found

•

radial decrease of the canalicular density from the Haversian canal towards the cement line

Immobilization and long-term recovery results in large changes in bone structure and strength but no corresponding alterations of osteocyte lacunar properties

The ability of osteocytes to demineralize the perilacunar matrix, osteocytic osteolysis, and thereby participate directly in bone metabolism, is an aspect of osteocyte biology that has received increasing attention during the last couple of years. The aim of the present work was to investigate whether osteocyte lacunar properties change during immobilization and subsequent recovery. A rat cortical bone model with negligible Haversian remodeling effects was used, with temporary immobilization of one hindlimb induced by botulinum toxin. Several complementary techniques covering multiple length scales enabled correlation of osteocyte lacunar properties to changes observed on the organ and tissue level of femoral bone. Bone structural parameters measured by μCT and mechanical properties were compared to sub-micrometer resolution SR μCT data mapping an unprecedented number (1.85 million) of osteocyte lacunae. Immobilization induced a significant reduction in aBMD, bone volume, tissue volume, and load to fracture, as well as the muscle mass of rectus femoris. During the subsequent recovery period, the bone structural and mechanical properties were only partly regained in spite of a long-term (28weeks) study period. No significant changes in osteocyte lacunar volume, density, oblateness, stretch, or orientation were detected upon immobilization or subsequent recovery. In conclusion, the bone architecture and not osteocyte lacunar properties or bone material characteristics dominate the immobilization response as well as the subsequent recovery.

Intraskeletal variation in human cortical osteocyte lacunar density: Implications for bone quality assessment

Osteocytes and their lacunocanalicular network have been identified as the regulator of bone quality and function by exerting extensive influence over metabolic processes, mechanical adaptation, and mineral homeostasis. Recent research has shown that osteocyte apoptosis leads to a decrease in bone quality and increase in bone fragility mediated through its effects on remodeling. The purpose of this study is to investigate variation in cortical bone osteocyte lacunar density with respect to major factors including sex, age, and intracortical porosity to establish both regional and systemic trends. Samples from the midshaft femur, midshaft rib and distal one-third diaphysis of the radius were recovered from 30 modern cadaveric individuals (15 males and 15 females) ranging from 49 to 100 years old. Thick ground undecalcified histological (80 μm) cross-sections were made and imaged under bright field microscopy. Osteocyte lacunar density (Ot.Lc.N/B.Ar) and intracortical porosity (%Po.Ar) were quantified. No significant sex differences in Ot.Lc.N/B.Ar or %Po.Ar were found in any element. Linear regressions demonstrated a significant decrease in osteocyte lacunar density (Ot.Lc.N/B.Ar) and increase in intracortical porosity (%Po.Ar) with age for the sex-pooled sample in the femur (R² = 0.208, 0.297 respectively) and radius (R² = 0.108, 0.545 respectively). Age was unable to significantly predict osteocyte lacunar density or intracortical porosity in the rib (R² = 0.058, 0.114 respectively). Comparisons of regression coefficients demonstrated a systemic trend in the decrease in osteocyte lacunar density (Ot.Lc.N/B.Ar) and increase in intracortical porosity (%Po.Ar) with age. In each element, intracortical porosity was significantly negatively correlated with lacunar density for which the radius demonstrated the strongest relationship (r = - 0.746). Using pore number (Po.N) as a proxy for available vascularity to support the osteocyte population, Po.N was able to predict 61.8% of variation in osteocyte lacunar number (Ot.Lc.N) in the rib. The femur and radius also demonstrated significant relationships between these variables (R² = 0.560 and 0.397 respectively). The results from this study indicate that although the femur, radius and rib may be experiencing systemically influenced declines in osteocyte lacunar density, there may be differential effects at each anatomical site potentially due to age related changes in mechanical loading. With decreasing osteocyte lacunar density in each element, intracortical porosity increased with likely direct impacts on gross bone strength. This study provides a foundation upon which to build interpretations of osteocyte lacunar density values and their effect on differential fracture risk for aging individuals.

Impact of Non-Invasively Induced Motor Deficits on Tibial Cortical Properties in Mutant Lurcher Mice

It has been shown that Lurcher mutant mice have significantly altered motor abilities, regarding their motor coordination and muscular strength because of olivorecebellar degeneration. We assessed the response of the cross-sectional geometry and lacuno-canalicular network properties of the tibial mid-diaphyseal cortical bone to motor differences between Lurcher and wild-type (WT) male mice from the B6CBA strain. The first data set used in the cross-sectional geometry analysis consists of 16 mice of 4 months of age and 32 mice of 9 months of age. The second data set used in the lacunar-canalicular network analysis consists of 10 mice of 4 months of age. We compared two cross-sectional geometry and four lacunar-canalicular properties by I-region using the maximum and minimum second moment of area and anatomical orientation as well as H-regions using histological differences within a cross section. We identified inconsistent differences in the studied cross-sectional geometry properties between Lurcher and WT mice. The biggest significant difference between Lurcher and WT mice is found in the number of canaliculi, whereas in the other studied properties are only limited. Lurcher mice exhibit an increased number of canaliculi (p < 0.01) in all studied regions compared with the WT controls. The number of canaliculi is also negatively correlated with the distance from the centroid in the Lurcher and positively correlated in the WT mice. When the Lurcher and WT sample is pooled, the number of canaliculi and lacunar volume is increased in the posterior Imax region, and in addition, midcortical H-region exhibit lower number of canaliculi, lacuna to lacuna distance and increased lacunar volume. Our results indicate, that the importance of precise sample selection within cross sections in future studies is highlighted because of the histological heterogeneity of lacunar-canalicular network properties within the I-region and H-region in the mouse cortical bone.

A Single Chance to Contact Multiple Targets: Distinct Osteocyte Morphotypes Shed Light on the Cellular Mechanism Ensuring the Robust Formation of Osteocytic Networks

The formation of the complex osteocytic network relies on the emission of long cellular processes involved in communication, mechanical strain sensing, and bone turnover control. Newly deposited osteocytic processes rapidly become trapped within the calcifying matrix, and, therefore, they must adopt their definitive conformation and contact their targets in a single morphogenetic event. However, the cellular mechanisms ensuring the robustness of this unique mode of morphogenesis remain unknown. To address this issue, we examined the developing calvaria of the amphibian Xenopus tropicalis by confocal, two-photon, and super-resolution imaging, and described flattened osteocytes lying within a woven bone structured in lamellae of randomly oriented collagen fibers. While most cells emit peripheral and perpendicular processes, we report two osteocytes morphotypes, located at different depth within the bone matrix and exhibiting distinct number and orientation of perpendicular cell processes. We show that this pattern is conserved with the chick Gallus gallus and suggest that the cellular microenvironment, and more particularly cell-cell contact, plays a fundamental role in the induction and stabilization of osteocytic processes. We propose that this intrinsic property might have been evolutionarily selected for its ability to robustly generate self-organizing osteocytic networks harbored by the wide variety of bone shapes and architectures found in extant and extinct vertebrates.

Osteocyte-directed bone demineralization along canaliculi

The mammalian skeleton stores calcium and phosphate ions in bone matrix. Osteocytes in osteocyte lacunae extend numerous dendrites into canaliculi less than a micron in diameter and which are distributed throughout bone matrix. Although osteoclasts are the primary bone-resorbing cells, osteocytes also reportedly dissolve hydroxyapatite at peri-lacunar bone matrix. However, robust three-dimensional evidence for peri-canalicular bone mineral dissolution has been lacking. Here we applied a previously reported Talbot-defocus multiscan tomography method for synchrotron X-ray microscopy and analyzed the degree of bone mineralization in mouse cortical bone around the lacuno-canalicular network, which is connected both to blood vessels and the peri- and endosteum. We detected cylindrical low mineral density regions spreading around canaliculi derived from a subset of osteocytes. Transmission electron microscopy revealed both intact and demineralized bone matrix around the canaliculus. Peri-canalicular low mineral density regions were also observed in osteopetrotic mice lacking osteoclasts, indicating that osteoclasts are dispensable for peri-canalicular demineralization. These data suggest demineralization can occur from within bone through the canalicular system, and that peri-canalicular demineralization occurs not uniformly but directed by individual osteocytes. Blockade of peri-canalicular demineralization may be a therapeutic strategy to increase bone mass and quality.

3D X-ray ultra-microscopy of bone tissue

We review the current X-ray techniques with 3D imaging capability at the nano-scale: transmission X-ray microscopy, ptychography and in-line phase nano-tomography. We further review the different ultra-structural features that have so far been resolved: the lacuno-canalicular network, collagen orientation, nano-scale mineralization and their use as basis for mechanical simulations. X-ray computed tomography at the micro-metric scale is increasingly considered as the reference technique in imaging of bone micro-structure. The trend has been to push towards increasingly higher resolution. Due to the difficulty of realizing optics in the hard X-ray regime, the magnification has mainly been due to the use of visible light optics and indirect detection of the X-rays, which limits the attainable resolution with respect to the wavelength of the visible light used in detection. Recent developments in X-ray optics and instrumentation have allowed to implement several types of methods that achieve imaging that is limited in resolution by the X-ray wavelength, thus enabling computed tomography at the nano-scale. We review here the X-ray techniques with 3D imaging capability at the nano-scale: transmission X-ray microscopy, ptychography and in-line phase nano-tomography. Further, we review the different ultra-structural features that have so far been resolved and the applications that have been reported: imaging of the lacuno-canalicular network, direct analysis of collagen orientation, analysis of mineralization on the nano-scale and use of 3D images at the nano-scale to drive mechanical simulations. Finally, we discuss the issue of going beyond qualitative description to quantification of ultra-structural features.

Organ and tissue level properties are more sensitive to age than osteocyte lacunar characteristics in rat cortical bone

Modeling and remodeling induce significant changes of bone structure and mechanical properties with age. Therefore, it is important to gain knowledge of the processes taking place in bone over time. The rat is a widely used animal model, where much data has been accumulated on age-related changes of bone on the organ and tissue level, whereas features on the nano- and micrometer scale are much less explored. We investigated the age-related development of organ and tissue level bone properties such as bone volume, bone mineral density, and load to fracture and correlated these with osteocyte lacunar properties in rat cortical bone. Femora of 14 to 42-week-old female Wistar rats were investigated using multiple complementary techniques including X-ray micro-computed tomography and biomechanical testing. The body weight, femoral length, aBMD, load to fracture, tissue volume, bone volume, and tissue density were found to increase rapidly with age at 14–30 weeks. At the age of 30–42 weeks, the growth rate appeared to decrease. However, no accompanying changes were found in osteocyte lacunar properties such as lacunar volume, ellipsoidal radii, lacunar stretch, lacunar oblateness, or lacunar orientation with animal age. Hence, the evolution of organ and tissue level properties with age in rat cortical bone is not accompanied by related changes in osteocyte lacunar properties. This suggests that bone microstructure and bone matrix material properties and not the geometric properties of the osteocyte lacunar network are main determinants of the properties of the bone on larger length scales.

Long bone maturation is driven by pore closing: A quantitative tomography investigation of structural formation in young C57BL/6 mice

During mammalian growth, long bones undergo extensive structural reorganization, transforming from primitive shapes in the limb buds into mature bones. Here we shed light on the steps involved in structural formation of the mineralized tissue in midshafts of C57BL/6 femurs, shortly after birth. By combining 3D micrometer-resolution X-ray microtomography with 2D histology, we study the transformation of the tissue from a partially-mineralized scaffold into a compact bone structure. We identify three growth phases that take place during murine long bone maturation: During a patterning phase (I) mineralized struts form a loosely connected foam-like cortical network. During a transitioning phase (II), the extensive non-mineralized tracts vanish, transforming the foam into a fully continuous mass, by 14 days of age. Concomitantly, closed porosity increases to about ∼ 1.4%, and stays at this level, also found in maturity. During a shaping phase (III), the bones gradually attain their characteristic intricate adult form. Architectured mineral depositioning--first in open foamy scaffolds, and later into solid bone material--is presumably a compromise between the mechanical needs of providing support to the body, and the biological requirements of vascularization and extensive nutritional needs in the early stages of bone formation.

High-throughput phenotyping and genetic linkage of cortical bone microstructure in the mouse

Background

Understanding cellular structure and organization, which plays an important role in biological systems ranging from mechanosensation to neural organization, is a complicated multifactorial problem depending on genetics, environmental factors, and stochastic processes. Isolating these factors necessitates the measurement and sensitive quantification of many samples in a reliable, high-throughput, unbiased manner. In this manuscript we present a pipelined approach using a fully automated framework based on Synchrotron-based X-ray Tomographic Microscopy (SRXTM) for performing a full 3D characterization of millions of substructures.

Results

We demonstrate the framework on a genetic study on the femur bones of in-bred mice. We measured 1300 femurs from a F2 cross experiment in mice without the growth hormone (which can confound many of the smaller structural differences between strains) and characterized more than 50 million osteocyte lacunae (cell-sized hollows in the bone). The results were then correlated with genetic markers in a process called quantitative trait localization (QTL). Our findings provide a mapping between regions of the genome (all 19 autosomes) and observable phenotypes which could explain between 8–40 % of the variance using between 2–10 loci for each trait. This map shows 4 areas of overlap with previous studies looking at bone strength and 3 areas not previously associated with bone.

Conclusions

The mapping of microstructural phenotypes provides a starting point for both structure-function and genetic studies on murine bone structure and the specific loci can be investigated in more detail to identify single gene candidates which can then be translated to human investigations. The flexible infrastructure offers a full spectrum of shape, distribution, and connectivity metrics for cellular networks and can be adapted to a wide variety of materials ranging from plant roots to lung tissue in studies requiring high sample counts and sensitive metrics such as the drug-gene interactions and high-throughput screening.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1617-y) contains supplementary material, which is available to authorized users.

Establishing biomechanical mechanisms in mouse models: practical guidelines for systematically evaluating phenotypic changes in the diaphyses of long bones

Mice are widely used in studies of skeletal biology, and assessment of their bones by mechanical testing is a critical step when evaluating the functional effects of an experimental perturbation. For example, a gene knockout may target a pathway important in bone formation and result in a "low bone mass" phenotype. But how well does the skeleton bear functional loads; eg, how much do bones deform during loading and how resistant are bones to fracture? By systematic evaluation of bone morphological, densitometric, and mechanical properties, investigators can establish the "biomechanical mechanisms" whereby an experimental perturbation alters whole-bone mechanical function. The goal of this review is to clarify these biomechanical mechanisms and to make recommendations for systematically evaluating phenotypic changes in mouse bones, with a focus on long-bone diaphyses and cortical bone. Further, minimum reportable standards for testing conditions and outcome variables are suggested that will improve the comparison of data across studies. Basic biomechanical principles are reviewed, followed by a description of the cross-sectional morphological properties that best inform the net cellular effects of a given experimental perturbation and are most relevant to biomechanical function. Although morphology is critical, whole-bone mechanical properties can only be determined accurately by a mechanical test. The functional importance of stiffness, maximum load, postyield displacement, and work-to-fracture are reviewed. Because bone and body size are often strongly related, strategies to adjust whole-bone properties for body mass are detailed. Finally, a comprehensive framework is presented using real data, and several examples from the literature are reviewed to illustrate how to synthesize morphological, tissue-level, and whole-bone mechanical properties of mouse long bones.