Transcriptional changes induced by bevacizumab combination therapy in responding and non-responding recurrent glioblastoma patients

Early 2-Factor Transcription Factors Associated with Progression and Recurrence in Bevacizumab-Responsive Subtypes of Glioblastoma

The early 2-factor (E2F) family of transcription factors, including E2F1 through 8, plays a critical role in apoptosis, metabolism, proliferation, and angiogenesis within glioblastoma (GBM). However, the specific functions of E2F transcription factors (E2Fs) and their impact on the malignancy of Bevacizumab (BVZ)-responsive GBM subtypes remain unclear. This study used data from The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), European Molecular Biology Laboratory's European Bioinformatics Institute (EMBL-EBI), and Gene Expression Omnibus (GEO) to explore the impact of eight E2F family members on the clinical characteristics of BVZ-responsive GBM subtypes and possible mechanisms of recurrence after BVZ treatment. Using machine learning algorithms, including TreeBagger and deep neural networks, we systematically predicted and validated GBM patient survival terms based on the expression profiles of E2Fs across BVZ-responsive GBM subtypes. Our bioinformatics analyses suggested that a significant increase in E2F8 post-BVZ treatment may enhance the function of angiogenesis and stem cell proliferation, implicating this factor as a candidate mechanism of GBM recurrence after treatment. In addition, BVZ treatment in unresponsive GBM patients may potentially worsen disease progression. These insights underscore that E2F family members play important roles in GBM malignancy and BVZ treatment response, highlighting their potential as prognostic biomarkers, therapeutic targets, and recommending precision BVZ treatment to individual GBM patients.

Single-cell RNA-seq reveals T cell exhaustion and immune response landscape in osteosarcoma

Background

T cell exhaustion in the tumor microenvironment has been demonstrated as a substantial contributor to tumor immunosuppression and progression. However, the correlation between T cell exhaustion and osteosarcoma (OS) remains unclear.

Methods

In our present study, single-cell RNA-seq data for OS from the GEO database was analysed to identify CD8+ T cells and discern CD8+ T cell subsets objectively. Subgroup differentiation trajectory was then used to pinpoint genes altered in response to T cell exhaustion. Subsequently, six machine learning algorithms were applied to develop a prognostic model linked with T cell exhaustion. This model was subsequently validated in the TARGETs and Meta cohorts. Finally, we examined disparities in immune cell infiltration, immune checkpoints, immune-related pathways, and the efficacy of immunotherapy between high and low TEX score groups.

Results

The findings unveiled differential exhaustion in CD8+ T cells within the OS microenvironment. Three genes related to T cell exhaustion (RAD23A, SAC3D1, PSIP1) were identified and employed to formulate a T cell exhaustion model. This model exhibited robust predictive capabilities for OS prognosis, with patients in the low TEX score group demonstrating a more favorable prognosis, increased immune cell infiltration, and heightened responsiveness to treatment compared to those in the high TEX score group.

Conclusion

In summary, our research elucidates the role of T cell exhaustion in the immunotherapy and progression of OS, the prognostic model constructed based on T cell exhaustion-related genes holds promise as a potential method for prognostication in the management and treatment of OS patients.

Overexpression of CD99 is associated with tumor adaptiveness and indicates the tumor recurrence and therapeutic responses in gliomas

