Madhava P. Aryal

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

Pharmacokinetic analysis of dynamic contrast-enhanced (DCE) MRI data allows estimation of quantitative imaging biomarkers such as Ktrans (rate constant for plasma/interstitium contrast reagent (CR) transfer) and ve (extravascular and extracellular volume fraction). However, the use of quantitative DCE-MRI in clinical practice is limited with uncertainty in arterial input function (AIF) determination being one of the primary reasons. In this multicenter study to assess the effects of AIF variations on pharmacokinetic parameter estimation, DCE-MRI data acquired at one center from 11 prostate cancer patients were shared among nine centers. Individual AIF from each data set was determined by each center and submitted to the managing center. These AIFs, along with a literature population averaged AIF, and their reference-tissue-adjusted variants were used by the managing center to perform pharmacokinetic data analysis using the Tofts model (TM). All other variables, including tumor region of interest (ROI) definition and pre-contrast T1, were kept constant to evaluate parameter variations caused solely by AIF discrepancies. Considerable parameter variations were observed with the within-subject coefficient of variation (wCV) of Ktrans obtained with unadjusted AIFs being as high as 0.74. AIF-caused variations were larger in Ktrans than ve and both were reduced when reference-tissue-adjusted AIFs were used. These variations were largely systematic, resulting in nearly unchanged parametric map patterns. The intravasation rate constant, kep (= Ktrans/ve), was less sensitive to AIF variation than Ktrans (wCV for unadjusted AIFs: 0.45 vs. 0.74), suggesting that it might be a more robust imaging biomarker of prostate microvasculature than Ktrans.

A Radiation-Induced Hippocampal Vascular Injury Surrogate Marker Predicts Late Neurocognitive Dysfunction

PURPOSE We aimed to develop a hippocampal vascular injury surrogate marker for early prediction of late neurocognitive dysfunction in patients receiving brain radiation therapy (RT). METHODS AND MATERIALS Twenty-seven patients (17 males and 10 females, 31-80 years of age) were enrolled in an institutional review board-approved prospective longitudinal study. Patients received diagnoses of low-grade glioma or benign tumor and were treated by (3D) conformal or intensity-modulated RT with a median dose of 54 Gy (50.4-59.4 Gy in 1.8-Gy fractions). Six dynamic-contrast enhanced MRI scans were performed from pre-RT to 18-month post-RT, and quantified for vascular parameters related to blood-brain barrier permeability, K(trans), and the fraction of blood plasma volume, Vp. The temporal changes in the means of hippocampal transfer constant K(trans) and Vp after starting RT were modeled by integrating the dose effects with age, sex, hippocampal laterality, and presence of tumor or edema near a hippocampus. Finally, the early vascular dose response in hippocampi was correlated with neurocognitive dysfunction at 6 and 18 months post-RT. RESULTS The mean K(trans) Increased significantly from pre-RT to 1-month post-RT (P<.0004), which significantly depended on sex (P<.0007) and age (P<.00004), with the dose response more pronounced in older females. Also, the vascular dose response in the left hippocampus of females correlated significantly with changes in memory function at 6 (r=-0.95, P<.0006) and 18-months (r=-0.88, P<.02) post-RT. CONCLUSIONS The early hippocampal vascular dose response could be a predictor of late neurocognitive dysfunction. A personalized hippocampus sparing strategy may be considered in the future.

Acute Temporal Changes of MRI-Tracked Tumor Vascular Parameters after Combined Anti-angiogenic and Radiation Treatments in a Rat Glioma Model: Identifying Signatures of Synergism

In this study we used magnetic resonance imaging (MRI) biomarkers to monitor the acute temporal changes in tumor vascular physiology with the aim of identifying the vascular signatures that predict response to combined anti-angiogenic and radiation treatments. Forty-three athymic rats implanted with orthotopic U-251 glioma cells were studied for approximately 21 days after implantation. Two MRI studies were performed on each animal, pre- and post-treatment, to measure tumor vascular parameters. Two animal groups received treatment comprised of Cilengitide, an anti-angiogenic agent and radiation. The first group received a subcurative regimen of Cilengitide 1 h before irradiation, while the second group received a curative regimen of Cilengitide 8 h before irradiation. Cilengitide was given as a single dose (4 mg/kg; intraperitoneal) after the pretreatment MRI study and before receiving a 20 Gy radiation dose. After irradiation, the post-treatment MRI study was performed at selected time points: 2, 4, 8 and 12 h (n = ≥5 per time point). Significant changes in vascular parameters were observed at early time points after combined treatments in both treatment groups (1 and 8 h). The temporal changes in vascular parameters in the first group (treated 1 h before exposure) resembled a previously reported pattern associated with radiation exposure alone. Conversely, in the second group (treated 8 h before exposure), all vascular parameters showed an initial response at 2–4 h postirradiation, followed by an apparent lack of response at later time points. The signature time point to define the “synergy” of Cilengitide and radiation was 4 h postirradiation. For example, 4 h after combined treatments using a 1 h separation (which followed the subcurative regimen), tumor blood flow was significantly decreased, nearly 50% below baseline (P = 0.007), whereas 4 h after combined treatments using an 8 h separation (which followed the curative regimen), tumor blood flow was only 10% less than baseline. Comparison between the first and second groups further revealed that most other vascular parameters were maximally different 4 h after combined treatments. In conclusion, the data are consistent with the assertion that the delivery of radiation at the vascular normalization time window of Cilengitide improves radiation treatment outcome. The different vascular responses after the different delivery times of combined treatments in light of the known tumor responses under similar conditions would indicate that timing has a crucial influence on treatment outcome and long-term survival. Tracking acute changes in tumor physiology after monotherapy or combined treatments appears to aid in identifying the beneficial timing for administration, and perhaps has predictive value. Therefore, judicial timing of treatments may result in optimal treatment response.

