Abstract
Purpose: The aim of this study was to investigate the effect of VEGF-targeted therapy (sunitinib) on molecular intratumoral heterogeneity (ITH) in metastatic clear cell renal cancer (mccRCC).
Experimental Design: Multiple tumor samples (n = 187 samples) were taken from the primary renal tumors of patients with mccRCC who were sunitinib treated (n = 23, SuMR clinical trial) or untreated (n = 23, SCOTRRCC study). ITH of pathologic grade, DNA (aCGH), mRNA (Illumina Beadarray) and candidate proteins (reverse phase protein array) were evaluated using unsupervised and supervised analyses (driver mutations, hypoxia, and stromal-related genes). ITH was analyzed using intratumoral protein variance distributions and distribution of individual patient aCGH and gene-expression clustering.
Results: Tumor grade heterogeneity was greater in treated compared with untreated tumors (P = 0.002). In unsupervised analysis, sunitinib therapy was not associated with increased ITH in DNA or mRNA. However, there was an increase in ITH for the driver mutation gene signature (DNA and mRNA) as well as increasing variability of protein expression with treatment (P < 0.05). Despite this variability, significant chromosomal and transcript changes to key targets of sunitinib, such as VHL, PBRM1, and CAIX, occurred in the treated samples.
Conclusions: These findings suggest that sunitinib treatment has significant effects on the expression and ITH of key tumor and treatment specific genes/proteins in mccRCC. The results, based on primary tumor analysis, do not support the hypothesis that resistant clones are selected and predominate following targeted therapy. Clin Cancer Res; 21(18); 4212–23. ©2015 AACR.
This article is featured in Highlights of This Issue, p. 4025
Translational Relevance
The effect of targeted therapy on intratumoral heterogeneity (ITH) has not been studied in renal cancer. We report a study using a unique set of sunitinib-naïve and -treated clear cell renal cancer samples from patients with metastatic disease enrolled in clinical trials. There was significant ITH at protein, mRNA, and DNA levels. Following sunitinib therapy, there was increased intratumor heterogeneity of tumor grade, renal cancer driver gene mutations (mRNA and DNA), and protein expression. Despite this, consistent genetic alterations to key renal cancer genes occurred following sunitinib therapy. These results suggest that the commonly proposed mechanism of drug treatment enforcing clonal selection of a preexisting clone should be reconsidered. Instead it appears that the tumor is becoming molecularly more complex, which concurs with the increased clinical difficulty in treating patients who have developed resistance to first-line targeted therapy.
Introduction
A number of targeted agents are available for the treatment of metastatic clear cell renal cell cancer (mccRCC), of which a VEGF tyrosine kinase inhibitor (TKI), sunitinib, is most prominent (1). Clinical benefit with sunitinib varies between patients with mccRCC; however, resistance usually develops within a year in all patients. Unlike some other tumors (2), the mechanism of TKI resistance in mccRCC is unclear.
Recent analysis of primary and metastatic RCC tissue showed significant DNA intratumoral heterogeneity (ITH), driving the hypothesis that a process of clonal evolution occurs within RCC (3, 4). Theoretically, the presence of ITH is relevant in the development of resistance to VEGF targeted therapy. It is hypothesized that one of two possible scenarios occurs in resistance to targeted therapy: (i) a single resistant clone, potentially present at baseline, predominates during resistance resulting in reduced ITH with therapy; or (ii) a multifactorial process of acquired polyclonal resistance occurs, in which there is an initial increase in ITH with the potential for development of multiple resistant clones (5).
To evaluate the effect of TKIs on ITH and provide further evidence regarding likely mechanisms of resistance to targeted therapy, multiregion tissue sampling, and DNA, RNA, protein expression analysis was performed on a cohort of primary ccRCC tumors from sunitinib-naïve and -treated patients with mccRCC to identify if: (i) molecular heterogeneity also occurred at a RNA and protein level; and (ii) a period of sunitinib treatment (18 weeks) exacerbated molecular heterogeneity. Ideally, to answer the latter question one would undertake multiregion tumor sampling both before and after the sunitinib treatment, but before nephrectomy. Such an experiment would not be ethical or permissible. As such, we used a comparator group in the form of patients with well-matched sunitinib-naïve tumors. Sunitinib has specific molecular targets and has an effect on stromal elements within the tumor (4, 6). It is conceivable that, although sunitinib has little effect on ITH of the whole genome it may have specific effects on drug targets and driver mutations. Also, sunitinib may have inconsistent effects on the ITH of different components of the tumor, for example, increased ITH may be solely due to increased variability of host response to the drug in stromal elements (such as variable lymphocyte infiltration; ref. 4), rather than direct changes to variability of tumor clones. As such, we hypothesized that sunitinib treatment results in an increase in ITH of treatment specific genes, due to treatment related pressure. Therefore, in addition to unsupervised analysis, supervised analysis was undertaken for ITH in genes associated with: (i) key driver mutations in renal cancer; (ii) stromal elements; and (iii) molecular (mainly hypoxia related) targets of sunitinib. Understanding the effect of therapy on ITH gives insight into mechanisms of drug resistance and facilitates future biomarker work and drug development.
Materials and Methods
Patients and samples
Adequate fresh-frozen primary ccRCC tissue for analysis was obtained from the cytoreductive nephrectomy samples of 23 sunitinib-naïve patients with mccRCC as part of the Scottish Collaboration On Translational Research into Renal Cell Cancer (SCOTRRCC) study (UK CRN ID: 12229). Fresh-frozen primary tumor tissue was also obtained from 27 patients with mccRCC treated with three cycles of sunitinib (18 weeks) followed by a cytoreductive nephrectomy after 2 weeks off sunitinib as part of the Upfront Sunitinib (SU011248) Therapy Followed by Surgery in Patients with Metastatic Renal Cancer: a Pilot Phase II Study (SuMR; ClinicalTrials.gov identifier: NCT01024205); 23 of these patients had adequate tissue for analysis. Investigations were approved by institutional review boards and written, informed consent was obtained from each patient.
Tissue mapping
Tissue was obtained from at least three spatially separate regions of the kidney tumor (more regions if the size of the tumor allowed). Tissue processing was then performed as previously described (7). Briefly, if required each piece of tumor tissue was mapped and divided into smaller pieces (∼1cm3) in preparation for lysate creation. A cryostat section was performed on each 1 cm3 piece of tissue to confirm the presence of viable ccRCC and for grading (low grade, high grade, or mixed low/high grade) by two independent histopathologists (M. O'Donnell and D. Berney). A minimum of four protein/DNA/RNA lysates were aimed for per patient, from each morphologically differing region within each tumor.
Nucleic acid/protein extraction from tissue
Nucleic acid.
DNA extraction from fresh-frozen tissue was carried out using the DNeasy Blood and Tissue Kit (Qiagen) according to the manufacturer's instructions. The only significant alternation to the standard protocol was the addition of a PBS washing step to remove OCT compound carried over from the process of producing cryostat sections. Samples were fragmented before proteinase K digestion using a Qiagen Tissue Lyser. aCGH hybridization was carried out as described previously (8).
RNA extraction was carried out using the miRNEasy Kit (Qiagen) according the manufacturer's protocols. Amplification and analysis of the RNA was as previously described (9).
Protein.
Fifty to 75 mg of tissue was placed into 2-mL tubes with 990 μL of lysis buffer 50 mmol/L Tris–HCl (pH 7.5), 5 mmol/L EGTA (pH 8.5), 150 mmol/L NaCl) supplemented with aprotinin (10 mg/mL), phosphatase inhibitor cocktail A (Sigma), phosphatase inhibitor cocktail B (Sigma), and a protease inhibitor cocktail (Roche). A single 5-mm steel ball was added to each tube and the samples were homogenized at 50 Hz twice for 5 minutes using a TissueLyser (Qiagen). Triton X-100 (10 μL) was added to each sample before centrifuging at 13,000 × g for 30 minutes at 4°C after which supernatants were transferred to fresh microcentrifuge tubes. Total protein concentrations were determined by BCA assay (Sigma) and normalized to 1 mg/mL.
Array-based comparative genomic hybridization
aCGH hybridization and analysis was carried out as recently described (8) using Roche UK Nimblegen 12 × 135 k whole-genome array. The CGH-segMNT module of NimbleScan was used for the analysis with a minimum segment length of five probes and an averaging window of 130 kb. Nimblegen arrays were positionally annotated based upon hg19 genomic coordinates and log ratio data were preprocessed in R as previously described (10). Briefly, array data were normalized with print tip Loess from the limma package to produce normalized log ratios, filtered to remove outliers based upon a 1 MAD deviation of each probe from its immediate genomic neighbors and smoothed with a circular binary segmentation algorithm from the DNACopy package. Smoothed log ratios were then thresholded for gain/loss (±0.1) and amplification/deletion (±0.45) to identify contiguous copy-number aberrations containing at least 15 consecutive probes. For clustering categorical aCGH states representing gain, loss, amplification, and deletion were clustered based upon Euclidean distance using a Ward's algorithm using the segmented values from the whole genome in each case. aCGH data available via GEO (Accession number GSE67818).
VHL genotyping
VHL genotyping (East of Scotland Regional Molecular Genetics Service, Ninewells Hospital, Dundee, UK) included sequence analysis of exons 1-3 of the VHL gene to detect mutations in the coding region plus multiplex ligation-dependent probe amplification for large deletions.
Gene-expression analysis
mRNA was amplified using the WT-Ovation FFPE System Version 2 (NuGEN), purified using the Qiaquick PCR Purification Kit (Qiagen), biotinylated using the Encore BiotinIL ModuleIL (NuGEN), purified using the minElute Reaction Cleanup Kit (Qiagen) and quantified using a Bioanalyser 2100 with the RNA 6000 Nano Kit (Agilent). cRNA was then hybridized to Human HT-12v3 expression Beadarrays (Illumina) according to the standard protocol for NuGEN-amplified samples. Gene-expression data were read and normalized with the lump package in R using variance stabilization and robust spline normalization. Illumina expression data were annotated based upon ensembl gene annotation (hg19, release 61). Processed data were filtered to remove probes flagged as undetected in more than 50% of the 192 samples analyzed and probes with replicate gene identifiers were removed selecting those with the highest variance for each gene. For unsupervised clustering, the most informative probes were selected as those with a median absolute deviation across the dataset of greater than 0.6 leaving an intrinsic dataset representing 1,787 probes from unique genes. Filtered expression data were then clustered based upon Euclidean distance using a Ward's algorithm. Gene-expression data available via GEO (Accession number GSE65615).
Supervised CNV and gene-expression analysis
Supervised analysis was performed to assess clustering of individual patient tumor DNA and mRNA samples and copy-number variation (CNV) or differential gene expression with respect to molecules of relevance in renal cell cancer: (i) key molecular targets of sunitinib and hypoxia related genes; (ii) stromal associated genes (extracellular matrix, stromal remodeling, and adhesion molecules) were used as previously described (11); and (iii) a driver mutation gene set which was derived from those genes mutated in greater than 3% of the cases in the TCGA renal clear cell dataset. For ordered clustering, treatment groups were ordered according to the number of copy-number changes in each case from left to right in the aCGH figures and for the mean total expression of all the genes in the list in the expression heat maps.
Reverse phase protein arrays
Reverse phase protein array (RPPA) samples protein levels were determined and slides spotted as described previously (7, 12). Briefly, batch effects across the three slides per marker were mitigated using ComBat (13) and normalized (VSN; 14). Fifty-eight proteins making up key functional groups of proteins in RCC pathogenesis or sunitinib response were selected for study (Supplementary Tables S1 and S2). There was no signal detected for three of the proteins (Ki67, FLT3, and phospho-Jak2), leaving 55 proteins for analysis (antibodies detailed in Supplementary Table S1).
Statistical analysis
As described previously, aCGH (15) and gene-expression (10) cluster dendrograms were created and cut into 30 branches, those patients whose tissue samples all fell within the same branch were defined as clustered. Heterogeneity was assigned where an individual patient's samples were split between two or more clusters within the dendrogram.
Intratumoral protein variance was calculated separately for untreated and sunitinib-treated patients in an ANOVA framework. Untreated and treated variance distributions combining data from all 55 proteins were compared by a Mann–Whitney test (MW). Intratumoral variances for individual proteins were compared by an F test where assumptions of normality and homoscedasticity held, respectively, assessed using the Lillefours and Fligner tests; FDR correction was applied (16). Differential protein expression between untreated and treated patient samples was tested for each protein using a Student t test where normality and homoscedasticity assumptions were met, otherwise MW was performed; FDR was applied over combined t test and MW values (16).
The MW/t-test and χ2 test/Fisher exact test were used to test differences in continuous variables and categorical variables between the two groups, respectively. Pearson's correlation coefficient was used to calculate correlation using nonparametric data. P values of <0.05 were taken as significant and all P values were two-tailed. SPSS version 21 (IBM Corporation) or the R package (version 2.11.1) was used for all statistical analyses.
Results
Patient demographics
Of the 50 patients included in this study there was adequate tissue for analysis from 46 patients. The characteristics and clinical outcomes of the two patient cohorts (sunitinib-treated and untreated) are given in Table 1; the two groups were comparable other than for tumor grade (see below). Multiple tissue samples were obtained for all patients (total = 187 samples). The precise number of replicate samples for analysis from each patient was driven by the specific availability of tissue. Multiple protein lysates were obtained from 44 patients (median 4 regions; range, 2–10 regions). Multiple DNA lysates were obtained from 42 patients (median 3 regions; range, 2–5 regions). Multiple RNA lysates were obtained from 36 patients (median 3 regions; range, 2–8 regions).
Patient demographics, pathology details, and clinical outcomes of patients donating tumor samples to this study
Morphologic heterogeneity of sunitinib-naïve and -treated mccRCC samples
Tumor grade was assessed in multiple regions from each tumor to assess ITH (Table 1). A homogeneous grade (tumor contained either low or high grade, but not both grades) occurred in 26 (57.8%) samples whereas heterogeneous grading (tumors containing both low- and high-grade elements) occurred in 19 (42.2%). Significantly, more treated tumors expressed morphologic heterogeneity compared with untreated tumors (65.2% vs. 18.2%; P = 0.002).
Unsupervised aCGH analysis to assess global genetic instability and ITH
Comparison of aCGH results, from multiple DNA samples, taken from spatially different areas of each tumor showed variable intratumoral clustering (Table 2). In sunitinib-naïve tumors, 23 of 67 (34.3%) patient DNA samples clustered compared with DNA from 25 of 61 (41.0%) sunitinib-treated samples (P = 0.44, Fig. 1A), suggesting substantial underlying genetic ITH in mccRCC primary tumors, which was not significantly altered by sunitinib. Although the extent of ITH did not alter with sunitinib, there were significantly greater levels of chromosomal gains, amplifications and deletions in treated compared with untreated samples (P < 0.03; Fig. 1B).
Unsupervised aCGH and gene-expression clustering results. A, heatmap demonstrating clustering of aCGH profiles of sunitinib-naïve and -treated samples. Categorical aCGH states (i.e., gains, losses, and amplifications) were used for clustering, using Ward's clustering algorithm based upon Euclidean distance. No significant difference in aCGH clustering was seen between treated and untreated patients (P = 0.44; Fisher test). B, comparisons of the number of aCGH aberrations in treated and untreated samples showed significantly more copy-number gains, amplifications, and deletions in treated patients than in untreated patients (t test). C, heatmap demonstrating unsupervised clustering of gene-expression profiles of sunitinib-naïve and -treated samples. Filtered gene-expression data were clustered based upon Euclidean distance using a Ward's algorithm. There were no significant difference in sample clustering between treated and untreated patients (P = 0.69; Fisher test). Heatmap is a representation of all probes, created using clusters containing treatment-related genes and those driving the dendrogram.
Details of individual patient clustering of DNA CNV and RNA expression
Unsupervised gene-expression analysis to determine global transcriptomic ITH
Hierarchical clustering of gene-expression results from multiple RNA samples from each tumor showed that there was intratumor clustering in 24 of 36 patients (Fig. 1C and Table 2). There was a greater proportion of patients whose RNA clustered compared with DNA clustering (Table 2) 24 of 36 (67%) versus 17 of 42 (40.5%) patients, respectively (P = 0.025). As with the DNA analysis, there was comparable intratumoral clustering of untreated patient RNA samples, 26 of 43 (60.5%) RNA samples, compared with, 48 of 73 (65.8%) RNA samples from treated patients, P = 0.69.
Correlation of sample clustering for DNA and gene expression
There were 33 patients with clustering results available for both DNA and RNA. Eight patients (24.2%) had both DNA and RNA samples that clustered, 17 patients (51.5%) had either RNA or DNA that clustered and in 8 patients (24.2%) neither RNA nor DNA was determined to have clustered. Clustering of either DNA or RNA made clustering in the other nucleic acid more likely (P = 0.015).
Effect of sunitinib treatment on intratumor protein heterogeneity and protein expression
Intratumoral variance in expression of candidate proteins, measured using RPPA, was marked in untreated patient samples (Fig. 2A). But protein variance was even higher in the sunitinib-treated samples compared with the untreated samples (40/55, P = 0.001; Fig. 2B). Despite exacerbation of protein variance by sunitinib, there were significant differences in median expression of 30 of the 55 proteins evaluated (Supplementary Table S2). There was no correlation of the effect of sunitinib on differential RNA and protein expression (Supplementary Fig. S1). Consistent, significant changes for both mRNA and protein were only seen for 1 of the 55 proteins evaluated (Cyclin D1).
Protein results for sunitinib-naïve and -treated samples. A, histogram demonstrating the intratumoral variance for each of the 55 proteins analyzed with RPPA in sunitinib-naïve patients. B, log2 within tumor variance ratio of sunitinib-naïve and -treated samples. The number of proteins analyzed for which the median intratumoral variance is greater in the sunitinib-naïve group (n = 15); greater in the sunitinib-treated group (n = 40). There are significantly more proteins with higher median variance in the treated group than the untreated group (P = 0.001, binomial test). Proteins with variance below zero on the y-axis are greater for the untreated primary tumor, those above the y-axis zero were greater in the treated primary tumor. Supplementary Table S2 details the differential protein expression seen following sunitinib therapy.
Supervised aCGH and gene-expression analysis to assess ITH of relevant pathways
Three key groups of genes, relevant in renal cancer pathogenesis and TKI treatment were assessed to observe for ITH and to assess which tumor component may be driving protein ITH following sunitinib therapy. Key driver mutations, hypoxia-associated genes, and stroma-related genes were assessed for aCGH-derived CNVs (Table 2; Supplementary Fig. S2) and gene-expression analysis (Table 2; Supplementary Fig. S2). Considering CNVs, there was a significant reduction in individual patient sample clustering of driver mutation (P = 0.032) and hypoxia genes (P < 0.001) in treated patients, suggesting increased ITH.
For gene-expression analysis clustering, driver mutation genes, hypoxia, and stroma-associated genes had less clustering in treated samples than untreated patient samples (Table 2; Supplementary Fig. S2); this reduction in clustering was significant for driver mutation (P = 0.006) and stroma-associated genes (P = 0.016) but not hypoxia genes (P = 0.29).
Differential CNV and gene expression between untreated and treated patients
Comparison of untreated and treated DNA samples revealed the significant gains or losses in the region of the following key genes (Fig. 3). In the treated samples there were gains of VEGFR3 and PDGFRβ and losses for VCAM1, HIF1α, CAIX, mTOR, BAP1, PBRM1, SETD2, and VHL. Of note, no significant CNVs were seen to VEGFR2 (KDR).
Supervised analysis of driver mutation, hypoxia, and stromal molecules CNV between untreated and treated patients. Heatmaps of driver mutations (A), hypoxia (B), and stromal-related genes (C), plotting gains/amplifications (red) and losses/deletions (blue) in corresponding gene regions. The right-hand bar represents the lod score (−log10) of the adjusted P value (Fisher test), dashed line represents P = 0.05, bars are colored when significant (red, significant gain in treated patient samples; blue, significant loss in treated patient samples). The genes regions are ordered according by level of significance.
Similarly, gene-expression profiles of untreated and treated patient tumor RNA samples were compared (Fig. 4). Vimentin, clusterin, PDGFα, PDGFβ, VEGFA, VEGFR2 (KDR), HIF1α, and CAIX all had significantly greater expression in treated tumor samples. PBRM1, KDM5C, VHL, and PDGFC had significantly lower levels of expression in the treated samples.
Supervised analysis of driver mutation, hypoxia, and stromal molecules gene expression between untreated and treated patients. Heatmaps for gene-expression levels of driver mutations (A), hypoxia (B), and stromal-related genes (C). The right-hand bar represents the lod score (−log10) of the adjusted P value (t test), dashed line represents P = 0.05, bars are colored when significant (red, greater expression in treated patient samples; blue, lower expression in treated patient samples). The genes are ordered according to the level of significance.
Discussion
Initially, we showed that primary ccRCC tissue from patients treated with sunitinib expressed greater morphologic heterogeneity (of tumor grade) compared with tumors from untreated patients. This was intriguing and, in view of the similar characteristics of the patient cohorts, was postulated to be treatment rather than patient related. These findings built on our previous findings from paired, sequential biopsy, and nephrectomy samples, which showed that sunitinib therapy was associated with increased tumor grade, lymphocyte infiltration, and proliferation (Ki-67 expression; ref. 4). This increase in pathologic ITH prompted further investigation at the molecular level to determine whether sunitinib treatment was associated with dynamic molecular changes resulting in a more aggressive tumor phenotype. There has been extensive work on DNA heterogeneity at a single time point, resulting in the hypothesis of specific somatic mutations, following a pattern of clonal evolution (3). In the work, presented above it was hypothesized that sunitinib treatment results in an increase in ITH of treatment-specific genes, due to treatment related pressure; such clonal divergence may facilitate the evolution of resistant phenotypes.
To test the hypothesis, we performed chromosomal, RNA, and protein-based analyses. Results identified ITH at every molecular level investigated, in both treated and untreated samples. A comparison of mRNA and DNA results showed less ITH at a transcript level, suggesting ITH was predominantly genetic. Overall, approximately only half of multiple DNA or mRNA samples from individual tumors clustered together demonstrating the extent of the ITH. As has been previously demonstrated (17), there was remarkable inconsistency regarding mRNA and protein expression in both treated and untreated samples. Cyclin D1 was the only molecule for which an mRNA–protein correlation was found, potentially suggesting it for future multilevel molecular analysis. Despite this, there was a positive correlation between tumors with both DNA clustering and mRNA clustering, suggesting some consistency within tumors at different levels of molecular analysis. Therefore, the process of ITH was not purely random or molecular technique specific, but individual tumors appeared to have different degrees of ITH, which was consistent across DNA and mRNA.
The effect of VEGF treatment on ITH at a molecular level is unknown in renal cancer. Indeed, to our knowledge there is an absence in the literature for any tumor type. Unsupervised analysis did not reveal any significant effect of sunitinib therapy on ITH at a chromosomal or mRNA level. This finding may have been expected as sunitinib has specific effects, and is therefore unlikely to affect the whole-genome significantly. However, results from supervised analyses showed a consistent increase in ITH for both DNA and mRNA of the renal cancer driver mutation genes. The degree of DNA and RNA ITH of stromal- and hypoxia-related genes following sunitinib therapy was inconsistent. The consistent increase in ITH in the driver mutations and the lack of consistent results in the stromal genes suggest that the findings are, at least in part, tumor related rather than purely stromal-related reactions to therapy. Our results suggest dynamic changes to heterogeneity in the tumor, vascular and stromal components may be distinct with sunitinib. Whether changes to specific genetic alterations/expression are due to direct effects of treatment on the tumor or indirect effects on the tumor/stroma ratio remain uncertain. These findings were supported by protein expression results where an increase in ITH was observed for selected tumor specific proteins. It is conceivable that sunitinib causes diversification of protein expression within tumors as a reaction to treatment associated factors such as hypoxia, as postulated previously for other tumors (18).
The relatively short period of targeted treatment (18 weeks), consistent time between cessation of sunitinib and nephrectomy and the significant changes to selected relevant targets indicate that the exacerbated ITH in treated samples is likely to be treatment related rather than simply time related.
These results allow the development of further hypotheses regarding the mechanism of acquired resistant to VEGF-targeted therapy in mccRCC. The findings presented here do not support the hypothesis that a single resistant clone predominates from a number of clones after treatment commences, as treatment was associated with more rather than less ITH. As such, it can be speculated that there may be an expansion of clones, facilitating the development of resistant clones, that is polyclonal drug resistance (5). Further work using circulating samples or multiple biopsies at multiple time points throughout drug treatment would be required to prove the development of polyclonal drug resistance and allow the development of clinically relevant strategies to overcome drug resistance.
These findings may be relevant in optimizing the timing of immune checkpoint inhibitors in metastatic RCC. Such immunotherapies may work best in tumors with a high mutational burden (19). Therefore, starting checkpoint inhibitor therapy after a short period of VEGF targeted therapy, when mutational burden may be higher is attractive. This approach should be considered for future clinical trials.
The final finding was that despite ITH, consistent chromosomal, RNA and protein changes to the targets of sunitinib were seen with treatment, in keeping with our previous findings regarding consistent changes in CAIX protein expression following sunitinib treatment (12). Notable CNVs identified were losses of the renal cancer driver genes SETD2, BAP1, VHL, and PBRM1 (latter two genes also had a significantly lower level of gene expression) following sunitinib therapy (20). BAP1 and PBRM1 loss have previously been shown to be associated with a poor prognosis (21). It should be noted that, together with greater levels of 3p loss in the treated patient cohort, enrichment of the 3p tumor suppressors maybe due to intrinsic differences in the unmatched patient cohorts rather than treatment; although this is thought to be unlikely due to the comparable VHL mutation rate of 78% in the untreated patients and 65% in the treated patients. Overall, these genetic alterations suggest that sunitinib is having direct effects on the tumor and further support the hypothesis that the changes seen with sunitinib are not purely a time related phenomena.
We recognize the limitations of this study. First, although the characteristics of the two groups were similar, the treated and untreated samples frozen were not matched from the same patients. However, obtaining multiple frozen samples from nephrectomy patients is challenging, requiring close interaction between urologists and pathologists (22); the multiregion sample set described here is to our knowledge the largest available. The ideal validation set of matched pre- and posttreatment multiregion-sampled tumors is not feasible or ethical and a similar validation cohort to ours was unavailable. However, the multilevel analysis performed in this study, with all analyses confirming increasing ITH with treatment, acts as internal validation of the conclusions. Second, there was a 2-week break between sunitinib completion and nephrectomy, which may have influenced the results in the treated patients. Finally, sunitinib was given for a fixed period and responses to sunitinib were variable. Therefore, treated samples may have been at different stages of clonal evolution. Finally, there is a justified concern that identification of specific tumor related results might be contaminated by treatment-induced changes to the proportion of the stromal component within the tumor. The supervised analysis presented in this study attempted to address this concern; however, this is a consideration when interpreting the data presented.
The results of this study demonstrate a more complex tumor appears to be developing following sunitinib therapy. The resulting tumor is likely to be difficult to target successfully using a “one size fits all approach” as is currently standard of care in VEGF-targeted therapy resistant disease. These findings may in part explain why treatment in relapsed disease is particularly challenging (23).
Disclosure of Potential Conflicts of Interest
G.D. Stewart, A.L.R. Lubbock, I.M. Overton, D.J. Harrison and T. Powles are listed as co-inventors on a pending patent, owned by The University Court of The University of Edinburgh and Queen Mary University of London, on molecular methods for stratifying renal cell carcinoma patients for predicting response to therapy and/or selecting treatment. A. Bex is a consultant/advisory board member for Pfizer. T. Powles reports receiving speakers bureau honoraria from and is a consultant/advisory board member for GlaxoSmithKline and Pfizer, and reports receiving commercial research grants from Pfizer. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: G.D. Stewart, I.M. Overton, D.J. Harrison, T. Powles
Development of methodology: G.D. Stewart, L. Eory, A.L.R. Lubbock, P. Mullen, I.M. Overton, D.J. Harrison, T. Powles
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): G.D. Stewart, F.C. O'Mahony, A. Laird, M. O'Donnell, P. Mullen, S.A. McNeill, A.C.P. Riddick, D. Berney, A. Bex, M. Aitchison, D.J. Harrison, T. Powles
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): G.D. Stewart, F.C. O'Mahony, A. Laird, L. Eory, A.L.R. Lubbock, A. Mackay, J. Nanda, P. Mullen, I.M. Overton, D.J. Harrison, T. Powles
Writing, review, and/or revision of the manuscript: G.D. Stewart, F.C. O'Mahony, A. Laird, A.L.R. Lubbock, M. O'Donnell, P. Mullen, S.A. McNeill, A.C.P. Riddick, A. Bex, M. Aitchison, I.M. Overton, D.J. Harrison, T. Powles
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): G.D. Stewart, A. Laird, A.L.R. Lubbock, J. Nanda, D. Berney, A. Bex, D.J. Harrison, T. Powles
Study supervision: G.D. Stewart, I.M. Overton, T. Powles
Other (review of pathologic material): M. O'Donnell
Grant Support
This work was supported by the Chief Scientist Office, Scotland (ETM37; to G.D. Stewart, A.C.P. Riddick, M. Aitchison, and D.J. Harrison), Cancer Research UK (Experimental Cancer Medicine Centre; to T. Powles, London and D.J. Harrison, Edinburgh), Medical Research Council (to A. Laird and D.J. Harrison), Royal College of Surgeons of Edinburgh (to A. Laird), Melville Trust (to A. Laird), Medical Research Council (MC_UU_12018/25; to I.M. Overton), Royal Society of Edinburgh Scottish Government Fellowship cofunded by Marie Curie Actions (to I.M. Overton), Renal Cancer Research Fund (to G.D. Stewart), Kidney Cancer Scotland (to G.D. Stewart) and an educational grant from Pfizer (to T. Powles).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Acknowledgments
The authors thank the Wellcome Trust Clinical Research Facility Genetics Core, Edinburgh for assistance with undertaking gene-expression analysis, and Morad Ansari for performing aCGH.
Footnotes
Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).
- Received February 4, 2015.
- Revision received April 15, 2015.
- Accepted May 3, 2015.
- ©2015 American Association for Cancer Research.