Quantitation of the susceptibility difference between trabecular bone and bone marrow: experimental studies

More accurate trabecular bone imaging using UTE MRI at the resonance frequency of fat

High-resolution magnetic resonance imaging (HR-MRI) has been increasingly used to assess the trabecular bone structure. High susceptibility at the marrow/bone interface may significantly reduce the marrow's apparent transverse relaxation time (T2*), overestimating trabecular bone thickness. Ultrashort echo time MRI (UTE-MRI) can minimize the signal loss caused by susceptibility-induced T2* shortening. However, UTE-MRI is sensitive to chemical shift artifacts, which manifest as spatial blurring and ringing artifacts partially due to non-Cartesian sampling. In this study, we proposed UTE-MRI at the resonance frequency of fat to minimize marrow-related chemical shift artifacts and the overestimation of trabecular thickness. Cubes of trabecular bone from six donors (75 ± 4 years old) were scanned using a 3 T clinical scanner at the resonance frequencies of fat and water, respectively, using 3D UTE sequences with five TEs (0.032, 1.1, 2.2, 3.3, and 4.4 ms) and a clinical 3D gradient echo (GRE) sequence at 0.2 × 0.2 × 0.4 mm3 voxel size. Trabecular bone thickness was measured in 30 regions of interest (ROIs) per sample. MRI results were compared with thicknesses obtained from micro-computed tomography (μCT) at 50 μm3 voxel size. Linear regression models were used to calculate the coefficient of determination between MRI- and μCT-based trabecular thickness. All MRI-based trabecular thicknesses showed significant correlations with μCT measurements. The correlations were higher (examined with paired Student's t-test, P < 0.01) for 3D UTE images performed at the fat frequency (R2 = 0.59-0.74, P < 0.01) than those at the water frequency (R2 = 0.18-0.52, P < 0.01) and clinical GRE images (R2 = 0.39-0.47, P < 0.01). Significantly reduced correlations were observed with longer TEs. This study highlighted the feasibility of UTE-MRI at the fat frequency for a more accurate assessment of trabecular bone thickness.

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.

Compositional evaluation of lesion and parent bone in patients with juvenile osteochondritis dissecans of the knee using T2 * mapping

Juvenile osteochondritis dissecans (JOCD) lesions contain cartilaginous, fibrous and osseous tissues which are difficult to distinguish with clinical, morphological magnetic resonance imaging (MRI). Quantitative T₂ * mapping has earlier been used to evaluate microstructure and composition of all aforementioned tissues as well as bone mineral density. However, the ability of T₂ * mapping to detect changes in tissue composition between different JOCD lesion regions, different disease stages, and between stable and unstable lesions has not been demonstrated. This study analyzed morphological and T₂ * MRI data from 25 patients (median age, 12.1 years) with 34 JOCD-affected and 13 healthy knees. Each lesion was assigned a stage reflecting the natural history of JOCD, with stages I and IV representing early and healed lesion, respectively. T₂ * values were evaluated within the progeny lesion, interface and parent bone of each lesion and in the control bone region. T₂ * was negatively correlated with JOCD stage in progeny lesion (ρ = -0.871; p < 0.001) and interface regions (ρ = -0.649; p < 0.001). Stage IV progeny showed significantly lower T₂ * than control bone (p = 0.028). T₂ * was significantly lower in parent bone than in control bone of patients with stable lesions (p = 0.009), but not in patients with unstable lesions (p = 0.14). Clinical significance: T₂ * mapping enables differentiation between different stages of JOCD and quantitative measurement of the ossification degree in progeny lesion and interface. The observed T₂ * decrease in healed and stable lesions may indicate increased bone density as a result of the active repair process. T₂ * mapping provides quantitative information about JOCD lesion composition.

Imaging of the Osteoporotic Spine - Quantitative Approaches in Diagnostics and for the Prediction of the Individual Fracture Risk

Osteoporosis is a highly prevalent systemic skeletal disease that is characterized by low bone mass and microarchitectural bone deterioration. It predisposes to fragility fractures that can occur at various sites of the skeleton, but vertebral fractures (VFs) have been shown to be particularly common. Prevention strategies and timely intervention depend on reliable diagnosis and prediction of the individual fracture risk, and dual-energy X-ray absorptiometry (DXA) has been the reference standard for decades. Yet, DXA has its inherent limitations, and other techniques have shown potential as viable add-on or even stand-alone options. Specifically, three-dimensional (3 D) imaging modalities, such as computed tomography (CT) and magnetic resonance imaging (MRI), are playing an increasing role. For CT, recent advances in medical image analysis now allow automatic vertebral segmentation and value extraction from single vertebral bodies using a deep-learning-based architecture that can be implemented in clinical practice. Regarding MRI, a variety of methods have been developed over recent years, including magnetic resonance spectroscopy (MRS) and chemical shift encoding-based water-fat MRI (CSE-MRI) that enable the extraction of a vertebral body's proton density fat fraction (PDFF) as a promising surrogate biomarker of bone health. Yet, imaging data from CT or MRI may be more efficiently used when combined with advanced analysis techniques such as texture analysis (TA; to provide spatially resolved assessments of vertebral body composition) or finite element analysis (FEA; to provide estimates of bone strength) to further improve fracture prediction. However, distinct and experimentally validated diagnostic criteria for osteoporosis based on CT- and MRI-derived measures have not yet been achieved, limiting broad transfer to clinical practice for these novel approaches. KEY POINTS::   · DXA is the reference standard for diagnosis and fracture prediction in osteoporosis, but it has important limitations.. · CT- and MRI-based methods are increasingly used as (opportunistic) approaches.. · For CT, particularly deep-learning-based automatic vertebral segmentation and value extraction seem promising.. · For MRI, multiple techniques including spectroscopy and chemical shift imaging are available to extract fat fractions.. · Texture and finite element analyses can provide additional measures for vertebral body composition and bone strength.. CITATION FORMAT: · Sollmann N, Kirschke JS, Kronthaler S et al. Imaging of the Osteoporotic Spine - Quantitative Approaches in Diagnostics and for the Prediction of the Individual Fracture Risk. Fortschr Röntgenstr 2022; DOI: 10.1055/a-1770-4626.

On quantification errors of R 2 * $$ {R}_2^{\ast } $$ and proton density fat fraction mapping in trabecularized bone marrow in the static dephasing regime

PURPOSE

To study the effect of field inhomogeneity distributions in trabecularized bone regions on the gradient echo (GRE) signal with short TEs and to characterize quantification errors on R 2 * $$ {R}_2^{\ast } $$ and proton density fat fraction (PDFF) maps when using a water-fat model with an exponential R 2 * $$ {R}_2^{\ast } $$ decay model at short TEs.

METHODS

Field distortions were simulated based on a trabecular bone micro CT dataset. Simulations were performed for different bone volume fractions (BV/TV) and for different bone-fat composition values. A multi-TE UTE acquisition was developed to acquire multiple UTEs with random order to minimize eddy currents. The acquisition was validated in phantoms and applied in vivo in a volunteer's ankle and knee. Chemical shift encoded MRI (CSE-MRI) based on a Cartesian multi-TE GRE scan was acquired in the spine of patients with metastatic bone disease.

RESULTS

Simulations showed that signal deviations from the exponential signal decay at short TEs were more prominent for a higher BV/TV. UTE multi-TE measurements reproduced in vivo the simulation-based predicted behavior. In regions with high BV/TV, the presence of field inhomogeneities induced an R 2 * $$ {R}_2^{\ast } $$ underestimation in trabecularized bone marrow when using CSE-MRI at 3T with a short TE.

CONCLUSION

R 2 * $$ {R}_2^{\ast } $$ can be underestimated when using short TEs (<2 ms at 3 T) and a water-fat model with an exponential R 2 * $$ {R}_2^{\ast } $$ decay model in multi-echo GRE acquisitions of trabecularized bone marrow.

Quantitative assessment of lumbar spine bone marrow in patients with different severity of CKD by IDEAL-IQ magnetic resonance sequence

