Next Article in Journal
Efficacy of Oncolytic Herpes Simplex Virus T-VEC Combined with BET Inhibitors as an Innovative Therapy Approach for NUT Carcinoma
Next Article in Special Issue
In Vivo Quantitative Imaging of Glioma Heterogeneity Employing Positron Emission Tomography
Previous Article in Journal
Mindfulness in Adolescent and Young Adult (AYA) Patients Undergoing Hematopoietic Stem Cell Transplantation (HSCT): A Qualitative Study
Previous Article in Special Issue
A Systematic Review of the Current Status and Quality of Radiomics for Glioma Differential Diagnosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Validation Study for Non-Invasive Prediction of IDH Mutation Status in Patients with Glioma Using In Vivo 1H-Magnetic Resonance Spectroscopy and Machine Learning

by
Elisabeth Bumes
1,*,
Claudia Fellner
2,
Franz A. Fellner
3,
Karin Fleischanderl
4,
Martina Häckl
5,
Stefan Lenz
4,
Ralf Linker
1,
Tim Mirus
5,
Peter J. Oefner
5,
Christian Paar
6,
Martin Andreas Proescholdt
7,
Markus J. Riemenschneider
8,
Katharina Rosengarth
7,
Serge Weis
9,
Christina Wendl
2,
Sibylle Wimmer
10,
Peter Hau
1,
Wolfram Gronwald
5,† and
Markus Hutterer
1,11,†
1
Department of Neurology and Wilhelm Sander-NeuroOncology Unit, Regensburg University Hospital, 93055 Regensburg, Germany
2
Department of Radiology and Division of Neuroradiology, Regensburg University Hospital, 93055 Regensburg, Germany
3
Central Institute of Radiology, Kepler University Hospital, 4021 Linz, Austria
4
Division of Molecular Pathology, Neuromed Campus, Kepler University Hospital, 4020 Linz, Austria
5
Institute of Functional Genomics, University of Regensburg, 93053 Regensburg, Germany
6
Institute of Laboratory Medicine, Kepler University Hospital, 4021 Linz, Austria
7
Department of Neurosurgery, Regensburg University Hospital, 93053 Regensburg, Germany
8
Department of Neuropathology, Regensburg University Hospital, 93053 Regensburg, Germany
9
Division of Neuropathology, Neuromed Campus, Kepler University Hospital, 4020 Linz, Austria
10
Institute of Neuroradiology, Neuromed Campus, Kepler University Hospital, 4020 Linz, Austria
11
Department of Neurology with Acute Geriatrics, Saint John of God Hospital Linz, 4021 Linz, Austria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2022, 14(11), 2762; https://doi.org/10.3390/cancers14112762
Submission received: 30 April 2022 / Revised: 24 May 2022 / Accepted: 31 May 2022 / Published: 2 June 2022
(This article belongs to the Special Issue In Vivo Quantitative Imaging of Gliomas)

Abstract

:

Simple Summary

The enzyme isocitrate dehydrogenase (IDH) affects glioma cell metabolism in multiple ways. Mutation of IDH is not only indicative of the presence of astrocytoma or oligodendroglioma but it also comes with a better prognosis and constitutes a promising therapeutic target. Therefore, determination of IDH mutation status is essential in clinical practice. In most patients, tissue can be obtained by resection or biopsy to determine IDH status histologically. However, in some cases, this is not possible for technical reasons. We recently showed in a small cohort of patients that non-invasive determination of IDH mutation status using proton magnetic resonance spectroscopy (1H-MRS) at 3.0 Tesla (T) together with machine learning techniques is feasible in a standard clinical setting and with acceptable effort. Here, we demonstrate that our approach showed comparably good results in sensitivity (82.6%) and specificity (72.7%) in a larger validation cohort employing 1H-MRS at 1.5 T in a retrospective, distinct setting. We concluded that our method works well regardless of the magnetic field strength and scanner used, and thus, may improve patient care.

Abstract

The isocitrate dehydrogenase (IDH) mutation status is an indispensable prerequisite for diagnosis of glioma (astrocytoma and oligodendroglioma) according to the WHO classification of brain tumors 2021 and is a potential therapeutic target. Usually, immunohistochemistry followed by sequencing of tumor tissue is performed for this purpose. In clinical routine, however, non-invasive determination of IDH mutation status is desirable in cases where tumor biopsy is not possible and for monitoring neuro-oncological therapies. In a previous publication, we presented reliable prediction of IDH mutation status employing proton magnetic resonance spectroscopy (1H-MRS) on a 3.0 Tesla (T) scanner and machine learning in a prospective cohort of 34 glioma patients. Here, we validated this approach in an independent cohort of 67 patients, for which 1H-MR spectra were acquired at 1.5 T between 2002 and 2007, using the same data analysis approach. Despite different technical conditions, a sensitivity of 82.6% (95% CI, 61.2–95.1%) and a specificity of 72.7% (95% CI, 57.2–85.0%) could be achieved. We concluded that our 1H-MRS based approach can be established in a routine clinical setting with affordable effort and time, independent of technical conditions employed. Therefore, the method provides a non-invasive tool for determining IDH status that is well-applicable in an everyday clinical setting.

Graphical Abstract

1. Introduction

