Imaging intratumor heterogeneity: role in therapy response, resistance, and clinical outcome

Intratumor partitioning and texture analysis of dynamic contrast-enhanced (DCE)-MRI identifies relevant tumor subregions to predict pathological response of breast cancer to neoadjuvant chemotherapy

PURPOSE

To predict pathological response of breast cancer to neoadjuvant chemotherapy (NAC) based on quantitative, multiregion analysis of dynamic contrast enhancement magnetic resonance imaging (DCE-MRI).

MATERIALS AND METHODS

In this Institutional Review Board-approved study, 35 patients diagnosed with stage II/III breast cancer were retrospectively investigated using 3T DCE-MR images acquired before and after the first cycle of NAC. First, principal component analysis (PCA) was used to reduce the dimensionality of the DCE-MRI data with high temporal resolution. We then partitioned the whole tumor into multiple subregions using k-means clustering based on the PCA-defined eigenmaps. Within each tumor subregion, we extracted four quantitative Haralick texture features based on the gray-level co-occurrence matrix (GLCM). The change in texture features in each tumor subregion between pre- and during-NAC was used to predict pathological complete response after NAC.

RESULTS

Three tumor subregions were identified through clustering, each with distinct enhancement characteristics. In univariate analysis, all imaging predictors except one extracted from the tumor subregion associated with fast washout were statistically significant (P < 0.05) after correcting for multiple testing, with area under the receiver operating characteristic (ROC) curve (AUC) or AUCs between 0.75 and 0.80. In multivariate analysis, the proposed imaging predictors achieved an AUC of 0.79 (P = 0.002) in leave-one-out cross-validation. This improved upon conventional imaging predictors such as tumor volume (AUC = 0.53) and texture features based on whole-tumor analysis (AUC = 0.65).

CONCLUSION

The heterogeneity of the tumor subregion associated with fast washout on DCE-MRI predicted pathological response to NAC in breast cancer. J. Magn. Reson. Imaging 2016;44:1107-1115.

Hypoxia-activated prodrugs: paths forward in the era of personalised medicine

Tumour hypoxia has been pursued as a cancer drug target for over 30 years, most notably using bioreductive (hypoxia-activated) prodrugs that target antineoplastic agents to low-oxygen tumour compartments. Despite compelling evidence linking hypoxia with treatment resistance and adverse prognosis, a number of such prodrugs have recently failed to demonstrate efficacy in pivotal clinical trials; an outcome that demands reflection on the discovery and development of these compounds. In this review, we discuss a clear disconnect between the pathobiology of tumour hypoxia, the pharmacology of hypoxia-activated prodrugs and the manner in which they have been taken into clinical development. Hypoxia-activated prodrugs have been evaluated in the manner of broad-spectrum cytotoxic agents, yet a growing body of evidence suggests that their activity is likely to be dependent on the coincidence of tumour hypoxia, expression of specific prodrug-activating reductases and intrinsic sensitivity of malignant clones to the cytotoxic effector. Hypoxia itself is highly variable between and within individual tumours and is not treatment-limiting in all cancer subtypes. Defining predictive biomarkers for hypoxia-activated prodrugs and overcoming the technical challenges of assaying them in clinical settings will be essential to deploying these agents in the era of personalised cancer medicine.

Decline of Tumor Vascular Function as Assessed by Dynamic Contrast-Enhanced Magnetic Resonance Imaging Is Associated With Poor Responses to Radiation Therapy and Chemotherapy

PURPOSE

To investigate whether changes in the volume transfer coefficient (K(trans)) in a growing tumor could be used as a surrogate marker for predicting tumor responses to radiation therapy (RT) and chemotherapy (CT).

METHODS AND MATERIALS

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) was consecutively performed on tumor-bearing mice, and temporal and spatial changes of K(trans) values were measured along with tumor growth. Tumor responses to RT and CT were studied before and after observed changes in K(trans) values with time.

RESULTS

Dynamic changes with an initial increase and subsequent decline in K(trans) values were found to be associated with tumor growth. When each tumor was divided into core and peripheral regions, the K(trans) decline was greater in core, although neither vascular structure or necrosis could be linked to this spatial difference. Tumor responses to RT were worse if applied after the decline of K(trans), and there was less drug distribution and cell death in the tumor core after CT.

CONCLUSION

The K(trans) value in growing tumors, reflecting the changes of tumor microenvironment and vascular function, is strongly associated with tumor responses to RT and CT and could be a potential surrogate marker for predicting the tumor response to these treatments.

