3307

CUTE: Compressed Sensing UTE for Multi-Echo T2* Mapping
Stefan Sommer1,2, Tom Hilbert3,4,5, Constantin von Deuster1,2, Natalie Hinterholzer2, Markus Klarhöfer1, and Daniel Nanz2,6
1Siemens Healthcare, Zurich, Switzerland, 2Swiss Center for Musculoskeletal Imaging (SCMI), Balgrist Campus, Zurich, Switzerland, 3Advanced Clinical Imaging Technology (ACIT), Siemens Healthcare, Lausanne, Switzerland, 4Department of Radiology, Lausanne University Hospital (CHUV), Lausanne, Switzerland, 5LTS5, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland, 6University of Zurich, Zurich, Switzerland

Synopsis

Clinical application of UTE (ultrashort echo time) sequences in musculoskeletal studies remains limited, despite their potential to unveil quantitative short-T2* information of e.g. tendons or ligaments. The major shortcoming of quantitative mapping is the considerably longer scan time compared to conventional, anatomical imaging sequences. We investigated the feasibility and use of an isotropic 3D radial compressed sensing UTE (CUTE) prototype sequence for quantitative mapping of short-T2* components in the knee. Similar image quality and T2* values were obtained for non-accelerated and accelerated acquisitions. Thus, CUTE sequences show great promise for fast and clinically feasible UTE T2* mapping.

Introduction

Transverse water-proton magnetization in densely packed tissues, such as ligament, tendon, and cortical bone, exhibits fast T2* decay due to slow molecular re-orientation and long correlation times. Sequences with center-out sampling trajectories and ultra-short echo times (UTE) can still detect a signal from rapidly decaying T2* tissue. Acquisition of multiple echoes allows for quantification of short and ultra-short relaxation times. Thus, UTE-type sequences potentially allow for exploration and investigation of early changes in tissue structure or composition, which might be indicative of pathologies such as tendon or meniscal degeneration1 or osteoarthritis2. UTE sequences have also been used to monitor treatment progression after anterior cruciate ligament repair or replacement3.
However, the validation of T2* mapping sequences that are also able to quantify short-T2* times remains limited, possibly due to the increased scan time in comparison to conventional imaging sequences and the need for gradient calibration4 if non-cartesian k-space sampling trajectories are used.
The aim of this study was to explore the potential of an isotropic 3D radial UTE prototype sequence in combination with a compressed sensing reconstruction for the quantification of short-T2* values of the knee within a clinically acceptable scan time.

Materials and Methods

Experiments were performed on a 3T whole-body MR scanner (MAGNETOM Prisma, Siemens Healthcare, Erlangen, Germany) equipped with a 1Tx/15Rx-channel knee coil (QED, Quality Electrodynamics, Mayfield Village, OH, USA).
Isotropic 3D UTE imaging was performed with a prototype sequence that acquired data in 5 consecutive radial center-out spokes (echoes), each followed by a rewinder gradient. Each spoke was acquired with a unique, slightly rotated orientation in k-space. A schematic of the sequence diagram and readout trajectory is shown in Figure 1. Multi-echo UTE data was acquired in the knee of a healthy volunteer with 3 different sets of 24’000 (R=1), 12’000 (R=2) and 6’000 (R=4) spokes and with corresponding acquisition times of 5:24, 2:42 and 1:21 [min:sec], respectively. The remaining parameters were identical for the three acquisitions: TR/TE 13.5/[0.04, 2.46, 4.91, 7.38, 9.84] ms, flip angle: 5°, matrix: 256, FoV: 160x160x160mm3, isotropic voxel size: 0.6mm, bandwidth: 1220Hz/px, RF-pulse duration: 40µs, dead time between RF and ADC: 10µs.
Prior to data gridding, the k-space trajectory was corrected for gradient imperfections using a gradient impulse response model4. The accelerated datasets (R=2, R=4) were reconstructed by an iterative compressed sensing approach5,6. The algorithm splits the problem to alternately solve data-consistency and enforce sparsity by soft-thresholding after a 3D wavelet transformation7. Regularization parameters were tuned manually.
Finally, R2* and M0 maps were estimated by a mono-exponential, nonlinear least‐squares fit for all datasets. Mean and standard deviation of T2* were calculated in ROIs (region of interest) placed in the tibial and femoral bone marrow, vastus lateralis and medialis muscle, and patellar tendon, respectively.

Results

