The Impact of Arterial Input Function Determination Variations on Prostate Dynamic Contrast-Enhanced Magnetic Resonance Imaging Pharmacokinetic Modeling: A Multicenter Data Analysis Challenge

Dynamic contrast-enhanced magnetic resonance imaging for head and neck cancers

Dynamic myraidpro contrast-enhanced magnetic resonance imaging (DCE-MRI) has been correlated with prognosis in head and neck squamous cell carcinoma as well as with changes in normal tissues. These studies implement different software, either commercial or in-house, and different scan protocols. Thus, the generalizability of the results is not confirmed. To assist in the standardization of quantitative metrics to confirm the generalizability of these previous studies, this data descriptor delineates in detail the DCE-MRI digital imaging and communications in medicine (DICOM) files with DICOM radiation therapy (RT) structure sets and digital reference objects (DROs), as well as, relevant clinical data that encompass a data set that can be used by all software for comparing quantitative metrics. Variable flip angle (VFA) with six flip angles and DCE-MRI scans with a temporal resolution of 5.5 s were acquired in the axial direction on a 3T MR scanner with a field of view of 25.6 cm, slice thickness of 4 mm, and 256×256 matrix size.

Robust and efficient pharmacokinetic parameter non-linear least squares estimation for dynamic contrast enhanced MRI of the prostate

PURPOSE

To describe an efficient numerical optimization technique using non-linear least squares to estimate perfusion parameters for the Tofts and extended Tofts models from dynamic contrast enhanced (DCE) MRI data and apply the technique to prostate cancer.

METHODS

Parameters were estimated by fitting the two Tofts-based perfusion models to the acquired data via non-linear least squares. We apply Variable Projection (VP) to convert the fitting problem from a multi-dimensional to a one-dimensional line search to improve computational efficiency and robustness. Using simulation and DCE-MRI studies in twenty patients with suspected prostate cancer, the VP-based solver was compared against the traditional Levenberg-Marquardt (LM) strategy for accuracy, noise amplification, robustness to converge, and computation time.

RESULTS

The simulation demonstrated that VP and LM were both accurate in that the medians closely matched assumed values across typical signal to noise ratio (SNR) levels for both Tofts models. VP and LM showed similar noise sensitivity. Studies using the patient data showed that the VP method reliably converged and matched results from LM with approximate 3× and 2× reductions in computation time for the standard (two-parameter) and extended (three-parameter) Tofts models. While LM failed to converge in 14% of the patient data, VP converged in the ideal 100%.

CONCLUSION

The VP-based method for non-linear least squares estimation of perfusion parameters for prostate MRI is equivalent in accuracy and robustness to noise, while being more reliably (100%) convergent and computationally about 3× (TM) and 2× (ETM) faster than the LM-based method.

Assessment of a statistical AIF extraction method for dynamic PET studies with 15O water and 18F fluorodeoxyglucose in locally advanced breast cancer patients

Blood flow-metabolism mismatch from dynamic positron emission tomography (PET) studies with O15-labeled water (H2O) and F18-labeled fluorodeoxyglucose (FDG) has been shown to be a promising diagnostic for locally advanced breast cancer (LABCa) patients. The mismatch measurement involves kinetic analysis with the arterial blood time course (AIF) as an input function. We evaluate the use of a statistical method for AIF extraction (SAIF) in these studies. Fifty three LABCa patients had dynamic PET studies with H2O and FDG. For each PET study, two AIFs were recovered, an SAIF extraction and also a manual extraction based on a region of interest placed over the left ventricle (LV-ROI). Blood flow-metabolism mismatch was obtained with each AIF, and kinetic and prognostic reliability comparisons were made. Strong correlations were found between kinetic assessments produced by both AIFs. SAIF AIFs retained the full prognostic value, for pathologic response and overall survival, of LV-ROI AIFs.

Linearization improves the repeatability of quantitative dynamic contrast-enhanced MRI

PURPOSE

The purpose of this study was to compare the repeatabilities of the linear and nonlinear Tofts and reference region models (RRM) for dynamic contrast-enhanced MRI (DCE-MRI).

MATERIALS AND METHODS