Glioma undergoes adaptive changes, leading to poor prognosis and resistance to treatment. CD99 influences the migration and invasion of glioma cells and plays an oncogene role. However, whether CD99 can affect the adaptiveness of gliomas is still lacking in research, making its clinical value underestimated. Here, we enrolled our in-house and public multiomics datasets for bioinformatic analysis and conducted immunohistochemistry staining to investigate the role of CD99 in glioma adaptive response and its clinical implications. CD99 is expressed in more adaptative glioma subtypes and cell states. Under hypoxic conditions, CD99 is upregulated in glioma cells and is associated with angiogenesis and metabolic adaptations. Gliomas with over-expressed CD99 also increased the immunosuppressive tumor-associated macrophages. The relevance with tumor adaptiveness of CD99 presented clinical significance. We discovered that CD99 overexpression is associated with short-time recurrence and validated its prognostic value. Additionally, Glioma patients with high expression of CD99 were resistant to chemotherapy and radiotherapy. The CD99 expression was also related to anti-angiogenic and immune checkpoint inhibitor therapy response. Inhibitors of the PI3K-AKT pathway have therapeutic potential against CD99-overexpressing gliomas. Our study identified CD99 as a biomarker characterizing the adaptive response in glioma. Gliomas with high CD99 expression are highly tolerant to stress conditions such as hypoxia and antitumor immunity, making treatment responses dimmer and tumor progression. Therefore, for patients with CD99-overexpressing gliomas, tumor adaptiveness should be fully considered during treatment to avoid drug resistance, and closer clinical monitoring should be carried out to improve the prognosis.

Polymorphism at codon 31 of CDKN1A (p21) as a predictive factor for bevacizumab therapy in glioblastoma multiforme

Glioblastoma (GBM), a prevalent and malignant brain tumor, poses a challenge in surgical resection due to its invasive nature within the brain parenchyma. CDKN1A (p21, Waf-1), a cyclin-dependent kinase inhibitor, plays a pivotal role in regulating cell growth arrest, terminal differentiation, and apoptosis. The existence of natural variants of CDKN1A has been associated with specific cancer types. In this retrospective study, our objective was to identify polymorphic variants of CDKN1A, specifically c.93C > A (codon 31 Ser31Arg), and investigate its potential impact within the scope of bevacizumab therapy for glioblastoma multiforme. This study involved a cohort of 139 unrelated adult Chinese GBM patients in Taiwan. Genomic DNA extracted from tumor samples was utilized for genotyping using the polymerase chain reaction (PCR) restriction fragment length polymorphism method (PCR-RFLP analysis). Through unconditional logistic regression analysis, odds ratios (ORs) with corresponding 95% confidence intervals (CIs) were calculated. Our findings unveiled that among these GBM patients, the distribution of codon 31 polymorphisms was as follows: 23.02% were Serine homozygotes (Ser/Ser), 27.34% were Arginine homozygotes (Arg/Arg), and 49.64% were Serine/Arginine heterozygotes (Ser/Arg). While CDKN1A c.93C > A polymorphisms did not exhibit a direct association with overall survival in GBM patients, noteworthy survival benefits emerged among individuals with Arg/Arg and Arg/Ser genotypes who received combined concurrent chemoradiotherapy (CCRT) and bevacizumab treatment compared to those who underwent CCRT alone. Our findings indicate a significant involvement of the CDKN1A c.93C > A polymorphism in the development and onset of GBM, offering potential implications for the early prognostication of bevacizumab therapy outcomes.

A programmed cell death-related model based on machine learning for predicting prognosis and immunotherapy responses in patients with lung adenocarcinoma

Background

lung adenocarcinoma (LUAD) remains one of the most common and lethal malignancies with poor prognosis. Programmed cell death (PCD) is an evolutionarily conserved cell suicide process that regulates tumorigenesis, progression, and metastasis of cancer cells. However, a comprehensive analysis of the role of PCD in LUAD is still unavailable.

Methods

We analyzed multi-omic variations in PCD-related genes (PCDRGs) for LUAD. We used cross-validation of 10 machine learning algorithms (101 combinations) to synthetically develop and validate an optimal prognostic cell death score (CDS) model based on the PCDRGs expression profile. Patients were classified based on their median CDS values into the high and low-CDS groups. Next, we compared the differences in the genomics, biological functions, and tumor microenvironment of patients between both groups. In addition, we assessed the ability of CDS for predicting the response of patients from the immunotherapy cohort to immunotherapy. Finally, functional validation of key genes in CDS was performed.

