Automated compartmental analysis for high-throughput skeletal phenotyping in femora of genetic mouse models

A new method for segmentation and analysis of bone callus in rodent fracture models using micro-CT

Fracture burden has created a need to better understand bone repair processes under different pathophysiological states. Evaluation of structural and material properties of the mineralized callus, which is integral to restoring biomechanical stability is, therefore, vital. Microcomputed tomography (micro-CT) can facilitate noninvasive imaging of fracture repair, however, current methods for callus segmentation are only semiautomated, restricted to defined regions, time/labor intensive, and prone to user variation. Herein, we share a new automatic method for segmenting callus in micro-CT tomograms that will allow for objective, quantitative analysis of the bone fracture microarchitecture. Fractured and nonfractured mouse femurs were scanned and processed by both manual and automated segmentation of fracture callus from cortical bone after which microarchitectural parameters were analyzed. All segmentation and analysis steps were performed using CTAn (Bruker) with automatic segmentation performed using the software's image-processing plugins. Results showed automatic segmentation reliably and consistently segmented callus from cortical bone, demonstrating good agreement with manual methods with low bias: tissue volume (TV): -0.320 mm3 , bone volume (BV): 0.0358 mm3 , and bone volume/tissue volume (BV/TV): -3.52%, and was faster and eliminated user-bias and variation. Method scalability and translatability across rodent models were verified in scans of fractured rat femora showing good agreement with manual methods with low bias: TV: -3.654 mm3 , BV: 0.830 mm3 , and BV/TV: 7.81%. Together, these data validate a new automated method for segmentation of callus and cortical bone in micro-CT tomograms that we share as a fast, reliable, and less user-dependent tool for application to study bone callus in fracture, and potentially elsewhere.

Mechanics and differential healing outcomes of small and large defect injuries of the tendon-bone attachment in the rat rotator cuff

INTRODUCTION

Rotator cuff tear size affects clinical outcomes following rotator cuff repair and is correlated with the risk of recurrent tendon defects. This study aimed to understand if and how the initial defect size influences the structural and mechanical outcomes of the injured rotator cuff attachment in vivo.

METHODS

Full-thickness punch injuries of the infraspinatus tendon-bone attachment in Long Evans rats were created to compare differences in healing outcomes between small and large defects. Biomechanical properties, gross morphology, bone remodeling, and cell and tissue morphology were assessed at both 3- and 8-weeks of healing.

RESULTS

At the time of injury (no healing), large defects had decreased mechanical properties compared to small defects, and both defect sizes had decreased mechanical properties compared to intact attachments. However, the mechanical properties of the two defect groups were not significantly different from each other after 8-weeks of healing and significantly improved compared to no healing but failed to return to intact levels. Local bone volume at the defect site was higher in large compared to small defects on average and increased from 3- to 8-weeks. In contrast, bone quality decreased from 3- to 8-weeks of healing and these changes were not dependent on defect size. Qualitatively, large defects had increased collagen disorganization and neovascularization compared to small defects.

DISCUSSION

In this study, we showed that both large and small defects did not regenerate the mechanical and structural integrity of the intact rat rotator cuff attachment following healing in vivo after 8 weeks of healing.

Micro-computed tomography assessment of bone structure in aging mice

High-resolution computed tomography (CT) is widely used to assess bone structure under physiological and pathological conditions. Although the analytic protocols and parameters for micro-CT (μCT) analyses in mice are standardized for long bones, vertebrae, and the palms in aging mice, they have not yet been established for craniofacial bones. In this study, we conducted a morphometric assessment of craniofacial bones, in comparison with long bones, in aging mice. Although age-related changes were observed in the microarchitecture of the femur, tibia, vertebra, and basisphenoid bone, and were more pronounced in females than in males, the microarchitecture of both the interparietal bone and body of the mandible, which develop by intramembranous ossification, was less affected by age and sex. By contrast, the condyle of the mandible was more affected by aging in males compared to females. Taken together, our results indicate that mouse craniofacial bones are uniquely affected by age and sex.

Growth and development of trabecular structure in the calcaneus of Japanese macaques (Macaca fuscata) reflects locomotor behavior, life history, and neuromuscular development

Bone structure dynamically adapts to its mechanical environment throughout ontogeny by altering the structure of trabecular bone, the three-dimensional mesh-like structure found underneath joint surfaces. Trabecular structure, then, can provide a record of variation in loading directions and magnitude; and in ontogenetic samples, it can potentially be used to track developmental shifts in limb posture. We aim to broaden the analysis of trabecular bone ontogeny by incorporating interactions between ontogenetic variation in locomotor repertoire, neuromuscular maturation, and life history. We examine the associations between these variables and age-related variation in trabecular structure in the calcaneus of Japanese macaques (Macaca fuscata). We used high-resolution micro-computed tomography scanning to image the calcaneus in a cross-sectional sample of 34 juvenile M. fuscata aged between 0 and 7 years old at the Primate Research Institute, Japan. We calculated whole bone averages of standard trabecular properties and generated whole-bone morphometric maps of bone volume fraction and Young's modulus. Trabecular structure becomes increasingly heterogeneous in older individuals. Bone volume fraction (BV/total volume [TV]) decreases during the first month of life and increases afterward, coinciding with the onset of independent locomotion in M. fuscata. At birth, primary Young's modulus is oriented orthogonal to the ossification center, but after locomotor onset bone structure becomes stiffest in the direction of joint surfaces and muscle attachments. Age-related variation in bone volume fraction is best predicted by an interaction between the estimated percentage of adult brain size, body mass, and locomotor onset. To explain our findings, we propose a model where interactions between age-related increases in body weight and maturation of the neuromuscular system alter the loading environment of the calcaneus, to which the internal trabecular structure dynamically adapts. This model cannot be directly tested based on our cross-sectional data. However, confirmation of the model by longitudinal experiments and in multiple species would show that trabecular structure can be used both to infer behavior from fossil morphology and serve as a valuable proxy for neuromuscular maturation and life history events like locomotor onset and the achievement of an adult-like gait. This approach could significantly expand our knowledge of the biology and behavior of fossil species.

Bone Mechanoregulation Allows Subject-Specific Load Estimation Based on Time-Lapsed Micro-CT and HR-pQCT in Vivo

Patients at high risk of fracture due to metabolic diseases frequently undergo long-term antiresorptive therapy. However, in some patients, treatment is unsuccessful in preventing fractures or causes severe adverse health outcomes. Understanding load-driven bone remodelling, i.e., mechanoregulation, is critical to understand which patients are at risk for progressive bone degeneration and may enable better patient selection or adaptive therapeutic intervention strategies. Bone microarchitecture assessment using high-resolution peripheral quantitative computed tomography (HR-pQCT) combined with computed mechanical loads has successfully been used to investigate bone mechanoregulation at the trabecular level. To obtain the required mechanical loads that induce local variances in mechanical strain and cause bone remodelling, estimation of physiological loading is essential. Current models homogenise strain patterns throughout the bone to estimate load distribution in vivo, assuming that the bone structure is in biomechanical homoeostasis. Yet, this assumption may be flawed for investigating alterations in bone mechanoregulation. By further utilising available spatiotemporal information of time-lapsed bone imaging studies, we developed a mechanoregulation-based load estimation (MR) algorithm. MR calculates organ-scale loads by scaling and superimposing a set of predefined independent unit loads to optimise measured bone formation in high-, quiescence in medium-, and resorption in low-strain regions. We benchmarked our algorithm against a previously published load history (LH) algorithm using synthetic data, micro-CT images of murine vertebrae under defined experimental in vivo loadings, and HR-pQCT images from seven patients. Our algorithm consistently outperformed LH in all three datasets. In silico-generated time evolutions of distal radius geometries (n = 5) indicated significantly higher sensitivity, specificity, and accuracy for MR than LH (p < 0.01). This increased performance led to substantially better discrimination between physiological and extra-physiological loading in mice (n = 8). Moreover, a significantly (p < 0.01) higher association between remodelling events and computed local mechanical signals was found using MR [correct classification rate (CCR) = 0.42] than LH (CCR = 0.38) to estimate human distal radius loading. Future applications of MR may enable clinicians to link subtle changes in bone strength to changes in day-to-day loading, identifying weak spots in the bone microstructure for local intervention and personalised treatment approaches.

Locomotory behaviour of early tetrapods from Blue Beach, Nova Scotia, revealed by novel microanatomical analysis

Evidence for terrestriality in early tetrapods is fundamentally contradictory. Fossil trackways attributed to early terrestrial tetrapods long predate the first body fossils from the Late Devonian. However, the Devonian body fossils demonstrate an obligatorily aquatic lifestyle. Complicating our understanding of the transition from water to land is a pronounced gap in the fossil record between the aquatic Devonian taxa and presumably terrestrial tetrapods from the later Early Carboniferous. Recent work suggests that an obligatorily aquatic habit persists much higher in the tetrapod tree than previously recognized. Here, we present independent microanatomical data of locomotor capability from the earliest Carboniferous of Blue Beach, Nova Scotia. The site preserves limb bones from taxa representative of Late Devonian to mid-Carboniferous faunas as well as a rich trackway record. Given that bone remodels in response to functional stresses including gravity and ground reaction forces, we analysed both the midshaft compactness profiles and trabecular anisotropy, the latter using a new whole bone approach. Our findings suggest that early tetrapods retained an aquatic lifestyle despite varied limb morphologies, prior to their emergence onto land. These results suggest that trackways attributed to early tetrapods be closely scrutinized for additional information regarding their creation conditions, and demand an expansion of sampling to better identify the first terrestrial tetrapods.

