Abnormal vascular structure and function within brain metastases is linked to pembrolizumab resistance.

BACKGROUND: We recently conducted a phase 2 trial (NCT028865685) evaluating intracranial efficacy of pembrolizumab for brain metastases (BM) of diverse histologies. Our study met its primary efficacy endpoint and illustrates that pembrolizumab exerts promising activity in a select group of patients with BM. Given the importance of aberrant vasculature in mediating immunosuppression, we explored the relationship between checkpoint inhibitor (ICI) efficacy and vascular architecture in the hopes of identifying potential mechanisms of intracranial ICI response or resistance for BM. METHODS: Using Vessel Architectural Imaging (VAI), a histologically validated quantitative metric for in vivo tumor vascular physiology, we analyzed dual echo DSC/DCE MRI for 44 patients on trial. Tumor and peri-tumor cerebral blood volume/flow, vessel size, arterial- and venous-dominance, and vascular permeability were measured before and after treatment with pembrolizumab. RESULTS: BM that progressed on ICI were characterized by a highly aberrant vasculature dominated by large-caliber vessels. In contrast, ICI-responsive BM possessed a more structurally balanced vasculature consisting of both small and large vessels, and there was a trend towards a decrease in under-perfused tissue, suggesting a reversal of the negative effects of hypoxia. In the peri-tumor region, development of smaller blood vessels, consistent with neo-angiogenesis, was associated with tumor growth before radiographic evidence of contrast enhancement on anatomical MRI. CONCLUSIONS: This study, one of the largest functional imaging studies for BM, suggests that vascular architecture is linked with ICI efficacy. Studies identifying modulators of vascular architecture, and effects on immune activity, are warranted and may inform future combination treatments.

Interpreting deep learning models for glioma survival classification using visualization and textual explanations.

BACKGROUND: Saliency-based algorithms are able to explain the relationship between input image pixels and deep-learning model predictions. However, it may be difficult to assess the clinical value of the most important image features and the model predictions derived from the raw saliency map. This study proposes to enhance the interpretability of saliency-based deep learning model for survival classification of patients with gliomas, by extracting domain knowledge-based information from the raw saliency maps. MATERIALS AND METHODS: Our study includes presurgical T1-weighted (pre- and post-contrast), T2-weighted and T2-FLAIR MRIs of 147 glioma patients from the BraTs 2020 challenge dataset aligned to the SRI 24 anatomical atlas. Each image exam includes a segmentation mask and the information of overall survival (OS) from time of diagnosis (in days). This dataset was divided into training ([Formula: see text]) and validation ([Formula: see text]) datasets. The extent of surgical resection for all patients was gross total resection. We categorized the data into 42 short (mean [Formula: see text] days), 30 medium ([Formula: see text] days), and 46 long ([Formula: see text] days) survivors. A 3D convolutional neural network (CNN) trained on brain tumour MRI volumes classified all patients based on expected prognosis of either short-term, medium-term, or long-term survival. We extend the popular 2D Gradient-weighted Class Activation Mapping (Grad-CAM), for the generation of saliency map, to 3D and combined it with the anatomical atlas, to extract brain regions, brain volume and probability map that reveal domain knowledge-based information. RESULTS: For each OS class, a larger tumor volume was associated with a shorter OS. There were 10, 7 and 27 tumor locations in brain regions that uniquely associate with the short-term, medium-term, and long-term survival, respectively. Tumors located in the transverse temporal gyrus, fusiform, and palladium are associated with short, medium and long-term survival, respectively. The visual and textual information displayed during OS prediction highlights tumor location and the contribution of different brain regions to the prediction of OS. This algorithm design feature assists the physician in analyzing and understanding different model prediction stages. CONCLUSIONS: Domain knowledge-based information extracted from the saliency map can enhance the interpretability of deep learning models. Our findings show that tumors overlapping eloquent brain regions are associated with short patient survival.

Structural and functional vascular dysfunction within brain metastases is linked to pembrolizumab inefficacy.

