Feasibility of assessing bone matrix and mineral properties in vivo by combined solid-state 1H and 31P MRI

Quantitative 1H Magnetic Resonance Imaging on Normal and Pathologic Rat Bones by Solid-State 1H ZTE Sequence with Water and Fat Suppression

BACKGROUND

Osteoporosis (OP) and osteomalacia (OM) are metabolic bone diseases characterized by mineral and matrix density changes. Quantitative bone matrix density differentiates OM from OP. MRI is a noninvasive and nonionizing imaging technique that can measure bone matrix density quantitatively in ex vivo and in vivo.

PURPOSE

To demonstrate water + fat suppressed 1H MRI to compute bone matrix density in ex vivo rat femurs in the preclinical model.

STUDY TYPE

Prospective.

ANIMAL MODEL

Fifteen skeletally mature female Sprague-Dawley rats, five per group (normal, ovariectomized (OVX), partially nephrectomized/vitamin D (Vit-D) deficient), 250-275 g, ∼15 weeks old.

FIELD STRENGTH/SEQUENCE

7T, zero echo time sequence with water + fat (VAPOR) suppression capability, μCT imaging, and gravimetric measurements.

ASSESSMENT

Cortical and trabecular bone segments from normal and disease models were scanned in the same coil along with a dual calibration phantom for quantitative assessment of bone matrix density.

STATISTICAL TESTS

ANOVA and linear regression were used for data analysis, with P-values <0.05 statistically significant.

RESULTS

The MRI-derived three-density PEG pellet densities have a strong linear relationship with physical density measures (r2 = 0.99). The Vit-D group had the lowest bone matrix density for cortical bone (0.47 ± 0.16 g cm-3), whereas the OVX had the lowest bone matrix density for trabecular bone (0.26 ± 0.04 g cm-3). Gravimetry results confirmed these MRI-based observations for Vit-D cortical (0.51 ± 0.07 g cm-3) and OVX trabecular (0.26 ± 0.03 g cm-3) bone groups.

DATA CONCLUSION

Rat femur images were obtained using a modified pulse sequence and a custom-designed double-tuned (1H/31P) transmit-receive solenoid-coil on a 7T preclinical MRI scanner. Phantom experiments confirmed a strong linear relation between MRI-derived and physical density measures and quantitative bone matrix densities in rat femurs from normal, OVX, and Vit-D deficient/partially nephrectomized animals were computed.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Quantitative 31P magnetic resonance imaging on pathologic rat bones by ZTE at 7T

BACKGROUND

Osteoporosis is characterized by low bone mineral density (BMD), which predisposes individuals to frequent fragility fractures. Quantitative BMD measurements can potentially help distinguish bone pathologies and allow clinicians to provide disease-relieving therapies. Our group has developed non-invasive and non-ionizing magnetic resonance imaging (MRI) techniques to measure bone mineral density quantitatively. Dual-energy X-ray Absorptiometry (DXA) is a clinically approved non-invasive modality to diagnose osteoporosis but has associated disadvantages and limitations.

PURPOSE

Evaluate the clinical feasibility of phosphorus (31P) MRI as a non-invasive and non-ionizing medical diagnostic tool to compute bone mineral density to help differentiate between different metabolic bone diseases.

MATERIALS AND METHODS

