Mechanical implications of estrogen supplementation in early postmenopausal women

More accurate trabecular bone imaging using UTE MRI at the resonance frequency of fat

High-resolution magnetic resonance imaging (HR-MRI) has been increasingly used to assess the trabecular bone structure. High susceptibility at the marrow/bone interface may significantly reduce the marrow's apparent transverse relaxation time (T2*), overestimating trabecular bone thickness. Ultrashort echo time MRI (UTE-MRI) can minimize the signal loss caused by susceptibility-induced T2* shortening. However, UTE-MRI is sensitive to chemical shift artifacts, which manifest as spatial blurring and ringing artifacts partially due to non-Cartesian sampling. In this study, we proposed UTE-MRI at the resonance frequency of fat to minimize marrow-related chemical shift artifacts and the overestimation of trabecular thickness. Cubes of trabecular bone from six donors (75 ± 4 years old) were scanned using a 3 T clinical scanner at the resonance frequencies of fat and water, respectively, using 3D UTE sequences with five TEs (0.032, 1.1, 2.2, 3.3, and 4.4 ms) and a clinical 3D gradient echo (GRE) sequence at 0.2 × 0.2 × 0.4 mm3 voxel size. Trabecular bone thickness was measured in 30 regions of interest (ROIs) per sample. MRI results were compared with thicknesses obtained from micro-computed tomography (μCT) at 50 μm3 voxel size. Linear regression models were used to calculate the coefficient of determination between MRI- and μCT-based trabecular thickness. All MRI-based trabecular thicknesses showed significant correlations with μCT measurements. The correlations were higher (examined with paired Student's t-test, P < 0.01) for 3D UTE images performed at the fat frequency (R2 = 0.59-0.74, P < 0.01) than those at the water frequency (R2 = 0.18-0.52, P < 0.01) and clinical GRE images (R2 = 0.39-0.47, P < 0.01). Significantly reduced correlations were observed with longer TEs. This study highlighted the feasibility of UTE-MRI at the fat frequency for a more accurate assessment of trabecular bone thickness.

Six-minute, in vivo MRI quantification of proximal femur trabecular bone 3D microstructure

BACKGROUND

Assessment of proximal femur trabecular bone microstructure in vivo by magnetic resonance imaging has recently been validated for acquiring information independent of bone mineral density in osteoporotic patients. However, the requisite signal-to-noise ratio (SNR) and resolution for interrogation of the trabecular microstructure at this anatomical location prolongs the scan duration and renders the imaging protocol clinically infeasible. Parallel imaging and compressed sensing (PICS) techniques can reduce the scan duration of the imaging protocol without substantially compromising image quality. The present work investigates the limits of acceleration for a commonly used PICS technique, ℓ1-ESPIRiT, for the purpose of quantifying measures of trabecular bone microarchitecture. Based on a desired error tolerance, a six-minute, prospectively accelerated variant of the imaging protocol was developed and assessed for intersession reproducibility and agreement with the longer reference scan.

PURPOSE

To investigate the limits of acceleration for MRI-based trabecular bone quantification by parallel imaging and compressed sensing reconstruction, and to develop a prototypical imaging protocol for assessing the proximal femur microstructure in a clinically practical scan time.

METHODS

Healthy participants (n = 11) were scanned by a 3D balanced steady-state free precession (bSSFP) sequence satisfying the Nyquist criterion with a scan duration of about 18 min. The raw data were retrospectively undersampled and reconstructed to mimic various acceleration factors ranging from 2 to 6. Trabecular volumes-of-interest in four major femoral regions (greater trochanter, intertrochanteric region, femoral neck, and femoral head) were analyzed and six relevant measures of trabecular bone microarchitecture (bone volume fraction, surface-to-curve ratio, erosion index, elastic modulus, trabecular thickness, plates-to-rods ratio) were obtained for images of all accelerations. To assess agreement, median percent error and intraclass correlation coefficients (ICCs) were computed using the fully-sampled data as reference. Based on this analysis, a prospectively 3-fold accelerated sequence with a duration of about 6 min was developed and the analysis was repeated.

RESULTS

A prospective acceleration factor of 3 demonstrated comparable performance in reproducibility and absolute agreement to the fully-sampled scan. The median CoV over all image-derived metrics was generally <6 % and ICCs >0.70. Also, measurements from prospectively 3-fold accelerated scans demonstrated in general median percent errors of <7 % and ICCs >0.70.

CONCLUSION

The present work proposes a method to make in vivo quantitative assessment of proximal femur trabecular microstructure with a clinically practical scan duration of about 6 min.

Exploring the Possibility of Measuring Vertebrae Bone Structure Metrics Using MDCT Images: An Unpaired Image-to-Image Translation Method

Bone structure metrics are vital for the evaluation of vertebral bone strength. However, the gold standard for measuring bone structure metrics, micro-Computed Tomography (micro-CT), cannot be used in vivo, which hinders the early diagnosis of fragility fractures. This paper used an unpaired image-to-image translation method to capture the mapping between clinical multidetector computed tomography (MDCT) and micro-CT images and then generated micro-CT-like images to measure bone structure metrics. MDCT and micro-CT images were scanned from 75 human lumbar spine specimens and formed training and testing sets. The generator in the model focused on learning both the structure and detailed pattern of bone trabeculae and generating micro-CT-like images, and the discriminator determined whether the generated images were micro-CT images or not. Based on similarity metrics (i.e., SSIM and FID) and bone structure metrics (i.e., bone volume fraction, trabecular separation and trabecular thickness), a set of comparisons were performed. The results show that the proposed method can perform better in terms of both similarity metrics and bone structure metrics and the improvement is statistically significant. In particular, we compared the proposed method with the paired image-to-image method and analyzed the pros and cons of the method used.

CT ‐based Stiffness Measures of Trabecular Bone Microstructure — Cadaveric Validation and In Vivo Application