The application of micro-CT in egg-laying hen bone analysis: introducing an automated bone separation algorithm

The application of micro-CT in small animal research, especially on bone health, has risen exponentially in recent years. However, its application in egg-laying hen bone analysis was still limited. This review introduces the technical aspects of micro-CT in egg-laying hen bone analysis, especially with the medullary bones presented in the cavity. In order to acquate application of micro-CT for laying hen bone research, image acquisition, reconstruction, and analysis settings need to be adjusted properly. The key difference regarding the application of micro-CT in laying hen bone compared to other small animals such as mice and rats was the larger bone size and more complex structures of medullary and trabecular bones. In order to analyze the details of laying hen bone structures, the volume of interest for laying hen should be selected at a region where all 3 bones are present (critical, trabecular, and medullary bone). Owing to the complexity of bone structures, the conventional techniques are not useful to distinguish the trabecular bone and medullary bone in laying hens accurately. In the current review, an automated segmentation algorithm is described to allow researchers to segment bone compartments without human bias. The algorithm is designed according to the morphology difference of medullary bones compared to trabecular and cortical bones. In this procedure, the loosely woven bones were separated by applying dual thresholds. The medullary calcium chunks were separated by opening or closing procedures, where we defined the diameter of medullary chunks being higher than the trabecular bone thickness as a separation trait. The application of micro-CT in laying hen bone health assessment will significantly expand our understanding of chicken bone physiology and osteoporosis, contributing to improve welfare in laying hens.

Bone and joint enhancement filtering: Application to proximal femur segmentation from uncalibrated computed tomography datasets

Methods for reliable femur segmentation enable the execution of quality retrospective studies and building of robust screening tools for bone and joint disease. An enhance-and-segment pipeline is proposed for proximal femur segmentation from computed tomography datasets. The filter is based on a scale-space model of cortical bone with properties including edge localization, invariance to density calibration, rotation invariance, and stability to noise. The filter is integrated with a graph cut segmentation technique guided through user provided sparse labels for rapid segmentation. Analysis is performed on 20 independent femurs. Rater proximal femur segmentation agreement was 0.21 mm (average surface distance), 0.98 (Dice similarity coefficient), and 2.34 mm (Hausdorff distance). Manual segmentation added considerable variability to measured failure load and volume (CV_RMS > 5%) but not density. The proposed algorithm considerably improved inter-rater reproducibility for all three outcomes (CV_RMS < 0.5%). The algorithm localized the periosteal surface accurately compared to manual segmentation but with a slight bias towards a smaller volume. Hessian-based filtering and graph cut segmentation localizes the periosteal surface of the proximal femur with comparable accuracy and improved precision compared to manual segmentation.

Hallmarks of frailty and osteosarcopenia in prematurely aged PolgA(D257A/D257A) mice

BACKGROUND

Frailty is a geriatric syndrome characterized by increased susceptibility to adverse health outcomes. One major determinant thereof is the gradual weakening of the musculoskeletal system and the associated osteosarcopenia. To improve our understanding of the underlying pathophysiology and, more importantly, to test potential interventions aimed at counteracting frailty, suitable animal models are needed.

METHODS

To evaluate the relevance of prematurely aged PolgA(D257A/D257A) mice as a model for frailty and osteosarcopenia, we quantified the clinical mouse frailty index in PolgA(D257A/D257A) and wild-type littermates (PolgA(+/+) , WT) with age and concertedly assessed the quantity and quality of bone and muscle tissue. Lastly, the anabolic responsiveness of skeletal muscle, muscle progenitors, and bone was assessed.

RESULTS

PolgA(D257A/D257A) accumulated health deficits at a higher rate compared with WT, resulting in a higher frailty index at 40 and 46 weeks of age (+166%, +278%, P < 0.0001), respectively, with no differences between genotypes at 34 weeks. Concomitantly, PolgA(D257A/D257A) displayed progressive musculoskeletal deterioration such as reduced bone and muscle mass as well as impaired functionality thereof. In addition to lower muscle weights (-14%, P < 0.05, -23%, P < 0.0001) and fibre area (-20%, P < 0.05, -22%, P < 0.0001) at 40 and 46 weeks, respectively, PolgA(D257A/D257A) showed impairments in grip strength and concentric muscle forces (P < 0.05). PolgA(D257A/D257A) mutation altered the acute response to various anabolic stimuli in skeletal muscle and muscle progenitors. While PolgA(D257A/D257A) muscles were hypersensitive to eccentric contractions as well as leucine administration, shown by larger downstream signalling response of the mechanistic target of rapamycin complex 1, myogenic progenitors cultured in vitro showed severe anabolic resistance to leucine and robust impairments in cell proliferation. Longitudinal micro-computed tomography analysis of the sixth caudal vertebrae showed that PolgA(D257A/D257A) had lower bone morphometric parameters (e.g. bone volume fraction, trabecular, and cortical thickness, P < 0.05) as well as reduced remodelling activities (e.g. bone formation and resorption rate, P < 0.05) compared with WT. When subjected to 4 weeks of cyclic loading, young but not aged PolgA(D257A/D257A) caudal vertebrae showed load-induced bone adaptation, suggesting reduced mechanosensitivity with age.

CONCLUSIONS

PolgA(D257A/D257A) mutation leads to hallmarks of age-related frailty and osteosarcopenia and provides a powerful model to better understand the relationship between frailty and the aging musculoskeletal system.

Mechanical regulation of bone formation and resorption around implants in a mouse model of osteopenic bone

Although mechanical stimulation is considered a promising approach to accelerate implant integration, our understanding of load-driven bone formation and resorption around implants is still limited. This lack of knowledge may delay the development of effective loading protocols to prevent implant loosening, especially in osteoporosis. In healthy bone, formation and resorption are mechanoregulated processes. In the intricate context of peri-implant bone regeneration, it is not clear whether bone (re)modelling can still be load-driven. Here, we investigated the mechanical control of peri-implant bone (re)modelling with a well-controlled mechanobiological experiment. We applied cyclic mechanical loading after implant insertion in tail vertebrae of oestrogen depleted mice and we monitored peri-implant bone response by in vivo micro-CT. Experimental data were combined with micro-finite element simulations to estimate local tissue strains in (re)modelling locations. We demonstrated that a substantial increase in bone mass around the implant could be obtained by loading the entire bone. This augmentation could be attributed to a large reduction in bone resorption rather than to an increase in bone formation. We also showed that following implantation, mechanical regulation of bone (re)modelling was transiently lost. Our findings should help to clarify the role of mechanical stimulation on the maintenance of peri-implant bone mass.

Consensus approach for 3D joint space width of metacarpophalangeal joints of rheumatoid arthritis patients using high-resolution peripheral quantitative computed tomography

Background

Joint space assessment for rheumatoid arthritis (RA) by ordinal conventional radiographic scales is susceptible to floor and ceiling effects. High-resolution peripheral quantitative computed tomography (HR-pQCT) provides superior resolution, and may detect earlier changes. The goal of this work was to compare existing 3D methods to calculate joint space width (JSW) metrics in human metacarpophalangeal (MCP) joints with HR-pQCT and reach consensus for future studies. Using the consensus method, we established reproducibility with repositioning as well as feasibility for use in second-generation HR-pQCT scanners.

Methods

Three published JSW methods were compared using datasets from individuals with RA from three research centers. A SPECTRA consensus method was developed to take advantage of strengths of the individual methods. Using the SPECTRA method, reproducibility after repositioning was tested and agreement between scanner generations was also established.

Results

When comparing existing JSW methods, excellent agreement was shown for JSW minimum and mean (ICC 0.987-0.996) but not maximum and volume (ICC 0.000-0.897). Differences were identified as variations in volume definitions and algorithmic differences that generated high sensitivity to boundary conditions. The SPECTRA consensus method reduced this sensitivity, demonstrating good scan-rescan reliability (ICC >0.911) except for minimum JSW (ICC 0.656). There was strong agreement between results from first- and second-generation HR-pQCT (ICC >0.833).

Conclusions

The SPECTRA consensus method combines unique strengths of three independently-developed algorithms and leverages underlying software updates to provide a mature analysis to measure 3D JSW. This method is robust with respect to repositioning and scanner generations, suggesting its suitability for detecting change.

The association between mineralised tissue formation and the mechanical local in vivo environment: Time-lapsed quantification of a mouse defect healing model