Structurally and functionally aberrant vasculature is a hallmark of tumor angiogenesis and treatment resistance. Given the synergistic link between aberrant tumor vasculature and immunosuppression, we analyzed perfusion MRI for 44 patients with brain metastases (BM) undergoing treatment with pembrolizumab. To date, vascular-immune communication, or the relationship between immune checkpoint inhibitor (ICI) efficacy and vascular architecture, has not been well-characterized in human imaging studies. We found that ICI-responsive BM possessed a structurally balanced vascular makeup, which was linked to improved vascular efficiency and an immune-stimulatory microenvironment. In contrast, ICI-resistant BM were characterized by a lack of immune cell infiltration and a highly aberrant vasculature dominated by large-caliber vessels. Peri-tumor region analysis revealed early functional changes predictive of ICI resistance before radiographic evidence on conventional MRI. This study was one of the largest functional imaging studies for BM and establishes a foundation for functional studies that illuminate the mechanisms linking patterns of vascular architecture with immunosuppression, as targeting these aspects of cancer biology may serve as the basis for future combination treatments.

Vascular differences between IDH-wildtype glioblastoma and astrocytoma IDH-mutant grade 4 at imaging and transcriptomic levels.

Global agreement in central nervous system (CNS) tumor classification is essential for predicting patient prognosis and determining the correct course of treatment, as well as for stratifying patients for clinical trials at international level. The last update by the World Health Organization of CNS tumor classification and grading in 2021 considered, for the first time, IDH-wildtype glioblastoma and astrocytoma IDH-mutant grade 4 as different tumors. Mutations in the genes isocitrate dehydrogenase (IDH) 1 and 2 occur early and, importantly, contribute to gliomagenesis. IDH mutation produces a metabolic reprogramming of tumor cells, thus affecting the processes of hypoxia and vascularity, resulting in a clear advantage for those patients who present with IDH-mutated astrocytomas. Despite the clinical relevance of IDH mutation, current protocols do not include full sequencing for every patient. Alternative biomarkers could be useful and complementary to obtain a more reliable classification. In this sense, magnetic resonance imaging (MRI)-perfusion biomarkers, such as relative cerebral blood volume and flow, could be useful from the moment of presurgery, without incurring additional financial costs or requiring extra effort. The main purpose of this work is to analyze the vascular and hemodynamic differences between IDH-wildtype glioblastoma and IDH-mutant astrocytoma. To achieve this, we evaluate and validate the association between dynamic susceptibility contrast-MRI perfusion biomarkers and IDH mutation status. In addition, to gain a deeper understanding of the vascular differences in astrocytomas depending on the IDH mutation, we analyze the transcriptomic bases of the vascular differences.

MR elastography identifies regions of extracellular matrix reorganization associated with shorter survival in glioblastoma patients.

Background: Biomechanical tissue properties of glioblastoma tumors are heterogeneous, but the molecular mechanisms involved and the biological implications are poorly understood. Here, we combine magnetic resonance elastography (MRE) measurement of tissue stiffness with RNA sequencing of tissue biopsies to explore the molecular characteristics of the stiffness signal. Methods: MRE was performed preoperatively in 13 patients with glioblastoma. Navigated biopsies were harvested during surgery and classified as "stiff" or "soft" according to MRE stiffness measurements (|G*|norm). Twenty-two biopsies from eight patients were analyzed by RNA sequencing. Results: The mean whole-tumor stiffness was lower than normal-appearing white matter. The surgeon's stiffness evaluation did not correlate with the MRE measurements, which suggests that these measures assess different physiological properties. Pathway analysis of the differentially expressed genes between "stiff" and "soft" biopsies showed that genes involved in extracellular matrix reorganization and cellular adhesion were overexpressed in "stiff" biopsies. Supervised dimensionality reduction identified a gene expression signal separating "stiff" and "soft" biopsies. Using the NIH Genomic Data Portal, 265 glioblastoma patients were divided into those with (n = 63) and without (n = 202) this gene expression signal. The median survival time of patients with tumors expressing the gene signal associated with "stiff" biopsies was 100 days shorter than that of patients not expressing it (360 versus 460 days, hazard ratio: 1.45, P < .05). Conclusion: MRE imaging of glioblastoma can provide noninvasive information on intratumoral heterogeneity. Regions of increased stiffness were associated with extracellular matrix reorganization. An expression signal associated with "stiff" biopsies correlated with shorter survival of glioblastoma patients.