Results

We constructed CDS based on four PCDRGs, which could effectively and consistently stratify patients with LUAD (patients with high CDS had poor prognoses). The performance of our CDS was superior compared to 77 LUAD signatures that have been previously published. The results revealed significant genetic alterations like mutation count, TMB, and CNV were observed in patients with high CDS. Furthermore, we observed an association of CDS with immune cell infiltration, microsatellite instability, SNV neoantigens. The immune status of patients with low CDS was more active. In addition, CDS could be reliable to predict therapeutic response in multiple immunotherapy cohorts. In vitro experiments revealed that high DNA damage inducible transcript 4 (DDIT4) expression in LUAD cells mediated protumor effects.

Conclusion

CDS was constructed based on PCDRGs using machine learning. This model could accurately predict patients' prognoses and their responses to therapy. These results provide new promising tools for clinical management and aid in designing personalized treatment strategies for patients with LUAD.

Circulating galectin-1 delineates response to bevacizumab in melanoma patients and reprograms endothelial cell biology

Blockade of vascular endothelial growth factor (VEGF) signaling with bevacizumab, a humanized anti-VEGF monoclonal antibody (mAb), or with receptor tyrosine kinase inhibitors, has improved progression-free survival and, in some indications, overall survival across several types of cancers by interrupting tumor angiogenesis. However, the clinical benefit conferred by these therapies is variable, and tumors from treated patients eventually reinitiate growth. Previously we demonstrated, in mouse tumor models, that galectin-1 (Gal1), an endogenous glycan-binding protein, preserves angiogenesis in anti-VEGF-resistant tumors by co-opting the VEGF receptor (VEGFR)2 signaling pathway in the absence of VEGF. However, the relevance of these findings in clinical settings is uncertain. Here, we explored, in a cohort of melanoma patients from AVAST-M, a multicenter, open-label, randomized controlled phase 3 trial of adjuvant bevacizumab versus standard surveillance, the role of circulating plasma Gal1 as part of a compensatory mechanism that orchestrates endothelial cell programs in bevacizumab-treated melanoma patients. We found that increasing Gal1 levels over time in patients in the bevacizumab arm, but not in the observation arm, significantly increased their risks of recurrence and death. Remarkably, plasma Gal1 was functionally active as it was able to reprogram endothelial cell biology, promoting migration, tubulogenesis, and VEGFR2 phosphorylation. These effects were prevented by blockade of Gal1 using a newly developed fully human anti-Gal1 neutralizing mAb. Thus, using samples from a large-scale clinical trial from stage II and III melanoma patients, we validated the clinical relevance of Gal1 as a potential mechanism of resistance to bevacizumab treatment.

ITGA5 Predicts Dual-Drug Resistance to Temozolomide and Bevacizumab in Glioma

AIMS

Anti-angiotherapy (Bevacizumab) is currently regarded as a promising option for glioma patients who are resistant to temozolomide (TMZ) treatment. But ongoing clinical research failed to meet therapeutic expectations. This study aimed to explore the pivotal genetic feature responsible for TMZ and Bevacizumab resistance in glioma patients.

METHODS

We downloaded the transcriptomic and methylation data of glioma patients from The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), and Gene Expression Omnibus (GEO) databases and grouped these patients into resistant and non-resistant groups based on their clinical profiles. Differentially expressed genes and pathways were identified and exhibited with software in R platform. A TMZ-resistant cell line was constructed for validating the expression change of the candidate gene, ITGA5. An ITGA5-overexpressing cell line was also constructed to investigate its biological function using the CCK8 assay, Western blot, periodic acid-Schiff (PAS) staining, and transcriptional sequencing.

RESULTS

