A comprehensive set of ultrashort echo time magnetic resonance imaging biomarkers to assess cortical bone health: A feasibility study at clinical field strength

INTRODUCTION

Conventional bone imaging methods primarily use X-ray techniques to assess bone mineral density (BMD), focusing exclusively on the mineral phase. This approach lacks information about the organic phase and bone water content, resulting in an incomplete evaluation of bone health. Recent research highlights the potential of ultrashort echo time magnetic resonance imaging (UTE MRI) to measure cortical porosity and estimate BMD based on signal intensity. UTE MRI also provides insights into bone water distribution and matrix organization, enabling a comprehensive bone assessment with a single imaging technique. Our study aimed to establish quantifiable UTE MRI-based biomarkers at clinical field strength to estimate BMD and microarchitecture while quantifying bound water content and matrix organization.

METHODS

Femoral bones from 11 cadaveric specimens (n = 4 males 67-92 yrs of age, n = 7 females 70-95 yrs of age) underwent dual-echo UTE MRI (3.0 T, 0.45 mm resolution) with different echo times and high resolution peripheral quantitative computed tomography (HR-pQCT) imaging (60.7 μm voxel size). Following registration, a 4.5 mm HR-pQCT region of interest was divided into four quadrants and used across the multi-modal images. Statistical analysis involved Pearson correlation between UTE MRI porosity index and a signal-intensity technique used to estimate BMD with corresponding HR-pQCT measures. UTE MRI was used to calculate T1 relaxation time and a novel bound water index (BWI), compared across subregions using repeated measures ANOVA.

RESULTS

The UTE MRI-derived porosity index and signal-intensity-based estimated BMD correlated with the HR-pQCT variables (porosity: r = 0.73, p = 0.006; BMD: r = 0.79, p = 0.002). However, these correlations varied in strength when we examined each of the four quadrants (subregions, r = 0.11-0.71). T1 relaxometry and the BWI exhibited variations across the four subregions, though these differences were not statistically significant. Notably, we observed a strong negative correlation between T1 relaxation time and the BWI (r = -0.87, p = 0.0006).

CONCLUSION

UTE MRI shows promise for being an innocuous method for estimating cortical porosity and BMD parameters while also giving insight into bone hydration and matrix organization. This method offers the potential to equip clinicians with a more comprehensive array of imaging biomarkers to assess bone health without the need for invasive or ionizing procedures.

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.

A preparation pulse for fast steady state approach in Actual Flip angle Imaging

BACKGROUND

Actual Flip angle Imaging (AFI) is a sequence used for B1 mapping, also embedded in the Variable flip angle with AFI for simultaneous estimation of T1 , B1 and equilibrium magnetization.

PURPOSE

To investigate the design of a preparation module for AFI to allow a fast approach to steady state (SS) without requiring the use of dummy acquisitions.

METHODS

The features of a preparation module with a B1 insensitive adiabatic pulse, spoiler gradients, and a recovery timeT r e c $T_{rec}$ were studied with simulations and validated via experiments and acquired with different k-space traveling strategies. The robustness of the flip angle of the preparation pulse on the acquired signal is studied.

RESULTS

When a 90° adiabatic pulse is used, the forthcomingT r e c $T_{rec}$ can be expressed as a function of repetition times and AFI flip angle only asTR 1 ( n + cos α ) / ( 1 - cos 2 α ) $\mathrm{TR_1}(n+\cos \alpha )/(1-\cos ^2\alpha )$ , where n represents the ratio between the two repetition times of AFI. The robustness of the method is demonstrated by showing that using the values further away from 90° still allows for a faster approach to SS than the use of dummy pulses.

CONCLUSIONS

The preparation module is particularly advantageous for low flip angles, as well as for AFI sequences that sample the center of the k-space early in the sequence, such as centric ordering acquisitions, and for ultrafast EPI-based AFI methods, thus allowing to reduce scanner overhead time.

Design and development of a novel flexible ultra-short echo time (FUSE) sequence

PURPOSE

To present the validation of a new Flexible Ultra-Short Echo time (FUSE) pulse sequence using a short-T2 phantom.

METHODS

FUSE was developed to include a range of RF excitation pulses, trajectories, dimensionalities, and long-T2 suppression techniques, enabling real-time interchangeability of acquisition parameters. Additionally, we developed an improved 3D deblurring algorithm to correct for off-resonance artifacts. Several experiments were conducted to validate the efficacy of FUSE, by comparing different approaches for off-resonance artifact correction, variations in RF pulse and trajectory combinations, and long-T2 suppression techniques. All scans were performed on a 3 T system using an in-house short-T2 phantom. The evaluation of results included qualitative comparisons and quantitative assessments of the SNR and contrast-to-noise ratio.

