Compartmental bone morphometry in the mouse femur: reproducibility and resolution dependence of microtomographic measurements

Loss of the glycosyltransferase Galnt11 affects vitamin D homeostasis and bone composition

O-glycosylation is a conserved posttranslational modification that impacts many aspects of organismal viability and function. Recent studies examining the glycosyltransferase Galnt11 demonstrated that it glycosylates the endocytic receptor megalin in the kidneys, enabling proper binding and reabsorption of ligands, including vitamin D-binding protein (DBP). Galnt11-deficient mice were unable to properly reabsorb DBP from the urine. Vitamin D plays an essential role in mineral homeostasis and its deficiency is associated with bone diseases such as rickets, osteomalacia, and osteoporosis. We therefore set out to examine the effects of the loss of Galnt11 on vitamin D homeostasis and bone composition. We found significantly decreased levels of serum 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D, consistent with decreased reabsorption of DBP. This was accompanied by a significant reduction in blood calcium levels and a physiologic increase in parathyroid hormone (PTH) in Galnt11-deficient mice. Bones in Galnt11-deficient mice were smaller and displayed a decrease in cortical bone accompanied by an increase in trabecular bone and an increase in a marker of bone formation, consistent with PTH-mediated effects on bone. These results support a unified model for the role of Galnt11 in bone and mineral homeostasis, wherein loss of Galnt11 leads to decreased reabsorption of DBP by megalin, resulting in a cascade of disrupted mineral and bone homeostasis including decreased circulating vitamin D and calcium levels, a physiological increase in PTH, an overall loss of cortical bone, and an increase in trabecular bone. Our study elucidates how defects in O-glycosylation can influence vitamin D and mineral homeostasis and the integrity of the skeletal system.

A microCT imaging protocol for reproducible and efficient quantitative morphometric analysis (QMA) of joint structures of the in situ mouse tibio-femoral joint

Micro-computed tomography (microCT) offers a three-dimensional (3D), high-resolution technique for the visualisation and analysis of bone microstructure. Using contrast-enhanced microCT, this capability has been expanded in recent studies to include cartilage morphometry and whole joint measures, known together as quantitative morphometric analysis (QMA). However, one of the main challenges in quantitative analysis of joint images is sensitivity to joint pose and alignment, which may influence measures related to both joint space and joint biomechanics. Thus, this study proposes a novel microCT imaging protocol for reproducible and efficient QMA of in situ mouse tibio-femoral joint. This work consists of two parts: an in situ diffusion kinetics study for a known cationic iodinated contrast agent (CA4+) for QMA of the cartilage, and a joint positioning and image processing workflow for whole joint QMA. In the diffusion kinetics study, 8 mice were injected at both of their tibio-femoral joints with distinct CA4+ concentrations and diffusion times. The mice were scanned at different time points after injection, and evaluated using attenuation and cartilage QMA measures. Results show that cartilage segmentation and QMA could be performed for CA4+ solution at a concentration of 48 mg/ml, and that reliable measurement and quantification of cartilage were achieved after 5 min of diffusion following contrast agent injection. We established the joint positioning and image processing workflow by developing a novel positioning device to control joint pose during scanning, and a spherical harmonics-based image processing workflow to ensure consistent alignment during image processing. Both legs of seven mice were scanned 10 times, 5 prior to receiving CA4+ and 5 after, and evaluated using whole joint QMA parameters. Joint QMA evaluation of the workflow showed excellent reproducibility; intraclass correlation coefficients ranged from 0.794 to 0.930, confirming that the imaging protocol enables reproducible and efficient QMA of joint structures in preclinical models, and that contrast agent injection did not cause significant alteration to the measured parameters.

Nuclear factor erythroid 2-related factor 2 (Nrf2) deficiency causes age-dependent progression of female osteoporosis

BACKGROUND

Nuclear factor erythroid 2-related factor 2 (Nrf2) is a crucial transcription factor for cellular redox homeostasis. The association of Nrf2 with elderly female osteoporotic has yet to be fully described. The aim was to elucidate a potential age-dependent Nrf2 contribution to female osteoporosis in mice.

METHODS

Eighteen female wild type (WT) and 16 Nrf2-knockout (KO) mice were sacrificed at different ages (12 weeks = young mature adult and 90 weeks = old) to analyze their femurs. The morphological properties (trabecular and cortical) were evaluated by micro-computed tomography (μCT) and compared to gold standard histochemistry analysis. The quasi-static compression tests were performed to calculate the mechanical properties of bones. Additionally, the population of bone resorbing cells and aromatase expression by osteocytes was immunohistochemically evaluated and empty osteocyte lacunae was counted in cortical bone.

RESULTS

Old Nrf2-KO mice revealed a significantly reduced trabecular bone mineral density (BMD), cortical thickness, cortical area, and bone fraction compared to old WT mice, regardless of no significant difference in skeletally mature young adult mice between WT and KO. Specifically, while all old WT mice showed thin metaphyseal trabeculae, trabecular bone was completely absent in 60% of old KO mice. Additionally, old KO mice showed significantly more osteoclast-like cells and fewer aromatase-positive osteocytes than WT mice, whereas the occurrence of empty osteocyte lacunae did not differ between both groups. Nrf2-KO mice further showed an age-dependently reduced fracture resilience compared to age-matched WT mice.

CONCLUSION

Our results suggest that chronic Nrf2 loss can lead to age-dependent progression of female osteoporosis.

Reproducibility of Densitometric and Biomechanical Assessment of the Mouse Tibia From In Vivo Micro-CT Images

Interventions for bone diseases (e.g. osteoporosis) require testing in animal models before clinical translation and the mouse tibia is among the most common tested anatomical sites. In vivo micro-Computed Tomography (microCT) based measurements of the geometrical and densitometric properties are non-invasive and therefore constitute an important tool in preclinical studies. Moreover, validated micro-Finite Element (microFE) models can be used for predicting the bone mechanical properties non-invasively. However, considering that the image processing pipeline requires operator-dependant steps, the reproducibility of these measurements has to be assessed. The aim of this study was to evaluate the intra- and inter-operator reproducibility of several bone parameters measured from microCT images. Ten in vivo microCT images of the right tibia of five mice (at 18 and 22 weeks of age) were processed. One experienced operator (intra-operator analysis) and three different operators (inter-operator) aligned each image to a reference through a rigid registration and selected a volume of interest below the growth plate. From each image the following parameters were measured: total bone mineral content (BMC) and density (BMD), BMC in 40 subregions (ten longitudinal sections, four quadrants), microFE-based stiffness and failure load. Intra-operator reproducibility was acceptable for all parameters (precision error, PE < 3.71%), with lowest reproducibility for stiffness (3.06% at week 18, 3.71% at week 22). The inter-operator reproducibility was slightly lower (PE < 4.25%), although still acceptable for assessing the properties of most interventions. The lowest reproducibility was found for BMC in the lateral sector at the midshaft (PE = 4.25%). Densitometric parameters were more reproducible than most standard morphometric parameters calculated in the proximal trabecular bone. In conclusion, microCT and microFE models provide reproducible measurements for non-invasive assessment of the mouse tibia properties.