Robust Intratumor Partitioning to Identify High-Risk Subregions in Lung Cancer: A Pilot Study

PURPOSE

To develop an intratumor partitioning framework for identifying high-risk subregions from (18)F-fluorodeoxyglucose positron emission tomography (FDG-PET) and computed tomography (CT) imaging and to test whether tumor burden associated with the high-risk subregions is prognostic of outcomes in lung cancer.

METHODS AND MATERIALS

In this institutional review board-approved retrospective study, we analyzed the pretreatment FDG-PET and CT scans of 44 lung cancer patients treated with radiation therapy. A novel, intratumor partitioning method was developed, based on a 2-stage clustering process: first at the patient level, each tumor was over-segmented into many superpixels by k-means clustering of integrated PET and CT images; next, tumor subregions were identified by merging previously defined superpixels via population-level hierarchical clustering. The volume associated with each of the subregions was evaluated using Kaplan-Meier analysis regarding its prognostic capability in predicting overall survival (OS) and out-of-field progression (OFP).

RESULTS

Three spatially distinct subregions were identified within each tumor that were highly robust to uncertainty in PET/CT co-registration. Among these, the volume of the most metabolically active and metabolically heterogeneous solid component of the tumor was predictive of OS and OFP on the entire cohort, with a concordance index or CI of 0.66-0.67. When restricting the analysis to patients with stage III disease (n=32), the same subregion achieved an even higher CI of 0.75 (hazard ratio 3.93, log-rank P=.002) for predicting OS, and a CI of 0.76 (hazard ratio 4.84, log-rank P=.002) for predicting OFP. In comparison, conventional imaging markers, including tumor volume, maximum standardized uptake value, and metabolic tumor volume using threshold of 50% standardized uptake value maximum, were not predictive of OS or OFP, with CI mostly below 0.60 (log-rank P>.05).

CONCLUSION

We propose a robust intratumor partitioning method to identify clinically relevant, high-risk subregions in lung cancer. We envision that this approach will be applicable to identifying useful imaging biomarkers in many cancer types.

Development of a nomogram combining clinical staging with (18)F-FDG PET/CT image features in non-small-cell lung cancer stage I-III

PURPOSE

Our goal was to develop a nomogram by exploiting intratumour heterogeneity on CT and PET images from routine (18)F-FDG PET/CT acquisitions to identify patients with the poorest prognosis.

METHODS

This retrospective study included 116 patients with NSCLC stage I, II or III and with staging (18)F-FDG PET/CT imaging. Primary tumour volumes were delineated using the FLAB algorithm and 3D Slicer™ on PET and CT images, respectively. PET and CT heterogeneities were quantified using texture analysis. The reproducibility of the CT features was assessed on a separate test-retest dataset. The stratification power of the PET/CT features was evaluated using the Kaplan-Meier method and the log-rank test. The best standard metric (functional volume) was combined with the least redundant and most prognostic PET/CT heterogeneity features to build the nomogram.

RESULTS

PET entropy and CT zone percentage had the highest complementary values with clinical stage and functional volume. The nomogram improved stratification amongst patients with stage II and III disease, allowing identification of patients with the poorest prognosis (clinical stage III, large tumour volume, high PET heterogeneity and low CT heterogeneity).

CONCLUSION

Intratumour heterogeneity quantified using textural features on both CT and PET images from routine staging (18)F-FDG PET/CT acquisitions can be used to create a nomogram with higher stratification power than staging alone.

Oxygen-Enhanced MRI Accurately Identifies, Quantifies, and Maps Tumor Hypoxia in Preclinical Cancer Models

There is a clinical need for noninvasive biomarkers of tumor hypoxia for prognostic and predictive studies, radiotherapy planning, and therapy monitoring. Oxygen-enhanced MRI (OE-MRI) is an emerging imaging technique for quantifying the spatial distribution and extent of tumor oxygen delivery in vivo. In OE-MRI, the longitudinal relaxation rate of protons (ΔR1) changes in proportion to the concentration of molecular oxygen dissolved in plasma or interstitial tissue fluid. Therefore, well-oxygenated tissues show positive ΔR1. We hypothesized that the fraction of tumor tissue refractory to oxygen challenge (lack of positive ΔR1, termed "Oxy-R fraction") would be a robust biomarker of hypoxia in models with varying vascular and hypoxic features. Here, we demonstrate that OE-MRI signals are accurate, precise, and sensitive to changes in tumor pO2 in highly vascular 786-0 renal cancer xenografts. Furthermore, we show that Oxy-R fraction can quantify the hypoxic fraction in multiple models with differing hypoxic and vascular phenotypes, when used in combination with measurements of tumor perfusion. Finally, Oxy-R fraction can detect dynamic changes in hypoxia induced by the vasomodulator agent hydralazine. In contrast, more conventional biomarkers of hypoxia (derived from blood oxygenation-level dependent MRI and dynamic contrast-enhanced MRI) did not relate to tumor hypoxia consistently. Our results show that the Oxy-R fraction accurately quantifies tumor hypoxia noninvasively and is immediately translatable to the clinic.

