Prediction of survival and progression in glioblastoma patients using temporal perfusion changes during radiochemotherapy

Cluster-based prognostication in glioblastoma: Unveiling heterogeneity based on diffusion and perfusion similarities

BACKGROUND

While the association between diffusion and perfusion magnetic resonance imaging (MRI) and survival in glioblastoma is established, prognostic models for patients are lacking. This study employed clustering of functional imaging to identify distinct functional phenotypes in untreated glioblastomas, assessing their prognostic significance for overall survival.

METHODS

A total of 289 patients with glioblastoma who underwent preoperative multimodal MR imaging were included. Mean values of apparent diffusion coefficient normalized relative cerebral blood volume and relative cerebral blood flow were calculated for different tumor compartments and the entire tumor. Distinct imaging patterns were identified using partition around medoids (PAM) clustering on the training dataset, and their ability to predict overall survival was assessed. Additionally, tree-based machine-learning models were trained to ascertain the significance of features pertaining to cluster membership.

RESULTS

Using the training dataset (231/289) we identified 2 stable imaging phenotypes through PAM clustering with significantly different overall survival (OS). Validation in an independent test set revealed a high-risk group with a median OS of 10.2 months and a low-risk group with a median OS of 26.6 months (P = 0.012). Patients in the low-risk cluster had high diffusion and low perfusion values throughout, while the high-risk cluster displayed the reverse pattern. Including cluster membership in all multivariate Cox regression analyses improved performance (P ≤ 0.004 each).

CONCLUSIONS

Our research demonstrates that data-driven clustering can identify clinically relevant, distinct imaging phenotypes, highlighting the potential role of diffusion, and perfusion MRI in predicting survival rates of glioblastoma patients.

Diffusion- and Perfusion-Weighted MRI Radiomics for Survival Prediction in Patients with Lower-Grade Gliomas

PURPOSE

Lower-grade gliomas of histologic grades 2 and 3 follow heterogenous clinical outcomes, which necessitates risk stratification. This study aimed to evaluate whether diffusion-weighted and perfusion-weighted MRI radiomics allow overall survival (OS) prediction in patients with lower-grade gliomas and investigate its prognostic value.

MATERIALS AND METHODS

In this retrospective study, radiomic features were extracted from apparent diffusion coefficient, relative cerebral blood volume map, and Ktrans map in patients with pathologically confirmed lower-grade gliomas (January 2012-February 2019). The radiomics risk score (RRS) calculated from selected features constituted a radiomics model. Multivariable Cox regression analysis, including clinical features and RRS, was performed. The models' integrated area under the receiver operating characteristic curves (iAUCs) were compared. The radiomics model combined with clinical features was presented as a nomogram.

RESULTS

The study included 129 patients (median age, 44 years; interquartile range, 37-57 years; 63 female): 90 patients for training set and 39 patients for test set. The RRS was an independent risk factor for OS with a hazard ratio of 6.01. The combined clinical and radiomics model achieved superior performance for OS prediction compared to the clinical model in both training (iAUC, 0.82 vs. 0.72, p=0.002) and test sets (0.88 vs. 0.76, p=0.04). The radiomics nomogram combined with clinical features exhibited good agreement between the actual and predicted OS with C-index of 0.83 and 0.87 in the training and test sets, respectively.

CONCLUSION

Adding diffusion- and perfusion-weighted MRI radiomics to clinical features improved survival prediction in lower-grade glioma.

Imaging Role in Diagnosis, Prognosis, and Treatment Response Prediction Associated with High-grade Glioma

Background

Glioma is one of the most drug and radiation-resistant tumors. Gliomas suffer from inter- and intratumor heterogeneity which makes the outcome of similar treatment protocols vary from patient to patient. This article is aimed to overview the potential imaging markers for individual diagnosis, prognosis, and treatment response prediction in malignant glioma. Furthermore, the correlation between imaging findings and biological and clinical information of glioma patients is reviewed.

Materials and Methods

