1095

Vessel architectural imaging in the human heart using heartbeat-to-heartbeat GESE-EPI
Maaike van den Boomen1,2, Mary Kate Manhard1,3, Kyrre E. Emblem4, David E. Sosnovik1,5,6, Niek H.J. Prakken2, Christopher Nguyen1,5,6, Kawin Setsompop1,3,7, and Ronald J.H. Borra2,8
1A.A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, United States, 2Department of Radiology, University Medical Center Groningen, Groningen, Netherlands, 3Department of Radiology, Harvard Medical School, Boston, MA, United States, 4Department of Diagnostic Physics, Oslo University Hospital, Oslo, Norway, 5Cardiovascular Research Center, Massachusetts General Hospital, Boston, MA, United States, 6Department of Medicine, Harvard Medical School, Boston, MA, United States, 7Division of Health Sciences and Technology, Harvard-MIT, Cambridge, MA, United States, 8Department of Nuclear Medicine and Molecular Imaging, University Medical Center Groningen, Groningen, Netherlands

Synopsis

Vessel architectural imaging (VAI) is explored in the heart by using a heartbeat-to-heartbeat GESE-EPI sequence upon injection of Gd-DTPA. Cardiac VAI can provide the vascular type, caliber, density and blood volume fraction indices in the myocardium, in line with previous work performed in the brain. Further histological validation of these indices is needed, but our initial results demonstrates the feasibility of this technique to advance cardiovascular research into cardiac microvascular dysfunction.

Introduction

Microvascular dysfunction is an increasing problem in cardiovascular risk population resulting in unexpected heart failure1. For cardiac MR applications the main vascular imaging technique is first-pass perfusion2, which does not provide specific microvascular information. For brain applications a vessel architectural imaging (VAI) approach has been demonstrated which allows for vascular volume, type, architecture and caliber calculations3–6. This technique is based on detecting the simultaneous dynamic T2- and T2*-changes upon injection of a contrast agent5-8. Cardiovascular research in microvascular disease could significantly benefit from the specific readouts that VAI provides. Recently, our group has shown that a GESE-EPI sequence gives stable heartbeat-to-heartbeat T2- and T2*-readouts9,10. In this work we evaluate this sequence in combination with a Gd-DTPA administration to determine the initial feasibility of cardiac VAI.

Methods

A cardiac GESE-EPI sequence was applied in 10 healthy participants on a 3T-Skyra (Siemens) with a Body-18/Spine-32 coil to acquire dynamic T2- and T2*-maps in a short-axis direction. The first dynamic acquisition was performed during a maximum achievable breath-hold to determine the oxygenation response9,10, while a second acquisition was done during a Gd-DTPA (0.2mmol/kg,4ml/s) injection over a second maximum achievable breath-hold for cardiac VAI. The sequence parameters were previously described9,10 including: 2-GE(TE)=9.8/23.6ms,3-mixed-SE(TE)=38.8/51.9/66.1ms,readout=120ms,TR=RR-interval. T2- and T2*-maps per heartbeat were created using an iterative parameter fitting technique11. The first TE of each heartbeat was used for segmentation and manually drawing a septal ROI12.

A dual-Gaussian membership function was fit to the dynamic $$$R_2=\frac{1}{T_2}$$$ and $$$R_2^*=\frac{1}{T_2^*}$$$ time-plots using nonlinear least-squares estimates. These estimated R2- and R2*-curves were combined in a time-parameterized “vessel vortex”3,6,13,14. According to previously published simulations3,5,6 the direction of the vortex indicates the predominant vessel type, which is either arterioles-dominated (clockwise), venule-dominated (counter-clockwise), or a mixture and/or capillaries (straight). Furthermore, the long axes of the vortex (Vindex) is proportional to the blood volume fraction (fBV)3,4, the slope of the vortex can be used to determine the relative vessel caliber index (Cindex)8,15 and also the relative vascular density index (Qindex) can be calculated16. $$V_{index}=\sqrt{\triangle {R_2^*}^2+\triangle {R_2}^2}$$ $$C_{index}=\frac{\triangle R_2^*}{\triangle R_2}$$ $$Q_{index}=\frac{\triangle R_2^*}{{\triangle R_2}^\frac{2}{3}}$$ Furthermore, a linear fit of the non-contrast-enhanced breath-hold data was used to determine the oxygenation response as a percentage in T2- and T2*-change9,10.

Results

The T2- and T2*-change due to the Gd-DTPA was seen over the time of the breath-hold(Fig1) and the voxel-by-voxel maximum $$$\triangle$$$R2- and $$$\triangle$$$R2*-maps are shown in Fig2a-b. These dynamic $$$\triangle$$$R2- and $$$\triangle$$$R2*-changes in an arterial(Fig2c) and a venule(Fig2d) ROI are shown to result in a clockwise and counter-clockwise vessel vortex, respectively. This vortex directional sensitivity results from the fact that temporal T2- and T2*-changes from Gd-DTPA are slightly shifted in time3,6(Fig3). Furthermore, all septal vortices from the 8 healthy participants that held their breath longer than the Gd-DTPA transit time, showed a predominant capillary architecture3(Fig4). In addition, the septal mean Vindex of 37.2±4.1s-1, mean Cindex of 2.9±0.3 and mean Qindex of 1.2±0.1s-1/3 were extracted(Fig5a-c). The oxygenation responses was determined as a T2- and T2*-change of 7.0±2.0% and 7.0±1.3%, respectively(Fig5d).