Simulated and experimental DCE-MRI data from 12 rats with a flank tumor of C6 glioma acquired over three consecutive days were analyzed using four quantitative and semi-quantitative DCE-MRI metrics. The quantitative methods used were: 1) linear Tofts model (LTM), 2) non-linear Tofts model (NTM), 3) linear RRM (LRRM), and 4) non-linear RRM (NRRM). The following semi-quantitative metrics were used: 1) maximum enhancement ratio (MER), 2) time to peak (TTP), 3) initial area under the curve (iauc64), and 4) slope. LTM and NTM were used to estimate K^trans, while LRRM and NRRM were used to estimate K^trans relative to muscle (R^Ktrans). Repeatability was assessed by calculating the within-subject coefficient of variation (wSCV) and the percent intra-subject variation (iSV) determined with the Gage R&R analysis.

RESULTS

The iSV for R^Ktrans using LRRM was two-fold lower compared to NRRM at all simulated and experimental conditions. A similar trend was observed for the Tofts model, where LTM was at least 50% more repeatable than the NTM under all experimental and simulated conditions. The semi-quantitative metrics iauc64 and MER were as equally repeatable as K^trans and R^Ktrans estimated by LTM and LRRM respectively. The iSV for iauc64 and MER were significantly lower than the iSV for slope and TTP.

CONCLUSION

In simulations and experimental results, linearization improves the repeatability of quantitative DCE-MRI by at least 30%, making it as repeatable as semi-quantitative metrics.

Automatic Detection and Quantitative DCE-MRI Scoring of Prostate Cancer Aggressiveness

Purpose

To develop a robust and clinically applicable automated method for analyzing Dynamic Contrast Enhanced (DCE-) MRI of the prostate as a guide for targeted biopsies and treatments.

Materials and methods

An unsupervised pattern recognition (PR) method was used to analyze prostate DCE-MRI from 71 sequential radiotherapy patients. Identified regions of interest (ROIs) with increased perfusion were assigned either to the peripheral (PZ) or transition zone (TZ). Six quantitative features, associated with the washin and washout part of the weighted average DCE curve from the ROI, were calculated. The associations between the assigned DCE-scores and Gleason Score (GS) were investigated. A heatmap of tumor aggressiveness covering the entire prostate was generated and validated with histopathology from MRI-ultrasound fused (MRI-US) targeted biopsies.

Results

The volumes of the PR-identified ROI’s were significantly correlated with the highest GS from the biopsy session for each patient. Following normalization (and only after normalization) with gluteus maximus muscle’s DCE signal, the quantitative features in PZ were significantly correlated with GS. These correlations straightened in subset of patients with available MRI-US biopsies when GS from the individual biopsies were used. Area under the receiver operating characteristics curve for discrimination between indolent vs aggressive cancer for the significant quantitative features reached 0.88–0.95. When DCE-scores were calculated in normal appearing tissues, the features were highly discriminative for cancer vs no cancer both in PZ and TZ. The generated heatmap of tumor aggressiveness coincided with the location and GS of the MRI-US biopsies.

Conclusion

A quantitative approach for DCE-MRI analysis was developed. The resultant map of aggressiveness correlated well with tumor location and GS and is applicable for integration in radiotherapy/radiology imaging software for clinical translation.

A Multi-Institutional Comparison of Dynamic Contrast-Enhanced Magnetic Resonance Imaging Parameter Calculations

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides quantitative metrics (e.g. K^trans, v_e) via pharmacokinetic models. We tested inter-algorithm variability in these quantitative metrics with 11 published DCE-MRI algorithms, all implementing Tofts-Kermode or extended Tofts pharmacokinetic models. Digital reference objects (DROs) with known K^trans and v_e values were used to assess performance at varying noise levels. Additionally, DCE-MRI data from 15 head and neck squamous cell carcinoma patients over 3 time-points during chemoradiotherapy were used to ascertain K^trans and v_e kinetic trends across algorithms. Algorithms performed well (less than 3% average error) when no noise was present in the DRO. With noise, 87% of K^trans and 84% of v_e algorithm-DRO combinations were generally in the correct order. Low Krippendorff’s alpha values showed that algorithms could not consistently classify patients as above or below the median for a given algorithm at each time point or for differences in values between time points. A majority of the algorithms produced a significant Spearman correlation in v_e of the primary gross tumor volume with time. Algorithmic differences in K^trans and v_e values over time indicate limitations in combining/comparing data from distinct DCE-MRI model implementations. Careful cross-algorithm quality-assurance must be utilized as DCE-MRI results may not be interpretable using differing software.

