ONCOhabitats: A system for glioblastoma heterogeneity assessment through MRI

Preoperative vascular heterogeneity based on dynamic susceptibility contrast MRI in predicting spatial pattern of locally recurrent high-grade gliomas

OBJECTIVES

To investigate if spatial recurrence pattern is associated with patient prognosis, and whether MRI vascular habitats can predict spatial pattern.

METHODS

In this retrospective study, 69 patients with locally recurrent high-grade gliomas (HGGs) were included. The cohort was divided into intra-resection cavity recurrence (ICR) and extra-resection cavity recurrence (ECR) patterns, according to the distance between the location of the recurrent tumor and the resection cavity or surgical region. Four vascular habitats, high angiogenic tumor, low angiogenic tumor, infiltrated peripheral edema, and vasogenic peripheral edema, were segmented and vascular heterogeneity parameters were analyzed. The survival and diagnostic performance under different spatial recurrence patterns were analyzed by Kaplan-Meier and ROC. A nomogram model was constructed by regression analysis and validated by bootstrapping technique.

RESULTS

Progression-free survival (PFS) and overall survival (OS) were longer for ICR (n = 32) than those for ECR (n = 37) (median PFS: 8 vs. 5 months, median OS: 17 vs. 13 months, p < 0.05). MRI vascular habitat analyses showed ECR had higher median relative cerebral blood volume (rCBVmedian) at each habitat than ICR (all p < 0.01). The rCBVmedian at IPE had good diagnostic performance (AUC: 0.727, 95%CI: 0.607, 0.828). The AUC of the nomogram based on MRI vascular habitats and clinical factors was 0.834 (95%CI: 0.726, 0.913) and was confirmed as 0.833 (95%CI: 0.830, 0.836) by bootstrapping validation.

CONCLUSIONS

The spatial pattern of locally recurrent HGGs is associated with prognosis. MRI vascular heterogeneity parameter could be used as a non-invasive imaging marker to predict spatial recurrence pattern.

CLINICAL RELEVANCE STATEMENT

Vascular heterogeneity parameters based on MRI vascular habitat analyses can non-invasively predict the spatial patterns of locally recurrent high-grade gliomas, providing a new diagnostic basis for clinicians to develop the extent of surgical resection and postoperative radiotherapy planning.

KEY POINTS

• Intra-resection cavity pattern was associated with longer progression-free survival and overall survival in locally recurrent high-grade gliomas. • Higher vascular heterogeneities in extra-resection cavity recurrence than in intra-resection cavity recurrence and the vascular heterogeneity parameters had good diagnostic performance in discriminating spatial recurrence pattern. • A nomogram model based on MRI vascular habitats and clinical factors had good performance in predicting spatial recurrence pattern.

Androgen deficiency is associated with a better prognosis in glioblastoma

BACKGROUND

The androgen receptor (AR) has been demonstrated to play a role in the pathogenesis of glioblastoma; however, the implications of circulating testosterone levels in the biology of glioblastoma remain unknown.

AIM

This study aimed to analyze the association between circulating testosterone levels and the prognosis of patients with glioblastoma.

METHODS

Forty patients with primary glioblastoma were included in the study. The main prognostic endpoint was progression-free survival (PFS). Circulating testosterone levels were used to determine the state of androgen deficiency (AD). AR expression was analyzed by reverse-transcriptase polymerase chain reaction, Western blot, and immunofluorescence. Survival analysis was performed using the log-rank test and univariate and multivariate Cox regression analysis.

RESULTS

Most of the patients showed AR expression, and it was mainly located in the cytoplasm, as well as in the nucleus of tumor cells. Patients with AD presented a better PFS than those patients with normal levels (252.0 vs. 135.0 days; p = 0.041). Furthermore, normal androgenic status was an independent risk factor for progression in a multivariate regression model (hazard ratio = 6.346; p = 0.004).

CONCLUSION

Circulating testosterone levels are associated with the prognosis of glioblastoma because patients with AD show a better prognosis than those with normal androgenic status.

Diffuse glioma molecular profiling with arterial spin labeling and dynamic susceptibility contrast perfusion MRI: A comparative study

Background

Evaluation of molecular markers (IDH, pTERT, 1p/19q codeletion, and MGMT) in adult diffuse gliomas is crucial for accurate diagnosis and optimal treatment planning. Dynamic Susceptibility Contrast (DSC) and Arterial Spin Labeling (ASL) perfusion MRI techniques have both shown good performance in classifying molecular markers, however, their performance has not been compared side-by-side.

Methods