Change of the cell morphology and polarity was closely associated with TMZ mono-resistance and TMZ/Bevacizumab dual resistance in glioma patients. The expression level of ITGA5 was effective in determining drug resistance and the outcome of glioma patients, which is regulated by methylation on two distinct sites. ITGA5 was augmented in TMZ-resistant glioma cells, while overexpressing ITGA5 altered the cell-promoted TMZ resistance through enhancing vascular mimicry (VM) formation correspondingly.

CONCLUSIONS

Both the epigenetic and transcriptional levels of ITGA5 are effective in predicting TMZ and Bevacizumab resistance, indicating that ITGA5 may serve as a predictor of the treatment outcomes of glioma patients.

A retrospective analysis of GSE84010: Cell adhesion molecules might contribute to bevacizumab resistance in glioblastoma

Bevacizumab (BEV) is an anti-angiogenesis antibody which has shown favorable therapeutic effects on some solid tumors. However, many clinical trials showed that BEV could only improve PFS instead of OS in glioblastoma (GBM) patients. However, some studies indicate that specific molecular subtypes of GBM could still benefit from combination treatment of BEV and Stupp protocol. Through the subgroup analysis of GSE84010 dataset, we found the neural and proneural subgroup can benefit from the administration of BEV in terms of OS, which is statistically significant. The further KEGG pathway enrichment analysis showed cell adhesion molecules (CAMs) pathway was enriched, and the expression of ITGAM has a predictive value for prognosis. These findings can provide some hints for future administration of BEV in newly diagnosed GBM patients.

Perspective: targeting VEGF-A and YKL-40 in glioblastoma - matter matters

Glioblastomas (GBM) are heterogeneous highly vascular brain tumors exploiting the unique microenvironment in the brain to resist treatment and anti-tumor responses. Anti-angiogenic agents, immunotherapy, and targeted therapy have been studied extensively in GBM patients over a number of decades with minimal success. Despite maximal efforts, prognosis remains dismal with an overall survival of approximately 15 months.

Bevacizumab, a humanized anti-vascular endothelial growth factor (VEGF) antibody, underwent accelerated approval by the U.S. Food and Drug Administration in 2009 for the treatment of recurrent GBM based on promising preclinical and early clinical studies. Unfortunately, subsequent clinical trials did not find overall survival benefit. Pursuing pleiotropic targets and leaning toward multitarget strategies may be a key to more effective therapeutic intervention in GBM, but preclinical evaluation requires careful consideration of model choices. In this study, we discuss bevacizumab resistance, dual targeting of pro-angiogenic modulators VEGF and YKL-40 in the context of brain tumor microenvironment, and how model choice impacts study conclusions and its translational significance.

A Six-lncRNA Signature for Immunophenotype Prediction of Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is the most aggressive primary tumor of the central nervous system. As biomedicine advances, the researcher has found the development of GBM is closely related to immunity. In this study, we evaluated the GBM tumor immunoreactivity and defined the Immune-High (IH) and Immune-Low (IL) immunophenotypes using transcriptome data from 144 tumors profiled by The Cancer Genome Atlas (TCGA) project based on the single-sample gene set enrichment analysis (ssGSEA) of five immune expression signatures (IFN-γ response, macrophages, lymphocyte infiltration, TGF-β response, and wound healing). Next, we identified six immunophenotype-related long non-coding RNA biomarkers (im-lncRNAs, USP30-AS1, HCP5, PSMB8-AS1, AL133264.2, LINC01684, and LINC01506) by employing a machine learning computational framework combining minimum redundancy maximum relevance algorithm (mRMR) and random forest model. Moreover, the expression level of identified im-lncRNAs was converted into an im-lncScore using the normalized principal component analysis. The im-lncScore showed a promising performance for distinguishing the GBM immunophenotypes with an area under the curve (AUC) of 0.928. Furthermore, the im-lncRNAs were also closely associated with the levels of tumor immune cell infiltration in GBM. In summary, the im-lncRNA signature had important clinical implications for tumor immunophenotyping and guiding immunotherapy in glioblastoma patients in future.