Effects of arterial input function selection on kinetic parameters in brain dynamic contrast-enhanced MRI

PURPOSE

Kinetic parameters derived from dynamic contrast-enhanced MRI (DCE-MRI) were suggested as a possible instrument for multi-parametric lesion characterization, but have not found their way into clinical practice yet due to inconsistent results. The quantification is heavily influenced by the definition of an appropriate arterial input functions (AIF). Regarding brain tumor DCE-MRI, there are currently several co-existing methods to determine the AIF frequently including different brain vessels as sources. This study quantitatively and qualitatively analyzes the impact of AIF source selection on kinetic parameters derived from commonly selected AIF source vessels compared to a population-based AIF model.

MATERIAL AND METHODS

74 patients with brain lesions underwent 3D DCE-MRI. Kinetic parameters [transfer constants of contrast agent efflux and reflux K^trans and k_ep and, their ratio, v_e, that is used to measure extravascular-extracellular volume fraction and plasma volume fraction v_p] were determined using extended Tofts model in 821 ROI from 4 AIF sources [the internal carotid artery (ICA), the closest artery to the lesion, the superior sagittal sinus (SSS), the population-based Parker model]. The effect of AIF source alteration on kinetic parameters was evaluated by tissue type selective intra-class correlation (ICC) and capacity to differentiate gliomas by WHO grade [area under the curve analysis (AUC)].

RESULTS

Arterial AIF more often led to implausible v_e >100% values (p<0.0001). AIF source alteration rendered different absolute kinetic parameters (p<0.0001), except for k_ep. ICC between kinetic parameters of different AIF sources and tissues were variable (0.08-0.87) and only consistent >0.5 between arterial AIF derived kinetic parameters. Differentiation between WHO III and II glioma was exclusively possible with v_p derived from an AIF in the SSS (p=0.03; AUC 0.74).

CONCLUSION

The AIF source has a significant impact on absolute kinetic parameters in DCE-MRI, which limits the comparability of kinetic parameters derived from different AIF sources. The effect is also tissue-dependent. The SSS appears to be the best choice for AIF source vessel selection in brain tumor DCE-MRI as it exclusively allowed for WHO grades II/III and III/IV glioma distinction (by v_p) and showed the least number of implausible v_e values.

DCE-MRI of the prostate using shutter-speed vs. Tofts model for tumor characterization and assessment of aggressiveness

PURPOSE

To quantify Tofts model (TM) and shutter-speed model (SSM) perfusion parameters in prostate cancer (PCa) and noncancerous peripheral zone (PZ) and to compare the diagnostic performance of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) to Prostate Imaging and Reporting and Data System (PI-RADS) classification for the assessment of PCa aggressiveness.

MATERIALS AND METHODS

Fifty PCa patients (mean age 60 years old) who underwent MRI at 3.0T followed by prostatectomy were included in this Institutional Review Board-approved retrospective study. DCE-MRI parameters (K^trans , v_e , k_ep [TM&SSM] and intracellular water molecule lifetime τ_i [SSM]) were determined in PCa and PZ. Differences in DCE-MRI parameters between PCa and PZ, and between models were assessed using Wilcoxon signed-rank tests. Receiver operating characteristic (ROC) analysis for differentiation between PCa and PZ was performed for individual and combined DCE-MRI parameters. Diagnostic performance of DCE-MRI parameters for identification of aggressive PCa (Gleason ≥8, grade group [GG] ≥3 or pathology stage pT3) was assessed using ROC analysis and compared with PI-RADSv2 scores.

RESULTS

