Implications of resolution and noise for in vivo micro-MRI of trabecular bone

A Novel Two-Compartment Model for Calculating Bone Volume Fractions and Bone Mineral Densities From Computed Tomography Images

Osteoporosis is a disease characterized by a degradation of bone structures. Various methods have been developed to diagnose osteoporosis by measuring bone mineral density (BMD) of patients. However, BMDs from these methods were not equivalent and were incomparable. In addition, partial volume effect introduces errors in estimating bone volume from computed tomography (CT) images using image segmentation. In this study, a two-compartment model (TCM) was proposed to calculate bone volume fraction (BV/TV) and BMD from CT images. The TCM considers bones to be composed of two sub-materials. Various equivalent BV/TV and BMD can be calculated by applying corresponding sub-material pairs in the TCM. In contrast to image segmentation, the TCM prevented the influence of the partial volume effect by calculating the volume percentage of sub-material in each image voxel. Validations of the TCM were performed using bone-equivalent uniform phantoms, a 3D-printed trabecular-structural phantom, a temporal bone flap, and abdominal CT images. By using the TCM, the calculated BV/TVs of the uniform phantoms were within percent errors of ±2%; the percent errors of the structural volumes with various CT slice thickness were below 9%; the volume of the temporal bone flap was close to that from micro-CT images with a percent error of 4.1%. No significant difference (p >0.01) was found between the areal BMD of lumbar vertebrae calculated using the TCM and measured using dual-energy X-ray absorptiometry. In conclusion, the proposed TCM could be applied to diagnose osteoporosis, while providing a basis for comparing various measurement methods.

Preclinical assessment of a new magnetic resonance-based technique for determining bone quality by characterization of trabecular microarchitecture

The utility of HR-CT to study longitudinal changes in bone microarchitecture is limited by subject radiation exposure. Although MR is not subject to this limitation, it is limited both by patient movement that occurs during prolonged scanning at distal sites, and by the signal-to-noise ratio that is achievable for high-resolution images in a reasonable scan time at proximal sites. Recently, a novel MR-based technique, fine structure analysis (FSA) (Chase et al. Localised one-dimensional magnetic resonance spatial frequency spectroscopy. PCT/US2012/068284 2012, James and Chase Magnetic field gradient structure characteristic assessment using one-dimensional (1D) spatial frequency distribution analysis. 7932720 B2, 2011) has been developed which provides both high-resolution and fast scan times, but which generates at a designated set of spatial positions (voxels) a one-dimensional signal of spatial frequencies. Appendix 1 provides a brief introduction to FSA. This article describes an initial exploration of FSA for the rapid, non-invasive characterization of trabecular microarchitecture in a preclinical setting. For L4 vertebrae of sham and ovariectomized (OVX) rats, we compared FSA-generated metrics with those from CT datasets and from CT-derived histomorphometry parameters, trabecular number (Tb.N), bone volume density (BV/TV), trabecular thickness (Tb.Th) and trabecular separation (Tb.Sp). OVX caused a reduction of the higher frequency structures that correspond to a denser trabecular lattice, while increasing the preponderance of lower frequency structures, which correspond to a more open lattice. As one example measure, the centroid of the FSA spectrum (which we refer to as fSAcB) showed strong correlation in the same region with CT-derived histomorphometry values: Tb.Sp: r -0.63, p < 0.001; Tb.N: r 0.71, p < 0.001; BV/TV: r 0.64, p < 0.001, Tb.Th: r 0.44, p < 0.05. Furthermore, we found a 17.5% reduction in fSAcB in OVX rats (p < 0.0001). In a longitudinal study, FSA showed that the age-related increase in higher frequency structures was abolished in OVX rats, being replaced with a 78-194% increase in lower frequency structures (2.4-2.8 objects/mm range), indicating a more sparse trabecular lattice (p < 0.05). The MR-based fine structure analysis enables high-resolution, radiation-free, rapid quantification of bone structures in one dimension (the specific point and direction being chosen by the clinician) of the spine.

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.