Radiomics: Images Are More than Pictures, They Are Data

In the past decade, the field of medical image analysis has grown exponentially, with an increased number of pattern recognition tools and an increase in data set sizes. These advances have facilitated the development of processes for high-throughput extraction of quantitative features that result in the conversion of images into mineable data and the subsequent analysis of these data for decision support; this practice is termed radiomics. This is in contrast to the traditional practice of treating medical images as pictures intended solely for visual interpretation. Radiomic data contain first-, second-, and higher-order statistics. These data are combined with other patient data and are mined with sophisticated bioinformatics tools to develop models that may potentially improve diagnostic, prognostic, and predictive accuracy. Because radiomics analyses are intended to be conducted with standard of care images, it is conceivable that conversion of digital images to mineable data will eventually become routine practice. This report describes the process of radiomics, its challenges, and its potential power to facilitate better clinical decision making, particularly in the care of patients with cancer.

Assessment of Tumor Angiogenesis: Dynamic Contrast-enhanced MR Imaging and Beyond

Dynamic contrast-enhanced (DCE) MR imaging is used increasingly often to evaluate tumor angiogenesis and the efficacy of antiangiogenic drugs. In clinical practice DCE-MR imaging applications are largely centered on lesion detection, characterization, and localization. In research, DCE-MR imaging helps inform decision making in early-phase clinical trials by showing efficacy and by selecting dose and schedule. However, the role of these techniques in patient selection is uncertain. Future research is required to optimize existing DCE-MR imaging methods and to fully validate these biomarkers for wider use in patient care and in drug development.

Early effects of low dose bevacizumab treatment assessed by magnetic resonance imaging

Background

Antiangiogenic treatments have been shown to increase blood perfusion and oxygenation in some experimental tumors, and to reduce blood perfusion and induce hypoxia in others. The purpose of this preclinical study was to investigate the potential of dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) and diffusion weighted MRI (DW-MRI) in assessing early effects of low dose bevacizumab treatment, and to investigate intratumor heterogeneity in this effect.

Methods

A-07 and R-18 human melanoma xenografts, showing high and low expression of VEGF-A, respectively, were used as tumor models. Untreated and bevacizumab-treated tumors were subjected to DCE-MRI and DW-MRI before treatment, and twice during a 7-days treatment period. Tumor images of K^trans (the volume transfer constant of Gd-DOTA) and v_e (the fractional distribution volume of Gd-DOTA) were produced by pharmacokinetic analysis of the DCE-MRI data, and tumor images of ADC (the apparent diffusion coefficient) were produced from DW-MRI data.

Results

Untreated A-07 tumors showed higher K^trans, v_e, and ADC values than untreated R-18 tumors. Untreated tumors showed radial heterogeneity in K^trans, i.e., K^trans was low in central tumor regions and increased gradually towards the tumor periphery. After the treatment, bevacizumab-treated A-07 tumors showed lower K^trans values than untreated A-07 tumors. Peripherial tumor regions showed substantial reductions in K^trans, whereas little or no effect was seen in central regions. Consequently, the treatment altered the radial heterogeneity in K^trans. In R-18 tumors, significant changes in K^trans were not observed. Treatment induced changes in tumor size, v_e, and ADC were not seen in any of the tumor lines.

Conclusions

Early effects of low dose bevacizumab treatment may be highly heterogeneous within tumors and can be detected with DCE-MRI.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1918-1) contains supplementary material, which is available to authorized users.

Association Between Computed Tomographic Features and Kirsten Rat Sarcoma Viral Oncogene Mutations in Patients With Stage I Lung Adenocarcinoma and Their Prognostic Value

BACKGROUND

We investigated the association between computed tomographic (CT) features and Kirsten rat sarcoma viral oncogene (KRAS) mutations in patients with stage I lung adenocarcinoma and their prognostic value.

PATIENTS AND METHODS

