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

Towards assessing and improving the reliability of ultrashort echo time quantitative magnetization transfer (UTE-qMT) MRI of cortical bone: In silico and ex vivo study

OBJECTIVE

To assess and improve the reliability of the ultrashort echo time quantitative magnetization transfer (UTE-qMT) modeling of the cortical bone.

MATERIALS AND METHODS

Simulation-based digital phantoms were created that mimic the UTE-qMT properties of cortical bones. A wide range of SNR from 25 to 200 was simulated by adding different levels of noise to the synthesized MT-weighted images to assess the effect of SNR on UTE-qMT fitting results. Tensor-based denoising algorithm was applied to improve the fitting results. These results from digital phantom studies were validated via ex vivo rat leg bone scans.

RESULTS

The selection of initial points for nonlinear fitting and the number of data points tested for qMT analysis have minimal effect on the fitting result. Magnetization exchange rate measurements are highly dependent on the SNR of raw images, which can be substantially improved with an appropriate denoising algorithm that gives similar fitting results from the raw images with an 8-fold higher SNR.

DISCUSSION

The digital phantom approach enables the assessment of the reliability of bone UTE-qMT fitting by providing the known ground truth. These findings can be utilized for optimizing the data acquisition and analysis pipeline for UTE-qMT imaging of cortical bones.

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.

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.

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.

Making the invisible visible-ultrashort echo time magnetic resonance imaging: Technical developments and applications

Magnetic resonance imaging (MRI) uses a large magnetic field and radio waves to generate images of tissues in the body. Conventional MRI techniques have been developed to image and quantify tissues and fluids with long transverse relaxation times (T2s), such as muscle, cartilage, liver, white matter, gray matter, spinal cord, and cerebrospinal fluid. However, the body also contains many tissues and tissue components such as the osteochondral junction, menisci, ligaments, tendons, bone, lung parenchyma, and myelin, which have short or ultrashort T2s. After radio frequency excitation, their transverse magnetizations typically decay to zero or near zero before the receiving mode is enabled for spatial encoding with conventional MR imaging. As a result, these tissues appear dark, and their MR properties are inaccessible. However, when ultrashort echo times (UTEs) are used, signals can be detected from these tissues before they decay to zero. This review summarizes recent technical developments in UTE MRI of tissues with short and ultrashort T2 relaxation times. A series of UTE MRI techniques for high-resolution morphological and quantitative imaging of these short-T2 tissues are discussed. Applications of UTE imaging in the musculoskeletal, nervous, respiratory, gastrointestinal, and cardiovascular systems of the body are included.

Ultrashort Echo Time Magnetic Resonance Imaging Techniques: Met and Unmet Needs in Musculoskeletal Imaging

This review article summarizes recent technical developments in ultrashort echo time (UTE) magnetic resonance imaging of musculoskeletal (MSK) tissues with short-T2 relaxation times. A series of contrast mechanisms are discussed for high-contrast morphological imaging of short-T2 MSK tissues including the osteochondral junction, menisci, ligaments, tendons, and bone. Quantitative UTE mapping of T1, T2*, T1ρ, adiabatic T1ρ, magnetization transfer ratio, MT modeling of macromolecular proton fraction, quantitative susceptibility mapping, and water content is also introduced. Met and unmet needs in MSK imaging are discussed. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 3.

Cortical Bone Mechanical Assessment via Free Water Relaxometry at 3 T

BACKGROUND

Investigation of cortical bone using magnetic resonance imaging is a developing field, which uses short/ultrashort echo time (TE) pulse sequences to quantify bone water content and to obtain indirect information about bone microstructure.

PURPOSE

To improve the accuracy of the previously proposed technique of free water T₁ quantification and to seek the relationship between cortical bone free water T₁ and its mechanical competence.

STUDY TYPE

Prospective.

SUBJECTS

Twenty samples of bovine tibia bone.

FIELD STRENGTH/SEQUENCES

3.0 T; ultra-fast two-dimensional gradient echo, Radio frequency-spoiled three-dimensional gradient echo.

ASSESSMENT

Cortical bone free water T₁ was quantified via three different methods: inversion recovery (IR), variable flip angle (VFA), and variable repetition time (VTR). Signal-to-noise ratio was measured by dividing the signal of each segmented sample to background noise. Segmentation was done manually. The effect of noise on T₁ quantification was evaluated. Then, the samples were subjected to mechanical compression test to measure the toughness, yield stress, ultimate stress, and Young modulus.

