Next Article in Journal
Dissecting the Methylomes of EGFR-Amplified Glioblastoma Reveals Altered DNA Replication and Packaging, and Chromatin and Gene Silencing Pathways
Next Article in Special Issue
The Impact of ioMRI on Glioblastoma Resection and Clinical Outcomes in a State-of-the-Art Neuro-Oncological Setup
Previous Article in Journal
Clinical Decision Support Systems for Brain Tumour Diagnosis and Prognosis: A Systematic Review
Previous Article in Special Issue
Simulated Adaptive Radiotherapy for Shrinking Glioblastoma Resection Cavities on a Hybrid MRI–Linear Accelerator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spectroscopic MRI-Based Biomarkers Predict Survival for Newly Diagnosed Glioblastoma in a Clinical Trial

by
Anuradha G. Trivedi
1,2,†,
Karthik K. Ramesh
1,2,†,
Vicki Huang
1,2,
Eric A. Mellon
3,
Peter B. Barker
4,
Lawrence R. Kleinberg
5,
Brent D. Weinberg
6,7,
Hui-Kuo G. Shu
1,7,* and
Hyunsuk Shim
1,2,6,7,*
1
Department of Radiation Oncology, Emory University School of Medicine, Atlanta, GA 30322, USA
2
Department of Biomedical Engineering, Emory University and Georgia Institute of Technology, Atlanta, GA 30332, USA
3
Department of Radiation Oncology, Sylvester Comprehensive Cancer Center, Miller School of Medicine, University of Miami, Miami, FL 45056, USA
4
Department of Radiology and Radiological Science, Johns Hopkins University, Baltimore, MD 21218, USA
5
Department of Radiation Oncology, Johns Hopkins University, Baltimore, MD 21218, USA
6
Department of Radiology and Imaging Sciences, Emory University School of Medicine, Atlanta, GA 30322, USA
7
Winship Cancer Institute, Emory University School of Medicine, Atlanta, GA 30322, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2023, 15(13), 3524; https://doi.org/10.3390/cancers15133524
Submission received: 26 May 2023 / Revised: 22 June 2023 / Accepted: 5 July 2023 / Published: 7 July 2023
(This article belongs to the Special Issue Advanced Imaging in Brain Tumor Patient Management)

Abstract

:

Simple Summary

Due to the infiltrative nature of glioblastoma, standard MRI techniques, such as T1-weighted contrast-enhanced (T1w-CE) and T2-weighted fluid-attenuated inversion recovery (FLAIR), imperfectly delineate radiation-targeted tumor volume. With spectroscopic MRI, the ratio of choline elevation, a tumor biomarker and N-acetylaspartate reduction, a healthy neuronal biomarker, can better determine the extent of the tumor. The aim of our secondary analysis was to determine if there was a relationship between survival outcomes and biomarkers identified by spectroscopic MRI for a cohort of 28 glioblastoma patients who received high-dose radiation guided by spectroscopic MRI. We determined that the volume of post-surgical spectroscopically abnormal tissue was a biomarker of overall and progression-free survival, whereas the volume of residual contrast enhancement, determined by T1w-CE MRI, was not. Our results suggest that accurate delineation and treatment of an infiltrative tumor not identified by contrast is a critical component of glioblastoma management and patient survival.

Abstract

Despite aggressive treatment, glioblastoma has a poor prognosis due to its infiltrative nature. Spectroscopic MRI-measured brain metabolites, particularly the choline to N-acetylaspartate ratio (Cho/NAA), better characterizes the extent of tumor infiltration. In a previous pilot trial (NCT03137888), brain regions with Cho/NAA ≥ 2x normal were treated with high-dose radiation for newly diagnosed glioblastoma patients. This report is a secondary analysis of that trial where spectroscopic MRI-based biomarkers are evaluated for how they correlate with progression-free and overall survival (PFS/OS). Subgroups were created within the cohort based on pre-radiation treatment (pre-RT) median cutoff volumes of residual enhancement (2.1 cc) and metabolically abnormal volumes used for treatment (19.2 cc). We generated Kaplan–Meier PFS/OS curves and compared these curves via the log-rank test between subgroups. For the subgroups stratified by metabolic abnormality, statistically significant differences were observed for PFS (p = 0.019) and OS (p = 0.020). Stratification by residual enhancement did not lead to observable differences in the OS (p = 0.373) or PFS (p = 0.286) curves. This retrospective analysis shows that patients with lower post-surgical Cho/NAA volumes had significantly superior survival outcomes, while residual enhancement, which guides high-dose radiation in standard treatment, had little significance in PFS/OS. This suggests that the infiltrating, non-enhancing component of glioblastoma is an important factor in patient outcomes and should be treated accordingly.

1. Introduction

