Recurrence of glioblastoma is associated with elevated microvascular transit time heterogeneity and increased hypoxia

Quantitative BOLD (qBOLD) imaging of oxygen metabolism and blood oxygenation in the human body: A scoping review

PURPOSE

There are many approaches to the quantitative BOLD (qBOLD) technique described in the literature, differing in pulse sequences, MRI parameters and data processing. Thus, in this review, we summarized the acquisition methods, approaches used for oxygenation quantification and clinical populations investigated.

METHODS

Three databases were systematically searched (Medline, Embase, and Web of Science) for published research that used qBOLD methods for quantification of oxygen metabolism. Data extraction and synthesis were performed by one author and reviewed by a second author.

RESULTS

A total of 93 relevant papers were identified. Acquisition strategies were summarized, and oxygenation parameters were found to have been investigated in many pathologies such as steno-occlusive diseases, stroke, glioma, and multiple sclerosis disease.

CONCLUSION

A summary of qBOLD approaches for oxygenation measurements and applications could help researchers to identify good practice and provide objective information to inform the development of future consensus recommendations.

Hemodynamic response function description in patients with glioma

INTRODUCTION

Functional magnetic resonance imaging is a powerful tool that has provided many insights into cognitive sciences. Yet, as its analysis is mostly based on the knowledge of an a priori canonical hemodynamic response function (HRF), its reliability in patients' applications has been questioned. There have been reports of neurovascular uncoupling in patients with glioma, but no specific description of the Hemodynamic Response Function (HRF) in glioma has been reported so far. The aim of this work is to describe the HRF in patients with glioma.

METHODS

Forty patients were included. MR images were acquired on a 1.5T scanner. Activated clusters were identified using a fuzzy general linear model; HRFs were adjusted with a double-gamma function. Analyses were undertaken considering the tumor grade, age, sex, tumor location, and activated location.

RESULTS

Differences are found in the occipital, limbic, insular, and sub-lobar areas, but not in the frontal, temporal, and parietal lobes. The presence of a glioma slows the time-to-peak and onset times by 5.2 and 3.8 % respectively; high-grade gliomas present 8.1 % smaller HRF widths than low-grade gliomas.

DISCUSSION AND CONCLUSION

There is significant HRF variation due to the presence of glioma, but the magnitudes of the observed differences are small. Most processing pipelines should be robust enough for this magnitude of variation and little if any impact should be visible on functional maps. The differences that have been observed in the literature between functional mapping obtained with magnetic resonance vs. that obtained with direct electrostimulation during awake surgery are more probably due to the intrinsic difference in the mapping process: fMRI mapping detects all recruited areas while intra-surgical mapping indicates only the areas indispensable for the realization of a certain task. Surgical mapping might not be the gold standard to use when trying to validate the fMRI mapping process.

Prediction of pseudoprogression in post-treatment glioblastoma using dynamic susceptibility contrast-derived oxygenation and microvascular transit time heterogeneity measures

OBJECTIVES

To evaluate the added value of MR dynamic susceptibility contrast (DSC)-perfusion-weighted imaging (PWI)-derived tumour microvascular and oxygenation information with cerebral blood volume (CBV) to distinguish pseudoprogression from true progression (TP) in post-treatment glioblastoma.

METHODS

This retrospective single-institution study included patients with isocitrate dehydrogenase (IDH) wild-type glioblastoma and a newly developed or enlarging measurable contrast-enhancing mass within 12 weeks after concurrent chemoradiotherapy. CBV, capillary transit time heterogeneity (CTH), oxygen extraction fraction (OEF), and cerebral metabolic rate of oxygen (CMRO2) were obtained from DSC-PWI. Predictors were selected using univariable logistic regression, and performance was measured with adjusted diagnostic odds with tumour volume and area under the curve (AUC) of receiver operating characteristics analysis.

RESULTS

A total of 103 patients were included (mean age, 59.6 years; 59 women), with 67 cases of TP and 36 cases of pseudoprogression. Pseudoprogression exhibited higher CTH (4.0 vs. 3.4, p = .019) and higher OEF (12.7 vs. 10.7, p = .014) than TP, but a similar CBV (1.48 vs. 1.53, p = .13) and CMRO2 (7.7 vs. 7.3s, p = .598). Independent of tumour volume, both high CTH (adjusted odds ratio [OR] 1.52; 95% confidence interval [CI]: 1.11-2.09, p = .009) and high OEF (adjusted OR 1.17; 95% CI:1.03-1.33, p = .016) were predictors of pseudoprogression. The combination of CTH, OEF, and CBV yielded higher diagnostic performance (AUC 0.71) than CBV alone (AUC 0.65).

CONCLUSION

High intratumoural capillary transit heterogeneity and high oxygen extraction fraction derived from DSC-PWI have enhanced the diagnostic value of CBV in pseudoprogression of post-treatment IDH-wild type glioblastoma.

CLINICAL RELEVANCE STATEMENT

In the early post-treatment stage of glioblastoma, pseudoprogression exhibited both high oxygen extraction fraction and high capillary transit heterogeneity and these dynamic susceptibility contrast-perfusion weighted imaging derived parameters have added value in cerebral blood volume-based noninvasive differentiation of pseudoprogression from true progression.

KEY POINTS

• Capillary transit time heterogeneity and oxygen extraction fraction can be measured noninvasively through processing of dynamic susceptibility contrast imaging. • Pseudoprogression exhibited higher capillary transit time heterogeneity and higher oxygen extraction fraction than true progression. • A combination of cerebral blood volume, capillary transit time heterogeneity, and oxygen extraction fraction yielded the highest diagnostic performance (area under the curve 0.71).

MRI-Based Assessment of Brain Tumor Hypoxia: Correlation with Histology