STATISTICAL TESTS

All the statistical analysis (Shapiro-Wilk, way analysis of variance, paired t test, Pearson correlation, and Bland-Altman plot) were done using SPSS.

RESULTS

Significant difference was found between T₁ quantification groups (P < 0.05). Average T₁ of each quantification method differed significantly after adding noise (P < 0.05). VFA-T₁ values significantly correlated with toughness (r = -0.68, P < 0.05), ultimate stress (r = -0.71, P < 0.05), and yield stress (r = -0.62, P < 0.05). No significant correlation was found between VTR-T₁ values and toughness (P = 0.07), ultimate stress (P = 0.47), yield stress (P = 0.30), and Young modulus (P = 0.39).

DATA CONCLUSION

Pore water T₁ value is associated with bone mechanical competence, and VFA method employing short-TE pulse sequence seems a superior technique to VTR method for this quantification.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: 1.

Association Between T2 * Relaxation Times Derived From Ultrashort Echo Time MRI and Symptoms During Exercise Therapy for Patellar Tendinopathy: A Large Prospective Study

Background

Exercise therapy is considered preferential treatment for patellar tendinopathy (PT). However, there is conflicting evidence for structural patellar tendon adaptation in response to exercise therapy and its association with symptoms is weak.

Purpose

To assess the association between 1) T₂^* relaxation times and symptom severity; 2) baseline T₂^* and clinical outcome; and 3) longitudinal T₂^* changes and clinical outcome in athletes with PT performing exercise therapy.

Study Type

Randomized controlled clinical trial.

Subjects

Seventy‐six athletes (18–35 years) with clinically diagnosed and ultrasound‐confirmed PT.

Field strength/Sequence

3D gradient echo sequence (3.0 T).

Assessment

Patients were enrolled in a randomized trial of progressive tendon‐loading exercises (PTLE) versus eccentric exercise therapy (EET). Symptoms were assessed using the Victorian Institute of Sports Assessment (VISA‐P) questionnaire. 3D‐Ultrashort echo time (UTE)‐MRI was acquired at baseline, 12 and 24 weeks. Voxel‐wise T₂^* relaxation times were quantified using mono‐exponential and bi‐exponential models. T₂^* analysis was performed in three patellar tendon tissue compartments representing: aligned collagen, degenerative tissue, and interface.

Statistical Tests

Adjusted general linear, mixed‐linear models, and generalized estimating equations.

Results

We included 76 patients with PT (58 men, mean age 24 ± 4 years); 38 in the PTLE‐group and 38 in the EET‐group, of which 57 subjects remained eligible for analysis. T₂^* relaxation times were significantly associated with VISA‐P in degenerative and interface tissues of the patellar tendon. No association was found between baseline T₂^* and VISA‐P after 24 weeks (P > 0.29). The estimated mean T₂^* in degenerative tissue decreased from 14 msec (95%CI: 12–16) at baseline to 13 msec (95%CI: 11–15) at 12 weeks and to 13 msec (95%CI: 10–15) at 24 weeks. The significant decrease in T₂^* from baseline to 24 weeks was associated with improved clinical outcome.

Data Conclusion

Tissue‐specific T₂^* relaxation times, identified with 3D‐UTE‐MRI, decreased significantly in athletes with patellar tendinopathy performing exercise therapy and this decrease was associated with improved clinical outcome.

Evidence Level

1

Technical Efficacy

Stage 4

In vivo quantitative imaging biomarkers of bone quality and mineral density using multi-band-SWIFT magnetic resonance imaging