SPHARM-PDM based image preprocessing pipeline for quantitative morphometric analysis (QMA) for in situ joint assessment in rabbit and rat models

The accessibility of quantitative measurements of joint morphometry depends on appropriate tibial alignment and volume of interest (VOI) selection of joint compartments; often a challenging and time-consuming manual task. In this work, we developed a novel automatic, efficient, and model-invariant image preprocessing pipeline that allows for highly reproducible 3D quantitative morphometric analysis (QMA) of the joint. The pipeline addresses the problem by deploying two modules: an alignment module and a subdivision module. Alignment is achieved by representing the tibia in its basic form using lower degree spherical harmonic basis functions and aligning using principal component analysis. The second module subdivides the joint into lateral and medial VOIs via a watershedding approach based on persistence homology. Multiple repeated micro-computed tomography scans of small (rat) and medium (rabbit) animal knees were processed using the pipeline to demonstrate model invariance. Existing QMA was performed to evaluate the pipeline’s ability to generate reproducible measurements. Intraclass correlation coefficient and mean-normalised root-mean-squared error of more than 0.75 and lower than 9.5%, respectively, were achieved for joint centre of mass, joint contact area under virtual loading, joint space width, and joint space volume. Processing time and technical requirements were reduced compared to manual processing in previous studies.

Quantitative measures of bone shape, cartilage morphometry and joint alignment are associated with disease in an ACLT and MMx rat model of osteoarthritis

Multi-scale, subject-specific quantitative methods to characterize and monitor osteoarthritis in animal models and therapeutic treatments could help reveal causal relationships in disease development and distinguish treatment strategies. In this work, we demonstrate a reproducible and sensitive quantitative image analysis to characterize bone, cartilage and joint measures describing a rat model of post-traumatic osteoarthritis. Eleven 3-month-old male Wistar rats underwent medial anterior cruciate ligament (ACL) transection and medial meniscectomy on the right knee to destabilise the right tibiofemoral joint. They were sacrificed 6 weeks post-surgery and a silicon-based micro-bead contrast agent was injected in the joint space, before scanning with micro-computed tomography (microCT). Subsequently, 3D quantitative morphometric analysis (QMA), previously developed for rabbit joints, was performed. This included cartilage, subchondral cortical and epiphyseal bone measures, as well as novel tibiofemoral joint metrics. Semi-quantitative evaluation was performed on matching two-dimensional (2D) histology and microCT images. Reproducibility of the QMA was tested on eleven age-matched additional joints. The results indicate the QMA method is accurate and reproducible and that microCT-derived cartilage measurements are valid for the analysis of rat joints. The pathologic changes caused by transection of the ACL and medial meniscectomy were reflected in measurements of bone shape, cartilage morphology, and joint alignment. Furthermore, we were able to identify model-specific predictive parameters based on morphometric parameters measured with the QMA.

OPG-Fc treatment partially rescues low bone mass phenotype in mature Bgn/Fmod deficient mice but is deleterious to the young mouse skeleton

Biglycan (Bgn) and Fibromodulin (Fmod) are small leucine rich proteoglycans (SLRPs) which are abundant in the extra-cellular matrix (ECM) of mineralized tissues. We have previously generated a Bgn/Fmod double knock-out (DKO) mouse model and found it has a 3-fold increase in osteoclastogenesis compared with Wild type (WT) controls, resulting in a markedly low bone mass (LBM) phenotype. To try and rescue/repair the LBM phenotype of Bgn/Fmod DKO mice by suppressing osteoclast formation and activity, 3- and 26-week-old Bgn/Fmod DKO mice and age/gender matched WT controls were treated with OPG-Fc for 6 weeks after which bone parameters were evaluated using DEXA, micro-computed tomography (μCT) and serum biomarkers analyses. In the appendicular skeleton, OPG-Fc treatment improved some morphometric and geometric parameters in both the trabecular and cortical compartments in Bgn/Fmod DKO female and male mice, especially in the repair module. For many of the skeletal parameters analyzed, the Bgn/Fmod DKO mice were more responsive to the treatment than their WT controls. In addition, we found that OPG-Fc treatment was not able to prevent or ameliorate the formation of ectopic ossification, which are common lesions seen in aged joints and are one of the phenotypical hallmarks of our Bgn/Fmod DKO model. Analysis of skull bones, specifically the occipital bone, showed the treatment recovered some parameters of LBM phenotype in the craniofacial skeleton, more so in the younger rescue module. Using OPG-Fc as treatment alleviated, yet did not completely restore, the severe osteopenia and mineralized tissue structural abnormalities that Bgn/Fmod DKO mice suffer from.

Effects of Early Life Stress on Bone Homeostasis in Mice and Humans

Bone pathology is frequent in stressed individuals. A comprehensive examination of mechanisms linking life stress, depression and disturbed bone homeostasis is missing. In this translational study, mice exposed to early life stress (MSUS) were examined for bone microarchitecture (μCT), metabolism (qPCR/ELISA), and neuronal stress mediator expression (qPCR) and compared with a sample of depressive patients with or without early life stress by analyzing bone mineral density (BMD) (DXA) and metabolic changes in serum (osteocalcin, PINP, CTX-I). MSUS mice showed a significant decrease in NGF, NPYR1, VIPR1 and TACR1 expression, higher innervation density in bone, and increased serum levels of CTX-I, suggesting a milieu in favor of catabolic bone turnover. MSUS mice had a significantly lower body weight compared to control mice, and this caused minor effects on bone microarchitecture. Depressive patients with experiences of childhood neglect also showed a catabolic pattern. A significant reduction in BMD was observed in depressive patients with childhood abuse and stressful life events during childhood. Therefore, future studies on prevention and treatment strategies for both mental and bone disease should consider early life stress as a risk factor for bone pathologies.

Collagen VIα2 chain deficiency causes trabecular bone loss by potentially promoting osteoclast differentiation through enhanced TNFα signaling

Type VI collagen is well known for its role in muscular disorders, however its function in bone is still not well understood. To examine its role in bone we analyzed femoral and vertebral bone mass by micro-computed tomography analysis, which showed lower bone volume/total volume and trabecular number in Col6α2-KO mice compared with WT. Dynamic histomorphometry showed no differences in trabecular bone formation between WT and Col6α2-KO mice based on the mineral appositional rate, bone formation rate, and mineralizing perimeter. Femoral sections were assessed for the abundance of Tartrate Resistant Acid Phosphatase-positive osteoclasts, which revealed that mutant mice had more osteoclasts compared with WT mice, indicating that the primary effect of Col6a2 deficiency is on osteoclastogenesis. When bone marrow stromal cells (BMSCs) from WT and Col6α2-KO mice were treated with rmTNFα protein, the Col6α2-KO cells expressed higher levels of TNFα mRNA compared with WT cells. This was accompanied by higher levels of p-p65, a down-stream target of TNFα, suggesting that BMSCs from Col6α2-KO mice are highly sensitive to TNFα signaling. Taken together, our data imply that Col6a2 deficiency causes trabecular bone loss by enhancing osteoclast differentiation through enhanced TNFα signaling.