Osteoporosis causes bone fragility and elevates fracture risk. Applications of finite element (FE) analysis (FEA) for assessment of trabecular bone (Tb) microstructural strength at whole‐body computed tomography (CT) imaging are limited due to challenges with Tb microstructural segmentation. We present a nonlinear FEA method for distal tibia CT scans evading binary segmentation of Tb microstructure, while accounting for bone microstructural distribution. First, the tibial axis in a CT scan was aligned with the FE loading axis. FE cubic mesh elements were modeled using image voxels, and CT intensity values were calibrated to ash density defining mechanical properties at individual elements. For FEA of an upright volume of interest (VOI), the bottom surface was fixed, and a constant displacement was applied at each vertex on the top surface simulating different loading conditions. The method was implemented and optimized using the ANSYS software. CT‐derived computational modulus values were repeat scan reproducible (intraclass correlation coefficient [ICC] ≥ 0.97) and highly correlated (r ≥ 0.86) with the micro‐CT (μCT)‐derived values. FEA‐derived von Mises stresses over the segmented Tb microregion were significantly higher (p < 1 × 10⁻¹¹) than that over the marrow space. In vivo results showed that both shear and compressive modulus for males were higher (p < 0.01) than for females. Effect sizes for different modulus measures between males and females were moderate‐to‐high (≥0.55) and reduced to small‐to‐negligible (<0.40) when adjusted for pure lean mass. Among body size and composition attributes, pure lean mass and height showed highest (r ∈ [0.45 0.56]) and lowest (r ∈ [0.25 0.39]) linear correlation, respectively, with FE‐derived modulus measures. In summary, CT‐based nonlinear FEA provides an effective surrogate measure of Tb microstructural stiffness, and the relaxation of binary segmentation will extend the scope for FEA in human studies using in vivo imaging at relatively low‐resolution. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Finite element analysis of trabecular bone microstructure using CT imaging and continuum mechanical modelling

Purpose

Osteoporosis is a bone disease associated with enhanced bone loss, microstructural degeneration, and fracture‐risk. Finite element (FE) modeling is used to estimate trabecular bone (Tb) modulus from high‐resolution three‐dimensional (3‐D) imaging modalities including micro‐computed tomography (CT), magnetic resonance imaging (MRI), and high‐resolution peripheral quantitative CT (HR‐pQCT). This paper validates an application of voxel‐based continuum finite element analysis (FEA) to predict Tb modulus from clinical CT imaging under a condition similar to in vivo imaging by comparing with measures derived by micro‐CT and experimental approaches.

Method

Voxel‐based continuum FEA methods for CT imaging were implemented using linear and nonlinear models and applied on distal tibial scans under a condition similar to in vivo imaging. First, tibial axis in a CT scan was aligned with the coordinate z‐axis at 150 μm isotropic voxels. FEA was applied on an upright cylindrical volume of interests (VOI) with its axis coinciding with the tibial bone axis. Voxel volume, edge, and vertex elements and their connectivity were defined as per the isotropic image grid. A calibration phantom was used to calibrate CT numbers in Hounsfield unit to bone mineral density (BMD) values, which was then converted into calcium hydroxyapatite (CHA) density. Mechanical properties at each voxel volume element was defined using its ash‐density defined on CT‐derived CHA density. For FEA, the bottom surface of the cylindrical VOI was fixed and a constant displacement was applied along the z‐direction at each vertex element on the top surface to simulate a physical axial compressive loading condition. Finally, a Poisson's ratio of 0.3 was applied, and Tb modulus (MPa) was computed as the ratio of average von Mises stress (MPa) of volume elements on the top surface and the applied displacement. FEA parameters including mesh element size, substep number, and different tolerance values were optimized.

Results

CT‐derived Tb modulus values using continuum FEA showed high linear correlation with the micro‐CT‐derived reference values (r ∈ [0.87 0.90]) as well as experimentally measured values (r ∈ [0.80 0.87]). Linear correlation of computed modulus with their reference values using continuum FEA with linear modeling was comparable with that obtained by nonlinear modeling. Nonlinear continuum FEA‐based modulus values (mean of 1087.2 MPa) showed greater difference from their reference values (mean of 1498.9 MPa using micro‐CT‐based FEA) as compared with linear continuum methods. High repeat CT scan reproducibility (intra‐class correlation [ICC] = 0.98) was observed for computed modulus values using both linear and nonlinear continuum FEA. It was observed that high stress regions coincide with Tb microstructure as fuzzily characterized by BMD values. Distributions of von Mises stress over Tb microstructure and marrow regions were significantly different (p < 10^–8).

Conclusion

Voxel‐based continuum FEA offers surrogate measures of Tb modulus from CT imaging under a condition similar to in vivo imaging that alleviates the need for segmentation of Tb and marrow regions, while accounting for bone distribution at the microstructural level. This relaxation of binary segmentation will extend the scope of FEA application to assess mechanical properties of bone microstructure at relatively low‐resolution imaging.

Effect of Low-Intensity Vibration on Bone Strength, Microstructure, and Adiposity in Pre-Osteoporotic Postmenopausal Women: A Randomized Placebo-Controlled Trial

There has been evidence that cyclical mechanical stimulation may be osteogenic, thus providing opportunities for nonpharmacological treatment of degenerative bone disease. Here, we applied this technology to a cohort of postmenopausal women with varying bone mineral density (BMD) T-scores at the total hip (-0.524 ± 0.843) and spine (-0.795 ± 1.03) to examine the response to intervention after 1 year of daily treatment with 10 minutes of vibration therapy in a randomized double-blinded trial. The device operates either in an active mode (30 Hz and 0.3 g) or placebo. Primary endpoints were changes in bone stiffness at the distal tibia and marrow adiposity of the vertebrae, based on 3 Tesla high-resolution MRI and spectroscopic imaging, respectively. Secondary outcome variables included distal tibial trabecular microstructural parameters and vertebral deformity determined by MRI, volumetric and areal bone densities derived using peripheral quantitative computed tomography (pQCT) of the tibia, and dual-energy X-ray absorptiometry (DXA)-based BMD of the hip and spine. Device adherence was 83% in the active group (n = 42) and 86% in the placebo group (n = 38) and did not differ between groups (p = .7). The mean 12-month changes in tibial stiffness in the treatment group and placebo group were +1.31 ± 6.05% and -2.55 ± 3.90%, respectively (group difference 3.86%, p = .0096). In the active group, marrow fat fraction significantly decreased after 12 months of intervention (p = .0003), whereas no significant change was observed in the placebo group (p = .7; group difference -1.59%, p = .029). Mean differences of the changes in trabecular bone volume fraction (p = .048) and erosion index (p = .044) were also significant, as was pQCT-derived trabecular volumetric BMD (vBMD; p = .016) at the tibia. The data are commensurate with the hypothesis that vibration therapy is protective against loss in mechanical strength and, further, that the intervention minimizes the shift from the osteoblastic to the adipocytic lineage of mesenchymal stem cells. © 2020 American Society for Bone and Mineral Research (ASBMR).

High resolution imaging in bone tissue research-review

