Lack of functional normalisation of tumour vessels following anti-angiogenic therapy in glioblastoma

Metabolic changes and anti-tumor effects of a ketogenic diet combined with anti-angiogenic therapy in a glioblastoma mouse model

The ketogenic diet (KD) is a high fat and low carbohydrate diet that produces ketone bodies through imitation of starvation. The combination of KD and Bevacizumab (Bev), a VEGF inhibitor, is considered to further reduce the supply of glucose to the tumor. The metabolite changes in U87 glioblastoma mouse models treated with KD and/or Bev were examined using gas chromatography-mass spectrometry. The combination therapy of KD and Bev showed a decrease in the rate of tumor growth and an increase in the survival time of mice, although KD alone did not have survival benefit. In the metabolome analysis, the pattern of changes for most amino acids are similar between tumor and brain tissues, however, some amino acids such as aspartic acid and glutamic acid were different between tumors and brain tissues. The KD enhanced the anti-tumor efficacy of Bev in a glioblastoma intracranial implantation mouse model, based on lowest levels of microvascular density (CD31) and cellular proliferation markers (Ki-67 and CCND1) in KD + Bev tumors compared to the other groups. These results suggested that KD combined with Bev may be a useful treatment strategy for patients with GBM.

Patient-derived organoids and orthotopic xenografts of primary and recurrent gliomas represent relevant patient avatars for precision oncology

Patient-based cancer models are essential tools for studying tumor biology and for the assessment of drug responses in a translational context. We report the establishment a large cohort of unique organoids and patient-derived orthotopic xenografts (PDOX) of various glioma subtypes, including gliomas with mutations in IDH1, and paired longitudinal PDOX from primary and recurrent tumors of the same patient. We show that glioma PDOXs enable long-term propagation of patient tumors and represent clinically relevant patient avatars that retain histopathological, genetic, epigenetic, and transcriptomic features of parental tumors. We find no evidence of mouse-specific clonal evolution in glioma PDOXs. Our cohort captures individual molecular genotypes for precision medicine including mutations in IDH1, ATRX, TP53, MDM2/4, amplification of EGFR, PDGFRA, MET, CDK4/6, MDM2/4, and deletion of CDKN2A/B, PTCH, and PTEN. Matched longitudinal PDOX recapitulate the limited genetic evolution of gliomas observed in patients following treatment. At the histological level, we observe increased vascularization in the rat host as compared to mice. PDOX-derived standardized glioma organoids are amenable to high-throughput drug screens that can be validated in mice. We show clinically relevant responses to temozolomide (TMZ) and to targeted treatments, such as EGFR and CDK4/6 inhibitors in (epi)genetically defined subgroups, according to MGMT promoter and EGFR/CDK status, respectively. Dianhydrogalactitol (VAL-083), a promising bifunctional alkylating agent in the current clinical trial, displayed high therapeutic efficacy, and was able to overcome TMZ resistance in glioblastoma. Our work underscores the clinical relevance of glioma organoids and PDOX models for translational research and personalized treatment studies and represents a unique publicly available resource for precision oncology.

Vascular dysfunction promotes regional hypoxia after bevacizumab therapy in recurrent glioblastoma patients

Background

Hypoxia is a driver of treatment resistance in glioblastoma. Antiangiogenic agents may transiently normalize blood vessels and decrease hypoxia before excessive pruning of vessels increases hypoxia. The time window of normalization is dose and time dependent. We sought to determine how VEGF blockade with bevacizumab modulates tumor vasculature and the impact that those vascular changes have on hypoxia in recurrent glioblastoma patients.

Methods

We measured tumor volume, vascular permeability (Ktrans), perfusion parameters (cerebral blood flow/volume, vessel caliber, and mean transit time), and regions of hypoxia in patients with recurrent glioblastoma before and after treatment with bevacizumab alone or with lomustine using [¹⁸F]FMISO PET-MRI. We also examined serial changes in plasma biomarkers of angiogenesis and inflammation.

Results

Eleven patients were studied. The magnitude of global tumor hypoxia was variable across these 11 patients prior to treatment and it did not significantly change after bevacizumab. The hypoxic regions had an inefficient vasculature characterized by elevated cerebral blood flow/volume and increased vessel caliber. In a subset of patients, there were tumor subregions with decreased mean transit times and a decrease in hypoxia, suggesting heterogeneous improvement in vascular efficiency. Bevacizumab significantly changed known pharmacodynamic biomarkers such as plasma VEGF and PlGF.