An improved understanding of how local mechanical stimuli guide the fracture healing process has the potential to enhance clinical treatment of bone injury. Recent preclinical studies of bone defect in animal models have used cross-sectional data to examine this phenomenon indirectly. In this study, a direct time-lapsed imaging approach was used to investigate the local mechanical strains that precede the formation of mineralised tissue at the tissue scale. The goal was to test two hypotheses: 1) the local mechanical signal that precedes the onset of tissue mineralisation is higher in areas which mineralise, and 2) this local mechanical signal is independent of the magnitude of global mechanical loading of the tissue in the defect. Two groups of mice with femoral defects of length 0.85 mm (n = 10) and 1.45 mm (n = 9) were studied, allowing for distinct distributions of tissue scale strains in the defects. The regeneration and (re)modelling of mineralised tissue was observed weekly using in vivo micro-computed tomography (micro-CT), which served as a ground truth for resolving areas of mineralised tissue formation. The mechanical environment was determined using micro-finite element analysis (micro-FE) on baseline images. The formation of mineralised tissue showed strong association with areas of higher mechanical strain (area-under-the-curve: 0.91 ± 0.04, true positive rate: 0.85 ± 0.05) while surface based strains could correctly classify 43% of remodelling events. These findings support our hypotheses by showing a direct association between the local mechanical strains and the formation of mineralised tissue.

Virtual supersampling as post-processing step preserves the trabecular bone morphometry in human peripheral quantitative computed tomography scans

In the clinical field of diagnosis and monitoring of bone diseases, high-resolution peripheral quantitative computed tomography (HR-pQCT) is an important imaging modality. It provides a resolution where quantitative bone morphometry can be extracted in vivo on patients. It is known that HR-pQCT provides slight differences in morphometric indices compared to the current standard approach micro-computed tomography (micro-CT). The most obvious reason for this is the restriction of the radiation dose and with this a lower image resolution. With advances in micro-CT evaluation techniques such as patient-specific remodeling simulations or dynamic bone morphometry, a higher image resolution would potentially also allow the application of such novel evaluation techniques to clinical HR-pQCT measurements. Virtual supersampling as post-processing step was considered to increase the image resolution of HR-pQCT scans. The hypothesis was that this technique preserves the structural bone morphometry. Supersampling from 82 μm to virtual 41 μm by trilinear interpolation of the grayscale values of 42 human cadaveric forearms resulted in strong correlations of structural parameters (R²: 0.96–1.00). BV/TV was slightly overestimated (4.3%, R²: 1.00) compared to the HR-pQCT resolution. Tb.N was overestimated (7.47%; R²: 0.99) and Tb.Th was slightly underestimated (-4.20%; R²: 0.98). The technique was reproducible with PE_%CV between 1.96% (SMI) and 7.88% (Conn.D). In a clinical setting with 205 human forearms with or without fracture measured at 82 μm resolution HR-pQCT, the technique was sensitive to changes between groups in all parameters (p < 0.05) except trabecular thickness. In conclusion, we demonstrated that supersampling preserves the bone morphometry from HR-pQCT scans and is reproducible and sensitive to changes between groups. Supersampling can be used to investigate on the resolution dependency of HR-pQCT images and gain more insight into this imaging modality.

Outer-Boundary Assisted Segmentation and Quantification of Trabecular Bones by an Imagej Plugin

Micro-computed tomography (micro-CT) is routinely used to assess bone quantity and trabecular microstructural properties in small animals under different bone loss conditions. However, the standard approach for trabecular analysis of micro-CT images is slice-by-slice semi-automatic hand-contouring, which is labor intensive and error prone. Described here is an efficient method for automatic segmentation of trabecular bones according to the bone's outer boundaries, where trabecular bones can be identified and segmented automatically with accuracy with less operator bias when appropriate segmentation parameters are set. To profile satisfactory segmentation parameters, an image stack of segmentation results is displayed, where all possible combinations of the segmentation parameters are changed one by one in sequence, and segmentation results with associated parameters can easily be visually checked. As a quality-control feature of the plugin, simulated standard objects are quantified where the measured quantities can be compared with theoretical values. Layer-by-layer quantification of trabecular properties and trabecular thicknesses are reported by such a plugin, and the distributions of such properties within the selected regions can be profiled easily. Although layer-by-layer quantification retains more information about trabecular bones and facilitates further statistical analysis of structural changes, such measures are unavailable from the output of current commercial software, where only a single quantified value for each parameter is reported for each sample. Therefore, the described workflows are better approaches for analyzing trabecular bones with accuracy and efficiency.

Impaired bone formation in ovariectomized mice reduces implant integration as indicated by longitudinal in vivo micro-computed tomography

Although osteoporotic bone, with low bone mass and deteriorated bone architecture, provides a less favorable mechanical environment than healthy bone for implant fixation, there is no general agreement on the impact of osteoporosis on peri-implant bone (re)modeling, which is ultimately responsible for the long term stability of the bone-implant system. Here, we inserted an implant in a mouse model mimicking estrogen deficiency-induced bone loss and we monitored with longitudinal in vivo micro-computed tomography the spatio-temporal changes in bone (re)modeling and architecture, considering the separate contributions of trabecular, endocortical and periosteal surfaces. Specifically, 12 week-old C57BL/6J mice underwent OVX/SHM surgery; 9 weeks after we inserted special metal-ceramics implants into the 6^th caudal vertebra and we measured bone response with in vivo micro-CT weekly for the following 6 weeks. Our results indicated that ovariectomized mice showed a reduced ability to increase the thickness of the cortical shell close to the implant because of impaired peri-implant bone formation, especially at the periosteal surface. Moreover, we observed that healthy mice had a significantly higher loss of trabecular bone far from the implant than estrogen depleted animals. Such behavior suggests that, in healthy mice, the substantial increase in peri-implant bone formation which rapidly thickened the cortex to secure the implant may raise bone resorption elsewhere and, specifically, in the trabecular network of the same bone but far from the implant. Considering the already deteriorated bone structure of estrogen depleted mice, further bone loss seemed to be hindered. The obtained knowledge on the dynamic response of diseased bone following implant insertion should provide useful guidelines to develop advanced treatments for osteoporotic fracture fixation based on local and selective manipulation of bone turnover in the peri-implant region.

Role of HTRA1 in bone formation and regeneration: In vitro and in vivo evaluation

The role of mammalian high temperature requirement protease A1 (HTRA1) in somatic stem cell differentiation and mineralized matrix formation remains controversial, having been demonstrated to impart either anti- or pro-osteogenic effects, depending on the in vitro cell model used. The aim of this study was therefore to further evaluate the role of HTRA1 in regulating the differentiation potential and lineage commitment of murine mesenchymal stem cells in vitro, and to assess its influence on bone structure and regeneration in vivo. Our results demonstrated that short hairpin RNA-mediated ablation of Htra1 in the murine mesenchymal cell line C3H10T1/2 increased the expression of several osteogenic gene markers, and significantly enhanced matrix mineralization in response to BMP-2 stimulation. These effects were concomitant with decreases in the expression of chondrogenic gene markers, and increases in adipogenic gene expression and lipid accrual. Despite the profound effects of loss-of-function of HTRA1 on this in vitro osteochondral model, these were not reproduced in vivo, where bone microarchitecture and regeneration in 16-week-old Htra1-knockout mice remained unaltered as compared to wild-type controls. By comparison, analysis of femurs from 52-week-old mice revealed that bone structure was better preserved in Htra1-knockout mice than age-matched wild-type controls. These findings therefore provide additional insights into the role played by HTRA1 in regulating mesenchymal stem cell differentiation, and offer opportunities for improving our understanding of how this multifunctional protease may act to influence bone quality.

Weak Epistasis Generally Stabilizes Phenotypes in a Mouse Intercross

The extent and strength of epistasis is commonly unresolved in genetic studies, and observed epistasis is often difficult to interpret in terms of biological consequences or overall genetic architecture. We investigated the prevalence and consequences of epistasis by analyzing four body composition phenotypes—body weight, body fat percentage, femoral density, and femoral circumference—in a large F2 intercross of B6-lit/lit and C3.B6-lit/lit mice. We used Combined Analysis of Pleiotropy and Epistasis (CAPE) to examine interactions for the four phenotypes simultaneously, which revealed an extensive directed network of genetic loci interacting with each other, circulating IGF1, and sex to influence these phenotypes. The majority of epistatic interactions had small effects relative to additive effects of individual loci, and tended to stabilize phenotypes towards the mean of the population rather than extremes. Interactive effects of two alleles inherited from one parental strain commonly resulted in phenotypes closer to the population mean than the additive effects from the two loci, and often much closer to the mean than either single-locus model. Alternatively, combinations of alleles inherited from different parent strains contribute to more extreme phenotypes not observed in either parental strain. This class of phenotype-stabilizing interactions has effects that are close to additive and are thus difficult to detect except in very large intercrosses. Nevertheless, we found these interactions to be useful in generating hypotheses for functional relationships between genetic loci. Our findings suggest that while epistasis is often weak and unlikely to account for a large proportion of heritable variance, even small-effect genetic interactions can facilitate hypotheses of underlying biology in well-powered studies.