The search strategy in this study is to select related studies from scientific websites such as PubMed, Scopus, Google Scholar, and Web of Science published until 2022. It comprised a combination of keywords such as Biomarkers, Diagnosis, Prognosis, Imaging techniques, and malignant glioma, according to Medical Subject Headings.

Results

Some imaging parameters that are effective in glioma management include: ADC, FA, Ktrans, regional cerebral blood volume (rCBV), cerebral blood flow (CBF), ve, Cho/NAA and lactate/lipid ratios, intratumoral uptake of 18F-FET (for diagnostic application), RD, ADC, ve, vp, Ktrans, CBFT1, rCBV, tumor blood flow, Cho/NAA, lactate/lipid, MI/Cho, uptakes of 18F-FET, 11C-MET, and 18F-FLT (for prognostic and predictive application). Cerebral blood volume and Ktrans are related to molecular markers such as vascular endothelial growth factor (VEGF). Preoperative ADCmin value of GBM tumors is associated with O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation status. 2-hydroxyglutarate metabolite and dynamic 18F-FDOPA positron emission tomography uptake are related to isocitrate dehydrogenase (IDH) mutations.

Conclusion

Parameters including ADC, RD, FA, rCBV, Ktrans, vp, and uptake of 18F-FET are useful for diagnosis, prognosis, and treatment response prediction in glioma. A significant correlation between molecular markers such as VEGF, MGMT, and IDH mutations with some diffusion and perfusion imaging parameters has been identified.

Implementation and evaluation of a dynamic contrast-enhanced MR perfusion protocol for glioblastoma using a 0.35 T MRI-Linac system

PURPOSE

MRI-linear accelerator (MRI-Linac) systems allow for daily tracking of MRI changes during radiotherapy (RT). Since one common MRI-Linac operates at 0.35 T, there are efforts towards developing protocols at that field strength. In this study we demonstrate the implementation of a post-contrast 3DT1-weighted (3D-T1w) and dynamic contrast-enhancement (DCE) protocol to assess glioblastoma response to RT using a 0.35 T MRI-Linac.

METHODS AND MATERIALS

The protocol implemented was used to acquire 3D-T1w and DCE data from a flow phantom and two patients with glioblastoma (a responder and a non-responder) who underwent RT on a 0.35 T MRI-Linac. The detection of post-contrast-enhanced volumes was evaluated by comparing the 3DT1w images from the 0.35 T MRI-Linac to images obtained using a 3 T scanner. The DCE data were tested temporally and spatially using data from a flow phantom and patients. Ktrans maps were derived from DCE at three time points (a week before treatment-Pre RT, four weeks through treatment-Mid RT, and three weeks after treatment-Post RT) and were validated with patients' treatment outcomes.

RESULTS

The 3D-T1w contrast-enhancement volumes were visually and volumetrically similar between 0.35 T MRI-Linac and 3 T. DCE images showed temporal stability, and associated Ktrans maps were consistent with patient response to treatment. On average, Ktrans values showed a 54 % decrease and 8.6 % increase for a responder and non-responder respectively when Pre RT and Mid RT images were compared.

CONCLUSION

Our findings support the feasibility of obtaining post-contrast 3D-T1w and DCE data from patients with glioblastoma using a 0.35 T MRI-Linac system.

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

INTRODUCTION

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

METHODS

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

RESULTS

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

CONCLUSION

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

Exploratory Analysis of Serial 18F-fluciclovine PET-CT and Multiparametric MRI during Chemoradiation for Glioblastoma