Cerebral hypoxia significantly impacts the progression of brain tumors and their resistance to radiotherapy. This study employed streamlined quantitative blood-oxygen-level-dependent (sqBOLD) MRI to assess the oxygen extraction fraction (OEF)-a measure of how much oxygen is being extracted from vessels, with higher OEF values indicating hypoxia. Simultaneously, we utilized vessel size imaging (VSI) to evaluate microvascular dimensions and blood volume. A cohort of ten patients, divided between those with glioma and those with brain metastases, underwent a 3 Tesla MRI scan. We generated OEF, cerebral blood volume (CBV), and vessel size maps, which guided 3-4 targeted biopsies per patient. Subsequent histological analyses of these biopsies used hypoxia-inducible factor 1-alpha (HIF-1α) for hypoxia and CD31 for microvasculature assessment, followed by a correlation analysis between MRI and histological data. The results showed that while the sqBOLD model was generally applicable to brain tumors, it demonstrated discrepancies in some metastatic tumors, highlighting the need for model adjustments in these cases. The OEF, CBV, and vessel size maps provided insights into the tumor's hypoxic condition, showing intertumoral and intratumoral heterogeneity. A significant relationship between MRI-derived measurements and histological data was only evident in the vessel size measurements (r = 0.68, p < 0.001).

Streamlined quantitative BOLD for detecting visual stimulus-induced changes in oxygen extraction fraction in healthy participants: toward clinical application in human glioma

OBJECTIVE

Monitoring brain oxygenation is critical in brain tumors, as low oxygenation influences tumor growth, pathological angiogenesis, and treatment resistance. This study examined the ability of the streamlined quantitative (sq)BOLD MRI technique to detect oxygenation changes in healthy individuals, as well as its potential application in a clinical setting.

METHODS