RESULTS

Using the capabilities of FUSE, we demonstrated that we could combine a shorter readout duration with our improved deblurring algorithm to effectively reduce off-resonance artifacts. Among the different RF and trajectory combinations, the spiral trajectory with the regular half-inc pulse achieves the highest SNRs. The dual-echo subtraction technique delivers better short-T2 contrast and superior suppression of water and agar signals, whereas the off-resonance saturation method successfully suppresses water and lipid signals simultaneously.

CONCLUSION

In this work, we have validated the use of our new FUSE sequence using a short T2 phantom, demonstrating that multiple UTE acquisitions can be achieved within a single sequence. This new sequence may be useful for acquiring improved UTE images and the development of UTE imaging protocols.

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.

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.

CT-like MRI: a qualitative assessment of ZTE sequences for knee osseous abnormalities

OBJECTIVE

To qualitatively evaluate the utility of zero echo-time (ZTE) MRI sequences in identifying osseous findings and to compare ZTE with optimized spoiled gradient echo (SPGR) sequences in detecting knee osseous abnormalities.

MATERIALS AND METHODS

ZTE and standard knee MRI sequences were acquired at 3T in 100 consecutive patients. Three radiologists rated confidence in evaluating osseous abnormalities and image quality on a 5-grade Likert scale in ZTE compared to standard sequences. In a subset of knees (n = 57) SPGR sequences were also obtained, and diagnostic confidence in identifying osseous structures was assessed, comparing ZTE and SPGR sequences. Statistical significance of using ZTE over SPGR was characterized with a paired t-test.

RESULTS

Image quality of the ZTE sequences was rated high by all reviewers with 278 out of 299 (100 studies, 3 radiologists) scores ≥ 4 on the Likert scale. Diagnostic confidence in using ZTE sequences was rated "very high confidence" in 97%, 85%, 71%, and 73% of the cases for osteophytosis, subchondral cysts, fractures, and soft tissue calcifications/ossifications, respectively. In 74% of cases with osseous findings, reviewer scores indicated confidence levels (score ≥ 3) that ZTE sequences improved diagnostic certainty over standard sequences. The diagnostic confidence in using ZTE over SPGR sequences for osseous structures as well as abnormalities was favorable and statistically significant (p < 0.01).

CONCLUSION

Incorporating ZTE sequences in the standard knee MRI protocol was technically feasible and improved diagnostic confidence for osseous findings in relation to standard MR sequences. In comparison to SPGR sequences, ZTE improved assessment of osseous abnormalities.

Zero echo time pediatric musculoskeletal magnetic resonance imaging: initial experience

BACKGROUND

Projection radiography (XR) is often supplemented by both CT (to evaluate osseous structures with ionizing radiation) and MRI (for marrow and soft-tissue assessment). Zero echo time (ZTE) MR imaging produces a "CT-like" osseous contrast that might obviate CT.

OBJECTIVE

This study investigated our institution's initial experience in implementing an isotropic ZTE MR imaging sequence for pediatric musculoskeletal examinations.

MATERIALS AND METHODS

Pediatric patients referred for extremity MRI at 3 tesla (T) underwent ZTE MR imaging to yield images with contrast similar to that of CT. A radiograph-like image was also created with ray-sum image processing. We assessed ZTE-CT/XR anatomical image quality (S_anat) from 0 (nondiagnostic) to 5 (outstanding). Further, we made image comparisons on a 5-point scale (S_comp) (range of -2 = conventional CT/XR greater anatomical delineation to +2 = ZTE-CT/XR greater anatomical delineation; 0=same) for three cohorts: (1) ZTE-XR to conventional radiography, (2) ZTE-CT to conventional CT and (3) pathological lesion assessment on ZTE-XR to conventional radiography. We measured cortical thickness of ZTE-XR and ZTE-CT and compared these with conventional imaging. We calculated confidence interval of proportions, Wilcoxon rank sum test and intraclass correlation coefficients for inter-reader agreement.

RESULTS

