Abstract:C. Yang, W. M. Stadler, M. Medved, G. S. Karczmar University of Chicago, Chicago, IL, United States Introduction: Accurate estimation of the arterial input function (AIF) is essential for quantitative analysis of DCE-MRI data. We have previously demonstrated that a common AIF can be estimated from the contrast agent concentration versus time curves (Ct(t)) in two reference tissues with different enhancement patterns (1). This double reference tissue method (DRTM) is based on the two assumptions that (a) each tissue’s AIF has the same shape, but may have a different arrival time; (b) we know a priori what kinetic model can adequately describe the Ct(t). Assumption (b) of the DRTM has limited its applications. For a particular tumor, it is usually unknown what kinetic model is adequate (2). This makes the tumors generally unfit as the reference tissues in the DRTM. Now we generalize the DRTM to a multiple reference tissue method (MRTM). Materials and Methods: In the MRTM, which uses two or more reference tissues with unknown kinetic parameters to estimate the AIF, we use an automatic model selection scheme to drop Assumption (b) in the DRTM. We start from a simple kinetic model, for example a two-compartment model (3) without inclusion of the blood plasma volume parameter, νp. With the simple model, we estimate the AIF and obtain the goodness of the fit to the data. Then we increase the model complexity for example by including the parameter νp. Using the new model, we will estimate the AIF again and obtain the goodness of the fit by the new AIF. We will stop increasing the model complexity when the improvement in the goodness of the fit is statistically insignificant. With this scheme, the adequate kinetic model for the reference tissues is determined from the data themselves. By using a non-parametric smoothing method (4) and modifying a Fast Fourier Transform algorithm (5), we are able to significantly increase the accuracy of the algorithm. We apply the MRTM to the DCE-MRI data of 8 renal cell cancer patients from an on-going clinical trial on the cancer drug Sorafenib. T1 weighted images of bone metastases in the pelvic region were scanned with 2 second resolution for about 7 minutes pre and post bolus injection of Gd-DTPA. Using a cluster analysis (6), voxels in the tumor ROI are divided into subsets such that all voxels within each subset have similar Ct(t). The tumor heterogeneity leads to many subsets having a significantly different mean Ct(t) (Fig. 1a) so that each subset can be considered as a reference tissue. Selecting 3-8 subsets as the “multiple reference tissues”, the AIF was then estimated with the MRTM using the two-compartment model. For each tumor the automatic model selection scheme was used to determine whether the inclusion of the blood plasma volume parameter was necessary. The adjusted Chi-square, defined as the sum of squared errors divided by the variance of the measurement noise and the degree of freedom, is used to indicate the goodness of the fit. For comparison, we also used a fixed bi-exponential AIF (3) to fit the tumor data. Results: One representative example is shown in Fig. 1. In 5 out of 8 patients, the contribution of the blood plasma to the measured contrast agent concentration in tumors was found significant. The AIFs estimated by the MRTM generally demonstrate both a first pass peak and a washout phase. Even the second pass of the contrast agent bolus can be shown with striking details due to the high signal to noise ratio (SNR) of the fitted AIF (Fig. 1b). The AIFs estimated by the MRTM all closely resemble published AIFs determined by blood sampling (7). The adjusted chi-square of the fit by the AIFs estimated with the MRTM is usually over 10 times less than that by the fixed bi-exponential AIF. Due to the improved fit by the estimated AIF, the kinetic parameters of the tumor estimated on a voxel-by-voxel basis were found to have high SNR (Fig 1c-1e). Discussions and Conclusions: The MRTM can be generally applied in DCE-MRI studies to estimate the local AIF of tumors using the heterogeneous tumor itself as a source of “multiple reference tissues”. The estimated local AIFs closely resembled published AIFs and significantly improved the goodness of the fit. Further studies are needed to investigate whether this method can provide more reproducible quantitative DCEMRI parameters to assess tumor perfusion and permeability, hence improving DCE-MRI based cancer diagnostics and therapeutic monitoring. References: 1. Yang C et al, MRM 2004; 52:1110. 2. Port RE et al, MRM 1999; 10:233. 3. Tofts PS et al, MRM 1991; 17:357. 4. Fan J and Gijbel I, Local polynomial modeling and its applications, 1996. 5. Riabkov DY and Di Bella EV, IEEE Trans Biomed Eng 2002; 49:1318. 6. Everitt BS, Cluster analysis, 1993. 7. Fritz-Hansen T et al, MRM 1996; 36:225.

Joint estimation of compartment-specific T2 relaxation and tumor microstructure using multi-TE IMPULSED MRI

Variations of quantitative perfusion measurement on dynamic contrast enhanced CT for colorectal cancer: implication of standardized image protocol.

Three-dimensional Simultaneous Brain Mapping of T1, T2, T2*and Magnetic Susceptibility with MR Multitasking

Three-dimensional simultaneous brain T 1 , T 2 , and ADC mapping with MR Multitasking.

Retrospective T2 Quantification from Conventional Weighted MRI of the Prostate Based on Deep Learning

Compartment-specific estimation of T2 and T2* with diffusion-PEPTIDE MRI

Assessing tumor microstructure with time‐dependent diffusion imaging: Considerations and feasibility on clinical MRI and MRI‐Linac

Three-dimensional simultaneous brain T<sub>1</sub>, T<sub>2</sub>, and ADC mapping with MR Multitasking

A multicompartment model for intratumor tissue-specific analysis of DCE-MRI using non-negative matrix factorization

Distortion-Free Diffusion MRI using an MRI-Guided Tri-Cobalt 60 Radiotherapy System: Sequence Verification and Preliminary Clinical Experience.

SU-D-303-02: Impact of Arterial Input Function Selection and T10 Correction On DCE-MRI Tumour Response Prediction Using Compared to Volumetric DCE CT

Partially Independent Component Analysis of Tumor Heterogeneities by DCE-MRI.

Consideration of transmembrane water exchange in pharmacokinetic model significantly improves the accuracy of DCE-MRI in estimating cellular density: A pilot study in glioblastoma multiforme

Quantification of the in vivo brain ultrashort‐T2* component in healthy volunteers

Comprehensive characterization of tumor therapeutic response with simultaneous mapping cell size, density, and transcytolemmal water exchange

Tumor Arterial Input Function ( AIF ) estimation in DCE-MRI studies using a multiple reference tissue method

Magnetic resonance imaging of mean cell size in human breast tumors

Voxel-wise mapping of DCE-MRI time-intensity-curve profiles enables visualizing and quantifying hemodynamic heterogeneity in breast lesions

T 2 ' Mapping of the Brain from Water-Unsuppressed 1 H-MRSI and Turbo Spin-Echo Data.

Model discovery approach enables noninvasive measurement of intra-tumoral fluid transport in dynamic MRI

Improving Microstructural Estimation in Time‐Dependent Diffusion MRI With a Bayesian Method