Pretreatment MRI data from 90 patients diagnosed with diffuse glioma (54 men/36 female, 53.1 ± 15.5 years, grades 2-4) were retrospectively analyzed. DSC-derived normalized cerebral blood flow/volume (nCBF/nCBV) and ASL-derived nCBF in tumor and perifocal edema were analyzed in patients with available IDH-mutation (n = 67), pTERT-mutation (n = 39), 1p/19q codeletion (n = 33), and MGMT promoter methylation (n = 31) status. Cross-validated uni- and multivariate logistic regression models assessed perfusion parameters' performance in molecular marker detection.

Results

ASL and DSC perfusion parameters in tumor and edema distinguished IDH-wildtype (wt) and pTERT-wt tumors from mutated ones. Univariate classification performance was comparable for ASL-nCBF and DSC-nCBV in IDH (maximum AUROCC 0.82 and 0.83, respectively) and pTERT (maximum AUROCC 0.70 and 0.81, respectively) status differentiation. The multivariate approach improved IDH (DSC-nCBV AUROCC 0.89) and pTERT (ASL-nCBF AUROCC 0.8 and DSC-nCBV AUROCC 0.86) classification. However, ASL and DSC parameters could not differentiate 1p/19q codeletion or MGMT promoter methylation status. Positive correlations were found between ASL-nCBF and DSC-nCBV/-nCBF in tumor and edema.

Conclusions

ASL is a viable gadolinium-free replacement for DSC for molecular characterization of adult diffuse gliomas.

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.

Preoperative neurocognitive function as an independent survival prognostic marker in primary glioblastoma

Background

Aim of the present study is to investigate whether preoperative neurocognitive status is prognostically associated with overall survival (OS) in newly diagnosed glioblastoma (GBM) patients.

Methods

Ninety patients with dominant-hemisphere IDH-wild-type GBM were assessed by Mini Mental Status Exam (MMSE), Trail Making Test (TMT) A and B parts, and Control Word Association Test (COWAT) phonemic and semantic subtests. Demographics, Karnofsky Performance Scale, tumor parameters, type of surgery, and adjuvant therapy data were available for patients.

Results

According to Cox proportional hazards model the neurocognitive variables of TMT B (P < .01), COWAT semantic subset (P < .05), and the MMSE (P < .01) were found significantly associated with survival prediction. From all other factors, only tumor volume and operation type (debulking vs biopsy) showed a statistical association (P < .05) with survival prediction. Kaplan Meier Long rank test showed statistical significance (P < .01) between unimpaired and impaired groups for TMT B, with median survival for the unimpaired group 26 months and 10 months for the impaired group, for COWAT semantic (P < .01) with median survival 23 months and 12 months, respectively and for MMSE (P < .01) with medial survival 19 and 12 months respectively.

Conclusions

Our study demonstrates that neurocognitive status at baseline-prior to treatment-is an independent prognostic factor for OS in wild-type GBM patients, adding another prognostic tool to assist physicians in selecting the best treatment plan.

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.

Visualising spatial heterogeneity in glioblastoma using imaging habitats

Glioblastoma is a high-grade aggressive neoplasm characterised by significant intra-tumoral spatial heterogeneity. Personalising therapy for this tumour requires non-invasive tools to visualise its heterogeneity to monitor treatment response on a regional level. To date, efforts to characterise glioblastoma's imaging features and heterogeneity have focussed on individual imaging biomarkers, or high-throughput radiomic approaches that consider a vast number of imaging variables across the tumour as a whole. Habitat imaging is a novel approach to cancer imaging that identifies tumour regions or 'habitats' based on shared imaging characteristics, usually defined using multiple imaging biomarkers. Habitat imaging reflects the evolution of imaging biomarkers and offers spatially preserved assessment of tumour physiological processes such perfusion and cellularity. This allows for regional assessment of treatment response to facilitate personalised therapy. In this review, we explore different methodologies to derive imaging habitats in glioblastoma, strategies to overcome its technical challenges, contrast experiences to other cancers, and describe potential clinical applications.

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.

Ellipsoid calculations versus manual tumor delineations for glioblastoma tumor volume evaluation

In glioblastoma, the response to treatment assessment is essentially based on the 2D tumor size evolution but remains disputable. Volumetric approaches were evaluated for a more accurate estimation of tumor size. This study included 57 patients and compared two volume measurement methods to determine the size of different glioblastoma regions of interest: the contrast-enhancing area, the necrotic area, the gross target volume and the volume of the edema area. The two methods, the ellipsoid formula (the calculated method) and the manual delineation (the measured method) showed a high correlation to determine glioblastoma volume and a high agreement to classify patients assessment response to treatment according to RANO criteria. This study revealed that calculated and measured methods could be used in clinical practice to estimate glioblastoma volume size and to evaluate tumor size evolution.

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.