This review article focuses on imaging of bone tissue to understand skeletal health with regards to bone quality. Skeletal fragility fractures are due to bone diseases such as osteoporosis which result in low bone mass and bone mineral density (BMD) leading to high risk of fragility fractures. Recent advances in imaging and analysis technologies have highly benefitted the field of biological sciences. In particular, their application in skeletal health has been of significant importance in understanding bone mechanical behavior (structure and properties) at the tissue level. While synchrotron based microCT technique has remained the gold standard for non-destructive evaluation of structure in material and biological sciences, several lab based microCT systems have been developed to provide high resolution imaging of specimens with greater access, and ease of use in laboratory settings. Lab based microCT scanners are widely used in the bone field as a standard tool to evaluate three-dimensional (3D) morphologies of bone structure at image resolutions appropriate for bone samples from small animals to bone biopsy specimens from humans. Both synchrotron and standard lab based microCT systems provide high resolution imaging ex vivo for a small sized specimen. A few X-ray based systems are also commercially available for in vivo scanning at relatively low image resolutions. Synchrotron-based CT microscopy is being used for various ultra-high-resolution image analyses using complex 3D software. However, the synchrotron-based CT technology is in high demand, allows only limited numbers of specimens, expensive, requires complex additional instrumentation, and is not easily available to researchers as it requires access to a synchrotron source which is always limited. Therefore, desktop laboratory scanners (microXCT, Zeiss/Xradia, Scanco, SkyScan. etc.), mimicking the synchrotron based CT technology or image resolution, have been developed to solve the accessibility issues. These lab based scanners have helped both material science, and the bone field to investigate bone tissue morphologies at submicron mage resolutions. Considerable progress has been made in both in vivo and ex vivo imaging towards providing high resolution images of bone tissue. Both clinical and research imaging technologies will continue to improve and help understand osteoporosis and other related skeletal issues in order to develop targeted treatments for bone fragility. This review summarizes the high resolution imaging work in bone research.

Quantitative Magnetic Resonance Imaging of Cortical and Trabecular Bone

Bone is a composite material consisting of mineral, organic matrix, and water. Water in bone can be categorized as bound water (BW), which is bound to bone mineral and organic matrix, or as pore water (PW), which resides in Haversian canals as well as in lacunae and canaliculi. Bone is generally classified into two types: cortical bone and trabecular bone. Cortical bone is much denser than trabecular bone that is surrounded by marrow and fat. Magnetic resonance (MR) imaging has been increasingly used for noninvasive assessment of both cortical bone and trabecular bone. Bone typically appears as a signal void with conventional MR sequences because of its short T2*. Ultrashort echo time (UTE) sequences with echo times 100 to 1,000 times shorter than those of conventional sequences allow direct imaging of BW and PW in bone. This article summarizes several quantitative MR techniques recently developed for bone evaluation. Specifically, we discuss the use of UTE and adiabatic inversion recovery prepared UTE sequences to quantify BW and PW, UTE magnetization transfer sequences to quantify collagen backbone protons, UTE quantitative susceptibility mapping sequences to assess bone mineral, and conventional sequences for high-resolution imaging of PW as well as the evaluation of trabecular bone architecture.

Developmental low-dose exposure to bisphenol A induces chronic inflammation, bone marrow fibrosis and reduces bone stiffness in female rat offspring only

BACKGROUND

Developmental exposure to low doses of the endocrine disruptor bisphenol A (BPA) is known to alter bone tissue in young rodents, although how bone tissue is affected in aged animals is not well known. We have recently shown that low-dose developmental exposure to BPA increases procollagen type I N-terminal propeptide (P1NP) levels, a peptide formed during type 1 collagen synthesis, in plasma of 5-week-old female rat offspring while male offspring showed reduced bone size.

OBJECTIVE

To analyze offspring bone phenotype at 52 weeks of age and clarify whether the BPA-induced increase in P1NP levels at 5 weeks is an early sign of bone marrow fibrosis development.

METHODS

As in our 5-week study, pregnant Fischer 344 rats were exposed to BPA via drinking water corresponding to 0.5 μg/kg BW/day (BPA0.5), which is in the range of human daily exposure, or 50 μg/kg BW/day (BPA50) from gestational day 3.5 until postnatal day 22. Controls were given only vehicle. The offspring were sacrificed at 52 weeks of age. Bone effects were analyzed using peripheral quantitative and micro-computed tomography (microCT), 3-point bending test, plasma markers and histological examination.

RESULTS

Compared to a smaller bone size at 5 weeks, at the age of 52 weeks, femur size in male offspring had been normalized in developmentally BPA-exposed rats. The 52-week-old female offspring showed, like the 5-week-old siblings, higher plasma P1NP levels compared to controls but no general increasing bone growth or strength. However, 2 out of 14 BPA-exposed female offspring bones developed extremely thick cortices later in life, discovered by systematic in vivo microCT scanning during the study. This was not observed in male offspring or in female controls. Biomechanical testing revealed that both doses of developmental BPA exposure reduced femur stiffness only in female offspring. In addition, histological analysis showed an increased number of fibrotic lesions only in the bone marrow of female rat offspring developmentally exposed to BPA. In line with this, plasma markers of inflammation, Tnf (in BPA0.5) and Timp1 (in BPA50) were increased exclusively in female offspring.

CONCLUSIONS

Developmental BPA exposure at an environmentally relevant concentration resulted in female-specific effects on bone as well as on plasma biomarkers of collagen synthesis and inflammation. Even a dose approximately eight times lower than the current temporary EFSA human tolerable daily intake of 4 μg/kg BW/day, appeared to induce bone stiffness reduction, bone marrow fibrosis and chronic inflammation in female rat offspring later in life.

Micro-Finite Element Analysis of the Proximal Femur on the Basis of High-Resolution Magnetic Resonance Images

PURPOSE OF REVIEW

Hip fractures have catastrophic consequences. The purpose of this article is to review recent developments in high-resolution magnetic resonance imaging (MRI)-guided finite element analysis (FEA) of the hip as a means to determine subject-specific bone strength.

RECENT FINDINGS

Despite the ability of DXA to predict hip fracture, the majority of fractures occur in patients who do not have BMD T scores less than - 2.5. Therefore, without other detection methods, these individuals go undetected and untreated. Of methods available to image the hip, MRI is currently the only one capable of depicting bone microstructure in vivo. Availability of microstructural MRI allows generation of patient-specific micro-finite element models that can be used to simulate real-life loading conditions and determine bone strength. MRI-based FEA enables radiation-free approach to assess hip fracture strength. With further validation, this technique could become a potential clinical tool in managing hip fracture risk.

Artificial Intelligence Applied to Osteoporosis: A Performance Comparison of Machine Learning Algorithms in Predicting Fragility Fractures From MRI Data

BACKGROUND

A current challenge in osteoporosis is identifying patients at risk of bone fracture.

PURPOSE

