In Vivo Quantitative MR Imaging of Bound and Pore Water in Cortical Bone

Ultrashort echo time MRI detects significantly lower collagen but higher pore water in the tibial cortex of female patients with osteopenia and osteoporosis

Ultrashort echo time (UTE) MRI can quantify the major proton pool densities in cortical bone, including total (TWPD), bound (BWPD), and pore water (PWPD) proton densities, as well as the macromolecular proton density (MMPD), associated with the collagen content, which is calculated using macromolecular fraction (MMF) from UTE magnetization transfer (UTE-MT) modeling. This study aimed to investigate the differences in water and collagen contents in tibial cortical bone, between female osteopenia (OPe) patients, osteoporosis (OPo) patients, and young participants (Young). Being postmenopausal and above 55 yr old were the inclusion criteria for OPe and OPo groups. The tibial shaft of 14 OPe (72.5 ± 6.8 yr old), 31 OPo (72.0 ± 6.4 yr old), and 31 young subjects (28.0 ± 6.1 yr old) were scanned using a knee coil on a clinical 3T scanner. Basic UTE, inversion recovery UTE, and UTE-MT sequences were performed. Investigated biomarkers were compared between groups using Kruskal-Wallis test. Spearman's correlation coefficients were calculated between the TH DXA T-score and UTE-MRI results. MMF, BWPD, and MMPD were significantly lower in OPo patients than in the young group, whereas T1, TWPD, and PWPD were significantly higher in OPo patients. The largest OPo/Young average percentage differences were found in MMF (41.9%), PWPD (103.5%), and MMPD (64.0%). PWPD was significantly higher (50.7%), while BWPD was significantly lower (16.4%) in OPe than the Young group on average. MMF was found to be significantly lower (27%) in OPo patients compared with OPe group. T1, MMF, TWPD, PWPD, and MMPD values significantly correlated with the TH DXA T-scores (provided by the patients and only available for OPe and OPo patients). DXA T-score showed the highest correlations with PWPD (R = 0.55) and MMF (R = 0.56) values. TWPD, PWPD, and MMF estimated using the UTE-MRI sequences were recommended to evaluate individuals with OPe and OPo.

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.

Toward the use of MRI measurements of bound and pore water in fracture risk assessment

The current clinical assessment of fracture risk lacks information about the inherent quality of a person's bone tissue. Working toward an imaging-based approach to quantify both a bone tissue quality marker (tissue hydration as water bound to the matrix) and a bone microstructure marker (porosity as water in pores), we hypothesized that the concentrations of bound water (Cbw) are lower and concentrations of pore water (Cpw) are higher in patients with osteoporosis (OP) than in age- and sex-matched adults without the disease. Using recent developments in ultrashort echo time (UTE) magnetic resonance imaging (MRI), maps of Cbw and Cpw were acquired from the uninjured distal third radius (Study 1) of 20 patients who experienced a fragility fracture of the distal radius (Fx) and 20 healthy controls (Non-Fx) and from the tibia mid-diaphysis (Study 2) of 30 women with clinical OP (low T-scores) and 15 women without OP (normal T-scores). In Study 1, Cbw was significantly lower (p = 0.0018) and Cpw was higher (p = 0.0022) in the Fx than in the Non-Fx group. In forward stepwise, logistic regression models using Bayesian Information Criterion for selecting the best set of predictors (from imaging parameters, age, BMI, and DXA scanner type), the area-under-the-receiver operator characteristics-curve (AUC with 95 % confidence intervals) was 0.73 (0.56, 0.86) for hip aBMD (best predictors without MRI) and 0.86 (0.70, 0.95) for the combination of Cbw and Cpw (best predictors overall). In Study 2, Cbw was significantly lower (p = 0.0005) in women with OP (23.8 ± 4.3 1H mol/L) than in women without OP (29.9 ± 6.4 1H mol/L); Cpw was significantly higher by estimate of 2.9 1H mol/L (p = 0.0298) with clinical OP, but only when accounting for the type of UTE-MRI scan with 3D providing higher values than 2D (p < 0.0001). Lastly, Cbw, but not Cpw, was sensitive to bone forming osteoporosis medications over 12-months. UTE-MRI-derived measurements of bound and pore water concentrations are potential, aBMD-independent predictors of fracture risk.

Assessing cortical bone porosity with MRI in an animal model of chronic kidney disease

Chronic kidney disease (CKD) is characterized by secondary hyperparathyroidism and an increased risk of hip fractures predominantly related to cortical porosity. Unfortunately, bone mineral density measurements and high-resolution peripheral computed tomography (HR-pQCT) imaging have shortcomings that limit their utility in these patients. Ultrashort echo time magnetic resonance imaging (UTE-MRI) has the potential to overcome these limitations by providing an alternative assessment of cortical porosity. The goal of the current study was to determine if UTE-MRI could detect changes in porosity in an established rat model of CKD. Cy/+ rats (n = 11), an established animal model of CKD-MBD, and their normal littermates (n = 12) were imaged using microcomputed tomography (microCT) and UTE-MRI at 30 and 35 weeks of age (which approximates late-stage kidney disease in humans). Images were obtained at the distal tibia and the proximal femur. Cortical porosity was assessed using the percent porosity (Pore%) calculated from microCT imaging and the porosity index (PI) calculated from UTE-MRI. Correlations between Pore% and PI were also calculated. Cy/+ rats had higher Pore% than normal rats at both skeletal sites at 35 weeks (tibia = 7.13 % +/- 5.59 % vs. 0.51 % +/- 0.09 %, femur = 19.99 % +/- 7.72 % vs. 2.72 % +/- 0.32 %). They also had greater PI at the distal tibia at 30 weeks of age (0.47 +/- 0.06 vs. 0.40 +/- 0.08). However, Pore% and PI were only correlated in the proximal femur at 35 weeks of age (ρ = 0.929, Spearman). These microCT results are consistent with prior studies in this animal model utilizing microCT imaging. The UTE-MRI results were inconsistent, resulting in variable correlations with microCT imaging, which may be related to suboptimal bound and pore water discrimination at higher magnetic field strengths. Nevertheless, UTE-MRI may still provide an additional clinical tool to assess fracture risk without using ionizing radiation in CKD patients.

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.