DCE-MRI parameters were significantly different between TM and SSM and between PZ and PCa (P < 0.037). Diagnostic performances of TM and SSM for differentiation of PCa from PZ were similar (highest AUC TM: K^trans +k_ep 0.76, SSM: τ_i +k_ep 0.80). PI-RADS outperformed TM and SSM DCE-MRI for identification of Gleason ≥8 lesions (AUC PI-RADS: 0.91, highest AUC DCE-MRI: K^trans +τ_i SSM 0.61, P = 0.002). The diagnostic performance of PI-RADS and DCE-MRI for identification of GG ≥3 and pT3 PCa was not significantly different (P > 0.213).

CONCLUSION

SSM DCE-MRI did not increase the diagnostic performance of DCE-MRI for PCa characterization. PI-RADS outperformed both TM and SSM DCE-MRI for identification of aggressive cancer.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 2 J. MAGN. RESON. IMAGING 2017;46:837-849.

Evaluation of Soft Tissue Sarcoma Response to Preoperative Chemoradiotherapy Using Dynamic Contrast-Enhanced Magnetic Resonance Imaging

This study aims to assess the utility of quantitative dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) parameters in comparison with imaging tumor size for early prediction and evaluation of soft tissue sarcoma response to preoperative chemoradiotherapy. In total, 20 patients with intermediate- to high-grade soft tissue sarcomas received either a phase I trial regimen of sorafenib + chemoradiotherapy (n = 8) or chemoradiotherapy only (n = 12), and underwent DCE-MRI at baseline, after 2 weeks of treatment with sorafenib or after the first chemotherapy cycle, and after therapy completion. MRI tumor size in the longest diameter (LD) was measured according to the RECIST (Response Evaluation Criteria In Solid Tumors) guidelines. Pharmacokinetic analyses of DCE-MRI data were performed using the Shutter-Speed model. After only 2 weeks of treatment with sorafenib or after 1 chemotherapy cycle, K^trans (rate constant for plasma/interstitium contrast agent transfer) and its percent change were good early predictors of optimal versus suboptimal pathological response with univariate logistic regression C statistics values of 0.90 and 0.80, respectively, whereas RECIST LD percent change was only a fair predictor (C = 0.72). Post-therapy K^trans, v_e (extravascular and extracellular volume fraction), and k_ep (intravasation rate constant), not RECIST LD, were excellent (C > 0.90) markers of therapy response. Several DCE-MRI parameters before, during, and after therapy showed significant (P < .05) correlations with percent necrosis of resected tumor specimens. In conclusion, absolute values and percent changes of quantitative DCE-MRI parameters provide better early prediction and evaluation of the pathological response of soft tissue sarcoma to preoperative chemoradiotherapy than the conventional measurement of imaging tumor size change.

Computational Challenges and Collaborative Projects in the NCI Quantitative Imaging Network

The Quantitative Imaging Network (QIN) of the National Cancer Institute (NCI) conducts research in development and validation of imaging tools and methods for predicting and evaluating clinical response to cancer therapy. Members of the network are involved in examining various imaging and image assessment parameters through network-wide cooperative projects. To more effectively use the cooperative power of the network in conducting computational challenges in benchmarking of tools and methods and collaborative projects in analytical assessment of imaging technologies, the QIN Challenge Task Force has developed policies and procedures to enhance the value of these activities by developing guidelines and leveraging NCI resources to help their administration and manage dissemination of results. Challenges and Collaborative Projects (CCPs) are further divided into technical and clinical CCPs. As the first NCI network to engage in CCPs, we anticipate a variety of CCPs to be conducted by QIN teams in the coming years. These will be aimed to benchmark advanced software tools for clinical decision support, explore new imaging biomarkers for therapeutic assessment, and establish consensus on a range of methods and protocols in support of the use of quantitative imaging to predict and assess response to cancer therapy.

Increased robustness in reference region model analysis of DCE MRI using two-step constrained approaches

PURPOSE

Reference region models (RRMs) can quantify tumor perfusion in dynamic contrast-enhanced MRI without an arterial input function. Inspection of the RRM reveals that one of the free parameters in the fit is uniquely linked to the reference region and is common to all voxels. A two-step approach is proposed that takes this constraint into account.

METHODS

Three constrained RRM (CRRM) approaches were devised and evaluated. Simulations were performed to compare their accuracy and precision over a range of noise and temporal resolutions. The CRRM was also applied on a virtual phantom that simulates different perfusion values. In vivo evaluation was performed on data from breast cancer and soft tissue sarcoma.