BACKGROUND

Chronic kidney disease (CKD) has a significant negative impact on bone health. Bone marrow is an essential component of bone, mainly composed of trabecular bone and fat. The IDEAL-IQ sequence of MRI allows indirect quantification of trabecular bone mass by R2* and direct quantification of bone marrow fat content by FF map, respectively.

OBJECTIVE

Our objective was to explore the association of CKD severity with bone marrow using IDEAL-IQ and whether mineral and bone metabolism markers alter this association.

METHOD

We recruited 68 CKD patients in this cross-sectional research (15 with CKD stages 3-4, 26 with stage 5, and 27 with stage 5d). All patients underwent lumbar spine IDEAL-IQ, BMD, and several bone metabolism markers (iPTH, 25-(OH)-VitD, calcium and phosphorus). Multiple linear regression analysis was used to examine the association of CKD severity with MRI measurements (R2* and FF).

RESULTS

More severe CKD was associated with a higher R2* value [CKD 5d versus 3-4: 30.077 s^-1 (95% CI: 12.937, 47.217), P for trend < 0.001], and this association was attenuated when iPTH was introduced [CKD 5d versus 3-4: 19.660 s^-1 (95% CI: 0.205, 39.114), P for trend = 0.042]. Furthermore, iPTH had an association with R2* value [iPTH (pg/mL): 0.033 s^-1 (95% CI: 0.001, 0.064), P = 0.041]. Besides, FF was mainly affected by age and BMI, but not CKD.

CONCLUSIONS

The bone marrow R2* value measured by IDEAL-IQ sequence is associated with CKD severity and iPTH. The R2* of IDEAL-IQ has the potential to reflect lumbar bone changes in patients with CKD.

Voxel-based mapping of five MR biomarkers in the wrist bone marrow

OBJECTIVE

MRI is a reliable and accurate technique to characterize rheumatoid arthritis. The aim of this study was to provide voxel-by-voxel 3D maps of the proton density fat fraction (PDFF), the T₁ of water (T_1W), the T₁ of fat (T_1F), the T₂* of water (T₂*_W), the T₂* of fat (T₂*_F) in the wrist bone marrow.

MATERIALS AND METHODS

The experiments were conducted on 14 healthy volunteers (mean age: 24 ± 4). The data were acquired at 1.5 T using two optimized four-echo 3D 1.2 × 1.2 × 1.2 mm³-isotropic spoiled gradient sequences. A repeatability study was carried out. The measurements were done using a homemade parametric viewer software.

RESULTS

The inter-volunteer results were, on average: PDFF = 86 ± 3%, T_1W = 441 ± 113 ms, T_1F = 245 ± 19 ms, T₂*_W = 6 ± 1 ms and T₂*_F = 16 ± 3 ms. The coefficients of variation were for fat based biomarkers CV_PDFF < 5%, CV_T1F < 15% and CV_T2*F < 10% in the repeatability study.

DISCUSSION

The protocol and quantification tool proposed in this study provide high-resolution voxel-by-voxel 3D maps of five biomarkers in the wrist in less than 4 min of acquisition.

Vertebral bone marrow T2* mapping using chemical shift encoding-based water-fat separation in the quantitative analysis of lumbar osteoporosis and osteoporotic fractures

BACKGROUND

Chemical shift encoding-based water-fat separation techniques have been used for fat quantification [proton density fat fraction (PDFF)], but they also enable the assessment of bone marrow T2*, which has previously been reported to be a potential biomarker for osteoporosis and may give insight into the cause of vertebral fractures (i.e., osteoporotic vs. traumatic) and the microstructure of the bone when applied to vertebral bone marrow.

METHODS

The 32 patients (78.1% with low-energy osteopenic/osteoporotic fractures, mean age 72.3±9.8 years, 76% women; 21.9% with high-energy traumatic fractures, 47.3±12.8 years, no women) were frequency-matched for age and sex to subjects without vertebral fractures (n=20). All study patients underwent 3T-MRI of the lumbar spine including sagittally acquired spoiled gradient echo sequences for chemical shift encoding-based water-fat separation, from which T2* values were obtained. Volumetric trabecular bone mineral density (BMD) and trabecular bone parameters describing the three-dimensional structural integrity of trabecular bone were derived from quantitative CT. Associations between T2* measurements, fracture status and trabecular bone parameters were assessed using multivariable linear regression models.

RESULTS

Mean T2* values of non fractured vertebrae in all patients showed a significant correlation with BMD (r=-0.65, P<0.001), trabecular number (TbN) (r=-0.56, P<0.001) and trabecular spacing (TbSp) (r=0.61, P<0.001); patients with low-energy osteoporotic vertebral fractures showed significantly higher mean T2* values than those with traumatic fractures (13.6±4.3 vs. 8.4±2.2 ms, P=0.01) as well as a significantly lower TbN (0.69±0.08 vs. 0.93±0.03 mm-1, P<0.01) and a significantly larger trabecular spacing (1.06±0.16 vs. 0.56±0.08 mm, P<0.01). Mean T2* values of osteoporotic patients with and without vertebral fracture showed no significant difference (13.5±3.4 vs. 15.6±3.5 ms, P=0.40). When comparing the mean T2* of the fractured vertebrae, no significant difference could be detected between low-energy osteoporotic fractures and high-energy traumatic fractures (12.6±5.4 vs. 8.1±2.4 ms, P=0.10).

CONCLUSIONS

T2* mapping of vertebral bone marrow using using chemical shift encoding-based water-fat separation allows for assessing osteoporosis as well as the trabecular microstructure and enables a radiation-free differentiation between patients with low-energy osteoporotic and high-energy traumatic vertebral fractures, suggesting its potential as a biomarker for bone fragility.

Texture Analysis Using CT and Chemical Shift Encoding-Based Water-Fat MRI Can Improve Differentiation Between Patients With and Without Osteoporotic Vertebral Fractures

PURPOSE

Osteoporosis is a highly prevalent skeletal disease that frequently entails vertebral fractures. Areal bone mineral density (BMD) derived from dual-energy X-ray absorptiometry (DXA) is the reference standard, but has well-known limitations. Texture analysis can provide surrogate markers of tissue microstructure based on computed tomography (CT) or magnetic resonance imaging (MRI) data of the spine, thus potentially improving fracture risk estimation beyond areal BMD. However, it is largely unknown whether MRI-derived texture analysis can predict volumetric BMD (vBMD), or whether a model incorporating texture analysis based on CT and MRI may be capable of differentiating between patients with and without osteoporotic vertebral fractures.

MATERIALS AND METHODS

Twenty-six patients (15 females, median age: 73 years, 11 patients showing at least one osteoporotic vertebral fracture) who had CT and 3-Tesla chemical shift encoding-based water-fat MRI (CSE-MRI) available were analyzed. In total, 171 vertebral bodies of the thoracolumbar spine were segmented using an automatic convolutional neural network (CNN)-based framework, followed by extraction of integral and trabecular vBMD using CT data. For CSE-MRI, manual segmentation of vertebral bodies and consecutive extraction of the mean proton density fat fraction (PDFF) and T2* was performed. First-order, second-order, and higher-order texture features were derived from texture analysis using CT and CSE-MRI data. Stepwise multivariate linear regression models were computed using integral vBMD and fracture status as dependent variables.

RESULTS

Patients with osteoporotic vertebral fractures showed significantly lower integral and trabecular vBMD when compared to patients without fractures (p<0.001). For the model with integral vBMD as the dependent variable, T2* combined with three PDFF-based texture features explained 40% of the variance (adjusted R²[Ra2] = 0.40; p<0.001). Furthermore, regarding the differentiation between patients with and without osteoporotic vertebral fractures, a model including texture features from CT and CSE-MRI data showed better performance than a model based on integral vBMD and PDFF only ( Ra2 = 0.47 vs. Ra2  = 0.81; included texture features in the final model: integral vBMD, CT_Short-run_emphasis, CT_Varianceglobal, and PDFF_Variance).