Reproducibility and Radiation Effect of High-Resolution In Vivo Micro Computed Tomography Imaging of the Mouse Lumbar Vertebra and Long Bone

A moderate radiation dose, in vivo µCT scanning protocol was developed and validated for long-term monitoring of multiple skeletal sites (femur, tibia, vertebra) in mice. A customized, 3D printed mouse holder was designed and utilized to minimize error associated with animal repositioning, resulting in good to excellent reproducibility in most cortical and trabecular bone microarchitecture and density parameters except for connectivity density. Repeated in vivo µCT scans of mice were performed at the right distal femur and the 4th lumbar vertebra every 3 weeks until euthanized at 9 weeks after the baseline scan. Comparing to the non-radiated counterparts, no radiation effect was found on trabecular bone volume fraction, osteoblast and osteoblast number/surface, or bone formation rate at any skeletal site. However, trabecular number, thickness, and separation, and structure model index were sensitive to ionizing radiation associated with the µCT scans, resulting in subtle but significant changes over multiple scans. Although the extent of radiation damage on most trabecular bone microarchitecture measures are comparable or far less than the age-related changes during the monitoring period, additional considerations need to be taken to minimize the confounding radiation factors when designing experiments using in vivo µCT imaging for long-term monitoring of mouse bone.

A new method to monitor bone geometry changes at different spatial scales in the longitudinal in vivo μCT studies of mice bones

Longitudinal studies of bone adaptation in mice using in vivo micro-computed tomography (μCT) have been commonly used for pre-clinical evaluation of physical and pharmacological interventions. The main advantage of this approach is to use each mouse as its own control, reducing considerably the sample size required by statistical power analysis. To date, multi-scale estimation of bone adaptations become essential since the bone activity that takes place at different scales may be associated with different bone mechanisms. Measures of bone adaptations at different time scales have been attempted in a previous study. This paper extends quantification of bone activity at different spatial scales with a proposition of a novel framework. The method involves applying level-set method (LSM) to track the geometric changes from the longitudinal in vivo μCT scans of mice tibia. Bone low- and high-spatial frequency patterns are then estimated using multi-resolution analysis. The accuracy of the framework is quantified by applying it to two times separated scanned images with synthetically manipulated global and/or local activity. The Root Mean Square Deviation (RMSD) was approximately 1.5 voxels or 0.7 voxels for the global low-spatial frequency or local high-spatial frequency changes, respectively. The framework is further applied to the study of bone changes in longitudinal datasets of wild-type mice tibiae over time and space. The results demonstrate the ability for the spatio-temporal quantification and visualisation of bone activity at different spatial scales in longitudinal studies thus providing further insight into bone adaptation mechanisms.

Effect of parathyroid hormone on the structural, densitometric and failure behaviors of mouse tibia in the spatiotemporal space

Parathyroid hormone (PTH) is an anabolic bone drug approved by the US Food and Drug Administration (FDA) to treat osteoporosis. However, previous studies using cross-sectional designs have reported variable and sometimes contradictory results. The aim of the present study was to quantify the localized effect of PTH on the structural and densitometric behaviors of mouse tibia and their links with the global mechanical behavior of bone using a novel spatiotemporal image analysis approach and a finite element analysis technique. Twelve female C57BL/6J mice were divided into two groups: the control and PTH treated groups. The entire right tibiae were imaged using an in vivo micro-computed tomography (μCT) system eight consecutive times. Next, the in vivo longitudinal tibial μCT images were rigidly registered and divided into 10 compartments across the entire tibial space. The bone volume (BV), bone mineral content (BMC), bone tissue mineral density (TMD), and tibial endosteal and periosteal areas (TEA and TPA) were quantified in each compartment. Additionally, finite element models of all the tibiae were generated to analyze the failure behavior of the tibia. It was found that both the BMC and BV started to increase in the proximal tibial region, and then the increases extended to the entire tibial region after two weeks of treatment (p < 0.05). PTH intervention significantly reduced the TEA in most tibial compartments after two weeks of treatment, and the TPA increased in most tibial regions after four weeks of treatment (p < 0.05). Tibial failure loads significantly increased after three weeks of PTH treatment (p < 0.01). The present study provided the first evidence of the localized effect of PTH on bone structural and densitometric properties, as well as their links with the global mechanical behaviors of bone, which are important pieces of information for unveiling the mechanism of PTH intervention.

A novel approach for computing 3D mice distal femur properties using high-resolution micro-computed tomography scanning

One of the most-scanned joints in preclinical animal models dealing with musculoskeletal pathologies is the mouse knee. While three-dimensional (3D) characterization of bone tissue porosity have previously been performed on cortical bone, it has not yet been comprehensively performed for the subchondral bone (SB) and the calcified cartilage (CC), which compose the subchondral mineralized zone (SMZ). Thus, it remains challenging to assess changes that occur in the SMZ of the mouse knee during pathologies such as osteoarthritis. One of the keys to addressing this challenge is to segment each layer to measure their morphologies, material properties, and porosity. Our study presents a novel approach for computing Tissue Mineral Density, 3D porosity, and the thickness of SB and CC in a mouse distal femur using High-Resolution Micro-Computed Tomography (HR-μCT). We have segmented the Vascular Porosity network, the osteocytes' lacunae of the SB, and the chondrocytes of the CC by using multi-thresholding and the percentage of chondrocytes porosity. Our results show a low intra- and inter-observer coefficient of variability. Regarding porosity and geometrical properties of both CC and SB, our results are within the range of the literature. Our approach opens new avenues for assessing porosity and vascular changes in the distal femur of preclinical animal models dealing with musculoskeletal pathologies such as osteoarthritis.

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.

Use of Micro-Computed Tomography for Bone Evaluation in Dentistry

Micro computed tomography (µCT) follows the same principle of computed-tomography used for patients, however providing higher-resolution. Using a non-destructive approach, samples can be scanned, and each section obtained is used to build a volume using tridimensional reconstruction. For bone analysis, it is possible to obtain information about the tissue's microarchitecture and composition. According to the characteristics of the bone sample (e.g. human or animal origin, long or irregular shape, epiphysis or diaphysis region) the pre-scanning parameters must be defined. The resolution (i.e. voxel size) should be chosen taking into account the features that will be evaluated, and the necessity to identify inner structures (e.g. bone channels and osteocyte lacunae). The region of interest should be delimited, and the threshold that defines the bone tissue set in order to proceed with binarization to separate the voxels representing bone from the other structures (channels, resorption areas, and medullary space). Cancellous bone is evaluated by means of the trabeculae characteristics and their connectivity. The cortex is evaluated in relation to the thickness and porosity. Bone mineral density can also be measured, by the amount of hydroxyapatite. Other parameters such as structure-model-index, anisotropy, and fractal dimension can be assessed. In conclusion, intrinsic and extrinsic determinants of bone quality can be assessed by µCT. In dentistry, this method can be used for evaluating bone loss, alterations in bone metabolism, or the effects of using drugs that impair bone remodeling, and also to assess the success rate of bone repair or surgical procedures.