A total of 79 patients with pathologic stage I lung adenocarcinoma, available KRAS mutational status, preoperative CT images available, and survival data were included in the present study. Seven CT features, including spiculation, concavity, ground-glass opacity, bubble-like lucency, air bronchogram, pleural retraction, and pleural attachment, were evaluated. The association among the clinical characteristics, CT features, and mutational status was analyzed using Student's t test, the χ(2) test or Fisher's exact test, and logistic regression. The association among CT features, mutational status, and overall survival was analyzed using Kaplan-Meier survival curves with the log-rank test and Cox proportional hazard regression.

RESULTS

The prevalence of KRAS mutations was 41.77%. Spiculation was significantly associated with the presence of KRAS mutations (odds ratio, 2.99; 95% confidence interval [CI], 1.16-7.68). Although KRAS mutational status was not significantly associated with overall survival, the presence of pleural attachment was associated with an increased risk of death (hazard ratio, 2.46; 95% CI, 1.09-5.53). When analyzing KRAS mutational status and pleural attachment combined, patients with wild-type KRAS and no pleural attachment had significantly better survival than did those with wild-type KRAS and pleural attachment (P = .014).

CONCLUSION

These data suggest that spiculation is associated with KRAS mutations and pleural attachment is associated with overall survival in patients with stage I lung adenocarcinoma. Combining the analysis of KRAS mutational status and CT features could better predict survival.

Improved accuracy and precision of tracer kinetic parameters by joint fitting to variable flip angle and dynamic contrast enhanced MRI data

PURPOSE

To improve the accuracy and precision of tracer kinetic model parameter estimates for use in dynamic contrast enhanced (DCE) MRI studies of solid tumors.

THEORY

Quantitative DCE-MRI requires an estimate of precontrast T1 , which is obtained prior to fitting a tracer kinetic model. As T1 mapping and tracer kinetic signal models are both a function of precontrast T1 it was hypothesized that its joint estimation would improve the accuracy and precision of both precontrast T1 and tracer kinetic model parameters.

METHODS

Accuracy and/or precision of two-compartment exchange model (2CXM) parameters were evaluated for standard and joint fitting methods in well-controlled synthetic data and for 36 bladder cancer patients. Methods were compared under a number of experimental conditions.

RESULTS

In synthetic data, joint estimation led to statistically significant improvements in the accuracy of estimated parameters in 30 of 42 conditions (improvements between 1.8% and 49%). Reduced accuracy was observed in 7 of the remaining 12 conditions. Significant improvements in precision were observed in 35 of 42 conditions (between 4.7% and 50%). In clinical data, significant improvements in precision were observed in 18 of 21 conditions (between 4.6% and 38%).

CONCLUSION

Accuracy and precision of DCE-MRI parameter estimates are improved when signal models are fit jointly rather than sequentially. Magn Reson Med 76:1270-1281, 2016. © 2015 Wiley Periodicals, Inc.

Prognostic Imaging Biomarkers in Glioblastoma: Development and Independent Validation on the Basis of Multiregion and Quantitative Analysis of MR Images

PURPOSE

To develop and independently validate prognostic imaging biomarkers for predicting survival in patients with glioblastoma on the basis of multiregion quantitative image analysis.

MATERIALS AND METHODS

This retrospective study was approved by the local institutional review board, and informed consent was waived. A total of 79 patients from two independent cohorts were included. The discovery and validation cohorts consisted of 46 and 33 patients with glioblastoma from the Cancer Imaging Archive (TCIA) and the local institution, respectively. Preoperative T1-weighted contrast material-enhanced and T2-weighted fluid-attenuation inversion recovery magnetic resonance (MR) images were analyzed. For each patient, we semiautomatically delineated the tumor and performed automated intratumor segmentation, dividing the tumor into spatially distinct subregions that demonstrate coherent intensity patterns across multiparametric MR imaging. Within each subregion and for the entire tumor, we extracted quantitative imaging features, including those that fully capture the differential contrast of multimodality MR imaging. A multivariate sparse Cox regression model was trained by using TCIA data and tested on the validation cohort.

RESULTS

The optimal prognostic model identified five imaging biomarkers that quantified tumor surface area and intensity distributions of the tumor and its subregions. In the validation cohort, our prognostic model achieved a concordance index of 0.67 and significant stratification of overall survival by using the log-rank test (P = .018), which outperformed conventional prognostic factors, such as age (concordance index, 0.57; P = .389) and tumor volume (concordance index, 0.59; P = .409).

CONCLUSION

The multiregion analysis presented here establishes a general strategy to effectively characterize intratumor heterogeneity manifested at multimodality imaging and has the potential to reveal useful prognostic imaging biomarkers in glioblastoma.