The role of statistical epistasis in the genetic architecture of complex traits has been of great interest to the genetics community since Fisher introduced the concept in 1918. However, assessing epistasis in human and model organism populations has been impeded by limited statistical power. To mitigate this limitation, we analyzed bone and body composition traits in an unusually large mouse intercross population of over 2000 mice, paired with a recently-developed computational approach that leverages information to detect interactions across multiple phenotypes. We discovered a large network of highly significant genetic interactions between variants that influence complex body composition traits. Although epistasis was abundant, the interaction network was dominated by epistasis that stabilizes phenotypes by reducing phenotypic deviation from the parent strains. Nevertheless, the observed network provides an overview of genetic architecture and specific hypotheses of how QTL combine to affect phenotypes. These findings suggest that epistatic effects are generally of lesser magnitude than main QTL effects, and therefore are unlikely to account for major components of variance, but also reinforce genetic interaction analysis as a potent tool for dissecting the biology of complex traits.

Three-Dimensional Quantitative Morphometric Analysis (QMA) for In Situ Joint and Tissue Assessment of Osteoarthritis in a Preclinical Rabbit Disease Model

This work utilises advances in multi-tissue imaging, and incorporates new metrics which define in situ joint changes and individual tissue changes in osteoarthritis (OA). The aims are to (1) demonstrate a protocol for processing intact animal joints for microCT to visualise relevant joint, bone and cartilage structures for understanding OA in a preclinical rabbit model, and (2) introduce a comprehensive three-dimensional (3D) quantitative morphometric analysis (QMA), including an assessment of reproducibility. Sixteen rabbit joints with and without transection of the anterior cruciate ligament were scanned with microCT and contrast agents, and processed for histology. Semi-quantitative evaluation was performed on matching two-dimensional (2D) histology and microCT images. Subsequently, 3D QMA was performed; including measures of cartilage, subchondral cortical and epiphyseal bone, and novel tibio-femoral joint metrics. Reproducibility of the QMA was tested on seven additional joints. A significant correlation was observed in cartilage thickness from matching histology-microCT pairs. The lateral compartment of operated joints had larger joint space width, thicker femoral cartilage and reduced bone volume, while osteophytes could be detected quantitatively. Measures between the in situ tibia and femur indicated an altered loading scenario. High measurement reproducibility was observed for all new parameters; with ICC ranging from 0.754 to 0.998. In conclusion, this study provides a novel 3D QMA to quantify macro and micro tissue measures in the joint of a rabbit OA model. New metrics were established consisting of: an angle to quantitatively measure osteophytes (σ), an angle to indicate erosion between the lateral and medial femoral condyles (ρ), a vector defining altered angulation (λ, α, β, γ) and a twist angle (τ) measuring instability and tissue degeneration between the femur and tibia, a length measure of joint space width (JSW), and a slope and intercept (m, Χ) of joint contact to demonstrate altered loading with disease progression, as well as traditional bone and cartilage and histo-morphometry measures. We demonstrate correlation of microCT and histology, sensitive discrimination of OA change and robust reproducibility.

Reproducibility of compartmental subchondral bone morphometry in the mouse tibiofemoral joint

AIM

Evidence suggests that subchondral bone can be used as a predictor for the onset of osteoarthritis. As such, there is a need to accurately and reproducibly quantify subchondral bone in areas where osteoarthritis develops. In this paper, we present a novel technique for the segmentation of subchondral bone in the tibiofemoral joint and assess the reproducibility of this method with multiple measures and users.

METHODS

The right hind leg of seven C57BL/6 mice were excised and imaged in μCT. The menisci and patella were manually segmented and the image data was Gaussian filtered and binarized. An in-house algorithm was used to generate cortical and epiphyseal volumes of interest and standard morphometric indices for bone were computed. The intraclass correlation coefficient (ICC), absolute precision error (PE(SD)), and precision error as a percentage of the coefficient of variation of the repeated measurements (PE(%CV)) were calculated for each index. Additionally, an inter-user study was performed using the same indices and statistics.

RESULTS

For repeated measures, ICC ranged from 0.869 (cortical bone volume fraction, femur) to 0.994 (degree of anisotropy, femur). Similarly, PE(%CV) ranged from 0.84% (cortical bone volume fraction, femur) to 5.11% (connectivity density, tibia). For repeated users, no effect was seen in the femur with a slight effect in the tibia.

CONCLUSIONS

A novel method for the automatic segmentation of cortical and epiphyseal bone is presented and is shown to be reproducible in C57BL/6 mice. This tool will allow for high-throughput studies of osteoarthritis in animal models.

Quantitative morphometric patterns in cartilage and bone from the humeral heads of end-stage osteoarthritis patients

OBJECTIVE

The purpose of this work is to investigate in a quantitative manner, the gross and regional structural patterns in cartilage and bone from the humeral head of end-stage OA patients, with the goal of identifying patterns of disease. Since the prevalence of primary OA of the shoulder is increasing as the population ages and the non-traumatic degenerative changes leading to this disease are poorly understood, a site-specific morphometric analysis speaks to the structure-function remodelling relationship of the pathological anatomy.

METHODS

Humeral heads were harvested from twenty-one patients undergoing shoulder arthroplasty for end-stage primary OA. The samples were scanned with micro-computed tomography and magnetic resonance imaging (MRI), and registered to provide reconstructed 3D datasets of the cartilage, cortical and trabecular bone tissues. Gross visual examination of the datasets allowed samples to be classified as OA-like, osteoporosis (OP)-like or OA/OP-like.

RESULTS

Volumes of interest (VOI) separating the OA-like samples into five distinct regions showed positive correlations between bone and cartilage morphometric parameters; specifically in areas where more cartilage has been lost, the underlying subchondral cortical bone was more porous and thicker, while the subchondral trabecular bone was more dense, including more connections and trabeculae. These differences were site-specific, where the central humeral head saw the greatest destruction of cartilage and bone sclerosis, followed by the anterior aspects.

CONCLUSION

The ability to correlate bone and cartilage changes is valuable, as these structural cues may allow a more targeted diagnostic approach and a better classification of the disease, leading to improved therapies.

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.

Development and optimization of a high-throughput micro-computed tomography imaging method incorporating a novel analysis technique to evaluate bone mineral density of arthritic joints in a rodent model of collagen induced arthritis

Rheumatoid arthritis (RA) is a chronic autoimmune disease resulting in joint inflammation, pain, and eventual bone loss. Bone loss and remodeling caused by symmetric polyarthritis, the hallmark of RA, is readily detectable by bone mineral density (BMD) measurement using micro-CT. Abnormalities in these measurements over time reflect the underlying pathophysiology of the bone. To evaluate the efficacy of anti-rheumatic agents in animal models of arthritis, we developed a high throughput knee and ankle joint imaging assay to measure BMD as a translational biomarker. A bone sample holder was custom designed for micro-CT scanning, which significantly increased assay throughput. Batch processing 3-dimensional image reconstruction, followed by automated image cropping, significantly reduced image processing time. In addition, we developed a novel, automated image analysis method to measure BMD and bone volume of knee and ankle joints. These improvements significantly increased the throughput of ex vivo bone sample analysis, reducing data turnaround from 5 days to 24 hours for a study with 200 rat hind limbs. Taken together, our data demonstrate that BMD, as quantified by micro-CT, is a robust efficacy biomarker with a high degree of sensitivity. Our innovative approach toward evaluation of BMD using optimized image acquisition and novel image processing techniques in preclinical models of RA enables high throughput assessment of anti-rheumatic agents offering a powerful tool for drug discovery.

Longitudinal assessment of oxytocin efficacy on bone and bone marrow fat masses in a rabbit osteoporosis model through 3.0-T magnetic resonance spectroscopy and micro-CT

This study aims to longitudinally assess the effect of oxytocin on bone and bone fat masses using micro-CT, in vivo magnetic resonance spectroscopy (MRS), and histopathological adipocyte quantification. Early in vivo oxytocin (OT) treatment to the osteoporosis (OP) rabbit model may reliably inhibit bone degeneration and reduce bone marrow fat accumulation by decreasing marrow adipocyte size and density.

INTRODUCTION

This study aims to longitudinally assess the effect of early OT treatment on bone and bone fat masses in a rabbit OP model by comparing the results of MRS and micro-CT with histopathological findings.

METHODS

Sixty 20-week-old female rabbits were randomly assigned into three groups. The control and OP groups were subjected to either sham surgery or bilateral ovariectomy (OVX). The OT group was subcutaneously injected with OT daily from the second week after OVX for 8 weeks. The left proximal femurs of the rabbits were evaluated through MRS, micro-CT, and histopathological examination at 0, 4, 8, 10, and 12 weeks after operation. Differences in fat fraction (FF) values, micro-CT parameters, and calculated pathological marrow adipocytes among three groups were analyzed.

RESULTS

The FF values of the OP group significantly increased (p = 0.019), but the tissue mineral density (TMD) decreased (p = 0.037) from eighth week compared with those of the control group. The FF values of the OT group significantly decreased (p = 0.044), but the TMD values increased (p = 0.042) from eighth week compared with those of the OP group. The adypocyte diameter of the OT group significantly decreased (p = 0.041) from eighth week and then adypocyte density did so too from tenth week, compared with those of the OP group at the same time point. No difference in adypocyte calculation was found between the OT and control groups until the 12th week after operation.

CONCLUSION

Early in vivo OT treatment slowed down bone deterioration and reduced bone marrow adiposity accumulation in a rabbit OP model, which is consistent with pathologic findings. OT treatment is a promising preventive OP therapy.