Eight Days of Earth Reambulation Worsen Bone Loss Induced by 1-Month Spaceflight in the Major Weight-Bearing Ankle Bones of Mature Mice

Spaceflight induces bone alterations with site-specific rates of bone loss according to the weight-bearing function of the bone. For the first time, this study aimed to characterize bone microarchitecture and density alterations of three ankle bones (calcaneus, navicular, and talus) of mice after spaceflight and to evaluate the impact of 8 days of Earth reambulation. Ten C57BL/6N male 4-month-old mice flew on the Bion-M1 biosatellite for 1 month; half were euthanized within 24-h of return and half after 8-days recovery on Earth. Bone microarchitecture and quality was assessed by microtomography (μCT). Whole calcaneus bone volume fraction decreased in Flight group (-6.4%, p < 0.05), and worsened in the Recovery group (-11.08%, p < 0.01), when compared to Control group. Navicular and talus trabecular bone volume fraction showed trends toward decrease in Flight and differences reached statistical significance in Recovery group (-8.16%; -8.87%, respectively; p < 0.05) when compared to Control group. At calcaneus, cortical thickness decreased in Recovery vs. Control groups (-11.69%; p < 0.01). Bone surface area, reflecting periosteal bone erosion, significantly increased in all bone sites analyzed. Qualitative analyses of 3-D bone reconstruction revealed local sites of cortical thinning and bone erosion, predominantly at articulations, muscle insertions, and ground contact bone sites. Overall, spaceflight-induced bone loss in ankle bones was site and compartment specific whilst the tissue mineral density of the remaining bone was preserved. Eight days after landing, bone status worsened as compared to immediate return.

Small leucine rich proteoglycans, a novel link to osteoclastogenesis

Biglycan (Bgn) and Fibromodulin (Fmod) are subtypes of the small leucine-rich family of proteoglycans (SLRP). In this study we examined the skeletal phenotype of BgnFmod double knockout (BgnFmod KO) mice and found they were smaller in size and have markedly reduced bone mass compared to WT. The low bone mass (LBM) phenotype is the result of both the osteoblasts and osteoclasts from BgnFmod KO mice having higher differentiation potential and being more active compared to WT mice. Using multiple approaches, we showed that both Bgn and Fmod directly bind TNFα as well as RANKL in a dose dependent manner and that despite expressing higher levels of both TNFα and RANKL, BgnFmod KO derived osteoblasts cannot retain these cytokines in the vicinity of the cells, which leads to elevated TNFα and RANKL signaling and enhanced osteoclastogenesis. Furthermore, adding either Bgn or Fmod to osteoclast precursor cultures significantly attenuated the cells ability to form TRAP positive, multinucleated giant cells. In summary, our data indicates that Bgn and Fmod expressed by the bone forming cells, are novel coupling ECM components that control bone mass through sequestration of TNFα and/or RANKL, thereby adjusting their bioavailability in order to regulate osteoclastogenesis.

Automated Quantitative Bone Analysis in In Vivo X-ray Micro-Computed Tomography

Measurement and analysis of bone morphometry in 3D micro-computed tomography volumes using automated image processing and analysis improve the accuracy, consistency, reproducibility, and speed of preclinical osteological research studies. Automating segmentation and separation of individual bones in 3D micro-computed tomography volumes of murine models presents significant challenges considering partial volume effects and joints with thin spacing, i.e., 50 to [Formula: see text]. In this paper, novel hybrid splitting filters are presented to overcome the challenge of automated bone separation. This is achieved by enhancing joint contrast using rotationally invariant second-derivative operators. These filters generate split components that seed marker-controlled watershed segmentation. In addition, these filters can be used to separate metaphysis and epiphysis in long bones, e.g., femur, and remove the metaphyseal growth plate from the detected bone mask in morphometric measurements. Moreover, for slice-by-slice stereological measurements of long bones, particularly curved bones, such as tibia, the accuracy of the analysis can be improved if the planar measurements are guided to follow the longitudinal direction of the bone. In this paper, an approach is presented for characterizing the bone medial axis using morphological thinning and centerline operations. Building upon the medial axis, a novel framework is presented to automatically guide stereological measurements of long bones and enhance measurement accuracy and consistency. These image processing and analysis approaches are combined in an automated streamlined software workflow and applied to a range of in vivo micro-computed tomography studies for validation.

Time course of peri-implant bone regeneration around loaded and unloaded implants in a rat model

The time-course of cancellous bone regeneration surrounding mechanically loaded implants affects implant fixation, and is relevant to determining optimal rehabilitation protocols following orthopaedic surgeries. We investigated the influence of controlled mechanical loading of titanium-coated polyether-ether ketone (PEEK) implants on osseointegration using time-lapsed, non-invasive, in vivo micro-computed tomography (micro-CT) scans. Implants were inserted into proximal tibial metaphyses of both limbs of eight female Sprague-Dawley rats. External cyclic loading (60 or 100 μm displacement, 1 Hz, 60 s) was applied every other day for 14 days to one implant in each rat, while implants in contralateral limbs served as the unloaded controls. Hind limbs were imaged with high-resolution micro-CT (12.5 μm voxel size) at 2, 5, 9, and 12 days post-surgery. Trabecular changes over time were detected by 3D image registration allowing for measurements of bone-formation rate (BFR) and bone-resorption rate (BRR). At day 9, mean %BV/TV for loaded and unloaded limbs were 35.5 ± 10.0% and 37.2 ± 10.0%, respectively, and demonstrated significant increases in bone volume compared to day 2. BRR increased significantly after day 9. No significant differences between bone volumes, BFR, and BRR were detected due to implant loading. Although not reaching significance (p = 0.16), an average 119% increase in pull-out strength was measured in the loaded implants. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:997-1006, 2017.

Exercise improves femoral whole-bone and tissue-level biomechanical properties in hyperphagic OLETF rats