Bone is a composite biomaterial of mineral crystals, organic matrix, and water. Each contributes to bone quality and strength and may change independently, or together, with disease progression and treatment. Even so, there is a near ubiquitous reliance on ionizing x-ray-based approaches to measure bone mineral density (BMD) which is unable to fully characterize bone strength and may not adequately predict fracture risk. Characterization of treatment efficacy in bone diseases of altered remodeling is complicated by the lack of imaging modality able to safely monitor material-level and biochemical changes in vivo. To improve upon the current state of bone imaging, we tested the efficacy of Multi Band SWeep Imaging with Fourier Transformation (MB-SWIFT) magnetic resonance imaging (MRI) as a readout of bone derangement in an estrogen deficient ovariectomized (OVX) rat model during growth. MB-SWIFT MRI-derived BMD correlated significantly with BMD measured using micro-computed tomography (μCT). In this rodent model, growth appeared to overcome estrogen deficiency as bone mass continued to increase longitudinally over the duration of the study. Nonetheless, after 10 weeks of intervention, MB-SWIFT detected significant changes consistent with estrogen deficiency in cortical water, cortical matrix organization (T₁), and marrow fat. Findings point to MB-SWIFT's ability to quantify BMD in good agreement with μCT while providing additive quantitative outcomes about bone quality in a manner consistent with estrogen deficiency. These results indicate MB-SWIFT as a non-ionizing imaging strategy with value for bone imaging and may be a promising technique to progress to the clinic for monitoring and clinical management of patients with bone diseases such as osteoporosis.

Single-sided NMR to estimate morphological parameters of the trabecular bone structure

PURPOSE

Single-sided ¹ H-NMR is proposed for the estimation of morphological parameters of trabecular bone, and potentially the detection of pathophysiological alterations of bone structure. In this study, a new methodology was used to estimate such parameters without using an external reference signal, and to study intratrabecular and intertrabecular porosities, with a view to eventually scanning patients.

METHODS

Animal trabecular bone samples were analyzed by a single-sided device. The Carr-Purcell-Meiboom-Gill sequence of ¹ H nuclei of fluids, including marrow, confined inside the bone, was analyzed by quasi-continuous T₂ distributions and separated into two ¹ H pools: short and long T₂ components. The NMR parameters were estimated using models of trabecular bone structure, and compared with the corresponding micro-CT.

RESULTS

Without any further assumptions, the internal reference parameter (short T₂ signal intensity fraction) enabled prediction of the micro-CT parameters BV/TV (volume of the trabeculae/total sample volume) and BS/TV (external surface of the trabeculae/total sample volume) with linear correlation coefficient >0.80. The assignment of the two pools to intratrabecular and intertrabecular components yielded an estimate of average intratrabecular porosity (33 ± 5)%. Using the proposed models, the NMR-estimated BV/TV and BS/TV were found to be linearly related to the corresponding micro-CT values with high correlation (>0.90 for BV/TV; >0.80 for BS/TV) and agreement coefficients.

CONCLUSION

Low-field, low-cost portable devices that rely on intrinsic magnetic field gradients and do not use ionizing radiation are viable tools for in vitro preclinical studies of pathophysiological structural alterations of trabecular bone.

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.

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.

Near infrared spectroscopic assessment of loosely and tightly bound cortical bone water

Water is an important component of bone and plays a key role in its mechanical and structural integrity. Water molecules in bone are present in different locations, including loosely or tightly bound to the matrix and/or mineral (biological apatite) phases. Identification of water location and interactions with matrix components impact bone function but have been challenging to assess. Here, we used near infrared (NIR) spectroscopy to identify loosely and tightly bound water present in cortical bone. In hydrated samples, NIR spectra have two primary water absorption bands at frequencies of ∼5200 and 7000 cm^-1. Using lyophilization and hydrogen-deuterium exchange assays, we showed that these absorption bands are primarily associated with loosely bound bone water. Using further demineralization assays, thermal denaturation, and comparison to standards, we found that these absorption bands have underlying components associated with water molecules tightly bound to bone. In dehydrated samples, the peak at ∼5200 cm^-1 was assigned to a combination of water tightly bound to collagen and to mineral, whereas the peak at 7000 cm^-1 was exclusively associated with tightly bound mineral water. We also found significant positive correlations between the NIR mineral absorption bands and the mineral content as determined by an established mid infrared spectroscopic parameter, phosphate/amide I. Moreover, the NIR water data showed correlation trends with tissue mineral density (TMD) in cortical bone tissues. These observations reveal the ability of NIR spectroscopy to non-destructively identify loosely and tightly bound water in bone, which could have further applications in biomineralization and biomedical studies.

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.

Tissue-Specific T2 * Biomarkers in Patellar Tendinopathy by Subregional Quantification Using 3D Ultrashort Echo Time MRI

Background

Quantitative MRI of patellar tendinopathy (PT) can be challenging due to spatial variation of T₂* relaxation times.

Purpose

1) To compare T₂* quantification using a standard approach with analysis in specific tissue compartments of the patellar tendon. 2) To evaluate test–retest reliability of different methods for fitting ultrashort echo time (UTE)‐relaxometry data.

