Chamith S. Rajapakse

Accuracy of High-Resolution In Vivo Micro Magnetic Resonance Imaging for Measurements of Microstructural and Mechanical Properties of Human Distal Tibial Bone

Micro magnetic resonance imaging (µMRI) is an in vivo imaging method that permits 3D quantification of cortical and trabecular bone microstructure. µMR images can also be used for building microstructural finite element (µFE) models to assess bone stiffness, which highly correlates with bones resistance to fractures. In order for µMRI‐based microstructural and µFE analyses to become standard clinical tools for assessing bone quality, validation with a current gold standard, namely, high‐resolution micro computed tomography (µCT), is required. Microstructural measurements of 25 human cadaveric distal tibias were performed for the registered µMR and µCT images, respectively. Next, whole bone stiffness, trabecular bone stiffness, and elastic moduli of cubic subvolumes of trabecular bone in both µMR and µCT images were determined by voxel‐based µFE analysis. The bone volume fraction (BV/TV), trabecular number (Tb.N*), trabecular spacing (Tb.Sp*), cortical thickness (Ct.Th), and structure model index (SMI) based on µMRI showed strong correlations with µCT measurements (r2 = 0.67 to 0.97), and bone surface‐to‐volume ratio (BS/BV), connectivity density (Conn.D), and degree of anisotropy (DA) had significant but moderate correlations (r2 = 0.33 to 0.51). Each of these measurements also contributed to one or many of the µFE‐predicted mechanical properties. However, model‐independent trabecular thickness (Tb.Th*) based on µMRI had no correlation with the µCT measurement and did not contribute to any mechanical measurement. Furthermore, the whole bone and trabecular bone stiffness based on µMRI were highly correlated with those of µCT images (r2 = 0.86 and 0.96), suggesting that µMRI‐based µFE analyses can directly and accurately quantify whole bone mechanical competence. In contrast, the elastic moduli of the µMRI trabecular bone subvolume had significant but only moderate correlations with their gold standards (r2 = 0.40 to 0.58). We conclude that most microstructural and mechanical properties of the distal tibia can be derived efficiently from µMR images and can provide additional information regarding bone quality.

In Vivo μMRI-Based Finite Element and Morphological Analyses of Tibial Trabecular Bone in Eugonadal and Hypogonadal Men Before and After Testosterone Treatment†‡

Osteoporosis is a major public health problem in men. Hypogonadal men have decreased BMD and deteriorated trabecular bone architecture compared with eugonadal men. Testosterone treatment improves their BMD and trabecular structure. We tested the hypothesis that testosterone replacement in hypogonadal men would also improve their bones mechanical properties. Ten untreated severely hypogonadal and 10 eugonadal men were selected. The hypogonadal men were treated with a testosterone gel for 24 mo to maintain their serum testosterone concentrations within the normal range. Each subject was assessed before and after 6, 12, and 24 mo of testosterone treatment by μMRI of the distal tibia. A subvolume of each μMR image was converted to a microfinite element (μFE) model, and six analyses were performed, representing three compression and three shear tests. The anisotropic stiffness tensor was calculated, from which the orthotropic elastic material constants were derived. Changes in microarchitecture were also quantified using newly developed individual trabeculae segmentation (ITS)‐based and standard morphological analyses. The accuracy of these techniques was examined with simulated μMR images. Significant differences in four estimated anisotropic elastic material constants and most morphological parameters were detected between the eugonadal and hypogonadal men. No significant change in estimated elastic moduli and morphological parameters was detected in the eugonadal group over 24 mo. After 24 mo of treatment, significant increases in estimated elastic moduli E22 (9.0%), E33 (5.1%), G23 (7.2%), and G12 (9.4%) of hypogonadal men were detected. These increases were accompanied by significant increases in trabecular plate thickness. These results suggest that 24 mo of testosterone treatment of hypogonadal men improves estimated elastic moduli of tibial trabecular bone by increased trabecular plate thickness.

Micro–MR Imaging–based Computational Biomechanics Demonstrates Reduction in Cortical and Trabecular Bone Strength after Renal Transplantation