Registration of 2D histological sections with 3D micro-CT datasets from small animal vertebrae and tibiae

The aim of this study was the registration of digitized thin 2D sections of mouse vertebrae and tibiae used for histomorphometry of trabecular bone structure into 3D micro computed tomography (μCT) datasets of the samples from which the sections were prepared. Intensity-based and segmentation-based registrations (SegRegs) of 2D sections and 3D μCT datasets were applied. As the 2D sections were deformed during their preparation, affine registration for the vertebrae was used instead of rigid registration. Tibiae sections were additionally cut on the distal end, which subsequently undergone more deformation so that elastic registration was necessary. The Jaccard distance was used as registration quality measure. The quality of intensity-based registrations and SegRegs was practically equal, although precision errors of the elastic registration of segmentation masks in tibiae were lower, while those in vertebrae were lower for the intensity-based registration. Results of SegReg significantly depended on the segmentation of the μCT datasets. Accuracy errors were reduced from approximately 64% to 42% when applying affine instead of rigid transformations for the vertebrae and from about 43% to 24% when using B-spline instead of rigid transformations for the tibiae. Accuracy errors can also be caused by the difference in spatial resolution between the thin sections (pixel size: 7.25 μm) and the μCT data (voxel size: 15 μm). In the vertebrae, average deformations amounted to a 6.7% shortening along the direction of sectioning and a 4% extension along the perpendicular direction corresponding to 0.13-0.17 mm. Maximum offsets in the mouse tibiae were 0.16 mm on average.

Therapeutic potential of adipose-derived stromal cells in age-related osteoporosis

Adipose-derived stromal cells (ASCs) are increasingly being used for orthopedic-based tissue engineering approaches due to their ability to readily undergo osteogenic differentiation. In the present study, we used in vitro and in vivo approaches to evaluate the use of ASCs as a treatment strategy for age-related osteoporosis. Molecular, histological and micro-computed tomography (micro-CT) based approaches confirmed that ASCs isolated from 18-week-old osteoporotic senescence-accelerated mice (SAMP6) were capable of undergoing osteogenesis when cultured in either silk fibroin (SF) scaffolds or scaffold-free microtissues (ASC-MT). A single intratibial injection of CM-Dil-labeled isogeneic ASCs or ASC-MT into SAMP6 recipients significantly improved trabecular bone quality after 6 weeks in comparison to untreated contralateral bones, as determined by micro-CT. Injected ASCs could be observed in paraffin wax bone sections at 24 h and 6 weeks post treatment and induced a significant increase in several molecular markers of bone turnover. Furthermore, a significant improvement in the osteogenic potential of osteoporotic patient-derived human bone marrow stromal cells (BMSCs) was observed when differentiated in conditioned culture media harvested from osteoporotic patient-derived human ASCs. These findings therefore support the use of ASCs as an autologous cell-based approach for the treatment of osteoporosis.

The Clinical Biomechanics Award 2012 - presented by the European Society of Biomechanics: large scale simulations of trabecular bone adaptation to loading and treatment

BACKGROUND

Microstructural simulations of bone remodeling are particularly relevant in the clinical management of osteoporosis. Before a model can be applied in the clinics, a validation against controlled in vivo data is crucial. Here we present a strain-adaptive feedback algorithm for the simulation of trabecular bone remodeling in response to loading and pharmaceutical treatment and report on the results of the large-scale validation against in vivo data.

METHODS

The algorithm follows the mechanostat principle and incorporates mechanical feedback, based on the local strain-energy density. For the validation, simulations of bone remodeling and adaptation in 180 osteopenic mice were performed. Permutations of the conditions for early (20th week) and late (26th week) loading of 8N or 0N, and treatments with bisphosphonates, or parathyroid hormone were simulated. Static and dynamic morphometry and local remodeling sites from in vivo and in silico studies were compared.

FINDINGS

For each study an individual set of model parameters was selected. Trabecular bone volume fraction was chosen as an indicator of the accuracy of the simulations. Overall errors for this parameter were 0.1-4.5%. Other morphometric indices were simulated with errors of less than 19%. Dynamic morphometry was more difficult to predict, which resulted in significant differences from the experimental data.

INTERPRETATION

We validated a new algorithm for the simulation of bone remodeling in trabecular bone. The results indicate that the simulations accurately reflect the effects of treatment and loading seen in respective experimental data, and, following adaptation to human data, could be transferred into clinics.

The role of the renal ammonia transporter Rhcg in metabolic responses to dietary protein

High dietary protein imposes a metabolic acid load requiring excretion and buffering by the kidney. Impaired acid excretion in CKD, with potential metabolic acidosis, may contribute to the progression of CKD. Here, we investigated the renal adaptive response of acid excretory pathways in mice to high-protein diets containing normal or low amounts of acid-producing sulfur amino acids (SAA) and examined how this adaption requires the RhCG ammonia transporter. Diets rich in SAA stimulated expression of enzymes and transporters involved in mediating NH4 (+) reabsorption in the thick ascending limb of the loop of Henle. The SAA-rich diet increased diuresis paralleled by downregulation of aquaporin-2 (AQP2) water channels. The absence of Rhcg transiently reduced NH4 (+) excretion, stimulated the ammoniagenic pathway more strongly, and further enhanced diuresis by exacerbating the downregulation of the Na(+)/K(+)/2Cl(-) cotransporter (NKCC2) and AQP2, with less phosphorylation of AQP2 at serine 256. The high protein acid load affected bone turnover, as indicated by higher Ca(2+) and deoxypyridinoline excretion, phenomena exaggerated in the absence of Rhcg. In animals receiving a high-protein diet with low SAA content, the kidney excreted alkaline urine, with low levels of NH4 (+) and no change in bone metabolism. Thus, the acid load associated with high-protein diets causes a concerted response of various nephron segments to excrete acid, mostly in the form of NH4 (+), that requires Rhcg. Furthermore, bone metabolism is altered by a high-protein acidogenic diet, presumably to buffer the acid load.

Mineralization kinetics in murine trabecular bone quantified by time-lapsed in vivo micro-computed tomography

Trabecular bone is a highly dynamic tissue due to bone remodeling, mineralization and demineralization. The mineral content and its spatial heterogeneity are main contributors to bone quality. Using time-lapsed in vivo micro-computed tomography (micro-CT), it is now possible to resolve in three dimensions where bone gets resorbed and deposited over several weeks. In addition, the gray values in the micro-CT images contain quantitative information about the local tissue mineral density (TMD). The aim of this study was to measure how TMD increases with time after new bone formation and how this mineralization kinetics is influenced by mechanical stimulation. Our analysis of changes in TMD was based on an already reported experiment on 15-week-old female mice (C57BL/6), where in one group the sixth caudal vertebra was mechanically loaded with 8N, while in the control group no loading was applied. Comparison of two consecutive images allows the categorization of bone into newly formed, resorbed, and quiescent bone for different time points. Gray values of bone in these categories were compared layer-wise to minimize the effects of beam hardening artifacts. Quiescent bone in the control group was found to mineralize with a rate of 8 ± 1 mgHA/cm(3) per week, which is about half as fast as observed for newly formed bone. Mechanical loading increased the rate of mineral incorporation by 63% in quiescent bone. The week before bone resorption, demineralization could be observed with a drop of TMD by 36 ± 4 mgHA/cm(3) in the control and 34 ± 3 mgHA/cm(3) in the loaded group. In conclusion, this study shows how time-lapsed in vivo micro-CT can be used to assess changes in TMD of bone with high spatial and temporal resolution. This will allow a quantification of how bone diseases and pharmaceutical interventions influence not only microarchitecture of trabecular bone, but also its material quality.

Osteocyte-derived insulin-like growth factor I is essential for determining bone mechanosensitivity

This study sought to determine whether deficient Igf1 expression in osteocytes would affect loading-induced osteogenic response. Tibias of osteocyte Igf1 conditional knockout (KO) mice (generated by cross-breeding Igf1 floxed mice with Dmp1-Cre transgenic mice) and wild-type (WT) littermates were subjected to four-point bending for 2 wk. Microcomputed tomography confirmed that the size of tibias of conditional mutants was smaller. Loading with an equivalent loading strain increased periosteal woven bone and endosteal lamellar bone formation in WT mice but not in conditional KO mice. Consistent with the lack of an osteogenic response, the loading failed to upregulate expression of early mechanoresponsive genes (Igf1, Cox-2, c-fos) or osteogenic genes (Cbfa-1, and osteocalcin) in conditional KO bones. The lack of osteogenic response was not due to reduced osteocyte density or insufficient loading strain. Deficient osteocyte Igf1 expression reduced the loading-induced upregulation of expression of canonical Wnt signaling genes (Wnt10b, Lrp5, Dkk1, sFrp2). The loading also reduced (by 40%) Sost expression in WT mice, but the loading not only did not reduce but upregulated (~1.5-fold) Sost expression in conditional KO mice. Conditional disruption of Igf1 in osteocytes also abolished the loading-induced increase in the bone β-catenin protein level. These findings suggest an impaired response in the loading-induced upregulation of the Wnt signaling in conditional KO mice. In summary, conditional disruption of Igf1 in osteocytes abolished the loading-induced activation of the Wnt signaling and the corresponding osteogenic response. In conclusion, osteocyte-derived IGF-I plays a key determining role in bone mechanosensitivity.