T1 relaxation of bound and pore water in cortical bone

MRI measures of bound and/or pore water concentration in cortical bone offer potential diagnostics of bone fracture risk. The transverse relaxation characteristics of both bound and pore water are relatively well understood and have been used to design clinical MRI pulse sequences to image each water pool quantitatively. However, these methods are also sensitive to longitudinal relaxation characteristics, which have been less well studied. Here, spectroscopic relaxometry measurements of 31 human cortical bone specimens provided a more detailed picture of T 1 of both bound and pore water. The results included mean, standard deviation, and range of T 1 spectra from both bound and pore water, as well as novel presentations of the 2D T 1 - T 2 distribution of pore water. Importantly, for each sample the pore water T 1 spectrum was found to span more than one order of magnitude and varied substantially across the 31 sample studies. Because many existing methods assume pore water T 1 to be mono-exponential and constant across individuals, the results were used to compute the potential effect neglecting this intra- and intersample T 1 variation on accurate MRI measurement of both bound and pore water concentrations. The greatest effect was found for adiabatic inversion recovery (AIR) based measurements of bound water concentration, which showed an average of 8.8% and as much as 37% error when using a common mono-exponential assumption of pore water T 1 . Despite these errors, the simulated AIR measurements were still moderately well correlated with the bound water concentrations derived from the spectroscopic data.

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.

Fast bound and pore water mapping of cortical bone with arbitrary slice oriented two-dimensional ultra-short echo time

PURPOSE

Extend fast, two-dimensional (2D) methods of bound and pore water mapping in bone to arbitrary slice orientation.

METHODS

To correct for slice profile artifacts caused by gradient errors of half pulse 2D ultra-short echo time (UTE), we developed a library of predistorted gradient waveforms that can be used to interpolate optimized gradient waveforms for 2D UTE slice selection. We also developed a method to estimate and correct for a bulk phase difference between the two half pulse excitations used for 2D UTE signal excitation. Bound water images were acquired in three healthy subjects with adiabatic inversion recovery prepared 2D UTE, while pore water images were acquired after short-T2 signals were suppressed with double adiabatic inversion recovery preparation. The repeatability of bound and pore water imaging with 2D UTE was tested by repeating acquisitions after repositioning.

RESULTS

The library-based interpolation of optimized slice select gradient waveforms combined with the method to estimate bulk phase between two excitations provided compact slice profiles for half pulse excited 2D UTE. Quantitative bound and pore water values were highly repeatable-the pooled SD of bound water across all three subjects was 0.38 mol1 $$ {}^1 $$ H/L, while pooled SD of pore water was 0.30 mol1 $$ {}^1 $$ H/L.

CONCLUSION

Fast, quantitative, 2D UTE-based bound and pore water images can be acquired at arbitrary oblique orientations after correcting for errors in the slice select gradient waveform and bulk phase shift between the two half acquisitions.

Synthetic CT in Musculoskeletal Disorders: A Systematic Review

ABSTRACT

Repeated computed tomography (CT) examinations increase patients' ionizing radiation exposure and health costs, making an alternative method desirable. Cortical and trabecular bone, however, have short T2 relaxation times, causing low signal intensity on conventional magnetic resonance (MR) sequences. Different techniques are available to create a "CT-like" contrast of bone, such as ultrashort echo time, zero echo time, gradient-echo, and susceptibility-weighted image MR sequences, and artificial intelligence. This systematic review summarizes the essential technical background and developments of ultrashort echo time, zero echo time, gradient-echo, susceptibility-weighted image MR imaging sequences and artificial intelligence; presents studies on research and clinical applications of "CT-like" MR imaging; and describes their main advantages and limitations. We also discuss future opportunities in research, which patients would benefit the most, the most appropriate situations for using the technique, and the potential to replace CT in the clinical workflow.

Comprehensive assessment of osteoporosis in lumbar spine using compositional MR imaging of trabecular bone

OBJECTIVES

To comprehensively assess osteoporosis in the lumbar spine, a compositional MR imaging technique is proposed to quantify proton fractions for all the water components as well as fat in lumbar vertebrae measured by a combination of a 3D short repetition time adiabatic inversion recovery prepared ultrashort echo time (STAIR-UTE) MRI and IDEAL-IQ.

METHODS