High-Grade Glioma Treatment Response Monitoring Biomarkers: A Position Statement on the Evidence Supporting the Use of Advanced MRI Techniques in the Clinic, and the Latest Bench-to-Bedside Developments. Part 1: Perfusion and Diffusion Techniques

Objective

Summarize evidence for use of advanced MRI techniques as monitoring biomarkers in the clinic, and highlight the latest bench-to-bedside developments.

Methods

Experts in advanced MRI techniques applied to high-grade glioma treatment response assessment convened through a European framework. Current evidence regarding the potential for monitoring biomarkers in adult high-grade glioma is reviewed, and individual modalities of perfusion, permeability, and microstructure imaging are discussed (in Part 1 of two). In Part 2, we discuss modalities related to metabolism and/or chemical composition, appraise the clinic readiness of the individual modalities, and consider post-processing methodologies involving the combination of MRI approaches (multiparametric imaging) or machine learning (radiomics).

Results

High-grade glioma vasculature exhibits increased perfusion, blood volume, and permeability compared with normal brain tissue. Measures of cerebral blood volume derived from dynamic susceptibility contrast-enhanced MRI have consistently provided information about brain tumor growth and response to treatment; it is the most clinically validated advanced technique. Clinical studies have proven the potential of dynamic contrast-enhanced MRI for distinguishing post-treatment related effects from recurrence, but the optimal acquisition protocol, mode of analysis, parameter of highest diagnostic value, and optimal cut-off points remain to be established. Arterial spin labeling techniques do not require the injection of a contrast agent, and repeated measurements of cerebral blood flow can be performed. The absence of potential gadolinium deposition effects allows widespread use in pediatric patients and those with impaired renal function. More data are necessary to establish clinical validity as monitoring biomarkers. Diffusion-weighted imaging, apparent diffusion coefficient analysis, diffusion tensor or kurtosis imaging, intravoxel incoherent motion, and other microstructural modeling approaches also allow treatment response assessment; more robust data are required to validate these alone or when applied to post-processing methodologies.

Conclusion

Considerable progress has been made in the development of these monitoring biomarkers. Many techniques are in their infancy, whereas others have generated a larger body of evidence for clinical application.

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.

Identifying clinically applicable machine learning algorithms for glioma segmentation: recent advances and discoveries

Background

While there are innumerable machine learning (ML) research algorithms used for segmentation of gliomas, there is yet to be a US FDA cleared product. The aim of this study is to explore the systemic limitations of research algorithms that have prevented translation from concept to product by a review of the current research literature.

Methods

We performed a systematic literature review on 4 databases. Of 11 727 articles, 58 articles met the inclusion criteria and were used for data extraction and screening using TRIPOD.

Results

We found that while many articles were published on ML-based glioma segmentation and report high accuracy results, there were substantial limitations in the methods and results portions of the papers that result in difficulty reproducing the methods and translation into clinical practice.

Conclusions

In addition, we identified that more than a third of the articles used the same publicly available BRaTS and TCIA datasets and are responsible for the majority of patient data on which ML algorithms were trained, which leads to limited generalizability and potential for overfitting and bias.

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.

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.

DCE-MRI in Glioma, Infiltration Zone and Healthy Brain to Assess Angiogenesis: A Biopsy Study

Purpose

To explore the focal predictability of vascular growth factor expression and neovascularization using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) in glioma.

Methods

120 brain biopsies were taken in vital tumor, infiltration zone and normal brain tissue of 30 glioma patients: 17 IDH(isocitrate dehydrogenase)-wildtype glioblastoma (GBM), 1 IDH-wildtype astrocytoma °III (together prognostic group 1), 3 IDH-mutated GBM (group 2), 3 anaplastic astrocytomas IDH-mutated (group 3), 4 anaplastic oligodendrogliomas and 2 low-grade oligodendrogliomas (together prognostic group 4). A mixed linear model evaluated the predictabilities of microvessel density (MVD), vascular area ratio (VAR), mean vessel size (MVS), vascular endothelial growth factor and receptors (VEGF-A, VEGFR‑2) and vascular endothelial-protein tyrosine phosphatase (VE-PTP) expression from Tofts model kinetic and model-free curve parameters.

Results