Anti-1-amino-3-18fluorine-fluorocyclobutane-1-carboxylic acid (18F-fluciclovine) positron emission tomography (PET) shows preferential glioma uptake but there is little data on how uptake correlates with post-contrast T1-weighted (Gd-T1) and dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) activity during adjuvant treatment. This pilot study aimed to compare 18F-fluciclovine PET, DCE-MRI and Gd-T1 in patients undergoing chemoradiotherapy for glioblastoma (GBM), and in a parallel pre-clinical GBM model, to investigate correlation between 18F-fluciclovine uptake, MRI findings, and tumour biology. 18F-fluciclovine-PET-computed tomography (PET-CT) and MRI including DCE-MRI were acquired before, during and after adjuvant chemoradiotherapy (60 Gy in 30 fractions with temozolomide) in GBM patients. MRI volumes were manually contoured; PET volumes were defined using semi-automatic thresholding. The similarity of the PET and DCE-MRI volumes outside the Gd-T1 volume boundary was measured using the Dice similarity coefficient (DSC). CT-2A tumour-bearing mice underwent MRI and 18F-fluciclovine PET-CT. Post-mortem mice brains underwent immunohistochemistry staining for ASCT2 (amino acid transporter), nestin (stemness) and Ki-67 (proliferation) to assess for biologically active tumour. 6 patients were recruited (GBM 1-6) and grouped according to overall survival (OS)-short survival (GBM-SS, median OS 249 days) and long survival (GBM-LS, median 903 days). For GBM-SS, PET tumour volumes were greater than DCE-MRI, in turn greater than Gd-T1. For GBM-LS, Gd-T1 and DCE-MRI were greater than PET. Tumour-specific 18F-fluciclovine uptake on pre-clinical PET-CT corresponded to immunostaining for Ki-67, nestin and ASCT2. Results suggest volumes of 18F-fluciclovine-PET activity beyond that depicted by DCE-MRI and Gd-T1 are associated with poorer prognosis in patients undergoing chemoradiotherapy for GBM. The pre-clinical model confirmed 18F-fluciclovine uptake reflected biologically active tumour.

Paradoxical Association Between Relative Cerebral Blood Volume Dynamics Following Chemoradiation and Increased Progression-Free Survival in Newly Diagnosed IDH Wild-Type MGMT Promoter Methylated Glioblastoma With Measurable Disease

Background and Purpose

While relative cerebral blood volume (rCBV) may be diagnostic and prognostic for survival in glioblastoma (GBM), changes in rCBV during chemoradiation in the subset of newly diagnosed GBM with subtotal resection and the impact of MGMT promoter methylation status on survival have not been explored. This study aimed to investigate the association between rCBV response, MGMT methylation status, and progression-free (PFS) and overall survival (OS) in newly diagnosed GBM with measurable enhancing lesions.

Methods

1,153 newly diagnosed IDH wild-type GBM patients were screened and 53 patients (4.6%) had measurable post-surgical tumor (>1mL). rCBV was measured before and after patients underwent chemoradiation. Patients with a decrease in rCBV >10% were considered rCBV Responders, while patients with an increase or a decrease in rCBV <10% were considered rCBV Non-Responders. The association between change in enhancing tumor volume, change in rCBV, MGMT promotor methylation status, and PFS or OS were explored.

Results

A decrease in tumor volume following chemoradiation trended towards longer OS (p=0.12; median OS=26.8 vs. 16.3 months). Paradoxically, rCBV Non-Responders had a significantly improved PFS compared to Responders (p=0.047; median PFS=9.6 vs. 7.2 months). MGMT methylated rCBV Non-Responders exhibited a significantly longer PFS compared to MGMT unmethylated rCBV Non-Responders (p<0.001; median PFS=0.5 vs. 7.1 months), and MGMT methylated rCBV Non-Responders trended towards longer PFS compared to methylated rCBV Responders (p=0.089; median PFS=20.5 vs. 13.8 months).

Conclusions

This preliminary report demonstrates that in newly diagnosed IDH wild-type GBM with measurable enhancing disease after surgery (5% of patients), an enigmatic non-response in rCBV was associated with longer PFS, particularly in MGMT methylated patients.

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.

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.

Neuro-Oncology Practice Clinical Debate: FDG PET to differentiate glioblastoma recurrence from treatment-related changes

The ability to accurately differentiate treatment-related changes (ie, pseudoprogression and radiation necrosis) from recurrent glioma remains a critical diagnostic problem in neuro-oncology. Because these entities are treated differently and have vastly different outcomes, accurate diagnosis is necessary to provide optimal patient care. In current practice, this diagnostic quandary commonly requires either serial imaging or histopathologic tissue confirmation. In this article, experts in the field debate the utility of 2-deoxy-2[18F]fluoro-d-glucose positron emission tomography (FDG PET) as an imaging tool to distinguish tumor recurrence from treatment-related changes in a patient with glioblastoma and progressive contrast enhancement on magnetic resonance (MR) following chemoradiotherapy.