Glioblastoma is the most common malignant primary brain tumor in need of novel treatment strategies due to its poor prognosis [1]. Despite improvements in the standard of care for glioblastoma patients, it is challenging to identify the extent of tumor infiltration with T1-weighted contrast-enhanced (T1w-CE) MRI and T2-weighted fluid-attenuated inversion recovery (FLAIR) MRI, which are currently used to determine the targets of radiation therapy [2]. Glioblastomas are characterized by a compromised blood–brain barrier and leaky vasculature, which can be identified by a gadolinium-based contrast agent in T1w-CE scans. FLAIR is used to identify tumor, edema, inflammation, and radiation effects but is not specific to tumor [3,4]. The combination of imaging methods may not identify the full extent of tumor infiltration, resulting in undertreatment during radiation therapy. The standard of care radiation therapy target is the residual contrast enhancing contour and resection cavity treated with 60 Gy with a lower dose to the surrounding FLAIR abnormality [5]. Temozolomide is administered concurrently, during and after radiation. This treatment results in a median overall survival (OS) of 16 months and median progression-free survival (PFS) of 4–7 months in historical cohorts [1,6,7].
To better identify infiltrative tumors for treatment targeting, members of our group have been developing a whole-brain, high-resolution, 3D echo-planar magnetic resonance spectroscopic imaging (MRSI) sequence that we have termed “spectroscopic MRI” [8,9,10,11]. Tumor cell metabolism differs from healthy brain tissue metabolism in that there are increased choline and reduced N-acetylaspartate levels in proliferating tumor cells. The ratio of choline to N-acetylaspartate (Cho/NAA) has been shown to be a highly specific tumor biomarker by image–histology correlation studies and predated tumor recurrence patterns [12].
We recently reported results from a multisite clinical study (NCT03137888) where 30 newly diagnosed glioblastoma patients were treated with escalated dose radiation (75 Gy), guided by the Cho/NAA abnormality on spectroscopic MRI. This treatment paradigm resulted in a median OS of 23.0 months and PFS of 16.6 months, a significant improvement from patients treated with standard of care radiation therapy dosing [13].
In this report, we perform a retrospective analysis of the data from our previously mentioned clinical study. We aim to assess spectroscopic MRI-based biomarkers of patient OS and PFS. Because historical results show that extent of resection (EOR) and residual contrast-enhancing volume are factors that correlate with survival, we hypothesize that incorporating metabolic information from spectroscopic MRI will better predict survival outcomes. We tested this using the Kaplan–Meier estimator for predicting PFS and OS and analyzed the difference in survival distributions based on the residual contrast-enhancing and spectroscopic MRI-based volumes.

2. Materials and Methods

2.1. Description of Tumor Volume Determination/Target Generation

Patients enrolled in this study had pathologically confirmed, newly diagnosed World Health Organization (WHO) grade IV glioblastoma, which at the time of enrollment, included isocitrate dehydrogenase (IDH) mutant patients. In 2021, the WHO definition of grade IV glioblastoma was changed to exclude IDH mutant patients [14]. In accordance with the updated definition, we performed our analysis excluding the two IDH mutant patients (n = 28). Further trial details, including patient selection and sample selection, as well as initial outcome results for all 30 enrolled patients, have been previously reported [13]. The patients in this study provided informed consent for participation according to the institutional review board at each participating institution, and each institutional review board approved this study. After surgical resection and prior to beginning radiation therapy, a spectroscopic MRI was acquired in addition to a standard MRI (T1w-CE and FLAIR MRIs). Spectroscopic MRIs were acquired on Siemens 3T scanners at each institution with an echo planar spectroscopic imaging pulse sequence with GRAPPA parallelization. Either a 20- or 32-channel head and neck coil was used with TE = 50 ms, TR = 1551 ms, and FA = 71° [13]. Scans were processed and registered to anatomic T1 MRIs for each patient before radiation treatment planning [15,16,17,18,19].
Pre-resection scans were obtained 0–18 days before surgery, and post-resection scans were obtained approximately 3–4 weeks after surgery and 1 week prior to the radiation therapy start date. For treatment planning, the contrast-enhancing regions of the pre-resection T1w-CE, including necrotic centers, were contoured. The residual post-contrast enhancement (rENH) tumor volume was semi-automatically generated on the Brain Imaging Collaboration Suite using the previously reported method by contouring only the contrast-enhancing regions of the T1w-CE scan for each patient inside of the FLAIR envelope, excluding the resection cavity [20,21]. All rENH volumes were checked and, if necessary, edited by a board certified neuroradiologist. The Brain Imaging Collaboration Suite also automatically generated a spectroscopic MRI contour for Cho/NAA ≥ 2x, which was then manually edited by site MRS experts based on spectral quality. The gross tumor volume 3 (GTV3) was calculated for each patient by combining the rENH volume with the Cho/NAA ≥ 2x volume to determine the escalated dose target. The GTV3 contours were inspected and approved by two board-certified radiation oncologists (one of them was the treating physician). This workflow was published in Gurbani et al. [21].
Briefly, for each patient enrolled, radiation treatment contours were guided by standard MRIs and spectroscopic MRI data as follows: FLAIR abnormal regions of the brain received 50.1 Gy of radiation, the resection cavity received 60 Gy, and the GTV3 contour received a boosted 75 Gy of radiation [21]. After completing radiation therapy, imaging was acquired every 2–3 months as per standard practice. This follow-up imaging was used to determine PFS, and OS was determined through chart review and communication between patients, their families, their oncology team, and published obituaries [13]. In this report, OS and PFS were calculated as of 17 months after completing enrollment of the last patient.
In Figure 1, an example patient from the completed clinical trial is shown. This patient had a gross total resection (GTR) with no residual enhancing tumor. For this patient, standard of care would have used the resection cavity with a margin for high-dose radiation (60 Gy). This figure demonstrates how, even when using a 20 mm margin, some abnormal tissue in the GTV3 volume detected by spectroscopic MRI is left out. Furthermore, most clinical practices use a margin between 5 and 10 mm, which, in this case, would have left out a significant portion of tumor that spectroscopic MRI was able to detect. It is important to note that the Dice similarity coefficients in the figure show how much of the cavity and margins overlap with the GTV3 contour. In addition, using large, homogeneous margins can lead to unnecessary treatment of healthy tissue with high-dose radiation. We calculated the percentages of volumes that would be outside of the GTV3 contour. For a 5, 10, and 20 mm expansion beyond the cavity, 41.3%, 49.7%, and 74.9% of each volume would fall outside of the GTV3 volume, respectively. While there is some overlap between standard of care high-dose targets and our GTV3 target, the target guided by spectroscopic MRI illustrates metabolically active tumors undetected by T1w-CE.
In Figure 2, two subjects are shown; Case 2 is the same patient shown in Figure 1. The vast difference in rENH and GTV3 contour volumes for the two subjects is shown here, particularly in Case 2 where the volume of residual enhancement was 0 cc. Both contours are also rendered in 3D to emphasize the spatial differences between the rENH and GTV3 contours.

