Main

The inhibition of vascular endothelial growth factor (VEGF) signalling can suppress tumour growth by suppressing tumour angiogenesis. Agents designed to target this pathway include the VEGF-A-neutralising antibody, bevacizumab, and a range of anti-angiogenic receptor tyrosine kinase inhibitors (TKIs), including sunitinib, pazopanib, cediranib and regorafenib. Vascular endothelial growth factor pathway-targeted agents have been tested in numerous malignancies and have been approved in several indications, including metastatic breast, colorectal, lung, ovarian, and renal cancer; hepatocellular carcinoma; glioblastoma; and neuroendocrine tumours (Jayson et al, 2012; Sennino and McDonald, 2012). However, individual responses are variable, with some patients demonstrating resistance from the outset (intrinsic resistance) and others going on to develop resistance following a brief period of disease control (acquired resistance). Unfortunately, our understanding of what determines response and resistance to these agents is limited, and the answers to many fundamental questions are still unresolved. For example, what mechanisms underlie intrinsic and acquired resistance to these agents? Can response and resistance be predicted? How do we overcome resistance in order to improve the efficacy of this therapeutic approach? (Bergers and Hanahan, 2008; Ebos et al, 2009; Rini and Atkins, 2009; Ebos and Kerbel, 2011; Bottsford-Miller et al, 2012; Jayson et al, 2012; Sennino and McDonald, 2012).

In patients with metastatic renal cell carcinoma (mRCC), the administration of sunitinib or pazopanib as a single-agent therapy has been shown to extend progression-free and overall survival, and this form of therapy is now the standard of care for mRCC (Motzer et al, 2007, 2009; Sternberg et al, 2010). Nevertheless, 20% of mRCC patients demonstrate intrinsic resistance and derive no benefit. Moreover, although the majority of patients do achieve a period of disease control with these drugs, resistance is invariably acquired, typically within months of starting the treatment (Motzer et al, 2007, 2009; Sternberg et al, 2010; Busch et al, 2011; Heng et al, 2012). The current standard practice to treat with single-agent TKIs in this setting, combined with the manifestation of both intrinsic and acquired resistance, makes mRCC an important setting in which to explore the mechanisms that underlie response and resistance to anti-angiogenic therapy in patients.

