B1-insensitive T2 mapping of healthy thigh muscles using a T2-prepared 3D TSE sequence

Insights into the Age Dependency of Compositional MR Biomarkers Quantifying the Health Status of Cartilage in Metacarpophalangeal Joints

(1) Background: We aim to investigate age-related changes in cartilage structure and composition in the metacarpophalangeal (MCP) joints using magnetic resonance (MR) biomarkers. (2) Methods: The cartilage tissue of 90 MCP joints from 30 volunteers without any signs of destruction or inflammation was examined using T1, T2, and T1ρ compositional MR imaging techniques on a 3 Tesla clinical scanner and correlated with age. (3) Results: The T1ρ and T2 relaxation times showed a significant correlation with age (T1ρ: Kendall-τ-b = 0.3, p < 0.001; T2: Kendall-τ-b = 0.2, p = 0.01). No significant correlation was observed for T1 as a function of age (T1: Kendall-τ-b = 0.12, p = 0.13). (4) Conclusions: Our data show an increase in T1ρ and T2 relaxation times with age. We hypothesize that this increase is due to age-related changes in cartilage structure and composition. In future examinations of cartilage using compositional MRI, especially T1ρ and T2 techniques, e.g., in patients with osteoarthritis or rheumatoid arthritis, the age of the patients should be taken into account.

Beyond mean value analysis - a voxel-based analysis of the quantitative MR biomarker water T2 in the presence of fatty infiltration in skeletal muscle tissue of patients with neuromuscular diseases

The main pathologies in the muscles of patients with neuromuscular diseases (NMD) are fatty infiltration and edema. Recently, quantitative magnetic resonance (MR) imaging for determination of the MR biomarkers proton density fat fraction (PDFF) and water T₂ (T_2w ) has been advanced. Biophysical effects or pathology can have different effects on MR biomarkers. Thus, for heterogeneously affected muscles, the routinely performed mean or median value analyses of MR biomarkers are questionable. Our work presents a voxel-based histogram analysis of PDFF and T_2w images to point out potential quantification errors. In 12 patients with NMD, chemical-shift encoding-based water-fat imaging for PDFF and T₂ mapping with spectral adiabatic inversion recovery (SPAIR) for T_2w determination was performed. Segmentation of nine thigh muscles was performed bilaterally (n = 216). PDFF and T₂ maps were coregistered. A voxel-based comparison of PDFF and T_2w showed a decreased T_2w with increasing PDFF. Mean T_2w and mean T_2w ^{without fatty voxels} (PDFF < 10%) show good agreement, whereas standard deviation (σ) T_2w and σ T_2w ^{without fatty voxels} show increasing difference with increasing values of σ. Thereby two subgroups can be observed, referring to muscles in which the exclusion of fatty voxels has a negligible influence versus muscles in which a strong dependency of the T_2w value distribution on the exclusion of fatty voxels is present. Because of the two opposite effects that influence T_2w in a voxel, namely, (i) a pathophysiologically increased water mobility leading to T_2w elevation, and (ii) a dependency of T_2w on the PDFF leading to decreased T_2w , the T_2w distribution within a muscle might be heterogenous and the routine mean or median analysis can lead to a misinterpretation of the muscle health. It was concluded that muscle T_2w mean values can wrongly suggest healthy muscle tissue. A deeper analysis of the underlying value distribution is necessary. Therefore, a quantitative analysis of T_2w histograms is a potential alternative.

Fast Open-Source Toolkit for Water T2 Mapping in the Presence of Fat From Multi-Echo Spin-Echo Acquisitions for Muscle MRI

Imaging has become a valuable tool in the assessment of neuromuscular diseases, and, specifically, quantitative MR imaging provides robust biomarkers for the monitoring of disease progression. Quantitative evaluation of fat infiltration and quantification of the T2 values of the muscular tissue's water component (wT2) are two of the most essential indicators currently used. As each voxel of the image can contain both water and fat, a two-component model for the estimation of wT2 must be used. In this work, we present a fast method for reconstructing wT2 maps obtained from conventional multi-echo spin-echo (MESE) acquisitions and released as Free Open Source Software. The proposed software is capable of fast reconstruction thanks to extended phase graphs (EPG) simulations and dictionary matching implemented on a general-purpose graphic processing unit. The program can also perform more conventional biexponential least-squares fitting of the data and incorporate information from an external water-fat acquisition to increase the accuracy of the results. The method was applied to the scans of four healthy volunteers and five subjects suffering from facioscapulohumeral muscular dystrophy (FSHD). Conventional multi-slice MESE acquisitions were performed with 17 echoes, and additionally, a 6-echo multi-echo gradient-echo (MEGE) sequence was used for an independent fat fraction calculation. The proposed reconstruction software was applied on the full datasets, and additionally to reduced number of echoes, respectively, to 8, 5, and 3, using EPG and biexponential least-squares fitting, with and without incorporating information from the MEGE acquisition. The incorporation of external fat fraction maps increased the robustness of the fitting with a reduced number of echoes per datasets, whereas with unconstrained fitting, the total of 17 echoes was necessary to retain an independence of wT2 from the level of fat infiltration. In conclusion, the proposed software can successfully be used to calculate wT2 maps from conventional MESE acquisition, allowing the usage of an optimized protocol with similar precision and accuracy as a 17-echo acquisition. As it is freely released to the community, it can be used as a reference for more extensive cohort studies.

