Automated simulation of areal bone mineral density assessment in the distal radius from high-resolution peripheral quantitative computed tomography

SummaryAn automated image processing method is presented for simulating areal bone mineral density measures using high-resolution peripheral quantitative computed tomography (HR-pQCT) in the ultra-distal radius. The accuracy of the method is validated against clinical dual X-ray absorptiometry (DXA). This technique represents a useful reference to gauge the utility of novel 3D quantification methods applied to HR-pQCT in multi-center clinical studies and potentially negates the need for separate forearm DXA measurements.IntroductionOsteoporotic status is primarily assessed by measuring areal bone mineral density (aBMD) using 2D dual X-ray absorptiometry (DXA). However, this technique does not sufficiently explain bone strength and fracture risk. High-resolution peripheral quantitative computed tomography (HR-pQCT) has been introduced as a method to quantify 3D bone microstructure and biomechanics. In this study, an automated method is proposed to simulate aBMD measures from HR-pQCT distal radius images.MethodsA total of 117 subject scans were retrospectively analyzed from two clinical bone quality studies. The distal radius was imaged by HR-pQCT and DXA on one of two devices (Hologic or Lunar). Areal BMD was calculated by simulation from HR-pQCT images (aBMD(sim)) and by standard DXA analysis (aBMD(dxa)).ResultsThe reproducibility of the simulation technique was 1.1% (root mean-squared coefficient of variation). HR-pQCT-based aBMD(sim) correlated strongly to aBMD(dxa) (Hologic: R (2) = 0.82, Lunar: R (2) = 0.87), though aBMD(sim) underestimated aBMD(dxa) for both DXA devices (p < 0.0001). Finally, aBMD(sim) predicted aBMD at the proximal femur and lumbar spine with equal power compared to aBMD(dxa).ConclusionThe results demonstrate that aBMD can be simulated from HR-pQCT images of the distal radius. This approach has the potential to serve as a surrogate forearm aBMD measure for clinical HR-pQCT studies when axial bone mineral density values are not required

HR‐pQCT measures of bone microarchitecture predict fracture : systematic review and meta‐analysis

HR‐pQCT is a non‐invasive imaging modality for assessing volumetric bone mineral density (vBMD) and microarchitecture of cancellous and cortical bone. The objective was to (i) assess fracture‐associated differences in HR‐pQCT bone parameters and (ii) to determine if HR‐pQCT is sufficiently precise to reliably detect these differences in individuals. We systematically identified 40 studies that used HR‐pQCT (39/40 used XtremeCT scanners) to assess 1291‐3253 and 3389‐10,687 individuals with and without fractures, respectively, ranging in age from 10.9 to 84.7 years with no comorbid conditions. Parameters describing radial and tibial bone density, microarchitecture, and strength were extracted and percentage differences between fracture and control subjects were estimated using a random effects meta‐analysis. An additional meta‐analysis of short‐term in vivo reproducibility of bone parameters assessed by XtremeCT was conducted to determine whether fracture‐associated differences exceeded the least significant change (LSC) required to discern measured differences from precision error. Radial and tibial HR‐pQCT parameters, including failure load, were significantly altered in fracture subjects, with differences ranging from −2.6% (95% CI: −3.4 to −1.9) in radial cortical vBMD to −12.6% (95% CI: −15.0 to −10.3) in radial trabecular vBMD. Fracture‐associated differences reported by prospective studies were consistent with those from retrospective studies, indicating that HR‐pQCT can predict incident fracture. Assessment of study quality, heterogeneity and publication biases verified the validity of these findings. Finally, we demonstrated that fracture‐associated deficits in total and trabecular vBMD, and certain tibial cortical parameters, can be reliably discerned from HR‐pQCT‐related precision error and can be used to detect fracture‐associated differences in individual patients. Although differences in other HR‐pQCT measures, including failure load, were significantly associated with fracture, improved reproducibility is needed to ensure reliable individual cross‐sectional screening and longitudinal monitoring. In conclusion, our study supports the use of HR‐pQCT in clinical fracture prediction

ORMIR_XCT: A Python package for high resolution peripheral quantitative computed tomography image processing

