4145

4D-MRI of Thorax and Abdomen with Variable Contrast at 1.5T
Orso Andrea Pusterla1,2,3, Francesco Santini2,3, Grzgorz Bauman2,3, Rahel Heule4, Alina Giger3, Philippe Claude Cattin3, Sairos Safai5, Sebastian Kozerke1, and Oliver Bieri2,3
1Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland, 2Department of Radiology, Division of Radiological Physics, University Hospital Basel, Basel, Switzerland, 3Department of Biomedical Engineering, University of Basel, Allschwil, Switzerland, 4High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany, 5Center for Proton Therapy, Paul Scherrer Institute, Villigen-PSI, Switzerland

Synopsis

Time-resolved volumetric imaging (4D-MRI) of moving organs is essential for several clinical applications, e.g., in interventional treatments, radiation therapies, and high-intensity focused ultrasound ablations. In this work, an interleaved time-resolved data-navigator acquisition strategy for 4D-MRI with flexible contrast is developed offering versatile settings for thoracic and abdominal imaging: for both, image and navigator, not only the flip angle but also the basic acquisition type (pulse sequence kernel) can be chosen independently from either balanced SSFP (i.e., T2/T1-weighted) or spoiled gradient echo (i.e., T1-weighted).

Introduction

Respiratory and organ motion challenge treatments of thoracic and abdominal targets (such as tumors), e.g., during radiotherapeutic interventions or high-intensity focused ultrasound ablations1,2. Hence, these clinical applications strongly rely on real-time and time-resolved volumetric MR imaging, occasionally paired with other imaging modalities such as ultrasonography3 or optical tracking devices4. For treatment planning, accurate organ motion modeling is necessary, while for therapy delivery, highly precise target localization and tracking is crucial.
One of the most sophisticated methods for capturing motions and drifts of targeted organs is based on continuous free-breathing acquisitions of 2D data-navigator slice pairs (with moving imaging data-slice and stationary navigator-slice), retrospectively stacked into time-resolved three-dimensional datasets (4D-MRI)5. While data and navigator images generally should have diametrical objectives, i.e., providing pathology-tuned image contrasts (data) and feature tracking for accurate motion modeling (navigator), current implementations do not offer sufficient flexibility since the imaging is performed with identical acquisition types and settings (i.e. pulse sequence kernels and flip angles). Moreover, the magnetization is in quasi steady-state for the navigator (continuously acquired), but in thermal equilibrium for the acquisition of the consecutive data slices, leading to different optimal acquisition parameter settings.
In this work, we thus propose Variable Contrast (VC) data-navigator imaging for 4D-MRI (VC-4DMRI). The implementation offers maximum flexibility, that is, for both, data and navigator images, not only different flip angles (FA) but also the sequence kernel can be set independently to offer either balanced steady-state free precession (bSSFP) or spoiled gradient echo (SPGR) contrast. Generally, VC-4DMRI is developed based on an ultra-fast SSFP kernel6, which offers substantially reduced repetition times and thus artifact-free thoracic bSSFP imaging.

Methods

VC-4DMRI was tested in healthy volunteers (in accordance with local ethical guidelines) on a 1.5T whole-body MR-system (MAGNETOM Avanto-Fit, Siemens Healthineers, Erlangen, Germany). A schematic illustration of the proposed VC-4DMRI scheme is shown in Figure 1, allowing selection of a specific contrast, that is T2/T1 (bSSFP) or T1 (SPGR), and a specific flip angle independently for both, data and navigator acquisitions. In addition, an interleaved slice-order acquisition scheme was employed to reduce partial saturation due to the overlapping slices (Figure 2).
Imaging was performed in the coronal orientation with a TE/TR = 0.78/1.88 ms, bandwidth = 1860 Hz/pixel, FA = 5-90°, GRAPPA factor 2, partial Fourier factor 7/8, field-of-view = 425×425 mm2, in-plane resolution = 2.2×2.2 mm2 (matrix size = 192×192), slice thickness = 9 mm, slice overlap = 35% (see Figure 2), acquisition time per image = 205 ms (i.e., sampling frequency = 5Hz).
Image stacking into 4D datasets was performed according to von Siebenthal et al.5, but by using a thorax-specific image registration approach7 for determining the respiratory state.

Results

Representative data-navigator acquisitions with VC-4DMRI using bSSFP-T2/T1 and SPGR-T1 contrasts are shown in Figure 3. Both pulse sequences (bSSFP and SPGR kernels) offer artifact-free thorax images. Vessels and fine vascular structures are well discernible in bSSFP, which clearly delivers a higher signal-to-noise-ratio (SNR) in comparison to SPGR (SNR bSSFP ≈ 35, SNR SPGR ≈ 11). In the spine and next to the ribcage, bSSFP suffers from signal drops. The interleaved slice-order acquisition scheme (see Figure 2) provides a signal increase of about 8% in the lung parenchyma and in the liver as compared to the sequential slice-order acquisition scheme.
The versatility of adapting the flip-angle for data and navigator imaging can be appreciated in Figure 4. VC-4DMRI offers for instance an improved discrimination of the hepatic vessels from the surrounding liver tissue using bSSFP with a very large FA (here 90°). Such a high FA for data-slices is only achievable combined with low-FA bSSFP or SPGR navigator-slices (here bSSFP with 10°) due to specific absorption rate (SAR) constraints (SAR$$$\propto$$$FA2).
Finally, a representative volume of a reconstructed VC-4DMRI dataset is shown in Figure 5. With the proposed setup, i.e. low-FA navigator combined with high-FA data-slices and both set to bSSFP contrast, time-resolved volumetric reconstructions are exceptionally improved since no saturation of signal is discernible around the navigator plane, as commonly observed with conventional approaches (i.e. using the same FA for data and navigator3).