CONCLUSION

Using texture analysis for spine CT and CSE-MRI can facilitate the differentiation between patients with and without osteoporotic vertebral fractures, implicating that future fracture prediction in osteoporosis may be improved.

Musculoskeletal MR Imaging Applications at Ultra-High (7T) Field Strength

Regulatory approval of ultrahigh field (UHF) MR imaging scanners for clinical use has opened new opportunities for musculoskeletal imaging applications. UHF MR imaging has unique advantages in terms of signal-to-noise ratio, contrast-to-noise ratio, spectral resolution, and multinuclear applications, thus providing unique information not available at lower field strengths. But UHF also comes with a set of technical challenges that are yet to be resolved and may not be suitable for all imaging applications. This review focuses on the latest research in musculoskeletal MR imaging applications at UHF including morphologic imaging, T₂, T₂∗, and T_1ρ mapping, chemical exchange saturation transfer, sodium imaging, and phosphorus spectroscopy imaging applications.

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.

MRI-Based Quantitative Osteoporosis Imaging at the Spine and Femur

Osteoporosis is a systemic skeletal disease with a high prevalence worldwide, characterized by low bone mass and microarchitectural deterioration, predisposing an individual to fragility fractures. Dual-energy X-ray absorptiometry (DXA) has been the clinical reference standard for diagnosing osteoporosis and for assessing fracture risk for decades. However, other imaging modalities are of increasing importance to investigate the etiology, treatment, and fracture risk. The purpose of this work is to review the available literature on quantitative magnetic resonance imaging (MRI) methods and related findings in osteoporosis at the spine and proximal femur as the clinically most important fracture sites. Trabecular bone microstructure analysis at the proximal femur based on high-resolution MRI allows for a better prediction of osteoporotic fracture risk than DXA-based bone mineral density (BMD) alone. In the 1990s, T₂ * mapping was shown to correlate with the density and orientation of the trabecular bone. Recently, quantitative susceptibility mapping (QSM), which overcomes some of the limitations of T₂ * mapping, has been applied for trabecular bone quantifications at the spine, whereas ultrashort echo time (UTE) imaging provides valuable surrogate markers of cortical bone quantity and quality. Magnetic resonance spectroscopy (MRS) and chemical shift encoding-based water-fat MRI (CSE-MRI) enable the quantitative assessment of the nonmineralized bone compartment through extraction of the bone marrow fat fraction (BMFF). Furthermore, CSE-MRI allows for the differentiation of osteoporotic vs. pathologic fractures, which is of high clinical relevance. Lastly, advanced postprocessing and image analysis tools, particularly considering statistical parametric mapping and region-specific BMFF distributions, have high potential to further improve MRI-based fracture risk assessments at the spine and hip. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 2.

Trabecular bone imaging using a 3D adiabatic inversion recovery prepared ultrashort TE Cones sequence at 3T

PURPOSE

To investigate direct imaging of trabecular bone using a 3D adiabatic inversion recovery prepared ultrashort TE cones (3D IR-UTE-Cones) sequence.

METHODS

The proposed 3D IR-UTE-Cones sequence used a broadband adiabatic inversion pulse together with a short TR/TI combination to suppress signals from long T2 tissues such as muscle and marrow fat, followed by multispoke UTE acquisition to detect signal from short T2 water components in trabecular bone. The feasibility of this technique for robust suppression of long T2 tissues was first demonstrated through numerical simulations. The proposed IR-UTE-Cones sequence was applied to a hip agarose bone phantom and to 6 healthy volunteers for morphologic and quantitative T 2 ∗ and proton density mapping of trabecular bone.

RESULTS

Numeric simulation suggests that the IR technique with a short TR/TI combination provides sufficient suppression of long T2 tissues with a wide range of T1 s. High contrast imaging of trabecular bone can be achieved ex vivo and in vivo, with fitted T 2 ∗ values of 0.3-0.45 ms and proton densities of 5-9 mol/L.

CONCLUSION

The 3D IR-UTE-Cones sequence with a short TR/TI combination provides robust suppression of long T2 tissues and allows both selective imaging and quantitative ( T 2 ∗ and proton density) assessment of short T2 water components in trabecular bone in vivo.

Diffusion MRI for Assessment of Bone Quality; A Review of Findings in Healthy Aging and Osteoporosis

Diffusion MRI (dMRI) is a growing imaging technique with the potential to provide biomarkers of tissue variation, such as cellular density, tissue anisotropy, and microvascular perfusion. However, the role of dMRI in characterizing different aspects of bone quality, especially in aging and osteoporosis, has not yet been fully established, particularly in clinical applications. The reason lies in the complications accompanied with implementation of dMRI in assessment of human bone structure, in terms of acquisition and quantification. Bone is a composite tissue comprising different elements, each contributing to the overall quality and functional competence of bone. As diffusion is a critical biophysical process in biological tissues, early changes of tissue microstructure and function can affect diffusive properties of the tissue. While there are multiple MRI methods to detect variations of individual properties of bone quality due to aging and osteoporosis, dMRI has potential to serve as a superior method for characterizing different aspects of bone quality within the same framework but with higher sensitivity to early alterations. This is mainly because several properties of the tissue including directionality and anisotropy of trabecular bone and cell density can be collected using only dMRI. In this review article, we first describe components of human bone that can be potentially detected by their diffusivity properties and contribute to variations in bone quality during aging and osteoporosis. Then we discuss considerations and challenges of dMRI in bone imaging, current status, and suggestions for development of dMRI in research studies and clinics to segregate different contributing components of bone quality in an integrated acquisition. Level of Evidence: 5 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2020;51:975-992.

On the sensitivity of quantitative susceptibility mapping for measuring trabecular bone density

Purpose

To develop a methodological framework to simultaneously measure R2* and magnetic susceptibility in trabecularized yellow bone marrow and to investigate the sensitivity of Quantitative Susceptibility Mapping (QSM) for measuring trabecular bone density using a non‐UTE multi‐gradient echo sequence.

Methods

The ankle of 16 healthy volunteers and two patients was scanned using a time‐interleaved multi‐gradient‐echo (TIMGRE) sequence. After field mapping based on water–fat separation methods and background field removal based on the Laplacian boundary value method, three different QSM dipole inversion schemes were implemented. Mean susceptibility values in regions of different trabecular bone density in the calcaneus were compared to the corresponding values in the R2* maps, bone volume to total volume ratios (BV/TV) estimated from high resolution imaging (in 14 subjects), and CT attenuation (in two subjects). In addition, numerical simulations were performed in a simplified trabecular bone model of randomly positioned spherical bone inclusions to verify and compare the scaling of R2* and susceptibility with BV/TV.

Results

Differences in calcaneus trabecularization were well depicted in susceptibility maps, in good agreement with high‐resolution MR and CT images. Simulations and in vivo scans showed a linear relationship of measured susceptibility with BV/TV and R2*. The ankle in vivo results showed a strong linear correlation between susceptibility and R2* (R² = 0.88, p < 0.001) with a slope and intercept of −0.004 and 0.2 ppm, respectively.

Conclusions

A method for multi‐paramteric mapping, including R2*‐mapping and QSM was developed for measuring trabecularized yellow bone marrow, showing good sensitivity of QSM for measuring trabecular bone density.

The robustness of T2 value as a trabecular structural index at multiple spatial resolutions of 7 Tesla MRI

PURPOSE

To evaluate the robustness of MR transverse relaxation times of trabecular bone from spin-echo and gradient-echo acquisitions at multiple spatial resolutions of 7 T.

METHODS

The effects of MRI resolutions to T₂ and T2* of trabecular bone were numerically evaluated by Monte Carlo simulations. T₂ , T2*, and trabecular structural indices from multislice multi-echo and UTE acquisitions were measured in defatted human distal femoral condyles on a 7 T scanner. Reference structural indices were extracted from high-resolution microcomputed tomography images. For bovine knee trabecular samples with intact bone marrow, T₂ and T2* were measured by degrading spatial resolutions on a 7 T system.