Advanced MR Techniques for Preoperative Glioma Characterization: Part 1.

Preoperative clinical magnetic resonance imaging (MRI) protocols for gliomas, brain tumors with dismal outcomes due to their infiltrative properties, still rely on conventional structural MRI, which does not deliver information on tumor genotype and is limited in the delineation of diffuse gliomas. The GliMR COST action wants to raise awareness about the state of the art of advanced MRI techniques in gliomas and their possible clinical translation or lack thereof. This review describes current methods, limits, and applications of advanced MRI for the preoperative assessment of glioma, summarizing the level of clinical validation of different techniques. In this first part, we discuss dynamic susceptibility contrast and dynamic contrast-enhanced MRI, arterial spin labeling, diffusion-weighted MRI, vessel imaging, and magnetic resonance fingerprinting. The second part of this review addresses magnetic resonance spectroscopy, chemical exchange saturation transfer, susceptibility-weighted imaging, MRI-PET, MR elastography, and MR-based radiomics applications. Evidence Level: 3 Technical Efficacy: Stage 2.

Advanced MR Techniques for Preoperative Glioma Characterization: Part 2.

Preoperative clinical MRI protocols for gliomas, brain tumors with dismal outcomes due to their infiltrative properties, still rely on conventional structural MRI, which does not deliver information on tumor genotype and is limited in the delineation of diffuse gliomas. The GliMR COST action wants to raise awareness about the state of the art of advanced MRI techniques in gliomas and their possible clinical translation. This review describes current methods, limits, and applications of advanced MRI for the preoperative assessment of glioma, summarizing the level of clinical validation of different techniques. In this second part, we review magnetic resonance spectroscopy (MRS), chemical exchange saturation transfer (CEST), susceptibility-weighted imaging (SWI), MRI-PET, MR elastography (MRE), and MR-based radiomics applications. The first part of this review addresses dynamic susceptibility contrast (DSC) and dynamic contrast-enhanced (DCE) MRI, arterial spin labeling (ASL), diffusion-weighted MRI, vessel imaging, and magnetic resonance fingerprinting (MRF). EVIDENCE LEVEL: 3. TECHNICAL EFFICACY: Stage 2.

Unlocking Bevacizumab's Potential: rCBVmax as a Predictive Biomarker for Enhanced Survival in Glioblastoma IDH-Wildtype Patients.

BACKGROUND: Aberrant vascular architecture and angiogenesis are hallmarks of glioblastoma IDH-wildtype, suggesting that these tumors are suitable for antiangiogenic therapy. Bevacizumab was FDA-approved in 2009 following promising results in two clinical trials. However, its use for recurrent glioblastomas remains a subject of debate, as it does not universally improve patient survival. PURPOSES: In this study, we aimed to analyze the influence of tumor vascularity on the benefit provided by BVZ and propose preoperative rCBVmax at the high angiogenic tumor habitat as a predictive biomarker to select patients who can benefit the most. METHODS: Clinical and MRI data from 106 patients with glioblastoma IDH-wildtype have been analyzed. Thirty-nine of them received BVZ, and the remaining sixty-seven did not receive a second-line treatment. The ONCOhabitats method was used to automatically calculate rCBV. RESULTS: We found a median survival from progression of 305 days longer for patients with moderate vascular tumors who received BVZ than those who did not receive any second-line treatment. This contrasts with patients with high-vascular tumors who only presented a median survival of 173 days longer when receiving BVZ. Furthermore, better responses to BVZ were found for the moderate-vascular group with a higher proportion of patients alive at 6, 12, 18, and 24 months after progression. CONCLUSIONS: We propose rCBVmax as a potential biomarker to select patients who can benefit more from BVZ after tumor progression. In addition, we propose a threshold of 7.5 to stratify patients into moderate- and high-vascular groups to select the optimal second-line treatment.

Quantification of Tissue Compression Identifies High-Grade Glioma Patients with Reduced Survival.