On the significance of motion degradation in high-resolution 3D μMRI of trabecular bone

RATIONALE AND OBJECTIVES

Subtle subject movement during high-resolution three-dimensional micro-magnetic resonance imaging of trabecular bone (TB) causes blurring, thereby rendering the data unreliable for quantitative analysis. In this work, the effects of translational and rotational motion displacements were evaluated qualitatively and quantitatively.

MATERIALS AND METHODS

In experiment 1, motion was induced by applying various simulated and previously observed in vivo trajectories as phase shifts to k-space or rotation angles to k-space segments of a virtually motion-free data set. In experiment 2, images that were visually free of motion artifacts from two groups of 10 healthy individuals, differing in age, were selected to probe the effects of motion on TB parameters. In both experiments, images were rated for motion severity, and the scores were compared to a focus criterion, the normalized gradient squared.

RESULTS

Strong correlations were observed between the motion quality scores and the corresponding normalized gradient squared values (R(2) = 0.52-0.64, P < .01). The results from experiment 1 demonstrated consistently lower image quality and alterations in structural parameters of 9% to 45% with increased amplitude of displacements. In experiment 2, the significant differences in structural parameter group means of the motion-free images were lost upon motion degradation. Autofocusing, a postprocessing correction method, partially recovered the sharpness of the original motion-free images in 13 of 20 subjects.

CONCLUSIONS

Quantitative TB structural measures are highly sensitive to subtle motion-induced degradation, which adversely affects precision and statistical power. The results underscore the influence of subject movement in high-resolution three-dimensional micro-magnetic resonance imaging and its correction for TB structure analysis.

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.

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.

Helmholtz-pair transmit coil with integrated receive array for high-resolution MRI of trabecular bone in the distal tibia at 7T

A Helmholtz-pair local transmit RF coil with an integrated four-element receive array RF coil and foot immobilization platform was designed and constructed for imaging the distal tibia in a whole-body 7T MRI scanner. Simulations and measurements of the B(1) field distribution of the transmit coil are described, along with SAR considerations for operation at 7T. Results of imaging the trabecular bone of three volunteers at 1.5T, 3T and 7T are presented, using identical 1.5T and 3T versions of the 7T four-element receive array. The spatially registered images reveal improved visibility for individual trabeculae and show average gains in SNR of 2.8× and 4.9× for imaging at 7T compared to 3T and 1.5T, respectively. The results thus display an approximately linear dependence of SNR with field strength and enable the practical utility of 7T scanners for micro-MRI of trabecular bone.

Radiological diagnostic progress in skeletal diseases

High-resolution bone imaging has made tremendous progress in the recent past. Both imaging modalities, computed tomography as well as MR imaging, have improved image quality. New developments such as HR-pQCT now make it possible to acquire in vivo images at peripheral sites with isotropic voxel size in a very short time. Further enhancements in the MR field have made it possible to image more central body sites such as the proximal femur with very high spatial resolution. New analysis methods can obtain direct estimates of biomechanical properties and important information related to bone's topology, as well as parameters of scale and orientation. These accomplishments will be essential in the noninvasive assessment of osteoporosis and fracture risk, will provide insight into the mechanisms behind bone loss, and will increasingly play a role as a tool for assessing treatment efficacy.

Non-destructive characterization of microdamage in cortical bone using low field pulsed NMR

The microcracking and damage accumulation process in human cortical bone was characterized by performing cyclic loading under four-point bending at ambient temperature. A non-destructive nuclear magnetic resonance (NMR) spin-spin (T(2)) relaxation technique was applied to quantify the apparent changes in bone porosity as a function of cyclic loading and prior damage accumulation, first to unloaded cortical bone to quantify the initial porosity and then to fatigued cortical bone that was subjected to cyclic loading to various levels of modulus degradation and microdamage in the form of microcracks. The NMR T(2) relaxation time and amplitude data of the fatigued bone were compared against the undamaged state. The difference in the T(2) relaxation time data was taken as a measure of the increase in pore size, bone porosity or microcrack density due to microdamage induced by cyclic loading. A procedure was developed to deduce the number and size distributions of microcracks formed in cortical bone. Serial sectioning of the fatigued bone showed the formation of microcracks along the cement lines or within the interstitial tissue. The results on the evolution of microdamage derived from NMR measurements were verified by independent experimental measurements of microcrack density using histological characterization techniques. The size distribution and population of the microcracks were then utilized in conjunction with an analytical model to predict the degradation of the elastic modulus of cortical bone as a function of damage accumulation.

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