Discussion

The first results of cardiac VAI with GESE-EPI show a relative minor decrease in $$$\triangle$$$R2 and $$$\triangle$$$R2* after the peak, resulting in incomplete vortex-loops compared to brain3,6,14. Nevertheless, the Vindex, Cindex, Qindex and vortex direction could still be determined(Fig2cd). These cardiac indices need further dedicated validation with histology in the heart, but as an initial step we compared them with histological validated brain-based VAI parameters. Firstly, in rat brain a Cindex of 4.8±0.38 and a Qindex of 0.72±0.21s-1/3 17 have been reported, compared to a 9.2±1.0 and 0.36±0.1521s-1/3 in tumorous tissue, which are in the same range as our cardiac indices. However, from these histology validation studies it has been recommend to use the ratio between indices of healthy and diseased tissue8,15. This emphasizes the need of clinical comparison of cardiac VAI in a cardiovascular disease model. Secondly, absolute values for the slope length in the brain were not reported previously, but it has been shown that these slopes correlate with fBV as $$$V_{index}=10s^{-1}\sim{fBV=3}$$$%3,13. This correlation would result in a mean fBV of 11.2±1.2% in our septal ROIs, which could be expected18. Thirdly, a complete curved vortex cannot be achieved due to leakage of the Gd-DTPA in cardiac tissue, which makes extraction of the oxygen saturation (area of the vortex)3 difficult. Usage of a contrast-agent that remains intravascular19 might solve this leakage problem, but as shown the oxygenation response can also be determined from a non-contrast-enhanced breath-hold acquisition9,10. Lastly, these first cardiac VAI results correspond with the expected high cardiac microvascular density and fBV and showed vortex shapes that are expected in healthy tissue3,13. We hypothesize that these cardiac VAI results using a GESE-EPI sequence, despite the ultimate need for histological validation, can be further evaluated in healthy and myocardial-infracted tissue, which would provide insights in the clinical relevance of cardiac VAI.

Conclusion

We successfully performed VAI in the heart for the first time using a GESE-EPI sequence providing dynamic contrast-enhanced myocardial T2*- and T2-maps per heartbeat over the time of a single breath-hold. Obtained values for the dominant vascular type, Vindex, Cindex and Qindex were proportional to previous results in the brain, but still need further histological and clinical validation in cardiac diseases.

Acknowledgements

This work was supported in part by the Dutch Heart Association (2016T042) and by NIH research grants (R01EB019437, R01EB020613, R01MH116173, and U01EB025162)

References

1. Petersen JW, Pepine CJ. Microvascular Coronary Dysfunction and Ischemic Heart Disease – Where Are We in 2014? Trends Cardiovasc Med 2016;25:98–103. doi:10.1016/j.tcm.2014.09.013.

2. Nagel E, Greenwood JP, McCann GP, Bettencourt N, Shah AM, Hussain ST, et al. Magnetic resonance perfusion or fractional flow reserve in coronary disease. N Engl J Med 2019;380:2418–28. doi:10.1056/NEJMoa1716734.

3. Emblem KE, Mouridsen K, Bjornerud A, Farrar CT, Borra RJH, Wen PY, et al. Vessel Architectural Imaging Identifies Cancer Patient Responders to Anti-angiogenic Therapy. Nat Med 2013;19:1178–83. doi:10.1038/nm.3289.

4. Emblem KE, Farrar CT, Gerstner ER, Batchelor TT, Borra RJH, Rosen BR, et al. Vessel calibre—a potential MRI biomarker of tumour response in clinical trials. Nat Rev Clin Oncol 2014;11:566–84. doi:10.1038/nrclinonc.2014.126.

5. Troprès I, Grimault S, Vaeth A, Grillon E, Julien C, Payen JF, et al. Vessel size imaging. Magn Reson Med 2001;45:397–408. doi:10.1002/1522-2594(200103)45:3<397::AID-MRM1052>3.0.CO;2-3.

6. Kiselev VG, Strecker R, Ziyeh S, Speck O, Hennig J. Vessel size imaging in humans. Magn Reson Med 2005;53:553–63. doi:10.1002/mrm.20383.

7. Prinster A, Pierpaoli C, Turner R, Jezzard P. Simultaneous Measurement of ΔR2 and ΔR2* in Cat Brain during Hypoxia and Hypercapnia. Neuroimage 1997;6:191–200. doi:10.1006/nimg.1997.0288.

8. Dennie J, Mandeville JB, Boxerman JL, Packard SD, Rosen BR, Weisskoff RM. NMR imaging of changes in vascular morphology due to tumor angiogenesis. Magn Reson Med 1998;40:793–9. doi:10.1002/mrm.1910400602.