The compression of peritumoral healthy tissue in brain tumor patients is considered a major cause of the life-threatening neurologic symptoms. Although significant deformations caused by the tumor growth can be observed radiologically, the quantification of minor tissue deformations have not been widely investigated. In this study, we propose a method to quantify subtle peritumoral deformations. A total of 127 MRI longitudinal studies from 23 patients with high-grade glioma were included. We estimate longitudinal displacement fields based on a symmetric normalization algorithm and we propose four biomarkers. We assess the interpatient and intrapatient association between proposed biomarkers and the survival based on Cox analyses, and the potential of the biomarkers to stratify patients according to their survival based on Kaplan−Meier analysis. Biomarkers show a significant intrapatient association with survival (p < 0.05); however, only compression biomarkers show the ability to stratify patients between those with higher and lower overall survival (AUC = 0.83, HR = 6.30, p < 0.05 for CompCH). The compression biomarkers present three times higher Hazard Ratios than those representing only displacement. Our study provides a robust and automated method for quantifying and delineating compression in the peritumoral area. Based on the proposed methodology, we found an association between lower compression in the peritumoral area and good prognosis in high-grade glial tumors.

Decreased tissue stiffness in glioblastoma by MR elastography is associated with increased cerebral blood flow.

PURPOSE: Understanding how mechanical properties relate to functional changes in glioblastomas may help explain different treatment response between patients. The aim of this study was to map differences in biomechanical and functional properties between tumor and healthy tissue, to assess any relationship between them and to study their spatial distribution. METHODS: Ten patients with glioblastoma and 17 healthy subjects were scanned using MR Elastography, perfusion and diffusion MRI. Stiffness and viscosity measurements G' and G'', cerebral blood flow (CBF), apparent diffusion coefficient (ADC) and fractional anisotropy (FA) were measured in patients' contrast-enhancing tumor, necrosis, edema, and gray and white matter, and in gray and white matter for healthy subjects. A regression analysis was used to predict CBF as a function of ADC, FA, G' and G''. RESULTS: Median G' and G'' in contrast-enhancing tumor were 13% and 37% lower than in normal-appearing white matter (P < 0.01), and 8% and 6% lower in necrosis than in contrast-enhancing tumor, respectively (P < 0.05). Tumors showed both inter-patient and intra-patient heterogeneity. Measurements approached values in normal-appearing tissue when moving outward from the tumor core, but abnormal tissue properties were still present in regions of normal-appearing tissue. Using both a linear and a random-forest model, prediction of CBF was improved by adding MRE measurements to the model (P < 0.01). CONCLUSIONS: The inclusion of MRE measurements in statistical models helped predict perfusion, with stiffer tissue associated with lower perfusion values.

Local detection of microvessels in IDH-wildtype glioblastoma using relative cerebral blood volume: an imaging marker useful for astrocytoma grade 4 classification.

BACKGROUND: The microvessels area (MVA), derived from microvascular proliferation, is a biomarker useful for high-grade glioma classification. Nevertheless, its measurement is costly, labor-intense, and invasive. Finding radiologic correlations with MVA could provide a complementary non-invasive approach without an extra cost and labor intensity and from the first stage. This study aims to correlate imaging markers, such as relative cerebral blood volume (rCBV), and local MVA in IDH-wildtype glioblastoma, and to propose this imaging marker as useful for astrocytoma grade 4 classification. METHODS: Data from 73 tissue blocks belonging to 17 IDH-wildtype glioblastomas and 7 blocks from 2 IDH-mutant astrocytomas were compiled from the Ivy GAP database. MRI processing and rCBV quantification were carried out using ONCOhabitats methodology. Histologic and MRI co-registration was done manually with experts' supervision, achieving an accuracy of 88.8% of overlay. Spearman's correlation was used to analyze the association between rCBV and microvessel area. Mann-Whitney test was used to study differences of rCBV between blocks with presence or absence of microvessels in IDH-wildtype glioblastoma, as well as to find differences with IDH-mutant astrocytoma samples. RESULTS: Significant positive correlations were found between rCBV and microvessel area in the IDH-wildtype blocks (p < 0.001), as well as significant differences in rCBV were found between blocks with microvascular proliferation and blocks without it (p < 0.0001). In addition, significant differences in rCBV were found between IDH-wildtype glioblastoma and IDH-mutant astrocytoma samples, being 2-2.5 times higher rCBV values in IDH-wildtype glioblastoma samples. CONCLUSIONS: The proposed rCBV marker, calculated from diagnostic MRIs, can detect in IDH-wildtype glioblastoma those regions with microvessels from those without it, and it is significantly correlated with local microvessels area. In addition, the proposed rCBV marker can differentiate the IDH mutation status, providing a complementary non-invasive method for high-grade glioma classification.