PURPOSE

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.

MATERIALS AND METHODS

The left distal tibial metaphysis of seven healthy volunteers was imaged at three time-points using a FLASE (fast large-angle spin-echo) pulse sequence at 137 x 137 x 410 mum(3) and (160 mum)(3) voxel sizes. Image artifacts, motion degradation, and serial image volume misalignments were controlled to maximize reproducibility of image-derived measures of scale, topology, orientation in terms of structural anisotropy, and finite-element derived Young's and shear moduli. Coefficients of variation (CV) and intraclass correlation coefficients (ICC) for structural and mechanical parameters were evaluated as measures of reproducibility and reliability. The ability of structural and mechanical parameters to distinguish between subjects was tested by analysis of variance.

RESULTS

Reproducibility was generally higher in the anisotropic data (CVs 1-5% versus 1-9% for isotropic images). Anisotropic voxel size yielded greater measurement reliability (ICCs 0.75-0.99, mean = 0.92 versus 0.62-0.99, mean = 0.86) and better discrimination of the seven subjects (75% versus 50% of the possible comparisons were significantly different [P < 0.05]) except for measures of structural anisotropy and topology. Isotropic resolution improved detection of structural orientation and permitted visualization of small perforations in longitudinal trabecular plates not detected at anisotropic resolution.

CONCLUSION

Improved image acquisition and processing techniques enhance reproducibility of structural and mechanical parameters derived from high-resolution MRI of metaphyseal bone in the distal tibia.

High-resolution imaging techniques for the assessment of osteoporosis

The importance of assessing the bone's microarchitectural make-up in addition to its mineral density in the context of osteoporosis has been emphasized in several publications. The high spatial resolution required to resolve the bone's microstructure in a clinically feasible scan time is challenging. At present, the best suited modalities meeting these requirements in vivo are high-resolution peripheral quantitative imaging (HR-pQCT) and magnetic resonance imaging (MRI). Whereas HR-pQCT is limited to peripheral skeleton regions like the wrist and ankle, MRI can also image other sites like the proximal femur but usually with lower spatial resolution. In addition, multidetector computed tomography has been used for high-resolution imaging of trabecular bone structure; however, the radiation dose is a limiting factor. This article provides an overview of the different modalities, technical requirements, and recent developments in this emerging field. Details regarding imaging protocols as well as image postprocessing methods for bone structure quantification are discussed.

Classification of trabeculae into three-dimensional rodlike and platelike structures via local inertial anisotropy

Trabecular bone microarchitecture is a significant determinant of the bone's mechanical properties and is thus of major clinical relevance in predicting fracture risk. The three-dimensional nature of trabecular bone is characterized by parameters describing scale, topology, and orientation of structural elements. However, none of the current methods calculates all three types of parameters simultaneously and in three dimensions. Here the authors present a method that produces a continuous classification of voxels as belonging to platelike or rodlike structures that determines their orientation and estimates their thickness. The method, dubbed local inertial anisotropy (LIA), treats the image as a distribution of mass density and the orientation of trabeculae is determined from a locally calculated tensor of inertia at each voxel. The orientation entropies of rods and plates are introduced, which can provide new information about microarchitecture not captured by existing parameters. The robustness of the method to noise corruption, resolution reduction, and image rotation is demonstrated. Further, the method is compared with established three-dimensional parameters including the structure-model index and topological surface-to-curve ratio. Finally, the method is applied to data acquired in a previous translational pilot study showing that the trabecular bone of untreated hypogonadal men is less platelike than that of their eugonadal peers.