A total of 182 participants underwent MRI, quantitative CT, and DXA. Lumbar collagen-bound water proton fraction (CBWPF), free water proton fraction (FWPF), total water proton fraction (TWPF), bone mineral density (BMD), and T-score were calculated in three vertebrae (L2-L4) for each subject. The correlations of the CBWPF, FWPF, and TWPF with BMD and T-score were investigated respectively. A comprehensive diagnostic model combining all the water components and clinical characteristics was established. The performances of all the water components and the comprehensive diagnostic model to discriminate between normal, osteopenia, and osteoporosis cohorts were also evaluated using receiver operator characteristic (ROC).

RESULTS

The CBWPF showed strong correlations with BMD (r = 0.85, p < 0.001) and T-score (r = 0.72, p < 0.001), while the FWPF and TWPF showed moderate correlations with BMD (r = 0.65 and 0.68, p < 0.001) and T-score (r = 0.47 and 0.49, p < 0.001). The high area under the curve values obtained from ROC analysis demonstrated that CBWPF, FWPF, and TWPF have the potential to differentiate the normal, osteopenia, and osteoporosis cohorts. At the same time, the comprehensive diagnostic model shows the best performance.

CONCLUSIONS

The compositional MRI technique, which quantifies CBWPF, FWPF, and TWPF in trabecular bone, is promising in the assessment of bone quality.

KEY POINTS

• Compositional MR imaging technique is able to quantify proton fractions for all the water components (i.e., collagen-bound water proton fraction (CBWPF), free water proton fraction (FWPF), and total water proton fraction (TWPF)) in the human lumbar spine. • The biomarkers derived from the compositional MR imaging technique showed moderate to high correlations with bone mineral density (BMD) and T-score and showed good performance in distinguishing people with different bone mass. • The comprehensive diagnostic model incorporating CBWPF, FWPF, TWPF, and clinical characteristics showed the highest clinical diagnostic capability for the assessment of osteoporosis.

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.

MRI-based mechanical competence assessment of bone using micro finite element analysis (micro-FEA): Review

Areal bone mineral density (aBMD) from dual-energy x-ray absorptiometry (DEXA) and volumetric bone mineral density (vBMD) have demonstrated limited capabilities in the evaluation of bone mechanical competence and prediction of bone fracture. Predicting the macroscopic mechanical behavior of the bone structure has been challenging because of the heterogeneous and anisotropic nature of bone, such as the dependencies on loading direction, anatomical location, and sample dimensions. Magnetic resonance imaging (MRI) has been introduced as a promising modality that can be coupled with finite element analysis (FEA) for the assessment of bone mechanical competence. This review article describes studies investigating MRI-based micro-FEA as a potential non-invasive method to predict bone mechanical competence and facilitate bone fracture risk estimation without exposure to ionizing radiation. Specifically, the steps, applications, and future potential of FEA using indirect and direct bone imaging are discussed.

High-contrast osteochondral junction imaging using a 3D dual adiabatic inversion recovery-prepared ultrashort echo time cones sequence

While conventional MRI sequences cannot visualize tissues from the osteochondral junction (OCJ) due to these tissues' short transverse T₂ /T₂ * relaxations, ultrashort echo time (UTE) sequences can overcome this limitation. A 2D UTE sequence with a dual adiabatic inversion recovery preparation (DIR-UTE) for selective imaging of short T₂ tissues with high contrast has previously been developed, but high sensitivity to eddy currents and aliased out-of-slice excitation make it difficult to image the thin layer of the OCJ in vivo. Here, we combine the DIR scheme with a 3D UTE cones sequence for volumetric imaging of OCJ tissues in vivo, aiming to generate higher OCJ contrast compared with a recently developed single IR-prepared UTE sequence with a fat saturation module (IR-FS-UTE). All sequences were implemented on a 3-T clinical scanner. The DIR-UTE cones sequence combined a 3D UTE cones sequence with two narrow-band adiabatic IR preparation pulses centered on water and fat spectrum frequencies, respectively. The 3D DIR-UTE cones sequence was first applied to a phantom, then to the knees of four healthy volunteers and four patients diagnosed with osteoarthritis and compared with the IR-FS-UTE sequence. In both phantom and volunteer studies, the proposed DIR-UTE cones sequence showed much higher contrast for OCJ imaging than the IR-FS-UTE sequence did. The 3D DIR-UTE cones sequence showed a significantly higher contrast-to-noise ratio between the OCJ and subchondral bone fat (mean, standard deviation [SD]: 25.7 ± 2.3) and between the OCJ and superficial layers of cartilage (mean, SD: 22.2 ± 3.5) compared with the IR-FS-UTE sequence (mean, SD: 10.8 ± 2.5 and 16.3 ± 2.6, respectively). The 3D DIR-UTE cones sequence is feasible for imaging of the OCJ region of the knee in vivo and produces both high resolution and high contrast.

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.

Quantifying cortical bone free water using short echo time (STE-MRI) at 1.5 T

PURPOSE

The purpose of our study was to use Dual-TR STE-MR protocol as a clinical tool for cortical bone free water quantification at 1.5 T and validate it by comparing the obtained results (MR-derived results) with dehydration results.

METHODS

Human studies were compliant with HIPPA and were approved by the institutional review board. Short Echo Time (STE) MR imaging with different Repetition Times (TRs) was used for quantification of cortical bone free water T₁ (T₁free) and concentration (ρ_free). The proposed strategy was compared with the dehydration technique in seven bovine cortical bone samples. The agreement between the two methods was quantified by using Bland and Altman analysis. Then we applied the technique on a cross-sectional population of thirty healthy volunteers (18F/12M) and examined the association of the biomarkers with age.

