2096

Flow Effects in AIF Measurements for Quantitative Myocardial Perfusion
Ye Tian1, Qi Huang1, Jason Mendes1, Ganesh Adluru1, and Edward DiBella1
1Radiology and Imaging Sciences, University of Utah, Salt Lake City, UT, United States

Synopsis

We study the flow effects in AIF estimation for myocardial perfusion MRI. Simulations were performed and radially acquired datasets were retrospectively analyzed. We reconstruct images using different number of rays and estimated AIF from each reconstruction. The results show an increase in peak AIF values with a longer acquisition window. For acquisition window of 125ms and 187ms, we found an averaged 10.6% and 22.9% increase when compared with acquisition window of 62ms. A low acquisition window is preferred to acquire an accurate AIF for quantitative myocardial perfusion MRI.

Introduction

Myocardial first-pass perfusion MRI can provide non-invasive, radiation-free diagnoses for coronary artery diseases. There is also increasing interest in getting quantitative myocardial blood flow (MBF) values from the first-pass perfusion images. Quantitative analysis requires an accurate estimation of the arterial input function (AIF). Popular methods for AIF estimations include the dual-bolus (1) or the dual-sequence (2) methods. However, the AIF can change across a cardiac cycle, and the diastolic AIF peak has been reported to be significantly higher than the systolic AIF peak (3); this may be caused by flow effects during the data acquisition. Here we study the flow effects with simulations and by retrospectively analyzing two types of human subjects datasets acquired with radial trajectories (4,5).

Theory

A commonly used pulse sequence for quantitative myocardial perfusion is saturation recovery preparation followed by a FLASH readout. During a typical acquisition window of 70-150ms, ‘fresh’ blood flows into or out of the LV if not during the quiescent phase of diastole which is relatively short. The mixture of blood that has seen a different number of slice selective alpha pulses can result in a signal intensity that is higher than still blood, giving a higher AIF estimation.
Radial trajectories go through the k-space center at every ray, such that each ray reflects the image contrast at that temporal point. We reconstructed radially acquired k-space data using a different number of rays starting from the first ray after each SR pulse and estimated AIF from each reconstruction. The AIFs estimated from different reconstructions should reflect different flow effects, with increasing flow effects in longer acquisition window.

Methods

Simulations
An in-vivo measured AIF [Gd] curve (acquisition window 47ms) was used in Bloch simulations as the ground truth. Different acquisition windows (31,47,62,125 and 187ms) combined with different constant flow rates (0.2,0.4,0.6,0.8 and 1m/s) were simulated. We chose TR=2.6ms, SRT=10ms, flip angle (FA)=14°.

Saturation recovery 72-ray datasets
Eleven patients with 46 injections (25 at stress) were acquired with TR/TE=2.6/1.1ms, SRT=10ms, FA=14°, 72 rays/frame, 2-4 short-axis slices per injection on a 3T scanner (5). Proton density (PD) weighted images were also acquired for each patient.

Continuous radial SMS (cSMS) datasets
Ten injections (3 were stress) from 7 patients acquired with TR/TE=2.5/1.2ms, SRT=10ms, and FA=14° were analyzed. Three slices were acquired simultaneously with a CAIPIRINHA phase modulation (6,7). A saturation pulse was applied after each ECG trigger, and radial data acquisition continued until the next trigger (4).

AIF estimation
For each dataset, five sets of AIFs were estimated from images reconstructed using different number of rays per time frame. We chose 12,18,24,48 and 72 rays/frame, correspond to acquisition windows of 31, 47, 62, 125 and 187ms.
For each injection, a region of interest (ROI) was drawn in the LV and the pixel intensities in the ROI were averaged to be the signal intensity (SI) curve. A look-up-table generated from Bloch simulations was used to convert measured SI curves (normalized by PD images) into gadolinium concentration curves.

Results

Simulation
Figure 1 shows the simulation results with different acquisition windows combined with different flow rates. The flow effects can increase the AIF estimation up to 70% with the highest flowrate and the longest acquisition window. However, with a short acquisition window (31 or 47ms), the overestimation was less than 11%.

In-vivo measurements
Figure 2 shows four typical AIFs at different injection doses and different physiological states (stress or rest), where the longer acquisition led to a higher AIF. Figure 3 summarizes all the results, including 148 sets of AIFs (each set has 5 AIFs with different acquisition windows) from 56 injections (2-4 slices per injection) of 18 patients. We used the 62ms AIF as a reference to calculate the percentage difference:
$$\mathrm{Percentage\quad Difference}=\frac{[\mathrm{Gd}]_{\mathrm{peak}}-[\mathrm{Gd}]_{\mathrm{peak}}^{\mathrm{ref}}}{[\mathrm{Gd}]_{\mathrm{peak}}^{\mathrm{ref}}},$$
where $$$[\mathrm{Gd}]_{\mathrm{peak}}$$$ is averaged 3 highest $$$[\mathrm{Gd}]$$$ values in each AIF curve. The $$$[\mathrm{Gd}]_{\mathrm{peak}}$$$ shows a small trend of increase (<5%) at low acquisition windows (31-62ms). With a longer acquisition window, the increase is significant, at 10.6% for 125ms and 22.9% for 187ms. The increase does not have a significant difference between stress and rest.

Discussion and Conclusion