2.2. Surgical Resection Classification

EOR was calculated based on the percent change between the pre-operative contrast-enhancing tumor volume and the rENH volume. When EOR ≥ 95%, the surgery was classified as gross total resection (GTR), while an EOR < 95% was considered subtotal resection (STR), excluding biopsy. A secondary method of surgical resection classification was also used, in which resection was classified as GTR if the rENH volume was ≤1.5 cc or STR if >1.5 cc. Surgery classification by EOR and rENH volume was used to compare the 28 patients considered in this analysis to patients in large historical cohorts [22,23,24].

2.3. Statistical Analysis/Survival Analysis

The Kaplan–Meier estimator was used to calculate survival curves based on OS and PFS, with survival measured as the time between surgical resection for each patient and their recorded OS and PFS dates [25]. Patients were first stratified into two groups using median rENH and median GTV3 volumes. Each Kaplan–Meier curve was compared using a log-rank test to determine differences between patients with rENH volumes greater than and less than the median rENH volume, and patients with GTV3 volumes greater than and less than the median GTV3 volume with regard to OS and PFS. Linear regression was performed and the associated R2 value was used to determine any correlation between the rENH and GTV3 volumes, which would confound any comparisons. Statistical analysis was conducted using the Python lifelines analysis library [26]. The secondary analysis performed in this report was approved by the institutional review boards at each institution.

3. Results

Of the 30 patients included in the original study, 2 had IDH mutations (IDH1 R132H) and 9 were O6-methylguanine-DNA methyltransferase (MGMT) hypermethylated. Details including age, IDH and MGMT statuses, resection classification, and treatment volumes for each patient can be found in Supplemental Table S1; subjects 1 and 15 were excluded from the analysis performed in this report due to their IDH mutations. As reported previously, 16 patients had died with a median OS of 23.0 months [13]. Excluding the two IDH mutant patients, the median time to follow-up for censored patients was 22.4 months (5.7–30.1 months), and the median PFS was 16.2 months, as shown on a Kaplan–Meier curve in Figure 3. The median age for the 28 patients was 59.3 years with a range of 38.1–71.6 years (minimum–maximum) and there were 10 females and 18 males. Without the two IDH mutant patients, the median OS becomes 21.6 months and the median PFS becomes 16.2 months.
The volume of high-risk disease remaining should be a predictor of survival outcomes for glioblastoma patients. EOR or rENH are classic ways to identify the amount of high-risk disease and have been shown in multiple previous studies to be correlated with both OS and PFS. Twelve out of twenty-eight patients in our cohort had GTR (defined by either EOR ≥ 95% or rENH ≤ 1.5 cc), as shown in Table 1. Using this cutoff (GTR versus less than GTR), we did not find a significant difference in survival outcomes. However, pre-treatment spectroscopic MRI scans present the opportunity to apply additional criteria for identifying high-risk disease that may be superior to rENH alone.
On our pilot study, GTV3 volume was hypothesized to provide a better representation of residual high-risk disease than rENH alone. For our cohort, median rENH was 2.1 cc (range: 0.0–15.4 cc) and median GTV3 was 19.2 cc (range: 0.9–65.0 cc). The GTV3 volumes were always greater than rENH volumes (median of 8.2 times larger). In Figure 4, a scatterplot of rENH vs. GTV3 volumes for each patient indicates a low correlation between the volumes (R2 = 0.127, p = 0.063).
Figure 5 shows Kaplan–Meier curves for OS and PFS using the rENH and GTV3 volumes as predictive markers. Figure 5A,B shows Kaplan–Meier curves for OS. In Figure 5A, the 28 patients are stratified by those with an rENH volume below and above the median (2.1 cc). Similarly, in Figure 5B, the 28 patients are stratified by those with a GTV3 volume below and above the median (19.2 cc). The difference in medians for the groups stratified by rENH is 1.7 months (p = 0.373), while the difference in medians for the groups stratified by GTV3 is 10.3 months (p = 0.020), indicating that GTV3 volume is a significant biomarker of OS. For patients with rENH < 2.1 cc, the median OS was 23.0 months; for rENH > 2.1 cc, the median OS was 21.3 months. For patients with GTV3 < 19.2 cc, the median OS was 23.0 months; for GTV3 > 19.2 cc, the median OS was 17.1 months.
Figure 5C,D shows Kaplan–Meier curves for PFS using the same predictive markers as the OS graphs. The difference in medians for the groups stratified by rENH is 6.2 months shown in Figure 5C (p = 0.286), while the difference in medians for the groups stratified by GTV3 is 11.4 months shown in Figure 5D (p = 0.019), indicating that GTV3 volume is a better predictive marker for PFS. For patients with rENH < 2.1 cc, the median PFS was 19.0 months; for rENH > 2.1 cc, the median PFS was 12.8 months. For patients with GTV3 < 19.2 cc, median PFS was 24.0 months; for GTV3 > 19.2 cc, the median PFS was 12.6 months.