Accuracy, repeatability, and interplatform reproducibility of T1 quantification methods used for DCE‐MRI: Results from a multicenter phantom study

To determine the in vitro accuracy, test‐retest repeatability, and interplatform reproducibility of T1 quantification protocols used for dynamic contrast‐enhanced MRI at 1.5 and 3 T.

Peritumoral tissue compression is predictive of exudate flux in a rat model of cerebral tumor: An MRI study in an embedded tumor

MRI estimates of extracellular volume and tumor exudate flux in peritumoral tissue are demonstrated in an experimental model of cerebral tumor. Peritumoral extracellular volume predicted the tumor exudate flux.

Multisite Concordance of DSC-MRI Analysis for Brain Tumors: Results of a National Cancer Institute Quantitative Imaging Network Collaborative Project

DSC-MR imaging data were collected after a preload and during a bolus injection of gadolinium contrast agent using a gradient recalled-echo-EPI sequence. Forty-nine low-grade and high-grade glioma datasets were uploaded to The Cancer Imaging Archive. Datasets included a predetermined arterial input function, enhancing tumor ROIs, and ROIs necessary to create normalized relative CBV and CBF maps. Seven sites computed 20 different perfusion metrics. For normalized relative CBV and normalized CBF, 93% and 94% of entries showed good or excellent cross-site agreement. All metrics could distinguish low- from high-grade tumors. BACKGROUND AND PURPOSE: Standard assessment criteria for brain tumors that only include anatomic imaging continue to be insufficient. While numerous studies have demonstrated the value of DSC-MR imaging perfusion metrics for this purpose, they have not been incorporated due to a lack of confidence in the consistency of DSC-MR imaging metrics across sites and platforms. This study addresses this limitation with a comparison of multisite/multiplatform analyses of shared DSC-MR imaging datasets of patients with brain tumors. MATERIALS AND METHODS: DSC-MR imaging data were collected after a preload and during a bolus injection of gadolinium contrast agent using a gradient recalled-echo–EPI sequence (TE/TR = 30/1200 ms; flip angle = 72°). Forty-nine low-grade (n = 13) and high-grade (n = 36) glioma datasets were uploaded to The Cancer Imaging Archive. Datasets included a predetermined arterial input function, enhancing tumor ROIs, and ROIs necessary to create normalized relative CBV and CBF maps. Seven sites computed 20 different perfusion metrics. Pair-wise agreement among sites was assessed with the Lin concordance correlation coefficient. Distinction of low- from high-grade tumors was evaluated with the Wilcoxon rank sum test followed by receiver operating characteristic analysis to identify the optimal thresholds based on sensitivity and specificity. RESULTS: For normalized relative CBV and normalized CBF, 93% and 94% of entries showed good or excellent cross-site agreement (0.8 ≤ Lin concordance correlation coefficient ≤ 1.0). All metrics could distinguish low- from high-grade tumors. Optimum thresholds were determined for pooled data (normalized relative CBV = 1.4, sensitivity/specificity = 90%:77%; normalized CBF = 1.58, sensitivity/specificity = 86%:77%). CONCLUSIONS: By means of DSC-MR imaging data obtained after a preload of contrast agent, substantial consistency resulted across sites for brain tumor perfusion metrics with a common threshold discoverable for distinguishing low- from high-grade tumors.

Multisite concordance of apparent diffusion coefficient measurements across the NCI quantitative imaging network

Abstract. Diffusion weighted MRI has become ubiquitous in many areas of medicine, including cancer diagnosis and treatment response monitoring. Reproducibility of diffusion metrics is essential for their acceptance as quantitative biomarkers in these areas. We examined the variability in the apparent diffusion coefficient (ADC) obtained from both postprocessing software implementations utilized by the NCI Quantitative Imaging Network and online scan time-generated ADC maps. Phantom and in vivo breast studies were evaluated for two (ADC2) and four (ADC4) b-value diffusion metrics. Concordance of the majority of implementations was excellent for both phantom ADC measures and in vivo ADC2, with relative biases <0.1% (ADC2) and <0.5% (phantom ADC4) but with higher deviations in ADC at the lowest phantom ADC values. In vivo ADC4 concordance was good, with typical biases of ±2% to 3% but higher for online maps. Multiple b-value ADC implementations were separated into two groups determined by the fitting algorithm. Intergroup mean ADC differences ranged from negligible for phantom data to 2.8% for ADC4 in vivo data. Some higher deviations were found for individual implementations and online parametric maps. Despite generally good concordance, implementation biases in ADC measures are sometimes significant and may be large enough to be of concern in multisite studies.