High-throughput quantification of the mechanical competence of murine femora--a highly automated approach for large-scale genetic studies

Animal models are widely used to gain insight into the role of genetics on bone structure and function. One of the main strategies to map the genes regulating specific traits is called quantitative trait loci (QTL) analysis, which generally requires a very large number of animals (often more than 1000) to reach statistical significance. QTL analysis for mechanical traits has been mainly based on experimental mechanical testing, which, in view of the large number of animals, is time consuming. Hence, the goal of the present work was to introduce an automated method for large-scale high-throughput quantification of the mechanical properties of murine femora. Specifically, our aims were, first, to develop and validate an automated method to quantify murine femoral bone stiffness. Second, to test its high-throughput capabilities on murine femora from a large genetic study, more specifically, femora from two growth hormone (GH) deficient inbred strains of mice (B6-lit/lit and C3.B6-lit/lit) and their first (F1) and second (F2) filial offsprings. Automated routines were developed to convert micro-computed tomography (micro-CT) images of femora into micro-finite element (micro-FE) models. The method was experimentally validated on femora from C57BL/6J and C3H/HeJ mice: for both inbred strains the micro-FE models closely matched the experimentally measured bone stiffness when using a single tissue modulus of 13.06 GPa. The mechanical analysis of the entire dataset (n=1990) took approximately 44 CPU hours on a supercomputer. In conclusion, our approach, in combination with QTL analysis could help to locate genes directly involved in controlling bone mechanical competence.

Decreased stability of erythroblastic islands in integrin β3-deficient mice

Erythroblasts proliferate and differentiate in hematopoietic organs within erythroblastic islands (EI) composed of erythropoietic progenitor cells attached to a central macrophage. This cellular interaction crucially involves the erythroid intercellular adhesion molecule-4 (ICAM-4) and αv integrin. Because integrins are biologically active as α/β heterodimers, we asked whether β3 could be a heterodimerization partner of αv integrin in EIs. To this end we compared stress erythropoiesis driven by two different mechanisms, namely that of integrin β3-deficient (β3(-/-)) mice that exhibit impaired hemostasis due to platelet dysfunction with that of systemically erythropoietin-overexpressing (tg6) mice. While compared to the respective wild type (wt) controls β3(-/-) mice had much less erythropoietic stimulation than tg6 mice β3(-/-) blood contained more erythrocytes of a lower maturity stage. Unexpectedly, membranes of peripheral erythrocytes from β3(-/-) mice (but not those from either wt control or from tg6 mice) contained calnexin, a chaperone that is normally completely lost during terminal differentiation of reticulocytes prior to their release into the circulation. In contrast to erythropoietin-overexpressing mice, the erythropoietic subpopulations representing orthochromatic erythroblasts and premature reticulocytes as well as the number of cells per EI were reduced in β3(-/-) bone marrow. In conclusion, absence of integrin β3 impairs adhesion of the latest erythroid developmental stage to the central macrophage of EIs resulting in preterm release of abnormally immature erythrocytes into the circulation.

Strain-adaptive in silico modeling of bone adaptation--a computer simulation validated by in vivo micro-computed tomography data

Computational models are an invaluable tool to test different mechanobiological theories and, if validated properly, for predicting changes in individuals over time. Concise validation of in silico models, however, has been a bottleneck in the past due to a lack of appropriate reference data. Here, we present a strain-adaptive in silico algorithm which is validated by means of experimental in vivo loading data as well as by an in vivo ovariectomy experiment in the mouse. The maximum prediction error following four weeks of loading resulted in 2.4% in bone volume fraction (BV/TV) and 8.4% in other bone structural parameters. Bone formation and resorption rate did not differ significantly between experiment and simulation. The spatial distribution of formation and resorption sites matched in 55.4% of the surface voxels. Bone loss was simulated with a maximum prediction error of 12.1% in BV/TV and other bone morphometric indices, including a saturation level after a few weeks. Dynamic rates were more difficult to be accurately predicted, showing evidence for significant differences between simulation and experiment (p<0.05). The spatial agreement still amounted to 47.6%. In conclusion, we propose a computational model which was validated by means of experimental in vivo data. The predictive value of an in silico model may become of major importance if the computational model should be applied in clinical settings to predict bone changes due to disease and test the efficacy of potential pharmacological interventions.

Textural characteristics of model and natural bone tissues and interfacial behavior of bound water

Water, as a probe liquid bound in model systems (highly disperse hydroxyapatite - protein composites as a model of the main components of bones) and rat bone tissues healthy and affected by osteoporosis occurred due to experimental Alzheimer's disease (EAD), has been investigated using low-temperature (1)H NMR spectroscopy, NMR cryoporometry, TG/DTG/DTA, DSC, and TG and DSC thermoporometry. The textural characteristics of these intact systems cannot be studied using the standard adsorption methods, but the cryoporometry and thermoporometry methods give these characteristics. The (1)H NMR spectra of water bound in model and natural bone tissues include signals, which can be assigned to strongly associated (typical) water (SAW, chemical shift of proton resonance δ(H)=5-6 ppm) and weakly associated (atypical) water (WAW) at δ(H)=1-2 ppm. Contributions of SAW and WAW give information on textural organization of both model and natural bones. The influence of such co-adsorbates as HCl, CDCl(3), CD(3)CN, C(6)D(6), and (CD(3))(2)SO on the interfacial behavior and clustering of bound water depends on their polarity, amounts of components, and textural and structural features of the materials analyzed with the (1)H NMR spectroscopy and cryoporometry methods. According to the NMR cryoporometry data, the EAD causes an increase in nanoporosity of the bone tissues. The total porosity and the specific surface area of biostructures (accessible for water molecules and estimated using NMR cryoporometry and TG thermoporometry methods with a model of cylindrical pores) are larger for the EAD sample. Weakly polar chloroform-d has a significant influence on the organization of water in the bone tissue, and this effect is greater for the EAD sample as more porous material.

Image interpolation allows accurate quantitative bone morphometry in registered micro-computed tomography scans

Time-lapsed in vivo micro-computed tomography is a powerful tool to analyse longitudinal changes in the bone micro-architecture. Registration can overcome problems associated with spatial misalignment between scans; however, it requires image interpolation which might affect the outcome of a subsequent bone morphometric analysis. The impact of the interpolation error itself, though, has not been quantified to date. Therefore, the purpose of this ex vivo study was to elaborate the effect of different interpolator schemes [nearest neighbour, tri-linear and B-spline (BSP)] on bone morphometric indices. None of the interpolator schemes led to significant differences between interpolated and non-interpolated images, with the lowest interpolation error found for BSPs (1.4%). Furthermore, depending on the interpolator, the processing order of registration, Gaussian filtration and binarisation played a role. Independent from the interpolator, the present findings suggest that the evaluation of bone morphometry should be done with images registered using greyscale information.

White spotting variant mouse as an experimental model for ovarian aging and menopausal biology

OBJECTIVE

Menopause is a unique phenomenon in modern women, as most mammalian species possess a reproductive period comparable with their life span. Menopause is caused by the depletion of germ cell-containing ovarian follicles and in laboratory studies is usually modeled in animals in which the ovarian function is removed through ovariectomy or chemical poisoning of the germ cells. Our objective was to explore and characterize the white spotting variant (Wv) mice that have reduced ovarian germ cell abundance, a result of a point mutation in the c-kit gene that decreases kinase activity, as a genetic model for use in menopause studies.

METHODS

Physiological and morphological features associated with menopause were determined in female Wv/Wv mice compared with age-matched wildtype controls. Immunohistochemistry was used to evaluate the presence and number of follicles in paraffin-embedded ovaries. Bone density and body composition were evaluated using the PIXImus x-ray densitometer, and lipids, calcium, and hormone levels were determined in serum using antigen-specific enzyme immunoassays. Heart and body weight were measured, and cardiac function was evaluated using transthoracic echocardiography.

RESULTS

The ovaries of the Wv/Wv females have a greatly reduced number of normal germ cells at birth compared with wildtype mice. The remaining follicles are depleted by around 2 months, and the ovaries develop benign epithelial lesions that resemble morphological changes that occur during ovarian aging, whereas a normal mouse ovary has numerous follicles at all stages of development and retains some follicles even in advanced age. Wv mice have elevated plasma gonadotropins and reduced estrogen and progesterone levels, a significant reduction in bone mass density, and elevated serum cholesterol and lipoprotein levels. Moreover, the Wv female mice have enlarged hearts and reduced cardiac function.

CONCLUSIONS

The reduction of c-kit activity in Wv mice leads to a substantially diminished follicular endowment in newborn mice and premature depletion of follicles in young mice, although mutant females have a normal life span after cessation of ovarian function. The Wv female mice exhibit consistent physiological changes that resemble common features of postmenopausal women. These alterations include follicle depletion, morphological aging of the ovary, altered serum levels of cholesterol, gonadotropins and steroid hormones, decreased bone density, and reduced cardiac function. These changes were not observed in male mice, either age-matched male Wv/Wv or wildtype mice, and are improbably caused by global loss of c-kit function. The Wv mouse may be a genetic, intact-ovary model that mimics closely the phenotypes of human menopause to be used for further studies to understand the mechanisms of menopausal biology.