Conclusions

The vascular signature in hypoxic tumor regions indicates a disorganized vasculature which, in most tumors, does not significantly change after bevacizumab treatment. While some tumor regions showed improved vascular efficiency following treatment, bevacizumab did not globally alter hypoxia or normalize tumor vasculature in glioblastoma.

Resolution of Spatial and Temporal Heterogeneity in Bevacizumab-Treated Breast Tumors by Eigenspectra Multispectral Optoacoustic Tomography

Understanding temporal and spatial hemodynamic heterogeneity as a function of tumor growth or therapy affects the development of novel therapeutic strategies. In this study, we employed eigenspectra multispectral optoacoustic tomography (eMSOT) as a next-generation optoacoustic method to impart high accuracy in resolving tumor hemodynamics during bevacizumab therapy in two types of breast cancer xenografts (KPL-4 and MDA-MB-468). Patterns of tumor total hemoglobin concentration (THb) and oxygen saturation (sO₂) were imaged in control and bevacizumab-treated tumors over the course of 58 days (KPL-4) and 16 days (MDA-MB-468), and the evolution of functional vasculature "normalization" was resolved macroscopically. An initial sharp drop in tumor sO₂ and THb content shortly after the initiation of bevacizumab treatment was followed by a recovery in oxygenation levels. Rim-core subregion analysis revealed steep spatial oxygenation gradients in growing tumors that were reduced after bevacizumab treatment. Critically, eMSOT imaging findings were validated directly by histopathologic assessment of hypoxia (pimonidazole) and vascularity (CD31). These data demonstrate how eMSOT brings new abilities for accurate observation of entire tumor responses to challenges at spatial and temporal dimensions not available by other techniques today. SIGNIFICANCE: Accurate assessment of hypoxia and vascularization over space and time is critical for understanding tumor development and the role of spatial heterogeneity in tumor aggressiveness, metastasis, and response to treatment.

Evaluation of a tumor electric field treatment system in a rat model of glioma

Objective

Glioma is a devastating disease lacking effective treatment. Tumor electric field therapy is emerging as a novel non‐invasive therapy. The current study evaluates the efficacy and safety of a self‐designed tumor electric field therapy system (TEFTS ASCLU‐300) in a rat orthotopic transplantation model of glioma.

Methods

A model of intracranial orthotopic transplantation was established in rats using glioma C6 cells. For electric field therapy, glioma‐bearing rats were exposed to alternating electric fields generated by a self‐developed TEFTS starting on either 1st (Group 2) or 3rd (Group 3) day after transplantation, while other conditions were maintained the same as non‐treated rats (Group 1). Glioma size, body weight, and overall survival (OS) were compared between groups. Immunohistochemical staining was applied to access tumor cell death and microvessel density within the tumor. In addition, the systemic effects of TEFTS on blood cells, vital organs, and hepatorenal functions were evaluated.

Results

TEFTS treatment significantly elongated the OS of tumor‐bearing rats compared with non‐treated rats (non‐treated vs treated: 24.77 ± 7.08 days vs 40.31 ± 19.11 days, P = .0031). Continuous TEFTS treatment starting on 1st or 3rd day significantly reduced glioma size at 2 and 3 weeks after tumor cell inoculation (Week 2: Group 1:289.95 ± 101.69 mm³; Group 2:70.45 ± 17.79 mm³; Group 3:73.88 ± 33.21 mm³, P < .0001. Week 3: Group 1:544.096 ± 78.53 mm³; Group 2:187.58 ± 78.44 mm³; Group 3:167.14 ± 109.96 mm³, P = .0005). Continuous treatment for more than 4 weeks inhibited tumor growth. The TEFTS treatment promoted tumor cell death, as demonstrated by increased number of Caspase 3⁺ cells within the tumor (non‐treated vs treated: 38.06 ± 10.04 vs 68.57 ± 8.09 cells/field, P = .0007), but had minimal effect on microvessel density, as shown by CD31 expression (non‐treated vs treated: 1.63 ± 0.09 vs 1.57 ± 0.13% of positively stained areas, P > .05). No remarkable differences were observed in hepatorenal function, blood cell counts, or other vital organs between non‐treated and treated groups.

Conclusion

The TEFTS developed by our research team was proved to be effective and safe to inhibit tumor growth and improve general outcomes in a rat model of brain glioma.