PURPOSE To examine the ability of three-dimensional micro-magnetic resonance (MR) imaging-based computational biomechanics to detect mechanical alterations in trabecular bone and cortical bone in the distal tibia of incident renal transplant recipients 6 months after renal transplantation and compare them with bone mineral density (BMD) outcomes. MATERIALS AND METHODS The study was approved by the institutional review board and complied with HIPAA guidelines. Written informed consent was obtained from all subjects. Micro-MR imaging of distal tibial metaphysis was performed within 2 weeks after renal transplantation (baseline) and 6 months later in 49 participants (24 female; median age, 44 years; range, 19-61 years) with a clinical 1.5-T whole-body imager using a modified three-dimensional fast large-angle spin-echo pulse sequence. Micro-finite-element models for cortical bone, trabecular bone, and whole-bone section were generated from each image by delineating the endosteal and periosteal boundaries. Mechanical parameters (stiffness and failure load) were estimated with simulated uniaxial compression tests on the micro-finite-element models. Structural parameters (trabecular bone volume fraction [BV/TV, bone volume to total volume ratio], trabecular thickness [TbTh], and cortical thickness [CtTh]) were computed from micro-MR images. Total hip and spine areal BMD were determined with dual-energy x-ray absorptiometry (DXA). Parameters obtained at the follow-up were compared with the baseline values by using parametric or nonparametric tests depending on the normality of data. RESULTS All mechanical parameters were significantly lower at 6 months compared with baseline. Decreases in cortical bone, trabecular bone, and whole-bone stiffness were 3.7% (P = .03), 4.9% (P = .03), and 4.3% (P = .003), respectively. Decreases in cortical bone, trabecular bone, and whole-bone failure strength were 7.6% (P = .0003), 6.0% (P = .004), and 5.6% (P = .0004), respectively. Conventional structural measures, BV/TV, TbTh, and CtTh, did not change significantly. Spine BMD decreased by 2.9% (P < .0001), while hip BMD did not change significantly at DXA. CONCLUSION MR imaging-based micro-finite-element analysis suggests that stiffness and failure strength of the distal tibia decrease over a 6-month interval after renal transplantation.

Computational biomechanics of the distal tibia from high-resolution MR and micro-CT images

The mechanical properties of bone estimated by micro-finite element (microFE) analysis on the basis of in vivo micro-MR images (microMRIs) of the distal extremities provide a new tool for direct assessment of the mechanical consequences of intervention. However, the accuracy of the method has not previously been investigated. Here, we compared microFE-derived mechanical parameters obtained from microMRIs at 160 microm isotropic voxel size now achievable in vivo with those derived from 25 microm isotropic (reference) microCT images of 30 cadaveric tibiae from 15 donors (4 females and 11 males, aged 55-84 years). Elastic and shear moduli estimated from 5mm(3) subvolumes of trabecular bone (TB) derived from microMRIs were significantly correlated with those derived from volume-matched reference microCT images (R(2)=0.60-0.67). Axial stiffness of whole-bone sections (including both cortical and trabecular compartments) derived from microMR-based models were highly correlated (R(2)=0.85) with those from high-resolution reference images. Further, microFE models generated from microCT images after downsampling to lower resolutions relevant to in vivo microMRI (100-160 microm) showed mechanical parameters to be strongly correlated (R(2)>0.93) with those derived at reference resolution (25 microm). Incorporation of grayscale image information into the microMR-based microFE model yielded slopes closer to unity than binarized models (1.07+/-0.15 vs. 0.71+/-0.11) when correlated with reference subregional elastic and shear moduli. This work suggests that elastic properties of distal tibia can be reliably estimated by microFE analysis from microMRIs obtainable at in vivo resolution.

Adverse Fat Depots and Marrow Adiposity Are Associated With Skeletal Deficits and Insulin Resistance in Long-Term Survivors of Pediatric Hematopoietic Stem Cell Transplantation.

Allogeneic hematopoietic stem‐cell transplantation (alloHSCT) survivors treated with total body irradiation (TBI) exhibit bone deficits and excess adiposity, potentially related to altered mesenchymal stem cell differentiation into osteoblasts or adipocytes. We examined associations among fat distribution, bone microarchitecture, and insulin resistance in alloHSCT survivors after TBI. This was a cross‐sectional observational study of 25 alloHSCT survivors (aged 12 to 25 years) a median of 9.7 (4.3 to 19.3) years after alloHSCT compared to 25 age‐, race‐, and sex‐matched healthy controls. Vertebral MR spectroscopic imaging and tibia micro‐MRI were used to quantify marrow adipose tissue (MAT) and trabecular microarchitecture. Additional measures included DXA whole‐body fat mass (WB‐FM), leg lean mass (Leg‐LM), trunk visceral adipose tissue (VAT), and CT calf muscle density. Insulin resistance in alloHSCT survivors was estimated by HOMA‐IR. AlloHSCT survivors had lower Leg‐LM (p < 0.001) and greater VAT (p < 0.01), MAT (p < 0.001), and fat infiltration of muscle (p = 0.04) independent of WB‐FM, versus matched controls; BMI did not differ. Survivors had lower bone volume fraction and abnormal microarchitecture including greater erosion and more rod‐like structure versus controls (all p = 0.04); 14 had vertebral deformities and two had compression fractures. Greater WB‐FM, VAT, MAT, and muscle fat infiltration were associated with abnormal trabecular microarchitecture (p < 0.04 for all). AlloHSCT HOMA‐IR was elevated, associated with younger age at transplantation (p < 0.01), and positively correlated with WB‐FM and VAT (both p < 0.01). In conclusion, the markedly increased marrow adiposity, abnormal bone microarchitecture, and abnormal fat distribution highlight the risks of long‐term treatment‐related morbidity and mortality in alloHSCT recipients after TBI. Trabecular deterioration was associated with marrow and visceral adiposity. Furthermore, long‐term survivors demonstrated sarcopenic obesity, insulin resistance, and vertebral deformities. Future studies are needed to identify strategies to prevent and treat metabolic and skeletal complications in this growing population of childhood alloHSCT survivors.

