Skip to main content
  • Research article
  • Open access
  • Published:

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

Abstract

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.

Peer Review reports

Background

To grow beyond a few millimeters in size, solid tumors need to establish vascular networks that can supply the tumor cells with oxygen and other nutrients [1]. The tumor cells produce and secrete several proteins that stimulate or inhibit angiogenesis, and the rate of angiogenesis is governed by the balance between these pro- and antiangiogenic factors [2]. Several antiangiogenic strategies are being investigated, including treatments with endogenous antiangiogenic facors or small peptides that mimic these factors [3, 4], monoclonal antibodies against proangiogenic factors or their receptors [5, 6], and tyrosine kinase inhibitors which may target multiple proangiogenic receptors [7]. The antiangiogenic treatments are generally not cytotoxic, and treatment-induced reductions in tumor volume often appear late compared to vascular effects [8]. It is therefore recognized that assessment of functional parameters are needed to detect early effects of antiangiogenic treatment.

Although antiangiogenic treatments may inhibit tumor growth when used alone, the therapeutic benefit may be even greater when used in combination with conventional therapies such as radiation and chemotherapy [9]. The effect of radiation and chemotherapy can be significantly affected by the tumor microenvironment, thus tumors with extensive hypoxia are more resistant to radiation and some forms of chemotherapy, and poor blood perfusion may reduce the uptake of chemotherapeutic drugs [10]. Antiangiogenic treatments have been reported to reduce blood perfusion and induce hypoxia in some experimental tumors [6, 11], and to increase blood perfusion and oxygenation in others [5, 12]. The reasons for these different effects are not well understood but may have significant impact on combination therapy [9]. It has been suggested that the effect of antiangiogenic treatment may vary with time after treatment, and that low doses of the antiangiogenic agent are required to increase blood perfusion and oxygenation [13]. It is also possible that the effect of antiangiogenic treatment may vary within tumors, although studies investigating this possibility are sparse.

Dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) and diffusion-weighted MRI (DW-MRI) have been used to evaluate the effect of antiangiogenic treatment [14]. In DCE-MRI, pharmacokinetic models are used to describe the tumor uptake of an intravenously administered contrast agent. The most common model is the generalized pharmacokinetic model of Tofts et al. [15]. In this model, the transfer rate constant, K trans, and the fractional distribution volume, v e, are estimated. K trans generally reflects blood perfusion and the vessel permeability - vessel surface area product, and v e reflects the extravascular extracellular volume fraction [15]. In DW-MRI, the apperant diffusion coefficient (ADC) is estimated. This parameter has been shown to reflect cell density and to be sensitive to necrotic tissue in untreated tumors [16, 17]. Reductions in K trans or K trans related parameters have been reported in most studies evaluating the effect of antiangiogetic treatment with DCE-MRI [14, 18], whereas both reductions and increases in ADC have been reported in studies evaluating the effect of antiangiogetic treatment with DW-MRI [19, 20]. In most of these studies, high doses of antiangiogenic agents have been used, and intratumor heterogeneity in the treatment-induced effects has not been investigated.

We have previously shown that DCE-MRI and DW-MRI are sensitive to effects of sunitinib treatment in human melanoma xenografts [21]. Sunitinib is a tyrosine kinase inhibitor which targets several receptors including vascular endothelial growth factor receptors 1-3 (VEGFR-1, -2, and -3), platelet-derived growth factor receptors α-β (PDGFR-α and PDGFR-β), stem cell growth factor receptor (c-KIT), and fms-like tyrosine kinase receptor 3 (FLT 3) [7]. In the previous study, we used a relatively high sunitinib dose which reduced microvascular density, increased hypoxic fractions, and induced necrosis. Moreover, the effect of treatment was evaluated once in xenografts from one melanoma line. In the current study, we evaluated the effect of low dose bevacizumab treatment with the same MR-techniques. Bevacizumab is a humanized monoclonal antibody that targets VEGF-A, and thus inhibits the VEGF-A pathway specifically [22]. Xenografts from two melanoma lines with different VEGF-A expression were used, and the tumors were subjected to DCE-MRI and DW-MRI before the treatment started and twice during a 7-days treatment period. We report that low dose bevacizumab treatment reduced K trans in the high VEGF-A expressing tumors, and that the effect was more pronounced in peripherial than in central tumor regions.

Methods

Mice and tumors

Adult (8–12 weeks of age) female BALB/c-nu/nu mice, bred at our research institute, were used as host animals for xenografted tumors. Animal care and experimental procedures were approved by the Institutional Committee on Research Animal Care and were performed in accordance with the Interdisciplinary Principles and Guidelines for the Use of Animals in Research, Marketing, and Education (New York Academy of Sciences, New York, NY, USA). The experiments were performed with tumors of the amelanotic human melanomas A-07 and R-18, established and characterized as described previously [23]. A-07 and R-18 cells were obtained from our frozen stock and were cultured in RPMI-1640 medium (25 mM HEPES and L-glutamine) supplemented with 13 % bovine calf serum, 250 mg/l penicillin, and 50 mg/l streptomycin. Approximately 3.5 × 105 cells in 10 Î¼l of Hanks’ balanced salt solution (HBSS) were inoculated intradermally in the hind leg by using a 100-μl Hamilton syringe.

