DCE-MRI High-resolution Metabolic Prostate Imaging is Insensitive to AIF Uncertainty
Xin Li1, Mark G. Garzotto2,3, Fergus V. Coakley4, Brendan Moloney1, William J. Woodward1, Yiyi Chen5, Wei Huang1, William D. Rooney1, and Charles S. Springer, Jr.1

1Advanced Imaging Research Center, Oregon Health & Science University, Portland, OR, United States, 2Portland VA Medical Center, Portland, OR, United States, 3Urology, Oregon Health & Science University, Portland, OR, United States, 4Department of Diagnostic Radiology, Oregon Health & Science University, Portland, OR, United States, 5Division of Biostatistics, Dept. of Public Health and Preventive Medicine, Knight Cancer Institute, Oregon Health and Science University, Portland, OR, United States

Synopsis

Accurate arterial input function (AIF) measurement in Dynamic Contrast Enhanced MRI (DCE-MRI) remains challenging. This hinders DCE-MRI’s wider adoption. Since the contrast reagent (CR) is detected indirectly through water proton R1 relaxation rate constant change, DCE-MRI intrinsically works as a dual-probe (CR and water) method. In this study, we demonstrate that while the common pharmacokinetic parameters associated with CR extravasation are highly sensitive to AIF accuracy, the transcytolemmal water exchange parameter is not. With the recent correlation of water exchange kinetics and cellular metabolic activity, this current work demonstrates the practicability of high-resolution metabolic imaging of the prostate.

Purpose

Accurate arterial input function (AIF) measurement in prostate Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI) remains challenging. This hinders DCE-MRI’s wider clinical adoption. Since Gadolinium based contrast reagent (CR) is detected indirectly through water proton R1 relaxation rate constant change, DCE-MRI is intrinsically a dual-molecular probe (CR and water) technique. In this study, we demonstrate that while the most common pharmacokinetic model parameters like Ktrans (CR transfer constant) and ve (extravascular, extracellular, volume fraction) are highly AIF-dependent for prostate DCE‑MRI, the transcytolemmal water exchange parameter is not. With the recent correlation of water exchange kinetics and cellular metabolic activity, this work demonstrates the practicability of a high-resolution metabolic imaging approach for the prostate.

Methods

Prostate DCE-MRI data were acquired on 13 subjects with a Siemens TIM Trio (3T) system under an IRB-approved protocol. RF receiving used a combination of Spine Matrix and flexible Body Matrix coil arrays. The DCE-MRI acquisition employed a 3D TurboFLASH pulse sequence with a 256*144*16 matrix size and a 360*203 mm2 FOV; slice thickness 3 or 3.2 mm. Imaging parameters are: TR/TE/FA, 5.0 ms/1.57 ms/15º; and image frame sampling interval is 6.3 s. Other details are reported previously.1 All subjects subsequently underwent standard ultrasound-guided systematic ten-core prostate biopsies. Malignancy was found in 5 subjects. One region of interest (ROI) was selected for each subject, resulting in 5 malignant and 8 benign ROI time-courses.1 The water exchange pharmacokinetic model assuming a single compartment of water signal2 extracts Ktrans, ve, and τi (mean intracellular water lifetime) parameters. Monte Carlo simulations were used to investigate parameter changes associated with AIF magnitude variation. The “true” AIF (AIF0) was that determined directly from the data1. A noiseless tissue time-course was first generated based on data acquisition detail and pharmacokinetic equations. Simulations then fitted pharmacokinetic parameters to the time-course with various AIF amplitudes scaled from 0.8 to 2.0 of AIF0 (step-size 0.1). For each amplitude-adjusted AIF (e.g., AIF1 = 130%∙AIF0), four hundred fittings were performed with random Gaussian noise added to the simulated time-course and a different initial guess parameter value set each time.

Results

The Figure 1 bar graph shows the Ktrans, ve, and τi accuracy (mean) and precision (SD) changes responding to a 30% AIF amplitude increase (AIF1). Ktrans = 0.4 min-1, ve = 0.25, and τi = 0.45 s base values were used to generate the simulated tissue DCE time-course. Ktrans and ve numerical values decreased nearly 30% with a 30% AIF magnitude increase, while the change for τi remains negligibly small (< 3%). The τi reciprocal (τi-1) is the unidirectional rate constant for equilibrium cellular water efflux (kio).4 Figure 2a shows a portion of an axial T2-weighted pelvic image (anterior up) of a 66 year old subject with left prostate malignancy. Figures 2b – 2d show zoomed parametric maps of a tumor-encompassing ROI in the left prostate indicated by the 2a yellow rectangle. In Figure 3 lesion ROI-averaged in vivo kio values are plotted against subsequently obtained ex vivo Gleason scores, for the five malignant cases. A clear trend of increasing kio with Gleason score (Spearman’s correlation coefficient, r = 0.908; p = 0.0333) is observed (the blue dashed curve is intended only to guide the eye). The error bars show kio standard errors and estimated Gleason score errors.3

Discussion

Although the AIF is hard to measure accurately from the CR bolus first-pass, an AIF error of less than 30% can be readily managed from the washout phase, taking into account the total blood volume and CR elimination rate. Furthermore, if kio increases are dominated by increases in the mean Na+,K+-ATPase [NKA] turnover for the cells in the voxel,4 this work suggests significant impact and unprecedented insight into prostate tumor progression in vivo. Active trans-membrane water cycling accompanies the continuous, homeostatic active trans-membrane ion and osmolyte cycling enzymatic activity, driven by the membrane P-type ATPase ion pump, constitutively characteristic of a living cell.4 Fig. 3 indicates a greater Gleason score correlates with greater mean NKA activity. Taken together, the kio map shows in vivo high-resolution metabolic imaging that is insensitive to AIF uncertainty (thus easy to quantify in that sense). It is important to point out that precision in the kio map is greatest when CR extravasation is extensive (large Ktrans),4 as always seen in prostate.1

Acknowledgements

Grant Support: NIH: RO1-NS040801, UO1-CA154602. R44 CA180425.

References

1. Li, et al. Magn Reson Med., 69:171-178, 2013. 2. Yankeelov, et al. Magn Reson Med ., 50:1151-1169, 2003. 3. Goodman et al. The Prostate 72:1389-1398, 2012. 4. Springer, et al. NMR Biomed. 27:760-73, 2014.

Figures

Figure 1. Bar graphs show the Ktrans, ve, and τi accuracy (mean) and precision (SD) changes responding to a 30% AIF amplitude increase (AIF1). A 30% magnitude increase of the AIF0 resulted in nearly 30% Ktrans and ve numerical value reductions, while the change for τi remains negligibly small.

Figure 2. a shows an axial T2-weighted pelvic image of a subject with left prostate malignancy. 2b – 2d show zoomed parametric maps of a tumor-encompassing ROI within the yellow rectangle of a. kep=Ktrans/ve and vi=1-ve. Dark frames surrounding the lesion ROI color maps are merely to provide image contrast.

Figure 3. Lesion ROI-averaged in vivo kio (1/τi) values are plotted against subsequently obtained ex vivo Gleason scores for the five malignant cases. The blue dashed curve is intended only to guide the eye. The error bars show kio standard errors and estimated Gleason score uncertainties.



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