T2 relaxation-time mapping in healthy and diseased skeletal muscle using extended phase graph algorithms

PURPOSE

Multi-echo spin-echo (MSE) transverse relaxometry mapping using multi-component models is used to study disease activity in neuromuscular disease by assessing the T₂ of the myocytic component (T_2water ). Current extended phase graph algorithms are not optimized for fat fractions above 50% and the effects of inaccuracies in the T_2fat calibration remain unexplored. Hence, we aimed to improve the performance of extended phase graph fitting methods over a large range of fat fractions, by including the slice-selection flip angle profile, a through-plane chemical-shift displacement correction, and optimized calibration of T_2fat .

METHODS

Simulation experiments were used to study the influence of the slice flip-angle profile with chemical-shift and T_2fat estimations. Next, in vivo data from four neuromuscular disease cohorts were studied for different T_2fat calibration methods and T_2water estimations.

RESULTS

Excluding slice flip-angle profiles or chemical-shift displacement resulted in a bias in T_2water up to 10 ms. Furthermore, a wrongly calibrated T_2fat caused a bias of up to 4 ms in T_2water . For the in vivo data, one-component calibration led to a lower T_2fat compared with a two-component method, and T_2water decreased with increasing fat fractions.

CONCLUSION

In vivo data showed a decline in T_2water for increasing fat fractions, which has important implications for clinical studies, especially in multicenter settings. We recommend using an extended phase graph-based model for fitting T_2water from MSE sequences with two-component T_2fat calibration. Moreover, we recommend including the slice flip-angle profile in the model with correction for through-plane chemical-shift displacements.

T2 mapping of the distal sciatic nerve in healthy subjects and patients suffering from lumbar disc herniation with nerve compression

OBJECTIVE

To measure T2 values for magnetic resonance neurography (MRN) of the healthy distal sciatic nerve and compare those to T2 changes in patients with nerve compression.

MATERIALS AND METHODS

Twenty-one healthy subjects and five patients with sciatica due to disc herniation underwent MRN using a T2-prepared turbo spin echo (TSE) sequence of the distal sciatic nerve bilaterally. Six and one of those healthy subjects further underwent a commonly used multi-echo spin-echo (MESE) sequence and magnetic resonance spectroscopy (MRS), respectively.

RESULTS

T2 values derived from the T2-prepared TSE sequence were 44.6 ± 3.0 ms (left) and 44.5 ± 2.6 ms (right) in healthy subjects and showed good inter-reader reliability. In patients, T2 values of 61.5 ± 6.2 ms (affected side) versus 43.3 ± 2.4 ms (unaffected side) were obtained. T2 values of MRS were in good agreement with measurements from the T2-prepared TSE, but not with those of the MESE sequence.

DISCUSSION

A T2-prepared TSE sequence enables precise determination of T2 values of the distal sciatic nerve in agreement with MRS. A MESE sequence tends to overestimate nerve T2 compared to T2 from MRS due to the influence of residual fat on T2 quantification. Our approach may enable to quantitatively assess direct nerve affection related to nerve compression.

Exploration of New Contrasts, Targets, and MR Imaging and Spectroscopy Techniques for Neuromuscular Disease - A Workshop Report of Working Group 3 of the Biomedicine and Molecular Biosciences COST Action BM1304 MYO-MRI