Bevacizumab treatment

Mice were given two intraperitoneal doses of 5 mg/kg bevacizumab (Avastin, F. Hoffman-La Roche, Basel, Switzerland) or vehicle (saline), with 3 days between the doses.

Anesthesia

MRI was carried out with anesthetized mice. Mice were given 0.5 L/min O2 containing ~4.0 % Sevofluran (Baxter, Illinois, USA) during MRI. Respiration rate and body core temperature were monitored continuously by using an abdominal pressure sensitive probe and a rectal temperature probe (Small Animal Instruments, New York, USA). The body core temperature of the mice was kept at 37 Â°C by adjusting the hot air flow automatically, and the sevofluran dose was adjusted to maintain a stable respiration rate.

MR scanner and coil

MRI was performed by using a Bruker Biospec 7.05 T bore magnet with a mouse quadrature volume coil (Bruker Biospin, Ettlingen, Germany). The tumors were positioned in the isocenter of the magnet and were imaged with axial slices covering the entire tumor volume.

DCE-MRI

A fast spin echo pulse sequence (RARE) with varying repetition time (TR = 200, 400, 800, 1500, and 5000 ms), an echo time (TE) of 8.5 ms, an image matrix of 128 × 128, a field of view (FOV) of 3 × 3 cm2, a slice thickness of 0.7 mm, and a slice gap of 0.3 mm was applied to measure precontrast T1-values (T10-map). Gd-DOTA (Dotarem, Guerbet, Paris, France), diluted to a final concentration of 0.06 M, was administered in the tail vein of mice in a bolus dose of 5.0 ml/kg during a period of 5 s by using an automated infusion pump (Harvard Apparatus, Holliston, MA, USA). A 3-dimensional SPGR pulse sequence (3D-FLASH) with a TR of 10 ms, a TE of 2.07 ms, a flip angle (α) of 20°, an image matrix of 128 × 128 × 10, and a FOV of 3 × 3 × 1 cm3 was applied to produce T1-weighted images with a spatial resolution of 0.23 × 0.23 × 1.0 mm3, and a temporal resolution of 14.8 s. T1-weighthed images were recorded before Gd-DOTA injection, and every 14.8 s for 15 min after the injection (6 precontrast, and 59 postcontrast images). According to the theoretical equation for SPGR pulse sequences [24, 25],

$$ S = {S}_0\cdot \frac{sin\alpha \cdot \left(1-{e}^{-TR/{T}_1}\right)}{1- cos\alpha \cdot {e}^{-TR/{T}_1}}\cdot {e}^{-TE/{T}_2^{*}}\approx {S}_0\cdot \frac{sin\alpha \cdot \left(1-{e}^{-TR/{T}_1}\right)}{1- cos\alpha \cdot {e}^{-TR/{T}_1}} $$

where S is the signal intensity, and S 0 is a constant depending on scanner gain and proton density. The approximation \( {e}^{-TE/{T}_2^{*}}=1 \) is valid when T *2  ≫ TE, which was verified to be the case in our experiments. Images of phantoms with different Gd-DOTA concentration showed that the signal intensities produced by the 3D-FLASH followed the theoretical equation, confirming that the pulse sequence was appropriate for measurement of contrast agent concentration (Additional file 1). In contrast, the 2-dimensional SPGR pulse sequences available on the 7.05 T Bruker scanner (2D-FLASH) produced signal intensities that deviated substantially from the theoretical equation, and were thus inappropriate for measurement of contrast agent concentration (Additional file 1). Concentration of Gd-DOTA was calculated from the T1-weighted images in three steps. First, the constant S 0 was calculated for each voxel by using the precontrast images and the T 10 -map. Seccond, T 1 -values were calculated for the postcontrast images. Third, the changes in T 1 -values were converted to Gd-DOTA concentrations (C) by using the equation [25]:

$$ C\cdot {r}_{Gd- DOTA} = \frac{1}{T_1}-\frac{1}{T_{10}} $$

r Gd-DOTA is the relaxivity of Gd-DOTA which was measured to be 3.70 mM−1s−1 for the 7.05 T Bruker scanner. The DCE-MRI series were analyzed on a voxel-by-voxel basis by using the pharmacokinetic model described by Tofts et al. [15], and the arterial input function of Benjaminsen et al. [26]:

$$ {C}_t(T)=\frac{K^{trans}}{\left(1-Hct\right)}\cdot {\displaystyle {\int}_0^T{C}_a}(t)\cdot {e}^{- \frac{K^{trans}\cdot \left(T-t\right)}{v_e}}\ dt, $$