All kinetic parameters were associated with VEGF‑A (all p < 0.001) expression. K^trans, k_ep and v_e were associated with VAR (p = 0.006, 0.004 and 0.01, respectively) and MVS (p = 0.0001, 0.02 and 0.003, respectively) but not MVD (p = 0.84, 0.74 and 0.73, respectively). Prognostic groups differed in K^trans (p = 0.007) and v_e (p = 0.004) values measured in the infiltration zone. Despite significant differences of VAR, MVS, VEGF‑A, VEGFR‑2, and VE-PTP in vital tumor tissue and the infiltration zone (p = 0.0001 for all), there was no significant difference between kinetic parameters measured in these zones.

Conclusion

The DCE-MRI kinetic parameters show correlations with microvascular parameters in vital tissue and also reveal blood-brain barrier abnormalities in the infiltration zones adequate to differentiate glioma prognostic groups.

Supplementary Information

The online version of this article (10.1007/s00062-021-01015-3) contains supplementary material, which is available to authorized users.

Noise dependency in vascular parameters from combined gradient-echo and spin-echo DSC MRI

Dynamic susceptibility contrast (DSC) imaging is a widely used technique for assessment of cerebral blood volume (CBV). With combined gradient-echo and spin-echo DSC techniques, measures of the underlying vessel size and vessel architecture can be obtained from the vessel size index (VSI) and vortex area, respectively. However, how noise, and specifically the contrast-to-noise ratio (CNR), affect the estimations of these parameters has largely been overlooked. In order to address this issue, we have performed simulations to generate DSC signals with varying levels of CNR, defined by the peak of relaxation rate curve divided by the standard deviation of the baseline. Moreover, DSC data from 59 brain cancer patients were acquired at two different 3 T-scanners (N = 29 and N = 30, respectively), where CNR and relative parameter maps were obtained. Our simulations showed that the measured parameters were affected by CNR in different ways, where low CNR led to overestimations of CBV and underestimations of VSI and vortex area. In addition, a higher noise-sensitivity was found in vortex area than in CBV and VSI. Results from clinical data were consistent with simulations, and indicated that CNR < 4 gives highly unreliable measurements. Moreover, we have shown that the distribution of values in the tumour regions could change considerably when voxels with CNR below a given cut off are excluded when generating the relative parameter maps. The widespread use of CBV and attractive potential of VSI and vortex area, makes the noise-sensitivity of these parameters found in our study relevant for further use and development of the DSC imaging technique. Our results suggest that the CNR has considerable impact on the measured parameters, with the potential to affect the clinical interpretation of DSC-MRI, and should therefore be taken into account in the clinical decision-making process.

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.

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.

Electronic supplementary material

The online version of this article (10.1007/s00330-020-07297-4) contains supplementary material, which is available to authorized users.

Multi-Parametric Deep Learning Model for Prediction of Overall Survival after Postoperative Concurrent Chemoradiotherapy in Glioblastoma Patients

This study aimed to investigate the performance of a deep learning-based survival-prediction model, which predicts the overall survival (OS) time of glioblastoma patients who have received surgery followed by concurrent chemoradiotherapy (CCRT). The medical records of glioblastoma patients who had received surgery and CCRT between January 2011 and December 2017 were retrospectively reviewed. Based on our inclusion criteria, 118 patients were selected and semi-randomly allocated to training and test datasets (3:1 ratio, respectively). A convolutional neural network–based deep learning model was trained with magnetic resonance imaging (MRI) data and clinical profiles to predict OS. The MRI was reconstructed by using four pulse sequences (22 slices) and nine images were selected based on the longest slice of glioblastoma by a physician for each pulse sequence. The clinical profiles consist of personal, genetic, and treatment factors. The concordance index (C-index) and integrated area under the curve (iAUC) of the time-dependent area-under-the-curve curves of each model were calculated to evaluate the performance of the survival-prediction models. The model that incorporated clinical and radiomic features showed a higher C-index (0.768 (95% confidence interval (CI): 0.759, 0.776)) and iAUC (0.790 (95% CI: 0.783, 0.797)) than the model using clinical features alone (C-index = 0.693 (95% CI: 0.685, 0.701); iAUC = 0.723 (95% CI: 0.716, 0.731)) and the model using radiomic features alone (C-index = 0.590 (95% CI: 0.579, 0.600); iAUC = 0.614 (95% CI: 0.607, 0.621)). These improvements to the C-indexes and iAUCs were validated using the 1000-times bootstrapping method; all were statistically significant (p < 0.001). This study suggests the synergistic benefits of using both clinical and radiomic parameters. Furthermore, it indicates the potential of multi-parametric deep learning models for the survival prediction of glioblastoma patients.