Mechanical Implications of Estrogen Supplementation in Early Postmenopausal Women

Whereas the structural implications of drug intervention are well established, there are few data on the possible mechanical consequences of treatment. In this work we examined the changes in elastic and shear moduli (EM and SM) in a region of trabecular bone in the distal radius and distal tibia of early postmenopausal women on the basis of MRI‐based micro‐finite‐element (µFE) analysis. Whole‐section axial stiffness (AS) encompassing both trabecular and cortical compartments was evaluated as well. The study was conducted on previously acquired high‐resolution images at the two anatomic sites. Images were processed to yield a 3D voxel array of bone‐volume fraction (BVF), which was converted to a µFE model of hexahedral elements in which tissue modulus was set proportional to voxel BVF. The study comprised 65 early postmenopausal women (age range 45 to 55 years), of whom 32 had chosen estrogen supplementation (estradiol group); the remainder had not (control group). Subjects had been scanned at baseline and 12 and 24 months thereafter. At the distal tibia, EM and SM were reduced by 2.9% to 5.5% in the control group (p < .05 to <.005), but there was no change in the estradiol subjects. AS decreased 3.9% (4.0%) in controls (p < .005) and increased by 5.8% (6.2%) in estradiol group subjects (p < .05) at 12 (24) months. At the distal radius, EM and SM changes from baseline were not significant, but at both time points AS was increased in estradiol group subjects and decreased in controls (p < .005 to <.05), albeit by a smaller margin than at the tibia. EM and SM were strongly correlated with BV/TV (r2 = 0.44 to 0.92) as well as with topologic parameters expressing the ratio of plates to rods (r2 = 0.45 to 0.82), jointly explaining up to 96% of the variation in the mechanical parameters. Finally, baseline AS was strongly correlated between the two anatomic sites (r2 = 0.58), suggesting that intersubject variations in the bones mechanical competence follows similar mechanisms. In conclusion, the results demonstrate that micro‐MRI‐based µFE models are suited for the study of the mechanical implications of antiresorptive treatment. The data further highlight the anabolic effect of short‐term estrogen supplementation.

Structural and mechanical parameters of trabecular bone estimated from in vivo high-resolution magnetic resonance images at 3 tesla field strength

To evaluate the performance of a new 3 Tesla (T) high‐resolution trabecular bone (TB) imaging technique at two resolution regimens in terms of serial reproducibility and sensitivity.

Cortical Bone Water Concentration: Dependence of MR Imaging Measures on Age and Pore Volume Fraction

PURPOSE To quantify bulk bone water to test the hypothesis that bone water concentration (BWC) is negatively correlated with bone mineral density (BMD) and is positively correlated with age, and to propose the suppression ratio (SR) (the ratio of signal amplitude without to that with long-T2 suppression) as a potentially stronger surrogate measure of porosity, which is evaluated ex vivo and in vivo. MATERIALS AND METHODS Human subject studies were conducted in compliance with institutional review board and HIPAA regulations. Healthy men and women (n = 72; age range, 20-80 years) were examined with a hybrid radial ultrashort echo time magnetic resonance (MR) imaging sequence at 3.0 T, and BWC was determined in the tibial midshaft. In a subset of 40 female subjects, the SR was measured with a similar sequence. Cortical volumetric BMD (vBMD) was measured by means of peripheral quantitative computed tomography (CT). The method was validated against micro-CT-derived porosity in 13 donor human cortical bone specimens. Associations among parameters were evaluated by using standard statistical tools. RESULTS BWC was positively correlated with age (r = 0.52; 95% confidence interval [CI]: 0.22, 0.73; P = .002) and negatively correlated with vBMD at the same location (r = -0.57; 95% CI: -0.76, -0.29; P < .001). Data were suggestive of stronger associations with SR (r = 0.64, 95% CI: 0.39, 0.81, P < .001 for age; r = -0.67, 95% CI: -0.82, -0.43, P < .001 for vBMD; P < .001 for both), indicating that SR may be a more direct measure of porosity. This interpretation was supported by ex vivo measurements showing SR to be strongly positively correlated with micro-CT porosity (r = 0.88; 95% CI: 0.64, 0.96; P < .001) and with age (r = 0.87; 95% CI: 0.62, 0.96; P < .001). CONCLUSION The MR imaging-derived SR may serve as a biomarker for cortical bone porosity that is potentially superior to BWC, but corroboration in larger cohorts is indicated.