RESULTS

The mean values of ρ_free for bovine cortical bone specimens were quantified as 4.37% and 5.34% by using STE-MR and dehydration techniques, respectively. The Bland and Altman analysis showed good agreement between the two methods along with the suggestion of 0.99% bias between them. Strong correlations were also reported between ρ_free (r² = 0.62) and T_1free and age (r² = 0.8). The reproducibility of the method, evaluated in eight subjects, yielded an intra-class correlation of 0.95.

CONCLUSION

STE-MR imaging with dual-TR strategy is a clinical solution for quantifying cortical bone ρ_free and T_1free.

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.

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.

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.

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.

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.

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.

Applications of PET-MRI in musculoskeletal disease

New integrated PET-MRI systems potentially provide a complete imaging modality for diagnosis and evaluation of musculoskeletal disease. MRI is able to provide excellent high-resolution morphologic information with multiple contrast mechanisms that has made it the imaging modality of choice in evaluation of many musculoskeletal disorders. PET offers incomparable abilities to provide quantitative information about molecular and physiologic changes that often precede structural and biochemical changes. In combination, hybrid PET-MRI can enhance imaging of musculoskeletal disorders through early detection of disease as well as improved diagnostic sensitivity and specificity. The purpose of this article is to review emerging applications of PET-MRI in musculoskeletal disease. Both clinical applications of malignant musculoskeletal disease as well as new opportunities to incorporate the molecular capabilities of nuclear imaging into studies of nononcologic musculoskeletal disease are discussed. Lastly, we discuss some of the technical considerations and challenges of PET-MRI as they specifically relate to musculoskeletal disease.

LEVEL OF EVIDENCE

5 TECHNICAL EFFICACY: Stage 3 J. Magn. Reson. Imaging 2018;48:27-47.

Imaging of the region of the osteochondral junction (OCJ) using a 3D adiabatic inversion recovery prepared ultrashort echo time cones (3D IR-UTE-cones) sequence at 3 T

The purpose of this study is to develop a 3D adiabatic inversion recovery prepared ultrashort echo time Cones (3D IR-UTE-Cones) sequence for high resolution and contrast imaging of the region of osteochondral junction (OCJ) of human knee joint using a clinical 3 T scanner. A feasibility study on direct imaging of the OCJ region was performed on a human patellar cartilage sample and on eight cadaveric knee joints using T₁ -weighted, proton density (PD)-weighted and short-T₂ -weighted 3D IR-UTE-Cones sequences. Contrast to noise ratio was measured to evaluate the effectiveness of the 3D IR-UTE-Cones sequences for selective imaging of the OCJ region. Computed tomography imaging was performed in parallel for the cadaveric knee joints. The optimized T₁ -weighted 3D IR-UTE-Cones sequence was used to image the knee joints of eight healthy volunteers and six patients with osteoarthritis (OA) to evaluate morphological changes in the OCJ region. Clinical PD- and T₂ -weighted FSE sequences were also performed for comparison. The T₁ -weighted 3D IR-UTE-Cones sequence showed high resolution and contrast bright band of the normal OCJ region in the cadaveric joints. Normal OCJ appearances were also seen in healthy volunteers. Abnormal OCJ regions, manifested as ill-defined, focal loss or non-visualization of the high intensity band adjacent to the subchondral bone plate, were observed in the knee joints of both ex vivo and in vivo OA patients. The 3D IR-UTE-Cones sequence can image OCJ regions ex vivo and in vivo, with abnormalities depicted with high resolution and contrast. The technique may be useful for demonstrating involvement of OCJ regions in early OA.

Quantitative MRI Musculoskeletal Techniques: An Update

OBJECTIVE. For many years, MRI of the musculoskeletal system has relied mostly on conventional sequences with qualitative analysis. More recently, using quantitative MRI applications to complement qualitative imaging has gained increasing interest in the MRI community, providing more detailed physiologic or anatomic information. CONCLUSION. In this article, we review the current state of quantitative MRI, technical and software advances, and the most relevant clinical and research musculoskeletal applications of quantitative MRI.

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.

MRI assessment of bone structure and microarchitecture

Osteoporosis is a disease of weak bone and increased fracture risk caused by low bone mass and microarchitectural deterioration of bone tissue. The standard-of-care test used to diagnose osteoporosis, dual-energy x-ray absorptiometry (DXA) estimation of areal bone mineral density (BMD), has limitations as a tool to identify patients at risk for fracture and as a tool to monitor therapy response. Magnetic resonance imaging (MRI) assessment of bone structure and microarchitecture has been proposed as another method to assess bone quality and fracture risk in vivo. MRI is advantageous because it is noninvasive, does not require ionizing radiation, and can evaluate both cortical and trabecular bone. In this review article, we summarize and discuss research progress on MRI of bone structure and microarchitecture over the last decade, focusing on in vivo translational studies. Single-center, in vivo studies have provided some evidence for the added value of MRI as a biomarker of fracture risk or treatment response. Larger, prospective, multicenter studies are needed in the future to validate the results of these initial translational studies.

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 5 J. MAGN. RESON. IMAGING 2017;46:323-337.

Three-dimensional adiabatic inversion recovery prepared ultrashort echo time cones (3D IR-UTE-Cones) imaging of cortical bone in the hip