Endothelial Progenitors in the Tumor Microenvironment

Tumor vascularization refers to the formation of new blood vessels within a tumor and is considered one of the hallmarks of cancer. Tumor vessels supply the tumor with oxygen and nutrients, required to sustain tumor growth and progression, and provide a gateway for tumor metastasis through the blood or lymphatic vasculature. Blood vessels display an angiocrine capacity of supporting the survival and proliferation of tumor cells through the production of growth factors and cytokines. Although tumor vasculature plays an essential role in sustaining tumor growth, it represents at the same time an essential way to deliver drugs and immune cells to the tumor. However, tumor vasculature exhibits many morphological and functional abnormalities, thus resulting in the formation of hypoxic areas within tumors, believed to represent a mechanism to maintain tumor cells in an invasive state.Tumors are vascularized through a variety of modalities, mainly represented by angiogenesis, where VEGF and other members of the VEGF family play a key role. This has represented the basis for the development of anti-VEGF blocking agents and their use in cancer therapy: however, these agents failed to induce significant therapeutic effects.Much less is known about the cellular origin of vessel network in tumors. Various cell types may contribute to tumor vasculature in different tumors or in the same tumor, such as mature endothelial cells, endothelial progenitor cells (EPCs), or the same tumor cells through a process of transdifferentiation. Early studies have suggested a role for bone marrow-derived EPCs; these cells do not are true EPCs but myeloid progenitors differentiating into monocytic cells, exerting a proangiogenic effect through a paracrine mechanism. More recent studies have shown the existence of tissue-resident endothelial vascular progenitors (EVPs) present at the level of vessel endothelium and their possible involvement as cells of origin of tumor vasculature.

The invasive proteome of glioblastoma revealed by laser-capture microdissection

Background

Glioblastomas are heterogeneous tumors composed of a necrotic and tumor core and an invasive periphery.

Methods

Here, we performed a proteomics analysis of laser-capture micro-dissected glioblastoma core and invasive areas of patient-derived xenografts.

Results

Bioinformatics analysis identified enriched proteins in central and invasive tumor areas. Novel markers of invasion were identified, the genes proteolipid protein 1 (PLP1) and Dynamin-1 (DNM1), which were subsequently validated in tumors and by functional assays.

Conclusions

In summary, our results identify new networks and molecules that may play an important role in glioblastoma development and may constitute potential novel therapeutic targets.

Spatially regularized estimation of the tissue homogeneity model parameters in DCE-MRI using proximal minimization

PURPOSE

The Tofts and the extended Tofts models are the pharmacokinetic models commonly used in dynamic contrast-enhanced MRI (DCE-MRI) perfusion analysis, although they do not provide two important biological markers, namely, the plasma flow and the permeability-surface area product. Estimates of such markers are possible using advanced pharmacokinetic models describing the vascular distribution phase, such as the tissue homogeneity model. However, the disadvantage of the advanced models lies in biased and uncertain estimates, especially when the estimates are computed voxelwise. The goal of this work is to improve the reliability of the estimates by including information from neighboring voxels.

THEORY AND METHODS

Information from the neighboring voxels is incorporated in the estimation process through spatial regularization in the form of total variation. The spatial regularization is applied on five maps of perfusion parameters estimated using the tissue homogeneity model. Since the total variation is not differentiable, two proximal techniques of convex optimization are used to solve the problem numerically.

RESULTS

The proposed algorithm helps to reduce noise in the estimated perfusion-parameter maps together with improving accuracy of the estimates. These conclusions are proved using a numerical phantom. In addition, experiments on real data show improved spatial consistency and readability of perfusion maps without considerable lowering of the quality of fit.

CONCLUSION

The reliability of the DCE-MRI perfusion analysis using the tissue homogeneity model can be improved by employing spatial regularization. The proposed utilization of modern optimization techniques implies only slightly higher computational costs compared to the standard approach without spatial regularization.

Blind deconvolution estimation of an arterial input function for small animal DCE-MRI

PURPOSE

One of the main obstacles for reliable quantitative dynamic contrast-enhanced (DCE) MRI is the need for accurate knowledge of the arterial input function (AIF). This is a special challenge for preclinical small animal applications where it is very difficult to measure the AIF without partial volume and flow artifacts. Furthermore, using advanced pharmacokinetic models (allowing estimation of blood flow and permeability-surface area product in addition to the classical perfusion parameters) poses stricter requirements on the accuracy and precision of AIF estimation. This paper addresses small animal DCE-MRI with advanced pharmacokinetic models and presents a method for estimation of the AIF based on blind deconvolution.