Gliomas are among the most common primary brain tumors in adults [1] and are mostly not curable. Despite a range of treatment options, such as tumor resection, radiotherapy, systemic therapies (e.g., alkylating agents), experimental targeted therapies, and tumor treating fields, the prognosis is generally poor [2,3,4].
According to the WHO classification of brain tumors 2021, the essential discriminating feature of all gliomas is the isocitrate dehydrogenase (IDH) mutation status [5]. The presence of an IDH mutated (IDHmut) status indicates a CNS WHO grade 2–4 astrocytoma or oligodendroglioma, whereas IDH wildtype (IDHwt) status in most cases implies the diagnosis of glioblastoma CNS WHO grade 4. Therefore, IDH status has a high value for the classification of gliomas. In addition, it is highly relevant for prognostication. While patients suffering from IDHwt glioblastoma receiving standard treatment have an average overall survival of 15–18 months [6,7], average overall survival can be up to 14.7 years in patients with oligodendroglioma in the presence of an IDH mutation and combined loss of heterozygosity (LOH) 1p/19q under sufficient treatment [8,9]. In addition, IDHmut is a putative, specific treatment target and has been investigated in several clinical trials using peptide vaccination and small-molecule inhibitor approaches [10,11]. Finally, repetitive evaluation of IDHmut might be a valuable tool for monitoring tumor response and relapse. However, sequential biopsy is not feasible in patients with glioma.
Non-invasive serial determination of IDH status by proton magnetic resonance spectroscopy (1H-MRS) can, therefore, be a valuable additional tool, supplementing standard immunohistochemical diagnosis for IDH status and treatment monitoring [3].
For these reasons, much effort has been put in the development of imaging methods that can reliably predict IDH mutation status in a clinically feasible manner. One focus is on the development of specialized applications of proton magnetic resonance imaging (1H-MRI) [12]. Good prediction values are reported in the analysis of standard 1H-MRI sequences using convoluted neuronal networks (CNNs). Employing regional cerebral blood volume (rCBV) maps, an accuracy of 83% in predicting IDH status in glioblastoma was obtained [13]. However, processing of radiomic features using artificial intelligence and machine learning techniques often ends up with inconclusive results, requires high technical expertise, and is time consuming in manual preparation [14,15].
In contrast, 1H-MRS is a well-established method to gain information about the metabolic composition of a brain lesion and is also suitable for detection of the IDHmut-derived metabolite 2-hydroxyglutarate (2-HG) at 2.24 ppm [16]. Current applications of 1H-MRS for the 2-HG based prediction of the IDH mutation status achieve a pooled sensitivity and specificity of 95% and 91%, respectively [16]. However, due to personnel, technical, and time requirements, these procedures are usually implemented only in a research setting. Usually, scanners with 3.0 Tesla (T) or 7.0 T and long examination times, which must be tolerated by the patient, are required [17,18]. This also means a high personnel effort both in acquisition and evaluation of the spectra. Furthermore, detection of the 2-HG peak at 2.24 ppm using a 3.0 T 1H-MRS is complicated by overlap with other metabolites, such as glutamine and glutamate or lipids [19]. These issues may be resolved with high technical effort such as two-dimensional measurements [20,21] or analysis of several metabolites, e.g., 2-HG and glutamate [22].
Oligodendroglioma and most astrocytoma belong to the group of IDHmut glioma and therefore show the signals of 2-HG in their spectra. Furthermore, LOH 1p/19q is a characteristic signature of oligodendroglioma that generally does not differ markedly in metabolism from IDHmut, non-co-deleted astrocytoma, except for elevated cystathionine levels that may be helpful in distinguishing IDH mutated gliomas with and without LOH 1p/19q [23]. Generally, 1H-MRS data for oligodendroglioma are sparse, but also suggest that 1H-MRS may be helpful in grading [24].
In our previous prospective study cohort (original cohort), we focused on a full spectrum analysis of in vivo 1H-MRS at 3.0 T in combination with a support vector machine (SVM) in 34 glioma patients. With this advanced diagnostic approach, we were able to achieve a specificity of 75.0% and a sensitivity of 95.5% for non-invasive determination of IDH mutation status [25]. This method demonstrated good technical feasibility and was well-tolerated by the patients with an additional examination period of approximately 20 min, making its use feasible in a standard clinical setting.
For further validation, we applied our approach of a full spectrum analysis in the present study to a larger retrospective cohort of patients with gliomas assessed by 1H-MRS at 1.5 T (validation cohort). The aim of this study was to investigate whether our established approach can be transferred to another environment with comparable specificity and sensitivity.

2. Materials and Methods

2.1. Validation Cohort

For this retrospective analysis, 100 preoperative spectra of patients with suspected glioma were available. 1H-MRS had been performed at 1.5 T at a single neuro-oncological site (Neuromed Campus, Kepler University Hospital, Linz, Austria) between 2002 and 2007. Spectra were excluded from further analysis if no tumor was diagnosed (n = 1) or IDH status could not be determined (n = 6). In these six cases, no IDH status was known beforehand, and no tissue was available for neuropathological reassessment. Additionally, quality of spectra was insufficient in n = 26 (Figure 1).
The mandatory 1H-MRS quality requirements were predefined for the amplitudes of the following metabolites: choline (Cho) > 0.2, creatine (Cr) > 0.1, myo-inositol (M-Ins) > 0, and creatine 2 (Cr2) > 0. Note that amplitudes are given in arbitrary units as defined by the vendor, and therefore, thresholds of amplitudes might not be applicable to other cohorts. Details and examples of sufficient and insufficient spectrum quality are shown in the supplement (Figure S1). Ultimately, 67 spectra were included from 44 patients with IDHwt and 23 patients with IDHmut status.
For 32 samples, IDH status had been assessed previously in clinical routine diagnostics by immunohistochemistry for the most common IDH mutation IDH1 R132H [26]. Neuropathological re-evaluation with next-generation sequencing (NGS) for this study was performed in 35 cases described in detail in Section 2.2. Thirteen patients with immunohistochemically confirmed IDHwt status were less than 55 years old at initial diagnosis. According to the current diagnostic guidelines [5,27], additional IDH sequencing was performed to detect rare IDH mutations particularly in these cases. In nine samples, IDHwt status was confirmed, whereas in four cases, DNA sequencing was not possible (no tissue available: n = 3, DNA concentration too low: n = 1).
From these four patients, three had a predefined, histologically confirmed glioblastoma, WHO grade IV, without notice of a lower-grade precursor lesion. Thus, these cases could be most likely considered primary IDHwt glioblastomas. The remaining sample was a histologically classical malignant diffusely infiltrating astrocytic glioma, most likely corresponding to a molecularly defined IDHwt glioblastoma according to the WHO 2021 classification.

2.2. Neuropathological Assessment of Glioma Tissue

2.2.1. Immunohistochemistry