RESULTS

In the defatted trabecular experiment, both T₂ and T2* values showed strong ( |r| > 0.80) correlations with trabecular spacing and number, at a high spatial resolution of 125 µm³ . The correlations for MR image-segmentation-derived structural indices were significantly degraded ( |r| < 0.50) at spatial resolutions of 250 and 500 µm³ . The correlations for T2* rapidly dropped ( |r| < 0.50) at a spatial resolution of 500 µm³ , whereas those for T₂ remained consistently high ( |r| > 0.85). In the bovine trabecular experiments with intact marrow, low-resolution (approximately 1 mm³ , 2 minutes) T₂ values did not shorten ( |r| > 0.95 with respect to approximately 0.4 mm³ , 11 minutes) and maintained consistent correlations ( |r| > 0.70) with respect to trabecular spacing (turbo spin echo, 22.5 minutes).

CONCLUSION

T₂ measurements of trabeculae at 7 T are robust with degrading spatial resolution and may be preferable in assessing trabecular spacing index with reduced scan time, when high-resolution 3D micro-MRI is difficult to obtain.

Bone mineral density in people living with HIV: a narrative review of the literature

Bone health status is largely absent in South Africa, the main reasons being the absence and cost-effectiveness of specific screening equipment for assessing bone mineral density (BMD). Various risk factors seem to play a role, some of which can be modified to change bone health status. Urbanisation is also a public health concern. Changing nutritional, as well as social behaviour, play integral roles in the prevalence and incidence of decreased BMD. Furthermore, human immunodeficiency virus (HIV) specifically, has a negative impact on BMD and although highly active antiretroviral therapy increases the prognosis for HIV-infected individuals, BMD still seem to decrease further. Dual energy X-ray absorptiometry is considered the gold standard for BMD assessment; however, recent developments have provided more cost-effective screening methods, among which heel quantitative ultrasound appears to be the most widely used in resource limited countries such as South Africa.

Simultaneous Quantification of Bone Edema/Adiposity and Structure in Inflamed Bone Using Chemical Shift-Encoded MRI in Spondyloarthritis

Purpose

To evaluate proton density fat fraction (PDFF) and R2* as markers of bone marrow composition and structure in inflamed bone in patients with spondyloarthritis.

Methods

Phantoms containing fat, water, and trabecular bone were constructed with proton density fat fraction (PDFF) and bone mineral density (BMD) values matching those expected in healthy bone marrow and disease states, and scanned using chemical shift‐encoded MRI (CSE‐MRI) at 3T. Measured PDFF and R2* values in phantoms were compared with reference FF and BMD values. Eight spondyloarthritis patients and 10 controls underwent CSE‐MRI of the sacroiliac joints. PDFF and R2* in areas of inflamed bone and fat metaplasia in patients were compared with normal bone marrow in controls.

Results

In phantoms, PDFF measurements were accurate over the full range of PDFF and BMD values. R2* measurements were positively associated with BMD but also were influenced by variations in PDFF. In patients, PDFF was reduced in areas of inflammation and increased in fat metaplasia compared to normal marrow. R2* measurements were significantly reduced in areas of fat metaplasia.

Conclusion

PDFF measurements reflect changes in marrow composition in areas of active inflammation and structural damage and could be used for disease monitoring in spondyloarthritis. R2* measurements may provide additional information bone mineral density but also are influenced by fat content. Magn Reson Med 79:1031–1042, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Evaluation of Spinal Bone Changes in Patients with Chronic Renal Failure by CT and Mr Imaging with Pathologic Correlation

To investigate bone changes in patients with chronic renal failure (CRF), bone mineral density (BMD) and T1 relaxation times were measured with CT and MR imaging and the results were correlated to histology. Excised lumbar vertebrae from 25 autopsy cases of CRF (18 males and 7 females), including 12 cases in which the patients had been receiving hemodialysis were examined. BMD and T1 relaxation time values were associated with specific histologic findings for cellularity, trabeculae, and peritrabecular fibrosis. Three vertebrae with low BMD showed increased hematopoietic marrow content, a finding not observed in primary osteoporosis. The vertebrae with osteosclerosis showed prolonged T1 relaxation time, which was due to increased amount of hematopoietic marrow, and the presence of thickened or many small irregular trabeculae or peritrabecular fibrosis. These findings may be useful in the evaluation of bone changes in patients with CRF.

Bone marrow segmentation based on a combined consideration of transverse relaxation processes and Dixon oscillations

The aim of this study was to demonstrate that gradient-echo sampling of single spin echoes can be used to isolate the signal from trabecular bone marrow, with high-quality segmentation and surface reconstructions resulting from the application of simple post-processing strategies. Theoretical expressions of the time-domain single-spin-echo signal were used to simulate signals from bone marrow, non-bone fatty deposits and muscle. These simulations were compared with and used to interpret signals obtained by the application of the gradient-echo sampling of a spin-echo sequence to image the knee and surrounding tissues at 1.5 T. Trabecular bone marrow has a much higher reversible transverse relaxation rate than surrounding non-bone fatty deposits and other musculoskeletal tissues. This observation, combined with a choice of gradient-echo spacing that accentuates Dixon-type oscillations from chemical-shift interference effects, enabled the isolation of bone marrow signal from surrounding tissues through the use of simple image subtraction and thresholding. Three-dimensional renderings of the marrow surface were then readily generated with this approach - renderings that may prove useful for bone morphology assessment, e.g. for the measurement of femoral anteversion. In conclusion, understanding the behavior of signals from bone marrow and surrounding tissue as a function of time through a spin echo facilitates the segmentation and reconstruction of bone marrow surfaces using straightforward post-processing strategies that are typically available on modern radiology workstations.

Focus on diffusion MR investigations of musculoskeletal tissue to improve osteoporosis diagnosis: a brief practical review

Nowadays, a huge number of papers have documented the ability of diffusion magnetic resonance imaging (D-MRI) to highlight normal and pathological conditions in a variety of cerebral, abdominal, and cardiovascular applications. To date, however, the role of D-MRI to investigate musculoskeletal tissue, specifically the cancellous bone, has not been extensively explored. In order to determine potentially useful applications of diffusion techniques in musculoskeletal investigation, D-MRI applications to detect osteoporosis disease were reviewed and further explained.

Validation of bone marrow fat quantification in the presence of trabecular bone using MRI

BACKGROUND

To validate six-echo, chemical-shift based MRI with T2 * correction for the quantification of bone marrow fat content in the presence of trabecular bone.

METHODS

Ten bone phantoms were made using trabecular bone cores extracted from the distal femur and proximal tibia of 20 human cadaveric knees. Bone marrow was removed from the cores and the marrow spaces were filled with water-fat gelatin to mimic bone marrow of known fat fractions. A chemical-shift based water-fat separation method with T2 * correction was used to generate fat fraction maps. The proton density fat fractions (PDFF) between marrow regions with and without bone were compared with the reference standard of known fat fraction using the squared Pearson correlation coefficient and unpaired t-test.

RESULTS

Strong correlations were found between the known fat fraction and measured PDFF in marrow without trabecular bone (R(2) = 0.99; slope = 0.99, intercept = 0.94) as well as in marrow with trabecular bone (R(2) = 0.97; slope = 1.0, intercept = -3.58). Measured PDFF between regions with and without bone were not significantly different (P = 0.5). However, PDFF was systematically underestimated by -3.2% fat fraction in regions containing trabecular bone.

CONCLUSION

Our implementation of a six-echo chemical-shift based MRI pulse sequence with T2 * correction provided an accurate means of determining fat content in bone marrow in the presence of trabecular bone.

T1/T2*-weighted MRI provides clinically relevant pseudo-CT density data for the pelvic bones in MRI-only based radiotherapy treatment planning

BACKGROUND AND PURPOSE