4. Discussion

In this paper, we have used data from a clinical trial using spectroscopic MRI to guide dose-escalation therapy to glioblastoma to assess the relationship between the rENH tumor and metabolically active tumor, as detected by spectroscopic MRI, and assess survival outcomes (PFS and OS) to determine if residual spectroscopically abnormal tumor (GTV3) was a biomarker for survival. For every subject, the rENH volume was smaller than residual non-enhancing tumor that was detected using spectroscopy, since GTV3 contains both residual enhancing tumor and non-enhancing metabolically active tumor. Additionally, there was low correlation between GTV3 and rENH volumes (R2 = 0.127). The lack of a linear relationship between the volumes is further supported by the insignificant p-value (p = 0.063).
The amount of residual spectroscopically abnormal tumor, as measured by the GTV3 volume, was a statistically significant biomarker for both OS (p = 0.020) and PFS (p = 0.019). We hypothesize that these results are attributed to the metabolically active tumor within the GTV3 volume that is not accounted for in the rENH volume. The large gap and lack of crossover between the GTV3-stratified Kaplan–Meier curves suggest that GTV3 is a better indicator of OS and PFS than rENH in this patient cohort. For this reason, GTV3 may be a better predictor of PFS and OS because it more accurately reflects how much non-infiltrating tumor is left post-resection. It also provides more precise targeting compared to adding a 1–2 cm margin around the resection cavity for the high-dose target, allowing higher radiation doses to be used with a decreased risk of damaging healthy tissue. Unlike prior studies, there was no relationship between rENH and survival or EOR and survival [27]. This result may be attributed to the relatively low amounts of residual contrast-enhancing tumor on most patients in this study, due to aggressive resection or higher radiation doses given to these regions.
For this study, the EOR classified using rENH volume (using a threshold of 1.5 cc) or EOR percentage (using a threshold of 95%) was similar to large historical cohorts, such as RTOG 0825 (637 patients), AVAGlio (921 patients), and EF-14 (695 patients) [22,23,24]. A 2005–2010 study of 128 patients from the Cleveland Clinic reported a similar distribution based on rENH and EOR [28]. Their study resulted in a median survival of 16 months for patients with rENH ≤ 1.5 cc and a median survival between 14 and 15 months for patients with EOR ≥ 95%. Spectroscopic MRI-guided dose-escalation resulted in a median survival of 23 months for patients with rENH ≤ 1.5 cc and a median survival of 27 months for patients with EOR ≥ 95%, which are more favorable outcomes. Similar studies have been previously reported, often with factors such as GTR versus STR, diffusion, or other volumes. Haj et al. performed a survival analysis between patients that received GTR versus STR [29]. They identified a 4-month difference in OS (GTR = 18.5 months, STR = 14.5 months) and a 1.6-month difference in PFS (GTR = 9.4 months, STR = 7.8 months). Awad et al. performed a survival analysis using only EOR and a combination of EOR and pre-resection contrast-enhancing volume [30]. Their EOR-only comparison failed to identify a relationship with survival outcomes; however, their EOR/pre-resection contrast-enhancing volume comparison was identified as a significant predictor of patient survival. Compared to the studies mentioned above, the number of patients in our trial is smaller and a larger study is warranted for fairer comparisons.
Limitations of this analysis include the relatively small cohort size, a treatment paradigm which was not a clinical standard of care, and data from only three clinical sites. More data is needed to determine if the amount of spectroscopically abnormal residual tumor after resection continues to have this predictive effect in larger patient groups, but this preliminary data is supportive of further study of using the extent of metabolically abnormal tissue to predict patient outcomes.

5. Conclusions

Due to the infiltrative nature of glioblastoma, T1w-CE and FLAIR MRIs may not optimally characterize the extent of tumor, resulting in undertreatment during radiation therapy. Our goal was to use data from a previously performed clinical trial to determine how spectroscopic MRI-based biomarkers predicted PFS and OS for patients in the trial. In this group, the volume of spectroscopically abnormal tissue (GTV3) was a statistically significant biomarker of both OS and PFS, and performed superior to the extent of resection and residual enhancing tumor. The results of this analysis support using spectroscopic MRI to guide high-dose radiation and suggest that the non-enhancing component of glioblastoma is an important biomarker of survival and should be treated accordingly.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15133524/s1. Table S1: Details about all 30 patients included in the original study.

Author Contributions