We used the asymmetric spin echo (ASE) technique with FLAIR preparation, along with model-based Bayesian inference to quantify the reversible transverse relaxation rate (R2') and oxygen extraction fraction (OEF) across the brain at baseline and during visual stimulation in eight healthy participants at 3T; and two patients with glioma at rest only.

RESULTS

Comparing sqBOLD-derived parameters between baseline and visual stimulation revealed a decrease in OEF from 0.56 ± 0.09 at baseline to 0.54 ± 0.07 at the activated state (p = 0.04, paired t test) within a functional localizer-defined volume of interest, and a decline in R2' from 6.5 ± 1.3s-1 at baseline to 6.2 ± 1.4s-1 at the activated state (p = 0.006, paired t test) in the visual cortex.

CONCLUSION

The sqBOLD technique is sensitive enough to detect and quantify changes in oxygenation in the healthy brain and shows potential for integration into clinical settings to provide valuable information on oxygenation in glioma.

Glioblastomas with and without peritumoral fluid-attenuated inversion recovery (FLAIR) hyperintensity present morphological and microstructural differences on conventional MR images

OBJECTIVES

Glioblastoma (GB) without peritumoral fluid-attenuated inversion recovery (FLAIR) hyperintensity is atypical and its characteristics are barely known. The aim of this study was to explore the differences in pathological and MRI-based intrinsic features (including morphologic and first-order features) between GBs with peritumoral FLAIR hyperintensity (PFH-bearing GBs) and GBs without peritumoral FLAIR hyperintensity (PFH-free GBs).

METHODS

In total, 155 patients with pathologically diagnosed GBs were retrospectively collected, which included 110 PFH-bearing GBs and 45 PFH-free GBs. The pathological and imaging data were collected. The Visually AcceSAble Rembrandt Images (VASARI) features were carefully evaluated. The first-order radiomics features from the tumor region were extracted from FLAIR, apparent diffusion coefficient (ADC), and T1CE (T1-contrast enhanced) images. All parameters were compared between the two groups of GBs.

RESULTS

The pathological data showed more alpha thalassemia/mental retardation syndrome X-linked (ATRX)-loss in PFH-free GBs compared to PFH-bearing ones (p < 0.001). Based on VASARI evaluation, PFH-free GBs had larger intra-tumoral enhancing proportion and smaller necrotic proportion (both, p < 0.001), more common non-enhancing tumor (p < 0.001), mild/minimal enhancement (p = 0.003), expansive T1/FLAIR ratio (p < 0.001) and solid enhancement (p = 0.009), and less pial invasion (p = 0.010). Moreover, multiple ADC- and T1CE-based first-order radiomics features demonstrated differences, especially the lower intensity heterogeneity in PFH-free GBs (for all, adjusted p < 0.05).

CONCLUSIONS

Compared to PFH-bearing GBs, PFH-free ones demonstrated less immature neovascularization and lower intra-tumoral heterogeneity, which would be helpful in clinical treatment stratification.

CLINICAL RELEVANCE STATEMENT

Glioblastomas without peritumoral FLAIR hyperintensity show less immature neovascularization and lower heterogeneity leading to potential higher treatment benefits due to less drug resistance and treatment failure.

KEY POINTS

• The study explored the differences between glioblastomas with and without peritumoral FLAIR hyperintensity. • Glioblastomas without peritumoral FLAIR hyperintensity showed less necrosis and contrast enhancement and lower intensity heterogeneity. • Glioblastomas without peritumoral FLAIR hyperintensity had less immature neovascularization and lower tumor heterogeneity.

Can Perfusion-Based Brain Tissue Oxygenation MRI Support the Understanding of Cerebral Abscesses In Vivo?

PURPOSE

The clinical condition of a brain abscess is a potentially life-threatening disease. The combination of MRI-based imaging, surgical therapy and microbiological analysis is critical for the treatment and convalescence of the individual patient. The aim of this study was to evaluate brain tissue oxygenation measured with dynamic susceptibility contrast perfusion weighted imaging (DSC-PWI) in patients with brain abscess and its potential benefit for a better understanding of the environment in and around brain abscesses.

METHODS

Using a local database, 34 patients (with 45 abscesses) with brain abscesses treated between January 2013 and March 2021 were retrospectively included in this study. DSC-PWI imaging and microbiological work-up were key inclusion criteria. These data were analysed regarding a correlation between DSC-PWI and microbiological result by quantifying brain tissue oxygenation in the abscess itself, the abscess capsula and the surrounding oedema and by using six different parameters (CBF, CBV, CMRO2, COV, CTH and OEF).

RESULTS

Relative cerebral blood flow (0.335 [0.18-0.613] vs. 0.81 [0.49-1.08], p = 0.015), relative cerebral blood volume (0.44 [0.203-0.72] vs. 0.87 [0.67-1.2], p = 0.018) and regional cerebral metabolic rate for oxygen (0.37 [0.208-0.695] vs. 0.82 [0.55-1.19], p = 0.022) were significantly lower in the oedema around abscesses without microbiological evidence of a specific bacteria in comparison with microbiological positive lesions.

CONCLUSIONS

The results of this study indicate a relationship between brain tissue oxygenation status in DSC-PWI and microbiological/inflammatory status. These results may help to better understand the in vivo environment of brain abscesses and support future therapeutic decisions.

Longitudinal study on MRI and neuropathological findings: Neither DSC-perfusion derived rCBVmax nor vessel densities correlate between newly diagnosed and progressive glioblastoma

INTRODUCTION

When evaluating MRIs for glioblastoma progression, previous scans are usually included into the review. Nowadays dynamic susceptibility contrast (DSC)-perfusion is an essential component in MR-diagnostics of gliomas, since the extent of hyperperfusion upon first diagnosis correlates with gene expression and survival. We aimed to investigate if this initial perfusion signature also characterizes the glioblastoma at time of progression. If so, DSC-perfusion data from the initial diagnosis could be of diagnostic benefit in follow-up assessments.

METHODS

We retrospectively identified 65 patients with isocitrate dehydrogenase wildtype glioblastoma who had received technically identical DSC-perfusion measurements at initial diagnosis and at time of first progression. We determined maximum relative cerebral blood volume values (rCBVmax) by standardized re-evaluation of the data including leakage correction. In addition, the corresponding tissue samples from 24 patients were examined histologically for the maximum vessel density within the tumor. Differences (paired t-test/ Wilcoxon matched pairs test) and correlations (Spearman) between the measurements at both timepoints were calculated.

RESULTS

The rCBVmax was consistently lower at time of progression compared to rCBVmax at time of first diagnosis (p < .001). There was no correlation between the rCBVmax values at both timepoints (r = .12). These findings were reflected in the histological examination, with a lower vessel density in progressive glioblastoma (p = .01) and no correlation between the two timepoints (r = -.07).

CONCLUSION

Our results suggest that the extent of hyperperfusion in glioblastoma at first diagnosis is not a sustaining tumor characteristic. Hence, the rCBVmax at initial diagnosis should be disregarded when reviewing MRIs for glioblastoma progression.

Differentiation of Glioblastoma and Brain Metastases by MRI-Based Oxygen Metabolomic Radiomics and Deep Learning

Glioblastoma (GB) and brain metastasis (BM) are the most frequent types of brain tumors in adults. Their therapeutic management is quite different and a quick and reliable initial characterization has a significant impact on clinical outcomes. However, the differentiation of GB and BM remains a major challenge in today's clinical neurooncology due to their very similar appearance in conventional magnetic resonance imaging (MRI). Novel metabolic neuroimaging has proven useful for improving diagnostic performance but requires artificial intelligence for implementation in clinical routines. Here; we investigated whether the combination of radiomic features from MR-based oxygen metabolism ("oxygen metabolic radiomics") and deep convolutional neural networks (CNNs) can support reliably pre-therapeutic differentiation of GB and BM in a clinical setting. A self-developed one-dimensional CNN combined with radiomic features from the cerebral metabolic rate of oxygen (CMRO₂) was clearly superior to human reading in all parameters for classification performance. The radiomic features for tissue oxygen saturation (mitoPO₂; i.e., tissue hypoxia) also showed better diagnostic performance compared to the radiologists. Interestingly, both the mean and median values for quantitative CMRO₂ and mitoPO₂ values did not differ significantly between GB and BM. This demonstrates that the combination of radiomic features and DL algorithms is more efficient for class differentiation than the comparison of mean or median values. Oxygen metabolic radiomics and deep neural networks provide insights into brain tumor phenotype that may have important diagnostic implications and helpful in clinical routine diagnosis.

Impact of molecular and clinical variables on survival outcome with immunotherapy for glioblastoma patients: A systematic review and meta-analysis

BACKGROUND

Given that only a subset of patients with glioblastoma multiforme (GBM) responds to immuno-oncology, this study aimed to assess the impact of multiple factors on GBM immunotherapy prognosis and investigate the potential predictors.

METHODS

A quantitative meta-analysis was conducted using the random-effects model. Several potential factors were also reviewed qualitatively.

RESULTS

A total of 39 clinical trials were included after screening 1317 papers. Patients with O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation [hazard ratio (HR) for overall survival (OS) = 2.30, p < 0.0001; HR for progression-free survival (PFS) = 2.10, p < 0.0001], gross total resection (HR for OS = 0.70, p = 0.02; HR for PFS = 0.56, p = 0.004), and no baseline steroid use (HR for OS = 0.52, p = 0.0002; HR for PFS = 0.61, p = 0.02) had a relatively significant favorable OS and PFS following immunotherapy. Patients with a Karnofsky Performance Status score < 80 (HR = 1.73, p = 0.0007) and undergoing two prior relapses (HR = 2.08, p = 0.003) were associated with worse OS. Age, gender, tumor programmed death-ligand 1 expression, and history of chemotherapy were not associated with survival outcomes. Notably, immunotherapy significantly improved the OS among patients undergoing two prior recurrences (HR = 0.40, p = 0.008) but not among patients in any other subgroups, as opposed to non-immunotherapy.

CONCLUSION

Several factors were associated with prognostic outcomes of GBM patients receiving immunotherapy; multiple recurrences might be a candidate predictor. More marker-driven prospective studies are warranted.

Decoding the Heterogeneity of Malignant Gliomas by PET and MRI for Spatial Habitat Analysis of Hypoxia, Perfusion, and Diffusion Imaging: A Preliminary Study

Background

Tumor heterogeneity poses major clinical challenges in high-grade gliomas (HGGs). Quantitative radiomic analysis with spatial tumor habitat clustering represents an innovative, non-invasive approach to represent and quantify tumor microenvironment heterogeneity. To date, habitat imaging has been applied mainly on conventional magnetic resonance imaging (MRI), although virtually extendible to any imaging modality, including advanced MRI techniques such as perfusion and diffusion MRI as well as positron emission tomography (PET) imaging.

Objectives

This study aims to evaluate an innovative PET and MRI approach for assessing hypoxia, perfusion, and tissue diffusion in HGGs and derive a combined map for clustering of intra-tumor heterogeneity.

Materials and Methods

Seventeen patients harboring HGGs underwent a pre-operative acquisition of MR perfusion (PWI), Diffusion (dMRI) and 18F-labeled fluoroazomycinarabinoside (18F-FAZA) PET imaging to evaluate tumor vascularization, cellularity, and hypoxia, respectively. Tumor volumes were segmented on fluid-attenuated inversion recovery (FLAIR) and T1 post-contrast images, and voxel-wise clustering of each quantitative imaging map identified eight combined PET and physiologic MRI habitats. Habitats’ spatial distribution, quantitative features and histopathological characteristics were analyzed.

Results

A highly reproducible distribution pattern of the clusters was observed among different cases, particularly with respect to morphological landmarks as the necrotic core, contrast-enhancing vital tumor, and peritumoral infiltration and edema, providing valuable supplementary information to conventional imaging. A preliminary analysis, performed on stereotactic bioptic samples where exact intracranial coordinates were available, identified a reliable correlation between the expected microenvironment of the different spatial habitats and the actual histopathological features. A trend toward a higher representation of the most aggressive clusters in WHO (World Health Organization) grade IV compared to WHO III was observed.

Conclusion

Preliminary findings demonstrated high reproducibility of the PET and MRI hypoxia, perfusion, and tissue diffusion spatial habitat maps and correlation with disease-specific histopathological features.

Cerebral oxygen extraction fraction MRI: Techniques and applications

The human brain constitutes 2% of the body's total mass but uses 20% of the oxygen. The rate of the brain's oxygen utilization can be derived from a knowledge of cerebral blood flow and the oxygen extraction fraction (OEF). Therefore, OEF is a key physiological parameter of the brain's function and metabolism. OEF has been suggested to be a useful biomarker in a number of brain diseases. With recent advances in MRI techniques, several MRI-based methods have been developed to measure OEF in the human brain. These MRI OEF techniques are based on the T₂ of blood, the blood signal phase, the magnetic susceptibility of blood-containing voxels, the effect of deoxyhemoglobin on signal behavior in extravascular tissue, and the calibration of the BOLD signal using gas inhalation. Compared to ¹⁵ O PET, which is considered the "gold standard" for OEF measurement, MRI-based techniques are non-invasive, radiation-free, and are more widely available. This article provides a review of these emerging MRI-based OEF techniques. We first briefly introduce the role of OEF in brain oxygen homeostasis. We then review the methodological aspects of different categories of MRI OEF techniques, including their signal mechanisms, acquisition methods, and data analyses. The strengths and limitations of the techniques are discussed. Finally, we review key applications of these techniques in physiological and pathological conditions.

Association between tissue hypoxia, perfusion restrictions, and microvascular architecture alterations with lesion-induced impairment of neurovascular coupling

Functional magnetic resonance imaging (fMRI) has been mainly utilized for the preoperative localization of eloquent cortical areas. However, lesion-induced impairment of neurovascular coupling (NVC) in the lesion border zone may lead to false-negative fMRI results. The purpose of this study was to determine physiological factors impacting the NVC. Twenty patients suffering from brain lesions were preoperatively examined using multimodal neuroimaging including fMRI, magnetoencephalography (MEG) during language or sensorimotor tasks (depending on lesion location), and a novel physiologic MRI approach for the combined quantification of oxygen metabolism, perfusion state, and microvascular architecture. Congruence of brain activity patterns between fMRI and MEG were found in 13 patients. In contrast, we observed missing fMRI activity in perilesional cortex that demonstrated MEG activity in seven patients, which was interpreted as lesion-induced impairment of NVC. In these brain regions with impaired NVC, physiologic MRI revealed significant brain tissue hypoxia, as well as significantly decreased macro- and microvascular perfusion and microvascular architecture. We demonstrated that perilesional hypoxia with reduced vascular perfusion and architecture is associated with lesion-induced impairment of NVC. Our physiologic MRI approach is a clinically applicable method for preoperative risk assessment for the presence of false-negative fMRI results and may prevent severe postoperative functional deficits.

Imaging effective oxygen diffusivity in the human brain with multiparametric magnetic resonance imaging

Cerebrovascular diseases can impair blood circulation and oxygen extraction from the blood. The effective oxygen diffusivity (EOD) of the capillary bed is a potential biomarker of microvascular function that has gained increasing interest, both for clinical diagnosis and for elucidating oxygen transport mechanisms. Models of capillary oxygen transport link EOD to measurable oxygen extraction fraction (OEF) and cerebral blood flow (CBF). In this work, we confirm that two well established mathematical models of oxygen transport yield nearly equivalent EOD maps. Furthermore, we propose an easy-to-implement and clinically applicable multiparametric magnetic resonance imaging (MRI) protocol for quantitative EOD mapping. Our approach is based on imaging OEF and CBF with multiparametric quantitative blood oxygenation level dependent (mq-BOLD) MRI and pseudo-continuous arterial spin labeling (pCASL), respectively. We evaluated the imaging protocol by comparing MRI-EOD maps of 12 young healthy volunteers to PET data from a published study in different individuals. Our results show comparably good correlation between MRI- and PET-derived cortical EOD, OEF and CBF. Importantly, absolute values of MRI and PET showed high accordance for all three parameters. In conclusion, our data indicates feasibility of the proposed MRI protocol for EOD mapping, rendering the method promising for future clinical evaluation of patients with cerebrovascular diseases.

Exploring the efficacy of tumor electric field therapy against glioblastoma: An in vivo and in vitro study

AIMS

Tumor electric fields therapy (TTFields) is emerging as a novel anti-cancer physiotherapy. Despite recent breakthroughs of TTFields in glioma treatment, the average survival time for glioblastoma patients with TTFields is <2 years, even when used in conjugation with traditional anti-cancer therapies. To optimize TTFields-afforded efficacy against glioblastoma, we investigated the cancer cell-killing effects of various TTFields paradigms using in vitro and in vivo models of glioblastoma.

METHODS

For in vitro studies, the U251 glioma cell line or primary cell cultures prepared from 20 glioblastoma patients were treated with the tumor electric field treatment (TEFT) system. Cell number, volume, and proliferation were measured after TEFT at different frequencies (100, 150, 180, 200, or 220 kHz), durations (24, 48, or 72 h), field strengths (1.0, 1.5, or 2.2V/cm), and output modes (fixed or random sequence output). A transwell system was used to evaluate the influence of TEFT on the invasiveness of primary glioblastoma cells. For in vivo studies, the therapeutic effect and safety profiles of random sequence electric field therapy in glioblastoma-transplanted rats were assessed by calculating tumor size and survival time and evaluating peripheral immunobiological and blood parameters, respectively.

RESULTS

In the in vitro settings, TEFT was robustly effective in suppressing cell proliferation of both the U251 glioma cell line and primary glioblastoma cell cultures. The anti-proliferation effects of TEFT were frequency- and "dose" (field strength and duration)-dependent, and contingent on the field sequence output mode, with the random sequence mode (TEFT-R) being more effective than the fixed sequence mode (TEFT-F). Genetic tests were performed in 11 of 20 primary glioblastoma cultures, and 6 different genetic traits were identified them. However, TEFT exhibited comparable anti-proliferation effects in all primary cultures regardless of their genetic traits. TEFT also inhibited the invasiveness of primary glioblastoma cells in transwell experiments. In the in vivo rat model of glioblastoma brain transplantation, treatment with TEFT-F or TEFT-R at frequency of 200 kHz and field strength of 2.2V/cm for 14 days significantly reduced tumor volume by 42.63% (TEFT-F vs. control, p = 0.0002) and 63.60% (TEFT-R vs. control, p < 0.0001), and prolonged animal survival time by 30.15% (TEFT-F vs. control, p = 0.0415) and 69.85% (TEFT-R vs. control, p = 0.0064), respectively. The tumor-bearing rats appeared to be well tolerable to TEFT therapies, showing only moderate increases in blood levels of creatine and red blood cells. Adverse skin reactions were common for TEFT-treated rats; however, skin reactions were curable by local treatment.

CONCLUSION

Tumor electric field treatment at optimal frequency, strength, and output mode markedly inhibits the cell viability, proliferation, and invasiveness of primary glioblastoma cells in vitro independent of different genetic traits of the cells. Moreover, a random sequence electric field output confers considerable anti-cancer effects against glioblastoma in vivo. Thus, TTFields are a promising physiotherapy for glioblastoma and warrants further investigation.

Metabolic Tumor Microenvironment Characterization of Contrast Enhancing Brain Tumors Using Physiologic MRI

The tumor microenvironment is a critical regulator of cancer development and progression as well as treatment response and resistance in brain neoplasms. The available techniques for investigation, however, are not well suited for noninvasive in vivo characterization in humans. A total of 120 patients (59 females; 61 males) with newly diagnosed contrast-enhancing brain tumors (64 glioblastoma, 20 brain metastases, 15 primary central nervous system (CNS) lymphomas (PCNSLs), and 21 meningiomas) were examined with a previously established physiological MRI protocol including quantitative blood-oxygen-level-dependent imaging and vascular architecture mapping. Six MRI biomarker maps for oxygen metabolism and neovascularization were fused for classification of five different tumor microenvironments: glycolysis, oxidative phosphorylation (OxPhos), hypoxia with/without neovascularization, and necrosis. Glioblastoma showed the highest metabolic heterogeneity followed by brain metastasis with a glycolysis-to-OxPhos ratio of approximately 2:1 in both tumor entities. In addition, glioblastoma revealed a significant higher percentage of hypoxia (24%) compared to all three other brain tumor entities: brain metastasis (7%; p < 0.001), PCNSL (8%; p = 0.001), and meningioma (8%; p = 0.003). A more aggressive biological brain tumor behavior was associated with a higher percentage of hypoxia and necrosis and a lower percentage of remaining vital tumor tissue and aerobic glycolysis. The proportion of oxidative phosphorylation, however, was rather similar (17–26%) for all four brain tumor entities. Tumor microenvironment (TME) mapping provides insights into neurobiological differences of contrast-enhancing brain tumors and deserves further clinical cancer research attention. Although there is a long roadmap ahead, TME mapping may become useful in order to develop new diagnostic and therapeutic approaches.

Physiological MRI of microvascular architecture, neovascularization activity, and oxygen metabolism facilitate early recurrence detection in patients with IDH-mutant WHO grade 3 glioma

PURPOSE

This study aimed to determine the diagnostic performance of physiological MRI biomarkers including microvascular perfusion and architecture, neovascularization activity, tissue oxygen metabolism, and tension for recurrence detection of IDH-mutant WHO grade 3 glioma.

METHODS

Sixty patients with IDH-mutant WHO grade 3 glioma who received overall 288 follow-up MRI examinations at 3 Tesla after standard treatment were retrospectively evaluated. A conventional MRI protocol was extended with a physiological MRI approach including vascular architecture mapping and quantitative blood-oxygen-level-dependent imaging which required 7 min extra data acquisition time. Custom-made MATLAB software was used for the calculation of MRI biomarker maps of microvascular perfusion and architecture, neovascularization activity, tissue oxygen metabolism, and tension. Statistical procedures included receiver operating characteristic analysis.

RESULTS

Overall, 34 patients showed recurrence of the WHO grade 3 glioma; of these, in 15 patients, recurrence was detected one follow-up examination (averaged 160 days) earlier by physiological MRI data than by conventional MRI. During this time period, the tumor volume increased significantly (P = 0.001) on average 7.4-fold from 1.5 to 11.1 cm³. Quantitative analysis of MRI biomarkers demonstrated microvascular but no macrovascular hyperperfusion in early recurrence. Neovascularization activity (AUC = 0.833), microvascular perfusion (0.682), and oxygen metabolism (0.661) showed higher diagnostic performance for early recurrence detection of WHO grade 3 glioma compared to conventional MRI including cerebral blood volume (0.649).

CONCLUSION

This study demonstrated that the targeted assessment of microvascular features and tissue oxygen tension as an early sign of neovascularization activity provided valuable information for recurrence diagnostic of WHO grade 3 glioma.

Cerebral and tumoral blood flow in adult gliomas: a systematic review of results from magnetic resonance imaging

Objective:

Blood flow is the rate of blood movement and relevant to numerous processes, though understudied in gliomas. The aim of this review was to pool blood flow metrics obtained from MRI modalities in adult supratentorial gliomas.

Methods:

MEDLINE, EMBASE and the Cochrane database were queried 01/01/2000–31/12/2019. Studies measuring blood flow in adult Grade II–IV supratentorial gliomas using dynamic susceptibility contrast (DSC) MRI, dynamic contrast enhanced MRI (DCE-MRI) or arterial spin labelling (ASL) were included. Absolute and relative cerebral blood flow (CBF), peritumoral blood flow and tumoral blood flow (TBF) were reported.

Results:

34 studies were included with 1415 patients and 1460 scans. The mean age was 52.4 ± 7.3 years. Most patients had glioblastoma (n = 880, 64.6%). The most common imaging modality was ASL (n = 765, 52.4%) followed by DSC (n = 538, 36.8%). Most studies were performed pre-operatively (n = 1268, 86.8%). With increasing glioma grade (II vs IV), TBF increased (70.8 vs 145.5 ml/100 g/min, p < 0.001) and CBF decreased (85.3 vs 49.6 ml/100 g/min, p < 0.001). In Grade IV gliomas, following treatment, CBF increased in ipsilateral (24.9 ± 1.2 vs 26.1 ± 0.0 ml/100 g/min, p < 0.001) and contralateral white matter (25.6 ± 0.2 vs 26.0± 0.0 ml/100 g/min, p < 0.001).

Conclusion:

Our findings demonstrate that increased mass effect from high-grade gliomas impairs blood flow within the surrounding brain that can improve with surgery.

Advances in knowledge:

This systematic review demonstrates how mass effect from brain tumours impairs blood flow in the surrounding brain parenchyma that can improve with treatment.

Physiological MRI Biomarkers in the Differentiation Between Glioblastomas and Solitary Brain Metastases

PURPOSE

Glioblastomas (GB) and solitary brain metastases (BM) are the most common brain tumors in adults. GB and BM may appear similar in conventional magnetic resonance imaging (cMRI). Their management strategies, however, are quite different with significant consequences on clinical outcome. The aim of this study was to evaluate the usefulness of a previously presented physiological MRI approach scoping to obtain quantitative information about microvascular architecture and perfusion, neovascularization activity, and oxygen metabolism to differentiate GB from BM.

PROCEDURES

Thirty-three consecutive patients with newly diagnosed, untreated, and histopathologically confirmed GB or BM were preoperatively examined with our physiological MRI approach as part of the cMRI protocol.

RESULTS

Physiological MRI biomarker maps revealed several significant differences in the pathophysiology of GB and BM: Central necrosis was more hypoxic in GB than in BM (30 %; P = 0.036), which was associated with higher neovascularization activity (65 %; P = 0.043) and metabolic rate of oxygen (48 %; P = 0.004) in the adjacent contrast-enhancing viable tumor parts of GB. In peritumoral edema, GB infiltration caused neovascularization activity (93 %; P = 0.018) and higher microvascular perfusion (30 %; P = 0.022) associated with higher tissue oxygen tension (33 %; P = 0.020) and lower oxygen extraction from vasculature (32 %; P = 0.040).

CONCLUSION

Our physiological MRI approach, which requires only 7 min of extra data acquisition time, might be helpful to noninvasively distinguish GB and BM based on pathophysiological differences. However, further studies including more patients are required.

Hypoxia and Microvascular Alterations Are Early Predictors of IDH-Mutated Anaplastic Glioma Recurrence

Simple Summary

Anaplastic gliomas (AGs) are considered the most common and aggressive primary brain tumors of young adults with inevitable recurrence and treatment failure. The aim of this study was to investigate whether the imaging biomarkers of hypoxia, microvascular architecture and neovascularization activity can be of assistance to detect pathophysiological changes in the early developmental stages of isocitrate-dehydrogenase (IDH) mutated AG recurrence. We evaluated 142 physiological magnetic resonance imaging follow-up examinations as a part of the conventional magnetic resonance imaging (MRI) protocol in 60 AG patients after standard therapy. Physiological MRI biomarkers showed intensifying local tissue hypoxia 250 days prior to radiological recurrence with following upregulation of neovascularization activity 50 to 70 days later. Integration of physiological MRI in the monitoring of AG patients may be of clinical significance to make personalized decision of early tumor recurrence without an additional delay for multimodal therapy.

Abstract

Anaplastic gliomas (AG) represents aggressive brain tumors that often affect young adults. Although isocitrate-dehydrogenase (IDH) gene mutation has been identified as a more favorable prognostic factor, most IDH-mutated AG patients are confronted with tumor recurrence. Hence, increased knowledge about pathophysiological precursors of AG recurrence is urgently needed in order to develop precise diagnostic monitoring and tailored therapeutic approaches. In this study, 142 physiological magnetic resonance imaging (phyMRI) follow-up examinations in 60 AG patients after standard therapy were evaluated and magnetic resonance imaging (MRI) biomarker maps for microvascular architecture and perfusion, neovascularization activity, oxygen metabolism, and hypoxia calculated. From these 60 patients, 34 patients developed recurrence of the AG, and 26 patients showed no signs for AG recurrence during the study period. The time courses of MRI biomarker changes were analyzed regarding early pathophysiological alterations over a one-year period before radiological AG recurrence or a one-year period of stable disease for patients without recurrence, respectively. We detected intensifying local tissue hypoxia 250 days prior to radiological recurrence which initiated upregulation of neovascularization activity 50 to 70 days later. These changes were associated with a switch from an avascular infiltrative to a vascularized proliferative phenotype of the tumor cells another 30 days later. The dynamic changes of blood perfusion, microvessel density, neovascularization activity, and oxygen metabolism showed a close physiological interplay in the one-year period prior to radiological recurrence of IDH-mutated AG. These findings may path the wave for implementing both new MR-based imaging modalities for routine follow-up monitoring of AG patients after standard therapy and furthermore may support the development of novel, tailored therapy options in recurrent AG.

Expression Profiling of Primary and Recurrent Glioblastomas Reveals a Reduced Level of Pentraxin 3 in Recurrent Glioblastomas

Glioblastomas (GBM) are highly infiltrative tumors and despite intensive treatment tumor recurrence is inevitable. The immune microenvironment in recurrent GBM is poorly characterized, but it is potentially influenced by therapeutic interventions with surgery, radiotherapy, and chemotherapy. The aim of this study was to obtain a deeper insight in the immune microenvironment in primary and recurrent GBM. Primary and recurrent glioblastoma samples from 18 patients were identified and expression profiling of 770 myeloid innate immune-related markers was performed. Leukemia inhibitory factor and pentraxin 3 were expressed at lower levels in recurrent tumors. Using in silico data and immunohistochemical staining, this was validated for pentraxin 3. Both high leukemia inhibitory factor and pentraxin 3 expression appeared to be associated with shorter survival in primary and recurrent GBM using in silico data. In primary GBM, gene set analysis also showed higher expression of genes involved in metabolism, extracellular matrix remodeling and complement activation, whereas genes involved in T cell activation and checkpoint signaling were expressed at higher levels in recurrent GBM. The reduced level of pentraxin 3 in recurrent glioblastomas and the gene set analysis results suggest an altered microenvironment in recurrent GBM that might be more active.

Advancements in Neuroimaging to Unravel Biological and Molecular Features of Brain Tumors

Simple Summary

Advanced neuroimaging is gaining increasing relevance for the characterization and the molecular profiling of brain tumor tissue. On one hand, for some tumor types, the most widespread advanced techniques, investigating diffusion and perfusion features, have been proven clinically feasible and rather robust for diagnosis and prognosis stratification. In addition, 2-hydroxyglutarate spectroscopy, for the first time, offers the possibility to directly measure a crucial molecular marker. On the other hand, numerous innovative approaches have been explored for a refined evaluation of tumor microenvironments, particularly assessing microstructural and microvascular properties, and the potential applications of these techniques are vast and still to be fully explored.

Abstract

In recent years, the clinical assessment of primary brain tumors has been increasingly dependent on advanced magnetic resonance imaging (MRI) techniques in order to infer tumor pathophysiological characteristics, such as hemodynamics, metabolism, and microstructure. Quantitative radiomic data extracted from advanced MRI have risen as potential in vivo noninvasive biomarkers for predicting tumor grades and molecular subtypes, opening the era of “molecular imaging” and radiogenomics. This review presents the most relevant advancements in quantitative neuroimaging of advanced MRI techniques, by means of radiomics analysis, applied to primary brain tumors, including lower-grade glioma and glioblastoma, with a special focus on peculiar oncologic entities of current interest. Novel findings from diffusion MRI (dMRI), perfusion-weighted imaging (PWI), and MR spectroscopy (MRS) are hereby sifted in order to evaluate the role of quantitative imaging in neuro-oncology as a tool for predicting molecular profiles, stratifying prognosis, and characterizing tumor tissue microenvironments. Furthermore, innovative technological approaches are briefly addressed, including artificial intelligence contributions and ultra-high-field imaging new techniques. Lastly, after providing an overview of the advancements, we illustrate current clinical applications and future perspectives.

Tissue Hypoxia and Alterations in Microvascular Architecture Predict Glioblastoma Recurrence in Humans

PURPOSE

Insufficient control of infiltrative glioblastoma (GBM) cells is a major cause of treatment failure and tumor recurrence. Hence, detailed insights into pathophysiologic changes that precede GBM recurrence are needed to develop more precise neuroimaging modalities for tailored diagnostic monitoring and therapeutic approaches.

EXPERIMENTAL DESIGN

Overall, 168 physiologic MRI follow-up examinations of 56 patients with GBM who developed recurrence after standard therapy were retrospectively evaluated, that is, two post-standard-therapeutic follow-ups before and one at radiological recurrence. MRI biomarkers for microvascular architecture and perfusion, neovascularization activity, oxygen metabolism, and hypoxia were determined for brain areas that developed in the further course into recurrence and for the recurrent GBM itself. The temporal pattern of biomarker changes was fitted with locally estimated scatterplot smoothing functions and analyzed for pathophysiologic changes preceding radiological GBM recurrence.

RESULTS

Our MRI approach demonstrated early pathophysiologic changes prior to radiological GBM recurrence in all patients. Analysis of the time courses revealed a model for the pathophysiology of GBM recurrence: 190 days prior to radiological recurrence, vascular cooption by GBM cells induced vessel regression, detected as decreasing vessel density/perfusion and increasing hypoxia. Seventy days later, neovascularization activity was upregulated, which reincreased vessel density and perfusion. Hypoxia, however, continued to intensify for 30 days and peaked 90 days before radiological recurrence.

CONCLUSIONS

Hypoxia may represent an early sign for GBM recurrence. This might become useful in the development of new combined diagnostic-therapeutic approaches for tailored clinical management of recurrent GBM. Further preclinical and in-human studies are required for validation and evaluation.

Emerging Techniques in Imaging of Glioma Microenvironment

Magnetic resonance imaging (MRI) has been the cornerstone of imaging of brain tumors in the past 4 decades. Conventional MRI remains the workhorse for neuro-oncologic imaging, not only for basic information such as location, extent, and navigation but also able to provide information regarding proliferation and infiltration, angiogenesis, hemorrhage, and more. More sophisticated MRI sequences have extended the ability to assess and quantify these features; for example, permeability and perfusion acquisitions can assess blood-brain barrier disruption and angiogenesis, diffusion techniques can assess cellularity and infiltration, and spectroscopy can address metabolism. Techniques such as fMRI and diffusion fiber tracking can be helpful in diagnostic planning for resection and radiation therapy, and more sophisticated iterations of these techniques can extend our understanding of neurocognitive effects of these tumors and associated treatment responses and effects. More recently, MRI has been used to go beyond such morphological, physiological, and functional characteristics to assess the tumor microenvironment. The current review highlights multiple recent and emerging approaches in MRI to characterize the tumor microenvironment.

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.

Development of a Non-invasive Assessment of Hypoxia and Neovascularization with Magnetic Resonance Imaging in Benign and Malignant Breast Tumors: Initial Results

PURPOSE

To develop a novel magnetic resonance imaging (MRI) approach for the noninvasive assessment of hypoxia and neovascularization in breast tumors.

PROCEDURES

In this IRB-approved prospective study, 20 patients with suspicious breast lesions (BI-RADS 4/5) underwent multiparametric breast MRI including quantitative BOLD (qBOLD) and vascular architecture mapping (VAM). Custom-made in-house MatLab software was used for qBOLD and VAM data postprocessing and calculation of quantitative MRI biomarker maps of oxygen extraction fraction (OEF), metabolic rate of oxygen (MRO₂), and mitochondrial oxygen tension (mitoPO₂) to measure tissue hypoxia and neovascularization including vascular architecture including microvessel radius (VSI), density (MVD), and type (MTI). Histopathology was used as standard of reference. Appropriate statistics were performed to assess and compare correlations between MRI biomarkers for hypoxia and neovascularization.

RESULTS

qBOLD and VAM data with good quality were obtained from all patients with 13 invasive ductal carcinoma (IDC) and 7 benign breast tumors with a lesion diameter of at least 10 mm in all spatial directions. MRI biomarker maps of oxygen metabolism and neovascularization demonstrated intratumoral spatial heterogeneity with a broad range of biomarker values. Bulk tumor neovasculature consisted of draining venous microvasculature with slow flowing blood. High OEF and low mitoPO₂ were associated with low MVD and vice versa. The heterogeneous pattern of MRO₂ values showed spatial congruence with VSI. IDCs showed significantly higher MRO₂ (P = 0.007), lower mitoPO₂ (P = 0.021), higher MVD (P = 0.005), and lower (i.e., more pathologic) MTI (P = 0.001) compared with benign breast tumors. These results indicate that IDCs consume more oxygen and are more hypoxic and neovascularized than benign tumors.

CONCLUSIONS

We developed a novel MRI approach for the noninvasive assessment of hypoxia and neovascularization in benign and malignant breast tumors that can be easily integrated in a diagnostic MRI protocol and provides insight into intratumoral heterogeneity.

Intratumoral heterogeneity of oxygen metabolism and neovascularization uncovers 2 survival-relevant subgroups of IDH1 wild-type glioblastoma

Background

The intratumoral heterogeneity of oxygen metabolism in combination with variable patterns of neovascularization (NV) as well as reprogramming of energy metabolism affects the landscape of tumor microenvironments (TMEs) in glioblastoma. Knowledge of the hypoxic and perivascular niches within the TME is essential for understanding treatment failure.

Methods

Fifty-two patients with untreated glioblastoma (isocitrate dehydrogenase 1 wild type [IDH1wt]) were examined with a physiological MRI protocol including a multiparametric quantitative blood oxygen level dependent (qBOLD) approach and vascular architecture mapping (VAM). Imaging biomarker information about oxygen metabolism (mitochondrial oxygen tension) and neovascularization (microvascular density and type) were fused for classification of 6 different TMEs: necrosis, hypoxia with/without neovascularization, oxidative phosphorylation (OxPhos), and glycolysis with/without neovascularization. Association of the different TME volume fractions with progression-free survival (PFS) was assessed using Kaplan-Meier analysis and Cox proportional hazards models.

Results

A common spatial structure of TMEs was detected: central necrosis surrounded by tumor hypoxia (with defective and functional neovasculature) and different TMEs with a predominance of OxPhos and glycolysis for energy production, respectively. The percentage of the different TMEs on the total tumor volume uncovered 2 clearly different subtypes of glioblastoma IDH1wt: a glycolytic dominated phenotype with predominantly functional neovasculature and a necrotic/hypoxic dominated phenotype with approximately 50% of defective neovasculature. Patients with a necrotic/hypoxic dominated phenotype showed significantly shorter PFS (P = 0.035).

Conclusions

Our non-invasive mapping approach allows for classification of the TME and detection of tumor-supportive niches in glioblastoma which may be helpful for both clinical patient management and research.

Gene Suppression of Transketolase-Like Protein 1 (TKTL1) Sensitizes Glioma Cells to Hypoxia and Ionizing Radiation

In several tumor entities, transketolase-like protein 1 (TKTL1) has been suggested to promote the nonoxidative part of the pentose phosphate pathway (PPP) and thereby to contribute to a malignant phenotype. However, its role in glioma biology has only been sparsely documented. In the present in vitro study using LNT-229 glioma cells, we analyzed the impact of TKTL1 gene suppression on basic metabolic parameters and on survival following oxygen restriction and ionizing radiation. TKTL1 was induced by hypoxia and by hypoxia-inducible factor-1α (HIF-1α). Knockdown of TKTL1 via shRNA increased the cells’ demand for glucose, decreased flux through the PPP and promoted cell death under hypoxic conditions. Following irradiation, suppression of TKTL1 expression resulted in elevated levels of reactive oxygen species (ROS) and reduced clonogenic survival. In summary, our results indicate a role of TKTL1 in the adaptation of tumor cells to oxygen deprivation and in the acquisition of radioresistance. Further studies are necessary to examine whether strategies that antagonize TKTL1 function will be able to restore the sensitivity of glioma cells towards irradiation and antiangiogenic therapies in the more complex in vivo environment.

Vascular Hysteresis Loops and Vascular Architecture Mapping in Patients with Glioblastoma treated with Antiangiogenic Therapy

In this study, we investigated the variability of vascular hysteresis loop (VHL) shapes and the spatial heterogeneity of neovascularization and microvascular alterations using vascular architecture mapping (VAM) in patients with recurrent glioblastoma during bevacizumab mono-therapy. VAM data were acquired in 13 patients suffering from recurrent glioblastoma prior to and 3 months after bevacizumab treatment onset using a dual contrast agent injections approach as part of routine MRI. Two patients were additionally examined after the first cycle of bevacizumab to check for early treatment response. VHLs were evaluated as biomarker maps of neovascularization activity: microvessel type indicator (MTI) and curvature (Curv) of the VHL-long-axis. Early response to bevacizumab was dominated by reduction of smaller microvasculature (around 10 µm). In the 3-month follow-up, responding tumors additionally showed a reduction in larger microvasculature (>20 µm). VAM biomarker images revealed spatially heterogeneous microvascular alterations during bevacizumab treatment. Responding, non-responding, progressive, and remote-progressive tumor areas were observed. MTI may be useful to predict responding and non-responding tumor regions, and Curv to assess severity of vasogenic edema. Analysis of VHLs in combination with VAM biomarkers may lead to a new perspective on investigating the spatial heterogeneity of neovascularization and microvascular alterations in glioblastoma during antiangiogenic therapy.