Multiparametric imaging of patient and tumour heterogeneity in non-small-cell lung cancer: quantification of tumour hypoxia, metabolism and perfusion

Purpose

Multiple imaging techniques are nowadays available for clinical in-vivo visualization of tumour biology. FDG PET/CT identifies increased tumour metabolism, hypoxia PET visualizes tumour oxygenation and dynamic contrast-enhanced (DCE) CT characterizes vasculature and morphology. We explored the relationships among these biological features in patients with non-small-cell lung cancer (NSCLC) at both the patient level and the tumour subvolume level.

Methods

A group of 14 NSCLC patients from two ongoing clinical trials (NCT01024829 and NCT01210378) were scanned using FDG PET/CT, HX4 PET/CT and DCE CT prior to chemoradiotherapy. Standardized uptake values (SUV) in the primary tumour were calculated for the FDG and hypoxia HX4 PET/CT scans. For hypoxia imaging, the hypoxic volume, fraction and tumour-to-blood ratio (TBR) were also defined. Blood flow and blood volume were obtained from DCE CT imaging. A tumour subvolume analysis was used to quantify the spatial overlap between subvolumes.

Results

At the patient level, negative correlations were observed between blood flow and the hypoxia parameters (TBR >1.2): hypoxic volume (−0.65, p = 0.014), hypoxic fraction (−0.60, p = 0.025) and TBR (−0.56, p = 0.042). At the tumour subvolume level, hypoxic and metabolically active subvolumes showed an overlap of 53 ± 36 %. Overlap between hypoxic sub-volumes and those with high blood flow and blood volume was smaller: 15 ± 17 % and 28 ± 28 %, respectively. Half of the patients showed a spatial mismatch (overlap <5 %) between increased blood flow and hypoxia.

Conclusion

The biological imaging features defined in NSCLC tumours showed large interpatient and intratumour variability. There was overlap between hypoxic and metabolically active subvolumes in the majority of tumours, there was spatial mismatch between regions with high blood flow and those with increased hypoxia.

Electronic supplementary material

The online version of this article (doi:10.1007/s00259-015-3169-4) contains supplementary material, which is available to authorized users.

Toward understanding and exploiting tumor heterogeneity

The extent of tumor heterogeneity is an emerging theme that researchers are only beginning to understand. How genetic and epigenetic heterogeneity affects tumor evolution and clinical progression is unknown. The precise nature of the environmental factors that influence this heterogeneity is also yet to be characterized. Nature Medicine, Nature Biotechnology and the Volkswagen Foundation organized a meeting focused on identifying the obstacles that need to be overcome to advance translational research in and tumor heterogeneity. Once these key questions were established, the attendees devised potential solutions. Their ideas are presented here.

Cluster analysis of quantitative parametric maps from DCE-MRI: application in evaluating heterogeneity of tumor response to antiangiogenic treatment

PURPOSE

The objective of this study was to compare a clustering approach to conventional analysis methods for assessing changes in pharmacokinetic parameters obtained from dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) during antiangiogenic treatment in a breast cancer model.

MATERIALS AND METHODS

BALB/c mice bearing established transplantable her2+ tumors were treated with a DNA-based antiangiogenic vaccine or with an empty plasmid (untreated group). DCE-MRI was carried out by administering a dose of 0.05 mmol/kg of Gadocoletic acid trisodium salt, a Gd-based blood pool contrast agent (CA) at 1T. Changes in pharmacokinetic estimates (K(trans) and vp) in a nine-day interval were compared between treated and untreated groups on a voxel-by-voxel analysis. The tumor response to therapy was assessed by a clustering approach and compared with conventional summary statistics, with sub-regions analysis and with histogram analysis.

RESULTS

Both the K(trans) and vp estimates, following blood-pool CA injection, showed marked and spatial heterogeneous changes with antiangiogenic treatment. Averaged values for the whole tumor region, as well as from the rim/core sub-regions analysis were unable to assess the antiangiogenic response. Histogram analysis resulted in significant changes only in the vp estimates (p<0.05). The proposed clustering approach depicted marked changes in both the K(trans) and vp estimates, with significant spatial heterogeneity in vp maps in response to treatment (p<0.05), provided that DCE-MRI data are properly clustered in three or four sub-regions.

CONCLUSIONS

This study demonstrated the value of cluster analysis applied to pharmacokinetic DCE-MRI parametric maps for assessing tumor response to antiangiogenic therapy.