High resolution peripheral quantitative computed tomography (HR-pQCT) is an imaging technique capable of imaging trabecular bone in-vivo. HR-pQCT has a wide range of applications, primarily focused on bone to improve our understanding of musculoskeletal diseases, assess epidemiological associations, and evaluate the effects of pharmaceutical interventions. Processing HR-pQCT images has largely been supported using the scanner manufacturer scripting language (Image Processing Language, IPL, Scanco Medical). However, by expanding image processing workflows outside of the scanner manufacturer software environment, users have the flexibility to apply more advanced mathematical techniques and leverage modern software packages to improve image processing. The ORMIR_XCT Python package was developed to reimplement some existing IPL workflows and provide an open and reproducible package allowing for the development of advanced HR-pQCT data processing workflows

Development of a novel method to measure bone marrow fat fraction in older women using high-resolution peripheral quantitative computed tomography

Bone marrow adipose tissue (BMAT) has been implicated in a number of conditions associated with bone deterioration and osteoporosis. Several studies have found an inverse relationship between BMAT and bone mineral density (BMD), and higher levels of BMAT in those with prevalent fracture. Magnetic resonance imaging (MRI) is the gold standard for measuring BMAT, but its use is limited by high costs and low availability. We hypothesized that BMAT could also be accurately quantified using high-resolution peripheral quantitative computed tomography (HR-pQCT). Methods: In the present study, a novel method to quantify the tibia bone marrow fat fraction, defined by MRI, using HR-pQCT was developed. In total, 38 postmenopausal women (mean [standard deviation] age 75.9 [3.1] years) were included and measured at the same site at the distal (n = 38) and ultradistal (n = 18) tibia using both MRI and HR-pQCT. To adjust for partial volume effects, the HR-pQCT images underwent 0 to 10 layers of voxel peeling to remove voxels adjacent to the bone. Linear regression equations were then tested for different degrees of voxel peeling, using the MRI-derived fat fractions as the dependent variable and the HR-pQCT-derived radiodensity as the independent variables. Results: The most optimal HR-pQCT derived model, which applied a minimum of 4 layers of peeled voxel and with more than 1% remaining marrow volume, was able to explain 76% of the variation in the ultradistal tibia bone marrow fat fraction, measured with MRI (p < 0.001). Conclusion: The novel HR-pQCT method, developed to estimate BMAT, was able to explain a substantial part of the variation in the bone marrow fat fraction and can be used in future studies investigating the role of BMAT in osteoporosis and fracture prediction

High-Impact Mechanical Loading Increases Bone Material Strength in Postmenopausal Women-A 3-Month Intervention Study.

Bone adapts to loading in several ways, including redistributing bone mass and altered geometry and microarchitecture. Because of previous methodological limitations, it is not known how the bone material strength is affected by mechanical loading in humans. The aim of this study was to investigate the effect of a 3-month unilateral high-impact exercise program on bone material properties and microarchitecture in healthy postmenopausal women. A total of 20 healthy and inactive postmenopausal women (aged 55.6 ± 2.3 years [mean ± SD]) were included and asked to perform an exercise program of daily one-legged jumps (with incremental number, from 3×10 to 4×20 jumps/d) during 3 months. All participants were asked to register their performed jumps in a structured daily diary. The participants chose one leg as the intervention leg and the other leg was used as control. The operators were blinded to the participant's choice of leg for intervention. The predefined primary outcome was change in bone material strength index (BMSi), measured at the mid tibia with a handheld reference probe indentation instrument (OsteoProbe). Bone microstructure, geometry, and density were measured with high-resolution peripheral quantitative computed tomography (XtremeCT) at the ultradistal and at 14% of the tibia bone length (distal). Differences were analyzed by related samples Wilcoxon signed rank test. The overall compliance to the jumping program was 93.6%. Relative to the control leg, BMSi of the intervention leg increased 7% or 0.89 SD (p = 0.046), but no differences were found for any of the XtremeCT-derived bone parameters. In conclusion, a unilateral high-impact loading program increased BMSi in postmenopausal women rapidly without affecting bone microstructure, geometry, or density, indicating that intense mechanical loading has the ability to rapidly improve bone material properties before changes in bone mass or structure. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals Inc