Neuromuscular diseases are characterized by progressive muscle degeneration and muscle weakness resulting in functional disabilities. While each of these diseases is individually rare, they are common as a group, and a large majority lacks effective treatment with fully market approved drugs. Magnetic resonance imaging and spectroscopy techniques (MRI and MRS) are showing increasing promise as an outcome measure in clinical trials for these diseases. In 2013, the European Union funded the COST (co-operation in science and technology) action BM1304 called MYO-MRI (www.myo-mri.eu), with the overall aim to advance novel MRI and MRS techniques for both diagnosis and quantitative monitoring of neuromuscular diseases through sharing of expertise and data, joint development of protocols, opportunities for young researchers and creation of an online atlas of muscle MRI and MRS. In this report, the topics that were discussed in the framework of working group 3, which had the objective to: Explore new contrasts, new targets and new imaging techniques for NMD are described. The report is written by the scientists who attended the meetings and presented their data. An overview is given on the different contrasts that MRI can generate and their application, clinical needs and desired readouts, and emerging methods.

Water T2 Mapping in Fatty Infiltrated Thigh Muscles of Patients With Neuromuscular Diseases Using a T2 -Prepared 3D Turbo Spin Echo With SPAIR

BACKGROUND

Muscle water T₂ (T_2w ) has been proposed as a biomarker to monitor disease activity and therapy effectiveness in patients with neuromuscular diseases (NMD). Multi-echo spin-echo (MESE) is known to be affected by fatty infiltration. A T₂ -prepared 3D turbo spin echo (TSE) is an alternative method for T₂ mapping, but has been only applied in healthy muscles.

PURPOSE

To examine the performance of T₂ -prepared 3D TSE in combination with spectral adiabatic inversion recovery (SPAIR) in measuring T_2w in fatty infiltrated muscles based on simulations and in vivo measurements in thigh muscles of patients with NMD.

STUDY TYPE

Prospective.

SUBJECTS

One healthy volunteer, 34 NMD patients.

FIELD STRENGTH/SEQUENCE

T₂ -prepared stimulated echo acquisition mode (STEAM) magnetic resonance spectroscopy (MRS), SPAIR STEAM MRS, and SPAIR T₂ -prepared STEAM MRS were performed in the subcutaneous fat of a healthy volunteer's thigh. T₂ mapping based on SPAIR 2D MESE and SPAIR T₂ -prepared 3D TSE was performed in the NMD patients' thigh region. Multi-TE STEAM MRS was performed for measuring a reference T_2w at different thigh locations.

ASSESSMENT

The behavior of the fat spectrum in the SPAIR T₂ -prepared 3D TSE was simulated using Bloch simulations. The in vivo T₂ results of the imaging methods were compared to the in vivo T_2w MRS results.

STATISTICAL TESTS

Pearson correlation coefficient with slope and intercept, relative error.

RESULTS

The simulated T₂ for the SPAIR T₂ -prepared 3D TSE sequence remained constant within a relative error of not more than 4% up to a fat fraction of 80%. In vivo T₂ values of SPAIR T₂ -prepared 3D TSE were in good agreement with the T_2w values of STEAM MRS (R = 0.86; slope = 1.12; intercept = -1.41 ms). In vivo T₂ values of SPAIR 2D MESE showed large deviations from the T_2w values of STEAM MRS (R = 0.14; slope = 0.32; intercept = 38.83 ms).

DATA CONCLUSION

The proposed SPAIR T₂ -prepared 3D TSE shows reduced sensitivity to fatty infiltration for T_2w mapping in the thigh muscles of NMD patients.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2020;51:1727-1736.

MR fingerprinting for water T1 and fat fraction quantification in fat infiltrated skeletal muscles

PURPOSE

To develop a fast MR fingerprinting (MRF) sequence for simultaneous estimation of water T1 (T1_H2O ) and fat fraction (FF) in fat infiltrated skeletal muscles.

METHODS

The MRF sequence for T1_H2O and FF quantification (MRF T1-FF) comprises a 1400 radial spokes echo train, following nonselective inversion, with varying echo and repetition time, as well as prescribed flip angle. Undersampled frames were reconstructed at different acquisition time-points by nonuniform Fourier transform, and a bi-component model based on Bloch simulations applied to adjust the signal evolution and extract T1_H2O and FF. The sequence was validated on a multi-vial phantom, in three healthy volunteers and five patients with neuromuscular diseases. We evaluated the agreement between MRF T1-FF parameters and reference values and confounding effects due to B0 and B1 inhomogeneities.

RESULTS

In phantom, T1_H2O and FF were highly correlated with references values measured with multi-inversion time inversion recovery-stimulated echo acquisition mode and Dixon, respectively (R² > 0.99). In vivo, T1_H2O and FF determined by the MRF T1-FF sequence were also correlated with reference values (R² = 0.98 and 0.97, respectively). The precision on T1_H2O was better than 5% for muscles where FF was less than 0.4. Both T1_H2O and FF values were not confounded by B0 nor B1 inhomogeneities.