PURPOSE

We present three-dimensional adiabatic inversion recovery prepared ultrashort echo time Cones (3D IR-UTE-Cones) imaging of cortical bone in the hip of healthy volunteers using a clinical 3T scanner.

METHODS

A 3D IR-UTE-Cones sequence, based on a short pulse excitation followed by a 3D Cones trajectory, with a nominal TE of 32μs, was employed for high contrast morphological imaging of cortical bone in the hip of heathy volunteers. Signals from soft tissues such as muscle and marrow fat were suppressed via adiabatic inversion and signal nulling. T₂^⁎ value of the cortical bone was also calculated based on 3D IR-UTE-Cones acquisitions with a series of TEs ranging from 0.032 to 0.8ms. A total of four healthy volunteers were recruited for this study. Average T₂^⁎ values and the standard deviation for four regions of interests (ROIs) at the greater trochanter, the femoral neck, the femoral head and the lesser trochanter were calculated.

RESULTS

The 3D IR-UTE-Cones sequence provided efficient suppression of soft tissues with excellent image contrast for cortical bone visualization in all volunteer hips. Exponential single component decay was observed for all ROIs, with averaged T₂^⁎ values ranging from 0.33 to 0.45ms, largely consistent with previously reported T₂^⁎ values of cortical bone in the tibial midshaft.

CONCLUSIONS

The 3D IR-UTE-Cones sequence allows in vivo volumetric imaging and quantitative T₂^⁎ measurement of cortical bone in the hip using a clinical 3T scanner.

Measurement of bound and pore water T1 relaxation times in cortical bone using three-dimensional ultrashort echo time cones sequences

PURPOSE

We present three-dimensional ultrashort echo time Cones (3D UTE Cones) techniques for quantification of total water T1 ( T1TW), bound water T1 ( T1BW), and pore water T1 ( T1PW) in vitro and in vivo using a 3 Tesla (T) scanner.

METHODS

T1TW, T1BW, and T1PW were measured with three-dimensional (3D) Cones and adiabatic inversion recovery Cone (IR-Cone) sequences. Two-dimensional (2D) nonselective ultrashort echo time (UTE) techniques, including saturation recovery, variable repetition times (TRs), and inversion recovery (IR) preparation approaches were compared with 3D-Cones techniques on bovine cortical bone samples (n = 8). The 3D Cones sequences were used to measure T1TW, T1BW, and T1PW in the tibial midshaft of healthy volunteers (n = 8).

RESULTS

Comparable T1 images were achieved for cortical bone between 3D Cones and 2D UTE techniques as well as those published in the literature. The 3D Cones sequences showed a mean T1TW of 208 ± 22 ms, a mean T1PW of 545 ± 28 ms, and a mean T1BW of 131 ± 12 ms for bovine cortical bone; and a mean T1TW of 246 ± 32 ms, a mean T1PW of 524 ± 46 ms, and a mean T1BW of 134 ± 11 ms for the tibial midshaft of healthy volunteers.

CONCLUSIONS

The 3D Cones sequences can be used for fast volumetric assessment of bound and pore water T1 images in vitro and in vivo. Magn Reson Med 77:2136-2145, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Bilateral Changes in Deep Tissue Environment After Manual Lymphatic Drainage in Patients with Breast Cancer Treatment-Related Lymphedema

BACKGROUND

Breast cancer treatment-related lymphedema (BCRL) arises from a mechanical insufficiency following cancer therapies. Early BCRL detection and personalized intervention require an improved understanding of the physiological processes that initiate lymphatic impairment. Here, internal magnetic resonance imaging (MRI) measures of the tissue microenvironment were paired with clinical measures of tissue structure to test fundamental hypotheses regarding structural tissue and muscle changes after the commonly used therapeutic intervention of manual lymphatic drainage (MLD).

METHODS AND RESULTS

Measurements to identify lymphatic dysfunction in healthy volunteers (n = 29) and patients with BCRL (n = 16) consisted of (1) limb volume, tissue dielectric constant, and bioelectrical impedance (i.e., non-MRI measures); (2) qualitative 3 Tesla diffusion-weighted, T₁-weighted and T₂-weighted MRI; and (3) quantitative multi-echo T₂ MRI of the axilla. Measurements were repeated in patients immediately following MLD. Normative control and BCRL T₂ values were quantified and a signed Wilcoxon Rank-Sum test was applied (significance: two-sided p < 0.05). Non-MRI measures yielded significant capacity for discriminating between arms with versus without clinical signs of BCRL, yet yielded no change in response to MLD. Alternatively, a significant increase in deep tissue T₂ on the involved (pre T₂ = 0.0371 ± 0.003 seconds; post T₂ = 0.0389 ± 0.003; p = 0.029) and contralateral (pre T₂ = 0.0365 ± 0.002; post T₂ = 0.0395 ± 0.002; p < 0.01) arms was observed. Trends for larger T₂ increases on the involved side after MLD in patients with stage 2 BCRL relative to earlier stages 0 and 1 BCRL were observed, consistent with tissue composition changes in later stages of BCRL manifesting as breakdown of fibrotic tissue after MLD in the involved arm. Contrast consistent with relocation of fluid to the contralateral quadrant was observed in all stages.

CONCLUSION