An often-overlooked comorbidity of type 2 diabetes (T2D) is increased fracture risk. Since traditional T2D therapies focus solely on glucose homeostasis, there is an increased need for effective treatment strategies that improve both metabolic and bone health. The current study evaluated if previously reported protective effects of exercise (EX) on cortical bone geometry and biomechanical properties are due to the prevention of obesity compared with obese Otsuka Long-Evans Tokushima Fatty (OLETF) rats or the unique skeletal effects of exercise. Four-week-old male OLETF rats were randomly assigned to 1 of 3 groups, each fed a standard rodent chow diet until 40 weeks of age: ad libitum-fed sedentary (O-SED), ad libitum-fed EX (O-EX), or a control group body-weight-matched to the O-EX group by caloric restriction (O-CR). Ad libitum-fed sedentary Long-Evans Tokushima Otsuka (L-SED) rats were used as a lean control. EX or CR prevented the excess body mass accumulation and elevated HbA1c observed in O-SED. Total-body bone mineral density was greater in O-EX than O-CR, but similar to L-SED and O-SED. O-EX and O-CR had lower osteocalcin and TRAP5b than O-SED and L-SED. EX or CR prevented the ∼3-fold increase in CTx in O-SED versus L-SED. EX increased femoral cortical mass accumulation and expansion at the mid-diaphysis compared with O-CR. EX or CR significantly increased tissue-level stiffness and strength compared with O-SED and L-SED, but O-EX had greater whole-bone stiffness than all other groups. In summary, EX has distinct advantages over CR for improving bone biomechanical properties in hyperphagic OLETF rats.

Lysine-specific demethylase 1 inhibitor rescues the osteogenic ability of mesenchymal stem cells under osteoporotic conditions by modulating H3K4 methylation

Bone tissue engineering may be hindered by underlying osteoporosis because of a decreased osteogenic ability of autologous seed cells and an unfavorably changed microenvironment in these patients. Epigenetic regulation plays an important role in the developmental origins of osteoporosis; however, few studies have investigated the potential of epigenetic therapy to improve or rescue the osteogenic ability of bone marrow mesenchymal stem cells (BMMSCs) under osteoporotic conditions. Here, we investigated pargyline, an inhibitor of lysine-specific demethylase 1 (LSD1), which mainly catalyzes the demethylation of the di- and mono-methylation of H3K4. We demonstrated that 1.5 mmol·L⁻¹ pargyline was the optimal concentration for the osteogenic differentiation of human BMMSCs. Pargyline rescued the osteogenic differentiation ability of mouse BMMSCs under osteoporotic conditions by enhancing the dimethylation level of H3K4 at the promoter regions of osteogenesis-related genes. Moreover, pargyline partially rescued or prevented the osteoporotic conditions in aged or ovariectomized mouse models, respectively. By introducing the concept of epigenetic therapy into the field of osteoporosis, this study demonstrated that LSD1 inhibitors could improve the clinical practice of MSC-based bone tissue engineering and proposes their novel use to treat osteoporosis.

Effect of combined treatment with zoledronic acid and parathyroid hormone on mouse bone callus structure and composition

In recent years, great interest in combined treatment of parathyroid hormone (PTH) with anti-resorptive therapy has emerged. PTH has been suggested to aid bridging of atrophic fractures and improve strength in closed fracture models. Bisphosphonate treatments typically result in a larger woven bone callus that is slower to remodel. The combination of both drugs has been demonstrated to be effective for the treatment of osteoporotic bone loss in many preclinical studies. However, the effect of combined treatment on fracture repair is still largely unexplored. In this study, we aimed to compare these drugs as single-agent and in combination in a murine closed fracture model. We wanted to assess potential differences in material properties, morphometry and in the development of the lacuno-canalicular network. A total of 40 female, 11-week-old wild type mice underwent a closed fracture on the midshaft of the tibia and were assigned to four groups (n=8-10 per group). Beginning on post-operative day 8, animals received different subcutaneous injections. Group 1 received a single injection of saline solution and Group 2 of zoledronic acid (ZA). Group 3 received daily dosing of PTH. Group 4 received a dual treatment, starting with a single dose of ZA followed by daily injection of PTH. Three weeks after fracture, all animals were euthanized and tibiae were assessed using micro-computed tomography (micro-CT), high-resolution micro-CT (HR micro-CT), Raman spectroscopy, quantitative histomorphometry, and deconvolution microscopy (DV microscopy). Combined treatment showed a significant increase of 41% in bone volume fraction and a significant decrease of 61% in the standard deviation of the trabecular spacing compared to vehicle, both known to be strong predictors of callus strength. An analysis via HR micro-CT showed similar results on all groups for lacunar numerical density, whereas mean lacuna volume was found to be higher compared to vehicle in treated groups, but only PTH mono-treatment showed a significant increase compared to vehicle (+45%). Raman spectroscopy did not reveal detectable changes in material properties of the bone calluses. Sclerostin staining, tartrate resistant acid phosphatase (TRAP) staining and canalicular analysis with DV microscopy on a subset of samples did not display distinctive difference in any of the treatments. We therefore consider PTH+ZA treatment beneficial for bone healing. No clear negative effect on bone quality was detected during this study.

Voxel size dependency, reproducibility and sensitivity of an in vivo bone loading estimation algorithm

A bone loading estimation algorithm was previously developed that provides in vivo loading conditions required for in vivo bone remodelling simulations. The algorithm derives a bone's loading history from its microstructure as assessed by high-resolution (HR) computed tomography (CT). This reverse engineering approach showed accurate and realistic results based on micro-CT and HR-peripheral quantitative CT images. However, its voxel size dependency, reproducibility and sensitivity still need to be investigated, which is the purpose of this study. Voxel size dependency was tested on cadaveric distal radii with micro-CT images scanned at 25 µm and downscaled to 50, 61, 75, 82, 100, 125 and 150 µm. Reproducibility was calculated with repeated in vitro as well as in vivo HR-pQCT measurements at 82 µm. Sensitivity was defined using HR-pQCT images from women with fracture versus non-fracture, and low versus high bone volume fraction, expecting similar and different loading histories, respectively. Our results indicate that the algorithm is voxel size independent within an average (maximum) error of 8.2% (32.9%) at 61 µm, but that the dependency increases considerably at voxel sizes bigger than 82 µm. In vitro and in vivo reproducibility are up to 4.5% and 10.2%, respectively, which is comparable to other in vitro studies and slightly higher than in other in vivo studies. Subjects with different bone volume fraction were clearly distinguished but not subjects with and without fracture. This is in agreement with bone adapting to customary loading but not to fall loads. We conclude that the in vivo bone loading estimation algorithm provides reproducible, sensitive and fairly voxel size independent results at up to 82 µm, but that smaller voxel sizes would be advantageous.

Dual-energy X-ray absorptiometry, peripheral quantitative computed tomography, and micro-computed tomography techniques are discordant for bone density and geometry measurements in the guinea pig