CONCLUSION

The MRF T1-FF sequence derived T1_H2O and FF values in voxels containing a mixture of water and fat protons. This method can be used to comprehend and characterize the effects of tissue water compartmentation and distribution on muscle T1 values in patients affected by chronic fat infiltration.

Multi-center evaluation of stability and reproducibility of quantitative MRI measures in healthy calf muscles

The purpose of this study was to evaluate temporal stability, multi-center reproducibility and the influence of covariates on a multimodal MR protocol for quantitative muscle imaging and to facilitate its use as a standardized protocol for evaluation of pathology in skeletal muscle. Quantitative T2, quantitative diffusion and four-point Dixon acquisitions of the calf muscles of both legs were repeated within one hour. Sixty-five healthy volunteers (31 females) were included in one of eight 3-T MR systems. Five traveling subjects were examined in six MR scanners. Average values over all slices of water-T2 relaxation time, proton density fat fraction (PDFF) and diffusion metrics were determined for seven muscles. Temporal stability was tested with repeated measured ANOVA and two-way random intraclass correlation coefficient (ICC). Multi-center reproducibility of traveling volunteers was assessed by a two-way mixed ICC. The factors age, body mass index, gender and muscle were tested for covariance. ICCs of temporal stability were between 0.963 and 0.999 for all parameters. Water-T2 relaxation decreased significantly (P < 10^-3 ) within one hour by ~ 1 ms. Multi-center reproducibility showed ICCs within 0.879-0.917 with the lowest ICC for mean diffusivity. Different muscles showed the highest covariance, explaining 20-40% of variance for observed parameters. Standardized acquisition and processing of quantitative muscle MRI data resulted in high comparability among centers. The imaging protocol exhibited high temporal stability over one hour except for water T2 relaxation times. These results show that data pooling is feasible and enables assembling data from patients with neuromuscular diseases, paving the way towards larger studies of rare muscle disorders.

T2 mapping of lumbosacral nerves in patients suffering from unilateral radicular pain due to degenerative disc disease

OBJECTIVE

Lumbosacral radicular syndrome (LRS) is a very common condition, often requiring diagnostic imaging with the aim of elucidating a structural cause when symptoms are longer lasting. However, findings on conventional anatomical MRI do not necessarily correlate with clinical symptoms, and it is primarily performed for the qualitative evaluation of surrounding compressive structures, such as herniated discs, instead of to evaluate the nerves directly. The present study investigated the performance of quantitative imaging by using magnetic resonance neurography (MRN) in patients with LRS.

METHODS

Eighteen patients (55.6% males, mean age 64.4 ± 10.2 years), with strict unilateral LRS matching at least one dermatome and suspected disc herniation, underwent high-resolution 3-T MRN using T2 mapping. On T2 maps, the presumably affected and contralateral unaffected nerves were identified; subsequent regions of interest (ROIs) were placed at preganglionic, ganglionic, and postganglionic sites; and T2 values were extracted. Patients then underwent an epidural steroid injection (ESI) with local anesthetic agents at the site of suspected nerve affection. T2 values of the affected nerves were compared against the contralateral nerves. Furthermore, receiver operating characteristics were calculated based on the measured T2 values and the responsiveness to ESI.

RESULTS

The mean T2 value was 77.3 ± 1.9 msec for affected nerves and 74.8 ± 1.4 msec for contralateral nerves (p < 0.0001). In relation to ESI performed at the site of suspected nerve affection, MRN with T2 mapping had a sensitivity/specificity of 76.9%/60.0% and a positive/negative predictive value of 83.3%/50.0%. Signal alterations in affected nerves according to qualitative visual inspection were present in only 22.2% of patients.

CONCLUSIONS

As one of the first of its kind, this study revealed elevated T2 values in patients suffering from LRS. T2 values of lumbosacral nerves might be used as more objective parameters to directly detect nerve affection in such patients.

Quantitative MRI Musculoskeletal Techniques: An Update

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

T2-Weighted Dixon Turbo Spin Echo for Accelerated Simultaneous Grading of Whole-Body Skeletal Muscle Fat Infiltration and Edema in Patients With Neuromuscular Diseases

OBJECTIVE

