2068

Respiratory motion-registered isotropic 3D whole-heart T2 mapping in patients with suspected inflammatory myocardial injury
Karolina Dorniak1, Lorenzo Di Sopra2, Agniezka Sabisz3, Anna Glinska3, Christopher W Roy2, Kamil Gorczewski4, Davide Piccini2,5, Jérôme Yerly2,6, Jadwiga Fijałkowska3, Edyta Szurowska3, Matthias Stuber2,6, and Ruud B van Heeswijk2
1Department of Noninvasive Cardiac Diagnostics, Medical University of Gdansk, Gdansk, Poland, 2Radiology, Lausanne University Hospital (CHUV) and University of Lausanne (UNIL), Lausanne, Switzerland, 32nd Department of Radiology, Medical University of Gdansk, Gdansk, Poland, 4Siemens Healthineers, Erlangen, Germany, 5Advanced Clinical Imaging Technology, Siemens Healthcare AG, Lausanne, Switzerland, 6Center for Biomedical Imaging (CIBM), Lausanne, Switzerland

Synopsis

T2 mapping can be used to effectively detect myocardial edema and inflammation. However, the focal nature of myocardial inflammation may render 2D approaches suboptimal and make whole-heart isotropic 3D mapping desirable. Unfortunately, at 1.5T, self-navigated 3D radial bSSFP results in too noisy images for adequate T2 mapping. In this study, we therefore used a respiratory motion-resolved reconstruction together with image registration to improve the 3D T2 mapping precision and accuracy at 1.5T in patients with inflammatory myocardial injury. The resulting myocardial T2 values matched those of the routine 2D T2 maps, with no discernible bias and slightly lower precision.

Background

Myocardial inflammation can be detected as elevated T2 relaxation times on T2 maps [1]. However, the focal nature of acute non-ischemic inflammatory myocardial injury in a spectrum of inflammatory conditions may render 2D T2 mapping of the heart suboptimal due to either insufficient spatial coverage or prolonged scan times. Hence, a non-invasive whole-heart isotropic 3D T2 quantification [2] is desirable in this setting, as was demonstrated previously in patients with acute graft rejection at 3T [3]. 3D techniques with radial [2], Cartesian [4,5], or hybrid [6] trajectories can be chosen, with intrinsic advantages such as inherent motion robustness for radial imaging and high signal-to-noise for Cartesian imaging. Unfortunately, at 1.5T, self-navigated 3D radial bSSFP results in too noisy images for adequate myocardial T2 mapping. To address this shortcoming, in this study we aimed to use respiratory motion-resolved reconstruction together with image registration of the resulting 4D images to decrease noise and motion artifacts, and thus to improve 3D T2 mapping precision and accuracy at 1.5T. These respiratory motion-resolved T2 maps were then compared to self-navigated 3D maps that were reconstructed from the same data, as well as to routine 2D maps.

Materials and Methods

This study was approved by the Institutional Review Board. Consecutive patients with suspected acute non-ischemic inflammatory myocardial injury (n=12, 6 cardiac sarcoidosis, 2 myocarditis, 2 acute graft rejection, 2 systemic sclerosis; 5 (42%) female, age 51±11y) were recruited. On a 1.5T clinical scanner (MAGNETOM Aera, Siemens Healthcare, Erlangen, Germany), routine 2D T2-prepared bSSFP T2 maps (pixel size 1.9x1.9mm2, slice thickness 8mm) [7] were acquired in a short-axis (SAX) orientation. Next, three 3D radial bSSFP volumes were acquired during free breathing (isotropic voxel size 1.6mm3, T2prep=0,30,60ms, flip angle=35°, interleaves of 50 lines acquired every third heartbeat) [2,8] using a prototype sequence.
For the respiratory motion-resolved reconstruction, a principal component analysis (PCA) was performed on the superior-inferior profiles that were acquired at the start of each interleave [9,10], and to partition the dataset into 4 different respiratory states. 4D (x-y-z-respiratory dimensions) images were then reconstructed with a compressed sensing algorithm that exploits sparsity along the respiratory dimension [11]. All respiratory bins were translationally and then non-rigidly registered to the end-expiration bin with Elastix [12], and were then averaged. After a second, similar, registration of the resulting three averaged T2-prepared images, pixel-wise T2 mapping [2] was performed. Since the motion is no longer resolved after these registrations, we named the resulting 3D map “motion-registered” T2 maps.
For the self-navigated reconstruction from the same 3D radial data, the 1D displacement of the left-ventricular blood pool along the superior-inferior readouts acquired at the start of each interleave was used to correct each interleave for respiratory motion in k-space prior to image reconstruction [8]. The resulting three 3D images were translationally and then non-rigidly registered, and pixel-wise T2 mapping [2] was performed.
The T2 values of the entire visible myocardium in the routine 2D maps and a matching single slice in both 3D volumes were then measured and compared using Matlab.

Results and Discussion

The respiratory motion-resolved reconstruction resulted in visibly well-separated motion states (Figure 1), and motion-registered isotropic 3D T2 maps of the heart were successfully obtained in all patients (Figure 2). Their myocardial T2 values matched those of the routine maps (51.1±3.5ms vs. 48.6±4.1ms, respectively, P=0.10), while the Bland-Altman analysis demonstrated that there was no discernible trend or bias when compared to the routine maps (Figure 3).
The motion-registered 3D T2 mapping precision was slightly lower than that of the routine technique, as evidenced by the myocardial T2 standard deviation (5.7±1.3ms vs. 4.5±1.9ms, P=0.03). This is most likely a combination of several factors: on the one hand, the compressed sensing reconstruction removes some noise, while on the other hand, residual 3D radial undersampling artifacts and local misregistration artifacts, as well as the much smaller voxel size will increase the T2 standard deviation. As expected, the self-navigated version resulted in significantly higher myocardial T2 values (58.2±9.2ms, P<0.001 versus both other techniques) and T2 standard deviation (9.7±5.3ms, P<0.001 versus both other techniques) compared to the routine maps, confirming the need for the motion-resolved reconstruction.