Glioblastoma versus solitary brain metastasis: MRI differentiation using the edema perfusion gradient.

BACKGROUND AND PURPOSE: Differentiation between glioblastoma multiforme (GBM) and solitary brain metastasis (SBM) remains a challenge in neuroradiology with up to 40% of the cases to be incorrectly classified using only conventional MRI. The inclusion of perfusion MRI parameters provides characteristic features that could support the distinction of these pathological entities. On these grounds, we aim to use a perfusion gradient in the peritumoral edema. METHODS: Twenty-four patients with GBM or an SBM underwent conventional and perfusion MR imaging sequences before tumors' surgical resection. After postprocessing of the images, quantification of dynamic susceptibility contrast (DSC) perfusion parameters was made. Three concentric areas around the tumor were defined in each case. The monocompartimental and pharmacokinetics parameters of perfusion MRI were analyzed in both series. RESULTS: DSC perfusion MRI models can provide useful information for the differentiation between GBM and SBM. It can be observed that most of the perfusion MR parameters (relative cerebral blood volume, relative cerebral blood flow, relative Ktrans, and relative volume fraction of the interstitial space) clearly show higher gradient for GBM than SBM. GBM also demonstrates higher heterogeneity in the peritumoral edema and most of the perfusion parameters demonstrate higher gradients in the area closest to the enhancing tumor. CONCLUSION: Our results show that there is a difference in the perfusion parameters of the edema between GBM and SBM demonstrating a vascularization gradient. This could help not only for the diagnosis, but also for planning surgical or radiotherapy treatments delineating the real extension of the tumor.

Lack of Benefit of Extending Temozolomide Treatment in Patients with High Vascular Glioblastoma with Methylated MGMT.

In this study, we evaluated the benefit on survival of the combination of methylation of O6-methylguanine-DNA methyltransferase (MGMT) promotor gene and moderate vascularity in glioblastoma using a retrospective dataset of 123 patients from a multicenter cohort. MRI processing and calculation of relative cerebral blood volume (rCBV), used to define moderate- and high-vascular groups, were performed with the automatic ONCOhabitats method. We assessed the previously proposed rCBV threshold (10.7) and the new calculated ones (9.1 and 9.8) to analyze the association with survival for different populations according to vascularity and MGMT methylation status. We found that patients included in the moderate-vascular group had longer survival when MGMT is methylated (significant median survival difference of 174 days, p = 0.0129*). However, we did not find significant differences depending on the MGMT methylation status for the high-vascular group (p = 0.9119). In addition, we investigated the combined correlation of MGMT methylation status and rCBV with the prognostic effect of the number of temozolomide cycles, and only significant results were found for the moderate-vascular group. In conclusion, there is a lack of benefit of extending temozolomide treatment for patients with high vascular glioblastomas, even presenting MGMT methylation. Preliminary results suggest that patients with moderate vascularity and methylated MGMT glioblastomas would benefit more from prolonged adjuvant chemotherapy.

Differential effect of vascularity between long- and short-term survivors with IDH1/2 wild-type glioblastoma.