Cohorts 1, 2 and 3 consisted of 40, 20 and 35 cases, respectively (age range 0.6-23.0 years). ZTE-CT versus CT and ZTE-XR versus radiography of cortical thicknesses were not significantly different (P=0.55 and P=0.31, respectively). Cortical delineation was rated diagnostic or better (score of 3, 4 or 5) in all cases (confidence interval of proportions = 100%) for ZTE-CT/XR. Similarly, intramedullary cavity delineation was rated diagnostic or better in all cases for ZTE-CT, and ZTE-XR was at least diagnostic in 58-63% of cases. For cohort 2, cortex and intramedullary cavity S_comp for ZTE-CT was comparable to those of conventional CT, with confidence interval of proportion (sum of score of -1 to +2) of 93-100% and 95%, respectively. Pathology visualized on ZTE-CT/XR was comparable; S_comp confidence interval of proportions was 95%/97-100%, with improved delineation of non-displaced fractures on ZTE-XR. Readers had moderate to near-perfect intraclass correlation coefficient (range=0.60-0.93).

CONCLUSION

Implementation of a diagnostic-quality ZTE MRI sequence in the pediatric population is feasible and can be performed as a complementary pulse sequence to enhance musculoskeletal MRI studies. Compared to conventional CT, ZTE has comparable cortical delineation, intramedullary cavity and pathology visualization. While not intended as a replacement for conventional radiography, ZTE-XR provides similar visualization of pathology.

Bone Phenotyping Approaches in Human, Mice and Zebrafish - Expert Overview of the EU Cost Action GEMSTONE ("GEnomics of MusculoSkeletal traits TranslatiOnal NEtwork")

A synoptic overview of scientific methods applied in bone and associated research fields across species has yet to be published. Experts from the EU Cost Action GEMSTONE ("GEnomics of MusculoSkeletal Traits translational Network") Working Group 2 present an overview of the routine techniques as well as clinical and research approaches employed to characterize bone phenotypes in humans and selected animal models (mice and zebrafish) of health and disease. The goal is consolidation of knowledge and a map for future research. This expert paper provides a comprehensive overview of state-of-the-art technologies to investigate bone properties in humans and animals - including their strengths and weaknesses. New research methodologies are outlined and future strategies are discussed to combine phenotypic with rapidly developing -omics data in order to advance musculoskeletal research and move towards "personalised medicine".

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.

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.

Characterization of hardware-related spatial distortions for IR-PETRA pulse sequence using a brain specific phantom

OBJECTIVE

Inversion recovery-pointwise encoding time reduction with radial acquisition (IR-PETRA) is an effective magnetic resonance (MR) pulse sequence in generating pseudo-CTs. The hardware-related spatial-distortion (HRSD) in MR images potentially deteriorates the accuracy of pseudo-CTs. Thus, we aimed at characterizing HRSD for IR-PETRA.

MATERIALS AND METHODS

gross-HRSD_overall (Euclidean-sum of gross-HRSD_i (i = x, y, z)) for IR-PETRA was assessed using a brain-specific phantom for two MR scanners (1.5 T-Aera and 3.0 T-Prisma). Moreover, hardware imperfections were analyzed by determining gradient-nonlinearity spatial-distortion (GNSD) and B₀-inhomogeneity spatial-distortion (B₀ISD) for magnetization-prepared rapid acquisition gradient-echo (MP-RAGE) which has well-known distortion characteristics.

RESULTS

In 3.0 T, maximum of gross-GNSD_overall (Euclidean-sum of gross-GNSD_i) and gross-B₀ISD for MP-RAGE was 2.77 mm and 0.57 mm, respectively. For this scanner, the mean and maximum of gross-HRSD_overall for IR-PETRA were 0.63 ± 0.38 mm and 1.91 mm, respectively. In 1.5 T, maximum of gross-GNSD_overall and gross-B₀ISD for MP-RAGE was 3.41 mm and 0.78 mm, respectively. The mean and maximum of gross-HRSD_overall for IR-PETRA were 1.02 ± 0.50 mm and 3.12 mm, respectively.

DISCUSSION

The spatial accuracy of MR images, besides being impacted by hardware performance, scanner capabilities, and imaging parameters, is mainly affected by its imaging strategy and data acquisition scheme. In 3.0 T, even without applying vendor correction algorithms, spatial accuracy of IR-PETRA image is sufficient for generating pseudo-CTs. In 1.5 T, distortion-correction is required to provide this accuracy.

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.

T1 measurement of bound water in cortical bone using 3D adiabatic inversion recovery ultrashort echo time (3D IR-UTE) Cones imaging

PURPOSE

We describe the measurement of bound water T₁ ( T 1 BW ) of cortical bone in vitro and in vivo with a 3D adiabatic inversion recovery ultrashort echo time (IR-UTE) Cones sequence using a clinical 3T scanner.

METHODS