Discussion and Conclusion

In this work, we have broadened and improved contemporary 4D-MRI with the possibility to acquire data and navigator slices with variable contrasts, allowing dedicated, versatile and optimized protocol settings for motion detection and target imaging with optimal contrast and improved SNR.
VC-4DMRI provides time-resolved volumetric imaging, e.g., of lung and liver with flexible contrast, high frame rate and high resolution (here: 2.5Hz for data slices, 2.2×2.2×5 mm3 after stacking reconstruction). In general, bSSFP acquisitions appear superior to SPGR, nevertheless the SPGR contrast might be advantageous for an improved differentiation and characterization of e.g. tumor tissues, or for imaging at field strength higher than 1.5T.
In summary, VC-4DMRI offers excellent prospects for time-resolved volumetric imaging of thoracic and abdominal organs during free-breathing with wide contrast flexibility for treatment planning and interventions.

Acknowledgements

This work was supported by the PHRT Initiative of the ETH Domain, Switzerland.

References

1. Bert C and Durante M. Motion in radiotherapy: particle therapy. Phys. Med. Biol. 56, 2011, R113.

2. Celicanin Z, Auboiroux V, Bieri O, et al. Real-time method for motion-compensated MR thermometry and MRgHIFU treatment in abdominal organs: MRgHIFU treatment method in abdominal organs. Magn Reson Med 2014;72:1087–95.

3. Giger A, Stadelmann M, Preiswerk F, et al. Ultrasound-driven 4D MRI. Phys Med Biol. 2018 Jul 16;63(14):145015.

4. Fattori G, Safai S, Carmona PF, et al. Monitoring of breathing motion in image-guided PBS proton therapy: comparative analysis of optical and electromagnetic technologies. Radiat Oncol 2017; 12:63. 10.1186/s13014-017-0797-9.

5. Von Siebenthal M, Székely G, Gamper U, et al. A 4D MR imaging of respiratory organ motion and its variability. Phys. Med. Biol., 2007, 52 1547.

6. Bieri O. Ultra-fast steady state free precession and its application to in vivo 1H morphological and functional lung imaging at 1.5 tesla. Magn Reson Med 2013; 70:657–663.

7. Sandkühler R, Pezold CS, Cattin P.C, et al. Adaptive graph diffusion regularisation for discontinuity preserving image registration. Int. Workshop on Biomedical Image Registration. Springer, Cham, 2018.

Figures

Figure 1. Schematic of the VC-4DMRI sequence for data-navigator acquisitions: the flip angles (marked in green) and the contrast (bSSFP or SPGR) can be selected independently for data and navigator. SPGR contrast is achieved by applying crusher gradients (colored in blue) and RF spoiling (not shown).

Figure 2. Exemplary order of data-navigator coronal acquisitions: slice-sequential (left), and slice-interleaved (right). Note the partial overlap between slices (both left and right). In the sequential scheme, data (D) and navigator (N) acquisitions are performed according to {D1,N,D2,N,D3,N,…} until the whole lung is sampled from anterior to posterior. In the interleaved scheme, the acquisition order is set up as {D1,N,D3,N,…,D2,N,D4,N,…} to avoid partial saturation caused by slice cross talk. On the right, data D2 and D4 are only partially shown for enhanced visualization.

Figure 3. VC-4DMRI data-navigator acquisitions using the four contrast possibilities offered by the sequence: bSSFP-bSSFP, bSSFP-SPGR, SPGR-SPGR, and SPGR-bSSFP. Here, the flip angles for bSSFP were set to 50°, and for SPGR to 5°.

Figure 4. VC-4DMRI acquired with bSSFP for both data slice and navigator, but using different FAs. Due to the mixed T2/T1 contrast of bSSFP and its non-linear behavior as a function of the FA, imaging using different FA allows for better localizing and identifying diverse anatomical structures: e.g., the hepatic vessels (marked by red arrowheads) using a FA of 90°, which are otherwise not well discernible from surrounding tissue using a FA of 50°. Due to specific absorption rate (SAR) constraints, a very large FA in the data slice can only be achieved in combination with a low-FA navigator.

Figure 5. Time-resolved volumetric data reconstructed from a 5min 30s long VC-4DMRI scan (28 data-slices measured 30 times): data and navigator were set to bSSFP contrast (FA data 90°, FA navigator 10°). A maximum intensity projection filter was applied for enhanced visualization. When the same FAs for data and navigator are used, as in classical approaches, 4D reconstructions are characterized by a signal drop (saturation) around the navigator plane. Here (FAs 90° and 10°), this signal drop is not discernible since the low-FA navigator saturates only marginally the data-slices.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
4145