9. van den Boomen M, Manhard MK, Nguyen C, Emblem K, Slart R, Sosnovik D, et al. Heartbeat-to-Heartbeat Quantitative Cardiac BOLD Imaging during a Single Breath-Hold using a combined Gradient Echo-Spin Echo EPI sequence. SCMR 22nd Annu Sci Sess Focus Sess 6; Seatlle, USA 2019:2.

10. van den Boomen M, Snel GJH, Nguyen C, Manhard MK, Sosnovik D, Rudi AJO, et al. Heartbeat-to-Heartbeat Quantitative Myocardial Oxygenation Imaging within a Single Breath- Hold using a Combined Gradient Echo-Spin Echo EPI ( GESE-EPI ) Sequence in Patients with Hypertension. ISMRM 27th Annu Meet Exhib 2019:585.

11. Schmiedeskamp H, Straka M, Newbould RD, Zaharchuk G, Andre JB, Olivot JM, et al. Combined spin- and gradient-echo perfusion-weighted imaging. Magn Reson Med n.d.;68:30–40. doi:10.1002/mrm.23195.

12. Heiberg E, Sjögren J, Ugander M, Carlsson M, Engblom H, Arheden H. Design and validation of Segment - freely available software for cardiovascular image analysis. BMC Med Imaging 2010;10. doi:10.1186/1471-2342-10-1.

13. Zhang K, Yun SD, Triphan SMF, Sturm VJ, Buschle LR, Hahn A, et al. Vessel architecture imaging using multiband gradient-echo/spin-echo EPI. PLoS One 2019;14:e0220939. doi:10.1371/journal.pone.0220939.

14. Xu C, Kiselev VG, Möller HE, Fiebach JB. Dynamic hysteresis between gradient echo and spin echo attenuations in dynamic susceptibility contrast imaging. Magn Reson Med 2013;69:981–91. doi:10.1002/mrm.24326.

15. Farrar CT, Kamoun WS, Ley CD, Kim YR, Kwon SJ, Dai G, et al. In vivo validation of MRI vessel caliber index measurement methods with intravital optical microscopy in a U87 mouse brain tumor model. Neuro Oncol 2010;12:341–50. doi:10.1093/neuonc/nop032.

16. Wu EX, Tang H, Jensen JH. High-resolution MR imaging of mouse brain microvasculature using the relaxation rate shift index Q. NMR Biomed 2004;17:507–12. doi:10.1002/nbm.921.

17. Jensen JH, Chandra R. MR imaging of microvasculature. Magn Reson Med 2000;44:224–30. doi:10.1002/1522-2594(200008)44:2<224::AID-MRM9>3.0.CO;2-M.

18. McCommis KS, Goldstein TA, Zhang H, Misselwitz B, Gropler RJ, Zheng J. Quantification of myocardial blood volume during dipyridamole and doubtamine stress: A perfusion CMR study. J Cardiovasc Magn Reson 2007;9:785–92. doi:10.1080/10976640701545206.

19. Gale EM, Wey HY, Ramsay I, Yen YF, Sosnovik DE, Caravan P. A manganese-based alternative to gadolinium: Contrast-enhanced MR angiography, excretion, pharmacokinetics, and metabolism. Radiology 2018;286:865–72. doi:10.1148/radiol.2017170977.

Figures

Gif of baseline and peak T2-and T2*-maps of the myocardium, pre- and post- the pass of a gadolinium-DTPA bolus. Here both the T2- and T2*-values over the whole myocardium are decreased at time of the bolus peak, as compared to the baseline.

The a) ΔR2-and b) ΔR2*-maps of the peak change after Gd-DTPA injection show variation over the myocardium, with higher ΔR2- and ΔR2*-changes near the LAD, RCA and LCX locations. The bolus Gd-DTPA causes the ΔR2* to peak before the ΔR2 in the fast flowing arteries c), resulting in a clockwise vortex, while in slow inflow venules the ΔR2 peak appears before the ΔR2* peak, resulting in a counterclockwise vortex d).

A cartoon of the Gd-DTPA in the vasculature shows the change in field inhomogeneity of the Gd present in the vasculature itself causing a T2*-relaxation change, while the effect of the Gd on the protons in the surrounding tissue and the vasculature itself causes the T2-relaxation change. In fast flowing artery the T2* effect occurs before the T2 effect, while in slow flowing venules this is the other way around.

Eight healthy participants showed a similar straight or slightly counterclockwise rotating vortex in their septal ROI, indicating (according to previously performed Monte-Carlo simulations3,19) predominant capillaries over macro-vasculature, which is also expected in healthy tissue14.

a) Mean vascular volume index (Vindex) of 37.2±4.1s-1, b) vessel caliber index (Cindex) of 2.9±0.3, c) vascular density index (Qindex) of 1.2±0.1s-1/3 and d) mean oxygenation response change of T2 and T2* of 7.0±2.0% and 7.0±1.3%, respectively, from all 8 healthy volunteers

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