To identify the machine learning classifiers that predict best osteoporotic bone fractures and, from the data, to highlight the imaging features and the anatomical regions that contribute most to prediction performance.

STUDY TYPE

Prospective (cross-sectional) case-control study.

POPULATION

Thirty-two women with prior fragility bone fractures, of mean age = 61.6 and body mass index (BMI) = 22.7 kg/m² , and 60 women without fractures, of mean age = 62.3 and BMI = 21.4 kg/m² . Field Strength/ Sequence: 3D FLASH at 3T.

ASSESSMENT

Quantitative MRI outcomes by software algorithms. Mechanical and topological microstructural parameters of the trabecular bone were calculated for five femoral regions, and added to the vector of features together with bone mineral density measurement, fracture risk assessment tool (FRAX) score, and personal characteristics such as age, weight, and height. We fitted 15 classifiers using 200 randomized cross-validation datasets. Statistical Tests: Data: Kolmogorov-Smirnov test for normality. Model Performance: sensitivity, specificity, precision, accuracy, F1-test, receiver operating characteristic curve (ROC). Two-sided t-test, with P < 0.05 for statistical significance.

RESULTS

The top three performing classifiers are RUS-boosted trees (in particular, performing best with head data, F1 = 0.64 ± 0.03), the logistic regression and the linear discriminant (both best with trochanteric datasets, F1 = 0.65 ± 0.03 and F1 = 0.67 ± 0.03, respectively). A permutation of these classifiers comprised the best three performers for four out of five anatomical datasets. After averaging across all the anatomical datasets, the score for the best performer, the boosted trees, was F1 = 0.63 ± 0.03 for All-features dataset, F1 = 0.52 ± 0.05 for the no-MRI dataset, and F1 = 0.48 ± 0.06 for the no-FRAX dataset. Data Conclusion: Of many classifiers, the RUS-boosted trees, the logistic regression, and the linear discriminant are best for predicting osteoporotic fracture. Both MRI and FRAX independently add value in identifying osteoporotic fractures. The femoral head, greater trochanter, and inter-trochanter anatomical regions within the proximal femur yielded better F1-scores for the best three classifiers.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:1029-1038.

The reproducibility of measuring trabecular bone parameters using a commercially available high-resolution magnetic resonance imaging approach: A pilot study

Bone imaging is currently the best non-invasive way to assess changes to bone associated with aging or chronic disease. However, common imaging techniques such as dual energy x-ray absorptiometry are associated with limitations. Magnetic resonance imaging (MRI) is a radiation-free technique that can measure bone microarchitecture. However, published MRI bone assessment protocols use specialized MRI coils and sequences and therefore have limited transferability across institutions. We developed a protocol on a Siemens 3 Tesla MRI machine, using a commercially available coil (Siemens 15 CH knee coil), and manufacturer supplied sequences to acquire images at the tibia. We tested the reproducibility of the FSE and the GE Axial sequences and hypothesized that both would generate reproducible trabecular bone parameters. Eight healthy adults (age 25.5 ± 5.4 years) completed three measurements of each MRI sequence at the tibia. Each of the images was processed for 8 different bone parameters (such as volumetric bone volume fraction). We computed the coefficient of variation (CV) and intraclass correlation coefficients (ICC) to assess reproducibility and reliability. Both sequences resulted in trabecular parameters that were reproducible (CV <5% for most) and reliable (ICC >80% for all). Our study is one of the first to report that a commercially available MRI protocol can result in reproducible data, and is significant as MRI may be an accessible method to measure bone microarchitecture in clinical or research environments. This technique requires further testing, including validation and evaluation in other populations.

•

Trabecular bone is difficult to measure using commercial MRI techniques

•

Reproducibility of a MRI protocol measuring trabecular bone was assessed

•

Tibia trabecular bone was reproducible using a knee coil and a FSE Axial sequence

•

Tibia trabecular bone was reproducible using a knee coil and a GE Axial sequence

Accuracy of MRI-based finite element assessment of distal tibia compared to mechanical testing

High-resolution MRI-derived finite element analysis (FEA) has been used in translational research to estimate the mechanical competence of human bone. However, this method has yet to be validated adequately under in vivo imaging spatial resolution or signal-to-noise conditions. We therefore compared MRI-based metrics of bone strength to those obtained from direct, mechanical testing. The study was conducted on tibiae from 17 human donors (12 males and five females, aged 33 to 88years) with no medical history of conditions affecting bone mineral homeostasis. A 25mm segment from each distal tibia underwent MR imaging in a clinical 3-Tesla scanner using a fast large-angle spin-echo (FLASE) sequence at 0.137mm×0.137mm×0.410mm voxel size, in accordance with in vivo scanning protocol. The resulting high-resolution MR images were processed and used to generate bone volume fraction maps, which served as input for the micro-level FEA model. Simulated compression was applied to compute stiffness, yield strength, ultimate strength, modulus of resilience, and toughness, which were then compared to metrics obtained from mechanical testing. Moderate to strong positive correlations were found between computationally and experimentally derived values of stiffness (R²=0.77, p<0.0001), yield strength (R²=0.38, p=0.0082), ultimate strength (R²=0.40, p=0.0067), and resilience (R²=0.46, p=0.0026), but only a weak, albeit significant, correlation was found for toughness (R²=0.26, p=0.036). Furthermore, experimentally derived yield strength and ultimate strength were moderately correlated with MRI-derived stiffness (R²=0.48, p=0.0022 and R²=0.58, p=0.0004, respectively). These results suggest that high-resolution MRI-based finite element (FE) models are effective in assessing mechanical parameters of distal skeletal extremities.

Advances in imaging approaches to fracture risk evaluation

Fragility fractures are a growing problem worldwide, and current methods for diagnosing osteoporosis do not always identify individuals who require treatment to prevent a fracture and may misidentify those not a risk. Traditionally, fracture risk is assessed using dual-energy X-ray absorptiometry, which provides measurements of areal bone mineral density at sites prone to fracture. Recent advances in imaging show promise in adding new information that could improve the prediction of fracture risk in the clinic. As reviewed herein, advances in quantitative computed tomography (QCT) predict hip and vertebral body strength; high-resolution HR-peripheral QCT (HR-pQCT) and micromagnetic resonance imaging assess the microarchitecture of trabecular bone; quantitative ultrasound measures the modulus or tissue stiffness of cortical bone; and quantitative ultrashort echo-time MRI methods quantify the concentrations of bound water and pore water in cortical bone, which reflect a variety of mechanical properties of bone. Each of these technologies provides unique characteristics of bone and may improve fracture risk diagnoses and reduce prevalence of fractures by helping to guide treatment decisions.

Clinical Evaluation of Bone Strength and Fracture Risk