Increased epithelial membrane protein 2 expression in glioblastoma after treatment with bevacizumab

Background

Antiangiogenic therapy with bevacizumab has failed to provide substantial gains in overall survival. Epithelial membrane protein 2 (EMP2) is a cell surface protein that has been previously shown to be expressed in glioblastoma, correlate with poor survival, and regulate neoangiogenesis in cell lines. Thus, the relationship between bevacizumab and EMP2 was investigated.

Methods

Tumor samples were obtained from 12 patients with newly diagnosed glioblastoma at 2 time points: (1) during the initial surgery and (2) during a subsequent surgery following disease recurrence post-bevacizumab treatment. Clinical characteristics and survival data from these patients were collected, and tumor samples were stained for EMP2 expression. The IVY Glioblastoma Atlas Project database was used to evaluate EMP2 expression levels in 270 samples by differing histological areas of the tumor.

Results

Patients with high EMP2 staining at initial diagnosis had decreased progression-free and overall survival after bevacizumab (median progression-free survival 4.6 months vs 5.9 months; log-rank P = .076 and overall survival 7.7 months vs 14.4 months; log-rank P = .011). There was increased EMP2 staining in samples obtained after bevacizumab treatment in both unpaired (mean H-score 2.31 vs 1.76; P = .006) and paired analyses (mean difference 0.571; P = .019). This expression increase correlated with length of bevacizumab therapy (R ²  = 0.449; Pearson P = .024).

Conclusions

Bevacizumab treatment increased EMP2 protein expression. This increase in EMP2 correlated with reduced mean survival time post-bevacizumab therapy. We hypothesize a role of EMP2 in clinical bevacizumab resistance and as a potential antiangiogenic therapeutic target in glioblastoma.

Angiotensinogen promoter methylation predicts bevacizumab treatment response of patients with recurrent glioblastoma

Patients with recurrent glioblastoma achieving response to bevacizumab combined with chemotherapy have clinical improvement and prolonged survival. High gene expression of angiotensinogen (AGT) is associated with a poor bevacizumab response. Because AGT expression is epigenetically regulated, we aimed to investigate whether AGT promoter methylation in tumor tissue predicts response to bevacizumab combination therapy in patients with recurrent glioblastoma. The study included 159 patients with recurrent glioblastoma, treated with bevacizumab combination treatment (training cohort, n = 77; validation cohort, n = 82). All patients could be evaluated for treatment response and biomarkers. DNA methylation of 4 CpG sites in the AGT promoter was measured using pyrosequencing. A model for nonresponse was established using logistic regression analysis. In the training cohort, lower methylation of each of the four CpG sites in the AGT promoter was significantly associated with nonresponse (all P < 0.05). Moreover, the mean methylation level of all four CpG sites was associated with an increased likelihood of not achieving response to bevacizumab combination therapy (twofold decrease: odds ratio = 3.01; 95% confidence interval: 1.41-6.44; P = 0.004). We developed a model for nonresponse in the training cohort, where a threshold of mean AGT promoter methylation levels was set to below 12%. The model could predict bevacizumab nonresponse with 96% specificity. Importantly, this predictor was also significantly associated with nonresponse in the validation cohort (P = 0.037). Taken together, our findings suggest that low AGT promoter methylation in tumor tissue predicts nonresponse to bevacizumab combination treatment in patients with recurrent glioblastoma. We have, thus, established and successfully validated a predictor for nonresponse that can be used to identify patients who will not benefit from bevacizumab combination therapy.

The Brain Penetrating and Dual TORC1/TORC2 Inhibitor, RES529, Elicits Anti-Glioma Activity and Enhances the Therapeutic Effects of Anti-Angiogenetic Compounds in Preclinical Murine Models

