Improving the estimation of prognosis for glioblastoma patients by MR based hemodynamic tissue signatures

The prognostic relevance of a gene expression signature in MRI-defined highly vascularized glioblastoma

Background

The vascular heterogeneity of glioblastomas (GB) remains an important area of research, since tumor progression and patient prognosis are closely tied to this feature. With this study, we aim to identify gene expression profiles associated with MRI-defined tumor vascularity and to investigate its relationship with patient prognosis.

Methods

The study employed MRI parameters calculated with DSC Perfusion Quantification of ONCOhabitats glioma analysis software and RNA-seq data from the TCGA-GBM project dataset. In our study, we had a total of 147 RNA-seq samples, which 15 of them also had MRI parameter information. We analyzed the gene expression profiles associated with MRI-defined tumor vascularity using differential gene expression analysis and performed Log-rank tests to assess the correlation between the identified genes and patient prognosis.

Results

The findings of our research reveal a set of 21 overexpressed genes associated with the high vascularity pattern. Notably, several of these overexpressed genes have been previously implicated in worse prognosis based on existing literature. Our log-rank test further validates that the collective upregulation of these genes is indeed correlated with an unfavorable prognosis. This set of genes includes a variety of molecules, such as cytokines, receptors, ligands, and other molecules with diverse functions.

Conclusions

Our findings suggest that the set of 21 overexpressed genes in the High Vascularity group could potentially serve as prognostic markers for GB patients. These results highlight the importance of further investigating the relationship between the molecules such as cytokines or receptors underlying the vascularity in GB and its observation through MRI and developing targeted therapies for this aggressive disease.

Tumour Size and Overall Survival in a Cohort of Patients with Unifocal Glioblastoma: A Uni- and Multivariable Prognostic Modelling and Resampling Study

Published models inconsistently associate glioblastoma size with overall survival (OS). This study aimed to investigate the prognostic effect of tumour size in a large cohort of patients diagnosed with GBM and interrogate how sample size and non-linear transformations may impact on the likelihood of finding a prognostic effect. In total, 279 patients with a IDH-wildtype unifocal WHO grade 4 GBM between 2014 and 2020 from a retrospective cohort were included. Uni-/multivariable association between core volume, whole volume (CV and WV), and diameter with OS was assessed with (1) Cox proportional hazard models +/- log transformation and (2) resampling with 1,000,000 repetitions and varying sample size to identify the percentage of models, which showed a significant effect of tumour size. Models adjusted for operation type and a diameter model adjusted for all clinical variables remained significant (p = 0.03). Multivariable resampling increased the significant effects (p < 0.05) of all size variables as sample size increased. Log transformation also had a large effect on the chances of a prognostic effect of WV. For models adjusted for operation type, 19.5% of WV vs. 26.3% log-WV (n = 50) and 69.9% WV and 89.9% log-WV (n = 279) were significant. In this large well-curated cohort, multivariable modelling and resampling suggest tumour volume is prognostic at larger sample sizes and with log transformation for WV.

Spectral Analysis Based on Hemodynamic Habitat Imaging Predicts Isocitrate Dehydrogenase Status and Prognosis in High-Grade Glioma

BACKGROUND

The intratumoral heterogeneity of high-grade gliomas (HGGs) is associated with isocitrate dehydrogenase (IDH) status and prognosis, which can be established by quantitative radioanalysis of spatial tumor habitats. Therefore, we designed a framework for tackling tumors based on spatial metabolism using the hemodynamic tissue signature (HTS), focusing on metabolic changes in tumor habitat to predict IDH status and assess prognosis in patients with HGG.

METHODS

Preoperative data for 121 patients with HGG with subsequent histologic confirmation of HGG were prospectively collected (January 2016 to December 2020). The HTS was mapped from the image data, chemical shift imaging voxels were selected from the HTS habitat as the region of interest, and the metabolic ratio of the HTS was calculated using weighted least square method fitting. The metabolic rate of the tumor enhancement area was used as a control to analyze the efficacy of each HTS metabolic rate in predicting the IDH status and prognosis of HGG.

RESULTS

Total choline (Cho)/total creatine and Cho/N-acetyl-aspartate showed significant differences between IDH-wildtype and IDH-mutant in high- and low-angiogenic enhanced tumor sites (P < 0.05); Cho/total creatine was an independent risk factor for prognosis of HGG patients in high-angiogenic enhanced tumor habitats, with significant differences in survival time between groups (P < 0.05). The metabolic ratio in the tumor enhanced area could not predict IDH status or evaluate prognosis.

CONCLUSIONS

Spectral analysis based on hemodynamic habitat imaging can clearly distinguish IDH mutations and the prognosis assessment is more accurate, rendering it superior to traditional spectral analysis in tumor enhancement areas.

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.

Hemodynamic Imaging in Cerebral Diffuse Glioma-Part A: Concept, Differential Diagnosis and Tumor Grading

Diffuse gliomas are the most common primary malignant intracranial neoplasms. Aside from the challenges pertaining to their treatment-glioblastomas, in particular, have a dismal prognosis and are currently incurable-their pre-operative assessment using standard neuroimaging has several drawbacks, including broad differentials diagnosis, imprecise characterization of tumor subtype and definition of its infiltration in the surrounding brain parenchyma for accurate resection planning. As the pathophysiological alterations of tumor tissue are tightly linked to an aberrant vascularization, advanced hemodynamic imaging, in addition to other innovative approaches, has attracted considerable interest as a means to improve diffuse glioma characterization. In the present part A of our two-review series, the fundamental concepts, techniques and parameters of hemodynamic imaging are discussed in conjunction with their potential role in the differential diagnosis and grading of diffuse gliomas. In particular, recent evidence on dynamic susceptibility contrast, dynamic contrast-enhanced and arterial spin labeling magnetic resonance imaging are reviewed together with perfusion-computed tomography. While these techniques have provided encouraging results in terms of their sensitivity and specificity, the limitations deriving from a lack of standardized acquisition and processing have prevented their widespread clinical adoption, with current efforts aimed at overcoming the existing barriers.