Preclinical models show that, as well as blocking the formation of new tumour vessels, TKIs can induce the regression of newly formed immature vessels (Mancuso et al, 2006; Welti et al, 2011). Studies on tissue derived from renal cancer patients provide histological evidence that treatment of renal cancer with TKIs can induce a reduction in tumour vessel density (Griffioen et al, 2012; Powles et al, 2013). Therefore, TKIs can suppress the tumour vascularisation process, and this constitutes an important mechanism of action for these drugs. Although histological studies are informative, non-invasive longitudinal measurements of vascular parameters in tumours treated with anti-angiogenic agents are also desirable. For example, the amplitude of enhancement measured using arterial phase contrast-enhanced computed tomography (CECT) can be utilised as a surrogate measure of tumour vessel density in renal cancer. In support of this, correlative radiological–pathological studies show that the amplitude of the enhancement directly correlates with the tumour vessel density in renal cancer (Jinzaki et al, 2000; Wang et al, 2006). Moreover, when the measured enhancement is normalised to the enhancement of the aorta, inherent variations in phase, delivery and distribution of contrast between scans can be removed, thereby facilitating the direct quantitative comparison of enhancement (and of tumour vessel density, indirectly) in different lesions and in the same lesion at several different time points. Furthermore, changes in tumour vascular parameters measured using advanced imaging techniques have prognostic significance for mRCC patients treated with TKIs (O'Connor and Jayson, 2012).

In the current study, we have performed a retrospective analysis of arterial phase CECT scans obtained from 56 lesions in 23 mRCC patients treated with TKIs. We have examined whether changes in the normalised mean enhancement (a surrogate marker of tumour vessel density) correlate with response and resistance to anti-angiogenic therapy within individual metastases. We have assessed both the initial response to treatment (that is, between the baseline and the first on-treatment scan) and the development of acquired resistance to TKI treatment.

Materials and methods

Patients and lesions

The study group consisted of patients with mRCC commencing treatment with a TKI between July 2006 and October 2009 at the Mount Vernon Cancer Centre, Northwood, UK. Only patients with clear cell mRCC were included, following histological confirmation using tissue from either the primary or a metastatic lesion. In all cases, patients received TKIs as first-line therapy for metastatic disease. The TKIs used were sunitinib, pazopanib, cediranib or regorafenib. Details of the patient and lesion characteristics, the treatment setting and dosing regimens used are given in Tables 1 and 2. The activity of the TKIs utilised in this study (IC50 for receptor inhibition with in vitro kinase assays) is shown in Table 3. Patients were required to have undergone at least two contrast-enhanced multidetector CT scans: one at baseline (pre-treatment) and one first on-treatment scan, typically following two cycles of therapy, using images acquired in the arterial phase of enhancement. Where available, lesions were also analysed longitudinally until disease progression, again using images acquired in the arterial phase of enhancement. Two of the patients included in the current study were also included in a previous retrospective analysis (Nathan et al, 2010).

Table 1 Patient and lesion characteristics
Table 2 Treatment setting and dosing regimen
Table 3 Potency of TKIs for inhibition of selected tyrosine kinase receptors

Imaging

The following protocol was used routinely for all patients and at all time points reported in this study. For oral contrast, 250 ml of EZ-CAT (EZEM, Bicester, UK) diluted in 750 ml of water, was ingested 1 h before the CT was performed. An iodinated contrast agent was administered by the injection of 100 ml of 350 mg ml−1 Optiray (Covidien Healthcare, Mansfield, MA, USA) at 4 ml s−1 via a pump injector, followed by a 50-ml saline chaser at the same rate. The thorax and upper abdomen (from the supraclavicular fossa to the iliac crest) were imaged in the arterial phase (25 s after intravenous contrast agent administration). The abdomen and pelvis (from the dome of the diaphragm to the pubis) were then imaged in the portal venous phase (70 s after intravenous contrast agent administration). Patients underwent contrast-enhanced multidetector CT (Somatom Definition or Sensation 16, Siemens Healthcare, Forchheim, Germany) using the following acquisition parameters: 120 kV, effective mAs 150 (the thorax) or 200 (the abdomen and the pelvis) with dose modulation, pitch 1.2, detector configuration 24 × 1.2 mm, slice collimation 3 mm, scan field of view 300–450 mm, matrix 512 mm.

Image analysis

The image analysis was performed using data acquired in the arterial phase only. Target lesions with long axis of 10 mm were identified, and automated volume of interest software (syngoCT Oncology, Siemens Healthcare) was used to calculate parameters of interest. We calculated the volume of the lesion (following manual correction of the tumour outline) and the enhancement at each voxel (measured in Hounsfield Units, HU). From these data, the mean enhancement (mE) for the whole tumour volume of interest was calculated. To take into account possible variation in the phase, delivery and distribution of the contrast agent between scans, the mE was normalised to the enhancement measured in the aorta (aE). The aE was calculated by placing a region of interest within the aorta in the same axial plane as the lesion. The normalised mean enhancement (nE) for each lesion was then calculated as follows:

nE=mE/aE

The software also calculated the maximum diameter of the lesion (measured in mm). All measurements were conducted by NV, a medical oncologist, and were confirmed by VG, an oncological radiologist with 15 years of experience in CT imaging. For the purpose of this study, no limit was placed on the number of target lesions included per patient or within an organ. The baseline enhancement and enhancement at the first on-treatment scan were expressed in terms of nE (normalised mean HU, nHU). The lesion size at the baseline and at the first on-treatment scan was expressed in terms of the longest diameter (mm). The change in lesion size was always expressed in terms of the percentage change in the longest diameter (%). Change in enhancement was always expressed in terms of a percentage change in nE (%). Time was expressed in terms of days.

Analysis of PFS

Progression-free survival (PFS) was defined as the number of days on therapy until therapy was discontinued because of disease progression, as documented in the patient’s medical records. For the purposes of correlating features of the response with PFS on a per-patient basis, the baseline nE, the percentage change in nE (between the baseline and the first on-treatment scan) and the percentage change in the lesion size (between the baseline and the first on-treatment scan) was averaged across all assessable lesions for each patient (n=1–7 lesions per patient). An assessment of the response using modified Choi criteria (Nathan et al, 2010) was performed as follows: partial response (PR; 10% reduction in the lesion size and 15% reduction in the nE between the baseline and the first on-treatment scan), stable disease (SD; <10% reduction in the lesion size and <10% increase in size between the baseline and the first on-treatment scan, irrespective of change in enhancement) and progressive disease (PD; 10% increase in the lesion size between the baseline and the first on-treatment scan, irrespective of change in enhancement). Response by RECIST v1.1 was quantified as described (Eisenhauer et al, 2009). Kaplan–Meier analysis was then used to examine for correlations between (a) the baseline nE and PFS, (b) changes in the nE and PFS, (c) the response according to modified Choi criteria and PFS or (d) response according to RECIST and PFS.

Generation of pseudo-coloured images

Grey scale axial images acquired by arterial phase CECT were pseudo-coloured using an image analysis software (Siemens Multimodality Work Platform, Siemens Healthcare). A colour look-up table called the ‘PET rainbow’ was applied so that each enhancement value (measured in HU) was assigned to an output colour. The dynamic range of the colour look-up table was set at a minimum of −954 HU and a maximum of +422 HU. The resulting pseudo-coloured images are therefore colour-coded maps demonstrating the level of enhancement measured in each pixel.

Statistical analysis

The Pearson rank test was used to test for correlations between: (a) the baseline enhancement and change in lesion size between the baseline and the first on-treatment scan, (b) change in enhancement and change in lesion size between the baseline and the first on-treatment scan, (c) the baseline size and baseline enhancement, (d) size at the first on-treatment scan and enhancement at the first on-treatment scan, (e) change in size and days elapsed between defined time points, (f) change in enhancement and days elapsed between defined time points. The log-rank test was used to test for significant differences in PFS when patients were stratified according to the baseline enhancement, change in enhancement, modified Choi response criteria and RECIST. For lesions that acquired resistance to therapy, the Mann–Whitney U-test was used to test whether the change in size and the change in enhancement was significant when comparing lesions that underwent PR as best response with lesions that underwent SD as best response. P-values0.05 were considered to be statistically significant.

Results

Details of the study population

A total of 56 lesions from 23 patients were analysed, with a median of 2 lesions (and a range of 1–7) per patient. The baseline CT scan was performed at a median of 14 (4–38) days before commencing the therapy. The first on-treatment scan was performed at a median of 72 (33–113) days after starting treatment. Longitudinal monitoring of lesions, where patients had a baseline scan, followed by two or more on-treatment scans, was available for 21 lesions amongst 13 patients. Amongst these 21 lesions, we identified 10 lesions amongst seven patients that could be followed until disease progression. It was not possible to follow lesions until disease progression in the remaining patients because of the following reasons: the patient stopped taking the drug for reasons other than disease progression; the lesion was no longer visible or measurable; subsequent scans were non-contrast (owing to worsening renal function precluding the IV contrast agent administration: renal impairment with a serum creatinine>120 μmol l−1 or eGFR<60 ml min−1); or the patient was lost to follow-up.

Correlation of imaging features during the initial response to therapy

To address the initial response to therapy, lesion size and lesion enhancement were measured at baseline and at the first on-treatment scan for all 56 lesions in the study. For measurement of the lesion size, we used the longest diameter of the lesion. For measurements of enhancement, the mean enhancement (mE) present across the entire volume of the lesion was calculated. This was then normalised to the enhancement in the aorta (aE), to give the normalised mean enhancement (nE) for the entire lesion. Change in lesion size and change in nE were both expressed in terms of percentage change from baseline.

We found no correlation between the baseline nE and reduction in lesion size (Figure 1A). However, we found a strong correlation between reduction in nE and reduction in lesion size (Figure 1B, r=0.654, P<0.0001). A possible explanation for this observation is that smaller tumours merely have lower enhancement. To address this, we examined whether the absolute lesion size (longest lesion diameter) correlated with the amplitude of nE. However, no correlation between the absolute lesion size and nE was found either at the baseline or at the first on-treatment scan (Supplementary Figure 1A and B). The correlation between reduction in nE and reduction in tumour size might also reflect differences between patients in terms of the time elapsed between scans. However, no correlations were found between either the change in size or the change in nE and the duration of time elapsed (a) from the baseline scan to the start of treatment, (b) from the start of treatment to the first on-treatment scan and (c) from the baseline scan to the first on-treatment scan (Supplementary Figures 2–4).

Figure 1
figure 1

Correlation of imaging features with change in lesion size and patient outcome. (A) Analysis of correlation between the baseline nE and the percentage change in the lesion size (baseline to the first on-treatment scan). Baseline nE is measured in normalised mean Hounsfield Units (nHU). (r=−0.142, P=0.3). (B) An analysis of the correlation between the change in nE (baseline to the first on-treatment scan) and percentage change in lesion size (baseline to the first on-treatment scan; r=0.654, P<0.0001). Change in nE is measured in percentage. (C) Kaplan–Meier analysis: PFS of patients with baseline nE greater than the median (>0.39 nHU) and baseline nE less than the median (<0.39 nHU), P=0.2. (D). Kaplan–Meier analysis: PFS of patients with the change in nE greater than the median (>51.1%) and change nE less than the median (<51.1%), P=0.9. (E) Kaplan–Meier analysis: PFS of patients with a partial response (PR), stable disease (SD) or progressive disease (PD) according to modified Choi criteria, P=0.05.

Correlation of imaging features with PFS

Data on PFS were available for 22 of the 23 patients in the study. The median PFS was 298 days (a range of 74–1450), which is comparable to that of previously published cohorts (Sternberg et al, 2010; Motzer et al, 2012a, 2012b). We examined whether the baseline nE or the change in nE (from baseline to the first on-treatment scan) correlated with PFS. No statistically significant difference in the median PFS was observed when patients with a baseline nE above the median value (>0.39 nHU) were compared with patients with a baseline nE below the median value (<0.39 nHU; Figure 1C). No statistically significant difference in the median PFS was observed when patients having a percentage change in nE that was above the median value (>51.1%) were compared with patients having a percentage change in nE that was below the median value (<51.1%; Figure 1D). We also applied the modified Choi criteria, which take into account both the change in size and enhancement, to categorise patient responses (Nathan et al, 2010). Using these criteria, 11 patients achieved a PR, 9 patients achieved SD and 2 patients had progressive disease. The difference in median PFS between patients experiencing SD (using modified Choi criteria) and those experiencing PR (using modified Choi criteria) was statistically significant (211 vs 508 days, P=0.05, Figure 1E). In contrast, when patients were stratified using RECIST criteria, no significant difference in median PFS between patients experiencing SD (n=17) and those experiencing PR (n=3 patients) was observed (P=0.9).

Subgrouping of individual lesions based on change in size and change in enhancement

We then examined whether individual lesions could be grouped together into categories based on both change in size and change in enhancement. To do this, we devised criteria that categorised lesions into one of four different ‘response subgroups’ (Figure 2 and Table 3).

Figure 2
figure 2

Categorisation of individual lesions into four response subgroups. Waterfall plot shows the percentage change in lesion size (black) and the percentage change in normalised mean enhancement (nE; red) between the baseline and the first on-treatment scan for each of the 56 lesions in the study. Criteria based on the change in size and the change in nE are used to classify response: response subgroup 1 (30% reduction in size and 50% reduction in enhancement), response subgroup 2 (<30% reduction in size, but <20% increase in size and 50% reduction in enhancement), response subgroup 3 (<30% reduction in size, but <20% increase in size and <50% reduction in enhancement) and response subgroup 4 (30% increase in size and <50% reduction in nE).

Seven lesions (12.5%) underwent what could be described as a PR according to conventional RECIST criteria, that is, a decrease in the length of the longest diameter of the individual lesion by 30%. All seven of these lesions also underwent a corresponding large drop in enhancement (range of −63.1 to −87.3%). Therefore, this represents a subgroup of lesions (response subgroup 1) that undergoes both a good vascular response (50% reduction in nE) and significant tumour shrinkage (30% decrease in size; Figure 2). These lesions were spread across a range of anatomical sites: the liver (two), lymph node (two), adrenal gland (one), lung (one) and pancreas (one). Imaging data captured before and after treatment from two lesions in response subgroup 1 is shown in Figure 3A and B.

Figure 3
figure 3

Representative axial images from individual lesions at the baseline and at the first on-treatment scan. Pseudo-coloured arterial phase contrast-enhanced CT images of individual lesions from response subgroups 1, 2, 3 and 4 are shown. In each panel, the left hand image is the baseline scan and the right hand image is the first on-treatment scan. (A) Hilar node metastasis from a 57-year-old man treated with pazopanib (response subgroup 1). (B) Lung metastasis from a 79-year-old woman treated with pazopanib (response subgroup 1). (C) Pancreatic metastasis from an 82-year-old man treated with pazopanib (response subgroup 2). (D) Pleural metastasis from a 52-year-old woman treated with sunitinib (response subgroup 3). (E) Pre-tracheal lymph node metastasis from a 64-year-old man treated with sunitinib (response subgroup 3). (F) Liver metastasis from a 65-year-old man treated with cediranib (response group 4). (G) Waterfall plot shows the percentage change in lesion size (black) and the percentage change in nE (red) between the baseline and first on-treatment scan for each of the lesions depicted in panels (AF).

Forty-three lesions (76.8%) could be classified as SD according to RECIST; that is, they exhibited <30% decrease in the length of the longest diameter, but <20% increase in the length of the longest diameter for the individual lesion. Here, the change in enhancement was more varied (a range of −79.5 to +8.8%). Interestingly, 14 of these lesions underwent a large drop in enhancement (a range of −52.3 to −79.5%), which is a similar range to that of the change observed in the response subgroup 1. Therefore, in this second subgroup of lesions (response subgroup 2), despite undergoing a good vascular response (50% reduction in nE), only disease stabilisation was achieved (<30% decrease in size, but <20% increase in size; Figure 2). The anatomical locations of these lesions were: the pancreas (four), lung/pleura (three), liver (two), lymph node (two), adrenal gland (one), primary (one) and renal bed (one). Imaging data captured before and after treatment from one lesion in response subgroup 2 are shown in Figure 3C (and a second example is shown in Figure 4F).

Figure 4
figure 4

Heterogeneity of response between different lesions in the same patient. (AE) The responses of different lesions within the same patient were compared. Waterfall plots showing the percentage change in lesion size (black) and the percentage change in nE (red) between the baseline and first on-treatment scan for five patients. Lesions are assigned identification numbers from 1 to 4 or 1 to 7, as appropriate, and are classified according to response subgroup 1, 2, 3 or 4 as appropriate. (A) A 69-year-old man treated with sunitinib: lesions 1–3 are lung metastases and lesion 4 is a mediastinal lymph node metastasis. (B) A 57-year-old man treated with sunitinib: lesions 1,3 and 4 are liver metastases and lesion 2 is a hilar lymph node metastasis. (C) An 82-year-old man treated with pazopanib: lesion 1 is a hilar lymph node metastasis, lesions 2 and 3 are head of pancreas metastases and lesion 4 is a tail of pancreas metastasis. (D) A 69-year-old woman treated with sunitinib: lesions 1–4 are all lung metastases. (E) A 52-year-old woman treated with sunitinib: lesions 1–6 are lesions present in either the pleura or lung, whereas lesion 7 is a liver metastasis. (F and G) Pseudo-coloured arterial phase contrast-enhanced CT images of two lesions from the patient featured in panel E are shown. The baseline and first on-treatment scans of lesion 1, a pleural metastasis (F), and of lesion 7, a liver metastasis (G), are shown.

In the remaining 29 lesions that could be classified as SD using RECIST, the change in enhancement was more modest (with a range of −49.3 to+8.8%). This response subgroup 3 is characterised by a poor vascular response (<50% reduction in nE) and disease stabilisation (<30% decrease in size, but <20% increase in size; Figure 2). The anatomical locations of these lesions were: lung/pleura (11), lymph node (7), adrenal gland (5), liver (5) and chest wall (1). Imaging data captured before and after treatment from two lesions in response subgroup 3 are shown in Figure 3D and E.

Six lesions (10.7%) were classed as progressive disease by RECIST (20% increase in the length of the longest diameter for the individual lesion). In these lesions, the change in enhancement was also variable (a range of −37.1 to+57.8), but there was significant overlap with the change observed in response subgroup 3. Therefore, this is a subgroup of lesions (response subgroup 4) characterised by a poor vascular response (<50% reduction in nE) and disease progression (20% increase in size; Figure 2). The anatomical locations of these lesions were: lung/pleura (three), liver (two) and adrenal gland (one). Imaging data captured before and after treatment from one lesion in response subgroup 4 are shown in Figure 3F (and a second example is shown in Figure 4G).

Heterogeneity of response between lesions in the same patient

We selected the five patients in our study who had four or more measurable lesions at baseline and examined for heterogeneity of response between lesions. To define heterogeneity of response, lesions were classified as belonging to one of the four previously defined response subgroups (see Table 4). Patients were defined as undergoing a homogenous response if all of their lesions could be classified as belonging to the same response subgroup. Patients were defined as undergoing a heterogeneous response if different lesions were classified as belonging to different response subgroups. In two out of the five patients, a homogeneous response to therapy was observed, with all four lesions belonging to response subgroup 3 (Figure 4A and B). However, a heterogeneous response was demonstrable in patients 3, 4 and 5 (Figure 4C–E). For example, in the fifth patient, lesion 1, undergoing a good vascular response (50% decrease in nE) and disease stabilisation (Figure 4E), was classified as belonging to subgroup 2. Imaging data obtained at the baseline and at the first on-treatment scan are presented for lesion 1 (Figure 4F). However, lesions 2–5, undergoing a poor vascular response (<50% decrease in nE) and disease stabilisation (Figure 4E), were classified in subgroup 3. Finally, lesions 6 and 7 were classified in subgroup 4, undergoing a poor vascular response (<50% reduction in enhancement) and disease progression (Figure 4E). Imaging data obtained at the baseline and at the first on-treatment scan are presented for lesion 7 (Figure 4G).

Table 4 Response subgroups 1, 2, 3 and 4: criteria for subgrouping of individual lesions based on change in size and change in enhancement

Longitudinal assessment of lesions until progression

A significant majority of renal cell carcinoma metastases acquire resistance to TKIs following a period of disease control. Where the data were available, we monitored individual lesions longitudinally until they progressed. A total of 10 lesions across seven patients met our criteria for progressive disease (20% increase in size from nadir). In five of these lesions, the best response before progression was classified as a PR (according to RECIST criteria), with the anatomical sites being the liver (two), lymph node (two) and pancreas (one). In the remaining five lesions, the best response before progression was classified as SD according to RECIST, with the anatomical sites being the lymph node (three), adrenal gland (one) and lung (one).

An example of longitudinal data from a lesion, with best response classified as PR, that subsequently progressed on therapy is presented in Figure 5A and B. The initial reduction in tumour size between the baseline and the first on-treatment scan was paralleled by a substantial reduction in enhancement (50% reduction in nE; Figure 5A). It is noteworthy that compared with the baseline scan, the lesion converts into a more homogenous, hypoattenuating appearance at 41 days after the commencement of the treatment (Figure 5B). The lesion’s nadir was recorded in a scan obtained after 249 days of therapy, at which time point the nE recorded still remained low compared with that of the baseline (Figure 5A). However, a further scan obtained after 326 days of therapy indicated that the lesion had progressed (20% increase in size from nadir). The progression was accompanied by a dramatic increase in the nE of the lesion (Figure 5A). Note the presence of hyperattenuating regions throughout the tumour at this time point (Figure 5B). A final scan obtained after 524 days of therapy confirmed that the lesion continued to progress and that the lesion enhancement remained elevated (Figure 5A and B). Therefore, the acquired resistance to therapy in this lesion was clearly accompanied by re-vascularisation.

Figure 5
figure 5

Longitudinal monitoring of lesions until progression. Where data were available, individual lesions were monitored longitudinally until they progressed. (A and B) Example of a lesion that underwent PR as the best response and then progressed (an 82-year-old man treated with pazopanib, hilar lymph node metastasis). (C and D) Example of a lesion that underwent SD as the best response and then progressed (61-year-old woman treated with regorafenib, supraclavicular lymph node). Graphs in A and C show percentage change in size (black) or percentage change in nE (red), relative to baseline, quantified at all time points for which data were available. In B and D, pseudo-coloured arterial phase contrast-enhanced CT images of the baseline scan (baseline), first on-treatment scan (treated), scan at nadir (nadir), scan at progression (PD) and/or the final available scan (final) are shown (number of days on treatment are indicated in parenthesis). (E) Percentage reductions in size from baseline to nadir in the lesions that underwent PR or SD as best response are compared (P=0.01). (F) Percentage reductions in nE between the baseline and nadir in the lesions that underwent PR or SD as best response are compared (P=0.06). (G) Percentage increases in size from nadir to progression in the lesions that underwent PR or SD as best response are compared (P=0.3). (H) Percentage increases in nE from nadir to progression in the lesions that underwent PR or SD as best response are compared (P=0.008).

An example of longitudinal data from a lesion, with best response classified as SD that subsequently progressed on therapy, is presented in Figure 5C and D. Here, despite a substantial reduction in enhancement between the baseline and the first on-treatment scan (50% reduction in nE), only a modest reduction in tumour size was observed. The lesion’s nadir was recorded in a scan obtained after 109 days of therapy, at which time point the nE also remained reduced relative to the baseline. Interestingly, although disease progression was observed on the final scan obtained after 249 days of therapy, the parallel increase in nE was only modest. Given this observation, we plotted the change in size and change in nE between the baseline and nadir and between the nadir and progression for all lesions that underwent PR as best response, and for all lesions that underwent SD as best response (Figure 5E–H). The change in size between the baseline and nadir was significantly greater in the PR group compared with the SD group (Figure 5E; P=0.01). There was a trend towards a greater reduction in nE in the PR group compared with the SD group (Figure 5F; P=0.06). During progression, the increase in size observed in the PR group was not statistically different from the SD group (Figures 5G, P=0.3). However, the increase in nE recorded during progression, even among this small number of lesions, was significantly greater in the PR group compared with the SD group (Figures 5H, P=0.008). Therefore, rebound re-vascularisation during progression was more significant in lesions that underwent PR as the best response than it was in lesions that underwent SD as best response. Importantly, we found no correlation between the increase in nE from nadir to progression and the increase in tumour size (Supplementary Figure 5A). Moreover, the time elapsed between scans did not correlate with the increase in nE from nadir to progression (Supplementary Figure 5B).

As we also recorded six lesions that progressed during the initial response to therapy, it was of interest to examine whether, as a group, these lesions also underwent an overall increase in enhancement during their progression. Interestingly, in contrast to the scenario of acquired resistance, where there was an overall increase in the enhancement during progression (mean change in nE between nadir to progression=+344±107% in lesions undergoing PR as the best response, mean change in nE between nadir to progression=+22±11% in lesions undergoing SD as best response, Figure 5H), there was no overall increase in enhancement at progression in lesions that underwent PD as the best response (mean change in nE from baseline to progression=−5.4±15%).

Discussion

Understanding the mechanisms that determine response and resistance to anti-angiogenic agents is a key goal, with potential relevance to a broad range of cancers. Here, we utilised the amplitude of nE on arterial phase CECT imaging as a surrogate marker of tumour vessel density in order to explore how changes in tumour vessel density correlate with the response and resistance to angiogenic therapy in mRCC patients. In the following discussion, we place these observations in context with regard to other published work in this area and with regard to understanding the mechanisms of response and resistance to anti-angiogenic agents.

We did not find any relationship between baseline lesion enhancement and change in lesion size or between baseline lesion enhancement and PFS, suggesting that baseline tumour vessel density does not determine response. As the mechanisms that determine response and resistance to anti-angiogenic agents are complex (Bergers and Hanahan, 2008; Ebos et al, 2009; Bottsford-Miller et al, 2012; Sennino and McDonald, 2012), it is perhaps not surprising that tumour vessel density at baseline does not predict response to anti-angiogenic therapy. However, in contrast to our findings, some other studies have reported an association between baseline vascular parameters and response to TKIs in mRCC. In patients treated with either sunitinib or sorafenib, Han et al (2010) reported a positive correlation between high pre-treatment enhancement on CECT and both tumour size reduction and PFS, whereas Fournier et al (2010) reported a positive correlation between high tumour perfusion, measured by dynamic contrast-enhanced CT, and objective response. In addition, two DCI–MRI studies showed that high Ktrans at baseline was positively correlated with longer PFS in patients treated with sorafenib (Flaherty et al, 2008; Hahn et al, 2008). However, there were technical differences between these studies and the current study that may explain the discrepancy. First, in three of these studies (Flaherty et al, 2008; Hahn et al, 2008; Fournier et al, 2010), the imaging modality employed was different to the method used here. Second, all four studies (Flaherty et al, 2008; Hahn et al, 2008; Fournier et al, 2010; Han et al, 2010) measured the enhancement within a single ROI on axial images without normalisation to the cardiac output, whereas we measured the enhancement from the entire tumour volume and normalised our measurements to the cardiac output.

We analysed the change in enhancement and change in size that occurred within individual lesions between the baseline and the first on-treatment scan. We found that reduction in tumour enhancement in response to TKI treatment was strongly correlated with a reduction in lesion size. Interestingly, two previously published imaging studies show that the reduction in enhancement measured using CECT is correlated with tumour shrinkage and longer PFS in mRCC patients (Cowey et al, 2010a; Smith et al, 2010a). Moreover, a reduction in equivalent tumour vascular parameters measured using other imaging modalities has been shown to correlate with PFS (Flaherty et al, 2008) or OS (Lassau et al, 2010) in mRCC patients. Taken together, these data suggest that the ability of these agents to induce a vascular response is associated with tumour shrinkage and clinical benefit. Here, we have also found that modified Choi response criteria, which take into account early changes in both enhancement and size, can discriminate patients with a longer PFS from patients with a shorter PFS. This finding is in agreement with four other imaging studies demonstrating that variants of the Choi response criteria can predict benefit in mRCC patients treated with TKIs (Nathan et al, 2010; van der Veldt et al, 2010; Smith et al, 2010b; Krajewski et al, 2011).

In our study, we also classified individual lesions into four different response subgroups, using criteria that were based on both the vascular response and the change in size observed after treatment (Table 4, Figure 6A–D). For example, lesions in response subgroup 1 underwent a good vascular response (50% reduction in nE) that was accompanied by significant tumour shrinkage (Figure 6A). However, in response subgroup 2, a good vascular response (50% reduction in nE) translated into only stabilisation of disease (Figure 6B). It is not currently clear why a good vascular response can result in significant tumour shrinkage or just disease stabilisation, but preclinical work may provide some clues. For example, there are mechanisms that permit tumour cells to survive even when the vascular supply is significantly reduced. These survival mechanisms include a reduced propensity for tumour cells to die under conditions of stress and may be driven by genetic aberrations such as the loss of p53 function (Yu et al, 2001, 2002). More recent work shows that tumours exposed to anti-VEGF therapy show signs of increased autophagy (Hu et al, 2012a, 2012b; Guo et al, 2013), altered metabolism (Nardo et al, 2011; Welti et al, 2012; Zulato et al, 2012; Xu et al, 2013) and adaptation to hypoxia (McIntyre et al, 2012; Rapisarda and Melillo, 2012). These may all constitute potential mechanisms via which tumours adapt and survive, despite a reduction in their vascular supply. Therefore, adaptation of tumour cells to loss of a vasculature may be an important mechanism of intrinsic resistance to anti-angiogenic therapy and may explain, at least in part, why a good vascular response does not always result in significant tumour shrinkage.

Figure 6
figure 6

Diagrams illustrating the initial response to therapy and acquired resistance to therapy. (AD) Diagram illustrating how individual lesions may be classified into response subgroups based on the vascular response and change in size after treatment with anti-angiogenic therapy. (A) Response subgroup 1: significant vascular response and significant tumour shrinkage. (B) Response subgroup 2: significant vascular response and stabilisation of disease. (C) Response subgroup 3: poor vascular response and stabilisation of disease. (D) Response subgroup 4: poor vascular response and progression of disease. (E and F) Diagram illustrating that acquired resistance to anti-angiogenic therapy is associated with rebound re-vascularisation. However, rebound re-vascularisation may be more aggressive in lesions having PR as best response (E) when compared with lesions having SD as best response (F).

Lesions in response subgroups 3 and 4 were characterised by a weaker vascular response (<50% reduction in nE) and only achieved disease stabilisation (subgroup 3, Figure 6C) or in fact progressed (subgroup 4, Figure 6D). This weaker vascular response suggests that a significant proportion of vessels survived TKI treatment in these lesions, which may explain why more significant tumour shrinkage was not observed. The TKIs utilised in this study are proposed to target tumour angiogenesis principally by blocking the pro-angiogenic signalling of VEGF and PDGF, growth factors that promote the growth and survival of new tumour vessels. However, there is preclinical and clinical evidence for various mechanisms that may limit the ability of TKIs to suppress tumour angiogenesis. For example, tumour vessels are known to undergo a maturation process that renders them insensitive to VEGF inhibition, and tumours that have a large fraction of mature vessels are therefore less probable to respond to anti-angiogenic TKIs (Bergers and Hanahan, 2008; Sitohy et al, 2011, 2012). Moreover, there is evidence that some metastatic tumours obtain access to a vasculature via alternative non-angiogenic mechanisms, including vessel co-option, vascular mimicry or vessel intussusception (Dome et al, 2007; Leite de Oliveira et al, 2011). Interestingly, vessel co-option has been described in brain, lung and liver metastases and is a process whereby, instead of inducing new blood vessel formation, tumours co-opt existing local blood vessels as they invade into the surrounding host tissue (Pezzella et al, 1996; Holash et al, 1999; Vermeulen et al, 2001; Leenders et al, 2004; Sardari Nia et al, 2007). Therefore, another important source of intrinsic resistance to anti-angiogenic therapy may be the presence of tumour vessels that are not affected by the treatment due to mechanisms such as those described above. This may explain the poor response observed in response subgroups 3 and 4.

Here, we also document that different responses to therapy can occur in different lesions within the same patient. This observation is important because it suggests that, when exposed to the same anti-angiogenic therapy, some metastases can respond well, whereas others may simultaneously progress. This poses challenges in terms of optimal patient management and continuation of therapy, and may herald early treatment failure. The mechanisms that underlie this heterogeneity of response are unclear. Deep sequencing studies have demonstrated the existence of genetic heterogeneity between different mRCC lesions within the same patient (Gerlinger et al, 2012), suggesting that tumours in different metastatic sites in the same patient may have different biology. If this includes (a) differences in their ability to adapt to loss of a vasculature, and/or (b) differences in the sensitivity of their vasculature to TKIs, this may give rise to heterogeneity of response to anti-angiogenic therapy between different lesions. It will be important to understand the mechanisms that underlie this heterogeneity in order to improve on the efficacy of current anti-angiogenic treatment approaches.

The development of acquired resistance to anti-angiogenic agents is universally observed. Here, we clearly demonstrate that acquired resistance in patients is accompanied by re-vascularisation. Smith et al (2010a) have reported a similar observation. The mechanism that permits this re-vascularisation has yet to be confirmed. However, upregulation of other secreted pro-angiogenic factors, such as FGF2, IL8 or HGF, has been reported in acquired resistance and may act to stimulate tumour angiogenesis even in the presence of ongoing VEGF receptor blockade (Casanovas et al, 2005; Huang et al, 2010; Kopetz et al, 2010; Shojaei et al, 2010; Welti et al, 2011; Porta et al, 2012). These data highlight the possibility that acquired resistance might be circumvented if a more sustained inhibition of angiogenesis could be achieved with appropriate therapy. Interestingly, re-establishment of a vasculature was more pronounced in lesions that had PR as the best response, when compared with lesions that had SD as the best response (Figure 6E and F). We believe that this has not been previously reported. One interpretation of these data is that the biology of lesions that undergo PR as the best response may be different to that of lesions that undergo SD as the best response. These PR lesions may be intrinsically more dependent on a vascular supply for survival and therefore more dependent on re-vascularisation in order to progress. Alternatively, in lesions that undergo PR as the best response, there may be a stronger selective pressure for resistant tumour subclones that go on to drive TKI-resistant tumour vascularisation. Importantly, if the biological mechanism mediating acquired resistance in PR lesions is different from that of SD lesions, then this may influence the type of therapeutic approach necessary to combat acquired resistance in these different scenarios. However, it should be noted that the number of lesions on which this observation is based is small (n=10) and so this finding does require validation in a larger set of lesions.

Although the four TKIs utilised in this study target similar signalling pathways, the agents are administered on different schedules (see Table 2) and they differ subtly in their potency against individual receptors (see Table 3). Of interest, a trial designed to examine the comparative efficacies of scheduling sunitinib at 50 mg daily for 28 days followed by 14 days off treatment vs 37.5 mg per day continuously (EFFECT trial) failed to demonstrate any significant difference in outcomes between treatment arms (Motzer et al, 2012b). Moreover, a randomised phase III trial designed to compare the efficacies of sunitinib and pazopanib in mRCC (COMPARZ trial) failed to demonstrate any significant difference in outcomes between patients treated with the two agents (Motzer et al, 2012a). Therefore, there is currently no clear evidence of differences in response depending on the schedule or type of TKI used in mRCC. Despite this, we cannot entirely rule out the possibility that differences in scheduling and potency between agents may have an influence on the results observed in our study.

In conclusion, here we utilised the amplitude of enhancement on arterial phase CECT imaging as a surrogate marker of tumour vessel density in order to explore how changes in tumour vessel density correlate with response and resistance to angiogenic therapy in mRCC patients. Tumour vessel density at baseline was not a predictor of response to treatment in our study, but we found a strong relationship between reduction in tumour vessel density and tumour shrinkage, supporting the hypothesis that tumour growth inhibition in mRCC is closely linked with the extent to which these agents suppress tumour vascularisation. By classifying the response of individual lesions, both in terms of their vascular response and change in size, we present evidence for two different forms of intrinsic resistance to anti-angiogenic therapy: (a) the adaptation of tumour cells to withdrawal of a vascular supply and (b) mechanisms that allow the tumour vasculature to survive the therapy. We also document rebound vascularisation in tumours that acquire resistance to TKIs and provide evidence for a relationship between the initial response to therapy and the amplitude of rebound re-vascularisation. Therefore, although this was a retrospective study with a relatively small sample size, we believe our data provide important insights. Moreover, although these data were acquired from patients with mRCC, we reason that the principles identified here may also have wider relevance to other cancer types in which anti-angiogenic agents are utilised.