Figure 2 shows sagittal magnitude and phase images for all 5 echoes from the fully sampled (R=1), twice (R=2) and four-fold (R=4) under-sampled data sets. While the twice under-sampled datasets were visually comparable to the full reference scans a lower SNR was observed for R=4, especially in the later echoes.
Figure 3 depicts the fitting results (M0 and R2*) in three orthogonal image planes for the three acquisitions (R=1, 2 and 4). All three data sets allowed for successful delineation of fine structures such as the medial collateral ligament. A prominent difference was an underestimation of cortical-bone R2* in both under-sampled datasets compared to the fully sampled dataset.
Figure 4 summarizes mean T2* value and standard deviation in different ROIs for R=1, 2 and 4. Relaxation times were in good agreement with literature and only minor discrepancies between different acquisitions were detected.

Discussion

Compressed sensing allows for under-sampling by a factor of 2-4 without detrimental loss in image quality. However, in areas with low signal intensities and very short-T2* (e.g. cortical bone), the SNR-loss introduced a bias into parts of the quantitative maps. A limitation of this study was the mono-exponential T2* decay model. Gained scan efficiency could be used to measure more echoes, which might allow for more sophisticated relaxation models such as multi-compartment fits.
It is also important to keep in mind that anisotropic structures like tendons and ligaments exhibit orientation-dependent T2* values due to the magic-angle effect, as presumably apparent in Fig 2., where the R2* values of the quadriceps tendon close to the patella are decreased in comparison to the rest of the tendon.

Conclusion

CUTE has the potential to bring a quantitative assessment of short-T2* species into the realm of clinical feasibility. It allows an efficient generation of UTE images, proton-density and T2* maps, which could enable a grading and assessment of early degenerative and other pathological processes in e.g. tendons, ligaments or menisci.

Acknowledgements

No acknowledgement found.

References

1. Breda, S. J. et al. Tissue-Specific T2* Biomarkers in Patellar Tendinopathy by Subregional Quantification Using 3D Ultrashort Echo Time MRI. J. Magn. Reson. Imaging 52, 420–430 (2020).

2. Titchenal, M. R. et al. Cartilage Subsurface Changes to Magnetic Resonance Imaging UTE-T2* 2 Years After Anterior Cruciate Ligament Reconstruction Correlate With Walking Mechanics Associated With Knee Osteoarthritis. Am. J. Sports Med. 46, 565–572 (2018).

3. Chu, C. R. & Williams, A. A. Quantitative MRI UTE-T2* and T2* Show Progressive and Continued Graft Maturation Over 2 Years in Human Patients After Anterior Cruciate Ligament Reconstruction. Orthop. J. Sport. Med. 7, (2019).

4. Vannesjo, S. J. et al. Image reconstruction using a gradient impulse response model for trajectory prediction. Magn. Reson. Med. 76, 45–58 (2016).

5. Lustig, M., Donoho, D. & Pauly, J. M. Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn. Reson. Med. 58, 1182–1195 (2007).

6. Donoho, D. L. Compressed sensing. IEEE Trans. Inf. Theory 52, 1289–1306 (2006).

7. Huang, F. et al. A rapid and robust numerical algorithm for sensitivity encoding with sparsity constraints: Self-feeding sparse SENSE. Magn. Reson. Med. 64, 1078–1088 (2010).

Figures

A schematic of the pulse diagram is depicted in (a). The TE is defined as the duration between the center of the radio frequency (RF) pulse and the start of the readout gradient. The dead time is defined as the time interval between the end of the RF-pulse and the beginning of the analog-to-digital converter (ADC). In (b) a simplified 2D representation of the readout directions in k-space is shown. Each spoke was rotated by a small angle for the joint sampling pattern to approach a uniform spherical sampling, while the angular difference between consecutive spokes was kept small.

Magnitude images and phase maps of a sagittal slice. Fully sampled (R=1, top), two-fold (R=2), and four-fold (R=4) under-sampled (bottom) datasets are shown for each of the 5 echoes (left to right), prior to R2* fitting.

Fitting results showing the initial signal strength (proportional to the transverse magnetization, M0) and R2* maps for the fully sampled (R=1), twice under-sampled (R=2) and four-fold under-sampled (R=4) acquisitions in three orthogonal planes. Identical scaling was used for all M0 [normalized signal intensity] and R2* [ms-1] maps.

T2* mean values and standard deviation for selected regions of interest (ROIs) and acceleration factors (R=1, R=2, R=4) are shown.

Proc. Intl. Soc. Mag. Reson. Med. 29 (2021)
3307