INTRODUCTION: IDH1/2 wt glioblastoma (GB) represents the most lethal tumour of the central nervous system. Tumour vascularity is associated with overall survival (OS), and the clinical relevance of vascular markers, such as rCBV, has already been validated. Nevertheless, molecular and clinical factors may have different influences on the beneficial effect of a favourable vascular signature. PURPOSE: To evaluate the association between the rCBV and OS of IDH1/2 wt GB patients for long-term survivors (LTSs) and short-term survivors (STSs). Given that initial high rCBV may affect the patient's OS in follow-up stages, we will assess whether a moderate vascularity is beneficial for OS in both groups of patients. MATERIALS AND METHODS: Ninety-nine IDH1/2 wt GB patients were divided into LTSs (OS ≥ 400 days) and STSs (OS < 400 days). Mann-Whitney and Fisher, uni- and multiparametric Cox, Aalen's additive regression and Kaplan-Meier tests were carried out. Tumour vascularity was represented by the mean rCBV of the high angiogenic tumour (HAT) habitat computed through the haemodynamic tissue signature methodology (available on the ONCOhabitats platform). RESULTS: For LTSs, we found a significant association between a moderate value of rCBVmean and higher OS (uni- and multiparametric Cox and Aalen's regression) (p = 0.0140, HR = 1.19; p = 0.0085, HR = 1.22) and significant stratification capability (p = 0.0343). For the STS group, no association between rCBVmean and survival was observed. Moreover, no significant differences (p > 0.05) in gender, age, resection status, chemoradiation, or MGMT methylation were observed between LTSs and STSs. CONCLUSION: We have found different prognostic and stratification effects of the vascular marker for the LTS and STS groups. We propose the use of rCBVmean at HAT as a vascular marker clinically relevant for LTSs with IDH1/2 wt GB and maybe as a potential target for randomized clinical trials focused on this group of patients.

MGMT methylation may benefit overall survival in patients with moderately vascularized glioblastomas.

OBJECTIVES: To assess the combined role of tumor vascularity, estimated from perfusion MRI, and MGMT methylation status on overall survival (OS) in patients with glioblastoma. METHODS: A multicentric international dataset including 96 patients from NCT03439332 clinical study were used to study the prognostic relationships between MGMT and perfusion markers. Relative cerebral blood volume (rCBV) in the most vascularized tumor regions was automatically obtained from preoperative MRIs using ONCOhabitats online analysis service. Cox survival regression models and stratification strategies were conducted to define a subpopulation that is particularly favored by MGMT methylation in terms of OS. RESULTS: rCBV distributions did not differ significantly (p > 0.05) in the methylated and the non-methylated subpopulations. In patients with moderately vascularized tumors (rCBV < 10.73), MGMT methylation was a positive predictive factor for OS (HR = 2.73, p = 0.003, AUC = 0.70). In patients with highly vascularized tumors (rCBV > 10.73), however, there was no significant effect of MGMT methylation (HR = 1.72, p = 0.10, AUC = 0.56). CONCLUSIONS: Our results indicate the existence of complementary prognostic information provided by MGMT methylation and rCBV. Perfusion markers could identify a subpopulation of patients who will benefit the most from MGMT methylation. Not considering this information may lead to bias in the interpretation of clinical studies. KEY POINTS: • MRI perfusion provides complementary prognostic information to MGMT methylation. • MGMT methylation improves prognosis in glioblastoma patients with moderate vascular profile. • Failure to consider these relations may lead to bias in the interpretation of clinical studies.

Higher vascularity at infiltrated peripheral edema differentiates proneural glioblastoma subtype.

BACKGROUND AND PURPOSE: Genetic classifications are crucial for understanding the heterogeneity of glioblastoma. Recently, perfusion MRI techniques have demonstrated associations molecular alterations. In this work, we investigated whether perfusion markers within infiltrated peripheral edema were associated with proneural, mesenchymal, classical and neural subtypes. MATERIALS AND METHODS: ONCOhabitats open web services were used to obtain the cerebral blood volume at the infiltrated peripheral edema for MRI studies of 50 glioblastoma patients from The Cancer Imaging Archive: TCGA-GBM. ANOVA and Kruskal-Wallis tests were carried out in order to assess the association between vascular features and the Verhaak subtypes. For assessing specific differences, Mann-Whitney U-test was conducted. Finally, the association of overall survival with molecular and vascular features was assessed using univariate and multivariate Cox models. RESULTS: ANOVA and Kruskal-Wallis tests for the maximum cerebral blood volume at the infiltrated peripheral edema between the four subclasses yielded false discovery rate corrected p-values of <0.001 and 0.02, respectively. This vascular feature was significantly higher (p = 0.0043) in proneural patients compared to the rest of the subtypes while conducting Mann-Whitney U-test. The multivariate Cox model pointed to redundant information provided by vascular features at the peripheral edema and proneural subtype when analyzing overall survival. CONCLUSIONS: Higher relative cerebral blood volume at infiltrated peripheral edema is associated with proneural glioblastoma subtype suggesting underlying vascular behavior related to molecular composition in that area.

The impact of EPI-based distortion correction of dynamic susceptibility contrast MRI on cerebral blood volume estimation in patients with glioblastoma.

PURPOSE: Relative cerebral blood volume (rCBV) from dynamic susceptibility contrast (DSC)-MRI is a valuable biomarker in patients with glioblastoma for assessing treatment response and predicting overall survival. DSC-MRI based on echo planar images (EPI) may possess severe geometric distortions from magnetic field inhomogeneities up to the order of centimeters. The aim of this study is to assess how much two readily available EPI-based geometric distortion correction methods, FSL TOPUP and EPIC, affect rCBV values from DSC-MRI in patients with confirmed glioblastoma. METHOD: We used a combined single-shot 2D gradient-echo (T2*), spin-echo (T2) EPI sequence to estimate both T2* and T2-weighted rCBV from the same contrast agent injection. Effects of distortion correction on the positive phase-encoded T2- and T2*-images were assessed in healthy anatomical brain regions in terms of Wilcoxon signed rank tests on median rCBV change and on Dice coefficients, as well as in tumor lesions in terms of Wilcoxon signed rank tests on median rCBV change. RESULTS: Our results show that following distortion correction, both gradient-echo and spin-echo rCBV increased in cortical areas of the frontal, temporal and occipital lobe, including the posterior orbital gyri in the frontal lobe and middle frontal gyri (p < 0.0008). Similar, improved Dice coefficients were observed for gradient-echo EPI in temporal, occipital and frontal lobe. Only spin-echo rCBV in enhancing lesion increased with correction (p = 0.0002). CONCLUSION: Our study sheds light on the importance of performing geometric distortion correction on EPI-based MRI data before assessing functional information such as rCBV values. Our findings may indicate that uncorrected rCBV values can be underestimated from positive phase-encoding EPI and that geometric distortion correction is warranted when comparing EPI-based data to conventional MRI.

Robust association between vascular habitats and patient prognosis in glioblastoma: An international multicenter study.

BACKGROUND: Glioblastoma (GBM) is the most aggressive primary brain tumor, characterized by a heterogeneous and abnormal vascularity. Subtypes of vascular habitats within the tumor and edema can be distinguished: high angiogenic tumor (HAT), low angiogenic tumor (LAT), infiltrated peripheral edema (IPE), and vasogenic peripheral edema (VPE). PURPOSE: To validate the association between hemodynamic markers from vascular habitats and overall survival (OS) in glioblastoma patients, considering the intercenter variability of acquisition protocols. STUDY TYPE: Multicenter retrospective study. POPULATION: In all, 184 glioblastoma patients from seven European centers participating in the NCT03439332 clinical study. FIELD STRENGTH/SEQUENCE: 1.5T (for 54 patients) or 3.0T (for 130 patients). Pregadolinium and postgadolinium-based contrast agent-enhanced T1 -weighted MRI, T2 - and FLAIR T2 -weighted, and dynamic susceptibility contrast (DSC) T2 * perfusion. ASSESSMENT: We analyzed preoperative MRIs to establish the association between the maximum relative cerebral blood volume (rCBVmax ) at each habitat with OS. Moreover, the stratification capabilities of the markers to divide patients into "vascular" groups were tested. The variability in the markers between individual centers was also assessed. STATISTICAL TESTS: Uniparametric Cox regression; Kaplan-Meier test; Mann-Whitney test. RESULTS: The rCBVmax derived from the HAT, LAT, and IPE habitats were significantly associated with patient OS (P < 0.05; hazard ratio [HR]: 1.05, 1.11, 1.28, respectively). Moreover, these markers can stratify patients into "moderate-" and "high-vascular" groups (P < 0.05). The Mann-Whitney test did not find significant differences among most of the centers in markers (HAT: P = 0.02-0.685; LAT: P = 0.010-0.769; IPE: P = 0.093-0.939; VPE: P = 0.016-1.000). DATA CONCLUSION: The rCBVmax calculated in HAT, LAT, and IPE habitats have been validated as clinically relevant prognostic biomarkers for glioblastoma patients in the pretreatment stage. This study demonstrates the robustness of the hemodynamic tissue signature (HTS) habitats to assess the GBM vascular heterogeneity and their association with patient prognosis independently of intercenter variability. LEVEL OF EVIDENCE: 3 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2020;51:1478-1486.

ONCOhabitats: A system for glioblastoma heterogeneity assessment through MRI.

BACKGROUND: Neuroimaging analysis is currently crucial for an early assessment of glioblastoma, to help improving treatment and tumor follow-up. To this end, multiple functional and morphological MRI sequences are usually employed, requiring the development of automated tools capable to extract the relevant information from these sources. In this work we present ONCOhabitats (https://www.oncohabitats.upv.es): an online open access system for glioblastoma analysis based on MRI data. METHODS: ONCOhabitats provides two main services for untreated glioblastomas: (1) malignant tissue segmentation, and (2) vascular heterogeneity assessment of the tumor. The segmentation service implements a deep patch-wise 3D Convolutional Neural Network with residual connections. The vascular heterogeneity assessment service implements the Hemodynamic Tissue Signature (HTS) method patented in P201431289, which aims to identify habitats within the tumor with early prognostic capabilities. RESULTS: The segmentation service was validated against the BRATS 2017 reference dataset, showing comparable results with current state-of-the-art methods (whole tumor Dice segmentation: 0.89). The vascular heterogeneity assessment service was validated in a retrospective cohort of 50 patients, in a study focused on predicting patient overall survival based on the HTS habitats. Cox proportional hazard regression analysis and Kaplan-Meier survival study showed significant positive correlations (p-value <.05) between the HTS habitats and patient overall survival. ONCOhabitats system also generates radiological reports for each service, including volumetries and perfusion measurements of the different regions of the lesion. CONCLUSION: ONCOhabitats system provides open-access services for glioblastoma heterogeneity assessment, implementing consolidated state-of-the-art techniques for medical image analysis. Additionally, we also give access to the scientific community to our computational resources, offering a computational capacity of about 300 cases per day.

Multi-parametric MR Imaging Biomarkers Associated to Clinical Outcomes in Gliomas: A Systematic Review.

PURPOSE: To systematically review evidence regarding the association of multiparametric biomarkers with clinical outcomes and their capacity to explain relevant subcompartments of gliomas. MATERIALS AND METHODS: Scopus database was searched for original journal papers from January 1st, 2007 to February 20th, 2017 according to PRISMA. Four hundred forty-nine abstracts of papers were reviewed and scored independently by two out of six authors. Based on those papers we analyzed associations between biomarkers, subcompartments within the tumor lesion, and clinical outcomes. From all the articles analyzed, the twenty-seven papers with the highest scores were highlighted to represent the evidence about MR imaging biomarkers associated with clinical outcomes. Similarly, eighteen studies defining subcompartments within the tumor region were also highlighted to represent the evidence of MR imaging biomarkers. Their reports were critically appraised according to the QUADAS-2 criteria. RESULTS: It has been demonstrated that multi-parametric biomarkers are prepared for surrogating diagnosis, grading, segmentation, overall survival, progression-free survival, recurrence, molecular profiling and response to treatment in gliomas. Quantifications and radiomics features obtained from morphological exams (T1, T2, FLAIR, T1c), PWI (including DSC and DCE), diffusion (DWI, DTI) and chemical shift imaging (CSI) are the preferred MR biomarkers associated to clinical outcomes. Subcompartments relative to the peritumoral region, invasion, infiltration, proliferation, mass effect and pseudo flush, relapse compartments, gross tumor volumes, and highrisk regions have been defined to characterize the heterogeneity. For the majority of pairwise cooccurrences, we found no evidence to assert that observed co-occurrences were significantly different from their expected co-occurrences (Binomial test with False Discovery Rate correction, α=0.05). The co-occurrence among terms in the studied papers was found to be driven by their individual prevalence and trends in the literature. CONCLUSION: Combinations of MR imaging biomarkers from morphological, PWI, DWI and CSI exams have demonstrated their capability to predict clinical outcomes in different management moments of gliomas. Whereas morphologic-derived compartments have been mostly studied during the last ten years, new multi-parametric MRI approaches have also been proposed to discover specific subcompartments of the tumors. MR biomarkers from those subcompartments show the local behavior within the heterogeneous tumor and may quantify the prognosis and response to treatment of gliomas.