Algorithmically Deduced FREM2 Molecular Pathway Is a Potent Grade and Survival Biomarker of Human Gliomas

Gliomas are the most common malignant brain tumors with high mortality rates. Recently we showed that the FREM2 gene has a role in glioblastoma progression. Here we reconstructed the FREM2 molecular pathway using the human interactome model. We assessed the biomarker capacity of FREM2 expression and its pathway as the overall survival (OS) and progression-free survival (PFS) biomarkers. To this end, we used three literature and one experimental RNA sequencing datasets collectively covering 566 glioblastomas (GBM) and 1097 low-grade gliomas (LGG). The activation level of deduced FREM2 pathway showed strong biomarker characteristics and significantly outperformed the FREM2 expression level itself. For all relevant datasets, it could robustly discriminate GBM and LGG (p < 1.63 × 10-13, AUC > 0.74). High FREM2 pathway activation level was associated with poor OS in LGG (p < 0.001), and low PFS in LGG (p < 0.001) and GBM (p < 0.05). FREM2 pathway activation level was poor prognosis biomarker for OS (p < 0.05) and PFS (p < 0.05) in LGG with IDH mutation, for PFS in LGG with wild type IDH (p < 0.001) and mutant IDH with 1p/19q codeletion(p < 0.05), in GBM with unmethylated MGMT (p < 0.05), and in GBM with wild type IDH (p < 0.05). Thus, we conclude that the activation level of the FREM2 pathway is a potent new-generation diagnostic and prognostic biomarker for multiple molecular subtypes of GBM and LGG.

Multi-parametric MRI (mpMRI) for treatment response assessment of radiation therapy

Magnetic resonance imaging (MRI) plays an important role in the modern radiation therapy (RT) workflow. In comparison with computed tomography (CT) imaging, which is the dominant imaging modality in RT, MRI possesses excellent soft-tissue contrast for radiographic evaluation. Based on quantitative models, MRI can be used to assess tissue functional and physiological information. With the developments of scanner design, acquisition strategy, advanced data analysis, and modeling, multiparametric MRI (mpMRI), a combination of morphologic and functional imaging modalities, has been increasingly adopted for disease detection, localization, and characterization. Integration of mpMRI techniques into RT enriches the opportunities to individualize RT. In particular, RT response assessment using mpMRI allows for accurate characterization of both tissue anatomical and biochemical changes to support decision-making in monotherapy of radiation treatment and/or systematic cancer management. In recent years, accumulating evidence have, indeed, demonstrated the potentials of mpMRI in RT response assessment regarding patient stratification, trial benchmarking, early treatment intervention, and outcome modeling. Clinical application of mpMRI for treatment response assessment in routine radiation oncology workflow, however, is more complex than implementing an additional imaging protocol; mpMRI requires additional focus on optimal study design, practice standardization, and unified statistical reporting strategy to realize its full potential in the context of RT. In this article, the mpMRI theories, including image mechanism, protocol design, and data analysis, will be reviewed with a focus on the radiation oncology field. Representative works will be discussed to demonstrate how mpMRI can be used for RT response assessment. Additionally, issues and limits of current works, as well as challenges and potential future research directions, will also be discussed.

Operator dependency of arterial input function in dynamic contrast-enhanced MRI

Objective

To investigate the effect of inter-operator variability in arterial input function (AIF) definition on kinetic parameter estimates (KPEs) from dynamic contrast-enhanced (DCE) MRI in patients with high-grade gliomas.

Methods

The study included 118 DCE series from 23 patients. AIFs were measured by three domain experts (DEs), and a population AIF (pop-AIF) was constructed from the measured AIFs. The DE-AIFs, pop-AIF and AUC-normalized DE-AIFs were used for pharmacokinetic analysis with the extended Tofts model. AIF-dependence of KPEs was assessed by intraclass correlation coefficient (ICC) analysis, and the impact on relative longitudinal change in K^trans was assessed by Fleiss’ kappa (κ).