Study Type

Prospective.

Subjects

Sixty‐five athletes with PT.

Field Strength/Sequence

3D UTE scans covering the patellar tendon were acquired using a 3.0T scanner and a 16‐channel surface coil.

Assessment

Voxelwise median T₂* was quantified with monoexponential, fractional‐order, and biexponential fitting. We applied two methods for T₂* analysis: first, a standard approach by analyzing all voxels covering the proximal patellar tendon. Second, within subregions of the patellar tendon, by using thresholds on biexponential fitting parameter percentage short T₂* (0–30% for mostly long T₂*, 30–60% for mixed T₂*, and 60–100% for mostly short T₂*).

Statistical Tests

Average test–retest reliability was assessed in three athletes using coefficients‐of‐variation (CV) and coefficients‐of‐repeatability (CR).

Results

With standard image analysis, we found a median [interquartile range, IQR] monoexponential T₂* of 6.43 msec [4.32–8.55] and fractional order T₂* 4.39 msec [3.06–5.78]. The percentage of short T₂* components was 52.9% [35.5–69.6]. Subregional monoexponential T₂* was 13.78 msec [12.11–16.46], 7.65 msec [6.49–8.61], and 3.05 msec [2.52–3.60] and fractional order T₂* 11.82 msec [10.09–14.44], 5.14 msec [4.25–5.96], and 2.19 msec [1.82–2.64] for 0–30%, 30–60%, and 60–100% short T₂*, respectively. Biexponential component short T₂* was 1.693 msec [1.417–2.003] for tissue with mostly short T₂* and long T₂* of 15.79 msec [13.47–18.61] for mostly long T₂*. The average CR (CV) was 2 msec (15%), 2 msec (19%) and 10% (22%) for monoexponential, fractional order and percentage short T₂*, respectively.

Data Conclusion

Patellar tendinopathy is characterized by regional variability in binding states of water. Quantitative multicompartment T₂* analysis in PT can be facilitated using a voxel selection method based on using biexponential fitting parameters.

Level of Evidence

1

Technical Efficacy Stage

1 J. Magn. Reson. Imaging 2020;52:420–430.

Evaluation of Osseous Morphology of the Hip Using Zero Echo Time Magnetic Resonance Imaging

BACKGROUND

Femoroacetabular impingement syndrome (FAIS) is a common disorder of the hip resulting in groin pain and ultimately osteoarthritis. Radiologic assessment of FAI morphologies, which may present with overlapping radiologic features of hip dysplasia, often requires the use of computed tomography (CT) for evaluation of osseous abnormality, owing to the difficulty of direct visualization of cortical and subchondral bone with conventional magnetic resonance imaging (MRI). The use of a zero echo time (ZTE) MRI pulse sequence may obviate the need for CT by rendering bone directly from MRI.

PURPOSE/HYPOTHESIS

The purpose was to explore the application of ZTE MRI to the assessment of osseous FAI and dysplasia morphologies of the hip. It was hypothesized that angular measurements from ZTE images would show significant agreement with measurements obtained from CT images.

STUDY DESIGN

Cohort study (diagnosis); Level of evidence, 2.

METHODS