Background. Glioblastoma multiforme (GBM) is a devastating disease showing a very poor prognosis. New therapeutic approaches are needed to improve survival and quality of life. GBM is a highly vascularized tumor and as such, chemotherapy and anti-angiogenic drugs have been combined for treatment. However, as treatment-induced resistance often develops, our goal was to identify and treat pathways involved in resistance to treatment to optimize the treatment strategies. Anti-angiogenetic compounds tested in preclinical and clinical settings demonstrated recurrence associated to secondary activation of the phosphatidylinositol 3-kinase (PI3K)/AKT/mTOR pathway. Aims. Here, we determined the sensitizing effects of the small molecule and oral available dual TORC1/TORC2 dissociative inhibitor, RES529, alone or in combination with the anti-VEGF blocking antibody, bevacizumab, or the tyrosine kinase inhibitor, sunitinib, in human GBM models. Results. We observed that RES529 effectively inhibited dose-dependently the growth of GBM cells in vitro counteracting the insurgence of recurrence after bevacizumab or sunitinib administration in vivo. Combination strategies were associated with reduced tumor progression as indicated by the analysis of Time to Tumor Progression (TTP) and disease-free survival (DSF) as well as increased overall survival (OS) of tumor bearing mice. RES529 was able to reduce the in vitro migration of tumor cells and tubule formation from both brain-derived endothelial cells (angiogenesis) and tumor cells (vasculogenic mimicry). Conclusions. In summary, RES529, the first dual TORC1/TORC2 dissociative inhibitor, lacking affinity for ABCB1/ABCG2 and having good brain penetration, was active in GBM preclinical/murine models giving credence to its use in clinical trial for patients with GBM treated in association with anti-angiogenetic compounds.

Hypoxia Regulated Gene Network in Glioblastoma Has Special Algebraic Topology Structures and Revealed Communications Involving Warburg Effect and Immune Regulation

Hypoxia regulated genes (HRGs) formed a complex molecular interaction network (MINW), contributing to many aspects of glioblastoma (GBM) tumor biology. However, little is known about the intrinsic structures of the HRGs-MINW, mainly due to a lack of analysis tools to decipher MINWs. By introducing general hyper-geometric distribution, we obtained a statistically reliable gene set of HRGs (SR-HRGs) from several datasets. Next, MINWs were reconstructed from several independent GBM expression datasets. Algebraic topological analysis was performed to quantitatively analyze the amount of equivalence classes of cycles in various dimensions by calculating the Betti numbers. Persistent homology analysis of a filtration of growing networks was further performed to examine robust topological structures in the network by investigating the Betti curves, life length of the cycles. Random networks with the same number of node and edge and degree distribution were produced as controls. As a result, GBM-HRGs-MINWs reconstructed from different datasets exhibited great consistent Betti curves to each other, which were significantly different from that of random networks. Furthermore, HRGs-MINWs reconstructed from normal brain expression datasets exhibited topological structures significantly different from that of GBM-HRGs-MINWs. Analysis of cycles in GBM-HRGs-MINWs revealed genes that had clinical implications, and key parts of the cycles were also identified in reconstructed protein-protein interaction networks. In addition, the cycles are composed by genes involved in the Warburg effect, immune regulation, and angiogenesis. In summary, GBM-HRGs-MINWs contained abundant molecular interacting cycles in different dimensions, which are composed by genes involved in multiple programs essential for the tumorigenesis of GBM, revealing novel interaction diagrams in GBM and providing novel potential therapeutic targets.

Neoadjuvant anti-PD-1 immunotherapy promotes a survival benefit with intratumoral and systemic immune responses in recurrent glioblastoma

