- Research article
- Open Access
Bevacizumab dose adjustment to improve clinical outcomes of glioblastoma
BMC Medicine volume 18, Article number: 142 (2020)
Glioblastoma (GBM) is one of the most aggressive and vascularized brain tumors in adults, with a median survival of 20.9 months. In newly diagnosed and recurrent GBM, bevacizumab demonstrated an increase in progression-free survival, but not in overall survival.
We conducted an in silico analysis of VEGF expression, in a cohort of 1082 glioma patients. Then, to determine whether appropriate bevacizumab dose adjustment could increase the anti-angiogenic response, we used in vitro and in vivo GBM models. Additionally, we analyzed VEGFA expression in tissue, serum, and plasma in a cohort of GBM patients before and during bevacizumab treatment.
We identified that 20% of primary GBM did not express VEGFA suggesting that these patients would probably not respond to bevacizumab therapy as we proved in vitro and in vivo. We found that a specific dose of bevacizumab calculated based on VEGFA expression levels increases the response to treatment in cell culture and serum samples from mice bearing GBM tumors. Additionally, in a cohort of GBM patients, we observed a correlation of VEGFA levels in serum, but not in plasma, with bevacizumab treatment performance.
Our data suggest that bevacizumab dose adjustment could improve clinical outcomes in Glioblastoma treatment.
Glioblastoma (GBM), a grade IV glioma according to the World Health Organization (WHO) classification, is one of the most common, aggressive, and highly vascularized brain tumors in adults, with a median survival of 20.9 months and an incidence rate of 3–4 newly diagnosed patients per 100,000 population [1,2,3]. The gold standard treatment consists of a maximal surgical resection, followed by radiotherapy with concomitant and adjuvant chemotherapy with temozolomide . Recently, it has been observed that using tumor-treating field (TTField) maintenance in combination with temozolomide has improved GBM overall survival (OS) with minor side effects . However, despite advances in surgical techniques and efforts in the development of new drugs, increasing patient survival rates is still a challenge.
GBMs are characterized by their low tumor cell differentiation, high proliferation and intracranial dissemination, elevated cellular heterogeneity, and abundant normal and aberrant vasculature [1, 5]. Pathological neovascularization is mostly triggered by hypoxia through the expression of vascular endothelial growth factor (VEGF) , which promotes angiogenesis, vasculogenesis, and vascular mimicry, cellular and molecular mechanisms that usually overlap each other .
Anti-angiogenic therapies have been shown to normalize the structure and function of tumor-associated blood vessels . Bevacizumab (Avastin®, Genentech/Roche), a recombinant humanized IgG1 monoclonal antibody that prevents the interaction of VEGFA with its receptors and inhibits downstream signaling pathways , was the first anti-angiogenic immunotherapy agent approved by the FDA to be used in GBM and in other cancer pathologies . Bevacizumab blocks the autocrine and paracrine VEGF signaling and reduces vessel leakiness and pressures within the brain tumor, triggering a vasogenic edema reduction and, consequently, lowering corticosteroid administration [11, 12]. Additionally, it may lead to vessel regression, thus depriving cancers of their nutrient source, which in turn may synergize the response to chemotherapy [13, 14]. However, some phase II and phase III clinical trials performed with bevacizumab in newly diagnosed (AVAglio-B021990 and RTOG-0825) and recurrent GBM patients failed to demonstrate a significant OS advantage, while revealing strong evidence of prolonged progression-free survival in GBM [15,16,17,18,19]. Indeed, the same results were observed in a recent meta-analysis carried out by Cochrane, in which authors emphasize that none of the clinical assays address the possibility that subsets of patients with GBM may benefit from bevacizumab therapy . Additionally, some GBM tumors grow in an angiogenesis-independent manner suggesting that these patients would probably respond poorly to anti-angiogenic therapy. Therefore, it is essential to identify a patient subpopulation that will benefit from the bevacizumab treatment . Although some predictive biomarkers have been suggested to be measured in plasma  and in solid biopsies , more robust and efficient biomarkers are needed to identify responder and non-responder patient subpopulations as well as to personalize the bevacizumab doses required to achieve antitumor effectiveness in the responder group.
We found that VEGFA expression levels could be used to calculate the specific dose of bevacizumab needed to increase the response to bevacizumab treatment. Additionally, GBM cell lines and patients with initial high VEGFA levels were more likely to respond positively to this anti-angiogenic therapy than other subpopulations of GBM with low VEGFA expression-secretion levels. Moreover, in a small cohort of GBM patients, we observed that VEGFA levels in serum, but not in plasma, correlated with response to bevacizumab treatment. Our data suggest that analyzing VEGFA levels in serum could be useful as a predictive biomarker to identify GBM patients that may benefit from bevacizumab treatment as well as to guide the recommended dosage.
Cell lines and culture
Three immortalized GBM cell lines (LN229, U87, and U373), human brain microvascular endothelial cells (HBMECs), and human mesenchymal stem cells (hMSCs) (a gift from Dr. Carmen Escobedo Lucea) were used. These cell lines were cultured in Dulbecco’s modified Eagles’ medium (DMEM) supplemented with 10% fetal bovine serum and 100 units/ml penicillin and 100 mg/ml streptomycin and were cultured in parallel at normoxia (21% O2) and hypoxia (1% O2). The GBM 12 OCT primary line was provided by Pilar Sanchez Gómez (Instituto de Salud Carlos III-UFIEC) , and HBMECs were isolated following the protocol published by Stins and colleagues .
Cancer stem cells (GBM27, GBM123, GBM128D, GBM38, GBM128B, and GBM18) were isolated from GBM surgical samples and cultured according to the protocol described previously . Tissue samples were obtained from patients operated at the Neurosurgery Department of Hospital Universitario la Fe (Valencia, Spain).
Bevacizumab treatment in vitro
Cells were grown for 72 h, and levels of VEGFA basal secretion were calculated by enzyme-linked immunosorbent assay (ELISA). Then, the amount of bevacizumab (Avastin®, Roche) needed to neutralize the VEGFA secretion was calculated in each case (Additional file: material and methods section).
Enzyme-linked immunosorbent assay (ELISA)
Basal levels of VEGFA in vitro or total serum levels isolated from mice and humans were calculated by ELISA, according to the VEGF Human ELISA Kit manufacturer’s instructions (Invitrogen). Absorbance detection was performed at 450 nm using a microplate reader.
RT and QRT-PCR
For RT-PCR, total RNA was isolated using RNeasy Mini or Micro kit (QIAGEN) following the manufacturer’s recommendations. One microgram of RNA was used for cDNA synthesis (High-Capacity cDNA Reverse Transcription Kit; Applied BioSystems). QRT-PCR reactions were performed in an optical 384-well plate equipped with an ABI PRISM 7900 HT sequence detection system (Applied Biosystems) using SYBR Green. Gene expression levels were quantified using the primers listed in Table 1, and two housekeeping genes were used to normalize the data.
Protein isolation and Western blotting
Cells and tumor specimens were lysed with radioimmunoprecipitation assay buffer (RIPA) containing a protease and phosphatase inhibitor cocktail (Roche). After centrifugation, the supernatant was collected and quantified using the BCA Protein Assay Kit (Calbiochem). Protein extracts were then separated by 12% SDS-PAGE and transferred to nitrocellulose membranes. After blocking for 1 h with 10% bovine serum albumin in T-TBS, the membranes were incubated with primary antibody against VEGFA (ab46154) overnight at 4 °C. After washing, the membranes were incubated with the corresponding peroxidase-labeled secondary antibody. Detection was performed using ECL reagents according to the manufacturer’s guidelines. Band intensities were quantified using Quantity One software (Bio-Rad Laboratories, Inc.) and were standardized to β-actin or GAPDH levels.
Cells were cultured in chamber slides and then fixed with 4% paraformaldehyde. After blocking for 1 h at room temperature (RT) with 10% BSA and 0.05% Triton X-100 in PBS, staining was performed using rabbit monoclonal anti-VEGFA (1:200, ab52917) followed by overnight incubation at 4 °C. Slides were washed and incubated for 1 h at RT with secondary antibody Alexa 555 goat anti-rabbit (1:200, Invitrogen) and Alexa Fluor 647 Phalloidin (1:100, Invitrogen). Then, the nuclei were stained with DAPI, (1:5000, Sigma-Aldrich), and slides were mounted with fluorescence mounting medium FluorSave™ reagent (Millipore). Fluorescence was observed under a Leica TCS SPE-inverted confocal microscope in the Confocal Unit of IiSGM. The experiment was performed twice, and in each event, three images of the different types of cells were taken. Negative controls with the secondary antibodies were carried out in all cases.
Tumor tissue sections were immersed in a 5% blocking solution of the specific serum and then incubated (overnight, 4 °C) in solutions containing the following primary antibodies: goat anti-mouse CD105 (R&D Systems) and mouse anti-human vimentin (Santa Cruz Biotechnology). Then, Alexa Fluor-conjugated secondary antibodies were used for 1 h (donkey anti-goat 555 and goat anti-mouse 488), and then the nuclei were counterstained with DAPI and coverslips were mounted using FluorSave™ reagent (Millipore). The immunofluorescence results were quantified using the Image-Pro Plus software where the fluorescence intensity of CD105 and VEGFA was calculated.
To produce conditioned media (CM), LN229, U87, and U373 were cultured for 72 h in the presence of IgG and standard or specific doses of bevacizumab. CM was sterile-filtered and stored at − 80 °C until use.
Matrigel™ Matrix High Concentration (BD Biosciences) was diluted with DMEM and used to coat 96-well plates at 37 °C for 30 min. Then, 2 × 104 HBMECs were resuspended in conditioned media DMEM 1:1 and seeded onto the Matrigel. The tube formation ability was measured after O/N incubation. The total length of the tubes and the number of branches were quantified using ImageJ.
Transwell migration assay
To evaluate cell migration, a 24-well plate with 8-μm pore size polycarbonate membrane inserts was used. Approximately 7.5 × 104 HBMECs were resuspended in 100 μl serum-free DMEM medium and seeded in the upper chamber, and 500 μl of conditioned media was added to the lower chamber. After 3 h, the migrated cells were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet. All experiments were performed in triplicate.
Wound healing assay
HBMECs were seeded in a 24-well plate overnight. The medium was removed, and the cells were washed twice with DMEM. A linear wound was created dragging a 200-μl pipette tip across the surface. Conditioned media from U87, U373, and LN229 cell lines with and without standard and specific doses of bevacizumab were added. The wound area was photographed with phase contrast at 0 h, and cells were allowed to migrate for 7 h. Scratched areas were measured with the ImageJ software, and the percent of wound closure was calculated.
xCELLigence real-time cell analysis (RTCA): migration
4 × 104 HBMECs were seeded per well into E16 plates (Roche) in DMEM, and conditioned media were loaded to the lower chamber. Then, the E16 plates were equilibrated to 37 °C and placed into the xCELLigence system. The impedance was recorded in 15 min intervals for 72 h, and the total number of cells migrating was quantified using a cell index.
In vivo Matrigel plug assay
Conditioned media of each cell line were concentrated using Amicon Ultra-15 (Millipore). Then, 150 μl of media concentrate was mixed with 350 μl Matrigel (BD Bioscience) and heparin. The mixture was injected subcutaneously into the flanks of 6-week-old immunodeficient mice (N = 5 mice per group). Mice were sacrificed after 7 days, and Matrigel plugs were removed.
Matrigel plugs and tumor tissues were fixed in 4% paraformaldehyde and embedded in paraffin, and 4-μm slices were obtained. The slides were stained with hematoxylin and eosin (H&E). Masson Trichrome Goldner staining was performed according to the manufacturer’s instructions (Bio Optica).
Heterotopic xenograft assay and bevacizumab treatment
1.5 × 106 LN229, U87, and U373 cells were resuspended in culture media, mixed with cold liquid Matrigel (BD), and then subcutaneously injected into athymic nude mice (N = 10 per group). Following Baumgarten’s recommendations for dose conversion between human/iv and mice/ip , specific and standard doses were calculated using in vitro hypoxic basal secretion. Final concentrations used were as follows: 60 mg/kg (U373-specific dose), 34 mg/kg (U87-specific dose), 10 mg/kg (LN229-specific dose), and 16.6 mg/kg bodyweight (standard dose). A control group (N = 10) received human polyclonal immunoglobulin (IgG, Intratect) at the same concentration as bevacizumab. Tumor volumes were measured with a caliper when they reached a visible size, and mice were euthanized when tumor burden became symptomatic.
Tissue microarray (TMA)
Eighty-seven tissue blocks were evaluated by a pathologist for the presence of malignant sections. Tumors were graded following the World Health Organization classification and divided into LGG and GBM (Additional file 1: Tables S1-2). Six TMAs were constructed from those formalin-fixed paraffin-embedded (FFPE) tissues using an arraying instrument (GALILEO CK 3500). Three tissue cores (0.6–1 mm diameter each) were made from each tissue block. Then, TMA blocks were cut at 4 μm and stained with hematoxylin & eosin. VEGFA signal was scored as follows: 0 = rare, 1 = localized, and 2 = diffuse.
FFPE and TMA sections were stained (as per the manufacturer’s staining protocol) with the Bond Polymer Refine Detection Kit on a Bond-max™ fully automated staining system (Leica Microsystems GmbH, Germany), using a rabbit polyclonal antibody against VEGFA (1:100, ab46154).
Solid surgical tissue samples and peripheral blood from patients and healthy donors were obtained from patients operated at HM Hospitales (HM), Madrid, Spain; Biobanco La Fe (PT17/0015/0043), Hospital Universitario la Fe (HUlaFe), Valencia, Spain; and the HGM BioBank, integrated in RETICS, National Network Biobanks, funding by Instituto de Salud Carlos III (Additional file 1: Tables S1-2). Peripheral blood samples from patients were collected prior to surgery, or prior and during Avastin treatment. These blood samples were left to clot for 30 min at room temperature, and the serum was isolated and stored at − 80 °C until use.
MR images were obtained using a 3-T MRI scanner (Achieva Intera 3T; Philips Healthcare, Best, the Netherlands) and 8-channel SENSE head coil. For DCE-MRI, baseline T1-weighted images were obtained with the following parameters: TR 3.5 ms, TE 1.7 ms, slice thickness 3.5 mm, field of view (FOV) 230 mm, matrix size of 116 × 117, 30 volumes, temporal resolution 5.7 s, and flip angles of 5° and 10° to create two pre-contrast datasets. Then, a DCE perfusion imaging dynamic series was performed using T1-weighted sequences with the same MR parameters except for an increased flip angle of 10°. At the end of the third volume acquisition, a bolus of 3 ml of gadolinium (Gd)-based contrast agent (gadobutrol 1 mmol/ml) was injected intravenously at a rate of 3 ml/s.
Structural contrast 3D T1 fast field echo (FFE) sequence was then performed, and the detail parameters were as follows: TR/TE = 5.8/2.7 ms, flip angle 8°, FOV 240 × 192 mm, matrix size 240 × 192, and reconstructed voxel size of 1 × 1 × 1 mm.
At that point, a T2-weighted perfusion MRI was performed. At the end of the third volume acquisition, a bolus of 6 ml of gadolinium was injected intravenously at a rate of 5 ml/s. The sequence was done with the following parameters: TR 17 ms, TE 25 ms, slice thickness 3.5 mm, field of view (FOV) 230 × 184 mm, matrix size of 80 × 60, 30 volumes, temporal resolution 2.6 s, and flip angles of 7°.
Parametric maps of cerebral blood volume (CBV) were acquired, and regional CBV values were calculated by a region of interest analysis between tumor and normal tissue. To observe BBB permeability, vascular constant transfer (Ktrans min− 1) values were calculated using syngo.via (Software Product Line for 3D Routine and Advanced Reading Siemens Healthcare). One region of interest was manually positioned on the solid tumoral area with maximal Ktrans value and compared with contralateral white matter.
Processed transcriptomic data form 981 glioma patients (151 astrocytoma, 109 oligoastrocytoma, 159 oligodendroglioma, and 562 glioblastoma) were obtained from the work published by Ceccarelli et al. in 2016 at https://tcga-data.nci.nih.gov/docs/publications/lgggbm_2015/  (Additional file 1: Tables S2-3). The histogram of the VEGFA expression of TCGA (Additional file 1: Figure S8 A) and our cohort (Additional file 1: Figure S8 B) was fitted with two Gaussian functions that represent two populations: high (green) and low (brown) VEGFA expression. The middle point between the positions of the two Gaussian functions was used to define the threshold between both populations.
High-volume data downloaded from the TCGA was managed using the R programming language.
For smaller datasets, statistical analysis was performed using a 2-tailed Student t test. Data are presented as means ± standard deviation and were calculated using the software package GraphPad Prism v. 5.0. Statistical values of P > 0.05 were not considered significant.
VEGFA expression level stratifies human glial tumors
To observe the clinical relevance of VEGFA, we first evaluated its mRNA expression in a cohort of 981 low- (LGG) and high-grade glioma (HGG) patients  (Additional file 1: Tables S2-3). Our results show that the level of VEGFA expression increased according to the histopathology degree. Distinct gene expression patterns were seen in the set of GBM samples, in which two different subgroups—high (green) and low (brown) expression—could be distinguished. Notably, there were 20% of GBM patients in the TCGA cohort and 25.9% in our cohort with low VEGFA levels. In the LGG population, VEGFA expression presented predominantly low values (Fig. 1a). Similar data were registered in our GBM cohort (N = 73); however, our LGG patients (N = 28) showed higher values than TCGA samples (Fig. 1b, Additional file 1: Table S1).
Results on VEGFA IHC staining are shown in Fig. 1c, d. 21.54% cases of de novo GBM and 52.94% of LGG patients showed rare or no expression of VEGFA protein (IHC score 0), 35.38% and 29.42% GBM and LGG respectively were scored as (1) localized and the rest of the patients and all the secondary GBMs had a very high score (2). Although the mean serum VEGFA concentrations (pg/ml) did not differ significantly between GBM and LGG, two secretion clusters could be distinguished in both groups. Moreover, VEGFA secretion values in both glioma groups were higher than in healthy donors (P < 0.01) and correlated with the tumor grade (Fig. 1e).
In concordance with the results observed in patients, distinct clusters based on VEGFA secretion and expression were detected in ten GBM cell lines, which showed a high correlation between supernatant VEGFA levels and mRNA values (Fig. 1f, g and Additional file 1: Figure S1). Remarkably, these patterns were reproduced in intracranial xenografts, showing a clear relation between in vitro and in vivo behavior (Fig. 1h). Together, these results demonstrate that glioma patients could be divided according to VEGFA expression and secretion levels.
Specific but not standard bevacizumab doses impair VEGFA expression and secretion
In view of the importance of VEGFA in the angiogenesis pathway and to better understand its role in GBM, we selected one representative cell line of each subgroup, U87 as high, U373 as medium, and LN229 as low expression.
The second set of experiments examined the ability of bevacizumab to inhibit VEGFA secretion in vitro at normoxic and hypoxic conditions. Although bevacizumab could compete with the detection antibody, ELISA assay is able to measure the total VEGFA present in the extracellular media . For this purpose, we used a standard dose (SD; 8.3 μg/ml) which mimics the clinical administration and a specific dose (Spe) calculated for each cell line (see Additional file 1: methods section). As can be seen in Fig. 2a and Additional file 1: Figure S2, the Spe calculated for each cell line (Additional file 1: Table S4) is able to neutralize VEGFA secretion in the three cell cultures studied. Once we evaluated both conditions, normoxia was established for the following experiments. Although VEGFA presence in the extracellular media was reduced using the SD, its complete inhibition was not reached. Next, we wondered if the inhibition of extracellular VEGFA could have an impact on VEGFA expression levels, as autocrine loop has been described . VEGFA transcript and protein levels were only decreased using Spe in U87 and U373 (Fig. 2b, c). Similarly, SD had only a small effect on cellular VEGFA expression in U87 and U373, but Spe significantly inhibited the VEGFA expressed in these cell lines (Fig. 2d). Bevacizumab Spe did not change the cell morphology (Additional file 1: Figure S3), and no effects were observed in hMSCs viability (Additional file 1: Figure S4) showing low toxicity levels. Our results prove that the neutralization of the secretion of VEGFA is bevacizumab dose-dependent.
Specific dosage of bevacizumab inhibited HBMEC migration and tube formation in tumor cells
In an effort to elucidate the underlying mechanism of bevacizumab in angiogenesis, we exposed HBMECs to conditioned media (CM) obtained by treating U87, U373, and LN229 cell lines with standard (SD) or specific doses (Spe) of bevacizumab. In each case, IgG were used as a control (C).
The antivascular potential of bevacizumab was first tested by tube formation assay after 16 h of CM incubation. As shown in Fig. 3a, the branch reduction number was only significant when U87 and U373-Spe-CM were used (P < 0.05). Moreover, U87-Spe-CM mediated the inhibition of tube formation compared with the U87-C-CM and SD-CM (P < 0.01, P < 0.05). In the case of U373-CM, the total network length was reduced using Spe-CM (P < 0.01). No effect was observed when LN229-CM was used. These data demonstrate that U87 and U373-Spe-CM can efficiently reduce HBMEC tube formation.
To study the cell migratory behavior, we carried out a transwell migration assay (Fig. 3b). After 3 h, the number of HBMECs that moved through the membrane with U87-C-CM was 6.7-fold higher than the number of cells incubated with U87-Spe-CM (P < 0.01). A similar percentage of migrated cells was observed within the U373 group, with statistically significant differences between the U373-C-CM versus SD and Spe (P < 0.05).
In the cells incubated with LN229-CM, only a slight reduction was observed when SD-CM was used (P < 0.05). These results suggest that the addition of a specific dose of bevacizumab inhibited the migration of endothelial cells.
Furthermore, a wound healing assay showed that the addition of bevacizumab reduced HBMEC motility in the U87 and U373-CM. No changes in cell migration were observed in the wells incubated with LN229-CM (Fig. 3c). Consistent with these observations, we noticed a marked decrease in the HBMEC migration activity using U87 and U373-Spe-CM by the xCELLigence Real Time Cell Analyzer (Fig. 3d) (controls without CM are shown in Additional file 1: Figure S5). Importantly, we found that Spe-CM causes the disruption of the autocrine VEGF/VEGFR2 signaling loop as shown in the dramatic decrease of VEGFR2 expression (Additional file 1: Figure S6).
Altogether, our results indicate a strong decrease in angiogenesis in vitro after the treatment of U87 and U373 cells with a specific dose of bevacizumab.
Specific dosage inhibits vascularization of Matrigel plugs from in vivo mouse models
Once we confirmed that CM from U87 and U373 treated with Spe of bevacizumab is sufficient to switch off the angiogenesis in vitro, we wondered whether the pathway would be limited in vivo. We concentrated conditioned media (c-CM) and mixed it with Matrigel and heparin. Then, it was injected into the mouse flanks. Macroscopic control plug images appeared red, whereas plugs treated with SD and Spe appeared pale in the three cell lines studied. The components of the cellular infiltrate could be observed in H&E and Masson images (Fig. 4a). The histological analysis demonstrated that the Spe markedly reduced blood vessel growth compared with the control group in the c-CM obtained from U87 and U373. Furthermore, the c-CM from U87 Spe and U373 SD and Spe, but not from LN229, presented a sustained decrease in neovascularization (Fig. 4b), which was confirmed in the significant reduction of specific markers for mouse endothelial cells such as CD31 and von Willebrand factor (VWF) (Fig. 4c). SD had no inhibition effect in U87. On the contrary, we observed an increase in endothelial marker expression (Fig. 4c, d). Analysis of angiogenesis-related factors (VEGFA, VEGFB, and VEGFR2) in Matrigel plugs obtained from c-CM-U87 and c-CM-U373 evidenced a high decrease in expression after the addition of the Spe of bevacizumab. We observed that only the addition of the Spe reduced VEGFA expression, while the use of both bevacizumab concentrations diminished VEGFB and VEGFR2 expression (Fig. 4d). These data confirm that the Spe of bevacizumab is needed to inhibit angiogenesis in vivo.
Specific dosages showed antitumor activity in angiogenic cell lines
To further confirm that the Spe of bevacizumab could only reduce tumor growth in high angiogenic cell lines, we implanted U87, U373, and LN229 subcutaneously into the flanks of nude mice. Bevacizumab treatment was administered twice a week following Baumgarten’s recommendations for dose extrapolations between human/iv and mice/ip . Remarkably, U87 and U373 tumor volume and growth were drastically suppressed by SD and in a more significant manner with Spe. However, bevacizumab treatment did not affect tumoral activity in the in vivo xenograft models formed by the low VEGFA-expressed cell line LN229 (Fig. 5a–c). More importantly, the amount of VEGFA secreted by only tumoral cells was quantified from mouse serum using a human-specific test. Values obtained demonstrate that SD could only reduce the amount of VEGFA ligand in the U87 xenograft model, whereas the Spe extended its capability to U87 and U373 models (P < 0.05). No effect was detected in LN229 tumors (Fig. 5d).
VEGFA protein expression in tissue was reduced in a dose-dependent manner in U87 and U373 (P < 0.01). Surprisingly, bevacizumab treatment increased VEGFA expression in LN229 tissue sections (Fig. 5e). The study of new blood vessel formation by CD105 confirms that bevacizumab treatment is only effective in cell lines with high VEGFA secretion (Fig. 5f), confirmed by the substantially reduced protein levels compared to the control (Fig. 5g). H&E staining revealed no changes in cell morphology; however, U373-xenografted mice treated with Spe had lower cellularity (Additional file 1: Figure S7). Taken together, our results suggest that two different groups of patients could be stratified depending on basal VEGFA secretion.
VEGFA serum concentration key to predict the response to bevacizumab therapy in glioma patients
In order to accomplish our study objective, we selected the patients treated with bevacizumab following the scheme depicted in Fig. 6a. The results obtained from surgical solid biopsy support the idea that some GBM patients have a tumor with a high angiogenic phenotype (patients 8 and 2), while in others, VEGFA protein expression is absent (patient 11) (Fig. 6b). Considering that there is no consensus for VEGFA measurement in the peripheral blood, we evaluated it in plasma and serum, and we observed that the serum mean concentration was almost 2.5 times the plasma value (Fig. 6c). Interestingly, the correlation between VEGFA expression in tumor tissue and serum (r = 0.87, P < 0.05) was much stronger than between the levels in tumor tissue and plasma (r = 0.75, P = 0.153) (Fig. 6d). To further confirm this observation, we collected blood prior to and after bevacizumab treatment. As expected, a strong association was shown between VEGFA serum levels, radiographic response, and RANO criteria (Additional file 1: Table S5). These values correlate with CBVr and Ktrans and could distinguish between potential bevacizumab responder (Fig. 6e–i) and non-responder patients (Fig. 6j–n and Additional file 1: Figure S9). No association was observed in the plasma (Additional file 1: Figure 1S0).
Glioma growth and progression are dependent on neovascularization , which involves endothelial cell proliferation, migration, and tube formation . Since VEGFA is the main pro-angiogenic factor, it has been used as a target in several pathologies . However, the clinical benefit in glioma patients has been questioned considering the results obtained in clinical trials , raising much controversy . Some possible explanations for the lack of therapeutic effects of bevacizumab are the acquired resistance observed in some patients and the fact that more than 30% of placebo arm patients received bevacizumab during disease progression . From our in silico analysis of TCGA data, we observed two different populations based on VEGFA expression which were confirmed in mRNA, protein, and serum levels in our own cohorts. The data suggest that there are more than 20% of GBM patients whose tumor growth is independent of the VEGFA pathway, so they would not benefit from an anti-angiogenic treatment.
Our expression results support the idea that GBM tumors are more angiogenic than LGG  and even more remarkably in recurrent GBM. These differences were not detected in the peripheral blood, where we observed similar VEGFA levels between HGG and LGG patients. However, significant differences were observed compared to healthy donors, confirming that brain tumor patients release more VEGFA into the bloodstream .
As bevacizumab is an antibody that neutralizes the cytokine VEGFA, we wondered if the benefit from this agent could correlate with its tumor expression pattern. Our results showed that there is a group of cell lines with low VEGFA expression/secretion (GBM27, LN229, 12OCT) that would not be candidates for this anti-angiogenic treatment as no response was observed with standard and specific bevacizumab doses in the cell line that represents that group (LN229). Moreover, VEGFA expression profiles from these cell lines are maintained in in vivo models, pointing to their potential use as pre-clinical models . Focusing on the cell lines that would respond to bevacizumab therapy, VEGFA transcript levels were only significantly reduced using a specific dose, which highlights the need to optimize antibody dosing [40, 41]. In this sense, a recent study has highlighted that adjusting doses of bevacizumab could increase concomitant chemotherapy delivery .
It is known that VEGFA directs binding to its receptor VEGFR2 and stimulates angiogenesis, endothelial cell migration, and cancer cell proliferation and survival [43, 44]. Our study confirms that in U87 and U373 GBM cells treated with the standard dose of bevacizumab, which mimics clinical administration, the migration of brain endothelial cells decreases in vitro and in vivo, with no effect observed in VEGF/VEGFR2 signaling. The inhibition of this switching loop together with the significant reduction in endothelial marker expression and HBMEC migration was only reached when high VEGFA cancer cell lines were treated with the specific dose. As no effect was detected in the LN229 cell line, we speculate that cell lines with a low VEGFA expression/secretion values could correspond with non-responder patients. These data suggest that VEGFA values could be used as a predictor of bevacizumab dose, increasing its therapeutic efficacy.
In this context, one of the challenges is to distinguish between VEGFA released from tumoral cells and other sources . The amount of hVEGFA measured in the serum from our xenograft mouse models corroborates that a specific dose was more effective than the standard dose in pro-angiogenic cell lines. We show that this treatment impairs tumor growth and results in less vasculature. In contrast, none of the bevacizumab dosages had an effect on the cell line with low VEGFA values, confirming our in vitro results. These data have several therapeutic implications, as some patients might not require bevacizumab treatment, while others need an adjusted dosage.
In this regard, several approaches have been studied to identify which patients would benefit from this treatment, for example, VEGFA plasma levels were measured prior to treatment in AVAGLIO cohort , and high VEGFA expression in tissue was correlated with the radiographic response in some patients . However, no association was observed with patient outcomes. In contrast, we and others have shown that VEGFA plasma values do not represent the real angiogenic tumor status . For this reason, as there are no established biomarkers to predict response to bevacizumab therapy, we propose the analysis of VEGFA in the original tissue and serum to help in patient stratification and to identify the mechanisms of resistance and escape routes [49, 50]. Moreover, the analysis of VEGFA in serum could help monitor disease progression and treatment efficacy, supporting MRI data.
Future experiments monitoring VEGFA serum levels with larger cohorts would be required to confirm these effects and to observe differences in the OS, and although our small cohort is highly representative, bevacizumab is not approved by the European Medicines Agency (EMA) . To overcome this issue, we analyzed the OS in the GBM-IDHwt TCGA cohort and divided the cohort based on the VEGFA expression into high and low clusters (Additional file 1: Table S6). Interestingly, when we divided the cohort according to treatment criteria, we did not see differences in OS in the full population (Additional file 1: Figure S11A). However, the median OS of patients with high VEGFA treated with bevacizumab increased in 8.5 months (Additional file 1: Figure S11C), and no effect was observed in patients with low VEGFA expression (Additional file 1: Figure S11B).
Although this is only a proof of concept, since those patients received a multimodal therapy including several drugs, this novel approach could indicate a change in the paradigm of GBM patients’ care.
Our results show that adjusting bevacizumab posology could be a successful strategy to further predict response stratification. Moreover, we propose to measure the amount of VEGFA in serum and not in plasma to enhance treatment efficacy.
Availability of data and materials
Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, et al. The 2016 World Health Organization classification of tumors of the central nervous system: a summary. Acta Neuropathol. 2016;131:803–20.
Ostrom QT, Cioffi G, Gittleman H, Patil N, Waite K, Kruchko C, et al. CBTRUS Statistical Report: primary brain and other central nervous system tumors diagnosed in the United States in 2012-2016. Neuro-Oncology. 2019;21:v1–100.
Stupp R, Taillibert S, Kanner A, Read W, Steinberg DM, Lhermitte B, et al. Effect of tumor-treating fields plus maintenance temozolomide vs maintenance temozolomide alone on survival in patients with glioblastoma a randomized clinical trial. JAMA - J Am Med Assoc. 2017;318:2306–16.
Stupp R, Hegi ME, Gilbert MR, Chakravarti A. Chemoradiotherapy in malignant glioma: standard of care and future directions. J Clin Oncol. 2007;25:4127–36.
Shergalis A, Bankhead A, Luesakul U, Muangsin N, Neamati N, Campus N, et al. Current challenges and opportunities in treating. Pharmacol Rev. 2018;70:412–45.
Krock BL, Skuli N, Simon MC. Hypoxia-induced angiogenesis: good and evil. Genes and Cancer. 2011;2:1117–33.
Angara K, Borin TF, Arbab AS. Vascular mimicry: a novel neovascularization mechanism driving anti-angiogenic therapy (AAT) resistance in glioblastoma. Transl Oncol. 2017;10:650–60.
Goel S, Wong AH, Jain RK. Vascular normalization as a therapeutic strategy. Cold Spring Harb Perspect Med. 2012;2:1–24.
Sullivan LA, Brekken RA. The VEGF family in cancer and antibody-based. 2010;2:165–75.
Mortimer J, Therapeutics E, Zonder H, Pal SK, Therapeutics E. HHS Public Access. Cancer Control. 2012;19:309–16.
Taylor Jennie GE. Anti-angiogenic therapy in high-grade glioma (treatment and toxicity). Curr Treat Options Neurol. 2014;15:328–37.
Hamerlik P, Lathia JD, Rasmussen R, Wu Q, Bartkova J, Lee M, et al. Autocrine VEGF-VEGFR2-Neuropilin-1 signaling promotes glioma stem-like cell viability and tumor growth. J Exp Med. 2012;209:507–20. https://doi.org/10.1084/jem.20111424.
Jakobsen JN, Urup T, Grunnet K, Toft A, Johansen MD, Poulsen SH, et al. Toxicity and efficacy of lomustine and bevacizumab in recurrent glioblastoma patients. J Neuro-Oncol. 2018;137:439–46.
Herrlinger U, Tzaridis T, Mack F, Steinbach JP, Schlegel U, Sabel M, et al. Lomustine-temozolomide combination therapy versus standard temozolomide therapy in patients with newly diagnosed glioblastoma with methylated MGMT promoter (CeTeG/NOA–09): a randomised, open-label, phase 3 trial. Lancet. 2019;393:678–88.
Chinot OL, Wick W, Mason W, Henriksson R, Saran F, Nishikawa R, et al. Bevacizumab plus radiotherapy-temozolomide for newly diagnosed glioblastoma. N Engl J Med. 2014;370:709–22.
Gilbert MR, Dignam JJ, Armstrong TS, Wefel JS, Blumenthal DT, Vogelbaum MA, et al. A randomized trial of bevacizumab for newly diagnosed glioblastoma. N Engl J Med. 2014;370:699–708.
Chinot OL, De La Motte RT, Moore N, Zeaiter A, Das A, Phillips H, et al. AVAglio: phase 3 trial of bevacizumab plus temozolomide and radiotherapy in newly diagnosed glioblastoma multiforme. Adv Ther. 2011;28:334–40.
Friedman HS, Prados MD, Wen PY, Mikkelsen T, Schiff D, Abrey LE, et al. Bevacizumab alone and in combination with irinotecan in recurrent glioblastoma. J Clin Oncol. 2009;27:4733–40.
Kreisl TN, Kim L, Moore K, Duic P, Royce C, Stroud I, et al. Phase II trial of single-agent bevacizumab followed by bevacizumab plus irinotecan at tumor progression in recurrent glioblastoma. J Clin Oncol. 2009;27:740–5.
Ameratunga M, Pavlakis N, Wheeler H, Grant R, Simes J, Khasraw M. Anti-angiogenic therapy for high-grade glioma (review); 2018.
Goyal P, Tenenbaum M, Gupta S, Kochar PS, Bhatt AA, Mangla M, et al. Survival prediction based on qualitative MRI diffusion signature in patients with recurrent high grade glioma treated with bevacizumab. Quant Imaging Med Surg. 2018;8:268–79.
Boisen MK, Holst CB, Consalvo N, Chinot OL, Johansen JS. Plasma YKL-40 as a biomarker for bevacizumab efficacy in patients with newly diagnosed glioblastoma in the phase 3 randomized AVAglio trial. Oncotarget. 2018;9:6752–62. https://doi.org/10.18632/oncotarget.22886.
Sandmann T, Bourgon R, Garcia J, Li C, Cloughesy T, Chinot OL, et al. Patients with proneural glioblastoma may derive overall survival benefit from the addition of bevacizumab to first-line radiotherapy and temozolomide: retrospective analysis of the AVAglio trial. J Clin Oncol. 2015;33:2735–44.
Zahonero C, Aguilera P, Ramírez-Castillejo C, Pajares M, Bolós MV, Cantero D, et al. Preclinical test of dacomitinib, an irreversible EGFR inhibitor, confirms its effectiveness for glioblastoma. Mol Cancer Ther. 2015;14:1548–58.
Stins MF, Badger J, Sik KK. Bacterial invasion and transcytosis in transfected human brain microvascular endothelial cells. Microb Pathog. 2001;30:19–28.
García-Romero N, González-Tejedo C, Carrión-Navarro J, Esteban-Rubio S, Rackov G, Rodríguez-Fanjul V, et al. Cancer stem cells from human glioblastoma resemble but do not mimic original tumors after in vitro passaging in serum-free media. Oncotarget. 2016;7:65888–90.
Von Baumgarten L, Brucker D, Tirniceru A, Kienast Y, Grau S, Burgold S, et al. Bevacizumab has differential and dose-dependent effects on glioma blood vessels and tumor cells. Clin Cancer Res. 2011;17:6192–205.
Ceccarelli M, Barthel FP, Malta TM, Sabedot TS, Salama SR, et al. Molecular profiling reveals biologically discrete subsets and pathways of progression in diffuse glioma. Cell. 2016;164(3):550–63.
Sumner G, Georgaros C, Rafique A, Dicioccio T, Martin J, Papadopoulos N, et al. Anti-VEGF drug interference with VEGF quantitation in the R&D systems human quantikine VEGF ELISA kit. Bioanalysis. 2019;11:381–92.
Arcondéguy T, Lacazette E, Millevoi S, Prats H, Touriol C. VEGF-A mRNA processing, stability and translation: a paradigm for intricate regulation of gene expression at the post-transcriptional level. Nucleic Acids Res. 2013;41:7997–8010.
Chi AS, Sorensen AG, Jain RK, Batchelor TT. Angiogenesis as a therapeutic target in malignant gliomas. Oncologist. 2009;14:621–36.
Norton K, Popel AS. Effects of endothelial cell proliferation and migration rates in a computational model of sprouting angiogenesis. Nat Publ Gr. 2016;6:36992.
Rosen L. VEGF-targeted therapy: therapeutic potential and recent advances 2005;10:382–391.
Wang N, Jain RK, Batchelor TT. New directions in anti-angiogenic therapy for glioblastoma; 2017. p. 321–32.
Hovey EJ, Hashemi-sadraei N, Ahluwalia MS. Bevacizumab in high-grade gliomas: a review of its uses, toxicity assessment, and future treatment challenges. 2013;:371–389.
Field KM, Jordan JT, Wen PY, Rosenthal MA, Reardon DA. Bevacizumab and glioblastoma: scientific review, newly reported updates, and ongoing controversies; 2015. p. 997–1007.
Plate KH, Breier G, Weich HA, Risau W. Vascular endothelial growth factor is a potential tumour angiogenesis factor in human gliomas in vivo. Nature. 1992;359:845–8.
Nowacka A. Serum VEGF-A concentrations in patients with central nervous system (CNS) tumors. Plos One. 2018;1:1–8.
Lenting K, Verhaak R, Pieter L, Leenders W. Glioma: experimental models and reality. Acta Neuropathol. 2017;133:263–82.
Falk AT, Barrière J, Franc E, Follana P. Bevacizumab: a dose review; 2015.
Ajlan A, Thomas P, Albakr A, Nagpal S. Optimizing bevacizumab dosing in glioblastoma: less is more. J Neuro-Oncol 2017;135:99–105.
Gerstner E, Emblem KE, Chang K, Vakulenko-Lagun B, Yen Y-F, Beers AL, et al. Bevacizumab reduces permeability and concurrent temozolomide delivery in a subset of patients with recurrent glioblastoma. Clin Cancer Res. 2019;:clincanres.1739.2019. doi:https://doi.org/10.1158/1078-0432.CCR-19-1739.
Suzuki Y, Montagne K, Nishihara A, Watabe T. BMPs promote proliferation and migration of endothelial cells via stimulation of VEGF-A/VEGFR2 and angiopoietin-1/Tie2 signalling. J Biochem. 2008;206:199–206.
Oommen S, Gupta SKVN. Vascular endothelial growth factor A (VEGF-A) induces endothelial and cancer cell migration through direct binding to integrin α9β1. J Biol Chem. 2011;286:1083–92.
Ferreira G, Galera BB, Duarte EC, Chen ES, Azis L, Damazo AS, et al. Prognostic significance of vascular endothelial growth factor polymorphisms in colorectal cancer patients. World J Gastrointest Oncol. 2017;9:78–86.
Nishikawa R, Saran F, Mason W, Wick W, Cloughesy TF, Henriksson R, et al. Biomarker (BM) evaluations in the phase III AVAglio study of bevacizumab (Bv) plus standard radiotherapy (RT) and temozolomide (T) for newly diagnosed glioblastoma (GBM). J Clin Oncol. 2013;31(15_suppl):2023.
Sathornsumetee S, Cao Y, Marcello JE, Ii JEH, Roger E, Desjardins A, et al. NIH Public Access. J Clin Oncol. 2014;26:271–8.
Hong Y, Hwang D, Hong S. Clinical usefulness of serum and plasma vascular endothelial growth factor in cancer patients: Which is the optimal specimen? 2000. p. 149–52.
Takano S, Ishikawa E, Matsuda M, Sakamoto N. The anti-angiogenic role of soluble-form VEGF receptor in malignant gliomas; 2016. p. 515–24.
de Groot JF, Piao Y, Tran H, Gilbert M, Wu HK, Liu J, et al. Myeloid biomarkers associated with glioblastoma response to anti-VEGF therapy with aflibercept. Clin Cancer Res. 2011;(14):4872–81.
Yu Z, Zhao G, Zhang Z, Li Y, Chen Y, Wang NAN, et al. Efficacy and safety of bevacizumab for the treatment of glioblastoma (review). Exp Ther Med. 2016;11(2):371–80.
We would like to acknowledge Julia Palacios García for the migration technical support. We also thank Victor González Rumayor for helping with the tissue microarray logistic.
We want to particularly acknowledge the patients in this study for their participation and the HGM BioBank integrated in RETICS, National Network Biobanks, and collaborating Centers for the generous gifts of clinical samples used in this work. The HGM BioBank integrated in National Network Biobanks is supported by the Instituto de Salud Carlos III, Spanish Science and Innovation Ministry (Grant no. PT17/0015/0042).
This work was supported by grants from the “Fondo de Investigaciones Sanitarias” (FIS) (PI17-01489), the Miguel Servet Program (CP11/00147) del Instituto de Salud Carlos III (AAS), and the Ministerio de Economía y Competitividad–FEDERER (RTC-2016-4990-1). IPA was supported by “Ayudas para la contratación de ayudantes de investigación cofinanciadas por el Fondo Social Europeo a través del Programa Operativo de Empleo Juvenil y la Iniciativa de Empleo Juvenil (YEI),” and SER was supported by FPI-CEU predoctoral fellowship.
Ethics approval and consent to participate
Procedures used on mice were approved by and performed according to the guidelines of the institutional animal care committee of the Instituto Carlos III in agreement with the European Union and National directives.
Samples from patients were kindly provided by the HGM BioBank, integrated in RETICS, National Network Biobanks, funding by Instituto de Salud Carlos III, from Biobanco La Fe (PT17/0015/0043), Valencia, Spain, and HM Hospitales, Madrid, Spain. Samples were processed following the current procedures and frozen immediately after their reception. All patients participating in the study gave their informed consent, and protocols were approved by institutional ethical committees.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
García-Romero, N., Palacín-Aliana, I., Madurga, R. et al. Bevacizumab dose adjustment to improve clinical outcomes of glioblastoma. BMC Med 18, 142 (2020). https://doi.org/10.1186/s12916-020-01610-0