Immunohistochemistry was performed on formalin-fixed and paraffin-embedded (FFPE) 5μM thick sections on Superfrost Plus Adhesion Microscope slides (J1800AMN2, Epredia, Kalamazoo, MI, USA). Deparaffinized, rehydrated sections underwent antigen retrieval in the PT Module Lab Vision (Thermo Scientific, Loughborough, UK) using the Dewax and HIER Buffer L (pH 6) (Epredia, DA Breda, The Netherlands). All subsequent steps were carried out using the Lab Vision Autostainer 480S (Thermo Scientific, Loughborough, UK) and the Ultravision Quanto Detection System HRP (TL-125-QHL, Epredia, DA Breda, The Netherlands). The reaction product was visualized using diaminobenzidine chromogen (DAB Quanto, TA-125-QDX, Epredia, DA Breda, The Netherlands). Sections were counterstained with Mayer’s hematoxylin (Thermo Scientific, Loughborough, UK). The following primary antibody was used: anti-IDH1 R132H (dilution 1:50, mouse monoclonal, clone H09, Dianova, Hamburg, Germany). As negative control, the primary antibody was omitted and was replaced by normal mouse serum.

2.2.2. DNA Extraction and Quantification

DNA extraction from FFPE tissue was performed with the blackPREP FFPE DNA Kit (Analytik Jena, Jena, Germany) and quantified with a Qubit 3.0 System using the Qubit™ dsDNA HS Assay Kit (Invitrogen, Dreieich, Germany).

2.2.3. Next-Generation Sequencing (NGS)

The IDH mutation status of samples with unclear immunohistochemistry and of patients younger than 55 years with negative IDH immunohistochemistry was determined by NGS using the Ion AmpliSeq™ Comprehensive Cancer Panel (Life Technologies, Darmstadt, Germany). NGS libraries were generated with the Ion AmpliSeqTM Library 2.0 Kit (Life Technologies, Darmstadt, Germany, User guide MAN0006735 Rev. B.0). DNA was pretreated for 30 min with 1 unit Uracil-DNA Glycosylase (Life Technologies, Darmstadt, Germany). Template preparation, chip loading, and sequencing were performed with the Ion Chef™ and the Ion GeneStudio™ S5 Plus Systems according to the Ion 510™ & Ion 520™ & Ion 530™ Kit—Chef user guide (MAN0016854 Rev. E.0, Life Technologies, Darmstadt, Germany). Data were analyzed with Torrent Suite™ software version 5.12.2 and Ion Reporter software version 5.16. Read alignments at the IDH1 and IDH2 hotspots were visually inspected with the Integrative Genomics Viewer (Broad Institute, Cambridge, MA, USA).

2.2.4. Multiplex Ligation-Dependent Probe Amplification (MLPA)

The copy number status of chromosomes 1p and 19q was determined by MLPA using the P088 Oligodendroglioma 1p-19q probe mix (MRC Holland, Amsterdam, Netherlands, Lot. D1-1119, COA Version 01) according to the MLPA General Protocol MDP version-007. Data were analyzed with Coffalyser.Net software version v.210604.1451 (MRC Holland, Amsterdam, Netherlands).

2.3. Data Acquisition of 1H-MRS

In the validation cohort, magnetic resonance imaging (MRI) and 1H-MRS were performed on clinical whole-body scanners, as performed in the original cohort. While the data for the original cohort were acquired on a state-of-the art 3.0 T scanner (MAGNETOM Skyra, Siemens Healthcare, Erlangen, Germany) between 2015 and 2019, data for the validation cohort had been measured between 2002 and 2007 on a 1.5 T scanner (MAGNETOM Symphony, Siemens Healthcare, Erlangen, Germany). MR spectra were acquired by a single-voxel technique employing a point-resolved spectroscopy (PRESS) sequence with a short echo time (TE) of 30 ms.
There were two major differences between the training and validation cohort: the voxel size for 1H-MRS was kept constant in the original cohort, while the voxel size was adapted to tumor size and shape in the validation cohort. In addition, positioning of the voxel into the most representative part of the tumor was guided exclusively by standard MRI for the validation cohort, while voxel positioning had been also supported by O-(2-18F-fluoroethyl)-L-tyrosine (18F-FET) positron emission tomography (PET) in the original cohort. Further technical details are given in the supplement (Table S1) and for the original cohort in Bumes et al., 2020 [25].

2.4. H-MRS Data Preprocessing

For the original cohort, in each spectrum, 1024 data points were acquired over a spectral width of 1200 Hz at 3.0 T, while for the validation cohort, 1024 data points were recorded over a spectral width of 1000 Hz at 1.5 T. Further technical details are given in the supplement (Text S1).
As the predictive model was trained on the original data, before application to the validation data (Figure 2), data sets had to be matched to each other. To this end, all spectra were aligned with respect to the choline peak while accounting for the differences in field strength and spectral width between the two data sets.

2.5. H-MRS Classification