Quantitative deep tissue T₂ MRI values yielded significant changes following MLD treatment, whereas non-MRI measurements did not vary. These findings highlight that internal imaging measures of tissue composition may be useful for evaluating how current and emerging therapies impact tissue function.

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

Purpose

To develop and evaluate an integrated imaging protocol for bone water and phosphorus quantification in vivo by solid-state ¹H and ³¹P MRI.

Materials and methods

All studies were HIPAA-compliant and were performed with institutional review board approval and written informed consent. Proton (¹H) ultra-short echo-time (UTE) and phosphorus (³¹P) zero echo-time (ZTE) sequences were designed and implemented on a 3 T clinical MR scanner to quantify bone water and mineral in vivo. The left tibia of ten healthy subjects (including both genders, 49±15 y/o) was examined with a custom-built ¹H/³¹P dual-frequency extremity RF coil. Total bone water (TW), water bound to the collagen matrix (BW) and bone ³¹P were quantified from MR images with respect to reference samples of known ¹H or ³¹P concentration, and pore water (PW) was subsequently determined from TW and BW. Porosity index (PI) was calculated as the ratio between UTE images acquired at two echo times. MRI parameters were compared with bone density measures obtained by high-resolution peripheral quantitative CT (HR-pQCT).

Results

The total scan time for the bone water and ³¹P quantification protocol was about 50 minutes. Average TW, BW, PW and ³¹P concentrations were 13.99±1.26, 10.39±0.80, 3.34±1.41 mol/L and 7.06±1.53 mol/L for the studied cohort, respectively, in good agreement with previous results conducted ex vivo. Average intra-subject coefficients of variation were 3.47%, 2.60% and 7.50% for TW, BW and PW and 5.60% for ³¹P. Negative correlations were observed between PW and vBMD (p<0.05) as well as between PI and ³¹P (p<0.05), while bone mineral content (BMC) estimated from ³¹P MRI and HR-pQCT were strongly positively correlated (p<0.0001).

Conclusion

This work demonstrates the feasibility of quantifying bone water and mineral phosphorus in human subjects in a single MRI session with a clinically practical imaging protocol.

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.

30-Second bound and pore water concentration mapping of cortical bone using 2D UTE with optimized half-pulses

PURPOSE

MRI of cortical bone has the potential to offer new information about fracture risk. Current methods are typically performed with 3D acquisitions, which suffer from long scan times and are generally limited to extremities. This work proposes using 2D UTE with half pulses for quantitatively mapping bound and pore water in cortical bone.

METHODS

Half-pulse 2D UTE methods were implemented on a 3T Philips Achieva scanner using an optimized slice-select gradient waveform, with preparation pulses to selectively image bound or pore water. The 2D methods were quantitatively compared with previously implemented 3D methods in the tibia in five volunteers.

RESULTS

The mean difference between bound and pore water concentration acquired from 3D and 2D sequences was 0.6 and 0.9 mol ¹ H/L_bone (3 and 12%, respectively). While 2D pore water methods tended to slightly overestimate concentrations relative to 3D methods, differences were less than scan-rescan uncertainty and expected differences between healthy and fracture-prone bones.

CONCLUSION

Quantitative bound and pore water concentration mapping in cortical bone can be accelerated by 2 orders of magnitude using 2D protocols with optimized half-pulse excitation. Magn Reson Med 77:945-950, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Potential of PET-MRI for imaging of non-oncologic musculoskeletal disease

Early detection of musculoskeletal disease leads to improved therapies and patient outcomes, and would benefit greatly from imaging at the cellular and molecular level. As it becomes clear that assessment of multiple tissues and functional processes are often necessary to study the complex pathogenesis of musculoskeletal disorders, the role of multi-modality molecular imaging becomes increasingly important. New positron emission tomography-magnetic resonance imaging (PET-MRI) systems offer to combine high-resolution MRI with simultaneous molecular information from PET to study the multifaceted processes involved in numerous musculoskeletal disorders. In this article, we aim to outline the potential clinical utility of hybrid PET-MRI to these non-oncologic musculoskeletal diseases. We summarize current applications of PET molecular imaging in osteoarthritis (OA), rheumatoid arthritis (RA), metabolic bone diseases and neuropathic peripheral pain. Advanced MRI approaches that reveal biochemical and functional information offer complementary assessment in soft tissues. Additionally, we discuss technical considerations for hybrid PET-MR imaging including MR attenuation correction, workflow, radiation dose, and quantification.

Fast volumetric imaging of bound and pore water in cortical bone using three-dimensional ultrashort-TE (UTE) and inversion recovery UTE sequences