A series IR-UTE data from 6 repetition times (TRs) with 5 inversion times (TIs) at each TR were acquired from 12 human tibial bone specimens, and data from 4 TRs with 5 TIs at each TR were acquired from the tibial midshafts of 8 healthy volunteers. The pore water nulling point was calculated from exponential fitting of the inversion recovery curve at each TR. Bone specimens and volunteers were then scanned again with the calculated nulling point at each TR. T 1 BW was derived through exponential fitting of data from IR-UTE images acquired at different TRs using the calculated pore water nulling point for each TR.

RESULTS

In vitro pore water nulling TIs were 141.3 ± 11.6, 123.4 ± 8.9, 101.3 ± 6.2, 88.9 ± 5.3, 74.8 ± 4.2, and 59.2 ± 3.9 ms for the 6 TRs of 500, 400, 300, 250, 200, and 150 ms, respectively. In vivo pore water nulling TIs were 132.8 ± 12.8, 110.3 ± 10.0, 80.0 ± 7.2, and 63.9 ± 5.4 ms for the 4 TRs of 400, 300, 200, and 150 ms, respectively. Excellent exponential fitting was achieved for IR-UTE imaging of bound water with pore water nulled at each TR. The mean T 1 BW was 106.9 ± 6.3 ms in vitro and 112.3 ± 16.4 ms in vivo.

CONCLUSION

Using the 3D IR-UTE Cones with a variable TR/TI approach, T 1 BW of cortical bone was calculated after complete nulling of pore water signals.

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.

Significant correlations between human cortical bone mineral density and quantitative susceptibility mapping (QSM) obtained with 3D Cones ultrashort echo time magnetic resonance imaging (UTE-MRI)

PURPOSE

Quantitative susceptibility mapping (QSM) MRI is a tool that can characterize changes in susceptibility, an intrinsic property which is associated with compositional changes in the tissue. Current QSM estimation of cortical bone is challenging because conventional clinical MRI cannot acquire signal in cortical bone. This study aimed to implement Cones 3D ultrashort echo time MRI (UTE-MRI) for ex vivo QSM measurements in human tibial cortical bone, investigating the correlations of QSM with volumetric intracortical bone mineral density (BMD).

MATERIALS AND METHODS

Nine tibial midshaft cortical bone specimens (25 mm long specimens cut at the mid-point of tibial shaft, 67 ± 20 years old, 5 women and 4 men) were scanned on a clinical 3 T MRI scanner for QSM measurement. The specimens were also scanned on a high-resolution micro-computed tomography (μCT) scanner for volumetric BMD estimation. QSM and μCT results were compared at approximately nine regions of interest (ROIs) per specimen.

RESULTS

Average 3D UTE-MRI QSM showed significantly strong correlation with volumetric BMD (R = -0.82, P < 0.01) and bone porosity (R = 0.72, P < 0.01). Combining all data points together (77 ROIs), QSM showed significant moderate to strong correlation with volumetric BMD after correction for interdependencies in specimens (R = -0.70, P < 0.01). The corrections were required because the data points were not independent in each specimen. Similarly, the correlation between QSM and porosity was significant (R = 0.68, P < 0.01).

CONCLUSIONS

These results suggest that the Cones 3D UTE-MRI QSM technique can potentially serve as a novel and accurate tool to assess intracortical bone mineral density whilst avoiding ionizing radiation.

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.

Ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) for simultaneous quantification of long and ultrashort T2 tissues

PURPOSE

To demonstrate an ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) method that allows quantifying relaxation times for muscle and bone in the musculoskeletal system and generating bone enhanced images that mimic CT scans.

METHODS

A fast imaging steady-state free precession MRF sequence with half pulse excitation and half projection readout was designed to sample fast T₂ decay signals. Varying echo time (TE) of a sinusoidal pattern was applied to enhance sensitivity for tissues with short and ultrashort T₂ values. The performance of UTE-MRF was evaluated via simulations, phantom, and in vivo experiments.

RESULTS

A minimal TE of 0.05 ms was achieved. Simulations indicated the sinusoidal TE sampling increased T₂ quantification accuracy in the cortical bone and tendon but had little impact on long T₂ muscle quantifications. For the rubber phantom, the averaged relaxometries from UTE-MRF (T₁ = 162 ms and T₂ = 1.07 ms) compared well with the gold standard (T₁ = 190 ms and T 2 ∗ = 1.03 ms). For the long T₂ agarose phantom, the linear regression slope between UTE-MRF and gold standard was 1.07 (R² = 0.991) for T₁ and 1.04 (R² = 0.994) for T₂ . In vivo experiments showed the detection of the cortical bone (averaged T₂ = 1.0 ms) and Achilles tendon (averaged T₂ = 15 ms). Scalp structures from the bone enhanced image show high similarity with CT.

CONCLUSION