Conceptualization, E.A.M., L.R.K., B.D.W., H.-K.G.S. and H.S.; Formal analysis, A.G.T. and K.K.R.; Funding acquisition, E.A.M. and H.S.; Methodology, B.D.W., H.-K.G.S. and H.S.; Software, A.G.T. and K.K.R.; Supervision, B.D.W., H.-K.G.S. and H.S.; Validation, B.D.W., H.-K.G.S. and H.S.; Visualization, A.G.T. and K.K.R.; Writing—original draft, A.G.T. and K.K.R.; Writing—review and editing, A.G.T., K.K.R., V.H., E.A.M., P.B.B., L.R.K., B.D.W., H.-K.G.S. and H.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NIH U01 CA264039 (E.A.M. and H.S.) and NIH F31 CA247564 (K.K.R.).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of Emory University (IRB00094188 on 14 April 2017).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank Eduard Schreibmann for his assistance in verifying treatment margins.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stupp, R.; Mason, W.P.; Van Den Bent, M.J.; Weller, M.; Fisher, B.; Taphoorn, M.J.B.; Belanger, K.; Brandes, A.A.; Marosi, C.; Bogdahn, U.; et al. Radiotherapy plus Concomitant and Adjuvant Temozolomide for Glioblastoma. N. Engl. J. Med. 2005, 352, 987–996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kelly, P.J.; Daumas-Duport, C.; Kispert, D.B.; Kall, B.A.; Scheithauer, B.W.; Illig, J.J. Imaging-based stereotaxic serial biopsies in untreated intracranial glial neoplasms. J. Neurosurg. 1987, 66, 865–874. [Google Scholar] [CrossRef] [PubMed]
  3. Pope, W.B.; Young, J.R.; Ellingson, B.M. Advances in MRI assessment of gliomas and response to anti-VEGF therapy. Curr. Neurol. Neurosci. Rep. 2011, 11, 336–344. [Google Scholar] [CrossRef] [Green Version]
  4. Tsuchiya, K.; Mizutani, Y.; Hachiya, J. Preliminary evaluation of fluid-attenuated inversion-recovery MR in the diagnosis of intracranial tumors. AJNR Am. J. Neuroradiol. 1996, 17, 1081–1086. [Google Scholar] [PubMed]
  5. Wernicke, A.G.; Smith, A.W.; Taube, S.; Mehta, M.P. Glioblastoma: Radiation treatment margins, how small is large enough? Pract. Radiat. Oncol. 2016, 6, 298–305. [Google Scholar] [CrossRef] [PubMed]
  6. Stupp, R.; Hegi, M.E.; Mason, W.P.; van den Bent, M.J.; Taphoorn, M.J.; Janzer, R.C.; Ludwin, S.K.; Allgeier, A.; Fisher, B.; Belanger, K.; et al. Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-year analysis of the EORTC-NCIC trial. Lancet Oncol. 2009, 10, 459–466. [Google Scholar] [CrossRef]
  7. Stupp, R.; Taillibert, S.; Kanner, A.A.; Kesari, S.; Steinberg, D.M.; Toms, S.A.; Taylor, L.P.; Lieberman, F.; Silvani, A.; Fink, K.L.; et al. Maintenance Therapy with Tumor-Treating Fields Plus Temozolomide vs Temozolomide Alone for Glioblastoma: A Randomized Clinical Trial. JAMA 2015, 314, 2535–2543. [Google Scholar] [CrossRef]
  8. Sabati, M.; Sheriff, S.; Gu, M.; Wei, J.; Zhu, H.; Barker, P.B.; Spielman, D.M.; Alger, J.R.; Maudsley, A.A. Multivendor implementation and comparison of volumetric whole-brain echo-planar MR spectroscopic imaging. Magn. Reson. Med. 2015, 74, 1209–1220. [Google Scholar] [CrossRef] [Green Version]
  9. Goryawala, M.; Saraf-Lavi, E.; Nagornaya, N.; Heros, D.; Komotar, R.; Maudsley, A.A. The Association between Whole-Brain MR Spectroscopy and IDH Mutation Status in Gliomas. J. Neuroimaging 2020, 30, 58–64. [Google Scholar] [CrossRef] [Green Version]
  10. Goryawala, M.Z.; Sheriff, S.; Stoyanova, R.; Maudsley, A.A. Spectral decomposition for resolving partial volume effects in MRSI. Magn. Reson. Med. 2018, 79, 2886–2895. [Google Scholar] [CrossRef]
  11. Goryawala, M.Z.; Sheriff, S.; Maudsley, A.A. Regional distributions of brain glutamate and glutamine in normal subjects. NMR Biomed. 2016, 29, 1108–1116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Cordova, J.S.; Shu, H.K.; Liang, Z.; Gurbani, S.S.; Cooper, L.A.; Holder, C.A.; Olson, J.J.; Kairdolf, B.; Schreibmann, E.; Neill, S.G.; et al. Whole-brain spectroscopic MRI biomarkers identify infiltrating margins in glioblastoma patients. Neuro-Oncology 2016, 18, 1180–1189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ramesh, K.; Mellon, E.A.; Gurbani, S.S.; Weinberg, B.D.; Schreibmann, E.; Sheriff, S.A.; Goryawala, M.; de le Fuente, M.; Eaton, B.R.; Zhong, J.; et al. A multi-institutional pilot clinical trial of spectroscopic MRI-guided radiation dose escalation for newly diagnosed glioblastoma. Neurooncol. Adv. 2022, 4, vdac006. [Google Scholar] [CrossRef]
  14. Louis, D.N.; Perry, A.; Wesseling, P.; Brat, D.J.; Cree, I.A.; Figarella-Branger, D.; Hawkins, C.; Ng, H.; Pfister, S.M.; Reifenberger, G. The 2021 WHO classification of tumors of the central nervous system: A summary. Neuro-Oncology 2021, 23, 1231–1251. [Google Scholar] [CrossRef] [PubMed]
  15. Gurbani, S.S.; Schreibmann, E.; Maudsley, A.A.; Cordova, J.S.; Soher, B.J.; Poptani, H.; Verma, G.; Barker, P.B.; Shim, H.; Cooper, L.A.D. A convolutional neural network to filter artifacts in spectroscopic MRI. Magn. Reson. Med. 2018, 80, 1765–1775. [Google Scholar] [CrossRef]
  16. Maudsley, A.A.; Darkazanli, A.; Alger, J.R.; Hall, L.O.; Schuff, N.; Studholme, C.; Yu, Y.; Ebel, A.; Frew, A.; Goldgof, D.; et al. Comprehensive processing, display and analysis for in vivo MR spectroscopic imaging. NMR Biomed. 2006, 19, 492–503. [Google Scholar] [CrossRef] [Green Version]
  17. Maudsley, A.A.; Domenig, C.; Govind, V.; Darkazanli, A.; Studholme, C.; Arheart, K.; Bloomer, C. Mapping of brain metabolite distributions by volumetric proton MR spectroscopic imaging (MRSI). Magn. Reason. Med. 2009, 61, 548–559. [Google Scholar] [CrossRef] [Green Version]
  18. Veenith, T.V.; Mada, M.; Carter, E.; Grossac, J.; Newcombe, V.; Outtrim, J.; Lupson, V.; Nallapareddy, S.; Williams, G.B.; Sheriff, S.; et al. Comparison of inter subject variability and reproducibility of whole brain proton spectroscopy. PLoS ONE 2014, 9, e115304. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Taub, E.; Salibi, N.; Uswatte, G.; Maudsley, A.A.; Sheriff, S.; Womble, B.; Mark, V.W.; Knight, D.C. Comparison of reproducibility of single voxel spectroscopy and whole-brain magnetic resonance spectroscopy imaging at 3T. NMR Biomed. 2018, 31, e3898. [Google Scholar] [CrossRef]
  20. Cordova, J.S.; Schreibmann, E.; Hadjipanayis, C.G.; Guo, Y.; Shu, H.K.; Shim, H.; Holder, C.A. Quantitative tumor segmentation for evaluation of extent of glioblastoma resection to facilitate multisite clinical trials. Transl. Oncol. 2014, 7, 40–47. [Google Scholar] [CrossRef] [Green Version]
  21. Gurbani, S.; Weinberg, B.; Cooper, L.; Mellon, E.; Schreibmann, E.; Sheriff, S.; Maudsley, A.; Goryawala, M.; Shu, H.K.; Shim, H. The Brain Imaging Collaboration Suite (BrICS): A Cloud Platform for Integrating Whole-Brain Spectroscopic MRI into the Radiation Therapy Planning Workflow. Tomography 2019, 5, 184–191. [Google Scholar] [CrossRef] [PubMed]
  22. Gilbert, M.R.; Dignam, J.J.; Armstrong, T.S.; Wefel, J.S.; Blumenthal, D.T.; Vogelbaum, M.A.; Colman, H.; Chakravarti, A.; Pugh, S.; Won, M.; et al. A Randomized Trial of Bevacizumab for Newly Diagnosed Glioblastoma. N. Engl. J. Med. 2014, 370, 699–708. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Chinot, O.L.; Wick, W.; Mason, W.; Henriksson, R.; Saran, F.; Nishikawa, R.; Carpentier, A.F.; Hoang-Xuan, K.; Kavan, P.; Cernea, D.; et al. Bevacizumab plus Radiotherapy–Temozolomide for Newly Diagnosed Glioblastoma. N. Engl. J. Med. 2014, 370, 709–722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Stupp, R.; Taillibert, S.; Kanner, A.; Read, W.; Steinberg, D.M.; Lhermitte, B.; Toms, S.; Idbaih, A.; Ahluwalia, M.S.; Fink, K.; et al. Effect of Tumor-Treating Fields Plus Maintenance Temozolomide vs Maintenance Temozolomide Alone on Survival in Patients With Glioblastoma. JAMA 2017, 318, 2306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Bland, J.M.; Altman, D.G. Statistics Notes: Survival probabilities (the Kaplan-Meier method). BMJ 1998, 317, 1572–1580. [Google Scholar] [CrossRef] [Green Version]
  26. Davidson-Pilon, C. Lifelines: Survival analysis in Python. J. Open Source Softw. 2019, 4, 1317. [Google Scholar] [CrossRef] [Green Version]
  27. Cordova, J.S.; Gurbani, S.S.; Holder, C.A.; Olson, J.J.; Schreibmann, E.; Shi, R.; Guo, Y.; Shu, H.-K.G.; Shim, H.; Hadjipanayis, C.G. Semi-Automated Volumetric and Morphological Assessment of Glioblastoma Resection with Fluorescence-Guided Surgery. Mol. Imaging Biol. 2016, 18, 454–462. [Google Scholar] [CrossRef] [Green Version]
  28. Grabowski, M.M.; Recinos, P.F.; Nowacki, A.S.; Schroeder, J.L.; Angelov, L.; Barnett, G.H.; Vogelbaum, M.A. Residual tumor volume versus extent of resection: Predictors of survival after surgery for glioblastoma. J. Neurosurg. 2014, 121, 1115–1123. [Google Scholar] [CrossRef] [Green Version]
  29. Haj, A.; Doenitz, C.; Schebesch, K.-M.; Ehrensberger, D.; Hau, P.; Putnik, K.; Riemenschneider, M.; Wendl, C.; Gerken, M.; Pukrop, T.; et al. Extent of Resection in Newly Diagnosed Glioblastoma: Impact of a Specialized Neuro-Oncology Care Center. Brain Sci. 2017, 8, 5. [Google Scholar] [CrossRef] [Green Version]
  30. Awad, A.-W.; Karsy, M.; Sanai, N.; Spetzler, R.; Zhang, Y.; Xu, Y.; Mahan, M.A. Impact of removed tumor volume and location on patient outcome in glioblastoma. J. Neuro-Oncol. 2017, 135, 161–171. [Google Scholar] [CrossRef]