Impact of uncertainty in longitudinal T1 measurements on quantification of dynamic contrast‐enhanced MRI

The objective of this study was to assess the uncertainty in T1 measurement, by estimating the repeatability coefficient (RC) from two repeated scans, in normal appearing brain tissues employing two different T1 mapping methods. All brain MRI scans were performed on a 3 T MR scanner in 10 patients who had low grade/benign tumors and partial brain radiation therapy (RT) without chemotherapy, at pre‐RT, 3 weeks into RT, end RT (6 weeks) and 11, 33, and 85 weeks after RT. T1‐weighted images were acquired using (1) a spoiled gradient echo sequence with two flip angles (2FA: 5° and 15°) and (2) a progressive saturation recovery sequence (pSR) with five different TR values (100–2000 ms). Manually drawn volumes of interest (VOIs) included left and right normal putamen and thalamus in gray matter, and frontal and parietal white matter, which were distant from tumors and received a total of accumulated radiation doses less than 5 Gy at 3 weeks. No significant changes or even trends in mean T1 from pre‐RT to 3 weeks into RT in these VOIs (p ≥ 0.11, Wilcoxon sign test) allowed us to calculate the repeatability statistics of between‐subject means of squares, within‐subject means of squares, F‐score, and RC. The 2FA method produced RCs in the range of (9.7–11.7)% in gray matter and (12.2–14.5)% in white matter; while the pSR method led to RCs ranging from 10.9 to 17.9% in gray matter and 7.5 to 10.3% in white matter. The overall mean (±SD) RCs produced by the two methods, 12.0 (±1.6)% for 2FA and 12.0 (±3.8)% for pSR, were not significantly different (p = 0.97). A similar repeatability in T1 measurement produced by the time efficient 2FA method compared with the time consuming pSR method demonstrates that the 2FA method is desirable to integrate into dynamic contrast‐enhanced MRI for rapid acquisition. Copyright

Tumor image signatures and habitats: a processing pipeline of multimodality metabolic and physiological images

Abstract. To create tumor “habitats” from the “signatures” discovered from multimodality metabolic and physiological images, we developed a framework of a processing pipeline. The processing pipeline consists of six major steps: (1) creating superpixels as a spatial unit in a tumor volume; (2) forming a data matrix {D} containing all multimodality image parameters at superpixels; (3) forming and clustering a covariance or correlation matrix {C} of the image parameters to discover major image “signatures;” (4) clustering the superpixels and organizing the parameter order of the {D} matrix according to the one found in step 3; (5) creating “habitats” in the image space from the superpixels associated with the “signatures;” and (6) pooling and clustering a matrix consisting of correlation coefficients of each pair of image parameters from all patients to discover subgroup patterns of the tumors. The pipeline was applied to a dataset of multimodality images in glioblastoma (GBM) first, which consisted of 10 image parameters. Three major image “signatures” were identified. The three major “habitats” plus their overlaps were created. To test generalizability of the processing pipeline, a second image dataset from GBM, acquired on the scanners different from the first one, was processed. Also, to demonstrate the clinical association of image-defined “signatures” and “habitats,” the patterns of recurrence of the patients were analyzed together with image parameters acquired prechemoradiation therapy. An association of the recurrence patterns with image-defined “signatures” and “habitats” was revealed. These image-defined “signatures” and “habitats” can be used to guide stereotactic tissue biopsy for genetic and mutation status analysis and to analyze for prediction of treatment outcomes, e.g., patterns of failure.

Temporal Feature Extraction from DCE-MRI to Identify Poorly Perfused Subvolumes of Tumors Related to Outcomes of Radiation Therapy in Head and Neck Cancer

This study aimed to develop an automated model to extract temporal features from DCE-MRI in head-and-neck (HN) cancers to localize significant tumor subvolumes having low blood volume (LBV) for predicting local and regional failure after chemoradiation therapy. Temporal features were extracted from time-intensity curves to build classification model for differentiating voxels with LBV from those with high BV. Support vector machine (SVM) classification was trained on the extracted features for voxel classification. Subvolumes with LBV were then assembled from the classified voxels with LBV. The model was trained and validated on independent datasets created from 456 873 DCE curves. The resultant subvolumes were compared to ones derived by a 2-step method via pharmacokinetic modeling of blood volume, and evaluated for classification accuracy and volumetric similarity by DSC. The proposed model achieved an average voxel-level classification accuracy and DSC of 82% and 0.72, respectively. Also, the model showed tolerance on different acquisition parameters of DCE-MRI. The model could be directly used for outcome prediction and therapy assessment in radiation therapy of HN cancers, or even supporting boost target definition in adaptive clinical trials with further validation. The model is fully automatable, extendable, and scalable to extract temporal features of DCE-MRI in other tumors.