where C t (T) is the Gd-DOTA concentration in the tumor tissue at time T, K trans is the transfer rate constant, Hct is the hematocrit, C a (t) is the arterial input function, and v e is the fractional distribution volume of Gd-DOTA. Numerical values of K trans and v e were determined for each voxel from the best curve fit. Unphysiological voxels (voxels with v e > 1) were excluded from the analysis. The number of unphysiological voxels did not differ between untreated and bevacizumab-treated tumors and were ~ 5 % for A-07 tumors, and ~ 2 % R-18 tumors. Calculation of Gd-DOTA concentrations and pharmacokinetic modeling were done with in-house-made software developed in Matlab (MathWorks, Natick, MA, USA).

DW-MRI

DW-MRI was carried out by applying a diffusion-weighted single-shot fast spin echo pulse sequence (RARE) with a TR of 1300 ms, a TE of 26 ms, an image matrix of 64 × 64, a FOV of 3 × 3 cm2, a slice thickness of 0.7 mm, and a slice gap of 0.3 mm. Four different diffusion-weightings with diffusion encoding constants (b) of 200, 400, 700 and 1000 s/mm2, a diffusion gradient duration of 7 ms, and a diffusion separation time of 14 ms were used. Diffusion sensitization gradients were applied in three orthogonal directions with the following physical gradient combinations: [1 0 0], [0 1 0], [0 0 1]. ADC maps were produced with in-house-made software developed in Matlab. Briefly, the directional diffusion images were averaged on a voxel-by-voxel basis to non-directional diffusion images. ADC values were calculated for each voxel by fitting signal intensities (S) to the mono-exponential model equation:

$$ \ln \left(S(b)\right)=-b\cdot \mathrm{A}\mathrm{D}\mathrm{C}+c $$

by using a linear least square fit algorithm. The signal decay of a large number of voxels was investigated to verify that the mono-exponential model gave good fits to the data. The fits generally had a correlation coefficient of 0.95–0.99. DW-MRI was performed before injection of contrast agent.

Statistical analysis

Statistical comparisons of data were carried out by the Student’s t test when the data complied with the conditions of normality and equal variance. Under other conditions, comparisons were done by nonparametric analysis using the Mann–Whitney rank sum test. The Kolmogorov-Smirnov method was used to test for normality, and the Levene’s test was used to test for equal variance. Probability values of P < 0.05, determined from two-sided tests, were considered significant. The statistical analysis was performed by using the SigmaStat statistical software (SPSS Science, Chicago, IL, USA).

Results

Untreated A-07 and R-18 tumors were subjected to DCE-MRI and DW-MRI to investigate whether the MR-techniques were sensitive to differences between these tumor lines. A-07 tumors generally showed higher uptake of Gd-DOTA than R-18 tumors, and the uptake differed substantially for individual voxels in both tumor lines. This is illustrated in Fig. 1 which shows plots of Gd-DOTA concentration versus time and the corresponding pharmacokinetic model fits for individual voxels in a representative A-07 and R-18 tumor. The signal-to-noise ratio was sufficiently high that well-defined pharmacokinetic model fits were produced for voxels with both high and low uptake of Gd-DOTA in both tumor lines. The K trans, v e, and ADC image and the K trans, v e, and ADC frequency distribution of a representative A-07 and R-18 tumor are presented in Fig. 2a-b. Untreated A-07 tumors showed significantly higher K trans, v e, and ADC values than untreated R-18 tumors (Fig. 2c-e; P < 0.001).

Fig. 1
figure 1

Uptake of Gd-DOTA in individual voxels. Plots of Gd-DOTA concentration versus time (symbols) and the corresponding pharmacokinetic model fits (solid lines) for individual voxels in a representative untreated A-07 tumor (a), and a representative untreated R-18 tumor (b). K trans and v e values for the individual voxels were determined by pharmacokinetic analysis and are shown in the panels

Fig. 2
figure 2

A-07 and R-18 tumors differed in K trans, v e, and ADC. a-b, K trans, v e, and ADC images, and K trans, v e, and ADC frequency distributions of a representative untreated A-07 tumor (a), and a representative untreated R-18 tumor (b). The images refer to the central axial section of the tumor, whereas the frequency distributions are based on the individual voxel values of all the sections of the tumor. Color bars show K trans, v e, or ADC scales, scale bars are 2 mm, and the vertical lines in the frequency distributions indicate median values. c-e, K trans (c), v e (d), and ADC values (e) in untreated A-07 and R-18 tumors. Colums, means of 9-13 tumors, bars, SEM. Statistical comparisons of the data were carried out by the Student’s t test or the Mann–Whitney rank sum test. Untreated A-07 tumors showed significantly higher, K trans, v e, and ADC values than untreated R-18 tumors (P < 0.001)

A-07 and R-18 tumors were divided in groups with matched tumor sizes to receive bevacizumab treatment or vehicle. The tumors were subjected to MRI before the treatment started (day 0), and twice during the treatment period (day 3 and day 7), allowing accurate measurement of tumor volume at these time points. All tumors grew during the 7-days treatment period, and the bevacizumab-treated tumors did not differ from the untreated tumors in size at any time point, regardless of whether A-07 (Fig. 3a; P > 0.05) or R-18 tumors (Fig. 3b; P > 0.05) were considered.

Fig. 3
figure 3