The study was done on radially acquired datasets. With Cartesian or spiral trajectories, the flow effects are different. With a centric Cartesian readout, arguably the flow effects can be minimal, however the image may suffer from an inherent high-pass filtering (8). With a linear Cartesian readout, the flow effects come from the k-space center readout, and may be similar with radial trajectories that average flow effects from all readouts. The 12-ray and 18-ray reconstructions suffer from artifacts as shown in Figure 4, but they may still provide useful AIF estimations, since averaging pixels in LV increase the signal-to-noise level. We also studied AIF difference across the apex, middle and base LV from 42 injections that have at least 3 slices and show the results in Figure 5.
The AIF estimation form the LV for quantitative myocardial perfusion is higher with an increased acquisition window; this is likely related to flow effects. A short acquisition window below or around 62ms and/or acquisition during a period of minimal LV blood flow is desired to acquire AIF with little flow effects when using a radial trajectory.

Acknowledgements

No acknowledgement found.

References

1. Christian TF, Aletras AH, Arai AE. Estimation of absolute myocardial blood flow during first-pass MR perfusion imaging using a dual-bolus injection technique: comparison to single-bolus injection method. J Magn Reson Imaging 2008;27(6):1271-1277.

2. Gatehouse PD, Elkington AG, Ablitt NA, Yang GZ, Pennell DJ, Firmin DN. Accurate assessment of the arterial input function during high-dose myocardial perfusion cardiovascular magnetic resonance. J Magn Reson Imaging 2004;20(1):39-45.

3. Radjenovic A, Biglands JD, Larghat A, Ridgway JP, Ball SG, Greenwood JP, Jerosch-Herold M, Plein S. Estimates of systolic and diastolic myocardial blood flow by dynamic contrast-enhanced MRI. Magn Reson Med 2010;64(6):1696-1703.

4. DiBella E, Likhite D, Adluru G, Welsh C, Wilson B. A New Hybrid Approach for Quantitative Multi-Slice Myocardial DCE Perfusion. 2016; Singapore. In the proceddings of ISMRM 2016 annual meeting. p 463.

5. Kim TH, Pack NA, Chen L, DiBella EV. Quantification of myocardial perfusion using CMR with a radial data acquisition: comparison with a dual-bolus method. J Cardiovasc Magn Reson 2010;12:45.

6. Wang H, Adluru G, Chen L, Kholmovski EG, Bangerter NK, DiBella EV. Radial simultaneous multi-slice CAIPI for ungated myocardial perfusion. Magn Reson Imaging 2016;34(9):1329-1336.

7. Yutzy SR, Seiberlich N, Duerk JL, Griswold MA. Improvements in multislice parallel imaging using radial CAIPIRINHA. Magn Reson Med 2011;65(6):1630-1637.

8. Kellman P, Hansen MS, Nielles-Vallespin S, Nickander J, Themudo R, Ugander M, Xue H. Myocardial perfusion cardiovascular magnetic resonance: optimized dual sequence and reconstruction for quantification. J Cardiovasc Magn Reson 2017;19(1):43.

Figures

Figure 1. Simulation results. The figure shows simulation results of different flow rates combined with different acquisition windows. Subfigures (a-d) show comparisons at the same acquisition window with different flow rates. Subfigures (e-g) show comparisons at the same flow rate but with different acquisition windows. Figure (h) illustrates the peak [Gd] percentage increase due to inflow with different flow rates and acquisition windows. Simulation results suggest that using a short acquisition window (31 or 47ms) can keep the AIF overestimation within 11%.

Figure 2. In-vivo AIF measurements. The figure shows 4 examples of the measured AIFs. Subfigure (a) shows an AIF of a low dose (0.004 mmol/kg) injection at stress. (b) shows an AIF of a standard dose (0.06 mmol/kg) injection at stress. (c) shows an AIF of a low dose (0.004 mmol/kg) injection at rest. (d) shows an AIF of a standard dose (0.02 mmol/kg) injection at stress. The AIFs from shorter acquisition windows (31, 47 and 62ms) share similar values, however the long acquisition window AIFs (125 and 187ms) are significantly higher. The AIFs of 31ms suffer from poor image quality.

Figure 3. The figure shows the peak [Gd] percentage difference when using AIF with a 62ms acquisition window as the reference. Subfigure (a) shows all AIFs, and (b) and (c) show rest and stress AIFs respectively. The dark lines illustrate the averaged values, and the dark areas are within one standard deviation. Each dot represents one AIF. The top numbers indicate the averaged percentage differences. The 125ms and 187ms acquisitions are higher than the 62ms acquisition by 10.6% and 22.9%, respectively. There is no significant difference between stress and rest scans.

Figure 4. Example images with different reconstruction windows. From top to bottom, the rows show a 72-ray dataset with dose 0.02mmol/kg (blue box), a 72-ray dataset with dose 0.004mmol/kg (orange box) and a cSMS dataset with dose 0.075mmol/kg (yellow box). The three slices shown in the cSMS dataset were acquired simultaneously. The shorter the reconstruction window, the poorer the image quality. However, given the large cavity size in the LV, it is possible to estimate an AIF out of each image.

Figure 5. T-test results. This table shows pairwise t-test results between different slices. The results suggest there is no significant difference between different slices. However, the reduced p-value at longer acquisition windows (125ms and 187ms) are suspected to be different flow patterns at different slice positions.

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