This study aims to examine agreement among bone mineral content (BMC) and density (BMD) estimates obtained using dual-energy X-ray absorptiometry (DXA), peripheral quantitative computed tomography (pQCT), and micro-computed tomography (μCT) against high-resolution μCT and bone ash of the guinea pig femur. Middle-aged (n = 40, 86 weeks) male guinea pigs underwent in vivo followed by ex vivo DXA (Hologic QDR 4500A) scanning for intact and excised femur BMC and areal density. To assess bone architecture and strength, excised femurs were scanned on pQCT (Stratec XCT 2000L) as well as on two μCT scanners (LaTheta LCT-200; Skyscan 1174), followed by three-point bending test. Reproducibility was determined using triplicate scans; and agreement assessed using Bland-Altman plots with reference methods being high-resolution μCT (Skyscan) for BMD and bone ashing for BMC. All techniques showed satisfactory ex vivo precision (CV 0.05-4.3 %). However, bias compared to the reference method was highest (207.5 %) in trabecular bone volume fraction (BV/TV) measured by LaTheta, and unacceptable in most total femur and cortical bone measurements. Volumetric BMD (vBMD) and BV/TV derived by LaTheta and pQCT at the distal metaphysis were biased from the Skyscan by an average of 49.3 and 207.5 %, respectively. Variability of vBMD, BV/TV and cross-sectional area at the diaphysis ranged from -5.5 to 30.8 %. LaTheta best quantified total femur BMC with an upper bias of 3.3 %. The observed differences among imaging techniques can be attributable to inherent dissimilarity in construction design, calibration, segmentation and scanning resolution used. These bone imaging tools are precise but are not comparable, at least when assessing guinea pig bones.

Minimizing Interpolation Bias and Precision Error in In Vivo µCT-Based Measurements of Bone Structure and Dynamics

In vivo µCT imaging allows for high-resolution, longitudinal evaluation of bone properties. Based on this technology, several recent studies have developed in vivo dynamic bone histomorphometry techniques that utilize registered µCT images to identify regions of bone formation and resorption, allowing for longitudinal assessment of bone remodeling. However, this analysis requires a direct voxel-by-voxel subtraction between image pairs, necessitating rotation of the images into the same coordinate system, which introduces interpolation errors. We developed a novel image transformation scheme, matched-angle transformation (MAT), whereby the interpolation errors are minimized by equally rotating both the follow-up and baseline images instead of the standard of rotating one image while the other remains fixed. This new method greatly reduced interpolation biases caused by the standard transformation. Additionally, our study evaluated the reproducibility and precision of bone remodeling measurements made via in vivo dynamic bone histomorphometry. Although bone remodeling measurements showed moderate baseline noise, precision was adequate to measure physiologically relevant changes in bone remodeling, and measurements had relatively good reproducibility, with intra-class correlation coefficients of 0.75-0.95. This indicates that, when used in conjunction with MAT, in vivo dynamic histomorphometry provides a reliable assessment of bone remodeling.

In vivo monitoring of bone architecture and remodeling after implant insertion: The different responses of cortical and trabecular bone

The mechanical integrity of the bone-implant system is maintained by the process of bone remodeling. Specifically, the interplay between bone resorption and bone formation is of paramount importance to fully understand the net changes in bone structure occurring in the peri-implant bone, which are eventually responsible for the mechanical stability of the bone-implant system. Using time-lapsed in vivo micro-computed tomography combined with new composite material implants, we were able to characterize the spatio-temporal changes of bone architecture and bone remodeling following implantation in living mice. After insertion, implant stability was attained by a quick and substantial thickening of the cortical shell which counteracted the observed loss of trabecular bone, probably due to the disruption of the trabecular network. Within the trabecular compartment, the rate of bone formation close to the implant was transiently higher than far from the implant mainly due to an increased mineral apposition rate which indicated a higher osteoblastic activity. Conversely, in cortical bone, the higher rate of bone formation close to the implant compared to far away was mostly related to the recruitment of new osteoblasts as indicated by a prevailing mineralizing surface. The behavior of bone resorption also showed dissimilarities between trabecular and cortical bone. In the former, the rate of bone resorption was higher in the peri-implant region and remained elevated during the entire monitoring period. In the latter, bone resorption rate had a bigger value away from the implant and decreased with time. Our approach may help to tune the development of smart implants that can attain a better long-term stability by a local and targeted manipulation of the remodeling process within the cortical and the trabecular compartments and, particularly, in bone of poor health.

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.

Bone adaptation to cyclic loading in murine caudal vertebrae is maintained with age and directly correlated to the local micromechanical environment

The ability of the skeleton to adapt to mechanical stimuli (mechanosensitivity) has most often been investigated at the whole-bone level, but less is known about the local mechanoregulation of bone remodeling at the bone surface, especially in context of the aging skeleton. The aim of this study was to determine the local and global mechanosensitivity of the sixth caudal vertebra during cyclic loading (8 N, three times per week, for six weeks) in mice aged 15, 52, and 82 weeks at the start of loading. Bone adaptation was monitored with in vivo micro-computed tomography. Strain energy density (SED), assumed as the mechanical stimulus for bone adaptation, was determined with micro-finite element models. Mechanical loading had a beneficial effect on the bone microstructure and bone stiffness in all age groups. Mineralizing surface was on average 13% greater (p<0.05) in loaded than control groups in 15- and 82-week-old mice, but not for 52-week-old mice. SED at the start of loading correlated to the change in bone volume fraction in the following 6 weeks for loaded groups (r(2)=0.69-0.85) but not control groups. At the local level, SED was 14-20% greater (p<0.01) at sites of bone formation, and 15-20% lower (p<0.01) at sites of bone resorption compared to quiescent bone surfaces for all age groups, indicating SED was a stimulus for bone adaptation. Taken together, these results support that mechanosensitivity is maintained with age in caudal vertebrae of mice at a local and global level. Since age-related bone loss was not observed in caudal vertebrae, results from the current study might not be translatable to aged humans.

Magnetic resonance imaging of the ear for patient-specific reconstructive surgery

Introduction

Like a fingerprint, ear shape is a unique personal feature that should be reconstructed with a high fidelity during reconstructive surgery. Ear cartilage tissue engineering (TE) advantageously offers the possibility to use novel 3D manufacturing techniques to reconstruct the ear, thus allowing for a detailed auricular shape. However it also requires detailed patient-specific images of the 3D cartilage structures of the patient’s intact contralateral ear (if available). Therefore the aim of this study was to develop and evaluate an imaging strategy for acquiring patient-specific ear cartilage shape, with sufficient precision and accuracy for use in a clinical setting.

Methods and Materials

Magnetic resonance imaging (MRI) was performed on 14 volunteer and six cadaveric auricles and manually segmented. Reproducibility of cartilage volume (Cg.V), surface (Cg.S) and thickness (Cg.Th) was assessed, to determine whether raters could repeatedly define the same volume of interest. Additionally, six cadaveric auricles were harvested, scanned and segmented using the same procedure, then dissected and scanned using high resolution micro-CT. Correlation between MR and micro-CT measurements was assessed to determine accuracy.

Results