The assessment of fatty infiltration and edema in the musculature of patients with neuromuscular diseases (NMDs) typically requires the separate performance of T1-weighted and fat-suppressed T2-weighted sequences. T2-weighted Dixon turbo spin echo (TSE) enables the generation of T2-weighted fat- and water-separated images, which can be used to assess both pathologies simultaneously. The present study examines the diagnostic performance of T2-weighted Dixon TSE compared with the standard sequences in 10 patients with NMDs and 10 healthy subjects.

METHODS

Whole-body magnetic resonance imaging was performed including T1-weighted Dixon fast field echo, T2-weighted short-tau inversion recovery, and T2-weighted Dixon TSE. Fatty infiltration and intramuscular edema were rated by 2 radiologists using visual semiquantitative rating scales. To assess intermethod and interrater agreement, weighted Cohen's κ coefficients were calculated.

RESULTS

The ratings of fatty infiltration showed high intermethod and high interrater agreement (T1-weighted Dixon fast field echo vs T2-weighted Dixon TSE fat image). The evaluation of edematous changes showed high intermethod and good interrater agreement (T2-weighted short-tau inversion recovery vs T2-weighted Dixon TSE water image).

CONCLUSIONS

T2-weighted Dixon TSE imaging is an alternative for accelerated simultaneous grading of whole-body skeletal muscle fat infiltration and edema in patients with NMDs.

High Isotropic Resolution T2 Mapping of the Lumbosacral Plexus with T2-Prepared 3D Turbo Spin Echo

PURPOSE

Isotropic high-resolution three-dimensional (3D) magnetic resonance neurography (MRN) is increasingly used to depict even small and highly oblique nerves of the lumbosacral plexus (LSP). The present study introduces a T2 mapping sequence (T2-prepared 3D turbo spin echo) that is B1-insensitive and enables quantitative assessment of LSP nerves.

METHODS

In this study 15 healthy subjects (mean age 28.5 ± 3.8 years) underwent 3 T MRN of the LSP area three times. The T2 values were calculated offline on a voxel-by-voxel basis and measured at three segments (preganglionic, ganglionic, postganglionic) of three LSP nerves (S1, L5, L4) by two independent investigators (experienced and novice). Normative data for the different nerves were extracted and intraclass correlation coefficients (ICCs) were calculated to assess reproducibility and interobserver reliability of T2 measurements.

RESULTS

The T2 mapping showed excellent reproducibility with ICCs ranging between 0.99 (S1 preganglionic) and 0.89 (L5 postganglionic). Interobserver reliability was less robust with ICCs ranging between 0.78 (S1 preganglionic) and 0.44 (L5 postganglionic) for S1 and L5. A mean T2 value of 74.6 ± 4.7 ms was registered for preganglionic segments, 84.7 ± 4.1 ms for ganglionic and 65.4 ± 2.5 ms for postganglionic segments, respectively. There was a statistically significant variation of T2 values across the nerve (preganglionic vs ganglionic vs postganglionic) for S1, L5, and L4.

CONCLUSION

Our approach enables isotropic high-resolution and B1-insensitive T2 mapping of LSP nerves with excellent reproducibility. It might reflect a robust and clinically useful method for future diagnostics of LSP pathologies.

T2 mapping with magnetization-prepared 3D TSE based on a modified BIR-4 T2 preparation

The purpose of this work was to investigate the performance of the modified BIR-4 T₂ preparation for T₂ mapping and propose a method to remove T₂ quantification errors in the presence of large B₁ and B₀ offsets. The theoretical investigation of the magnetization evolution during the T₂ preparation in the presence of B₁ and B₀ offsets showed deviations from a mono-exponential T₂ decay (two parameter fit). A three parameter fit was used to improve T₂ accuracy. Furthermore, a two parameter fit with an additional saturation preparation scan was proposed to improve T₂ accuracy and precision. These three fitting methods were compared based on simulations, phantom measurements and an in vivo healthy volunteer study of the neck musculature using a 3D TSE readout. The results based upon the pure two parameter fit overestimated T₂ in regions with high B₀ offsets (up to 40% in phantoms). The three parameter fit T₂ values were robust to B₀ offsets but with higher standard deviation (up to 40% in simulations). The two parameter fit with the saturation preparation yielded high robustness towards B₀ offsets with a noise performance comparable to that of the two parameter fit. In the volunteer study the T₂ values obtained by the pure two parameter fit showed a dependence on the field inhomogeneities, whereas the T₂ values from the proposed fitting approach were shown to be insensitive to B₀ offsets. The proposed method enabled accurate and precise T₂ mapping in the presence of large B₁ and B₀ offsets.