Multi-modality CT imaging of human bone for improved validation of subject-specific finite element analysis

INTRODUCTION Finite element analysis (FE) is a promising alternative to dual-energy x-ray absorptiometry for improved prediction of fracture risk as FE can incorporate bone density, geometry, and microarchitecture [1]. However, the accuracy of FE models is heavily influenced by the spatial resolution of the CT scan [2].  In this study we quantify the architectural differences in trabecular bone between different modalities, specifically high-resolution peripheral quantitative computed tomography (HR-pQCT) and micro-computed tomography (μCT), and measure their effect on the resulting FE models. METHODS Two cadaveric, human distal tibiae were imaged using HR-pQCT (XtremeCT, Scanco Medical, Switzerland), with a voxel size of 82μm. After imaging, four cubes, 10mm edge length, were cut and scanned using μCT (μCT-35, Scanco Medical, Switzerland.), with a voxel size of 20μm. During a sensitivity analysis, the μCT image data was segmented using a threshold based technique and resampled to voxel sizes of 40μm, 60μm, and 80μm to assess the effect of voxel size on the FE results. The FE results were statistical compared with a one-way repeated measures ANOVA and a Tukey’s post-hoc test. Based on the known location of each cube recorded during the cutting process, the μCT data was manually aligned to the larger HR-pQCT image, where mutual information registration was applied for accurate alignment.  Virtual cubes were extracted from the registered HR-pQCT data. All cube image data was rescaled to preserve the bone volume/total volume ratio. Subsequently, all image data was converted to hexahedron elements for FE analysis and subjected to 1% uniaxial compression (FAIM v6.0, Numerics88) in the x-, y-, and z-directions. The resulting FE data was compared with a two-way ANOVA, with Bonferroni multiple comparisons test. RESULTS The sensitivity analysis found no statistically significant differences in reaction force between any cubes with different resolutions when compressed in the z-direction. The mean percent error, when compared to the 20μm cube, for the 40μm, 60μm, and 80μm cubes was 0.06%, 0.45%, and 0.76%, respectively. In the y-direction, there were significant differences between the 20μm and the 80μm (p < 0.01; 62.51% error). When loaded in the x-direction, there were significant differences between the 20μm FE cube and the 60μm and 80μm FE cubes (p<0.05; respective mean percent errors of 121.27% and 182.83%). With increasing voxel size, the reaction force was overestimated. As shown in figure 1, the registration between μCT and HR-pQCT was successful. However, there were statistically significant differences found between the μCT versus HR-pQCT FE data in the y- and x-direction (p<0.01; 23.95% and 21.04% error, respectively) but none were found for the z-direction (p>0.05, 0.14% error.) Figure 1. Registered distal tibia (HR-pQCT - red) to the μCT acquired cubes (yellow). The registered image appears green in overlap. DISCUSSION AND CONCLUSIONS It was shown there is little difference in the FE results loaded in the z-direction upon rescaling μCT data (mean percent errors less than 0.8%). However, there are large changes in the non-axial (i.e. x and y) directions. It is recommended that μCT data be rescaled up to 40μm. In addition, it was shown that μCT data could be successfully registered to HR-pQCT (see figure 1). For a given resolution, μCT data is only comparable to HR-pQCT data in the z-direction. In conclusion, HR-pQCT may not be an ideal substitute for μCT micro-architectural data.  This study improves our understanding on the use of HR-pQCT imaging for quantification of bone micro-architecture

High-Resolution Cone-Beam Computed Tomography is a Fast and Promising Technique to Quantify Bone Microstructure and Mechanics of the Distal Radius