Thirty-eight hips from 23 patients were imaged with ZTE MRI and CT. Clinically relevant angular measurements of hip morphology were made in both modalities and compared to assess agreement. Measurements included coronal and sagittal center-edge angles, femoral neck-shaft angle, acetabular version (at 1-, 2-, and 3-o'clock positions), Tönnis angle, alpha angle, and modified-beta angle. Interrater agreement was assessed for a subset of 10 hips by 2 raters. Intermodal agreement was assessed on the complete cohort and a single rater.

RESULTS

Interrater agreement was demonstrated in both CT and ZTE, with intraclass correlation coefficient values ranging from 0.636 to 0.990 for ZTE and 0.747 to 0.983 for CT, indicating "good" to "excellent" agreement. Intermodal agreement was also shown to be significant, with intraclass correlation coefficients ranging from 0.618 to 0.904.

CONCLUSION

Significant agreement of angular measurements for hip morphology exists between ZTE MRI and CT imaging. ZTE MRI may be an effective method to quantitatively evaluate osseous hip morphology.

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.

Short-T2 MRI: Principles and recent advances

Among current modalities of biomedical and diagnostic imaging, MRI stands out by virtue of its versatile contrast obtained without ionizing radiation. However, in various cases, e.g., water protons in tissues such as bone, tendon, and lung, MRI performance is limited by the rapid decay of resonance signals associated with short transverse relaxation times T₂ or T₂*. Efforts to address this shortcoming have led to a variety of specialized short-T₂ techniques. Recent progress in this field expands the choice of methods and prompts fresh considerations with regard to instrumentation, data acquisition, and signal processing. In this review, the current status of short-T₂ MRI is surveyed. In an attempt to structure the growing range of techniques, the presentation highlights overarching concepts and basic methodological options. The most frequently used approaches are described in detail, including acquisition strategies, image reconstruction, hardware requirements, means of introducing contrast, sources of artifacts, limitations, and applications.

Environmentally-Controlled Near Infrared Spectroscopic Imaging of Bone Water

We have designed an environmentally-controlled chamber for near infrared spectroscopic imaging (NIRSI) to monitor changes in cortical bone water content, an emerging biomarker related to bone quality assessment. The chamber is required to ensure repeatable spectroscopic measurements of tissues without the influence of atmospheric moisture. A calibration curve to predict gravimetric water content from human cadaveric cortical bone was created using NIRSI data obtained at six different lyophilization time points. Partial least squares (PLS) models successfully predicted bone water content that ranged from 0–10% (R = 0.96, p < 0.05, root mean square error of prediction (RMSEP) = 7.39%), as well as in the physiologic range of 4–10% of wet tissue weight (R = 0.87, p < 0.05, RMSEP = 14.5%). Similar results were obtained with univariate and bivariate regression models for prediction of water in the 0–10% range. Further, we identified two new NIR bone absorbances, at 6560 cm⁻¹ and 6688 cm⁻¹, associated with water and collagen respectively. Such data will be useful in pre-clinical studies that investigate changes in bone quality with disease, aging and with therapeutic use.

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.

Fast quantitative 3D ultrashort echo time MRI of cortical bone using extended cones sampling

PURPOSE

To investigate the effect of stretching sampling window on quantitative 3D ultrashort TE (UTE) imaging of cortical bone at 3 T.

METHODS

Ten bovine cortical bone and 17 human tibial midshaft samples were imaged with a 3T clinical MRI scanner using an 8-channel knee coil. Quantitative 3D UTE imaging biomarkers, including T₁ , T 2 ∗ , magnetization transfer ratio and magnetization transfer modeling, were performed using radial or spiral Cones sampling trajectories with various durations. Errors in UTE-MRI biomarkers as a function of sampling time were evaluated using radial sampling as a reference standard.

RESULTS

For both bovine and human cortical bone samples, no significant differences were observed for all UTE biomarkers (single-component T 2 ∗ , bicomponent T 2 ∗ and relative fractions, T₁ , magnetization transfer ratio, and magnetization transfer modeling of macromolecular fraction) for spiral sampling windows of 992 µs to 1600 µs compared with a radial sampling window of 688 µs.

CONCLUSION

The total scan time can be reduced by 76% with quantification errors less than 5%. Quantitative UTE-MRI techniques can be greatly accelerated using longer sampling windows without significant quantification errors.

Detecting stress injury (fatigue fracture) in fibular cortical bone using quantitative ultrashort echo time-magnetization transfer (UTE-MT): An ex vivo study

Bone stress injury (BSI) incidents have been increasing amongst athletes in recent years as a result of more intense sporting activities. Cortical bone in the tibia and fibula is one of the most common BSI sites. Nowadays, clinical magnetic resonance imaging (MRI) is the recommended technique for BSI diagnosis at an early stage. However, clinical MRI focuses on edema observations in surrounding soft tissues, rather than the injured components of the bone. Specifically, both normal and injured bone are invisible in conventional clinical MRI. In contrast, ultrashort echo time (UTE)-MRI is able to detect the rapidly decaying signal from the bone. This study aimed to employ UTE-MRI for fatigue fracture detection in fibula cortical bone through an ex vivo investigation. Fourteen human fibular samples (47 ± 20 years old, four women) were subjected to cyclic loading on a four-point bending setup. The loading was displacement controlled to induce -5000 ± 1500 μ-strain at 4 Hz. Loading was stopped when bone stiffness was reduced by 20%. Fibula samples were imaged twice, using UTE-MRI and micro-computed tomography (μCT), first pre-loading and second post-loading. After loading, the macromolecular fraction (MMF) from UTE-MT modeling demonstrated a significant decrease (12% ± 20%, P = 0.02) on average. Single-component T2 * also decreased significantly by BSI (12% ± 11%, P = 0.01) on average. MMF reduction is hypothesized to be a result of collagenous matrix rupture and water increase. However, faster T2 * decay might be a result of water shifts towards newly developed microcracks with higher susceptibility. Despite this good sensitivity level of the UTE-MRI technique, the μCT-based porosity at a voxel size of 9 μm was not affected by loading. UTE-MRI shows promise as a new quantitative technique to detect BSI.

Visualization of calcium phosphate cement in teeth by zero echo time 1 H MRI at high field

¹ H magnetic resonance imaging (MRI) by a zero echo time (ZTE) sequence is an excellent method to image teeth. Calcium phosphate cement (CPC) materials are applied in the restoration of tooth lesions, but it has not yet been investigated whether they can be detected by computed tomography (CT) or MRI. The aim of this study was to optimize high-field ZTE imaging to enable the visualization of a new CPC formulation implanted in teeth and to apply this in the assessment of its decomposition in vivo. CPC was implanted in three human and three goat teeth ex vivo and in three goat teeth in vivo. An ultrashort echo time (UTE) sequence with multiple flip angles and echo times was applied at 11.7 T to measure T₁ and T₂ * values of CPC, enamel and dentin. Teeth with CPC were imaged with an optimized ZTE sequence. Goat teeth implanted with CPC in vivo were imaged after 7 weeks ex vivo. T₂ * relaxation of implanted CPC, dentin and enamel was better fitted by a model assuming a Gaussian rather than a Lorentzian distribution. For CPC and human enamel and dentin, the average T₂ * values were 273 ± 19, 562 ± 221 and 476 ± 147 μs, respectively, the average T₂ values were 1234 ± 27, 963 ± 151 and 577 ± 41 μs, respectively, and the average T₁ values were 1065 ± 45, 972 ± 40 and 903 ± 7 ms, respectively. In ZTE images, CPC had a higher signal-to-noise-ratio than dentin and enamel because of the higher water content. Seven weeks after in vivo implantation, the CPC-filled lesions showed less homogeneous structures, a lower T₁ value and T₂ * separated into two components. MRI by ZTE provides excellent contrast for CPC in teeth and allows its decomposition to be followed.

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.

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.

Cortical Bone Porosity: What Is It, Why Is It Important, and How Can We Detect It?

There is growing recognition of the role of micro-architecture in osteoporotic bone loss and fragility. This trend has been driven by advances in imaging technology, which have enabled a transition from measures of mass to micro-architecture. Imaging trabecular bone has been a key research focus, but advances in resolution have also enabled the detection of cortical bone micro-architecture, particularly the network of vascular canals, commonly referred to as 'cortical porosity.' This review aims to provide an overview of what this level of porosity is, why it is important, and how it can be characterized by imaging. Moving beyond a 'trabeculocentric' view of bone loss holds the potential to improve diagnosis and monitoring of interventions. Furthermore, cortical porosity is intimately linked to the remodeling process, which underpins bone loss, and thus a larger potential exists to improve our fundamental understanding of bone health through imaging of both humans and animal models.

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.

Can ultrashort-TE (UTE) MRI sequences on a 3-T clinical scanner detect signal directly from collagen protons: freeze-dry and D2 O exchange studies of cortical bone and Achilles tendon specimens

Ultrashort-TE (UTE) sequences can obtain signal directly from short-T2 , collagen-rich tissues. It is generally accepted that bound and free water can be detected with UTE techniques, but the ability to detect protons directly on the collagen molecule remains controversial. In this study, we investigated the potential of UTE sequences on a 3-T clinical scanner to detect collagen protons via freeze-drying and D2 O-H2 O exchange studies. Experiments were performed on bovine cortical bone and human Achilles tendon specimens, which were either subject to freeze-drying for over 66 h or D2 O-H2 O exchange for 6 days. Specimens were imaged using two- and three-dimensional UTE with Cones trajectory techniques with a minimum TE of 8 μs at 3 T. UTE images before treatment showed high signal from all specimens with bi-component T2 * behavior. Bovine cortical bone showed a shorter T2 * component of 0.36 ms and a longer T2 * component of 2.30 ms with fractions of 78.2% and 21.8% by volume, respectively. Achilles tendon showed a shorter T2 * component of 1.22 ms and a longer T2 * component of 15.1 ms with fractions of 81.1% and 18.9% by volume, respectively. Imaging after freeze-drying or D2 O-H2 O exchange resulted in either the absence or near-absence of signal. These results indicate that bound and free water are the sole sources of UTE signal in bovine cortical bone and human Achilles tendon samples on a clinical 3-T scanner. Protons on the native collagen molecule are not directly visible when imaged using UTE sequences. Copyright © 2016 John Wiley & Sons, Ltd.

Separation of collagen-bound and porous bone water transverse relaxation in mice: proposal of a multi-step approach

The separation and quantification of collagen-bound water (CBW) and pore water (PW) components of the cortical bone signal are important because of their different contribution to bone mechanical properties. Ultrashort TE (UTE) imaging can be used to exploit the transverse relaxation from CBW and PW, allowing their quantification. We tested, for the first time, the feasibility of UTE measurements in mice for the separation and quantification of the transverse relaxation of CBW and PW in vivo using three different approaches for T2 * determination. UTE sequences were acquired at 4.7 T in six mice with 10 different TEs (50-5000 μs). The transverse relaxation time T2 * of CBW (T2 *cbw ) and PW (T2 *pw ) and the CBW fraction (bwf) were computed using a mono-exponential (i), a standard bi-exponential (ii) and a new multi-step bi-exponential (iii) approach. Regions of interest were drawn at multiple levels of the femur and vertebral body cortical bone for each mouse. The sum of the normalized squared residuals (Res) and the homogeneity of variance were tested to compare the different methods. In the femur, approach (i) yielded mean T2 * ± standard deviation (SD) of 657 ± 234 μs. With approach (ii), T2 *cbw , T2 *pw and bwf were 464 ± 153 μs, 15 777 ± 10 864 μs and 57.6 ± 9.9%, respectively. For approach (iii), T2 *cbw , T2 *pw and bwf were 387 ± 108 μs, 7534 ± 2765 μs and 42.5 ± 6.2%, respectively. Similar values were obtained from vertebral bodies. Res with approach (ii) was lower than with the two other approaches (p < 0.007), but T2 *pw and bwf variance was lower with approach (iii) than with approach (ii) (p < 0.048). We demonstrated that the separation and quantification of cortical bone water components with UTE sequences is feasible in vivo in mouse models. The direct bi-exponential approach exhibited the best approximation to the measured signal curve with the lowest residuals; however, the newly proposed multi-step algorithm resulted in substantially lower variability of the computed parameters. Copyright © 2016 John Wiley & Sons, Ltd.

Separation of collagen-bound and porous bone-water longitudinal relaxation in mice using a segmented inversion recovery zero-echo-time sequence

PURPOSE

Cortical bone mechanical properties are related to the collagen-bound water (CBW) and pore water (PW) components of cortical bone. The study evaluates the feasibility of zero-echo-time imaging in mice in vivo for longitudinal relaxation time (T1) measurements in cortical bone and separation of CBW and PW components.

METHODS

Zero-echo-time data were acquired at 4.7 Tesla in six mice with 14 different inversion times (0-2,600 ms). Region-of-interest analysis was performed at level of femur diaphysis. The T1 of cortical bone and of CBW (T1cbw) and PW (T1pw) as well as the CBW fraction (cbwf) was computed using a mono-exponential and a bi-exponential fitting approach, respectively. The sum of the squared residuals (Res) to the fit was provided for both approaches.

RESULTS

For the mono-exponential model, mean T1 ± standard deviation (SD) was 1,057 ± 160 ms. The bi-exponential approach provided a reliable separation of two different bone-water components, with a mean T1cbw of 213 ± 95 ms, T1pw of 2,152 ± 894 ms, and cbwf of 7.4 ± 2.7 %. Lower Res was obtained with bi-exponential approach (P < 0.001), and Res mean values ± SD were 0.016 ± 0.007 (bi-exponential) and 0.033 ± 0.016 (mono-exponential).

CONCLUSION

Zero-echo-time imaging allows for longitudinal relaxation measurements of cortical bone in vivo in mice models, with a reliable separation of PW and CBW components using a bi-exponential curve fitting approach. Magn Reson Med 77:1909-1915, 2017. © 2016 International Society for Magnetic Resonance in Medicine.