RESULTS

In simulations, the CRRM consistently improved precision and had better accuracy at low signal-to-noise ratio (SNR). In virtual phantom, the CRRMs were able to fit voxels that had similar kinetics to the reference tissue, whereas the unconstrained models failed to accurately fit these voxels. In the in vivo data, the constrained approaches produced parameter maps that had less variability and were in better agreement with the Tofts model.

CONCLUSION

These findings indicate that the two-step fitting approach of the CRRM can reduce the variability of perfusion estimates for quantifying perfusion with dynamic contrast-enhanced (DCE) MRI. Magn Reson Med 78:1547-1557, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

DCE-MRI Perfusion and Permeability Parameters as predictors of tumor response to CCRT in Patients with locally advanced NSCLC

In this prospective study, 36 patients with stage III non-small cell lung cancers (NSCLC), who underwent dynamic contrast-enhanced MRI (DCE-MRI) before concurrent chemo-radiotherapy (CCRT) were enrolled. Pharmacokinetic analysis was carried out after non-rigid motion registration. The perfusion parameters [including Blood Flow (BF), Blood Volume (BV), Mean Transit Time (MTT)] and permeability parameters [including endothelial transfer constant (K^trans), reflux rate (K_ep), fractional extravascular extracellular space volume (V_e), fractional plasma volume (V_p)] were calculated, and their relationship with tumor regression was evaluated. The value of these parameters on predicting responders were calculated by receiver operating characteristic (ROC) curve. Multivariate logistic regression analysis was conducted to find the independent variables. Tumor regression rate is negatively correlated with V_e and its standard variation V_e_SD and positively correlated with K^trans and K_ep. Significant differences between responders and non-responders existed in K^trans, K_ep, V_e, V_e_SD, MTT, BV_SD and MTT_SD (P < 0.05). ROC indicated that V_e < 0.24 gave the largest area under curve of 0.865 to predict responders. Multivariate logistic regression analysis also showed V_e was a significant predictor. Baseline perfusion and permeability parameters calculated from DCE-MRI were seen to be a viable tool for predicting the early treatment response after CCRT of NSCLC.

Relative sensitivities of DCE-MRI pharmacokinetic parameters to arterial input function (AIF) scaling

Dynamic-Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) has been used widely for clinical applications. Pharmacokinetic modeling of DCE-MRI data that extracts quantitative contrast reagent/tissue-specific model parameters is the most investigated method. One of the primary challenges in pharmacokinetic analysis of DCE-MRI data is accurate and reliable measurement of the arterial input function (AIF), which is the driving force behind all pharmacokinetics. Because of effects such as inflow and partial volume averaging, AIF measured from individual arteries sometimes require amplitude scaling for better representation of the blood contrast reagent (CR) concentration time-courses. Empirical approaches like blinded AIF estimation or reference tissue AIF derivation can be useful and practical, especially when there is no clearly visible blood vessel within the imaging field-of-view (FOV). Similarly, these approaches generally also require magnitude scaling of the derived AIF time-courses. Since the AIF varies among individuals even with the same CR injection protocol and the perfect scaling factor for reconstructing the ground truth AIF often remains unknown, variations in estimated pharmacokinetic parameters due to varying AIF scaling factors are of special interest. In this work, using simulated and real prostate cancer DCE-MRI data, we examined parameter variations associated with AIF scaling. Our results show that, for both the fast-exchange-limit (FXL) Tofts model and the water exchange sensitized fast-exchange-regime (FXR) model, the commonly fitted CR transfer constant (K(trans)) and the extravascular, extracellular volume fraction (ve) scale nearly proportionally with the AIF, whereas the FXR-specific unidirectional cellular water efflux rate constant, kio, and the CR intravasation rate constant, kep, are both AIF scaling insensitive. This indicates that, for DCE-MRI of prostate cancer and possibly other cancers, kio and kep may be more suitable imaging biomarkers for cross-platform, multicenter applications. Data from our limited study cohort show that kio correlates with Gleason scores, suggesting that it may be a useful biomarker for prostate cancer disease progression monitoring.