Glioblastoma is the most common primary malignant brain tumor in adults and is associated with poor survival. The Ivy Foundation Early Phase Clinical Trials Consortium conducted a randomized, multi-institution clinical trial to evaluate immune responses and survival following neoadjuvant and/or adjuvant therapy with pembrolizumab in 35 patients with recurrent, surgically resectable glioblastoma. Patients who were randomized to receive neoadjuvant pembrolizumab, with continued adjuvant therapy following surgery, had significantly extended overall survival compared to patients that were randomized to receive adjuvant, post-surgical programmed cell death protein 1 (PD-1) blockade alone. Neoadjuvant PD-1 blockade was associated with upregulation of T cell- and interferon-γ-related gene expression, but downregulation of cell-cycle-related gene expression within the tumor, which was not seen in patients that received adjuvant therapy alone. Focal induction of programmed death-ligand 1 in the tumor microenvironment, enhanced clonal expansion of T cells, decreased PD-1 expression on peripheral blood T cells and a decreasing monocytic population was observed more frequently in the neoadjuvant group than in patients treated only in the adjuvant setting. These findings suggest that the neoadjuvant administration of PD-1 blockade enhances both the local and systemic antitumor immune response and may represent a more efficacious approach to the treatment of this uniformly lethal brain tumor.

Outcome of Bevacizumab Therapy in Patients with Recurrent Glioblastoma Treated with Angiotensin System Inhibitors

PURPOSE

Antihypertensive therapy may improve bevacizumab efficacy in cancer patients. We examined efficacy and toxicity of angiotensin system inhibitors (ASI) and other antihypertensive drugs in bevacizumab treated recurrent glioblastoma patients.

METHODS

We retrospectively combined a national prescription registry with a clinical database with recurrent glioblastoma patients (n = 243).

RESULTS

Patients who initiated ASI after bevacizumab (n = 26) showed a tendency towards improved progression-free survival and overall survival (OS) with hazard rate (HR) reductions (HR = 0.70 and HR = 0.79, respectively). Calcium antagonists during bevacizumab therapy significantly improved OS (HR = 0.57).

CONCLUSIONS

Overall the study supports a potential beneficial effect of antihypertensive treatment on prognosis of bevacizumab treated glioblastoma patients.

Glioblastoma endothelium drives bevacizumab-induced infiltrative growth via modulation of PLXDC1

Bevacizumab, a VEGF‐targeting monoclonal antibody, may trigger an infiltrative growth pattern in glioblastoma. We investigated this pattern using both a human specimen and rat models. In the human specimen, a substantial fraction of infiltrating tumor cells were located along perivascular spaces in close relationship with endothelial cells. Brain xenografts of U87MG cells treated with bevacizumab were smaller than controls (p = 0.0055; Student t‐test), however, bands of tumor cells spread through the brain farther than controls (p < 0.001; Student t‐test). Infiltrating tumor Cells exhibited tropism for vascular structures and propensity to form tubules and niches with endothelial cells. Molecularly, bevacizumab triggered an epithelial to mesenchymal transition with over‐expression of the receptor Plexin Domain Containing 1 (PLXDC1). These results were validated using brain xenografts of patient‐derived glioma stem‐like cells. Enforced expression of PLXDC1 in U87MG cells promoted brain infiltration along perivascular spaces. Importantly, PLXDC1 inhibition prevented perivascular infiltration and significantly increased the survival of bevacizumab‐treated rats. Our study indicates that bevacizumab‐induced brain infiltration is driven by vascular endothelium and depends on PLXDC1 activation of tumor cells.

What's new?

Bevacizumab, a VEGF‐targeting monoclonal antibody, has been observed to trigger an infiltrative growth pattern in glioblastoma as an escape mechanism. The mechanisms underlying this gliomatosis‐like growth pattern, however, remain unclear. Here, the authors found that the infiltrative growth pattern occurs mostly along perivascular spaces and relies on the over‐expression of PLXDC1 by tumor cells and on the restoration of the endothelial component of blood brain barrier. Altogether, the data show that the brain infiltration induced by bevacizumab is mainly driven by the vascular endothelium. Importantly, inhibition of PLXDC1 prevents bevacizumab‐induced infiltrative growth, resulting in significant increase of survival.