Figure 1. A comparison of an example gross tumor volume 3 (GTV3) treatment volume to the resection cavity with standard of care margins. Due to no residual enhancing tumor, standard of care would involve treating the cavity with a margin. (A) An axial view of the GTV3 contour (green) that extends past the resection cavity (blue), resection cavity with 5 mm margin (yellow), 10 mm margin (orange), and 20 mm margin (red). (B) This 3D view of the GTV3 contour shows a small amount of infiltrating tumor (arrow) that extends even beyond the 20 mm margin. (CE) 3D views of the resection cavity with various margins and volumes are calculated. A Dice overlap score was calculated between each margin and GTV3. For (CF), the Dice scores were 0.18, 0.47, 0.59, and 0.40, respectively. The percentages of volumes in (CF) not included in (B), were 37.4%, 41.3%, 49.7%, and 74.9%, respectively.
Figure 1. A comparison of an example gross tumor volume 3 (GTV3) treatment volume to the resection cavity with standard of care margins. Due to no residual enhancing tumor, standard of care would involve treating the cavity with a margin. (A) An axial view of the GTV3 contour (green) that extends past the resection cavity (blue), resection cavity with 5 mm margin (yellow), 10 mm margin (orange), and 20 mm margin (red). (B) This 3D view of the GTV3 contour shows a small amount of infiltrating tumor (arrow) that extends even beyond the 20 mm margin. (CE) 3D views of the resection cavity with various margins and volumes are calculated. A Dice overlap score was calculated between each margin and GTV3. For (CF), the Dice scores were 0.18, 0.47, 0.59, and 0.40, respectively. The percentages of volumes in (CF) not included in (B), were 37.4%, 41.3%, 49.7%, and 74.9%, respectively.
Cancers 15 03524 g001
Figure 2. Two sample cases are shown. In the first case, the first column shows the residual post-contrast enhancement (rENH) contour (pink outline) with a volume of 12.11 cc, overlaid on the T1-weighted contrast-enhanced (T1w-CE) MRI (first row) for a 54-year-old female glioblastoma patient. The second column shows the gross tumor volume 3 (GTV3) contour (pink outline), over 5 times larger than rENH with a volume of 63.37 cc. The second row is a three-dimensional rendering of both contour volumes overlaid on a T1 pre-contrast image. In the second case of a 56-year-old female glioblastoma patient, there is no residual contrast enhancement visible on the T1w-CE MRI, whereas the GTV3 volume is 18.45 cc.
Figure 2. Two sample cases are shown. In the first case, the first column shows the residual post-contrast enhancement (rENH) contour (pink outline) with a volume of 12.11 cc, overlaid on the T1-weighted contrast-enhanced (T1w-CE) MRI (first row) for a 54-year-old female glioblastoma patient. The second column shows the gross tumor volume 3 (GTV3) contour (pink outline), over 5 times larger than rENH with a volume of 63.37 cc. The second row is a three-dimensional rendering of both contour volumes overlaid on a T1 pre-contrast image. In the second case of a 56-year-old female glioblastoma patient, there is no residual contrast enhancement visible on the T1w-CE MRI, whereas the GTV3 volume is 18.45 cc.
Cancers 15 03524 g002
Figure 3. Kaplan–Meier survival curve of progression-free survival (PFS) for the 28 IDH wild-type patients included in this analysis. The median PFS is 16.2 months and median time to follow-up is 22.4 months.
Figure 3. Kaplan–Meier survival curve of progression-free survival (PFS) for the 28 IDH wild-type patients included in this analysis. The median PFS is 16.2 months and median time to follow-up is 22.4 months.
Cancers 15 03524 g003
Figure 4. Gross tumor volume 3 (GTV3) and residual post-contrast enhancement (rENH) scatter plot indicating low correlation between the two volumes based on R2 = 0.127 (p = 0.063). Each marker ‘x’ in the scatter plot corresponds to the GTV3 and rENH volume for each patient in the study. The horizontal blue line indicates the median GTV3 volume (19.2 cc), and the vertical orange line indicates the median rENH volume (2.1 cc).
Figure 4. Gross tumor volume 3 (GTV3) and residual post-contrast enhancement (rENH) scatter plot indicating low correlation between the two volumes based on R2 = 0.127 (p = 0.063). Each marker ‘x’ in the scatter plot corresponds to the GTV3 and rENH volume for each patient in the study. The horizontal blue line indicates the median GTV3 volume (19.2 cc), and the vertical orange line indicates the median rENH volume (2.1 cc).
Cancers 15 03524 g004
Figure 5. (A) Kaplan–Meier survival curve for overall survival (OS) stratified into equal groups by the median residual contrast enhancing (rENH) tumor volume (2.1 cc). The median OS for patients with rENH > 2.1 cc is 21.3 months, and the median for OS for patients with rENH < 2.1 cc is 23.0 months. (B) Kaplan–Meier survival curve for OS stratified into equal groups by the median gross tumor volume 3 (GTV3) (19.2 cc). The median OS for patients with GTV3 > 19.2 cc is 17.1 months, and the median OS for patients with GTV3 < 19.2 cc is 27.4 months. (C) Kaplan–Meier survival curve for progression-free survival (PFS) stratified into equal groups by the median rENH tumor volume (2.1 cc). The median PFS for patients with rENH > 2.1 cc is 12.8 months, and the median PFS for patients with rENH < 2.1 cc is 19.0 months. (D) Kaplan–Meier survival curve for PFS stratified into equal groups by the median GTV3 (19.2 cc). The median PFS for patients with GTV3 > 19.2 cc is 12.6 months, and the median PFS for patients with GTV3 < 19.2 cc is 24.0 months.
Figure 5. (A) Kaplan–Meier survival curve for overall survival (OS) stratified into equal groups by the median residual contrast enhancing (rENH) tumor volume (2.1 cc). The median OS for patients with rENH > 2.1 cc is 21.3 months, and the median for OS for patients with rENH < 2.1 cc is 23.0 months. (B) Kaplan–Meier survival curve for OS stratified into equal groups by the median gross tumor volume 3 (GTV3) (19.2 cc). The median OS for patients with GTV3 > 19.2 cc is 17.1 months, and the median OS for patients with GTV3 < 19.2 cc is 27.4 months. (C) Kaplan–Meier survival curve for progression-free survival (PFS) stratified into equal groups by the median rENH tumor volume (2.1 cc). The median PFS for patients with rENH > 2.1 cc is 12.8 months, and the median PFS for patients with rENH < 2.1 cc is 19.0 months. (D) Kaplan–Meier survival curve for PFS stratified into equal groups by the median GTV3 (19.2 cc). The median PFS for patients with GTV3 > 19.2 cc is 12.6 months, and the median PFS for patients with GTV3 < 19.2 cc is 24.0 months.
Cancers 15 03524 g005
Table 1. Classification of all 28 patients based on extent of resection (EOR) and residual post-contrast enhancement (rENH) volume. Patients with EOR ≥ 95% had gross total resection (GTR), EOR < 95% had subtotal resection (STR), and the remaining 5 patients had biopsy only. Patients with rENH ≤ 1.5 cc also had GTR, and rENH > 1.5 cc had STR.
Table 1. Classification of all 28 patients based on extent of resection (EOR) and residual post-contrast enhancement (rENH) volume. Patients with EOR ≥ 95% had gross total resection (GTR), EOR < 95% had subtotal resection (STR), and the remaining 5 patients had biopsy only. Patients with rENH ≤ 1.5 cc also had GTR, and rENH > 1.5 cc had STR.
ClassificationNumber (%) of Patients
EORGTR (EOR ≥ 95%)12 (42.9%)
STR (EOR < 95%)11 (39.3%)
rENHGTR (rENH ≤ 1.5 cc)13 (46.4%)
STR (rENH > 1.5 cc)15 (53.6%)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Trivedi, A.G.; Ramesh, K.K.; Huang, V.; Mellon, E.A.; Barker, P.B.; Kleinberg, L.R.; Weinberg, B.D.; Shu, H.-K.G.; Shim, H. Spectroscopic MRI-Based Biomarkers Predict Survival for Newly Diagnosed Glioblastoma in a Clinical Trial. Cancers 2023, 15, 3524. https://doi.org/10.3390/cancers15133524