In radiotherapy (RT), target soft tissues are best defined on MR images. In several cases, CT imaging is needed only for dose calculation and generation of digitally reconstructed radiographs (DRRs). Image co-registration errors between MRI and CT can be avoided by using MRI-only based treatment planning, especially in the pelvis. Since electron density information can not be directly derived from the MRI, a method is needed to convert MRI data into CT like data. We investigated whether there is a relationship between MRI intensity and Hounsfield unit (HU) values for the pelvic bones. The aim was to generate a method to convert bone MRI intensity into HU data surrogate for RT treatment planning.

MATERIAL AND METHODS

The HU conversion model was generated for 10 randomly chosen prostate cancer patients and independent validation was performed in another 10 patients. Data consisted of 800 image voxels chosen within the pelvic bones in both T1/T2*-weighted gradient echo and CT images. Relation between MRI intensity and electron density was derived from calibrated HU-values. The proposed method was tested by constructing five "pseudo"-CT series.

RESULTS

We found that the MRI intensity is related to the HU value within a HU range from 0 to 1400 within the pelvic bones. The mean prediction error of the conversion model was 135 HU. Dose calculation based on the pseudo-CT images was accurate and the generated DRRs were of good quality.

CONCLUSIONS

The proposed method enables generation of clinically relevant pseudo-CT data for the pelvic bones from one MRI series. It is simpler than previously reported approaches which require either acquisition of several MRI series or T2* maps with special imaging sequences. The method can be applied with commercial clinical image processing software. The application requires segmentation of the bones in the MR images.

Noninvasive imaging of bone microarchitecture

The noninvasive quantification of peripheral compartment-specific bone microarchitecture is feasible with high-resolution peripheral quantitative computed tomography (HR-pQCT) and high-resolution magnetic resonance imaging (HR-MRI). In addition to classic morphometric indices, both techniques provide a suitable basis for virtual biomechanical testing using finite element (FE) analyses. Methodical limitations, morphometric parameter definition, and motion artifacts have to be considered to achieve optimal data interpretation from imaging studies. With increasing availability of in vivo high-resolution bone imaging techniques, special emphasis should be put on quality control including multicenter, cross-site validations. Importantly, conclusions from interventional studies investigating the effects of antiosteoporotic drugs on bone microarchitecture should be drawn with care, ideally involving imaging scientists, translational researchers, and clinicians.

Quantitative analysis of the diffusion-weighted steady-state free precession signal in vertebral bone marrow lesions

OBJECTIVES

: Diffusion-weighted steady-state free precession (DW-SSFP) sequences have shown great potential for the differential diagnosis of benign osteoporotic and malignant neoplastic vertebral compression fractures, which appear hypo- to isointense or hyperintense in DW-SSFP magnetic resonance imaging, respectively. In contrast to other diffusion weighting sequences, the DW-SSFP signal depends not only on the apparent diffusion coefficient (ADC), but also on the tissue relaxation times and sequence parameters. The purpose of the present study was to provide a detailed analysis of the DW-SSFP signal in benign and malignant vertebral lesions (VLs) and in vertebral bone marrow (VBM) to understand the observed signal alterations and their dependence on tissue and sequence parameters.

MATERIALS AND METHODS

: Magnetic resonance imaging was performed in 40 patients with benign (n = 20) or malignant (n = 20) VLs to determine the fat fraction and tissue parameters (ADC, T1, T2, T2*) for both the water and fat signal. With these values, the DW-SSFP signal was simulated and compared with the measured signals for different diffusion gradients by determining the signal intensity ratio between the SSFP signals of the lesions and of normal-appearing VBM for both malignant and benign VLs.

RESULTS

: The simulated DW-SSFP contrast agreed well with the measured contrast and provided a very good differentiation between benign osteoporotic and malignant VLs. ADCs were significantly different in both lesion types (malignant 1.36 vs. osteoporotic 1.77 × 10 mm/s); however, the observed contrast differences were caused predominantly by an opposed-phase readout in combination with significantly different T2* values (malignant 22 vs. osteoporotic 14 ms) and fat fractions (malignant 3.9% vs. osteoporotic 12%) in the lesions as well as significantly different fat fractions in normal-appearing VBM (malignant 42% vs. osteoporotic 52%) of both patient groups.

CONCLUSIONS

: Although the ADCs of the evaluated malignant and benign VLs showed highly significant differences, the influence of diffusion on the DW-SSFP signal contrast is relatively low compared with other tissue parameters due to the very complex signal mechanism of the SSFP sequence. Thus, the observed DW-SSFP signal contrast of different VLs (hypo-/isointense vs. hyperintense signal) is rather fat- and T2*-weighted than diffusion-weighted. The intermediate diffusion weighting of the applied SSFP sequence, however, helps to shift the different contrasts into a signal range that is easily visually accessible.

Correlation of signal attenuation-based quantitative magnetic resonance imaging with quantitative computed tomographic measurements of subchondral bone mineral density in metacarpophalangeal joints of horses

OBJECTIVE

To evaluate the ability of signal attenuation-based quantitative magnetic resonance imaging (QMRI) to estimate subchondral bone mineral density (BMD) as assessed via quantitative computed tomography (QCT) in osteoarthritic joints of horses.

SAMPLE POPULATION

20 metacarpophalangeal joints from 10 horse cadavers.

PROCEDURES

Magnetic resonance (MR) images (dorsal and transverse T1-weighted gradient recalled echo [GRE] and dorsal T2*-weighted GRE fast imaging employing steady-state acquisition [T2*-FIESTA]) and transverse single-slice computed tomographic (CT) images of the joints were acquired. Magnetic resonance signal intensity (SI) and CT attenuation were quantified in 6 regions of interest (ROIs) in the subchondral bone of third metacarpal condyles. Separate ROIs were established in the air close to the joint and used to generate corrected ratios and SIs. Computed tomographic attenuation was corrected by use of a calibration phantom to obtain a K(2)HPO(4)-equivalent density of bone. Correlations between QMRI performed with different MR imaging sequences and QCT measurements were evaluated. The intraobserver repeatability of ROI measurements was tested for each modality.

RESULTS

Measurement repeatability was excellent for QCT (R(2) = 98.3%) and QMRI (R(2) = 98.8%). Transverse (R(2) = 77%) or dorsal (R(2) = 77%) T1-weighted GRE and QCT BMD measurements were negatively correlated, as were dorsal T2*-FIESTA and QCT (R(2) = 80%) measurements. Decreased bone SI during MR imaging linearly reflected increased BMD.

CONCLUSIONS AND CLINICAL RELEVANCE

Results of this ex vivo study suggested that signal attenuation-based QMRI was a reliable, clinically applicable method for indirect estimation of subchondral BMD in osteoarthritic metacarpophalangeal joints of horses.

Imaging of the musculoskeletal system in vivo using ultra-high field magnetic resonance at 7 T

Recently, great progress has been made in particularly in the imaging of cartilage and bone structure. Increased interest has focused on high-field (3 Tesla) imaging and more recently on ultra-high field (UHF) magnetic resonance imaging (MRI) at 7 T for in vivo imaging. Because the signal-to-noise ratio (SNR) scales linearly with field strength, a substantial increase in SNR is expected compared with lower field strengths. This gain in SNR can be used to increase spatial resolution or reduce imaging time. The goal of this review was to highlight recent developments and challenges in in vivo musculoskeletal (MSK) imaging using UHF-MRI at 7 T. One focus of this review is on the emerging methodology of quantitative MRI for the assessment of trabecular bone structure at the tibia, wrist, and knee. In particular for this application, susceptibility effects between the bone and bone marrow transitions that scale with field strength have to be considered. Another important MSK application is the characterization of knee cartilage morphology. The higher SNR provided by UHF-MRI is a potential advantage for visualizing, segmenting, and analyzing cartilage. Standard clinical MSK imaging relies heavily on T1, T2, and proton density weighted fast spin echo sequences. However, fast spin echo imaging has proven to be very challenging at higher fields because of very high specific absorption rates, using multiple pulses in a short time frame; thus the imaging protocols have to be adapted and gradient echo sequences may be more beneficial. Imaging of more central body parts such as the spine at 7 T is still in its infancy and dedicated coils have to be developed.