Low dose bevacizumab treatment did not affect tumor growth. Tumor size before treatment, and 3 and 7 days after the treatment started in untreated and bevacizumab treated A-07 (a) and R-18 (b) tumors. Colums, means of 4-7 tumors, bars, SEM. Statistical comparisons of the data were carried out by the Student’s t test or the Mann–Whitney rank sum test. Significant differences in tumor size were not found between untreated and bevacizumab-treated A-07 tumors (a; P > 0.05), or between untreated and bevacizumab-treated R-18 tumors (b; P > 0.05)

A-07 and R-18 tumors differed in their response to low dose bevacizumab treatment. This is illustrated qualitatively in Fig. 4 which shows the K trans images and the K trans frequency distributions of a representative untreated A-07 tumor (Fig. 4a), a representative bevacizumab-treated A-07 tumor (Fig. 4b), a representative untreated R-18 tumor (Fig. 4c), and a representative bevacizumab-treated R-18 tumor (Fig. 4d). The images were recorded before treatment, and 3 and 7 days after the treatment start. Quantitative studies showed that K trans values were significantly reduced during growth in A-07 tumors (Fig. 5a; day 7 vs day 0: P = 0.032 for untreated tumors, and P = 0.015 for bevacizumab-treated tumors). After the treatment period, K trans values were lower in bevacizumab-treated than in untreated A-07 tumors, suggesting that the treatment reduced K trans (Fig. 5a). This difference was borderline significant when the absolute values of K trans were considered (P = 0.053) and significant when the K trans values were normalized to the pretreatment values (P = 0.032). In R-18 tumors, changes in K trans during growth were small, and significant differences between untreated and bevacizumab-treated tumors were not found, regardless of whether the absolute or normalized K trans values were considered (Fig. 5b; P > 0.05). v e and ADC values did not change during the treatment period for any of the tumor lines. Thus neither growth-induced nor treatment-induced changes in these parameters were observed (Fig. 5a-b; P > 0.05).

Fig. 4
figure 4

Low dose bevacizumab treatment reduced K trans in A-07 tumors. K trans images and K trans frequency distributions recorded before treatment, and 3 and 7 days after the treatment started. The panels show a representative untreated A-07 tumor (a), a representative bevacizumab-treated A-07 tumor (b), a representative untreated R-18 tumor (c), and a representative bevacizumab-treated R-18 tumor (d). The images refer to the central axial section of the tumor, whereas the frequency distributions are based on the individual voxel values of all the sections of the tumor. Color bars show K trans scales, scale bars are 2 mm, and the vertical lines in the frequency distributions indicate median values

Fig. 5
figure 5

The effect of low dose bevacizumab treatment on K trans, v e, and ADC. K trans, normalized K trans, v e, and ADC before treatment, and 3 and 7 days after the treatment started, in untreated and bevacizumab treated A-07 (a) and R-18 (b) tumors. Colums, means of 4-7 tumors, bars, SEM. Statistical comparisons of the data were carried out by the Student’s t test or the Mann–Whitney rank sum test. P-values are indicated in the panels where the statistical tests revealed significant or borderline significant differences between untreated and bevacizumab-treated tumors

To investigate intratumor heterogeneity in treatment effects, tumor images were divided in 5 concentric regions of interest (ROIs) as illustrated in Fig. 6a-b. Before the treatment period, A-07 tumors showed radial heterogeneity in K trans, i.e., K trans values were low in the central ROIs and increased gradually towards the tumor periphery (Fig. 6c; day 0, ROI 1 vs ROI 5: P < 0.001). In untreated A-07 tumors, K trans was similarly reduced in all ROIs during the treatment period, and thus tumor growth did not alter the radial heterogeneity (Fig. 6c; day 0 vs vehicle day 7). Compared with untreated tumors, bevacizumab-treated A-07 tumors showed similar K trans values in the central ROIs and significantly lower K trans values in the tumor periphery after the treatment period (Fig. 6c; vehicle day 7 vs bevacizumab day 7: P > 0.05 for central ROIs, and P = 0.016 for peripherial ROI). This implies that the treatment was more effective for peripherial than for central tumor regions, and consequently the treatment altered the radial heterogeneity in A-07 tumors. Radial heterogeneity in K trans was also found in R-18 tumors before the treatment period (Fig. 6d; day 0, ROI 1 vs ROI 5: P = 0.048). After the treatment period, bevacizumab-treated R-18 tumors did not differ from untreated R-18 tumors in any ROI, indicating that the treatment did not affect any region in these tumors (Fig. 6d; vehicle day 7 vs bevacizumab day 7: P > 0.05).

Fig. 6
figure 6