First, all 1H-MR spectra of the original cohort were used for training of an SVM with linear kernel. To this end, the R library e1071 (http://cran.r-project.org/web/packages/e1071/(accessed on 22 February 2022)) was applied. Throughout training, the cost parameter was left on its default value of one. Next, the preprocessed and adjusted spectra of the validation cohort were presented to the trained model as test cases.
The complete data acquisition and analysis pipeline is summarized in Figure 2, and a detailed description is presented in the supplement (Text S2).
We have developed a user-friendly R-shiny app that allows application of the trained model to additional data. New data are automatically adapted with respect to field strength, spectral width, and number of acquired data points. The app including the fully trained model predicts IDH mutation status and is freely accessible for scientific purposes: https://www.uni-regensburg.de/medicine/functional-genomics/staff/prof-wolfram-gronwald/software/index.html (accessed on 29 April 2022) (Figure S2). In addition, the app can also be used as a web application: https://IDH-prediction.spang-lab.de (accessed on 29 April 2022).

3. Results

3.1. Patient Characteristics of the Validation Cohort

Among the 67 patients enrolled, 44 cases were IDHwt, whereas 21 patients carried an IDH1 R132H mutation (Figure 1). Additionally, two rare IDH mutations (R132S and R172M) were detected, bringing the total to 23 IDHmut patients.
Age at initial diagnosis was higher in patients with IDHwt glioma, with a median of 61 years (range 21 to 83 years), than in patients with IDHmut glioma, with a median of 37 years (range 18 to 69 years). The exact age distribution at first diagnosis is shown in the supplement (Figure S3).
As expected, patients with IDHwt status usually showed WHO grade IV in the initial histological assessment according to the WHO classification of brain tumors 2000 (Table 1).
In total, 30% of patients with IDHmut status (n = 7) presented with LOH 1p/19q, thus confirming an oligodendroglioma according to the WHO classifications 2016 and 2021. No LOH 1p/19q was present in ten samples (43%), while assessment was not possible in six cases (26%) due to technical limitations.
Since an IDHwt status excludes the diagnosis of oligodendroglioma, we did not determine LOH status in IDHwt samples in general.
Patient characteristics of the original cohort are published in Bumes et al. [25].

3.2. IDH Mutation Caused Specific Alterations in 1H-MR Spectra

Distinct alterations characteristic for the presence of an IDH mutation were apparent in the spectra beyond the 2-HG peak at 2.24 ppm, especially in the region between 4.1 ppm–3.5 ppm, as already reported for the original cohort [25]. In two exemplary spectra of the validation cohort (Figure 3A,B) acquired with a routine 1.5 T scanner and without standardized voxel sizes, this region is marked by a blue box including an increase in the region of the M-Ins signal at 3.53 ppm in the IDHmut spectrum (Figure 3A). In Figure 3C, the averaged spectra of the IDHmut and IDHwt group show the variation in the region of interest even more clearly.
Thus, the spectral differences observed in the original cohort held true for the validation cohort. Furthermore, clear differences between the IDHmut and IDHwt group were apparent in the spectral region corresponding to lactate and methylene (CH2) groups of lipids at around 1.3 ppm, showing a considerably increased intensity in the IDHwt group (Figure 3C).

3.3. A Linear SVM Provided High Sensitivity and Specificity in Detecting IDHmut

Next, a linear SVM was applied to the classification of 1H-MR spectra from the validation cohort. Employing the original cohort, the number of discriminatory features was optimized in a cross-validation for an optimal prediction accuracy, resulting in the selection of two predictive features. The final model was subsequently trained on the complete original cohort, using the two most differentially expressed features between the IDHwt and IDHmut groups. Following matching of the validation and original data that were independently acquired at 1.5 T and 3.0 T, respectively, the trained model was applied to achieve classification of the validation data in the next step, as described in the Materials and Methods section.
We showed an accuracy of 76.1% (95% CI, 64.1–85.7%), a sensitivity of 82.6% (95% CI, 61.2–95.1%), and a specificity of 72.7% (95% CI, 57.2–85.0%). The area under the ROC curve (AUC) amounted to 0.82 (Figure 4). A full description of the obtained statistical parameters for both the original and the validation cohort is given in Table 2.
Visual inspection of false positive and false negative predicted spectra revealed no evidence of systematic error. Quality of the spectra was sufficient in all false predicted samples; examples are shown in the supplement (Figure S4).

3.4. SVM Did Not Predict LOH 1p/19q Status

Subsequently, we investigated whether the non-invasive prediction of the LOH 1p/19q status was possible. To this end, an SVM was trained to separate the spectra of IDHmut subgroups with and without combined LOH 1p/19q. According to the current diagnostic standards, a combined LOH 1p/19q is present only in IDHmut gliomas, so the neuropathological re-evaluation of the LOH status was performed only in IDHmut samples as described in Section 2.2.4.
The small group sizes of the validation cohort (Figure 5A), due to technical reasons caused by the age of the samples, did not allow a statistically meaningful analysis of the spectra. With the same technical assessment in both cohorts, it was feasible to add the samples from the original cohort to the validation cohort to obtain a larger group (total cohort: IDHmut with LOH 1p/19q, n = 13; IDHmut without LOH 1p/19q, n = 26).
As demonstrated in the averaged spectra (Figure 5B), only minor differences in the full spectrum analysis of the total cohort were visible between IDHmut subgroups with and without LOH 1p/19q. Results of the used classification approach showed that the observed differences did not allow reliable prediction of the LOH 1p/19q status (not shown).

4. Discussion

Here, we showed that assessment of IDH mutation status using single-voxel 1H-MRS on a 1.5 T scanner produced acceptable sensitivity (82.6%), specificity (72.7%), and AUC (0.82), if compared to tissue-based methods. These findings corroborate data from our recently published original cohort, which yielded a higher sensitivity (95.5%) but similar specificity (75.0%) and AUC (0.83).
The data provided in Table 2 show that for both cohorts, reliable classification results were obtained. This is especially noteworthy as disease prevalence differed considerably between the two cohorts. In the original cohort used for training of the model, 22 out of 34 patients were bearing an IDH mutation (65%), while in the validation, this was true for only 23 out of 67 patients (34%). This further demonstrates the robustness of the proposed approach.
The good reproducibility is further remarkable given the different technical conditions. The single-voxel spectroscopy in our validation cohort, which used spectra generated from 2002 to 2007, was performed on a 1.5 T scanner, while the original cohort was examined on a modern 3.0 T scanner from 2015 to 2020. Although similar measurement parameters were used in both cohorts, the 1H-MRS at 1.5 T yielded, compared to 3.0 T, a reduced spectral resolution and a lower signal-to-noise ratio (SNR). Therefore, the technical conditions of the validation cohort were inferior to those of the original cohort.
Furthermore, differences in the positioning and configuration of the voxel were evident for both cohorts. In our original cohort, 18F-FET PET was used as an additional tool for voxel positioning, which was not available for the validation cohort. Therefore, voxel positioning in the validation cohort might not always have captured the region of highest biological activity, further impeding sample classification.
While in the original cohort, a constant isotropic voxel length of 14 mm was used to obtain comparable spectra, the voxel size in the validation cohort was adjusted to the region of interest (based on structural MRI) to achieve the best possible coverage of the tumor region. This approach causes minimal partial volume effects from surrounding tissue. However, for small voxel size, the number of averages was not increased due to time reasons, resulting in a relatively low SNR in some cases. To exclude spectra with too low SNR, quality criteria for amplitudes of metabolites were introduced (Figure 1 and Figure S1). On the other hand, due to the use of a constant isotropic voxel length of 14 mm, there was a higher risk of partial volume effects in the original cohort.
Despite these substantial differences, our approach of full spectrum analysis proved to be sufficiently robust to assess IDH mutation status. Hence, the approach can be applied easily to other data sets acquired with different technical settings. The technical requirements for determination of IDH mutation status were also met with 1.5 T scanners, which are frequently used in clinical practice. Furthermore, the additional examination time of approximately 20 min is tolerable for most patients. Note that even better results may be expected with the use of 3.0 T scanners, which are becoming increasingly available in clinical routine.
A meta-analysis of 14 studies comprising 460 glioma patients showed excellent diagnostic performance of 2-HG 1H-MRS for prediction of IDHmut glioma, with a pooled sensitivity and specificity of 95% (95% CI, 85–98%) and 91% (95% CI, 83–96%), respectively [16]. However, determination of the 2-HG peak is technically much more elaborate and difficult to implement in clinical routine. Usually, a scanner with at least 3.0 T, often 7.0 T, is used with long examination times that are not tolerable for neurologically impaired patients. Berrington et al. [17] showed that reliable detection of 2-HG is possible using an ultra-high field 7.0 T scanner. Similarly, Shen et al. [18] employed a 7.0 T scanner for the detection of 2-HG. Additionally, two-dimensional MR spectra have been successfully applied for the reliable detection of 2-HG [20,21].
Other approaches, such as the use of radiomic features, employ standard MRI sequences and therefore skip the additional examination time for 1H-MRS. However, the prediction quality depends very much on the algorithm used [14]. Furthermore, it is a technically complex and time-consuming process (e.g., tumor segmentation is needed) difficult to be implemented in clinical routine.
In comparison, despite its simplicity and lower technical requirements, our approach using full spectrum analysis still allowed reliable prediction of IDH mutation status with both good sensitivity and specificity. We encourage colleagues to use our free app to predict IDH mutation status for scientific purposes (https://www.uni-regensburg.de/medicine/functional-genomics/staff/prof-wolfram-gronwald/software/index.html (accessed on 29 April 2022)) (Figure S2), which is also available as a web application (https://IDH-prediction.spang-lab.de (accessed on 29 April 2022)).
Published data using a MEGA-PRESS sequence on a 3.0 T scanner showed that the spectra of IDH mutant gliomas with and without LOH 1p/19q are distinct, especially in cystathionine levels [23]. We could not replicate these alterations in our analysis, probably because of different technical conditions, as we used a PRESS sequence on a 1.5 T scanner. Furthermore, only a small number of histologically confirmed spectra with LOH 1p/19q (n = 13) could be included in our analysis.
This study has several weaknesses. First, both the original and the validation cohort are relatively small. We had to exclude some spectra of the validation cohort either due to unavailability of IDH status or poor spectra quality.
Neuropathological reassessment of IDH status by sequencing was not possible for the entire validation cohort. Previously, IDH status was determined by immunohistochemistry. Here, if feasible, IDH status was confirmed by sequencing, particularly in patients with IDHwt status who were younger than 55 years at initial diagnosis. However, due to the time passed, tissue was not always available or technical problems prevented neuropathological reassessment. We believe the risk of missing a rare IDH mutation in patients with IDHwt immunohistochemistry but no retroactive sequencing to be very low. As outlined above, the clinical information and/or the histological pattern in these cases were consistent with IDHwt glioblastoma and no suspicion of misclassification was found in the analysis of the spectra.
This study also has several strengths. Despite the different conditions used to acquire 1H-MR spectra in the two cohorts, diagnostic performance of our approach, in particular with regard to specificity, was sufficiently robust to predict IDH mutation status with acceptable effort of time and human resources.

5. Conclusions

We conclude that assessment of IDH mutation status is feasible in a standard clinical setting with little time and staff effort using 1H-MRS on routine 1.5 T and 3.0 T MR scanners. The method may, therefore, be of value for the assessment of IDH status in patients where biopsy or resection is not feasible, and in a setting where sequential performance is required. Thus, our approach works independently of various technical requirements and should be further developed to enable the non-invasive determination of IDH mutation status in clinical practice.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14112762/s1, Figure S1: Quality of 1H-MR spectra; Table S1: Technical parameters and conditions for the original and validation cohort; Text S1: Matching of the original and the validation cohort; Text S2: Workflow of data collection and analysis in the original and validation cohort; Figure S2: Screenshot of our App predicting IDH mutation status; Figure S3: Patient characteristics: Distribution of age at first diagnosis in the IDHwt group (n = 44) and the IDHmut group (n = 23); Figure S4: Examples of false positive and false negative predicted 1H-MR spectra.

Author Contributions

Conceptualization, W.G. and M.H. (Markus Hutterer); methodology, C.F., F.A.F., K.F., M.H. (Martina Häckl), M.H. (Markus Hutterer), S.L., T.M., C.P., K.R., S.W. (Serge Weis), C.W., S.W. (Sibylle Wimmer) and W.G.; software, M.H. (Martina Häckl) and W.G.; validation, E.B., W.G., P.H. and M.H. (Markus Hutterer); formal analysis, E.B., C.F., S.L., S.W. (Sibylle Wimmer) and W.G.; investigation, C.F., F.A.F., S.L., C.P., S.W. (Sibylle Wimmer) and W.G.; resources, F.A.F., C.W., R.L. and S.W. (Serge Weis); data curation, E.B., W.G. and M.H. (Markus Hutterer); writing—original draft preparation, E.B.; writing—review and editing, C.F., F.A.F., K.F., M.H. (Martina Häckl), M.H. (Markus Hutterer), S.L., R.L., T.M., P.J.O., C.P., M.A.P., M.J.R., K.R., S.W. (Serge Weis), C.W., S.W. (Sibylle Wimmer), P.H. and W.G.; visualization, E.B., C.F. and W.G.; supervision, P.H., W.G. and M.H. (Markus Hutterer); project administration, E.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German Research Foundation, KFO 262.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the local Ethics Committee (Upper Austria, protocol code EK 1076/2020, date of approval 20 May 2020).

Informed Consent Statement

Due to the retrospective nature of the study, informed consent of the patients was not performed.

Data Availability Statement

Data are contained within the article or Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ostrom, Q.T.; Cioffi, G.; Waite, K.; Kruchko, C.; Barnholtz-Sloan, J.S. CBTRUS Statistical Report: Primary Brain and Other Central Nervous System Tumors Diagnosed in the United States in 2014–2018. Neuro-Oncology 2021, 23, iii1–iii105. [Google Scholar] [CrossRef]
  2. Wen, P.Y.; Weller, M.; Lee, E.Q.; Alexander, B.M.; Barnholtz-Sloan, J.S.; Barthel, F.P.; Batchelor, T.T.; Bindra, R.S.; Chang, S.M.; Chiocca, E.A.; et al. Glioblastoma in adults: A Society for Neuro-Oncology (SNO) and European Society of Neuro-Oncology (EANO) consensus review on current management and future directions. Neuro-Oncology 2020, 22, 1073–1113. [Google Scholar] [CrossRef]
  3. Weller, M.; van den Bent, M.; Preusser, M.; Le Rhun, E.; Tonn, J.C.; Minniti, G.; Bendszus, M.; Balana, C.; Chinot, O.; Dirven, L.; et al. EANO guidelines on the diagnosis and treatment of diffuse gliomas of adulthood. Nat. Rev. Clin. Oncol. 2021, 18, 170–186. [Google Scholar] [CrossRef]
  4. Stupp, R.; Taillibert, S.; Kanner, A.; Read, W.; Steinberg, D.; 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: A Randomized Clinical Trial. JAMA 2017, 318, 2306–2316. [Google Scholar] [CrossRef] [Green Version]
  5. Louis, D.N.; Perry, A.; Wesseling, P.; Brat, D.J.; Cree, I.A.; Figarella-Branger, D.; Hawkins, C.; Ng, H.K.; Pfister, S.M.; Reifenberger, G.; et al. The 2021 WHO Classification of Tumors of the Central Nervous System: A summary. Neuro-Oncology 2021, 23, 1231–1251. [Google Scholar] [CrossRef]
  6. Stupp, R.; Mason, W.P.; van den Bent, M.J.; Weller, M.; Fisher, B.; Taphoorn, M.J.; 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]
  7. Gilbert, M.R.; Yuan, Y.; Wu, J.; Mendoza, T.; Vera, E.; Omuro, A.; Lieberman, F.; Robins, H.I.; Gerstner, E.R.; Wu, J.; et al. A phase II study of dose-dense temozolomide and lapatinib for recurrent low-grade and anaplastic supratentorial, infratentorial, and spinal cord ependymoma. Neuro-Oncology 2021, 23, 468–477. [Google Scholar] [CrossRef]
  8. Cairncross, G.; Wang, M.; Shaw, E.; Jenkins, R.; Brachman, D.; Buckner, J.; Fink, K.; Souhami, L.; Laperriere, N.; Curran, W.; et al. Phase III trial of chemoradiotherapy for anaplastic oligodendroglioma: Long-term results of RTOG 9402. J. Clin. Oncol. 2013, 31, 337–343. [Google Scholar] [CrossRef]
  9. Van den Bent, M.J.; Brandes, A.A.; Taphoorn, M.J.; Kros, J.M.; Kouwenhoven, M.C.; Delattre, J.Y.; Bernsen, H.J.; Frenay, M.; Tijssen, C.C.; Grisold, W.; et al. Adjuvant procarbazine, lomustine, and vincristine chemotherapy in newly diagnosed anaplastic oligodendroglioma: Long-term follow-up of EORTC brain tumor group study 26951. J. Clin. Oncol. 2013, 31, 344–350. [Google Scholar] [CrossRef]
  10. Platten, M.; Bunse, L.; Wick, A.; Bunse, T.; Le Cornet, L.; Harting, I.; Sahm, F.; Sanghvi, K.; Tan, C.L.; Poschke, I.; et al. A vaccine targeting mutant IDH1 in newly diagnosed glioma. Nature 2021, 592, 463–468. [Google Scholar] [CrossRef]
  11. Kaminska, B.; Czapski, B.; Guzik, R.; Król, S.K.; Gielniewski, B. Consequences of IDH1/2 Mutations in Gliomas and an Assessment of Inhibitors Targeting Mutated IDH Proteins. Molecules 2019, 24, 968. [Google Scholar] [CrossRef] [Green Version]
  12. Ellingson, B.M.; Wen, P.Y.; Cloughesy, T.F. Modified Criteria for Radiographic Response Assessment in Glioblastoma Clinical Trials. Neurotherapeutics 2017, 14, 307–320. [Google Scholar] [CrossRef] [Green Version]
  13. Pasquini, L.; Napolitano, A.; Tagliente, E.; Dellepiane, F.; Lucignani, M.; Vidiri, A.; Ranazzi, G.; Stoppacciaro, A.; Moltoni, G.; Nicolai, M.; et al. Deep Learning Can Differentiate IDH-Mutant from IDH-Wild GBM. J. Pers. Med. 2021, 11, 290. [Google Scholar] [CrossRef]
  14. Pasquini, L.; Napolitano, A.; Lucignani, M.; Tagliente, E.; Dellepiane, F.; Rossi-Espagnet, M.C.; Ritrovato, M.; Vidiri, A.; Villani, V.; Ranazzi, G.; et al. AI and High-Grade Glioma for Diagnosis and Outcome Prediction: Do All Machine Learning Models Perform Equally Well? Front. Oncol. 2021, 11, 601425. [Google Scholar] [CrossRef]
  15. Sotoudeh, H.; Shafaat, O.; Bernstock, J.D.; Brooks, M.D.; Elsayed, G.A.; Chen, J.A.; Szerip, P.; Chagoya, G.; Gessler, F.; Sotoudeh, E.; et al. Artificial Intelligence in the Management of Glioma: Era of Personalized Medicine. Front. Oncol. 2019, 9, 768. [Google Scholar] [CrossRef]
  16. Suh, C.H.; Kim, H.S.; Jung, S.C.; Choi, C.G.; Kim, S.J. 2-Hydroxyglutarate MR spectroscopy for prediction of isocitrate dehydrogenase mutant glioma: A systemic review and meta-analysis using individual patient data. Neuro-Oncology 2018, 20, 1573–1583. [Google Scholar] [CrossRef] [Green Version]
  17. Berrington, A.; Voets, N.L.; Larkin, S.J.; de Pennington, N.; McCullagh, J.; Stacey, R.; Schofield, C.J.; Jezzard, P.; Clare, S.; Cadoux-Hudson, T.; et al. A comparison of 2-hydroxyglutarate detection at 3 and 7 T with long-TE semi-LASER. NMR Biomed. 2018, 31, e3886. [Google Scholar] [CrossRef]
  18. Shen, X.; Voets, N.L.; Larkin, S.J.; de Pennington, N.; Plaha, P.; Stacey, R.; McCullagh, J.S.O.; Schofield, C.J.; Clare, S.; Jezzard, P.; et al. A Noninvasive Comparison Study between Human Gliomas with IDH1 and IDH2 Mutations by MR Spectroscopy. Metabolites 2019, 9, 35. [Google Scholar] [CrossRef] [Green Version]
  19. Askari, P.; Dimitrov, I.E.; Ganji, S.K.; Tiwari, V.; Levy, M.; Patel, T.R.; Pan, E.; Mickey, B.E.; Malloy, C.R.; Maher, E.A.; et al. Spectral fitting strategy to overcome the overlap between 2-hydroxyglutarate and lipid resonances at 2.25 ppm. Magn. Reson. Med. 2021, 86, 1818–1828. [Google Scholar] [CrossRef]
  20. Andronesi, O.C.; Kim, G.S.; Gerstner, E.; Batchelor, T.; Tzika, A.A.; Fantin, V.R.; Vander Heiden, M.G.; Sorensen, A.G. Detection of 2-hydroxyglutarate in IDH-mutated glioma patients by in vivo spectral-editing and 2D correlation magnetic resonance spectroscopy. Sci. Transl. Med. 2012, 4, 116ra114. [Google Scholar] [CrossRef] [Green Version]
  21. Verma, G.; Mohan, S.; Nasrallah, M.P.; Brem, S.; Lee, J.Y.; Chawla, S.; Wang, S.; Nagarajan, R.; Thomas, M.A.; Poptani, H. Non-invasive detection of 2-hydroxyglutarate in IDH-mutated gliomas using two-dimensional localized correlation spectroscopy (2D L-COSY) at 7 Tesla. J. Transl. Med. 2016, 14, 274. [Google Scholar] [CrossRef] [Green Version]
  22. Nagashima, H.; Tanaka, K.; Sasayama, T.; Irino, Y.; Sato, N.; Takeuchi, Y.; Kyotani, K.; Mukasa, A.; Mizukawa, K.; Sakata, J.; et al. Diagnostic value of glutamate with 2-hydroxyglutarate in magnetic resonance spectroscopy for IDH1 mutant glioma. Neuro-Oncology 2016, 18, 1559–1568. [Google Scholar] [CrossRef] [Green Version]
  23. Branzoli, F.; Pontoizeau, C.; Tchara, L.; Di Stefano, A.L.; Kamoun, A.; Deelchand, D.K.; Valabrègue, R.; Lehéricy, S.; Sanson, M.; Ottolenghi, C.; et al. Cystathionine as a marker for 1p/19q codeleted gliomas by in vivo magnetic resonance spectroscopy. Neuro-Oncology 2019, 21, 765–774. [Google Scholar] [CrossRef]
  24. Smits, M. Imaging of oligodendroglioma. Br. J. Radiol. 2016, 89, 20150857. [Google Scholar] [CrossRef]
  25. Bumes, E.; Wirtz, F.P.; Fellner, C.; Grosse, J.; Hellwig, D.; Oefner, P.J.; Häckl, M.; Linker, R.; Proescholdt, M.; Schmidt, N.O.; et al. Non-Invasive Prediction of IDH Mutation in Patients with Glioma WHO II/III/IV Based on F-18-FET PET-Guided In Vivo (1)H-Magnetic Resonance Spectroscopy and Machine Learning. Cancers 2020, 12, 3406. [Google Scholar] [CrossRef]
  26. Yan, H.; Parsons, D.W.; Jin, G.; McLendon, R.; Rasheed, B.A.; Yuan, W.; Kos, I.; Batinic-Haberle, I.; Jones, S.; Riggins, G.J.; et al. IDH1 and IDH2 mutations in gliomas. N. Engl. J. Med. 2009, 360, 765–773. [Google Scholar] [CrossRef]
  27. Wen, P.Y.; Huse, J.T. 2016 World Health Organization Classification of Central Nervous System Tumors. Continuum 2017, 23, 1531–1547. [Google Scholar] [CrossRef]
Figure 1. Overview of exclusion criteria and mutation status of 1H-MR spectra. 1H-MRS: proton magnetic resonance spectroscopy; IDH: isocitrate dehydrogenase; mut: mutated; wt: wildtype.
Figure 1. Overview of exclusion criteria and mutation status of 1H-MR spectra. 1H-MRS: proton magnetic resonance spectroscopy; IDH: isocitrate dehydrogenase; mut: mutated; wt: wildtype.
Cancers 14 02762 g001
Figure 2. Workflow of data collection and analysis in the original and validation cohort. 1H-MRS: 1H-magnetic resonance spectra; IDH: isocitrate dehydrogenase; mut: mutated; wt: wildtype; SVM: support vector machine.
Figure 2. Workflow of data collection and analysis in the original and validation cohort. 1H-MRS: 1H-magnetic resonance spectra; IDH: isocitrate dehydrogenase; mut: mutated; wt: wildtype; SVM: support vector machine.
Cancers 14 02762 g002
Figure 3. 1H-MR spectra obtained for an IDHmut patient with R132H mutation (A) and an IDHwt patient (B) of the validation cohort. (C) Averaged spectra of IDHmut group (marked in red, n = 23) and of IDHwt group (marked in green, n = 44) of the validation cohort. The region used for classification including M-Ins peak is indicated by the blue box. 1H-MR: 1H-magnetic resonance; IDH: isocitrate dehydrogenase; mut: mutated; wt: wildtype; M-Ins: myo-inositol.
Figure 3. 1H-MR spectra obtained for an IDHmut patient with R132H mutation (A) and an IDHwt patient (B) of the validation cohort. (C) Averaged spectra of IDHmut group (marked in red, n = 23) and of IDHwt group (marked in green, n = 44) of the validation cohort. The region used for classification including M-Ins peak is indicated by the blue box. 1H-MR: 1H-magnetic resonance; IDH: isocitrate dehydrogenase; mut: mutated; wt: wildtype; M-Ins: myo-inositol.
Cancers 14 02762 g003
Figure 4. Area under the ROC curve (AUC) in prediction of IDH status in the validation cohort. IDH: isocitrate dehydrogenase.
Figure 4. Area under the ROC curve (AUC) in prediction of IDH status in the validation cohort. IDH: isocitrate dehydrogenase.
Cancers 14 02762 g004
Figure 5. LOH 1p/19q status in the original and validation cohort in the IDHmut subgroup. (A) Distribution of LOH 1p/19q status in the original, validation, and total cohort. (B) Averaged spectra of IDHmut group with LOH 1p/19q (marked in red, n = 13) and of IDHmut group without LOH 1p/19q (marked in green, n = 26) in the total cohort. LOH: loss of heterozygosity; IDH: isocitrate dehydrogenase; mut: mutated.
Figure 5. LOH 1p/19q status in the original and validation cohort in the IDHmut subgroup. (A) Distribution of LOH 1p/19q status in the original, validation, and total cohort. (B) Averaged spectra of IDHmut group with LOH 1p/19q (marked in red, n = 13) and of IDHmut group without LOH 1p/19q (marked in green, n = 26) in the total cohort. LOH: loss of heterozygosity; IDH: isocitrate dehydrogenase; mut: mutated.
Cancers 14 02762 g005
Table 1. Characteristics of the validation cohort. WHO grade according to initial histology (WHO classification of brain tumors from 2000). IDH: isocitrate dehydrogenase; mut: mutated; wt: wildtype; LOH: loss of heterozygosity.
Table 1. Characteristics of the validation cohort. WHO grade according to initial histology (WHO classification of brain tumors from 2000). IDH: isocitrate dehydrogenase; mut: mutated; wt: wildtype; LOH: loss of heterozygosity.
IDH Mutation
IDHwtIDHmut
NumberIn %NumberIn %
Sex
men2659%1565%
women1841%835%
WHO
II37%1357%
III49%939%
IV3682%14%
unknown12%00%
LOH1p19q
no49%1043%
yes00%730%
unknown4091%626%
Total44 23
Table 2. Summary of classification results obtained for the original and the validation cohort. CI: confidence interval; AUC: area under the curve.
Table 2. Summary of classification results obtained for the original and the validation cohort. CI: confidence interval; AUC: area under the curve.
Original CohortValidation Cohort
Value (95% CI)Value (95% CI)
Accuracy (%)88.24 (72.55–96.70)76.12 (64.14–85.69)
Sensitivity (%)95.45 (77.16–99.88)82.61 (61.22–95.05)
Specificity (%)75.00 (42.81 94.51)72.73 (57.21–85.04)
Positive likelihood ratio3.82 (1.43–10.22)3.03 (1.80–5.08)
Negative likelihood ratio0.06 (0.01–0.42)0.24 (0.10–0.59)
Positive predictive value (%)87.50 (72.35–94.93)61.29 (48.55–72.66)
Negative predictive value (%)90.00 (56.33–98.43)88.89 (76.32–95.20)
Disease prevalence (%)64.71 (46.49–80.25)34.33 (23.15–46.94)
AUC0.830.82
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bumes, E.; Fellner, C.; Fellner, F.A.; Fleischanderl, K.; Häckl, M.; Lenz, S.; Linker, R.; Mirus, T.; Oefner, P.J.; Paar, C.; et al. Validation Study for Non-Invasive Prediction of IDH Mutation Status in Patients with Glioma Using In Vivo 1H-Magnetic Resonance Spectroscopy and Machine Learning. Cancers 2022, 14, 2762. https://doi.org/10.3390/cancers14112762

AMA Style

Bumes E, Fellner C, Fellner FA, Fleischanderl K, Häckl M, Lenz S, Linker R, Mirus T, Oefner PJ, Paar C, et al. Validation Study for Non-Invasive Prediction of IDH Mutation Status in Patients with Glioma Using In Vivo 1H-Magnetic Resonance Spectroscopy and Machine Learning. Cancers. 2022; 14(11):2762. https://doi.org/10.3390/cancers14112762

Chicago/Turabian Style

Bumes, Elisabeth, Claudia Fellner, Franz A. Fellner, Karin Fleischanderl, Martina Häckl, Stefan Lenz, Ralf Linker, Tim Mirus, Peter J. Oefner, Christian Paar, and et al. 2022. "Validation Study for Non-Invasive Prediction of IDH Mutation Status in Patients with Glioma Using In Vivo 1H-Magnetic Resonance Spectroscopy and Machine Learning" Cancers 14, no. 11: 2762. https://doi.org/10.3390/cancers14112762

APA Style

Bumes, E., Fellner, C., Fellner, F. A., Fleischanderl, K., Häckl, M., Lenz, S., Linker, R., Mirus, T., Oefner, P. J., Paar, C., Proescholdt, M. A., Riemenschneider, M. J., Rosengarth, K., Weis, S., Wendl, C., Wimmer, S., Hau, P., Gronwald, W., & Hutterer, M. (2022). Validation Study for Non-Invasive Prediction of IDH Mutation Status in Patients with Glioma Using In Vivo 1H-Magnetic Resonance Spectroscopy and Machine Learning. Cancers, 14(11), 2762. https://doi.org/10.3390/cancers14112762

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