Conclusions

Respiratory motion-resolved 3D radial imaging at 1.5T led to precise and accurate isotropic 3D whole-heart T2 maps in a small patient cohort. Future efforts will be directed towards studies in larger patient cohorts as well as accelerating the mapping.

Acknowledgements

No acknowledgement found.

References

1. Messroghli DR, Moon JC, Ferreira VM, Grosse-Wortmann L, He T, Kellman P, Mascherbauer J, Nezafat R, Salerno M, Schelbert EB, Taylor AJ, Thompson R, Ugander M, van Heeswijk RB, Friedrich MG. Clinical recommendations for cardiovascular magnetic resonance mapping of T1, T2, T2* and extracellular volume: A consensus statement by the Society for Cardiovascular Magnetic Resonance (SCMR) endorsed by the European Association for Cardiovascular Imaging (EACVI). J Cardiovasc Magn Reson. 2018 Feb 7;20(1):9. doi: 10.1186/s12968-017-0408-9

2. van Heeswijk RB, Piccini D, Feliciano H, Hullin R, Schwitter J, Stuber M. Self-navigated isotropic three-dimensional cardiac T2 mapping. Magn Reson Med. 2015 Apr;73(4):1549-54. doi: 10.1002/mrm.25258

3. van Heeswijk RB, Piccini D, Tozzi P, Rotman S, Meyer P, Schwitter J, Stuber M, Hullin R. Three-Dimensional Self-Navigated T2 Mapping for the Detection of Acute Cellular Rejection After Orthotopic Heart Transplantation. Transplant Direct. 2017 Mar 28;3(4):e149. doi: 10.1097/TXD.0000000000000635

4. Ding H, Fernandez-de-Manuel L, Schär M, Schuleri KH, Halperin H, He L, Zviman MM, Beinart R, Herzka DA. Three-dimensional whole-heart T2 mapping at 3T. Magn Reson Med. 2015 Sep;74(3):803-16. doi: 10.1002/mrm.25458

5. Milotta G, Ginami G, Bustin A, Neji R, Prieto C, Botnar RM. 3D Whole-heart free-breathing qBOOST-T2 mapping. Magn Reson Med. 2019 Oct 21. doi: 10.1002/mrm.28039. [Epub ahead of print] 6. Yang HJ, Sharif B, Pang J, Kali A, Bi X, Cokic I, Li D, Dharmakumar R. Free-breathing, motion-corrected, highly efficient whole heart T2 mapping at 3T with hybrid radial-cartesian trajectory. Magn Reson Med. 2016 Jan;75(1):126-36. doi: 10.1002/mrm.25576. Epub 2015 Mar 6.

7. Giri S, Chung YC, Merchant A, Mihai G, Rajagopalan S, Raman SV, Simonetti OP. T2 quantification for improved detection of myocardial edema. J Cardiovasc Magn Reson. 2009 Dec 30;11:56. doi: 10.1186/1532-429X-11-56

8. Piccini D, Littmann A, Nielles-Vallespin S, Zenge MO. Epub 2011 Dec 28. Respiratory self-navigation for whole-heart bright-blood coronary MRI: methods for robust isolation and automatic segmentation of the blood pool. Magn Reson Med. 2012 Aug;68(2):571-9. doi: 10.1002/mrm.23247.

9. Di Sopra L, Piccini D, Coppo S, Stuber M, Yerly J. An automated approach to fully self-gated free-running cardiac and respiratory motion-resolved 5D whole-heart MRI. Magn Reson Med. 2019 Dec;82(6):2118-2132. doi: 10.1002/mrm.27898

10. Feng L, Axel L, Chandarana H, Block KT, Sodickson DK, Otazo R. XD-GRASP: Golden-angle radial MRI with reconstruction of extra motion-state dimensions using compressed sensing. Magn Reson Med. 2016 Feb;75(2):775-88. doi: 10.1002/mrm.25665

11. Piccini D, Feng L, Bonanno G, Coppo S, Yerly J, Lim RP, Schwitter J, Sodickson DK, Otazo R, Stuber M. Four-dimensional respiratory motion-resolved whole heart coronary MR angiography. Magn Reson Med. 2017 Apr;77(4):1473-1484. doi: 10.1002/mrm.26221

12. Klein S, Staring M, Murphy K, Viergever MA, Pluim JP. Elastix: a toolbox for intensity based medical image registration. IEEE Trans Med Imaging. 2010 Jan;29(1):196-205. doi: 10.1109/TMI.2009.2035616

Figures

Figure 1. An animation of the respiratory-resolved reconstruction of a source image with T2prep duration of 30ms in a patient. The breathing motion is well-resolved, and results in a non-rigid translation of the heart.

Figure 2. Motion-registered 3D T2 maps of the heart of three patients. A-C) Routine 2D SAX T2 maps. D-F) Matching slices from the motion-registered 3D T2 maps; T2 values closely match those of the routine maps. G-L) Perpendicular long-axis (LAX) maps from the same 3D datasets as in D-F. Patient 1 (female, 37y.o.) and patient 2 (male, 44y.o.) had myocarditis, while patient 3 (male, 65y.o.) had sarcoidosis. Local minor T2 fluctuation can be observed in all three orientations of the 3D T2 maps and all three patients. However, global T2 values did not significantly differ in any of the patients.

Figure 3. Bland-Altman analysis of the 4D T2 mapping technique compared to the routine mapping technique. The bias is slightly but not significantly positive, while there is no discernible trend or outlying value.

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