Intratumor heterogeneity in the effect of low dose bevacizumab treatment. a-b, K trans image of a representative untreated A-07 tumor (a), and image illustrating how the tumor was divided in 5 concentric circular ROIs (b). The circular ROIs are bounded by lines drawn at distances of nR/5 from the tumor center, where R is tumor radius and n is ROI number. Color bar shows K trans scale, scale bars are 2 mm. c-d, K trans in 5 concentric circular ROIs before treatment, and in untreated and bevacizumab-treated tumors 7 days after the treatment started. The graphs refer to A-07 (c) and R-18 (d) tumors. Symbols, means of 4–9 tumors, bars, SEM. Statistical comparisons of the data were carried out by the Student’s t test or the Mann–Whitney rank sum test. After the treatment period (day 7), differences in K trans values between untreated and bevacizumab-treated A-07 tumors were not significant in ROI 1–3 (P > 0.05), were borderline significant in ROI 4 (P = 0.063), and were significant in ROI 5 (P = 0.016). For R-18 tumors, significant differences between untreated and bevacizumab-treated tumors were not found in any ROI (P > 0.05)

Discussion

A-07 and R-18 melanoma xenografts were used as tumor models in the current study. We have previously shown that A-07 and R-18 cells differ substantially in the expression and secretion of VEGF-A, and that A-07 tumors have higher microvascular density, higher blood perfusion, and lower cell density than R-18 tumors [27–29]. In the current study, we demonstrate that these tumor lines also differ in the MR-derived parameters K trans, v e, and ADC, and we demonstrate that the tumor lines differ in their response to low dose bevacizumab treatment.

K trans generally reflects blood perfusion and the vessel permeability - vessel surface area product [15]. However, in tumors with a high and uniform vessel permeability, the uptake of small-size contrast agents like Gd-DOTA is not limited by the vessel permeability, and K trans reflects blood perfusion [30]. We have previously shown that A-07 and R-18 tumors show high and similar permeability for macromolecules, and that K trans reflects blood perfusion in these tumor lines [31, 32]. Consequently, the difference in K trans values between A-07 and R-18 tumors reported here probably reflected a difference in blood perfusion between the tumor lines. v e reflects the extravascular extracellular volume fraction which is inversely correlated to the cell density [15], and ADC has been shown to reflect cell density and to be sensitive to necrosis [16, 17]. Untreated A-07 and R-18 tumors show insignificant necrotic fractions but differ substantially in cell density [29]. The difference in v e and ADC values between A-07 and R-18 tumors thus probably reflected a difference in cell density between the tumor lines.

We have previously shown that sunitinib treatment reduces microvascular density and K trans values in A-07 tumors, suggesting that the sunitinib-induced reduction in K trans reflected reduced blood perfusion in that study [21]. It is highly likely that the treatment-induced reduction in K trans observed in the current study also reflected reduced blood perfusion, because both sunitinib and bevacizumab treatment inhibit the VEGF-A pathway [7]. Inhibition of the VEGF-A pathway has also been shown to reduce vessel permeability in experimental tumors [5], and consequently, the bevacizumab-induced reduction in K trans reported here may also have been influenced by reduced vessel permeability. The bevacizumab treatment did not affect v e and ADC values in any of the tumor lines, implying that the treatment did not change cell density and did not induce necrosis. This observation confirms that the bevacizumab dose was low. In contrast, increased ADC values reflecting induction of tumor necrosis have been observed after sunitinib treatment in A-07 tumors [21].

The different effect of low dose bevacizumab treatment between A-07 and R-18 tumors was probably a consequence of a difference in the rate of VEGF-A induced angiogenesis. Thus A-07 tumors show high VEGF-A expression, high microvascular density, and high pretreatment K trans values, and respond to bevacizumab treatment with reduced K trans values, whereas R-18 tumors show low VEGF-A expression, low microvascular density, low pretreatment K trans values, and no change in K trans values after bevacizumab treatment.

In A-07 tumors, bevacizumab treatment reduced K trans in peripherial regions with high pretreatment K trans values and had little or no effect in central regions with low pretreatment K trans values. We have previously demonstrated that untreated A-07 tumors show radial heterogeneity in several vascular parameters including microvascular density and blood perfusion [26, 33], suggesting that these tumors show similar heterogeneity in the rate of angiogenesis. These observations suggest that bevacizumab treatment was most effective in tumor regions with a high angiogenic rate. This suggestion is consistent with several studies reporting that antiangiogenic agents selectively removes immature blood vessels, because the fraction of immature blood vessels is expected to be high in tumor regions with high angiogenic rates [6, 12, 34].

In most studies evaluating the effect of antiangiogenic treatment with MR techniques, intratumor heterogeneity in the treatment effect has not been investigated [35]. Our study demonstrates that the effect of antiangiogenic treatment may be highly heterogeneous within tumors, and that careful monitoring of intratumor heterogeneities may provide mechanistic information about treatment effects and may identify poorly responding tumor regions. Detection of poorly responding regions can be important in a therapeutic setting because these regions may repopulate the tumor even if the treatment completely eradicates the tumor mass in other regions. If the effect of antiangiogenic treatment is evaluated with average parameters, poorly responding tumor regions may be overlooked. In addition, our study suggests that treatment-induced effects may be separated from growth-induced effects by evaluating changes in intratumor heterogeneities.