METHODS

A parametric AIF model designed for small animal physiology and use of advanced pharmacokinetic models is proposed. The parameters of the AIF are estimated using multichannel blind deconvolution.

RESULTS

Evaluation on simulated data show that for realistic signal to noise ratios blind deconvolution AIF estimation leads to comparable results as the use of the true AIF. Evaluation on real data based on DCE-MRI with two contrast agents of different molecular weights showed a consistence with the known effects of the molecular weight.

CONCLUSION

Multi-channel blind deconvolution using the proposed AIF model specific for small animal DCE-MRI provides reliable perfusion parameter estimates under realistic signal to noise conditions.

Angiogenesis and radiological tumor growth in patients with glioblastoma

BACKGROUND

The preoperative growth of human glioblastomas (GBMs) has been shown to vary among patients. In animal studies, angiogenesis has been linked to hypoxia and faster growth of GBM, however, its relation to the growth of human GBMs is sparsely studied. We have therefore aimed to look for associations between radiological speed of growth and microvessel density (MVD) counts of the endothelial markers vWF (Factor VIII related antigen) and CD105 (endoglin).

METHODS

Preoperative growth was estimated from segmented tumor volumes of two preoperative T1-weighted postcontrast magnetic resonance imaging scans taken ≥14 days apart in patients with newly diagnosed GBMs. A Gompertzian growth curve was computed from the volume data and separated the patients into two groups of either faster or slower tumor growth than expected. MVD counts of the immunohistochemical markers von Willebrand factor (vWF) (a pan-endothelial marker) and CD105 (a marker of proliferating endothelial cells) were assessed for associations with fast-growing tumors using Mann-Whitney U tests and a multivariable binary logistic regression analysis.

RESULTS

We found that only CD105-MVD was significantly associated with faster growth in a univariable analysis (p = 0.049). However, CD105-MVD was no longer significant when corrected for the presence of thromboses and high cellular density in a multivariable model, where the latter features were significant independent predictors of faster growth with respective odds ratios 4.2 (95% confidence interval, 1.2, 14.3), p = 0.021 and 2.6 (95% confidence interval, 1.0, 6.5), p = 0.048.

CONCLUSIONS

MVDs of neither endothelial marker were independently associated with faster growth, suggesting angiogenesis-independent processes contribute to faster glioblastoma growth.

Spatiotemporal heterogeneity of tumor vasculature during tumor growth and antiangiogenic treatment: MRI assessment using permeability and blood volume parameters

Tumor heterogeneity is an important concept when assessing intratumoral variety in vascular phenotypes and responses to antiangiogenic treatment. This study explored spatiotemporal heterogeneity of vascular alterations in C6 glioma mice during tumor growth and antiangiogenic treatment on serial MR examinations (days 0, 4, and 7 from initiation of vehicle or multireceptor tyrosine kinase inhibitor administration). Transvascular permeability (TP) was quantified on dynamic‐contrast‐enhanced MRI (DCE‐MRI) using extravascular extracellular agent (Gd‐DOTA); blood volume (BV) was estimated using intravascular T₂ agent (SPION). With regard to region‐dependent variability in vascular phenotypes, the control group demonstrated higher TP in the tumor center than in the periphery, and greater BV in the tumor periphery than in the center. This distribution pattern became more apparent with tumor growth. Antiangiogenic treatment effect was regionally heterogeneous: in the tumor center, treatment significantly suppressed the increase in TP and decrease in BV (ie, typical temporal change in the control group); in the tumor periphery, treatment‐induced vascular alterations were insignificant and BV remained high. On histopathological examination, the control group showed greater CD31, VEGFR2, Ki67, and NG2 expression in the tumor periphery than in the center. After treatment, CD31 and Ki67 expression was significantly suppressed only in the tumor center, whereas VEGFR2 and α‐caspase 3 expression was decreased and NG2 expression was increased in the entire tumor. These results demonstrate that MRI can reliably depict spatial heterogeneity in tumor vascular phenotypes and antiangiogenic treatment effects. Preserved angiogenic activity (high BV on MRI and high CD31) and proliferation (high Ki67) in the tumor periphery after treatment may provide insights into the mechanism of tumor resistance to antiangiogenic treatment.