High-resolution MRI of internal field diffusion-weighting in trabecular bone

Trabecular bone structure is known to play a crucial role in the overall strength, and thus fracture risk, of such areas of the skeleton as the vertebrae, spine, femur, tibiae, or radius. Several MR methods devoted to probing this structure depend upon the susceptibility difference between the solid bone matrix and the intervening fluid/marrow/fat, usually in the context of a linewidth (1/T(2)') measurement or mapping technique. A recently demonstrated new approach to this system involves using internal gradients to encode diffusion weighting, and extracting structural information (e.g., surface-to-volume ratio) from the resulting signal decay. This contrast method has been demonstrated in bulk measurements on cleaned, water-saturated bovine trabecular bone samples. In the present work, microscopic imaging (0.156 mm in-plane resolution) is performed in order to spatially resolve this contrast on the trabecular level, and confirm its interpretation for the bulk measurements. It is found that the local rate of decay due to diffusion in the internal field (DDIF) is maximal close to the trabecular surfaces. The overall decay rate in a lower resolution scan probes the abundance of these surfaces, and provides contrast beyond that found in conventional proton density weighted or T(1)-weighted imaging. Furthermore, a microscopic calculation of internal field distributions shows a qualitative distinction between the structural sensitivities of DDIF and T(2)'. DDIF contrast is highly localized around trabecular walls than is the internal field itself, making it a less sensitive but more specific measure of such important properties as trabecular number.

Early detection and quantification of murine melanoma bone metastases with magnetic resonance imaging

OBJECTIVE

Bone metastases occur in approximately 80% of patients with advanced cancer and cause significant morbidity. There are currently no established means by which to identify the early growth of micro-metastatic cells or their effects on bone at a time when curative therapy might be initiated. We postulated that high-resolution magnetic resonance imaging (MRI) could detect and quantify the growth and destructive effects of bone micrometastases.

DESIGN

Using a mouse model for metastasis of malignant melanoma, we have examined the ability of MRI to quantify cortical bone destruction and the percentage of the medullary cavity occupied by tumour, trabecular bone, and marrow. The results from MRI were compared to histomorphometry (the reference standard) and to radiographs.

RESULTS

In vivo gradient-echo and spin-echo MRI demonstrated that metastatic melanoma replaced the marrow space but that the cortical bone integrity was preserved (P < or = 0.001). The smallest detectable micrometastasis had an area of 0.323 mm(2). In contrast, we observed no trends after quantifying the radiograph data.

CONCLUSION

These approaches delineated the limits of MRI in its ability to quantify tumour burden and the effect on bone in this model. Given the increasing use of MRI as a non-invasive clinical diagnostic method, the present findings may be applicable in detecting bone metastases in the clinical setting at an early and potentially treatable stage.

In vivo quantification of response to treatment in patients with multiple myeloma by 1H magnetic resonance spectroscopy of bone marrow

OBJECT

Magnetic resonance imaging (MRI) is the gold standard non-invasive technique to detect malignant disease in the bone marrow. Proton magnetic resonance spectroscopy (MRS) can be performed as a quick adjunct to routine spinal MRI. We performed proton MRS to patients with multiple myeloma (MM) at diagnosis and after treatment to investigate the possible correlation of MRS data with response to therapy.

PATIENTS AND METHODS

Twenty-one patients with newly diagnosed MM underwent combined MRI/MRS explorations of a transverse center section in the fifth lumbar vertebral body. MRS was acquired with STEAM and 40 ms TE. Areas of unsuppressed water and lipid resonances were used to calculate the lipid-to-water ratio (LWR).

RESULTS

No association was detected between initial LWRs and the clinical characteristics of patients. Post treatment MRS was available in 16 patients of whom 11 (69%) presented an LWR increase, this included all complete responders (8/8, 100%, P = 0.012). A post-treatment LWR value equal to or larger than one is proposed as a non-invasive marker of complete response to treatment.

CONCLUSION

Only patients responding to treatment presented a significant increase in bone marrow LWR after therapy. MRS may provide an adequate quantification of response to chemotherapy in patients with MM.

Structural and functional assessment of trabecular and cortical bone by micro magnetic resonance imaging

Osteoporosis is a multifactorial disorder of bone mineral homeostasis affecting the elderly. It is a major public health issue with significant socioeconomic consequences. Recent findings suggest that bone loss-the key manifestation of the disease-is accompanied by architectural deterioration, both affecting the bone's mechanical competence and susceptibility to fracture. This article reviews the potential of quantitative micro MRI (mu-MRI), including a discussion of the technical requirements for image acquisition, processing, and analysis for assessing the architectural implications of osteoporosis and as a means to monitor the response to treatment. With current technology, the resolution achievable in clinically acceptable scan times and necessary signal-to-noise ratio (SNR) is comparable to trabecular thickness. This limited spatial resolution regime demands processing and analysis algorithms designed to operate under such limiting conditions. It is shown that three different classes of structural parameters can be distinguished, characterizing scale, topology, and orientation. There is considerable evidence that osteoporotic bone loss affects all three classes but that topological changes, resulting from conversion of trabecular plates to rods, with the latter's eventual disconnection, are particularly prominent. Clinical applications discussed can be divided into those dealing with assessment of osteoporotic fracture risk as opposed to the study of the effect of disease progression and regression in response to treatment. Current data suggest that noninvasive assessment of cortical and trabecular bone (TB) architecture by mu-MRI may provide new surrogate endpoints to assess the efficacy of intervention in osteoporosis treatment and prevention.

Bone marrow changes in beta-thalassemia major: quantitative MR imaging findings and correlation with iron stores

The purpose of this study is to describe the MR imaging features of bone marrow in beta-thalassemia major and investigate their relation to ferritin, liver and spleen siderosis. Spinal bone marrow was prospectively assessed on abdominal MR studies of 40 transfused beta-thalassemic patients and 15 controls using T1-w, Pd, T2*-w Gradient Echo (GRE) and T1-w turbo Spin Echo (TSE) sequences. Signal intensity (SI) ratios of liver, spleen and bone marrow to paraspinous muscles (L/M, S/M, B/M respectively) and the respective T2 relaxation rates (1/T2) were calculated. Serum ferritin levels were recorded. Bone marrow hypointensity in at least T2*-w GRE sequence was noted in 29/40 (72.5%) patients. Eleven/40 patients exhibited normal B/M on all MR sequences. Five/40 patients had normal B/M and low L/M. B/M correlated with L/M in T1-w TSE sequence only (r = 0.471, p = 0.05). B/M correlated with S/M and mean ferritin values in all sequences (r > 0.489, p < 0.01 and r > - 0.496, p < 0.03 respectively). Marrow 1/T2 did not correlate with ferritin values or liver and spleen 1/T2. B/M in transfused beta-thalassemic patients is related to splenic siderosis and ferritin levels. Although marrow is usually hypointense, it may occasionally display normal SI coexisting with liver hypointensity, a pattern typical of primary hemochromatosis.

Prediction of mechanical properties of trabecular bone using quantitative MRI

Techniques for quantitative magnetic resonance imaging (MRI) have been developed for non-invasive estimation of the mineral density and structure of trabecular bone. The R*(2) relaxation rate (i.e. 1/T*(2)) is sensitive to bone mineral density (BMD) via susceptibility differences between trabeculae and bone marrow, and by binarizing MRI images, structural variables, such as apparent bone volume fraction, can be assessed. In the present study, trabecular bone samples of human patellae were investigated in vitro at 1.5 T to determine the ability of MRI-derived variables (R*(2) and bone volume fraction) to predict the mechanical properties (Young's modulus, yield stress and ultimate strength). Further, the MRI variables were correlated with reference measurements of volumetric BMD and bone area fraction as determined with a clinical pQCT system. The MRI variables correlated significantly (p < 0.01) with the mechanical variables (r = 0.32-0.46), BMD (r = 0.56) and bone structure (r = 0.51). A combination of R*(2) and MRI-derived bone volume fraction further improved the prediction of yield stress and ultimate strength. Although pQCT showed a trend towards better prediction of the mechanical properties, current results demonstrate the feasibility of combined MR imaging of marrow susceptibility and bone volume fraction in predicting the mechanical strength of trabecular bone and bone mineral density.

Quantitative MRI for the assessment of bone structure and function

Osteoporosis is the most common degenerative disease in the elderly. It is characterized by low bone mass and structural deterioration of bone tissue, leading to morbidity and increased fracture risk in the hip, spine and wrist-all sites of predominantly trabecular bone. Bone densitometry, currently the standard methodology for diagnosis and treatment monitoring, has significant limitations in that it cannot provide information on the structural manifestations of the disease. Recent advances in imaging, in particular MRI, can now provide detailed insight into the architectural consequences of disease progression and regression in response to treatment. The focus of this review is on the emerging methodology of quantitative MRI for the assessment of structure and function of trabecular bone. During the past 10 years, various approaches have been explored for obtaining image-based quantitative information on trabecular architecture. Indirect methods that do not require resolution on the scale of individual trabeculae and therefore can be practiced at any skeletal location, make use of the induced magnetic fields in the intertrabecular space. These fields, which have their origin in the greater diamagnetism of bone relative to surrounding marrow, can be measured in various ways, most typically in the form of R2', the recoverable component of the total transverse relaxation rate. Alternatively, the trabecular network can be quantified by high-resolution MRI (micro-MRI), which requires resolution adequate to at least partially resolve individual trabeculae. Micro-MRI-based structure analysis is therefore technically demanding in terms of image acquisition and algorithms needed to extract the structural information under conditions of limited signal-to-noise ratio and resolution. Other requirements that must be met include motion correction and image registration, both critical for achieving the reproducibility needed in repeat studies. Key clinical applications targeted involve fracture risk prediction and evaluation of the effect of therapeutic intervention.

Fully balanced steady-state 3D-spin-echo (bSSSE) imaging at 3 Tesla

In this work a 3D fully balanced steady-state spin-echo sequence (bSSSE) was designed and implemented on a 3-Tesla whole-body MR scanner. In contrast to previously introduced nonbalanced SSSE (nbSSSE) sequences, bSSSE features fully rewound imaging gradients. Additionally, all phase encodings and the readout prephasing are performed after the refocusing pulse. Its performance was compared with previously used gradient-echo (GE) and nbSSSE sequences. The signal response of each sequence was simulated by numerical solution of the Bloch equation. All sequences were then optimized based on these simulations for high-resolution (HR) imaging of tissues with short relaxation times and applied to in vivo HR-MRI of trabecular bone. bSSSE outperformed nbSSSE sequences in terms of SNR and SNR efficiency. Measurements of structural bone parameters revealed thinner trabeculae and a lower bone/bone-marrow fraction for the bSSSE sequence compared to the nbSSSE sequence. By applying a parallel imaging technique based on generalized autocalibrating partially parallel acquisition (GRAPPA) with a reduction factor of 2, we were able to maintain a clinically feasible scan time using bSSSE. Comparisons of structural bone parameters revealed a difference of less than 3% for all structural parameters between parallel and conventional imaging.

Implications of pulse sequence in structural imaging of trabecular bone

PURPOSE

To investigate the SNR and image properties of 3D steady-state free precession (SSFP), fast large-angle spin echo (FLASE), gradient-recalled acquisition in steady state (GRASS), and spoiled GRASS (SPGR) for structural imaging of trabecular bone (TB).

MATERIALS AND METHODS

SNR was examined theoretically and experimentally on phantoms, bone specimens, and in vivo. The bone volume fraction, TB thickness, and echo time (TE) dependence of the thickness were compared. The trabecula was modeled as a cylinder in simulations to examine the intra-voxel spin-dephasing in SSFP and GRASS. Images were acquired on a 1.5 T Siemens Sonata system (40 mT/m maximum gradient, 200 T/m/s peak slew rate).

RESULTS

Within the hardware and safety limit constraints, SNR of FLASE was superior, followed by SSFP, GRASS, and SPGR. The trabecular thickness derived from gradient-echo images was 10-45% greater than that obtained with FLASE. Conversely, SSFP images delineated partial volume trabeculae better than FLASE. Simulations indicated that the artifactual thickening was more severe in SSFP than in GRASS, which was attributed to off-resonance effects from the induced gradients at the bone/marrow interface.

CONCLUSION

FLASE had the highest SNR and was insensitive to susceptibility dephasing. Although SSFP has superior SNR compared to GRASS, off-resonance effects and duty cycle limitations may compromise its practicality in this application. Inc.

Correlation between the bone marrow blood perfusion and lipid water content on the lumbar spine in female subjects

PURPOSE

To assess the marrow lipid water ratio (LWR) and spectral line width (LW) of lumbar vertebra by proton MR spectroscopy ((1)H MRS) and correlate with the marrow blood perfusion (MBP) by contrast-enhanced dynamic MRI (dMRI) in female subjects.

MATERIALS AND METHODS

A total of 50 female subjects were included. Single-voxel (1)H MRS was measured at the L3 vertebrae. In the dMRI study, the maximum enhancement percentage (E(max)) and enhancement slope were calculated from time-intensity curves. Pearson's correlation was calculated to explore the association of (1)H MRS, age, and perfusion parameters. Partial correlation was then used to control confounder effect.

RESULTS

LWR was negatively correlated with E(max) (r = -0.72; P < 0.0001), but positively correlated with age (r = 0.63; P < 0.0001). Lipid LW was negatively correlated with age (r = -0.33; P <0.05), but positively correlated with E(max) (r = 0.4; P < 0.05). After adjusting for the age effect by partial correlation, marrow LWR still negatively correlated with E(max) (r = -0.63; P < 0.0001). After adjusting for the E(max) effect, marrow LWR still positively correlated with age (r = 0.37; P < 0.05).

CONCLUSION

After controlling the age factor, the LWR and water LW has significant negative correlation with vertebral BMP. Further investigation is needed to explore the relationship between bone marrow fat metabolism and angiogenesis.

Assessment of compressive modulus, hydraulic permeability and matrix content of trypsin-treated nucleus pulposus using quantitative MRI

A clinical strength MRI and intact bovine caudal intervertebral discs were used to test the hypotheses that (1) mechanical loading and trypsin treatment induce changes in NMR parameters, mechanical properties and biochemical contents; and (2) mechanical properties are quantitatively related to NMR parameters. MRI acquisitions, confined compression stress-relaxation experiments, and biochemical assays were applied to determine the NMR parameters (relaxation times T1 and T2, magnetization transfer ratio (MTR) and diffusion trace (TrD)), mechanical properties (compressive modulus H(A0) and hydraulic permeability k(0)), and biochemical contents (H(2)O, proteoglycan and total collagen) of nucleus pulposus tissue from bovine caudal discs subjected to one of two injections and one of two mechanical loading conditions. Significant correlations were found between k(0) and T1 (r=0.75,p=0.03), T2 (r=0.78, p=0.02), and TrD (r=0.85, p=0.007). A trend was found between H(A0) and TrD (r=0.56, p=0.12). However, loading decreased these correlations (r=0.4, p=0.2). The significant effect of trypsin treatment on mechanical properties, but not on NMR parameters, may suggest that mechanical properties are more sensitive to the structural changes induced by trypsin treatment. The significant effect of loading on T1 and T2, but not on H(A0) or k(0), may suggest that NMR parameters are more sensitive to the changes in water content enhanced by loading. We conclude that MRI offers promise as a sensitive and non-invasive technique for describing alterations in material properties of intervertebral disc nucleus, and our results demonstrate that the hydraulic permeability correlated more strongly to the quantitative NMR parameters than did the compressive modulus; however, more studies are necessary to more precisely characterize these relationships.

Trabecular bone volume fraction measurements of a large number of subjects using a compact MRI

Trabecular bone volume fraction (TBVF) and speed of sound (SOS) were measured for the right calcanei of 416 female volunteers. The TBVF was measured with a compact MRI developed in our laboratory. The SOS was measured with a commercial quantitative ultrasound system. It was observed that the correlation coefficient between TBVF and SOS and that between TBVF and age varied depending on the location of region of interest (ROI) in the calcaneus. As a result, an optimum circular ROI with a diameter of 20 mm was determined so that the correlation coefficients were maximized. In the optimum ROI, transverse relaxation time (T2) of the bone marrow protons of the calcaneus was found to be concentrated in a narrow range over the subjects. This result suggested that a 50% scan time reduction in the TBVF measurements could be made by skipping the T2 correction procedure.

Magnetic resonance imaging of trabecular and cortical bone in mice: comparison of high resolution in vivo and ex vivo MR images with corresponding histology

Measurements of bone morphometry and remodeling have been shown to reflect bone strength and can be used to diagnose degenerative bone disease. In this study, in vivo and ex vivo magnetic resonance imaging (MRI) techniques to assess trabecular and cortical bone properties have been compared to each other and to histology as a novel means for the quantification of bone. Femurs of C57Bl/6 mice were examined both in vivo and ex vivo on an 11.7 T MRI scanner, followed by histologic processing and morphometry. A thresholding analysis technique was applied to the MRI images to generate contour lines and to delineate the boundaries between bone and marrow. Using MRI, an optimal correlation with histology was obtained with an in vivo longitudinal sectioned short echo time gradient-echo versus an in vivo long echo time spin-echo sequence or an ex vivo pulse sequence. Gradient-echo images were acquired with a maximum in-plane resolution of 35 microm. Our results demonstrated that in both the in vivo and ex vivo data sets, the percent area of marrow increases and percent area of trabecular bone and cortical bone thickness decreases moving from the epiphyseal growth plate to the diaphysis. These changes, observed with MRI, correlate with the histological data. Investigations using in vivo MRI gradient-echo sequences consistently gave the best correlation with histology. Our quantitative evaluation using both ex vivo and in vivo MRI was found to be an effective means to visualize non-invasively the normal variation in trabecular and cortical bone as compared to a histological "gold standard" The experiments validated in vivo MRI as a potential high resolution technique for investigating both soft tissue, such as marrow, and bone without radiation exposure.

Correlation of bone marrow lipid water content with bone mineral density on the lumbar spine

STUDY DESIGN

Prospective study.

OBJECTIVES

To assess the proton MR spectroscopy (1H MRS) of vertebral bone marrow and correlate the lipid water ratio (LWR) and spectral line width (LW) with bone mineral density (BMD) in female subjects.

SUMMARY OF BACKGROUND DATA

The mechanism of bone marrow fat accumulation and bone mineral content is poorly understood. Proton MR spectroscopy was used to demonstrate the lipid and water spectra in the bone marrow. We try to assess the possible interaction between the bone marrow lipid content, aging, and BMD.

METHODS

Proton MRS and BMD of the lumbar spine were performed in 52 female subjects (mean age, 58 years; SD, 10 years). They were 13 premenopausal and 39 postmenopausal women. The BMD (g/cm2) was measured using dual energy radiograph absorptiometry at the lumbar spine. Single voxel 1H MRS was measured at L3 vertebral body by stimulated echo-acquisition mode (STEAM) sequence and demonstrated two major peaks (lipid and water). Comparisons of the differences between the two subgroups were made. Pearson's correlation was also calculated to explore the association of 1H MRS measurements with age and BMD. Partial correlation was further conducted when controlling the variable such as age or BMD.

RESULTS

BMD and LWR had statistically significant difference between the pre- and postmenopausal subgroups (P < 0.001), while lipid LW had a borderline difference and water LW had no difference. LWR was positively correlated with age (r = 0.52 and P < 0.0001) and negatively correlated with BMD (r = -0.40 and P = 0.003) for all the subjects. Lipid LW was negatively correlated with age (r = -0.32 and P = 0.0197) and positively correlated with BMD (r = 0.67 and P < 0.0001). When controlling for BMD effect, only LWR is statistically correlated with age (partial r = 0.39, P = 0.0045), while only the lipid LW is statistically correlated with BMD when controlling for age (partial r = 0.63, P < 0.0001). None of the correlations between water LW and age or BMD was significant. In the subgroups, only the lipid LW is significantly correlated with BMD (r = 0.78, P = 0.0016 in premenopausal women; r = 0.62, P < 0.0001 in postmenopausal women).

CONCLUSION

The LWR had a positive correlation with the age, while the lipid LW had a positive correlation with BMD, even after controlling the age factor. The bone marrow lipid water content and metabolism acted as important roles in the internal environment of bone and influenced bone mineralization.

Optimisation of T2 and M0 measurements of bi-exponential systems

Cramer-Rao theory and computer simulations were used to show that the errors involved in calculating the magnetization and relaxation parameters of a two-component system decrease with: (1) increasing SNR, (2) increasing number of echoes used in the fitting procedure, and (3) increasing ratio of the relaxation times of the two components, T(22)/T(21). Images of bi-compartmental phantoms of known T(2) values were acquired using an optimized imaging sequence, and an optimized fitting algorithm was used to calculate the T(2) values of the two components by fitting the resulting images to a bi-exponential decay model. Accuracy better than 6% was achieved in the calculations of the T(2) values of the two components, and region fitting provided better accuracy than pixel-by-pixel fitting. The procedures were used to calculate the T(2) and M(0) values of equine carpal bones with known degree of radiographic bone sclerosis. Although the T(2) and M(0) values of both water and fat components all decreased with the degree of radiographic bone sclerosis, the transverse relaxation of the water component, T(2W), showed a greater decrease with advanced stages of bone sclerosis.

In-vivo high resolution three-dimensional MRI studies of rat joints at 7 T

It is important to obtain high resolution images of joints for the study of disease, especially in rodent experimental models. We optimized (1)H magnetic resonance imaging three-dimensional sequences at 7 T, with lipid signal suppression, and T(1) and T(2) measurements for in-vivo experiments on rat joints, in order to assess the effectiveness of high-field MRI. The method was validated by applying it to the early diagnosis of arthritis. We studied the progress of rheumatoid arthritis in an arthritic rat model. We observed the rats' knees for 21 days after inducing arthritis. The images acquired over one hour had a high resolution of 1.75 x 10(-3) mm(3), (105 x 105 x 145 microm(3)) which allowed us to spot the early stages of joint degeneration, such as bone erosion, and to observe an apparent 'MRI' loss of cartilage thickness, attributed to dehydration of the cartilage tissue. The MR images obtained during the early stages of rheumatoid arthritis enabled us to study joint changes accurately before any histological signs of attack were visible.

Investigation of the microstructure of the isolated rat heart: A comparison betweenT*2- and diffusion-weighted MRI

Myocardial fiber structure can be determined with diffusion-weighted (DW) MRI as well as with high-resolution T*(2) imaging. The purpose of the present study was twofold: to provide a more quantitative description of T*(2)-based myocardial fiber contrast, and to compare the T*(2)-based fiber structure with high-resolution (78 microm in-plane, 1-mm slice thickness) DW images of the isolated rat heart at 11.75 T. This study demonstrates that the static dephasing regime is responsible for visualization of myocardial microstructure, and that the dynamic dephasing regime can be neglected. In comparison with DW experiments, T*(2) mapping and DW images yield almost equivalent information on myocardial fiber structure.

Osteoporosis imaging

Because osteoporotic fractures may be prevented, diagnostic techniques are essential in the assessment of osteoporosis. Conventional radiographs of the spine are not suited for diagnosing early osteoporosis, but they show fractures that may have no clinical symptoms. The radiologist should be aware of the enormous significance of these fractures for future osteoporotic fractures. Bone mass measurements are standard techniques in the diagnosis of osteoporosis, which are the basis of the WHO definition of osteoporosis. In this article the authors presented these standard techniques and newer diagnostic techniques that provide insights in the structure of trabecular bone.