Obtaining high-resolution scans of bones and joints for clinical applications is challenging. HR-pQCT is considered the best technology to acquire high-resolution images of the peripheral skeleton in vivo, but a breakthrough for widespread clinical applications is still lacking. Recently, we showed on trapezia that CBCT is a promising alternative providing a larger FOV at a shorter scanning time. The goals of this study were to evaluate the accuracy of CBCT in quantifying trabecular bone microstructural and predicted mechanical parameters of the distal radius, the most often investigated skeletal site with HR-pQCT, and to compare it with HR-pQCT. Nineteen radii were scanned with four scanners: (1) HR-pQCT (XtremeCT, Scanco Medical AG, @ (voxel size) 82 mu m), (2) HR-pQCT (XtremeCT-II, Scanco, @60.7 mu m), (3) CBCT (NewTom 5G, Cefla, @75 mu m) reconstructed and segmented using in-house developed software and (4) microCT (VivaCT40, Scanco, @19 mu m-gold standard). The following parameters were evaluated: predicted stiffness, strength, bone volume fraction (BV/TV) and trabecular thickness (Tb.Th), separation (Tb.Sp) and number (Tb.N). The overall accuracy of CBCT with in-house optimized algorithms in quantifying bone microstructural parameters was comparable (R-2 = 0.79) to XtremeCT (R-2 = 0.76) and slightly worse than XtremeCT-II (R-2 = 0.86) which were both processed with the standard manufacturer's technique. CBCT had higher accuracy for BV/TV and Tb.Th but lower for Tb.Sp and Tb.N compared to XtremeCT. Regarding the mechanical parameters, all scanners had high accuracy (R-2 >= 0.96). While HR-pQCT is optimized for research, the fast scanning time and good accuracy renders CBCT a promising technique for high-resolution clinical scanning

Clinical observation of diminished bone quality and quantity through longitudinal HR-pQCT-derived remodeling and mechanoregulation.

High resolution peripheral quantitative computed tomography (HR-pQCT) provides methods for quantifying volumetric bone mineral density and microarchitecture necessary for early diagnosis of bone disease. When combined with a longitudinal imaging protocol and finite element analysis, HR-pQCT can be used to assess bone formation and resorption (i.e., remodeling) and the relationship between this remodeling and mechanical loading (i.e., mechanoregulation) at the tissue level. Herein, 25 patients with a contralateral distal radius fracture were imaged with HR-pQCT at baseline and 9-12 months follow-up: 16 patients were prescribed vitamin D3 with/without calcium supplement based on a blood biomarker measures of bone metabolism and dual-energy X-ray absorptiometry image-based measures of normative bone quantity which indicated diminishing (n = 9) or poor (n = 7) bone quantity and 9 were not. To evaluate the sensitivity of this imaging protocol to microstructural changes, HR-pQCT images were registered for quantification of bone remodeling and image-based micro-finite element analysis was then used to predict local bone strains and derive rules for mechanoregulation. Remodeling volume fractions were predicted by both average values of trabecular and cortical thickness and bone mineral density (R2 > 0.8), whereas mechanoregulation was affected by dominance of the arm and group classification (p < 0.05). Overall, longitudinal, extended HR-pQCT analysis enabled the identification of changes in bone quantity and quality too subtle for traditional measures

In Vivo Evaluation of the Presence of Bone Marrow in Cortical Porosity in Postmenopausal Osteopenic Women

This is the first observational study examining cortical porosity in vivo in postmenopausal osteopenic women and to incorporate data from two different imaging modalities to further examine the nature of cortical porosity. The goal of this study was to combine high-resolution peripheral computed tomography (HR-pQCT) images, which contain high spatial resolution information of the cortical structure, and magnetic resonance (MR) images, which allow the visualization of soft tissues such as bone marrow, to observe the amount of cortical porosity that contains bone marrow in postmenopausal osteopenic women. The radius of 49 and the tibia of 51 postmenopausal osteopenic women (age 56 ± 3.7) were scanned using both HR-pQCT and MR imaging. A normalized mutual information registration algorithm was used to obtain a three-dimensional rigid transform which aligned the MR image to the HR-pQCT image. The aligned images allowed for the visualization of bone marrow in cortical pores. From the HR-pQCT image, the percent cortical porosity, the number of cortical pores, and the size of each cortical pore was determined. By overlaying the aligned MR and HR-pQCT images, the percent of cortical pores containing marrow, the number of cortical pores containing marrow, and the size of each cortical pore containing marrow were measured. While the amount of cortical porosity did not vary greatly between subjects, the type of cortical pore, containing marrow vs. not containing marrow, varied highly between subjects. The results suggest that cortical pore spaces contain components of varying composition, and that there may be more than one mechanism for the development of cortical porosity