Performance of the MRI-based virtual bone biopsy in the distal radius: Serial reproducibility and reliability of structural and mechanical parameters in women representative of osteoporosis study populations

Serial reproducibility and reliability critically determine sensitivity to detect changes in response to intervention and provide a basis for sample size estimates. Here, we evaluated the performance of the MRI-based virtual bone biopsy in terms of 26 structural and mechanical parameters in the distal radius of 20 women in the age range of 50 to 75 years (mean=62.0 years, S.D.=8.1 years), representative of typical study populations in drug intervention trials and fracture studies. Subjects were examined three times at average intervals of 20.2 days (S.D.=14.5 days) by MRI at 1.5 T field strength at a voxel size of 137×137×410 μm(3). Methods involved prospective and retrospective 3D image registration and auto-focus motion correction. Analyses were performed from a central 5×5×5 mm(3) cuboid subvolume and trabecular volume consisting of a 13 mm axial slab encompassing the entire medullary cavity. Whole-volume axial stiffness and sub-regional Youngs and shear moduli were computed by finite-element analysis. Whole-volume-derived aggregate mean coefficient of variation of all structural parameters was 4.4% (range 1.8% to 7.7%) and 4.0% for axial stiffness; corresponding data in the subvolume were 6.5% (range 1.6% to 13.0%) for structural, and 5.5% (range 4.6% to 6.5%) for mechanical parameters. Aggregate ICC was 0.976 (range 0.947 to 0.986) and 0.992 for whole-volume-derived structural parameters and axial stiffness, and 0.946 (range 0.752 to 0.991) and 0.974 (range 0.965 to 0.978) for subvolume-derived structural and mechanical parameters, respectively. The strongest predictors of whole-volume axial stiffness were BV/TV, junction density, skeleton density and Tb.N (R(2) 0.79-0.87). The same parameters were also highly predictive of sub-regional axial modulus (R(2) 0.88-0.91). The data suggest that the method is suited for longitudinal assessment of the response to therapy. The underlying technology is portable and should be compatible with all general-purpose MRI scanners, which is appealing considering the very large installed base of this modality.

Volumetric Cortical Bone Porosity Assessment with MR Imaging: Validation and Clinical Feasibility

PURPOSE To develop a method to assess volumetric cortical bone porosity in clinically practical acquisition times by measuring the signal decay at only two echo times (TEs) as part of a single three-dimensional ultrashort TE (UTE) magnetic resonance (MR) examination. MATERIALS AND METHODS The study was approved by the institutional review board and complied with HIPAA guidelines. Written informed consent was obtained from all subjects. A marker of cortical bone porosity called porosity index was defined as the ratio of UTE image intensities at a long and short TE, and the results were compared with biexponential analysis. Porosity index of midtibia cortical bone samples obtained from 16 donors was compared with ground-truth porosity by using micro-computed tomographic (CT) imaging and bone mineral density by peripheral quantitative CT scanner. Reproducibility of porosity index were tested in volunteers, and clinical feasibility was evaluated in postmenopausal women. Interparameter associations were assessed by using Pearson or Spearman correlation coefficient. RESULTS Bone specimen porosity index was correlated with micro-CT imaging porosity (R(2) = 0.79) and pore size (R(2) = 0.81); age (R(2) = 0.64); peripheral quantitative CT scanner density (R(2) = 0.49, negatively); and pore water fraction (R(2) = 0.62) and T2* (R(2) = 0.64) by biexponential analysis. The reproducibility study yielded a coefficient of variation of 2.2% and intraclass correlation coefficient of 0.97. The study that involved postmenopausal women showed a wide range of porosity index (15%-38%). CONCLUSION A two-point MR imaging method to assess cortical bone porosity in humans was conceived and validated. This approach has the potential for clinical use to assess changes in cortical bone porosity that result from disease or in response to therapy. (©) RSNA, 2015 Online supplemental material is available for this article.