Treatment-induced reductions in tumor size generally occur late after antiangiogenic treatment [8]. However, if non-responding patients could be identified shortly after treatment initiation, any ineffective treatment could be stopped to avoid toxicity, and other treatments could be considered. In the current study, low dose bevacizumab reduced K trans without affecting tumor growth, suggesting that DCE-MRI may be used to identify patients that respond to low dose bevacizumab treatment before treatment-induced reductions in tumor size can be detected.

It has been suggested that antiangiogenic agents including bevacizumab can selectively remove immature blood vessels, increase tumor perfusion, and increase oxygenation [5, 12, 13]. These effects have been labeled vascular normalization and have been reported to occur within a limited time period [36]. Vascular normalization may be reversed if the treatment is stopped, and tumors may switch to other angiogenesis pathways during treatment and become resistant to the antiangiogenic agents. Moreover, the beneficial effects of vascular normalization may be balanced by severe vascular regression after prolonged exposure to antiangiogenic agents, or if the dose of the antiangiogenic agent is too large [36]. Appropriate timing and low doses are thus required to induce vascular normalization. It has also been demonstrated that inhibition of the VEGF-A pathway fails to normalize tumor vasculature and induces hypoxia in some preclinical tumors, suggesting that vascular normalization cannot be induced in all tumor models [6, 34]. In the current study, low dose bevacizumab treatment did not increase blood perfusion in A-07 and R-18 human melanoma xenografts. It is unlikely that the lack of increased blood perfusion was due to inadequate observation time points or to too large bevacizumab dose, because both similar and higher bevacizumab doses have been shown to increase blood perfusion and oxygenation at similar time points in several tumor models, including human breast carcinoma xenografts, human ovarian carcinoma xenografts, human neuroblastoma xenografts, murine melanoma, and murine breast carcinoma [5, 12 37]. The effect of low dose bevacizumab treatment reported here is similar to our previous experience with sunitinib treatment. In those studies, sunitinib treatment did not improve vascular function but induced hypoxia in A-07 and R-18 tumors [11, 21]. Taken together, our current and previous studies suggest that inhibition of the VEGF-A pathway does not induce vascular normalization in these melanoma lines.

In tumors where antiangiogenic treatment induces hypoxia, neoadjuvant antiangiogenic therapy is expected to reduce the effect of radiation and chemotherapy [9, 10]. In contrast, neoadjuvant bevacizumab treatment has been shown to enhance the effect of radiation and chemotherapy in preclinical tumors where bevacizumab normalizes the vasculature and the microenvironment [5, 12]. We have previously shown that DCE-MRI and DW-MRI can be used to identify tumors where antiangiogenic treatment does not normalize the microenvironment [21]. In that study, sunitinib treatment reduced K trans values and increased ADC values reflecting reduced perfusion and induction of tumor necrosis. The current study suggests that DCE-MRI can be used to identify such tumors, also when the treatment does not induce necrosis. These tumors respond to antiangiogenic treatment with reduced K trans values and no change in ADC values. Others have reported that vascular normalization results in increased K trans values and reduced ADC values [19 38]. Taken together, these studies suggest that DCE-MRI and DW-MRI may be used to monitor the effect of antiangiogenic treatment to detect vascular normalization, and to identify tumors where such treatment does not induce vascular normalization. Importantly, the MR-techniques are able to identify tumors where antiangiogenic treatment does not normalize the vasculature also when the treatment effect is small and tumor necrosis is not induced.

Conclusion

A-07 and R-18 tumors differed in the response to low dose bevacizumab treatment, and the response was associated with the rate of VEGF-A induced angiogenesis. Effects of low dose bevacizumab treatment were detected by DCE-MRI before tumor growth was affected. Our study suggests that DCE-MRI may be used to identify tumors where antiangiogenic treatment does not induce vascular normalization, also when the treatment does not induce necrosis. Moreover, the effect of low dose bevacizumab treatment was highly heterogeneous within A-07 tumors. Our study demonstrates that careful monitoring of intratumor heterogeneity may identify poorly responding tumor regions, may provide mechanistic information about the treatment effect, and may be used to differentiate treatment-induced from growth-induced effects in tumors similar to A-07 tumors.