AMA Style

Trivedi AG, Ramesh KK, Huang V, Mellon EA, Barker PB, Kleinberg LR, Weinberg BD, Shu H-KG, Shim H. Spectroscopic MRI-Based Biomarkers Predict Survival for Newly Diagnosed Glioblastoma in a Clinical Trial. Cancers. 2023; 15(13):3524. https://doi.org/10.3390/cancers15133524

Chicago/Turabian Style

Trivedi, Anuradha G., Karthik K. Ramesh, Vicki Huang, Eric A. Mellon, Peter B. Barker, Lawrence R. Kleinberg, Brent D. Weinberg, Hui-Kuo G. Shu, and Hyunsuk Shim. 2023. "Spectroscopic MRI-Based Biomarkers Predict Survival for Newly Diagnosed Glioblastoma in a Clinical Trial" Cancers 15, no. 13: 3524. https://doi.org/10.3390/cancers15133524

APA Style

Trivedi, A. G., Ramesh, K. K., Huang, V., Mellon, E. A., Barker, P. B., Kleinberg, L. R., Weinberg, B. D., Shu, H. -K. G., & Shim, H. (2023). Spectroscopic MRI-Based Biomarkers Predict Survival for Newly Diagnosed Glioblastoma in a Clinical Trial. Cancers, 15(13), 3524. https://doi.org/10.3390/cancers15133524

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop