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 K
trans, v
e, and k
ep. 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 (K
trans, v
e, and k
ep) 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 T
1 were kept the same.
Mean
tumor K
trans, v
e, and k
ep (= K
trans/v
e)
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. K
trans 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 K
trans but not in the distribution pattern, can be seen
among the K
trans maps obtained with different unadjusted AIFs. These
differences are lessened when the K
trans parameter was derived with
reference-tissue-adjusted AIFs.
Figure 2
shows the column graph of wCV values for the K
trans, v
e
and k
ep 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 K
trans
and v
e having the largest and smallest variations, respectively. The
agreement in the K
trans and v
e 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 K
trans maps obtained with different AIFs are mostly
in voxel K
trans values. The
pattern of voxel K
trans 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. k
ep
variation is shown to be less sensitive to AIF uncertainty than K
trans,
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.