Good inter- and intra-rater reproducibility was observed (precision errors <4% for Cg.S and <9% for Cg.V and Cg.Th). Intraclass correlations were good for Cg.V and Cg.S (>0.82), but low for Cg.Th (<0.23) due to similar average Cg.Th between patients. However Pearson’s coefficients showed that the ability to detect local cartilage shape variations is unaffected. Good correlation between clinical MRI and micro-CT (r>0.95) demonstrated high accuracy.

Discussion and Conclusion

This study demonstrated that precision and accuracy of the proposed method was high enough to detect patient-specific variation in ear cartilage geometry. The present study provides a clinical strategy to access the necessary information required for the production of 3D ear scaffolds for TE purposes, including detailed patient-specific shape. Furthermore, the protocol is applicable in daily clinical practice with existing infrastructure.

Gene-by-diet interactions influence calcium absorption and bone density in mice

Dietary calcium (Ca) intake is needed to attain peak bone mineral density (BMD). Habitual low Ca intake increases intestinal Ca absorption efficiency to protect bone mass, but the mechanism controlling, and the impact of genetics on, this adaptive response is not clear. We fed 11 genetically diverse inbred mouse lines a normal (0.5%) or low (0.25%) Ca diet from 4 to 12 weeks of age (n = 8 per diet per line) and studied the independent and interacting effects of diet and genetics on Ca and bone metabolism. Significant genetic variation was observed in all bone, renal, and intestinal phenotypes measured including Ca absorption. Also, adaptation of Ca absorption and bone parameters to low dietary Ca was significantly different among the lines. Ca absorption was positively correlated to femur BMD (r = 0.17, p = 0.02), and distal femur bone volume/tissue volume (BV/TV) (r = 0.34, p < 0.0001). Although Ca absorption was correlated to 1,25 dihydroxyvitamin D (1,25(OH)2 D) (r = 0.35, p < 0.0001), the adaptation of Ca absorption to low Ca intake did not correlate to diet-induced adaptation of 1,25(OH)2 D across the 11 lines. Several intestinal proteins have been proposed to mediate Ca absorption: claudins 2 and 12, voltage gated Ca channel v1.3 (Cav1.3), plasma membrane Ca ATPase 1b (PMCA1b), transient receptor potential vanilloid member 6 (TRPV6), and calbindin D9k (CaBPD9k). Only the mRNA levels for TRPV6, CaBPD9k, and PMCA1b were related to Ca absorption (r = 0.42, 0.43, and 0.21, respectively). However, a significant amount of the variation in Ca absorption is not explained by the current model and suggests that novel mechanisms remain to be determined. These observations lay the groundwork for discovery-focused initiatives to identify novel genetic factors controlling gene-by-diet interactions affecting Ca/bone metabolism.

Three-dimensional image registration improves the long-term precision of in vivo micro-computed tomographic measurements in anabolic and catabolic mouse models

Micro-computed tomography (micro-CT) is a widely used technique to track bone structural and mineral changes in small animals in vivo. Precise definition of volumes of interest (VOIs) in follow-up scans is required to accurately quantify these changes. To improve precision, VOIs can be transferred from baseline images onto follow-ups using image registration. We studied the performance of a registration procedure applied to in vivo data sets of anabolic and osteoporotic bone changes in mice. Micro-CT image data from two separate CD1 mouse data sets were studied. The first included a group treated with parathyroid hormone (PTH) and control and the second, an ovariectomy (OVX) group and control. Micro-CT was performed once per week for 4 weeks at the proximal tibia starting at treatment onset (PTH data set) or after surgery (OVX data set). A series consisting entirely of user-defined VOIs and a registered series where VOIs defined at baseline were transferred to follow-ups were created. Standard bone structural and mineral measurements were calculated. Image registration resulted in a 13-56 % reduction in precision error. Significant effects of registration to detect PTH-induced changes in BV/TV and trabecular BMD were observed. When changes were very pronounced or small, the qualitative improvement observed for the registered data set did not reach statistical significance. This study documents an increase in long-term precision of micro-CT measurements with image registration. Sensitivity to detect changes was improved but not uniform for all parameters. Future study of this technique on images with a smaller voxel size (<19 μm) may capture the effect in greater detail, in particular for trabecular thickness, where changes may be too small to be observed with the voxel size used here. Our results document the value of registration and indicate that the magnitude of improvement depends on the model and treatment chosen.

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.

Trabecular bone adapts to long-term cyclic loading by increasing stiffness and normalization of dynamic morphometric rates

Bone has the ability to adapt to external loading conditions. Especially the beneficial effect of short-term cyclic loading has been investigated in a number of in vivo animal studies. The aim of this study was to assess the long-term effect (>10 weeks) of cyclic mechanical loading on the bone microstructure, bone stiffness, and bone remodeling rates. Mice were subjected to cyclic mechanical loading at the sixth caudal vertebra with 8N or 0N (control) three times per week for a total period of 14 weeks. Structural bone parameters were determined from in vivo micro-computed tomography (micro-CT) scans performed at week 0, 4, 6, 8, 10, 12, and 14. Mechanical parameters were derived from micro-finite element analysis. Dynamic bone morphometry was calculated using registration of serial micro-CT scans. Bone volume fraction and trabecular thickness increased significantly more for the loaded group than for the control group (p = 0.006 and p = 0.002 respectively). The trabecular bone microstructure adapted to the load of 8N in approximately ten weeks, indicated by the trabecular bone volume fraction, which increased from 16.7% at 0 weeks to 21.6% at week 10 and only showed little change afterwards (bone volume fraction of 21.5% at 14 weeks). Similarly bone stiffness - (at the start of the experiment 649N/mm) - reached 846N/mm at 10 weeks in the loaded group and was maintained to the end of the experiment (850N/mm). At 4 weeks the bone formation rate was 32% greater and the bone resorption rate 22% less for 8N compared to 0N. This difference was significantly reduced as the bone adapted to 8N, with 8N remodeling rates returning to the values of the 0N group at approximately 10 weeks. Together these data suggest that once bone has adapted to a new loading state, the remodeling rates reduce gradually while maintaining bone volume fraction and stiffness.

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.

Longitudinal in vivo imaging of bone formation and resorption using fluorescence molecular tomography

Bone research often focuses on anatomical imaging of the bone microstructure, but in order to gain better understanding in how bone remodeling is modulated through interventions also bone formation and resorption processes should be investigated. With this in mind, the purpose of this study was to establish a longitudinal in vivo imaging approach of bone formation and resorption using fluorescence molecular tomography (FMT). In this study the reproducibility, accuracy and sensitivity of FMT for bone imaging were assessed by performing longitudinal measurements with FMT and comparing it to in vivo micro-computed tomography on a set of control mice, and mice in which load-adaptation was induced in the sixth caudal vertebra. The precision error for FMT measurements, expressed as coefficient of variation, was smaller than 16%, indicating acceptable reproducibility. A correlation was found between bone resorption measured with FMT and bone resorption rate measured with in vivo micro-computed tomography only over the first 14days (R=0.81, p<0.01), but not between bone formation measured with FMT and bone formation rate measured with in vivo micro-CT. Bone formation measured by FMT was 89-109% greater (p<0.05) for mice subjected to mechanical loading than control mice. Bone resorption was 5-8% lower, but did not reach a significant difference between groups, indicating moderate sensitivity for FMT. In conclusion, in vivo FMT in mouse tail bones is feasible but needs to be optimized for monitoring load adaptation in living mice.

