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.

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

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

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

Restriction-induced time-dependent transcytolemmal water exchange: Revisiting the Kärger exchange model

Restriction-Induced Time-Dependent Transcytolemmal Water Exchange: Revisiting the Kӓrger Exchange Model

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

Water exchange detected by shutter speed dynamic contrast enhanced-MRI help distinguish solitary brain metastasis from glioblastoma

Correlative assessment of perfusion-related parameters of intravoxel incoherent motion with dynamic contrast ;enhancement in MR imaging of gliomas

Abstract A024: DCE-MRI-based biophysical simulation to forecast NAT response in HER2+ breast cancer patients, with glucose characterization and orthogonal validation using FDG-PET and 64Cu studies

Cellular EXchange Imaging (CEXI): Evaluation of a diffusion model including water exchange in cells using numerical phantoms of permeable spheres

Utilizing Dynamic Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) to Analyze Interstitial Fluid Flow and Transport in Glioblastoma and the Surrounding Parenchyma in Human Patients

The Role of T1-weighted Dynamic Contrast-Enhanced Perfusion Magnetic Resonance Imaging in Gross Target Volume Delineation of Glioma Patients

High-Resolution Diffusion-Weighted Imaging of C6 Glioma on a 7T BioSpec MRI Scanner: Correlation of Tumor Cellularity and Nuclear-to-Cytoplasmic Ratio with Apparent Diffusion Coefficient.

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

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

Quantitative Analysis of Permeability for Glioma Grading Using Dynamic Contrast-Enhanced Magnetic Resonance Imaging

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

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

Microstructure estimation from diffusion-MRI: Compartmentalized models in permeable cellular tissue

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

Impact of Arterial Input Function and Pharmacokinetic Models on DCE-MRI Biomarkers for Detection of Vascular Effect Induced by Stroma-Directed Drug in an Orthotopic Mouse Model of Pancreatic Cancer