PURPOSE OF REVIEW

This paper seeks to evaluate and compare recent advances in the clinical assessment of the changes in bone mechanical properties that take place as a result of osteoporosis and other metabolic bone diseases and their treatments.

RECENT FINDINGS

In addition to the standard of DXA-based areal bone mineral density (aBMD), a variety of methods, including imaging-based structural measurements, finite element analysis (FEA)-based techniques, and alternate methods including ultrasound, bone biopsy, reference point indentation, and statistical shape and density modeling, have been developed which allow for reliable prediction of bone strength and fracture risk. These methods have also shown promise in the evaluation of treatment-induced changes in bone mechanical properties. Continued technological advances allowing for increasingly high-resolution imaging with low radiation dose, together with the expanding adoption of DXA-based predictions of bone structure and mechanics, as well as the increasing awareness of the importance of bone material properties in determining whole-bone mechanics, lead us to anticipate substantial future advances in this field.

An update on bone imaging and markers in chronic kidney disease

Bone disorders in chronic kidney disease (CKD) are associated with heightened risks of fractures, vascular calcification, poor quality of life and mortality compared to the general population. However, diagnosis and management of these disorders in CKD are complex and appreciably limited by current diagnostic modalities. Areas covered: Bone histomorphometry remains the gold standard for diagnosis but is not widely utilised and lacks feasibility as a monitoring tool. In practice, non-invasive imaging and biochemical markers are preferred to guide therapeutic decisions. Expert commentary: This review aims to summarize the risk factors for, and spectrum of bone disease in CKD, as well as appraise the clinical utility of dual energy X-ray densitometry, peripheral quantitative computed tomography, high-resolution peripheral quantitative computed tomography, and bone turnover markers.

Bone quality: the determinants of bone strength and fragility

Bone fragility is a major health concern, as the increased risk of bone fractures has devastating outcomes in terms of mortality, decreased autonomy, and healthcare costs. Efforts made to address this problem have considerably increased our knowledge about the mechanisms that regulate bone formation and resorption. In particular, we now have a much better understanding of the cellular events that are triggered when bones are mechanically stimulated and how these events can lead to improvements in bone mass. Despite these findings at the molecular level, most exercise intervention studies reveal either no effects or only minor benefits of exercise programs in improving bone mineral density (BMD) in osteoporotic patients. Nevertheless, and despite that BMD is the gold standard for diagnosing osteoporosis, this measure is only able to provide insights regarding the quantity of bone tissue. In this article, we review the complex structure of bone tissue and highlight the concept that its mechanical strength stems from the interaction of several different features. We revisited the available data showing that bone mineralization degree, hydroxyapatite crystal size and heterogeneity, collagen properties, osteocyte density, trabecular and cortical microarchitecture, as well as whole bone geometry, are determinants of bone strength and that each one of these properties may independently contribute to the increased or decreased risk of fracture, even without meaningful changes in aBMD. Based on these findings, we emphasize that while osteoporosis (almost) always causes bone fragility, bone fragility is not always caused just by osteoporosis, as other important variables also play a major role in this etiology. Furthermore, the results of several studies showing compelling data that physical exercise has the potential to improve bone quality and to decrease fracture risk by influencing each one of these determinants are also reviewed. These findings have meaningful clinical repercussions as they emphasize the fact that, even without leading to improvements in BMD, exercise interventions in patients with osteoporosis may be beneficial by improving other determinants of bone strength.

Sustained delivery of rhBMP-2 by means of poly(lactic-co-glycolic acid) microspheres: cranial bone regeneration without heterotopic ossification or craniosynostosis

BACKGROUND

Commercially available recombinant human bone morphogenetic protein 2 (rhBMP2) has demonstrated efficacy in bone regeneration, but not without significant side effects. The authors used rhBMP2 encapsulated in poly(lactic-co-glycolic acid) (PLGA) microspheres placed in a rabbit cranial defect model to test whether low-dose, sustained delivery can effectively induce bone regeneration.

METHODS

The rhBMP2 was encapsulated in 15% PLGA using a double-emulsion, solvent extraction/evaporation technique, and its release kinetics and bioactivity were tested. Two critical-size defects (10 mm) were created in the calvaria of New Zealand white rabbits (5 to 7 months of age, male and female) and filled with a collagen scaffold containing either (1) no implant, (2) collagen scaffold only, (3) PLGA-rhBMP2 (0.1 μg per implant), or (4) free rhBMP2 (0.1 μg per implant). After 6 weeks, the rabbits were killed and defects were analyzed by micro-computed tomography, histology, and finite element analysis.

RESULTS

The rhBMP2 delivered by means of bioactive PLGA microspheres resulted in higher volumes and surface area coverage of new bone than an equal dose of free rhBMP2 by micro-computed tomography (p=0.025 and p=0.025). Finite element analysis indicated that the mechanical competence using the regional elastic modulus did not differ with rhBMP2 exposure (p=0.70). PLGA-rhBMP2 did not demonstrate heterotopic ossification, craniosynostosis, or seroma formation.

CONCLUSIONS

Sustained delivery by means of PLGA microspheres can significantly reduce the rhBMP2 dose required for de novo bone formation. Optimization of the delivery system may be a key to reducing the risk for recently reported rhBMP2-related adverse effects.

Novel assessment tools for osteoporosis diagnosis and treatment

This review describes new technologies for the diagnosis and treatment, including fracture risk prediction, of postmenopausal osteoporosis. Four promising technologies and their potential for clinical translation and basic science studies are discussed. These include reference point indentation (RPI), Fourier transform infrared spectroscopy (FTIR), Raman spectroscopy, and magnetic resonance imaging (MRI). While each modality exploits different physical principles, the commonality is that none of them require use of ionizing radiation. To provide context for the new developments, brief summaries are provided for the current state of biomarker assays, fracture risk assessment (FRAX), and other fracture risk prediction algorithms and quantitative ultrasound (QUS) measurements.

Registration-based autofocusing technique for automatic correction of motion artifacts in time-series studies of high-resolution bone MRI

PURPOSE

To develop a registration-based autofocusing (RAF) motion correction technique for high-resolution trabecular bone (TB) imaging and to evaluate its performance on in vivo MR data.

MATERIALS AND METHODS

The technique combines serial registration with a previously developed motion correction technique - autofocusing - for automatic correction of subject movement degradation of MR images acquired in longitudinal studies. The method was tested on in vivo images of the distal radius to measure improvements in serial reproducibility of parameters in 12 women (ages 50-75 years), and to compare with the navigator echo-based correction and autofocusing. Furthermore, the technique's ability to optimize the sensitivity to detect simulated bone loss was ascertained.