A comparison of the thresholding strategies of micro-CT for periodontal bone loss: a pilot study

OBJECTIVES

Micro-CT provides three-dimensional details and has been widely used for biomedical assessments. This study aimed to determine the most appropriate threshold method for quantitatively assessing the dynamics of periodontal destruction.

METHODS

Inflammation was induced by submerging a silk ligature in the sulcus of the maxillary second molars of rats, and the animals were killed prior to ligature placement and after 7 and 21 days. The maxillae were examined for the bone resorptive activities by micro-CT, histology and tartrate-resistant acid phosphatase staining. The imaging threshold was determined by CT phantom, global and local algorithms. A bone fraction measurement from each threshold-determining technique was compared with histomorphometry. The reliability and reproducibility were examined by the intraclass correlation coefficient (ICC) and the coefficient of variation.

RESULTS

Significant reduction of inflammatory infiltration (p < 0.01) and active osteoclastic resorption (p < 0.05) from Day 7 to Day 21 were noted. High inter- and intraexaminer agreement were demonstrated in both histomorphometric and micro-CT assessments (ICC > 0.98). The algorithm-based technique demonstrated stronger correlation to histomorphometry than phantom-based thresholds, and the highest agreement was presented by the local algorithm (ICC > 0.96). This, however, was considerably computationally expensive.

CONCLUSIONS

The local threshold-determining algorithm is suggested for examining inflammation-induced bone loss. Further investigation will be aimed at enhancing computational efficiency.

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.

MicroCT morphometry analysis of mouse cancellous bone: intra- and inter-system reproducibility

The agreement between measurements and the relative performance reproducibility among different microcomputed tomography (microCT) systems, especially at voxel sizes close to the limit of the instruments, is not known. To compare this reproducibility 3D morphometric analyses of mouse cancellous bone from distal femoral epiphyses were performed using three different ex vivo microCT systems: GE eXplore Locus SP, Scanco μCT35 and Skyscan 1172. Scans were completed in triplicate at 12 μm and 8 μm voxel sizes and morphometry measurements, from which relative values and dependence on voxel size were examined. Global and individual visually assessed thresholds were compared. Variability from repeated scans at 12 μm voxel size was also examined. Bone volume fraction and trabecular separation values were similar, while values for relative bone surface, trabecular thickness and number varied significantly across the three systems. The greatest differences were measured in trabecular thickness (up to 236%) and number (up to 218%). The relative dependence of measurements on voxel size was highly variable for the trabecular number (from 0% to 20% relative difference between measurements from 12 μm and 8 μm voxel size scans, depending on the system). The intra-system reproducibility of all trabecular measurements was also highly variable across the systems and improved for BV/TV in all the systems when a smaller voxel size was used. It improved using a smaller voxel size in all the other parameters examined for the Scanco system, but not consistently so for the GE or the Skyscan system. Our results indicate trabecular morphometry measurements should not be directly compared across microCT systems. In addition, the conditions, including voxel size, for trabecular morphometry studies in mouse bone should be chosen based on the specific microCT system and the measurements of main interest.

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.

Disruption of bone development and homeostasis by trisomy in Ts65Dn Down syndrome mice

Down syndrome (DS) is a genetic disorder resulting from trisomy 21 that causes cognitive impairment, low muscle tone and craniofacial alterations. Morphometric studies of the craniofacial and appendicular skeleton in individuals with DS suggest that bone development and homeostasis are affected by trisomy. The Ts65Dn mouse model has three copies of approximately half the genes found on human chromosome 21 and exhibits craniofacial skeletal and size differences similar to those observed in humans with DS. We hypothesized that Ts65Dn and euploid mice have distinct differences in bone development and homeostasis influencing both the craniofacial and appendicular skeletal phenotypes. Quantitative assessment of structural and mechanical properties of the femur in Ts65Dn and control mice at 6 and 16 weeks of age revealed significant deficiencies in trabecular and cortical bone architecture, bone mineral density, bone formation, and bone strength in trisomic bone. Furthermore, bone mineral density and dynamic dentin formation rate of the skull and incisor, respectively, were also reduced in Ts65Dn mice, demonstrating that trisomy significantly affects both the craniofacial and appendicular skeleton.

Towards quantitative 3D imaging of the osteocyte lacuno-canalicular network

Osteocytes are the most abundant cells in bone and the only cells embedded in the bone mineral matrix. They form an extended, three-dimensional (3D) network, whose processes interconnecting the cell bodies reside in thin canals, the canaliculi. Together with the osteocyte lacunae, the canaliculi form the lacuno-canalicular network (LCN). As the negative imprint of the cellular network within bone tissue, the LCN morphology is considered to play a central role for bone mechanosensation and mechanotransduction. However, the LCN has neither been visualized nor quantified in an adequate way up to now. On this account, this article summarizes the current state of knowledge of the LCN morphology and then reviews different imaging methods regarding the quantitative 3D assessment of bone tissue in general and of the LCN in particular. These imaging methods will provide new insights in the field of bone mechanosensation and mechanotransduction and thus, into processes of strain sensation and transduction, which are tightly associated with osteocyte viability and bone quality.

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.

Beam hardening artifacts in micro-computed tomography scanning can be reduced by X-ray beam filtration and the resulting images can be used to accurately measure BMD

Bone mineral density (BMD) measurements are critical in many research studies investigating skeletal integrity. For pre-clinical research, micro-computed tomography (microCT) has become an essential tool in these studies. However, the ability to measure the BMD directly from microCT images can be biased by artifacts, such as beam hardening, in the image. This three-part study was designed to understand how the image acquisition process can affect the resulting BMD measurements and to verify that the BMD measurements are accurate. In the first part of this study, the effect of beam hardening-induced cupping artifacts on BMD measurements was examined. In the second part of this study, the number of bones in the X-ray path and the sampling process during scanning was examined. In the third part of this study, microCT-based BMD measurements were compared with ash weights to verify the accuracy of the measurements. The results indicate that beam hardening artifacts of up to 32.6% can occur in sample sizes of interest in studies investigating mineralized tissue and affect mineral density measurements. Beam filtration can be used to minimize these artifacts. The results also indicate that, for murine femora, the scan setup can impact densitometry measurements for both cortical and trabecular bone and morphologic measurements of trabecular bone. Last, when a scan setup that minimized all of these artifacts was used, the microCT-based measurements correlated well with ash weight measurements (R(2)=0.983 when air was excluded), indicating that microCT can be an accurate tool for murine bone densitometry.

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.