We report the three-dimensional ultrashort-TE (3D UTE) and adiabatic inversion recovery UTE (IR-UTE) sequences employing a radial trajectory with conical view ordering for bi-component T2 * analysis of bound water (T2 *(BW) ) and pore water (T2 *(PW) ) in cortical bone. An interleaved dual-echo 3D UTE acquisition scheme was developed for fast bi-component analysis of bound and pore water in cortical bone. A 3D IR-UTE acquisition scheme employing multiple spokes per IR was developed for bound water imaging. Two-dimensional UTE (2D UTE) and IR-UTE sequences were employed for comparison. The sequences were applied to bovine bone samples (n = 6) and volunteers (n = 6) using a 3-T scanner. Bi-component fitting of 3D UTE images of bovine samples showed a mean T2 *(BW) of 0.26 ± 0.04 ms and T2 *(PW) of 4.16 ± 0.35 ms, with fractions of 21.5 ± 3.6% and 78.5 ± 3.6%, respectively. The 3D IR-UTE signal showed a single-component decay with a mean T2 *(BW) of 0.29 ± 0.05 ms, suggesting selective imaging of bound water. Similar results were achieved with the 2D UTE and IR-UTE sequences. Bi-component fitting of 3D UTE images of the tibial midshafts of healthy volunteers showed a mean T2 *(BW) of 0.32 ± 0.08 ms and T2 *(PW) of 5.78 ± 1.24 ms, with fractions of 34.2 ± 7.4% and 65.8 ± 7.4%, respectively. Single-component fitting of 3D IR-UTE images showed a mean T2 *(BW) of 0.35 ± 0.09 ms. The 3D UTE and 3D IR-UTE techniques allow fast volumetric mapping of bound and pore water in cortical bone. Copyright © 2016 John Wiley & Sons, Ltd.

Magnetization transfer imaging of cortical bone in vivo using a zero echo time sequence in mice at 4.7 T: a feasibility study

OBJECTIVE

To investigate the feasibility of magnetization transfer (MT) imaging in mice in vivo for the assessment of cortical bone.

MATERIALS AND METHODS

MT-zero echo time data were acquired at 4.7 T in six mice using MT preparation pulses with two different flip angles (FAs) and a series of ten different off-resonance frequencies (500-15000 Hz). Regions of interest were drawn at multiple levels of the femoral cortical bone. The MT ratio (MTR) was computed for each combination of FAs and off-resonance frequencies. T1 measurements were used to estimate the direct saturation (DS) using a Bloch equation simulation. Estimation of the absorption line width of cortical bone from T2* measurements was also performed.

RESULTS

MTR values were higher using 3000° FA than 1000° FA. MTR values decreased toward higher off-resonance frequencies. Maximum mean MTR ± standard deviation (SD) of 58.57 ± 5.22 (range 50.44-70.61) was measured with a preparation pulse of 3000° and off-resonance frequency of 500 Hz. Maximum "true" MT effect was estimated at around 2-3 and 5 kHz, respectively, for 1000° and 3000° FA. Mean full width at half maximum ± SD of 577 ± 91 Hz was calculated for the absorption spectral line of the cortical bone.

CONCLUSION

MT imaging can be used for the assessment of cortical bone in mice in vivo. DS effects are negligible using preparation pulses with off-resonance frequencies greater than 3 kHz.

Solid-State Quantitative (1)H and (31)P MRI of Cortical Bone in Humans

Magnetic resonance imaging (MRI) plays a pivotal role for assessment of the musculoskeletal system. It is currently the clinical modality of choice for evaluation of soft tissues including cartilage, ligaments, tendons, muscle, and bone marrow. By comparison, the study of calcified tissue by MRI is still in its infancy. In this article, we review the potential of the modality for assessment of cortical bone properties known to be affected in degenerative bone disease, with focus on parameters related to matrix and mineral densities, and porosity, by means of emerging solid-state (1)H and (31)P MRI techniques. In contrast to soft tissues, the MRI signal in calcified tissues has very short lifetime, on the order of 100 μs to a few milliseconds, demanding customized imaging approaches that allow capture of the signal almost immediately after excitation. The technologies described are suited for quantitatively imaging human cortical bone in specimens as well as in vivo in patients on standard clinical imagers, yielding either concentrations in absolute units when measured against a reference standard, or more simply, in the form of surrogate biomarkers. The two major water fractions in cortical bone are those of collagen-bound and pore water occurring at an approximately 3:1 ratio. Collagen-bound water density provides a direct quantitative measure of osteoid density. While at an earlier stage of development, quantification of mineral phosphorus by (31)P MRI yields mineral density and, together with knowledge of matrix density, should allow quantification of the degree of bone mineralization.

MRI-derived bound and pore water concentrations as predictors of fracture resistance

Accurately predicting fracture risk in the clinic is challenging because the determinants are multi-factorial. A common approach to fracture risk assessment is to combine X-ray-based imaging methods such as dual-energy X-ray absorptiometry (DXA) with an online Fracture Risk Assessment Tool (FRAX) that includes additional risk factors such as age, family history, and prior fracture incidents. This approach still does not adequately diagnose many individuals at risk, especially those with certain diseases like type 2 diabetes. As such, this study investigated bound water and pore water concentrations (Cbw and Cpw) from ultra-short echo time (UTE) magnetic resonance imaging (MRI) as new predictors of fracture risk. Ex vivo cadaveric arms were imaged with UTE MRI as well as with DXA and high-resolution micro-computed tomography (μCT), and imaging measures were compared to both whole-bone structural and material properties as determined by three-point bending tests of the distal-third radius. While DXA-derived areal bone mineral density (aBMD) and μCT-derived volumetric BMD correlated well with structural strength, they moderately correlated with the estimate material strength with gender being a significant covariate for aBMD. MRI-derived measures of Cbw and Cpw had a similar predictive ability of material strength as aBMD but did so independently of gender. In addition, Cbw was the only imaging parameter to significantly correlate with toughness, the energy dissipated during fracture. Notably, the strength of the correlations with the material properties of bone tended to be higher when a larger endosteal region was used to determine Cbw and Cpw. These results indicate that MRI measures of Cbw and Cpw have the ability to probe bone material properties independent of bone structure or subject gender. In particular, toughness is a property of fracture resistance that is not explained by X-ray based methods. Thus, these MRI-derived measures of Cbw and Cpw in cortical bone have the potential to be useful in clinical populations for evaluating fracture risk, especially involving diseases that affect material properties of the bone beyond its strength.