References

  1. Folkman J. What is the evidence that tumors are angiogenesis dependent? J Natl Cancer Inst. 1990;82(1):4–6.

    Article  CAS  PubMed  Google Scholar 

  2. Carmeliet P, Jain RK. Angiogenesis in cancer and other diseases. Nature. 2000;407(6801):249–57.

    Article  CAS  PubMed  Google Scholar 

  3. Jia Y, Liu M, Huang W, Wang Z, He Y, Wu J, et al. Recombinant human endostatin endostar inhibits tumor growth and metastasis in a mouse xenograft model of colon cancer. Pathol Oncol Res. 2012;18(2):315–23.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Gaustad JV, Simonsen TG, Andersen LM, Rofstad EK. Thrombospondin-1 domain-containing peptide properdistatin improves vascular function in human melanoma xenografts. Microvasc Res. 2015;98:159–65.

    Article  CAS  PubMed  Google Scholar 

  5. Dickson PV, Hamner JB, Sims TL, Fraga CH, Ng CY, Rajasekeran S, et al. Bevacizumab-induced transient remodeling of the vasculature in neuroblastoma xenografts results in improved delivery and efficacy of systemically administered chemotherapy. Clin Cancer Res. 2007;13(13):3942–50.

    Article  CAS  PubMed  Google Scholar 

  6. Franco M, Man S, Chen L, Emmenegger U, Shaked Y, Cheung AM, et al. Targeted anti-vascular endothelial growth factor receptor-2 therapy leads to short-term and long-term impairment of vascular function and increase in tumor hypoxia. Cancer Res. 2006;66(7):3639–48.

    Article  CAS  PubMed  Google Scholar 

  7. Roskoski Jr R. Sunitinib: a VEGF and PDGF receptor protein kinase and angiogenesis inhibitor. Biochem Biophys Res Commun. 2007;356(2):323–8.

    Article  CAS  PubMed  Google Scholar 

  8. Morgan B, Horsfield MA, Steward WP. The role of imaging in the clinical development of antiangiogenic agents. Hematol Oncol Clin North Am. 2004;18(5):1183–206.

    Article  PubMed  Google Scholar 

  9. Horsman MR, Siemann DW. Pathophysiologic effects of vascular-targeting agents and the implications for combination with conventional therapies. Cancer Res. 2006;66(24):11520–39.

    Article  CAS  PubMed  Google Scholar 

  10. Brown JM, Giaccia AJ. The unique physiology of solid tumors: opportunities (and problems) for cancer therapy. Cancer Res. 1998;58(7):1408–16.

    CAS  PubMed  Google Scholar 

  11. Gaustad JV, Simonsen TG, Leinaas MN, Rofstad EK. Sunitinib treatment does not improve blood supply but induces hypoxia in human melanoma xenografts. BMC Cancer. 2012;12(1):388.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Dings RP, Loren M, Heun H, McNiel E, Griffioen AW, Mayo KH, et al. Scheduling of radiation with angiogenesis inhibitors anginex and Avastin improves therapeutic outcome via vessel normalization. Clin Cancer Res. 2007;13(11):3395–402.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Winkler F, Kozin SV, Tong RT, Chae SS, Booth MF, Garkavtsev I, et al. Kinetics of vascular normalization by VEGFR2 blockade governs brain tumor response to radiation: role of oxygenation, angiopoietin-1, and matrix metalloproteinases. Cancer Cell. 2004;6(6):553–63.

    CAS  PubMed  Google Scholar 

  14. Li SP, Padhani AR. Tumor response assessments with diffusion and perfusion MRI. J Magn Reson Imaging. 2012;35(4):745–63.

    Article  PubMed  Google Scholar 

  15. Tofts PS, Brix G, Buckley DL, Evelhoch JL, Henderson E, Knopp MV, et al. Estimating kinetic parameters from dynamic contrast-enhanced T(1)-weighted MRI of a diffusable tracer: standardized quantities and symbols. J Magn Reson Imaging. 1999;10(3):223–32.

    Article  CAS  PubMed  Google Scholar 

  16. Yankeelov TE, Arlinghaus LR, Li X, Gore JC. The role of magnetic resonance imaging biomarkers in clinical trials of treatment response in cancer. Semin Oncol. 2011;38(1):16–25.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Padhani AR. Diffusion magnetic resonance imaging in cancer patient management. Semin Radiat Oncol. 2011;21(2):119–40.

    Article  PubMed  Google Scholar 

  18. O’Connor JP, Jackson A, Parker GJ, Roberts C, Jayson GC. Dynamic contrast-enhanced MRI in clinical trials of antivascular therapies. Nat Rev Clin Oncol. 2012;9(3):167–77.

    Article  PubMed  Google Scholar 

  19. Batchelor TT, Duda DG, di Tomaso E, Ancukiewicz M, Plotkin SR, Gerstner E, et al. Phase II study of cediranib, an oral pan-vascular endothelial growth factor receptor tyrosine kinase inhibitor, in patients with recurrent glioblastoma. J Clin Oncol. 2010;28(17):2817–23.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Loveless ME, Lawson D, Collins M, Nadella MV, Reimer C, Huszar D, et al. Comparisons of the efficacy of a Jak1/2 inhibitor (AZD1480) with a VEGF signaling inhibitor (cediranib) and sham treatments in mouse tumors using DCE-MRI, DW-MRI, and histology. Neoplasia. 2012;14(1):54–64.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Gaustad JV, Pozdniakova V, Hompland T, Simonsen TG, Rofstad EK. Magnetic resonance imaging identifies early effects of sunitinib treatment in human melanoma xenografts. J Exp Clin Cancer Res. 2013;32:93.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Hurwitz H, Fehrenbacher L, Novotny W, Cartwright T, Hainsworth J, Heim W, et al. Bevacizumab plus irinotecan, fluorouracil, and leucovorin for metastatic colorectal cancer. N Engl J Med. 2004;350(23):2335–42.

    Article  CAS  PubMed  Google Scholar 

  23. Rofstad EK. Orthotopic human melanoma xenograft model systems for studies of tumour angiogenesis, pathophysiology, treatment sensitivity and metastatic pattern. Br J Cancer. 1994;70(5):804–12.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Hittmair K, Gomiscek G, Langenberger K, Recht M, Imhof H, Kramer J. Method for the quantitative assessment of contrast agent uptake in dynamic contrast-enhanced MRI. Magn Reson Med. 1994;31(5):567–71.

    Article  CAS  PubMed  Google Scholar 

  25. Yankeelov TE, Gore JC. Dynamic contrast enhanced magnetic resonance imaging in oncology: theory, data acquisition, analysis, and examples. Curr Med Imaging Rev. 2009;3(2):91–107.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Benjaminsen IC, Graff BA, Brurberg KG, Rofstad EK. Assessment of tumor blood perfusion by high-resolution dynamic contrast-enhanced MRI: a preclinical study of human melanoma xenografts. Magn Reson Med. 2004;52(2):269–76.

    Article  PubMed  Google Scholar 

  27. Rofstad EK, Halsør EF. Vascular endothelial growth factor, interleukin 8, platelet-derived endothelial cell growth factor, and basic fibroblast growth factor promote angiogenesis and metastasis in human melanoma xenografts. Cancer Res. 2000;60(17):4932–8.

    CAS  PubMed  Google Scholar 

  28. Tufto I, Lyng H, Rofstad EK. Vascular density in human melanoma xenografts: relationship to angiogenesis, perfusion and necrosis. Cancer Lett. 1998;123(2):159–65.

    Article  CAS  PubMed  Google Scholar 

  29. Tufto I, Rofstad EK. Interstitial fluid pressure, fraction of necrotic tumor tissue, and tumor cell density in human melanoma xenografts. Acta Oncol. 1998;37(3):291–7.

    Article  CAS  PubMed  Google Scholar 

  30. Tofts PS. Modeling tracer kinetics in dynamic Gd-DTPA MR imaging. J Magn Reson Imaging. 1997;7(1):91–101.

    Article  CAS  PubMed  Google Scholar 

  31. Bjørnæs I, Rofstad EK. Microvascular permeability to macromolecules in human melanoma xenografts assessed by contrast-enhanced MRI--intertumor and intratumor heterogeneity. Magn Reson Imaging. 2001;19(5):723–30.

    Article  PubMed  Google Scholar 

  32. Graff BA, Benjaminsen IC, Brurberg KG, Ruud EB, Rofstad EK. Comparison of tumor blood perfusion assessed by dynamic contrast-enhanced MRI with tumor blood supply assessed by invasive imaging. J Magn Reson Imaging. 2005;21(3):272–81.

    Article  PubMed  Google Scholar 

  33. Gaustad JV, Brurberg KG, Simonsen TG, Mollatt CS, Rofstad EK. Tumor vascularity assessed by magnetic resonance imaging and intravital microscopy imaging. Neoplasia. 2008;10(4):354–62.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Fenton BM, Paoni SF. The addition of AG-013736 to fractionated radiation improves tumor response without functionally normalizing the tumor vasculature. Cancer Res. 2007;67(20):9921–8.

    Article  CAS  PubMed  Google Scholar 

  35. O’Connor JP, Rose CJ, Waterton JC, Carano RA, Parker GJ, Jackson A. Imaging intratumor heterogeneity: role in therapy response, resistance, and clinical outcome. Clin Cancer Res. 2015;21(2):249–57.

    Article  PubMed  Google Scholar 

  36. Jain RK. Normalization of tumor vasculature: an emerging concept in antiangiogenic therapy. Science. 2005;307(5706):58-62.

  37. Moestue SA, Huuse EM, Lindholm EM, Bofin A, Engebraaten O, Maelandsmo GM et al. Low-molecular contrast agent dynamic contrast-enhanced (DCE)-MRI and diffusion-weighted (DW)-MRI in early assessment of bevacizumab treatment in breast cancer xenografts. J Magn Reson Imaging. 2013;38(5):1043-53.

  38. Sorensen AG, Emblem KE, Polaskova P, Jennings D, Kim H, Ancukiewicz M et al. Increased survival of glioblastoma patients who respond to antiangiogenic therapy with elevated blood perfusion. Cancer Res. 2012;72(2):402-7.

Download references

Acknowledgements

Financial support was received from the Norwegian Cancer Society and the South-Eastern Norway Regional Health Authority.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jon-Vidar Gaustad.

Additional information

Competing interests

The authors declare no competing interests.

Authors’ contributions

JVG, TGS, and EKR conceived and designed the study. JVG, RS, and LMKA performed the experiments. JVG, TGS, RS, CSW, and EKR analyzed and interpreted the data. JVG wrote the manuscript. All authors read and approved the final manuscript.

Additional file

Additional file 1:

Signal intensities in phantoms measured with 2D-FLASH and 3D-FLASH. (PDF 73 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gaustad, JV., Simonsen, T.G., Smistad, R. et al. Early effects of low dose bevacizumab treatment assessed by magnetic resonance imaging. BMC Cancer 15, 900 (2015). https://doi.org/10.1186/s12885-015-1918-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-015-1918-1

Keywords