The different contributions of cortical and trabecular bone to implant anchorage in a human vertebra

The quality of the peri-implant bone and the strength of the bone-implant interface are important factors for implant anchorage. With regard to peri-implant bone, cortical and trabecular compartments both contribute to the load transfer from the implant to the surrounding bone but their relative roles have yet to be investigated in detail. However, this knowledge is crucial for the better understanding of implant failure and for the development of new implants. This is especially true for osteoporotic bone, which is characterized by a deterioration of the trabecular architecture and a thinning of the cortical shell, leading to a higher probability of implant loosening. The aim of this study was to investigate the relative biomechanical roles of cortical and trabecular bone on implant pull-out stiffness in human vertebrae. The starting point of our investigation was a micro-computed tomography scan of an adult human vertebra. The cortical shell was identified and an implant was digitally inserted into the vertebral body. Pull-out tests were simulated with micro-finite element analysis and the apparent stiffness of the system with various degrees of shell thickness and bone volume fraction was computed. Our computational models demonstrated that cortical bone, although being very thin, plays a major role in the mechanical competence of the bone-implant construct.

Mouse tail vertebrae adapt to cyclic mechanical loading by increasing bone formation rate and decreasing bone resorption rate as shown by time-lapsed in vivo imaging of dynamic bone morphometry

It is known that mechanical loading leads to an increase in bone mass through a positive shift in the balance between bone formation and bone resorption. How the remodeling sites change over time as an effect of loading remains, however, to be clarified. The purpose of this paper was to investigate how bone formation and resorption sites are modulated by mechanical loading over time by using a new imaging technique that extracts three dimensional formation and resorption parameters from time-lapsed in vivo micro-computed tomography images. To induce load adaptation, the sixth caudal vertebra of C57BL/6 mice was cyclically loaded through pins in the adjacent vertebrae at either 8 N or 0 N (control) three times a week for 5 min (3000 cycles) over a total of 4 weeks. The results showed that mechanical loading significantly increased trabecular bone volume fraction by 20% (p<0.001) and cortical area fraction by 6% (p<0.001). The bone formation rate was on average 23% greater (p<0.001) and the bone resorption rate was on average 25% smaller (p<0.001) for the 8 N group than for the 0 N group. The increase in bone formation rate for the 8 N group was mostly an effect of a significantly increased surface of bone formation sites (on average 16%, p<0.001), while the thickness of bone formation packages was less affected (on average 5% greater, p<0.05). At the same time the surface of bone resorption sites was significantly reduced (on average 15%, p<0.001), while the depth of resorption pits remained the same. For the 8 N group, the strength of the whole bone increased significantly by 24% (p<0.001) over the loading period, while the strain energy density in the trabecular bone decreased significantly by 24% (p<0.001). In conclusion, mouse tail vertebrae adapt to mechanical loading by increasing the surface of formation sites and decreasing the surface of resorption sites, leading to an overall increase in bone strength. This new imaging technique will provide opportunities to investigate in vivo bone remodeling in the context of disease and treatment options, with the added value that both bone formation and bone resorption parameters can be nondestructively calculated over time.

The importance of murine cortical bone microstructure for microcrack initiation and propagation

In order to better understand bone postyield behavior and consequently bone failure behavior, this study aimed first to investigate cortical bone microstructure and second, to relate cortical bone microstructure to microdamage initiation and propagation in C57BL/6 (B6) and C3H/He (C3H) mice; two murine inbred strains known for their differences in bone phenotype. Murine femora of B6 and C3H were loaded axially under compression in a stepwise manner. For each loading step, 3D data sets at a nominal resolution of 700 nm were acquired by means of synchrotron radiation-based computed tomography. Cortical bone microstructure was divided into three phases: the canal network, the osteocyte lacunar system, and microdamage. Canal volume density and canal unit volume both correlated highly to crack number density (canal volume density: R(2)=0.64, p<0.005 and canal unit volume: R(2)=0.75, p<0.001). Moreover, the large canal units in C3H bone were responsible for more microdamage accumulation compared to B6 bones. This more pronounced microdamage accumulation due to large intracortical bone voids, which eventually leads to a fatal macrocrack (fracture), represents a potential contributing factor to the higher incidence of bone fractures in the elderly.

In vivo micro-computed tomography allows direct three-dimensional quantification of both bone formation and bone resorption parameters using time-lapsed imaging

Bone is a living tissue able to adapt its structure to external influences such as altered mechanical loading. This adaptation process is governed by two distinct cell types: bone-forming cells called osteoblasts and bone-resorbing cells called osteoclasts. It is therefore of particular interest to have quantitative access to the outcomes of bone formation and resorption separately. This article presents a non-invasive three-dimensional technique to directly extract bone formation and resorption parameters from time-lapsed in vivo micro-computed tomography scans. This includes parameters such as Mineralizing Surface (MS), Mineral Apposition Rate (MAR), and Bone Formation Rate (BFR), which were defined in accordance to the current nomenclature of dynamic histomorphometry. Due to the time-lapsed and non-destructive nature of in vivo micro-computed tomography, not only formation but also resorption can now be assessed quantitatively and time-dependent parameters Eroded Surface (ES) as well as newly defined indices Mineral Resorption Rate (MRR) and Bone Resorption Rate (BRR) are introduced. For validation purposes, dynamic formation parameters were compared to the traditional quantitative measures of dynamic histomorphometry, where MAR correlated with R = 0.68 and MS with R = 0.78 (p < 0.05). Reproducibility was assessed in 8 samples that were scanned 5 times and errors ranged from 0.9% (MRR) to 6.6% (BRR). Furthermore, the new parameters were applied to a murine in vivo loading model. A comparison of directly extracted parameters between formation and resorption within each animal revealed that in the control group, i.e., during normal remodeling, MAR was significantly lower than MRR (p < 0.01), whereas MS compared to ES was significantly higher (p < 0.0001). This implies that normal remodeling seems to take place by many small formation packets and few but large resorption volumes. After 4 weeks of mechanical loading, newly extracted trabecular BFR and MS were significantly higher (p < 0.01) in the loading compared to the control group. At the same time, ES was significantly decreased (p < 0.01). This indicates that modeling induced by mechanical loading takes place primarily by increased area, not width of formation packets. With these results, we conclude that the non-invasive direct technique is well suited to extract dynamic bone morphometry parameters and eventually gain more insight into the processes of bone adaptation not only for formation but also resorption.

High-throughput full-automatic synchrotron-based tomographic microscopy

At the TOMCAT (TOmographic Microscopy and Coherent rAdiology experimenTs) beamline of the Swiss Light Source with an energy range of 8-45 keV and voxel size from 0.37 µm to 7.4 µm, full tomographic datasets are typically acquired in 5 to 10 min. To exploit the speed of the system and enable high-throughput studies to be performed in a fully automatic manner, a package of automation tools has been developed. The samples are automatically exchanged, aligned, moved to the correct region of interest, and scanned. This task is accomplished through the coordination of Python scripts, a robot-based sample-exchange system, sample positioning motors and a CCD camera. The tools are suited for any samples that can be mounted on a standard SEM stub, and require no specific environmental conditions. Up to 60 samples can be analyzed at a time without user intervention. The throughput of the system is dependent on resolution, energy and sample size, but rates of four samples per hour have been achieved with 0.74 µm voxel size at 17.5 keV. The maximum intervention-free scanning time is theoretically unlimited, and in practice experiments have been running unattended as long as 53 h (the average beam time allocation at TOMCAT is 48 h per user). The system is the first fully automated high-throughput tomography station: mounting samples, finding regions of interest, scanning and reconstructing can be performed without user intervention. The system also includes many features which accelerate and simplify the process of tomographic microscopy.

Genetic loci for bone architecture determined by three-dimensional CT in crosses with the diabetic GK rat

The F344 rat carries alleles contributing to bone fragility while the GK rat spontaneously develops type-2 diabetes. These characteristics make F344×GK crosses well suited for the identification of genes related to bone size and allow for future investigation on the association with type-2 diabetes. The aim of this study was to identify quantitative trait loci (QTLs) for bone size phenotypes measured by a new application of three-dimensional computed tomography (3DCT) and to investigate the effects of sex- and reciprocal cross. Tibia from male and female GK and F344 rats, representing the parental, F1 and F2 generations, were examined with 3DCT and analyzed for: total and cortical volumetric BMD, straight and curved length, peri- and endosteal area at mid-shaft. F2 progeny (108 male and 98 female) were genotyped with 192 genome-wide microsatellite markers (average distance 10 cM). Sex- and reciprocal cross-separated QTL analyses were performed for the identification of QTLs linked to 3DCT phenotypes and true interactions were confirmed by likelihood ratio analysis in all F2 animals. Several genome-wide significant QTLs were found in the sex- and reciprocal cross-separated progeny on chromosomes (chr) 1, 3, 4, 9, 10, 14, and 17. Overlapping QTLs for both males and females in the (GK×F344)F2 progeny were located on chr 1 (39-67 cM). This region confirms previously reported pQCT QTLs and overlaps loci for fasting glucose. Sex separated linkage analysis confirmed a male specific QTL on chr 9 (67-82 cM) for endosteal area at the fibula site. Analyses separating the F2 population both by sex and reciprocal cross identified cross specific QTLs on chr 14 (males) and chr 3 and 4 (females). Two loci, chr 4 and 6, are unique to 3DCT and separate from pQCT generated loci. The 3DCT method was highly reproducible and provided high precision measurements of bone size in the rat enabling identification of new sex- and cross-specific loci. The QTLs on chr 1 indicate potential genetic association between bone-related phenotypes and traits affecting type-2 diabetes. The results illustrate the complexity of the genetic architecture of bone size phenotypes and demonstrate the importance of complementary methods for bone analysis.