Reference point indentation is insufficient for detecting alterations in traditional mechanical properties of bone under common experimental conditions

Reference point indentation (RPI) was developed as a novel method to assess mechanical properties of bone in vivo, yet it remains unclear what aspects of bone dictate changes/differences in RPI-based parameters. The main RPI parameter, indentation distance increase (IDI), has been proposed to be inversely related to the ability of bone to form/tolerate damage. The goal of this work was to explore the relationshipre-intervention RPI measurebetween RPI parameters and traditional mechanical properties under varying experimental conditions (drying and ashing bones to increase brittleness, demineralizing bones and soaking in raloxifene to decrease brittleness). Beams were machined from cadaveric bone, pre-tested with RPI, subjected to experimental manipulation, post-tested with RPI, and then subjected to four-point bending to failure. Drying and ashing significantly reduced RPI's IDI, as well as ultimate load (UL), and energy absorption measured from bending tests. Demineralization increased IDI with minimal change to bending properties. Ex vivo soaking in raloxifene had no effect on IDI but tended to enhance post-yield behavior at the structural level. These data challenge the paradigm of an inverse relationship between IDI and bone toughness, both through correlation analyses and in the individual experiments where divergent patterns of altered IDI and mechanical properties were noted. Based on these results, we conclude that RPI measurements alone, as compared to bending tests, are insufficient to reach conclusions regarding mechanical properties of bone. This proves problematic for the potential clinical use of RPI measurements in determining fracture risk for a single patient, as it is not currently clear that there is an IDI, or even a trend of IDI, that can determine clinically relevant changes in tissue properties that may contribute to whole bone fracture resistance.

Age-related changes in the fracture resistance of male Fischer F344 rat bone

In addition to the loss in bone volume that occurs with age, there is a decline in material properties. To test new therapies or diagnostic tools that target such properties as material strength and toughness, a pre-clinical model of aging would be useful in which changes in bone are similar to those that occur with aging in humans. Toward that end, we hypothesized that similar to human bone, the estimated toughness and material strength of cortical bone at the apparent-level decreases with age in the male Fischer F344 rat. In addition, we tested whether the known decline in trabecular architecture in rats translated to an age-related decrease in vertebra (VB) strength and whether non-X-ray techniques could quantify tissue changes at micron and sub-micron length scales. Bones were harvested from 6-, 12-, and 24-month (mo.) old rats (n=12 per age). Despite a loss in trabecular bone with age, VB compressive strength was similar among the age groups. Similarly, whole-bone strength (peak force) in bending was maintained (femur) or increased (radius) with aging. There was though an age-related decrease in post-yield toughness (radius) and bending strength (femur). The ability to resist crack initiation was actually higher for the 12-mo. and 24-mo. than for 6-mo. rats (notch femur), but the estimated work to propagate the crack was less for the aged bone. For the femur diaphysis region, porosity increased while bound water decreased with age. For the radius diaphysis, there was an age-related increase in non-enzymatic and mature enzymatic collagen crosslinks. Raman spectroscopy analysis of embedded cross-sections of the tibia mid-shaft detected an increase in carbonate subsitution with advanced aging for both inner and outer tissue.

Evaluation of bound and pore water in cortical bone using ultrashort-TE MRI

Bone water exists in different states with the majority bound to the organic matrix and to mineral, and a smaller fraction in 'free' form in the pores of cortical bone. In this study, we aimed to develop and evaluate ultrashort-TE (UTE) MRI techniques for the assessment of T2*, T1 and concentration of collagen-bound and pore water in cortical bone using a 3-T clinical whole-body scanner. UTE MRI, together with an isotope study using tritiated and distilled water (THO-H2O) exchange, as well as gravimetric analysis, were performed on ten sectioned bovine bone samples. In addition, 32 human cortical bone samples were prepared for comparison between the pore water concentration measured with UTE MRI and the cortical porosity derived from micro-computed tomography (μCT). A short T2* of 0.27 ± 0.03 ms and T1 of 116 ± 6 ms were observed for collagen-bound water in bovine bone. A longer T2* of 1.84 ± 0.52 ms and T1 of 527 ± 28 ms were observed for pore water in bovine bone. UTE MRI measurements showed a pore water concentration of 4.7-5.3% by volume and collagen-bound water concentration of 15.7-17.9% in bovine bone. THO-H2O exchange studies showed a pore water concentration of 5.9 ± 0.6% and collagen-bound water concentration of 18.1 ± 2.1% in bovine bone. Gravimetric analysis showed a pore water concentration of 6.3 ± 0.8% and collagen-bound water concentration of 19.2 ± 3.6% in bovine bone. A mineral water concentration of 9.5 ± 0.6% was derived in bovine bone with the THO-H2O exchange study. UTE-measured pore water concentration is highly correlated (R(2) = 0.72, p < 0.0001) with μCT porosity in the human cortical bone study. Both bovine and human bone studies suggest that UTE sequences could reliably measure collagen-bound and pore water concentration in cortical bone using a clinical scanner.