Fifteen ex-vivo rat bones in three groups [control, ovariectomized (osteoporosis), and vitamin-D deficient (osteomalacia - hypo-mineralized) were scanned to compute BMD. A double-tuned (1H/31P) transmit-receive single RF coil was custom-designed and in-house-built with a better filling factor and strong radiofrequency (B1) field to acquire solid-state 31P MR images from rat femurs with an optimum signal-to-noise ratio (SNR). Micro-computed tomography (μCT) and gold-standard gravimetric analyses were performed to compare and validate MRI-derived bone mineral densities.

RESULTS

Three-dimensional 31P MR images of rat bones were obtained with a zero-echo-time (ZTE) sequence with 468 μm spatial resolution and 12-17 SNR on a Bruker 7 T Biospec having multinuclear capability. BMD was measured quantitatively on cortical and trabecular bones with a known standard reference. A strong positive correlation (R = 0.99) and a slope close to 1 in phantom measurements indicate that the densities measured by 31P ZTE MRI are close to the physical densities in computing quantitative BMD. The 31P NMR properties (resonance linewidth of 4 kHz and T1 of 67 s) of ex-vivo rat bones were measured, and 31P ZTE imaging parameters were optimized. The BMD results obtained from MRI are in good agreement with μCT and gravimetry results.

CONCLUSION

Quantitative measurements of BMD on ex-vivo rat femurs were successfully conducted on a 7 T preclinical scanner. This study suggests that quantitative measurements of BMD are feasible on humans in clinical MRI with suitable hardware, RF coils, and pulse sequences with optimized parameters within an acceptable scan time since human femurs are approximately ten times larger than rat femurs. As MRI provides quantitative in-vivo data, various systemic musculoskeletal conditions can be diagnosed potentially in humans.

Bone Biomarkers Based on Magnetic Resonance Imaging

Magnetic resonance imaging (MRI) is increasingly used to evaluate the microstructural and compositional properties of bone. MRI-based biomarkers can characterize all major compartments of bone: organic, water, fat, and mineral components. However, with a short apparent spin-spin relaxation time (T2*), bone is invisible to conventional MRI sequences that use long echo times. To address this shortcoming, ultrashort echo time MRI sequences have been developed to provide direct imaging of bone and establish a set of MRI-based biomarkers sensitive to the structural and compositional changes of bone. This review article describes the MRI-based bone biomarkers representing total water, pore water, bound water, fat fraction, macromolecular fraction in the organic matrix, and surrogates for mineral density. MRI-based morphological bone imaging techniques are also briefly described.

In Vivo Assessment of Bone Quality Without X-rays

PURPOSE OF REVIEW

This review summarizes recent advances in the assessment of bone quality using non-X-ray techniques.

RECENT FINDINGS

Quantitative ultrasound (QUS) provides multiple measurements of bone characteristics based on the propagation of sound through bone, the attenuation of that sound, and different processing techniques. QUS parameters and model predictions based on backscattered signals can discriminate non-fracture from fracture cases with accuracy comparable to standard bone mineral density (BMD). With advances in magnetic resonance imaging (MRI), bound water and pore water, or a porosity index, can be quantified in several long bones in vivo. Since such imaging-derived measurements correlate with the fracture resistance of bone, they potentially provide new BMD-independent predictors of fracture risk. While numerous measurements of mineral, organic matrix, and bound water by Raman spectroscopy correlate with the strength and toughness of cortical bone, the clinical assessment of person's bone quality using spatially offset Raman spectroscopy (SORS) requires advanced spectral processing techniques that minimize contaminating signals from fat, skin, and blood. Limiting exposure of patients to ionizing radiation, QUS, MRI, and SORS has the potential to improve the assessment of fracture risk and track changes of new therapies that target bone matrix and micro-structure.

Automated, calibration-free quantification of cortical bone porosity and geometry in postmenopausal osteoporosis from ultrashort echo time MRI and deep learning

BACKGROUND

Assessment of cortical bone porosity and geometry by imaging in vivo can provide useful information about bone quality that is independent of bone mineral density (BMD). Ultrashort echo time (UTE) MRI techniques of measuring cortical bone porosity and geometry have been extensively validated in preclinical studies and have recently been shown to detect impaired bone quality in vivo in patients with osteoporosis. However, these techniques rely on laborious image segmentation, which is clinically impractical. Additionally, UTE MRI porosity techniques typically require long scan times or external calibration samples and elaborate physics processing, which limit their translatability. To this end, the UTE MRI-derived Suppression Ratio has been proposed as a simple-to-calculate, reference-free biomarker of porosity which can be acquired in clinically feasible acquisition times.

PURPOSE

To explore whether a deep learning method can automate cortical bone segmentation and the corresponding analysis of cortical bone imaging biomarkers, and to investigate the Suppression Ratio as a fast, simple, and reference-free biomarker of cortical bone porosity.

METHODS

In this retrospective study, a deep learning 2D U-Net was trained to segment the tibial cortex from 48 individual image sets comprised of 46 slices each, corresponding to 2208 training slices. Network performance was validated through an external test dataset comprised of 28 scans from 3 groups: (1) 10 healthy, young participants, (2) 9 postmenopausal, non-osteoporotic women, and (3) 9 postmenopausal, osteoporotic women. The accuracy of automated porosity and geometry quantifications were assessed with the coefficient of determination and the intraclass correlation coefficient (ICC). Furthermore, automated MRI biomarkers were compared between groups and to dual energy X-ray absorptiometry (DXA)- and peripheral quantitative CT (pQCT)-derived BMD. Additionally, the Suppression Ratio was compared to UTE porosity techniques based on calibration samples.

RESULTS

The deep learning model provided accurate labeling (Dice score 0.93, intersection-over-union 0.88) and similar results to manual segmentation in quantifying cortical porosity (R2 ≥ 0.97, ICC ≥ 0.98) and geometry (R2 ≥ 0.82, ICC ≥ 0.75) parameters in vivo. Furthermore, the Suppression Ratio was validated compared to established porosity protocols (R2 ≥ 0.78). Automated parameters detected age- and osteoporosis-related impairments in cortical bone porosity (P ≤ .002) and geometry (P values ranging from <0.001 to 0.08). Finally, automated porosity markers showed strong, inverse Pearson's correlations with BMD measured by pQCT (|R| ≥ 0.88) and DXA (|R| ≥ 0.76) in postmenopausal women, confirming that lower mineral density corresponds to greater porosity.

CONCLUSION

This study demonstrated feasibility of a simple, automated, and ionizing-radiation-free protocol for quantifying cortical bone porosity and geometry in vivo from UTE MRI and deep learning.

MRI Quantification of Cortical Bone Porosity, Mineralization, and Morphologic Structure in Postmenopausal Osteoporosis

Background Preclinical studies have suggested that solid-state MRI markers of cortical bone porosity, morphologic structure, mineralization, and osteoid density are useful measures of bone health. Purpose To explore whether MRI markers of cortical bone porosity, morphologic structure, mineralization, and osteoid density are affected in postmenopausal osteoporosis (OP) and to examine associations between MRI markers and bone mineral density (BMD) in postmenopausal women. Materials and Methods In this single-center study, postmenopausal women were prospectively recruited from January 2019 to October 2020 into two groups: participants with OP who had not undergone treatment, defined as having any dual-energy x-ray absorptiometry (DXA) T-score of -2.5 or less, and age-matched control participants without OP (hereafter, non-OP). Participants underwent MRI in the midtibia, along with DXA in the hip and spine, and peripheral quantitative CT in the midtibia. Specifically, MRI measures of cortical bone porosity (pore water and total water), osteoid density (bound water [BW]), morphologic structure (cortical bone thickness), and mineralization (phosphorous [P] density [31P] and 31P-to-BW concentration ratio) were quantified at 3.0 T. MRI measures were compared between OP and non-OP groups and correlations with BMD were assessed. Results Fifteen participants with OP (mean age, 63 years ± 5 [SD]) and 19 participants without OP (mean age, 65 years ± 6) were evaluated. The OP group had elevated pore water (11.6 mol/L vs 9.5 mol/L; P = .007) and total water densities (21.2 mol/L vs 19.7 mol/L; P = .03), and had lower cortical bone thickness (4.8 mm vs 5.6 mm; P < .001) and 31P density (6.4 mol/L vs 7.5 mol/L; P = .01) than the non-OP group, respectively, although there was no evidence of a difference in BW or 31P-to-BW concentration ratio. Pore and total water densities were inversely associated with DXA and peripheral quantitative CT BMD (P < .001), whereas cortical bone thickness and 31P density were positively associated with DXA and peripheral quantitative CT BMD (P = .01). BW, 31P density, and 31P-to-BW concentration ratio were positively associated with DXA (P < .05), but not with peripheral quantitative CT. Conclusion Solid-state MRI of cortical bone was able to help detect potential impairments in parameters reflecting porosity, morphologic structure, and mineralization in postmenopausal osteoporosis. © RSNA, 2023 Supplemental material is available for this article. See also the editorial by Bae in this issue.

MRI-based porosity index (PI) and suppression ratio (SR) in the tibial cortex show significant differences between normal, osteopenic, and osteoporotic female subjects

Introduction

Ultrashort echo time (UTE) MRI enables quantitative assessment of cortical bone. The signal ratio in dual-echo UTE imaging, known as porosity index (PI), as well as the signal ratio between UTE and inversion recovery UTE (IR-UTE) imaging, known as the suppression ratio (SR), are two rapid UTE-based bone evaluation techniques developed to reduce the time demand and cost in future clinical studies. The goal of this study was to investigate the performance of PI and SR in detecting bone quality differences between subjects with osteoporosis (OPo), osteopenia (OPe), and normal bone (Normal).

Methods

Tibial midshaft of fourteen OPe (72 ± 6 years old), thirty-one OPo (72 ± 6 years old), and thirty-seven Normal (36 ± 19 years old) subjects were scanned using dual-echo UTE and IR-UTE sequences on a clinical 3T scanner. Measured PI, SR, and bone thickness were compared between OPo, OPe, and normal bone (Normal) subjects using the Kruskal-Wallis test by ranks. Spearman's rank correlation coefficients were calculated between dual-energy x-ray absorptiometry (DEXA) T-score and UTE-MRI results.

Results

PI was significantly higher in the OPo group compared with the Normal (24.1%) and OPe (16.3%) groups. SR was significantly higher in the OPo group compared with the Normal (41.5%) and OPe (21.8%) groups. SR differences between the OPe and Normal groups were also statistically significant (16.2%). Cortical bone was significantly thinner in the OPo group compared with the Normal (22.0%) and OPe (13.0%) groups. DEXA T-scores in subjects were significantly correlated with PI (R=-0.32), SR (R=-0.50), and bone thickness (R=0.51).

Discussion

PI and SR, as rapid UTE-MRI-based techniques, may be useful tools to detect and monitor bone quality changes, in addition to bone morphology, in individuals affected by osteoporosis.

Interleaved and simultaneous multi-nuclear magnetic resonance in vivo. Review of principles, applications and potential

Magnetic resonance signals from different nuclei can be excited or received at the same time,rendering simultaneous or rapidly interleaved multi-nuclear acquisitions feasible. The advan-tages are a reduction of total scan time compared to sequential multi-nuclear acquisitions or that additional information from heteronuclear data is obtained at thesame time and anatomical position. Information content can be qualitatively increased by delivering a more comprehensive MR-based picture of a transient state (such as an exercise bout). Also, combiningnon-proton MR acquisitions with ¹ Hinformation (e.g., dynamic shim updates and motion correction) can be used to improve data quality during long scans and benefits image coregistration. This work reviews the literature on interleaved and simultaneous multi-nuclear MRI and MRS in vivo. Prominent use cases for this methodology in clinical and research applications are brain and muscle, but studies have also been carried out in other targets, including the lung, knee, breast and heart. Simultaneous multi-nuclear measurements in the liver and kidney have also been performed, but exclusively in rodents. In this review, a consistent nomenclature is proposed, to help clarify the terminology used for this principle throughout the literature on in-vivo MR. An overview covers the basic principles, the technical requirements on the MR scanner and the implementations realised either by MR system vendors or research groups, from the early days until today. Considerations regarding the multi-tuned RF coils required and heteronuclear polarisation interactions are briefly discussed, and fields for future in-vivo applications for interleaved multi-nuclear MR pulse sequences are identified.

Bone hydration: How we can evaluate it, what can it tell us, and is it an effective therapeutic target?

Water constitutes roughly a quarter of the cortical bone by volume yet can greatly influence mechanical properties and tissue quality. There is a growing appreciation for how water can dynamically change due to age, disease, and treatment. A key emerging area related to bone mechanical and tissue properties lies in differentiating the role of water in its four different compartments, including free/pore water, water loosely bound at the collagen/mineral interfaces, water tightly bound within collagen triple helices, and structural water within the mineral. This review summarizes our current knowledge of bone water across the four functional compartments and discusses how alterations in each compartment relate to mechanical changes. It provides an overview on the advent of- and improvements to- imaging and spectroscopic techniques able to probe nano-and molecular scales of bone water. These technical advances have led to an emerging understanding of how bone water changes in various conditions, of which aging, chronic kidney disease, diabetes, osteoporosis, and osteogenesis imperfecta are reviewed. Finally, it summarizes work focused on therapeutically targeting water to improve mechanical properties.

Comparison of magnetic resonance imaging and computed tomography for bone assessment of neurogenic heterotopic ossification of the hip: a preliminary study

Background

Neurogenic heterotopic ossification (NHO) is a frequent complication, often involving the hip. The functional impact may require surgical management and pre-surgical imaging assessment is necessary, usually by computed tomography (CT). We aimed to compare the performances of magnetic resonance imaging (MRI) and CT for bone assessment on pre-surgical imaging of the heterotopic ossifications and their features in NHO of the hip.

Methods

This single-center prospective preliminary study included all patients who underwent surgery for NHO with joint limitation from July 2019 to March 2020. All patients had a CT after biphasic iodinated solution injection and an MRI including T1-weighted, STIR and ZTE sequences. Standardized reports were completed for both exams for each patient, evaluating location, implantation and fragmentation of NHO, relation to the joint capsule and bone mineralization, then were compared.

Results

Seven patients from 32 to 70 years old (mean = 50.2 ± 17.2 years) were evaluated. NHO were bilateral in 2 patients, for a total of nine hips: six right hips and three left hips. Observed concordance rates between MRI and CT were, respectively, 94.4% for location, 100% for circumferential extension, 87.3% for implantation 88.9% for fragmentation, 77.8% for relation to the joint capsule and 66.7% for bone mineralization. It was 100% for femoral neck fracture and osteonecrosis of the femoral head.

Conclusion

This preliminary study suggests that pre-surgical MRI imaging should be considered as effective as CT for bone assessment of NHO and their features.

Trial registration: ClinicalTrials.gov, NCT03832556. Registered February 6, 2019, https://clinicaltrials.gov/ct2/show/NCT03832556.

Advancements in Osteoporosis Imaging, Screening, and Study of Disease Etiology

PURPOSE OF REVIEW

The purpose of this review is to inform researchers and clinicians with the most recent imaging techniques that are employed (1) to opportunistically screen for osteoporosis and (2) to provide a better understanding into the disease etiology of osteoporosis.

RECENT FINDINGS

Phantomless calibration techniques for computed tomography (CT) may pave the way for better opportunistic osteoporosis screening and the retroactive analysis of imaging data. Additionally, hardware advances are enabling new applications of dual-energy CT and cone-beam CT to the study of bone. Advances in MRI sequences are also improving imaging evaluation of bone properties. Finally, the application of image registration techniques is enabling new uses of imaging to investigate soft tissue-bone interactions as well as bone turnover. While DXA remains the most prominent imaging tool for osteoporosis diagnosis, new imaging techniques are becoming more widely available and providing additional information to inform clinical decision-making.

Survey of MRI Usefulness for the Clinical Assessment of Bone Microstructure

Bone microarchitecture has been shown to provide useful information regarding the evaluation of skeleton quality with an added value to areal bone mineral density, which can be used for the diagnosis of several bone diseases. Bone mineral density estimated from dual-energy X-ray absorptiometry (DXA) has shown to be a limited tool to identify patients’ risk stratification and therapy delivery. Magnetic resonance imaging (MRI) has been proposed as another technique to assess bone quality and fracture risk by evaluating the bone structure and microarchitecture. To date, MRI is the only completely non-invasive and non-ionizing imaging modality that can assess both cortical and trabecular bone in vivo. In this review article, we reported a survey regarding the clinically relevant information MRI could provide for the assessment of the inner trabecular morphology of different bone segments. The last section will be devoted to the upcoming MRI applications (MR spectroscopy and chemical shift encoding MRI, solid state MRI and quantitative susceptibility mapping), which could provide additional biomarkers for the assessment of bone microarchitecture.

MRI-derived porosity index is associated with whole-bone stiffness and mineral density in human cadaveric femora

Ultrashort echo time (UTE) magnetic resonance imaging (MRI) measures proton signals in cortical bone from two distinct water pools, bound water, or water that is tightly bound to bone matrix, and pore water, or water that is freely moving in the pore spaces in bone. By isolating the signal contribution from the pore water pool, UTE biomarkers can directly quantify cortical bone porosity in vivo. The Porosity Index (PI) is one non-invasive, clinically viable UTE-derived technique that has shown strong associations in the tibia with μCT porosity and other UTE measures of bone water. However, the efficacy of the PI biomarker has never been examined in the proximal femur, which is the site of the most catastrophic osteoporotic fractures. Additionally, the loads experienced during a sideways fall are complex and the femoral neck is difficult to image with UTE, so the usefulness of the PI in the femur was unknown. Therefore, the aim of this study was to examine the relationships between the PI measure in the proximal cortical shaft of human cadaveric femora specimens compared to (1) QCT-derived bone mineral density (BMD) and (2) whole bone stiffness obtained from mechanical testing mimicking a sideways fall. Fifteen fresh, frozen whole cadaveric femora specimens (age 72.1 ± 15.0 years old, 10 male, 5 female) were scanned on a clinical 3-T MRI using a dual-echo UTE sequence. Specimens were then scanned on a clinical CT scanner to measure volumetric BMD (vBMD) and then non-destructively mechanically tested in a sideways fall configuration. The PI in the cortical shaft demonstrated strong correlations with bone stiffness (r = -0.82, P = 0.0014), CT-derived vBMD (r = -0.64, P = 0.0149), and with average cortical thickness (r = -0.60, P = 0.0180). Furthermore, a hierarchical regression showed that PI was a strong predictor of bone stiffness which was independent of the other parameters. The findings from this study validate the MRI-derived porosity index as a useful measure of whole-bone mechanical integrity and stiffness.

Quantitative Ultrashort Echo Time (UTE) Magnetic Resonance Imaging of Bone: An Update

Bone possesses a highly complex hierarchical structure comprised of mineral (~45% by volume), organic matrix (~35%) and water (~20%). Water exists in bone in two forms: 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. Magnetic resonance (MR) imaging has been increasingly used for assessment of cortical and trabecular bone. However, bone appears as a signal void on conventional MR sequences because of its short T2*. Ultrashort echo time (UTE) sequences with echo times (TEs) 100-1,000 times shorter than those of conventional sequences allow direct imaging of BW and PW in bone. A series of quantitative UTE MRI techniques has been developed for bone evaluation. UTE and adiabatic inversion recovery prepared UTE (IR-UTE) sequences have been developed to quantify BW and PW. UTE magnetization transfer (UTE-MT) sequences have been developed to quantify collagen backbone protons, and UTE quantitative susceptibility mapping (UTE-QSM) sequences have been developed to assess bone mineral.

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.

Magnetic Resonance Imaging of Hard Tissues and Hard Tissue Engineered Bio-substitutes

Magnetic resonance imaging (MRI) is a non-invasive diagnostic imaging tool based on the detection of protons into the tissues. This imaging technique is remarkable because of high spatial resolution, strong soft tissue contrast and specificity, and good depth penetration. However, MR imaging of hard tissues, such as bone and teeth, remains challenging due to low proton content in such tissues as well as to very short transverse relaxation times (T₂). To overcome these issues, new MRI techniques, such as sweep imaging with Fourier transformation (SWIFT), ultrashort echo time (UTE) imaging, and zero echo time (ZTE) imaging, have been developed for hard tissues imaging with promising results reported. Within this article, MRI techniques developed for the detection of hard tissues, such as bone and dental tissues, have been reviewed. The main goal was thus to give a comprehensive overview on the corresponding (pre-) clinical applications and on the potential future directions with such techniques applied. In addition, a section dedicated to MR imaging of novel biomaterials developed for hard tissue applications was given as well.

Impact of gradient imperfections on bone water quantification with UTE MRI

PURPOSE

The impact of gradient imperfections on UTE images and UTE image-derived bone water quantification was investigated at 3 T field strength.

METHODS

The effects of simple gradient time delays and eddy currents on UTE images, as well as the effects of gradient error corrections, were studied with simulation and phantom experiments. The k-space trajectory was mapped with a 2D sequence with phase encoding on both spatial axes by measuring the phase of the signal in small time increments during ramp-up of the read gradient. In vivo 3D UTE images were reconstructed with and without gradient error compensation to determine the bias in bone water quantification. Finally, imaging was performed on 2 equally configured Siemens TIM Trio systems (Siemens Medical Solutions, Erlangen, Germany) to investigate the impact of such gradient imperfections on inter-scanner measurement bias.

RESULTS

Compared to values derived from UTE images with full gradient error compensation, total bone water was found to deviate substantially with no (up to 17%) or partial (delay-only) compensation (up to 10.8%). Bound water, obtained with inversion recovery-prepared UTE, was somewhat less susceptible to gradient errors (up to 2.2% for both correction strategies). Inter-scanner comparison indicated a statistically significant bias between measurements from the 2 MR systems for both total and bound water, which either vanished or was substantially reduced following gradient error correction.

CONCLUSION

Gradient imperfections impose spatially dependent artifacts on UTE images, which compromise not only bone water quantification accuracy but also inter-scanner measurement agreement if left uncompensated.

MRI-based assessment of proximal femur strength compared to mechanical testing

Half of the women who sustain a hip fracture would not qualify for osteoporosis treatment based on current DXA-estimated bone mineral density criteria. Therefore, a better approach is needed to determine if an individual is at risk of hip fracture from a fall. The objective of this study was to determine the association between radiation-free MRI-derived bone strength and strain simulations compared to results from direct mechanical testing of cadaveric femora. Imaging was conducted on a 3-Tesla MRI scanner using two sequences: one balanced steady-state free precession sequence with 300 μm isotropic voxel size and one spoiled gradient echo with anisotropic voxel size of 234 × 234 × 1500 μm. Femora were dissected free of soft-tissue and 4350-ohm strain-gauges were securely applied to surfaces at the femoral shaft, inferior neck, greater trochanter, and superior neck. Cadavers were mechanically tested with a hydraulic universal test frame to simulate loading in a sideways fall orientation. Sideways fall forces were simulated on MRI-based finite element meshes and bone stiffness, failure force, and force for plastic deformation were computed. Simulated bone strength metrics from the 300 μm isotropic sequence showed strong agreement with experimentally obtained values of bone strength, with stiffness (r = 0.88, p = 0.0002), plastic deformation point (r = 0.89, p < 0.0001), and failure force (r = 0.92, p < 0.0001). The anisotropic sequence showed similar trends for stiffness, plastic deformation point, and failure force (r = 0.68, 0.70, 0.84; p = 0.02, 0.01, 0.0006, respectively). Surface strain-gauge measurements showed moderate to strong agreement with simulated magnitude strain values at the greater trochanter, superior neck, and inferior neck (r = -0.97, -0.86, 0.80; p ≤0.0001, 0.003, 0.03, respectively). The findings from this study support the use of MRI-based FE analysis of the hip to reliably predict the mechanical competence of the human femur in clinical settings.

Correlations of cortical bone microstructural and mechanical properties with water proton fractions obtained from ultrashort echo time (UTE) MRI tricomponent T2* model

Mechanical and microstructural evaluations of cortical bone using ultrashort echo time magnetic resonance imaging (UTE-MRI) have been performed increasingly in recent years. UTE-MRI acquires considerable signal from cortical bone and enables quantitative bone evaluations. Fitting bone apparent transverse magnetization (T2*) decay using a bicomponent model has been regularly performed to estimate bound water (BW) and pore water (PW) in the quantification of bone matrix and porosity, respectively. Human cortical bone possesses a considerable amount of fat, which appears as MRI T2* signal oscillation and can subsequently lead to BW overestimation when using a bicomponent model. Tricomponent T2* fitting model has been developed to improve BW and PW estimations by accounting for fat contribution in the MRI signal. This study aimed to investigate the correlations of microstructural and mechanical properties of human cortical bone with water pool fractions obtained from a tricomponent T2* model. 135 cortical bone strips (~4 × 2 × 40 mm³ ) from tibial and femoral midshafts of 37 donors (61 ± 24 years old) were scanned using ten sets of dual-echo 3D-UTE-Cones sequences (TE = 0.032-24.0 ms) on a 3 T MRI scanner for T2* fitting analyses. Average bone porosity and pore size were measured using microcomputed tomography (μCT) at 9 μm voxel size. Bone mechanical properties were measured using 4-point bending tests. Using a tricomponent model, bound water fraction (Frac_BW ) showed significant strong (R = 0.70, P < 0.01) and moderate (R = 0.58-0.62, P < 0.01) correlations with porosity and mechanical properties, respectively. Correlations of bone microstructural and mechanical properties with water pool fractions were higher for tricomponent model results compared with the bicomponent model. The tricomponent T2* fitting model is suggested as a useful technique for cortical bone evaluation where the MRI contribution of bone fat is accounted for.

Water proton density in human cortical bone obtained from ultrashort echo time (UTE) MRI predicts bone microstructural properties

PURPOSE

To investigate the correlations between cortical bone microstructural properties and total water proton density (TWPD) obtained from three-dimensional ultrashort echo time Cones (3D-UTE-Cones) magnetic resonance imaging techniques.

MATERIALS AND METHODS

135 cortical bone samples were harvested from human tibial and femoral midshafts of 37 donors (61 ± 24 years old). Samples were scanned using 3D-UTE-Cones sequences on a clinical 3T MRI and on a high-resolution micro-computed tomography (μCT) scanner. TWPD was measured using 3D-UTE-Cones MR images. Average bone porosity, pore size, and bone mineral density (BMD) were measured from μCT images at 9 μm voxel size. Pearson's correlation coefficients between TWPD and μCT-based measures were calculated.

RESULTS

TWPD showed significant moderate correlation with both average bone porosity (R = 0.66, p < 0.01) and pore size (R = 0.57, p < 0.01). TWPD also showed significant strong correction with BMD (R = 0.71, p < 0.01).

CONCLUSIONS

The presented 3D-UTE-Cones imaging technique allows assessment of TWPD in human cortical bone. This quick UTE-MRI-based technique was capable of predicting bone microstructure differences with significant correlations. Such correlations highlight the potential of UTE-MRI-based measurement of bone water proton density to assess bone microstructure.

An Update in Qualitative Imaging of Bone Using Ultrashort Echo Time Magnetic Resonance

Bone is comprised of mineral, collagenous organic matrix, and water. X-ray-based techniques are the standard approach for bone evaluation in clinics, but they are unable to detect the organic matrix and water components in bone. Magnetic resonance imaging (MRI) is being used increasingly for bone evaluation. While MRI can non-invasively assess the proton pools in soft tissues, cortical bone typically appears as a signal void with clinical MR techniques because of its short T2*. New MRI techniques have been recently developed to image bone while avoiding the ionizing radiation present in x-ray-based methods. Qualitative bone imaging can be achieved using ultrashort echo time (UTE), single inversion recovery UTE (IR-UTE), dual-inversion recovery UTE (Dual-IR-UTE), double-inversion recovery UTE (Double-IR-UTE), and zero echo time (ZTE) sequences. The contrast mechanisms as well as the advantages and disadvantages of each technique are discussed.

Volumetric mapping of bound and pore water as well as collagen protons in cortical bone using 3D ultrashort echo time cones MR imaging techniques

Cortical bone assessment using magnetic resonance imaging (MRI) has recently received great attention in an effort to avoid the potential harm associated with ionizing radiation-based techniques. Ultrashort echo time MRI (UTE-MRI) techniques can acquire signal from major hydrogen proton pools in cortical bone, including bound and pore water, as well as from the collagen matrix. This study aimed to develop and evaluate the feasibility of a technique for mapping bound water, pore water, and collagen proton densities in human cortical bone ex vivo and in vivo using three-dimensional UTE Cones (3D-UTE-Cones) MRI. Eight human tibial cortical bone specimens (63 ± 19 years old) were scanned using 3D-UTE-Cones sequences on a clinical 3 T MRI scanner and a micro-computed tomography (μCT) scanner. Total, bound, and pore water proton densities (TWPD, BWPD, and PWPD, respectively) were measured using UTE and inversion recovery UTE (IR-UTE) imaging techniques. Macromolecular proton density (MMPD), a collagen representation, was measured using TWPD and macromolecular fraction (MMF) obtained from two-pool UTE magnetization transfer (UTE-MT) modeling. The correlations between proton densities and μCT-based measures were investigated. The 3D-UTE-Cones techniques were further applied on ten young healthy (34 ± 3 years old) and five old (78 ± 6 years old) female volunteers to evaluate the techniques' feasibility for translational clinical applications. In the ex vivo study, PWPD showed the highest correlations with bone porosity and bone mineral density (BMD) (R = 0.79 and - 0.70, p < 0.01). MMPD demonstrated moderate to strong correlations with bone porosity and BMD (R = -0.67 and 0.65, p < 0.01). MMPD showed strong correlation with age in specimens from female donors (R = -0.91, p = 0.03, n = 5). The presented comprehensive 3D-UTE-Cones imaging protocol allows quantitative mapping of protons in major pools of cortical bone ex vivo and in vivo. PWPD and MMPD can serve as potential novel biomarkers to assess bone matrix and microstructure, as well as bone age- or injury-related variations.

Assessing cortical bone mechanical properties using collagen proton fraction from ultrashort echo time magnetization transfer (UTE-MT) MRI modeling

Cortical bone shows as a signal void when using conventional clinical magnetic resonance imaging (MRI). Ultrashort echo time MRI (UTE-MRI) can acquire high signal from cortical bone, thus enabling quantitative assessments. Magnetization transfer (MT) imaging combined with UTE-MRI can indirectly assess protons in the organic matrix of bone. This study aimed to examine UTE-MT MRI techniques to estimate the mechanical properties of cortical bone. A total of 156 rectangular human cortical bone strips were harvested from the tibial and femoral midshafts of 43 donors (62 ± 22 years old, 62 specimens from females, 94 specimens from males). Bone specimens were scanned using UTE-MT sequences on a clinical 3 T MRI scanner and on a micro-computed tomography (μCT) scanner. A series of MT pulse saturation powers (400°, 600°, 800°) and frequency offsets (2, 5, 10, 20, 50 kHz) was used to measure the macromolecular fraction (MMF) utilizing a two-pool MT model. Failure mechanical properties of the bone specimens were measured using 4-point bending tests. MMF from MRI results showed significant strong correlations with cortical bone porosity (R = -0.72, P < 0.01) and bone mineral density (BMD) (R = +0.71, P < 0.01). MMF demonstrated significant moderate correlations with Young modulus, yield stress, and ultimate stress (R = 0.60-0.61, P < 0.01). These results suggest that the two-pool UTE-MT model focusing on the organic matrix of bone can potentially serve as a novel tool to detect the variations of bone mechanical properties and intracortical porosity.

Collagen proton fraction from ultrashort echo time magnetization transfer (UTE-MT) MRI modelling correlates significantly with cortical bone porosity measured with micro-computed tomography (μCT)

Intracortical bone porosity is a key microstructural parameter that determines bone mechanical properties. While clinical MRI visualizes the cortical bone with a signal void, ultrashort echo time (UTE) MRI can acquire high signal from cortical bone, thus enabling quantitative assessments. Magnetization transfer (MT) imaging combined with UTE-MRI can indirectly assess protons in the bone collagenous matrix, which are inversely related to porosity. This study aimed to examine UTE-MT MRI techniques to evaluate intracortical bone porosity. Eighteen human cortical bone specimens from the tibial and fibular midshafts were scanned using UTE-MT sequences on a clinical 3 T MRI scanner and on a high-resolution micro-computed tomography (μCT) scanner. A series of MT pulse saturation powers (500°, 1000°, 1500°) and frequency offsets (2, 5, 10, 20, 50 kHz) were used to measure the macromolecular fraction (MMF) and macromolecular T2 (T2MM ) using a two-pool MT model. The measurements were made on 136 different regions of interest (ROIs). ROIs were selected at three cortical bone layers (from endosteum to periosteum) and four anatomical sites (anterior, mid-medial, mid-lateral, and posterior) to provide a wide range of porosity. MMF showed moderate to strong correlations with intracortical bone porosity (R = -0.67 to -0.73, p < 0.01) and bone mineral density (BMD) (R = +0.46 to +0.70, p < 0.01). Comparing the average MMF between cortical bone layers revealed a significant increase from the endosteum towards the periosteum. Such a pattern was in agreement with porosity reduction and BMD increase towards the periosteum. These results suggest that the two-pool UTE-MT technique can potentially serve as a novel and accurate tool to assess intracortical bone porosity.

In vivo bone 31 P relaxation times and their implications on mineral quantification

PURPOSE

The intersubject variations in bone phosphorus-31 (³¹ P) T₁ and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> , as well as the implications on in vivo ³¹ P MRI-based bone mineral quantification, were investigated at 3T field strength.

METHODS

A technique that isolates the bone signal from the composite in vivo ³¹ P spectrum was first evaluated via simulation and experiments ex vivo and subsequently applied to measure the T₁ of bone ³¹ P collectively with a spectroscopic saturation recovery sequence in a group of healthy subjects aged 26 to 76 years. <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> was derived from the bone signal linewidth. The density of bone ³¹ P was derived for all subjects from ³¹ P zero TE images acquired in the same scan session using the measured relaxation times. Test-retest experiments were also performed to evaluate repeatability of this in vivo MRI-based bone mineral quantification protocol.

RESULTS

The T₁ obtained in vivo using the proposed spectral separation method combined with saturation recovery sequence is 38.4 ± 1.5 s for the subjects studied. Average ³¹ P density found was 6.40 ± 0.58 mol/L (corresponding to 1072 ± 98 mg/cm³ mineral density), in good agreement with an earlier study in specimens from donors of similar age range. Neither the relaxation times (P = 0.18 for T₁ , P = 0.99 for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> ) nor ³¹ P density (P = 0.55) were found to correlate with subject age. Average coefficients of variation for the repeat study were 1.5%, 2.6%, and 4.4% for bone ³¹ P T₁ , <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> , and density, respectively.

CONCLUSION

Neither ³¹ P T₁ nor <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> varies significantly in healthy adults across a 50-year age range, therefore obviating the need for subject-specific measurements.

Correction: Feasibility of assessing bone matrix and mineral properties in vivo by combined solid-state 1H and 31P MRI

[This corrects the article DOI: 10.1371/journal.pone.0173995.].