Rapid in vivo whole body composition of rats using cone beam μCT

Precise, noninvasive analysis and quantification of in vivo body composition is essential for research involving longitudinal, small-animal disease models. We investigated the feasibility and precision of a rapid, flat-panel μCT scanner to report whole body adipose tissue volume (ATV), lean tissue volume (LTV), skeletal tissue volume (STV), and bone mineral content (BMC) in 25 postmortem female and 52 live male Sprague-Dawley rats. μCT images, acquired in three 90-mm segments and reconstructed with 308 μm of isotropic voxel spacing, formed contiguous image volumes of each entire rat specimen. Three signal-intensity thresholds (determined to be -186, 5, and 155 HU) were used to classify each voxel as adipose, lean, or skeletal tissue, respectively. Tissue masses from the volume fractions of ATV, LTV, and STV were calculated from assumed tissue densities of 0.95, 1.05, and 1.92 g/cm(-3), respectively. A CT-derived total mass was calculated for each rat and compared with the gravimetrically measured mass, which differed on average for the postmortem female and the live male group by 2.5 and 1.1%, respectively. To evaluate the accuracy of the CT-derived body composition technique, following the live male study excised muscle tissue in the lower right leg of all rats in group B were compared with the image-derived LT measurement of the same regional compartment and found to differ on average by 2.2%. Through repeated CT measurements of postmortem specimens, the whole body ATV, LTV, STV, and BMC measurement analysis gave a precision value of ±0.6, 1.9, 1.7, and 0.5% of the average value, respectively.

Guidelines for assessment of bone microstructure in rodents using micro-computed tomography

Use of high-resolution micro-computed tomography (microCT) imaging to assess trabecular and cortical bone morphology has grown immensely. There are several commercially available microCT systems, each with different approaches to image acquisition, evaluation, and reporting of outcomes. This lack of consistency makes it difficult to interpret reported results and to compare findings across different studies. This article addresses this critical need for standardized terminology and consistent reporting of parameters related to image acquisition and analysis, and key outcome assessments, particularly with respect to ex vivo analysis of rodent specimens. Thus the guidelines herein provide recommendations regarding (1) standardized terminology and units, (2) information to be included in describing the methods for a given experiment, and (3) a minimal set of outcome variables that should be reported. Whereas the specific research objective will determine the experimental design, these guidelines are intended to ensure accurate and consistent reporting of microCT-derived bone morphometry and density measurements. In particular, the methods section for papers that present microCT-based outcomes must include details of the following scan aspects: (1) image acquisition, including the scanning medium, X-ray tube potential, and voxel size, as well as clear descriptions of the size and location of the volume of interest and the method used to delineate trabecular and cortical bone regions, and (2) image processing, including the algorithms used for image filtration and the approach used for image segmentation. Morphometric analyses should be based on 3D algorithms that do not rely on assumptions about the underlying structure whenever possible. When reporting microCT results, the minimal set of variables that should be used to describe trabecular bone morphometry includes bone volume fraction and trabecular number, thickness, and separation. The minimal set of variables that should be used to describe cortical bone morphometry includes total cross-sectional area, cortical bone area, cortical bone area fraction, and cortical thickness. Other variables also may be appropriate depending on the research question and technical quality of the scan. Standard nomenclature, outlined in this article, should be followed for reporting of results.

Reproducibility of bone micro-architecture measurements in rodents by in vivo micro-computed tomography is maximized with three-dimensional image registration

In vivo micro-computed tomography (microCT) is a new method to monitor longitudinal changes of bone micro-architecture. Common animal models of bone diseases are mice and rats, and it is important to know the reproducibility of the bone measurements in order to correctly interpret results. When performing baseline and follow-up acquisitions, variation in the scan region will influence the parameters, and it has yet to be investigated if three-dimensional (3D) registration can improve the reproducibility. Two typical breeds of mice and one typical breed of rats were scanned four times each using microCT and standard bone morphological and density measurements were calculated. Image registration was used to find the overlapping regions between the scans within each series of animal data and only overlapping regions were analyzed for the bone parameters. Reproducibility was determined for each animal both pre- and post-registration. For the rats, results included a bone volume ratio (BV/TV) precision error of 5.46%, cortical thickness (Ct.Th) error of 1.97%, and tissue mineral density (TMD) of 2.00%. For the BL6 mice, precision errors were 3.00% (BV/TV), 0.95% (Ct.Th), and 0.94% (TMD), and for the C3H mice 2.68% (BV/TV), 1.52% (Ct.Th), and 1.72% (TMD). After image registration there was a significant improvement in reproducibility in most parameters for the rats. In general, metric parameters such as bone volume ratio had better reproducibility than the non-metric parameters connectivity density and structure model index. With 3D registration, reproducibility improved the results obtained by the experienced operators in this study. Registration may serve to equalize reproducibility of operators with different skill levels and across laboratories. It also improves efficiency by reducing the amount of hand-contouring required. This reproducibility data will be important for the interpretation of current and future longitudinal microCT studies.

Revealing the interplay of bone and cartilage in osteoarthritis through multimodal imaging of murine joints

Osteoarthritis (OA) affects both cartilage and bone tissues, and the subsequent breakdown of the two tissues appears to be interrelated. The interest in the role of subchondral bone changes with OA is growing, and one suggestion is that a simple inverse correlation exists between the cartilage loss and increased bone mineral density. In this work the STR/ort mouse is used as a model for human OA, in order to investigate disease progression. The aim of the work is to elucidate the tempero-spatial relationships between bone and cartilage architecture and determine whether a simple inverse correlation is satisfactory. We employ 3D whole joint quantitative imaging techniques for assessment of subchondral bone and articular cartilage. The knee joints of mice aged 3, 4, 7 and 10 months are scanned with muCT and then the tibial plateaus are scanned with CLSM. The results show that depending on site (medial and lateral), compartment (epiphyseal, metaphyseal, cortical), and age (3, 4, 7, 10 months), the subchondral bone undergoes changes that lead to an altered architecture. This is primarily seen as densification of the cortex and epiphysis in the STR/ort mice, with a significant change occurring between 7 and 10 months, while the medial cartilage thickness is significantly reduced after 7 months. Using a novel multimodal imaging approach, morphometric changes in the murine osteoarthritic knee joint are elucidated. It is seen that a complex interplay of events - both spatially and temporally - is involved in OA onset and progression. The initial measured differences between the two strains suggest a possible morphological phenotype involved in OA resistance/vulnerability. Temporally the changes have a strong strain:age dependence, although no separate timeline of events between the two tissues could be discerned. Spatially, the changes to medial and lateral morphometry across the cartilage and bone, indicate a relationship to altered joint mechanics.

Phenotypic integration among trabecular and cortical bone traits establishes mechanical functionality of inbred mouse vertebrae

Conventional approaches to identifying quantitative trait loci (QTLs) regulating bone mass and fragility are limited because they examine cortical and trabecular traits independently. Prior work examining long bones from young adult mice and humans indicated that skeletal traits are functionally related and that compensatory interactions among morphological and compositional traits are critical for establishing mechanical function. However, it is not known whether trait covariation (i.e., phenotypic integration) also is important for establishing mechanical function in more complex, corticocancellous structures. Covariation among trabecular, cortical, and compositional bone traits was examined in the context of mechanical functionality for L(4) vertebral bodies across a panel of 16-wk-old female AXB/BXA recombinant inbred (RI) mouse strains. The unique pattern of randomization of the A/J and C57BL/6J (B6) genome among the RI panel provides a powerful tool that can be used to measure the tendency for different traits to covary and to study the biology of complex traits. We tested the hypothesis that genetic variants affecting vertebral size and mass are buffered by changes in the relative amounts of cortical and trabecular bone and overall mineralization. Despite inheriting random sets of A/J and B6 genomes, the RI strains inherited nonrandom sets of cortical and trabecular bone traits. Path analysis, which is a multivariate analysis that shows how multiple traits covary simultaneously when confounding variables like body size are taken into consideration, showed that RI strains that tended to have smaller vertebrae relative to body size achieved mechanical functionality by increasing mineralization and the relative amounts of cortical and trabecular bone. The interdependence among corticocancellous traits in the vertebral body indicated that variation in trabecular bone traits among inbred mouse strains, which is often thought to arise from genetic factors, is also determined in part by the adaptive response to variation in traits describing the cortical shell. The covariation among corticocancellous traits has important implications for genetic analyses and for interpreting the response of bone to genetic and environmental perturbations.