The UTE-MRF with sinusoidal TEs can simultaneously quantify T₁ , T₂ , proton density, and B₀ in long, short, even ultrashort T₂ musculoskeletal structures. Bone enhanced images can be achieved in the brain with UTE-MRF.

Fat suppression for ultrashort echo time imaging using a single-point Dixon method

PURPOSE

In ultrashort echo time (UTE) imaging, fat suppression can improve short T₂ * contrast but can also reduce short T₂ * signals. The conventional two-point Dixon (2p-Dixon) method does not perform well due to short T₂ * decay. In this study, we propose a new method to suppress fat for high contrast UTE imaging of short T₂ tissues, utilizing a single-point Dixon (1p-Dixon) method.

METHODS

The proposed method utilizes dual-echo UTE imaging, where UTE is followed by the second TE, chosen flexibly. Fat is estimated by applying a 1p-Dixon method to the non-UTE image after correction of phase errors, which is used to suppress fat in the UTE image. In vivo ankle and knee imaging were performed at 3 T to evaluate the proposed method.

RESULT

It was observed that fat and water signals in tendons were misestimated by the 2p-Dixon method due to signal decay, while the 1p-Dixon method showed reliable fat and water separation not affected by the short T₂ * signal decay. Compared with the conventional chemical shift based fat saturation technique, the 1p-Dixon based approach showed much stronger signal intensities in the Achilles, quadriceps, and patellar tendons, with significantly improved contrast to noise ratios (CNRs) of 11.8 ± 2.2, 16.0 ± 1.6, and 26.8 ± 1.3 with the 1p-Dixon method and 0.6 ± 0.2, 4.6 ± 1.0, and 17.5 ± 1.4 with regular fat saturation, respectively.

CONCLUSION

The proposed 1p-Dixon based fat suppression allows more flexible selection of imaging parameters and more accurate fat and water separation over the conventional 2p-Dixon in UTE imaging. Moreover, the proposed method provides much improved CNR for short T₂ tissues over the conventional fat saturation method.

Recent Advances in the Genetics of Fractures in Osteoporosis

Genetic susceptibility, together with old age, female sex, and low bone mineral density (BMD) are amongst the strongest determinants of fracture risk. Tmost recent large-scale genome-wide association study (GWAS) meta-analysis has yielded fifteen loci. This review focuses on the advances in the research of genetic determinants of fracture risk. We first discuss the genetic architecture of fracture risk, touching upon different methods and overall findings. We then discuss in a second paragraph the most recent advances in the field and focus on the genetics of fracture risk and also of other endophenotypes closely related to fracture risk such as bone mineral density (BMD). Application of state-of-the-art methodology such as Mendelian randzation in fracture GWAS are reviewed. The final part of this review touches upon potential future directions in genetic research of osteoporotic fractures.

Accurate T1 mapping of short T2 tissues using a three-dimensional ultrashort echo time cones actual flip angle imaging-variable repetition time (3D UTE-Cones AFI-VTR) method

PURPOSE

To develop an accurate T1 measurement method for short T2 tissues using a combination of a 3-dimensional ultrashort echo time cones actual flip angle imaging technique and a variable repetition time technique (3D UTE-Cones AFI-VTR) on a clinical 3T scanner.

METHODS

First, the longitudinal magnetization mapping function of the excitation pulse was obtained with the 3D UTE-Cones AFI method, which provided information about excitation efficiency and B1 inhomogeneity. Then, the derived mapping function was substituted into the VTR fitting to generate accurate T1 maps. Numerical simulation and phantom studies were carried out to compare the AFI-VTR method with a B1 -uncorrected VTR method, a B1 -uncorrected variable flip angle (VFA) method, and a B1 -corrected VFA method. Finally, the 3D UTE-Cones AFI-VTR method was applied to bovine bone samples (N = 6) and healthy volunteers (N = 3) to quantify the T1 of cortical bone.

RESULTS

Numerical simulation and phantom studies showed that the 3D UTE-Cones AFI-VTR technique provides more accurate measurement of the T1 of short T2 tissues than the B1 -uncorrected VTR and VFA methods or the B1 -corrected VFA method. The proposed 3D UTE-Cones AFI-VTR method showed a mean T1 of 240 ± 25 ms for bovine cortical bone and 218 ± 10 ms for the tibial midshaft of human volunteers, respectively, at 3 T.

CONCLUSION

The 3D UTE-Cones AFI-VTR method can provide accurate T1 measurements of short T2 tissues such as cortical bone. Magn Reson Med 80:598-608, 2018. © 2018 International Society for Magnetic Resonance in Medicine.