Results

There was a moderate to substantial agreement (ICC 0.51–0.76) between KPEs when using DE-AIFs, while AUC-normalized AIFs yielded ICC 0.77–0.95 for K^trans, k_ep and v_e and ICC 0.70 for v_p. Inclusion of the pop-AIF did not reduce agreement. Agreement in relative longitudinal change in K^trans was moderate (κ = 0.591) using DE-AIFs, while AUC-normalized AIFs gave substantial (κ = 0.809) agreement.

Discussion

AUC-normalized AIFs can reduce the variation in kinetic parameter results originating from operator input. The pop-AIF presented in this work may be applied in absence of a satisfactory measurement.

Radiomics-based neural network predicts recurrence patterns in glioblastoma using dynamic susceptibility contrast-enhanced MRI

Glioblastoma remains the most devastating brain tumor despite optimal treatment, because of the high rate of recurrence. Distant recurrence has distinct genomic alterations compared to local recurrence, which requires different treatment planning both in clinical practice and trials. To date, perfusion-weighted MRI has revealed that perfusional characteristics of tumor are associated with prognosis. However, not much research has focused on recurrence patterns in glioblastoma: namely, local and distant recurrence. Here, we propose two different neural network models to predict the recurrence patterns in glioblastoma that utilizes high-dimensional radiomic profiles based on perfusion MRI: area under the curve (AUC) (95% confidence interval), 0.969 (0.903-1.000) for local recurrence; 0.864 (0.726-0.976) for distant recurrence for each patient in the validation set. This creates an opportunity to provide personalized medicine in contrast to studies investigating only group differences. Moreover, interpretable deep learning identified that salient radiomic features for each recurrence pattern are related to perfusional intratumoral heterogeneity. We also demonstrated that the combined salient radiomic features, or "radiomic risk score", increased risk of recurrence/progression (hazard ratio, 1.61; p = 0.03) in multivariate Cox regression on progression-free survival.

MR-Guided Radiotherapy for Brain and Spine Tumors

MRI is the standard modality to assess anatomy and response to treatment in brain and spine tumors given its superb anatomic soft tissue contrast (e.g., T1 and T2) and numerous additional intrinsic contrast mechanisms that can be used to investigate physiology (e.g., diffusion, perfusion, spectroscopy). As such, hybrid MRI and radiotherapy (RT) devices hold unique promise for Magnetic Resonance guided Radiation Therapy (MRgRT). In the brain, MRgRT provides daily visualizations of evolving tumors that are not seen with cone beam CT guidance and cannot be fully characterized with occasional standalone MRI scans. Significant evolving anatomic changes during radiotherapy can be observed in patients with glioblastoma during the 6-week fractionated MRIgRT course. In this review, a case of rapidly changing symptomatic tumor is demonstrated for possible therapy adaptation. For stereotactic body RT of the spine, MRgRT acquires clear isotropic images of tumor in relation to spinal cord, cerebral spinal fluid, and nearby moving organs at risk such as bowel. This visualization allows for setup reassurance and the possibility of adaptive radiotherapy based on anatomy in difficult cases. A review of the literature for MR relaxometry, diffusion, perfusion, and spectroscopy during RT is also presented. These techniques are known to correlate with physiologic changes in the tumor such as cellularity, necrosis, and metabolism, and serve as early biomarkers of chemotherapy and RT response correlating with patient survival. While physiologic tumor investigations during RT have been limited by the feasibility and cost of obtaining frequent standalone MRIs, MRIgRT systems have enabled daily and widespread physiologic measurements. We demonstrate an example case of a poorly responding tumor on the 0.35 T MRIgRT system with relaxometry and diffusion measured several times per week. Future studies must elucidate which changes in MR-based physiologic metrics and at which timepoints best predict patient outcomes. This will lead to early treatment intensification for tumors identified to have the worst physiologic responses during RT in efforts to improve glioblastoma survival.