RESULTS

The new technique yielded superior reproducibility of image-derived structural and mechanical parameters. Average coefficient of variation across all parameters improved by 12.5%, 27.0%, 33.5%, and 37.0%, respectively, following correction by navigator echoes, autofocusing, and the RAF technique (without and with correction for rotational motion); average intra-class correlation coefficient increased by 1.2%, 2.2%, 2.8%, and 3.2%, respectively. Furthermore, simulated bone loss (5%) was well recovered independent of the choice of reference image (4.71% or 4.86% with respect to using either the original or the image subjected to bone loss) in the time series.

CONCLUSION

The data suggest that our technique simultaneously corrects for intra-scan motion corruption while improving inter-scan registration. Furthermore, the technique is not biased by small changes in bone architecture between time-points.

Finite element analysis applied to 3-T MR imaging of proximal femur microarchitecture: lower bone strength in patients with fragility fractures compared with control subjects

PURPOSE

To determine the feasibility of using finite element analysis applied to 3-T magnetic resonance (MR) images of proximal femur microarchitecture for detection of lower bone strength in subjects with fragility fractures compared with control subjects without fractures.

MATERIALS AND METHODS

This prospective study was institutional review board approved and HIPAA compliant. Written informed consent was obtained. Postmenopausal women with (n = 22) and without (n = 22) fragility fractures were matched for age and body mass index. All subjects underwent standard dual-energy x-ray absorptiometry. Images of proximal femur microarchitecture were obtained by using a high-spatial-resolution three-dimensional fast low-angle shot sequence at 3 T. Finite element analysis was applied to compute elastic modulus as a measure of strength in the femoral head and neck, Ward triangle, greater trochanter, and intertrochanteric region. The Mann-Whitney test was used to compare bone mineral density T scores and elastic moduli between the groups. The relationship (R(2)) between elastic moduli and bone mineral density T scores was assessed.

RESULTS

Patients with fractures showed lower elastic modulus than did control subjects in all proximal femur regions (femoral head, 8.51-8.73 GPa vs 9.32-9.67 GPa; P = .04; femoral neck, 3.11-3.72 GPa vs 4.39-4.82 GPa; P = .04; Ward triangle, 1.85-2.21 GPa vs 3.98-4.13 GPa; P = .04; intertrochanteric region, 1.62-2.18 GPa vs 3.86-4.47 GPa; P = .006-.007; greater trochanter, 0.65-1.21 GPa vs 1.96-2.62 GPa; P = .01-.02), but no differences in bone mineral density T scores. There were weak relationships between elastic moduli and bone mineral density T scores in patients with fractures (R(2) = 0.25-0.31, P = .02-.04), but not in control subjects. CONCLUSION Finite element analysis applied to high-spatial-resolution 3-T MR images of proximal femur microarchitecture can allow detection of lower elastic modulus, a marker of bone strength, in subjects with fragility fractures compared with control subjects. MR assessment of proximal femur strength may provide information about bone quality that is not provided by dual-energy x-ray absorptiometry.

Effects of testosterone and growth hormone on the structural and mechanical properties of bone by micro-MRI in the distal tibia of men with hypopituitarism

CONTEXT

Severe deficiencies of testosterone (T) and GH are associated with low bone mineral density (BMD) and increased fracture risk. Replacement of T in hypogonadal men improves several bone parameters. Replacement of GH in GH-deficient men improves BMD.

OBJECTIVE

Our objective was to determine whether T and GH treatment together improves the structural and mechanical parameters of bone more than T alone in men with hypopituitarism.

DESIGN AND SUBJECTS

This randomized, prospective, 2-year study included 32 men with severe deficiencies of T and GH due to panhypopituitarism.

INTERVENTION

Subjects were randomized to receive T alone (n = 15) or T and GH (n = 17) for 2 years.

MAIN OUTCOME MEASURES

We evaluated magnetic resonance microimaging-derived structural (bone volume fraction [BVF] and trabecular thickness) and mechanical (axial stiffness [AS], a measure of bone strength) properties of the distal tibia at baseline and after 1 and 2 years of treatment.

RESULTS

Treatment with T and GH did not affect BVF, thickness, or AS differently from T alone. T treatment in all subjects for 2 years increased trabecular BVF by 9.6% (P < .0001), trabecular thickness by 2.6% (P < .001), and trabecular AS by 9.8% (P < .001). In contrast, testosterone treatment in all subjects significantly increased cortical thickness by 2.4% (P < .01) but decreased cortical BVF by -4.7% (P < .01) and cortical AS by -6.9% (P < .01).

CONCLUSION

Combined T and GH treatment of men with hypopituitarism for 2 years did not improve the measured structural or mechanical parameters of the distal tibia more than T alone. However, testosterone significantly increased the structural and mechanical properties of trabecular bone but decreased most of these properties of cortical bone, illustrating the potential importance of assessing trabecular and cortical bone separately in future studies of the effect of testosterone on bone.

Potential of in vivo MRI-based nonlinear finite-element analysis for the assessment of trabecular bone post-yield properties

PURPOSE

Bone strength is the key factor impacting fracture risk. Assessment of bone strength from high-resolution (HR) images have largely relied on linear micro-finite element analysis (μFEA) even though failure always occurs beyond the yield point, which is outside the linear regime. Nonlinear μFEA may therefore be more informative in predicting failure behavior. However, existing nonlinear models applied to trabecular bone (TB) have largely been confined to micro-computed tomography (μCT) and, more recently, HR peripheral quantitative computed tomography (HR-pQCT) images, and typically have ignored evaluation of the post-yield behavior. The primary purpose of this work was threefold: (1) to provide an improved algorithm and program to assess TB yield as well as post-yield properties; (2) to explore the potential benefits of nonlinear μFEA beyond its linear counterpart; and (3) to assess the feasibility and practicality of performing nonlinear analysis on desktop computers on the basis of micro-magnetic resonance (μMR) images obtained in vivo in patients.

METHODS

A method for nonlinear μFE modeling of TB yield as well as post-yield behavior has been designed where material nonlinearity is captured by adjusting the tissue modulus iteratively according to the tissue-level effective strain obtained from linear analysis using a computationally optimized algorithm. The software allows for images at in vivo μMRI resolution as input with retention of grayscale information. Associations between axial stiffness estimated from linear analysis and yield as well as post-yield parameters from nonlinear analysis were investigated from in vivo μMR images of the distal tibia (N = 20; ages: 58-84) and radius (N = 20; ages: 50-75).

RESULTS