Hemodynamic Imaging in Cerebral Diffuse Glioma—Part B: Molecular Correlates, Treatment Effect Monitoring, Prognosis, and Future Directions

Gliomas, and glioblastoma in particular, exhibit an extensive intra- and inter-tumoral molecular heterogeneity which represents complex biological features correlating to the efficacy of treatment response and survival. From a neuroimaging point of view, these specific molecular and histopathological features may be used to yield imaging biomarkers as surrogates for distinct tumor genotypes and phenotypes. The development of comprehensive glioma imaging markers has potential for improved glioma characterization that would assist in the clinical work-up of preoperative treatment planning and treatment effect monitoring. In particular, the differentiation of tumor recurrence or true progression from pseudoprogression, pseudoresponse, and radiation-induced necrosis can still not reliably be made through standard neuroimaging only. Given the abundant vascular and hemodynamic alterations present in diffuse glioma, advanced hemodynamic imaging approaches constitute an attractive area of clinical imaging development. In this context, the inclusion of objective measurable glioma imaging features may have the potential to enhance the individualized care of diffuse glioma patients, better informing of standard-of-care treatment efficacy and of novel therapies, such as the immunotherapies that are currently increasingly investigated. In Part B of this two-review series, we assess the available evidence pertaining to hemodynamic imaging for molecular feature prediction, in particular focusing on isocitrate dehydrogenase (IDH) mutation status, MGMT promoter methylation, 1p19q codeletion, and EGFR alterations. The results for the differentiation of tumor progression/recurrence from treatment effects have also been the focus of active research and are presented together with the prognostic correlations identified by advanced hemodynamic imaging studies. Finally, the state-of-the-art concepts and advancements of hemodynamic imaging modalities are reviewed together with the advantages derived from the implementation of radiomics and machine learning analyses pipelines.

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.

Survival prediction of glioblastoma patients-are we there yet? A systematic review of prognostic modeling for glioblastoma and its clinical potential

Glioblastoma is associated with a poor prognosis. Even though survival statistics are well-described at the population level, it remains challenging to predict the prognosis of an individual patient despite the increasing number of prognostic models. The aim of this study is to systematically review the literature on prognostic modeling in glioblastoma patients. A systematic literature search was performed to identify all relevant studies that developed a prognostic model for predicting overall survival in glioblastoma patients following the PRISMA guidelines. Participants, type of input, algorithm type, validation, and testing procedures were reviewed per prognostic model. Among 595 citations, 27 studies were included for qualitative review. The included studies developed and evaluated a total of 59 models, of which only seven were externally validated in a different patient cohort. The predictive performance among these studies varied widely according to the AUC (0.58-0.98), accuracy (0.69-0.98), and C-index (0.66-0.70). Three studies deployed their model as an online prediction tool, all of which were based on a statistical algorithm. The increasing performance of survival prediction models will aid personalized clinical decision-making in glioblastoma patients. The scientific realm is gravitating towards the use of machine learning models developed on high-dimensional data, often with promising results. However, none of these models has been implemented into clinical care. To facilitate the clinical implementation of high-performing survival prediction models, future efforts should focus on harmonizing data acquisition methods, improving model interpretability, and externally validating these models in multicentered, prospective fashion.

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.

Vascular habitat analysis based on dynamic susceptibility contrast perfusion MRI predicts IDH mutation status and prognosis in high-grade gliomas

OBJECTIVE

The current study aimed to evaluate the clinical practice for hemodynamic tissue signature (HTS) method in IDH genotype prediction in three groups derived from high-grade gliomas.

METHODS

Preoperative MRI examinations of 44 patients with known grade and IDH genotype were assigned into three study groups: glioblastoma multiforme, grade III, and high-grade gliomas. Perfusion parameters were analyzed and were used to automatically draw the four reproducible habitats (high-angiogenic enhancing tumor habitats, low-angiogenic enhancing tumor habitats, infiltrated peripheral edema habitats, vasogenic peripheral edema habitats) related to vascular heterogeneity. These four habitats were then compared between inter-patient with IDH mutation and their wild-type counterparts at these three groups, respectively. The discriminating potential for HTS in assessing IDH mutation status prediction was assessed by ROC curves.

RESULTS

Compared with IDH wild type, IDH mutation had significantly decreased relative cerebral blood volume (rCBV) at the high-angiogenic enhancing tumor habitats and low-angiogenic enhancing tumor habitats. ROC analysis revealed that the rCBVs in habitats had great ability to discriminate IDH mutation from their wild type in all groups. In addition, the Kaplan-Meier survival analysis yielded significant differences for the survival times observed from the populations dichotomized by low (< 4.31) and high (> 4.31) rCBV in the low-angiogenic enhancing tumor habitat.

CONCLUSIONS

The HTS method has been proven to have high prediction capabilities for IDH mutation status in high-grade glioma patients, providing a set of quantifiable habitats associated with tumor vascular heterogeneity.

KEY POINTS

• The HTS method has a high accuracy for molecular stratification prediction for all subsets of HGG. • The HTS method can give IDH mutation-related hemodynamic information of tumor-infiltrated and vasogenic edema. • IDH-relevant rCBV difference in habitats will be a great prognosis factor in HGG.