The Impact of Arterial Input Function Determination Variation on Prostate Dynamic Contrast-Enhanced Magnetic Resonance Imaging Pharmacokinetic Modeling: A Multicenter Data Analysis Challenge
Wei Huang1, Yiyi Chen1, Andriy Fedorov2, Xia Li3, Guido Jajamovich4, Dariya I Malyarenko5, Madhava Aryal5, Peter S LaViolette6, Matthew J Oborski7, Finbarr O’Sullivan8, Richard G Abramson9, Mark Muzi10, Kourosh Jafari-Khouzani 11, Aneela Afzal1, Alina Tudorica1, Brendan Moloney1, Cecilia Besa4, Jayashree Kalpathy-Cramer11, James M Mountz7, Charles M Laymon7, Kathleen Schmainda6, Yue Cao5, Thomas L Chenevert5, Bachir Taouli4, Thomas E Yankeelov9, Fiona Fennessy2, and Xin Li1

1Oregon Health & Science University, Portland, OR, United States, 2Brigham and Women’s Hospital and Harvard Medical School, Boston, MA, United States, 3General Electric Global Research, Niskayuna, NY, United States, 4Icahn School of Medicine at Mount Sinai, New York, NY, United States, 5University of Michigan, Ann Arbor, MI, United States, 6Medical College of Wisconsin, Milwaukee, WI, United States, 7University of Pittsburgh, Pittsburg, PA, United States, 8University College Cork, Cork, Ireland, 9Vanderbilt University, Nashville, TN, United States, 10University of Washington, Seattle, WA, United States, 11Massachusetts General Hospital and Harvard Medical School, Boston, MA, United States

Synopsis

Dynamic Contrast-Enhanced MRI (DCE-MRI) pharmacokinetic modeling is widely used to extract tissue specific quantitative parameters. However, the accuracy and precision of these parameters can be affected by many factors, with arterial input function (AIF) determination being a primary source of uncertainties. In this multicenter study, we sought to evaluate variations in DCE-MRI parameters estimated from shared prostate DCE-MRI data as a result of differences in AIFs.

Introduction

Pharmacokinetic analysis of Dynamic Contrast-Enhanced MRI (DCE-MRI) data allows extraction and mapping of quantitative parameters of tissue biology in vivo, such as Ktrans, ve, and kep. However, the accuracy and precision of these parameters can be affected by many factors, with arterial input function (AIF) determination being a primary element. In this multicenter study, we sought to evaluate variations in DCE-MRI parameters estimated from shared prostate DCE-MRI data as a result of differences in AIFs.

Methods

Prostate DCE-MRI data were collected on 11 subjects at one center and shared with eight other Quantitative Imaging Network (QIN) centers. The data acquisition details have been previously reported 1. Each of the 9 centers determined individual AIFs with its site-specific method (fully automated, semi-automated, manual, etc.). All individual AIF files were submitted to a managing center (one of the 9 centers), which performed DCE-MRI pharmacokinetic modeling with all AIFs using the Tofts model (TM) 2. A previously published population AIF 3 was also included in the analysis for comparison. In addition, a reference tissue (using obturator muscle in this study) approach 4 was also used for concentration adjustment of all individual AIFs. The same 11 data sets were then modeled with the reference-tissue-adjusted AIFs, resulting in a combined total of twenty sets of pharmacokinetic parameters (Ktrans, ve, and kep) for each subject, which were then classified into AIF-unadjusted (unadj.) and AIF-reference-tissue adjusted (adj.) groups. To insure variations in derived kinetic parameters were from AIF variations only, all other parameters for the model fitting including tumor ROI definition and pre-contrast T1 were kept the same. Mean tumor Ktrans, ve, and kep (= Ktrans/ve) values for each patient were summarized. To assess agreements and variations of the DCE-MRI parameters obtained with different AIFs within a group (unadj. and adj.), intra-class correlation coefficient (ICC) and within-subject coefficient of variation (wCV) were computed, respectively.

Results

The top image of Figure 1 zoomed to the prostate area shows a DCE-MRI slice of one subject. The cyan-colored ROI demarks the lesion area for subsequent TM modeling and parameter comparisons. Ktrans color maps generated by TM analysis using unadjusted AIFs from the 9 centers (numbered 1 to 9) are shown in the middle panels and those with reference-tissue-adjusted AIFs are shown at the bottom. Substantial variations, mostly in the magnitude of voxel Ktrans but not in the distribution pattern, can be seen among the Ktrans maps obtained with different unadjusted AIFs. These differences are lessened when the Ktrans parameter was derived with reference-tissue-adjusted AIFs. Figure 2 shows the column graph of wCV values for the Ktrans, ve and kep parameters obtained with the unadjusted (shaded light grey) and adjusted (dark grey) AIFs. The respective 95% confidence intervals (CI) are shown as error bars. Figure 3 shows the column graph of ICC values with the same layout as Fig. 2.

Discussion and Conclusion

Results from this multicenter study clearly show that variations in AIF quantification result in variations in pharmacokinetic parameter values estimated from prostate DCE-MRI data, with Ktrans and ve having the largest and smallest variations, respectively. The agreement in the Ktrans and ve parameters obtained with muscle-reference-region-adjusted AIFs is improved compared to that with unadjusted AIFs (Figs. 1-3). It is important to note that DCE-MRI parameter variations caused by AIF variations are mostly systematic. As shown in Fig. 1, the differences among the prostate tumor Ktrans maps obtained with different AIFs are mostly in voxel Ktrans values. The pattern of voxel Ktrans distribution largely remains similar. This suggests that assessment of tumor heterogeneity through texture analysis of DCE-MRI parametric maps may not be affected greatly by variations in AIF determination. Consequently, despite the significant challenge in accurate AIF determination, quantitative DCE-MRI remains a promising imaging tool for tumor characterization and therapeutic monitoring. This observation needs further validation in future studies. kep variation is shown to be less sensitive to AIF uncertainty than Ktrans, suggesting it may be a more robust pharmacokinetic parameter for characterizing prostate microvasculature. In addition, parameters obtained using individually determined AIF demonstrates better consistency (higher ICC values, etc.) than those obtained with a generic population AIF (often derived from a different acquisition protocol) for substantial contrast agent extravasation situations like that in the prostate.

Acknowledgements

Grant Support: NIH U01-CA154602, U01-CA151261, U01-CA183848, U01-CA154601, 5U01-CA148131, U01-CA176110, U01-CA172320, U01-CA142565, U01-CA166104, U01-CA140230.

References

1. Hegde et al. J Magn Reson Imaging. 2013;37:1035-54. 2. Tofts et al. J Magn Reson Imaging 1999;10:223-232. 3. Parker et al. Magn Reson Med. 2006 Nov;56(5):993-1000. 4. Li et al. Magn Reson Med. 2013 Mar 27;69:171-178.

Figures

Figure 1. a shows a cropped DCE-MRI image of the prostate area. The cyan-colored ROI demarks the lesion area for subsequent TM modeling. b and c show lesion Ktrans color maps by TM analysis of the data using unadjusted AIFs determined by the 9 centers and reference-tissue-adjusted AIFs, respectively.

Figure 2. Column graph of wCV for the Ktrans, ve and kep parameters obtained with the unadjusted (shaded light grey) and adjusted (dark grey) AIFs. The respective 95% confidence intervals (CI) are shown as error bars.

Figure 3. Column graph of ICC for the Ktrans, ve and kep parameters obtained with the unadjusted (shaded light grey) and adjusted (dark grey) AIFs. The respective 95% confidence intervals (CI) are shown as error bars.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0580