All simulations were completed in 1 h or less for 61 strain levels using a desktop computer (dual quad-core Xeon 3.16 GHz CPUs equipped with 40 GB of RAM). Although yield stress and ultimate stress correlated strongly (R(2) > 0.95, p < 0.001) with axial stiffness, toughness correlated moderately at the distal tibia (R(2) = 0.81, p < 0.001) and only weakly at the distal radius (R(2) = 0.34, p = 0.007). Further, toughness was found to vary by up to 16% for bone of very similar axial stiffness (<2%).

CONCLUSIONS

The work demonstrates the practicality of nonlinear μFE simulations at in vivo μMRI resolution, as well as its potential for providing additional information beyond that obtainable from linear analysis. The data suggest that a direct assessment of toughness may provide information not captured by stiffness.

Assessment of trabecular bone yield and post-yield behavior from high-resolution MRI-based nonlinear finite element analysis at the distal radius of premenopausal and postmenopausal women susceptible to osteoporosis

RATIONALE AND OBJECTIVES

To assess the performance of a nonlinear microfinite element model on predicting trabecular bone yield and post-yield behavior based on high-resolution in vivo magnetic resonance images via the serial reproducibility.

MATERIALS AND METHODS

The nonlinear model captures material nonlinearity by iteratively adjusting tissue-level modulus based on tissue-level effective strain. It enables simulations of trabecular bone yield and post-yield behavior from micro magnetic resonance images at in vivo resolution by solving a series of nonlinear systems via an iterative algorithm on a desktop computer. Measures of mechanical competence (yield strain/strength, ultimate strain/strength, modulus of resilience, and toughness) were estimated at the distal radius of premenopausal and postmenopausal women (N = 20, age range 50-75) in whom osteoporotic fractures typically occur. Each subject underwent three scans (20.2 ± 14.5 days). Serial reproducibility was evaluated via coefficient of variation (CV) and intraclass correlation coefficient (ICC).

RESULTS

Nonlinear simulations were completed in an average of 14 minutes per three-dimensional image data set involving analysis of 61 strain levels. The predicted yield strain/strength, ultimate strain/strength, modulus of resilience, and toughness had a mean value of 0.78%, 3.09 MPa, 1.35%, 3.48 MPa, 14.30 kPa, and 32.66 kPa, respectively, covering a substantial range by a factor of up to 4. Intraclass correlation coefficient ranged from 0.986 to 0.994 (average 0.991); CV ranged from 1.01% to 5.62% (average 3.6%), with yield strain and toughness having the lowest and highest CV values, respectively.

CONCLUSIONS

The data suggest that the yield and post-yield parameters have adequate reproducibility to evaluate treatment effects in interventional studies within short follow-up periods.

Micro-finite element analysis applied to high-resolution MRI reveals improved bone mechanical competence in the distal femur of female pre-professional dancers

Micro-finite element analysis applied to high-resolution (0.234-mm length scale) MRI reveals greater whole and cancellous bone stiffness, but not greater cortical bone stiffness, in the distal femur of female dancers compared to controls. Greater whole bone stiffness appears to be mediated by cancellous, rather than cortical bone adaptation.

INTRODUCTION

The purpose of this study was to compare bone mechanical competence (stiffness) in the distal femur of female dancers compared to healthy, relatively inactive female controls.

METHODS

This study had institutional review board approval. We recruited nine female modern dancers (25.7±5.8 years, 1.63±0.06 m, 57.1±4.6 kg) and ten relatively inactive, healthy female controls matched for age, height, and weight (32.1±4.8 years, 1.6±0.04 m, 55.8±5.9 kg). We scanned the distal femur using a 7-T MRI scanner and a three-dimensional fast low-angle shot sequence (TR/TE=31 ms/5.1 ms, 0.234 mm×0.234 mm×1 mm, 80 slices). We applied micro-finite element analysis to 10-mm-thick volumes of interest at the distal femoral diaphysis, metaphysis, and epiphysis to compute stiffness and cross-sectional area of whole, cortical, and cancellous bone, as well as cortical thickness. We applied two-tailed t-tests and ANCOVA to compare groups.

RESULTS

Dancers demonstrated greater whole and cancellous bone stiffness and cross-sectional area at all locations (p<0.05). Cortical bone stiffness, cross-sectional area, and thickness did not differ between groups (>0.08). At all locations, the percent of intact whole bone stiffness for cortical bone alone was lower in dancers (p<0.05). Adjustment for cancellous bone cross-sectional area eliminated significant differences in whole bone stiffness between groups (p>0.07), but adjustment for cortical bone cross-sectional area did not (p<0.03).

CONCLUSIONS

Modern dancers have greater whole and cancellous bone stiffness in the distal femur compared to controls. Elevated whole bone stiffness in dancers may be mediated via cancellous, rather than cortical bone adaptation.

Magnetic resonance of calcified tissues

MRI of the human body is largely made possible by the favorable relaxation properties of protons of water and triacyl glycerides prevalent in soft tissues. Hard tissues--key among them bone--are generally less amenable to measurement with in vivo MR imaging techniques, not so much as a result of the lower proton density but rather due to the extremely short life-times of the proton signal in water bound to solid-like entities, typically collagen, or being trapped in micro-pores. Either mechanism can enhance T2 relaxation by up to three orders of magnitude relative to their soft-tissue counterparts. Detection of these protons requires solid-state techniques that have emerged in recent years and that promise to add a new dimension to the study of hard tissues. Alternative approaches to probe calcified tissues exploit their characteristic magnetic properties. Bone, teeth and extra-osseous calcium-containing biomaterials are unique in that they are more diamagnetic than all other tissues and thus yield information indirectly by virtue of the induced magnetic fields present in their vicinity. Progress has also been made in methods allowing very high-resolution structural imaging of trabecular and cortical bone relying on detection of the surrounding soft-tissues. This brief review, much of it drawn from work conducted in the author's laboratory, seeks to highlight opportunities with focus on early-stage developments for image-based assessment of structure, function, physiology and mechanics of calcified tissues in humans via liquid and solid-state approaches, including proton, deuteron and phosphorus NMR and MRI.

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.

Computationally-optimized bone mechanical modeling from high-resolution structural images

Image-based mechanical modeling of the complex micro-structure of human bone has shown promise as a non-invasive method for characterizing bone strength and fracture risk in vivo. In particular, elastic moduli obtained from image-derived micro-finite element (μFE) simulations have been shown to correlate well with results obtained by mechanical testing of cadaveric bone. However, most existing large-scale finite-element simulation programs require significant computing resources, which hamper their use in common laboratory and clinical environments. In this work, we theoretically derive and computationally evaluate the resources needed to perform such simulations (in terms of computer memory and computation time), which are dependent on the number of finite elements in the image-derived bone model. A detailed description of our approach is provided, which is specifically optimized for μFE modeling of the complex three-dimensional architecture of trabecular bone. Our implementation includes domain decomposition for parallel computing, a novel stopping criterion, and a system for speeding up convergence by pre-iterating on coarser grids. The performance of the system is demonstrated on a dual quad-core Xeon 3.16 GHz CPUs equipped with 40 GB of RAM. Models of distal tibia derived from 3D in-vivo MR images in a patient comprising 200,000 elements required less than 30 seconds to converge (and 40 MB RAM). To illustrate the system's potential for large-scale μFE simulations, axial stiffness was estimated from high-resolution micro-CT images of a voxel array of 90 million elements comprising the human proximal femur in seven hours CPU time. In conclusion, the system described should enable image-based finite-element bone simulations in practical computation times on high-end desktop computers with applications to laboratory studies and clinical imaging.

In vivo estimation of bone stiffness at the distal femur and proximal tibia using ultra-high-field 7-Tesla magnetic resonance imaging and micro-finite element analysis

The goal of this study was to demonstrate the feasibility of using 7-Tesla (7T) magnetic resonance imaging (MRI) and micro-finite element analysis (µFEA) to evaluate mechanical and structural properties of whole, cortical, and trabecular bone at the distal femur and proximal tibia in vivo. 14 healthy subjects were recruited (age 40.7 ± 15.7 years). The right knee was scanned on a 7T MRI scanner using a 28 channel-receive knee coil and a three-dimensional fast low-angle shot sequence (TR/TE 20 ms/5.02 ms, 0.234 mm × 0.234 mm × 1 mm, 80 axial images, 7 min 9 s). Bone was analyzed at the distal femoral metaphysis, femoral condyles, and tibial plateau. Whole, cortical, and trabecular bone stiffness was computed using µFEA. Bone volume fraction (BVF), bone areas, and cortical thickness were measured. Trabecular bone stiffness (933.7 ± 433.3 MPa) was greater than cortical bone stiffness (216 ± 152 MPa) at all three locations (P < 0.05). Across locations, there were no differences in bone stiffness (whole, cortical, or trabecular). Whole, cortical, and trabecular bone stiffness correlated with BVF (R ≥ 0.69, P < 0.05) and inversely correlated with corresponding whole, cortical, and trabecular areas (R ≤ -0.54, P < 0.05), but not with cortical thickness (R < -0.11, P > 0.05). Whole, cortical, and trabecular stiffness correlated with body mass index (R ≥ 0.62, P < 0.05). In conclusion, at the distal femur and proximal tibia, trabecular bone contributes 66-74% of whole bone stiffness. 7T MRI and µFEA may be used as a method to provide insight into how structural properties of cortical or trabecular bone affect bone mechanical competence in vivo.

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 Young's 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.

Quantitative assessment of trabecular bone micro-architecture of the wrist via 7 Tesla MRI: preliminary results

OBJECT

The goal of this study was to determine the feasibility of performing quantitative 7 T magnetic resonance imaging (MRI) assessment of trabecular bone micro-architecture of the wrist, a common fracture site.

MATERIALS AND METHODS

The wrists of 4 healthy subjects (1 woman, 3 men, 28 ± 8.9 years) were scanned on a 7 T whole body MR scanner using a 3D fast low-angle shot (FLASH) sequence (TR/TE = 20/4.5 m s, 0.169 × 0.169 × 0.5 mm). Trabecular bone was segmented and divided into 4 or 8 angular subregions. Total bone volume (TBV), bone volume fraction (BVF), surface-curve ratio (SC), and erosion index (EI) were computed. Subjects were scanned twice to assess measurement reproducibility.

RESULTS

Group mean subregional values for TBV, BVF, SC, and EI (8 subregion analysis) were as follows: 8489 ± 3686, 0.27 ± 0.045, 9.61 ± 6.52; and 1.43 ± 1.25. Within each individual, there was subregional variation in TBV, SC, and EI (>5%), but not BVF (<5%). Intersubject variation (≥12%) existed for all parameters. Within-subject coefficients of variation were ≤10%.

CONCLUSION

This is the first study to perform quantitative 7T MRI assessment of trabecular bone micro-architecture of the wrist. This method could be utilized to study perturbations in bone structure in subjects with osteoporosis or other bone disorders.

Performance of μMRI-Based virtual bone biopsy for structural and mechanical analysis at the distal tibia at 7T field strength

PURPOSE

To assess the performance of a 3D fast spin echo (FSE) pulse sequence utilizing out-of-slab cancellation through phase alternation and micro-magnetic resonance imaging (μMRI)-based virtual bone biopsy processing methods to probe the serial reproducibility and sensitivity of structural and mechanical parameters of the distal tibia at 7.0T.

MATERIALS AND METHODS

The distal tibia of five healthy subjects was imaged at three timepoints with a 3D FSE sequence at 137 × 137 × 410 μm(3) voxel size. Follow-up images were retrospectively 3D registered to baseline images. Coefficients of variation (CV) and intraclass correlation coefficients (ICCs) for measures of scale and topology of the whole tibial trabecular bone (TB) cross-section as well as finite-element-derived Young's and shear moduli of central cuboidal TB subvolumes (8 × 8 × 5 mm(3) ) were evaluated as measures of reproducibility and reliability. Four additional cubic TB subregions (anterior, medial, lateral, and posterior) of similar dimensions were extracted and analyzed to determine associations between whole cross-section and subregional structural parameters.

RESULTS

The mean signal-to-noise ratio (SNR) over the 15 image acquisitions was 27.5 ± 2.1. Retrospective registration yielded an average common analysis volume of 67% across the three exams per subject. Reproducibility (mean CV = 3.6%; range, 1.5%-5%) and reliability (ICCs, 0.95-0.99) of all parameters permitted parameter-based discrimination of the five subjects in spite of the narrow age range (26-36 years) covered. Parameters characterizing topology were better able to distinguish two individuals who demonstrated similar values for scalar measurements (≈ 34% difference, P < 0.001). Whole-section axial stiffness encompassing the cortex was superior at distinguishing two individuals relative to its central subregional TB counterpart (≈ 8% difference; P < 0.05). Interregion comparisons showed that although all parameters were correlated (mean R(2) = 0.78; range 0.57-0.99), the strongest associations observed were those for the erosion index (mean R(2) = 0.95, P ≤ 0.01).

CONCLUSION

The reproducibility and structural and mechanical parameter-based discriminative ability achieved in five healthy subjects suggests that 7T-derived μMRI of TB can be applied towards serial patient studies of osteoporosis and may enable earlier detection of disease or treatment-based effects.

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.