1. Introduction
Gliomas represent the most clinically important group of primary brain tumors. Traditionally, they have been classified into WHO groups I—IV to evaluate their malignant potential by analysis of their morphological features. However, the past decades of research have led to the discovery of many molecular alterations in gliomas that have a great impact on the tumor’s malignancy and, accordingly, to the patient’s prognosis [
1]. Among such alterations, the mutation of isocitrate dehydrogenase (IDH) enzymes 1 or 2 is highly correlated with the patient’s overall survival, and the effect is present regardless of the tumor’s histopathological WHO grade [
2,
3,
4,
5]. IDH mutation also seems to play a pivotal role in the carcinogenesis of other solid tumors, such as cholangiocarcinoma, where it is also a major target for medical therapy [
6,
7,
8].
Normally, IDH enzymes catalyze the oxidative decarboxylation of isocitrate to form a-ketoglutarate (aKG) in the Krebs cycle. IDH1 and IDH2 localize differently in the cell but share the same function; hence, they are hereafter referred to collectively as IDH. The mutation of IDH confers a neomorphic enzyme activity that catalyzes the reduction of aKG into the putative oncometabolite D-2-hydroxyglutarate (D2HG) [
9]. The accumulation of D2HG is further associated with the hypermethylation of DNA and chromatin, which is thought to dysregulate cell epigenetics [
10,
11].
IDH mutation status is an important factor for surgical decision-making: patients with IDH-mutated tumors are more likely to have a good long-term prognosis, and thus favor aggressive gross total resection with more survival benefit to gain. Patients with IDH wild-type tumors have a generally poorer prognosis and, therefore, conservative resection to avoid neurological deficit is favored [
12,
13,
14]. The effect of gross total resection on survival remains also in recurrent diseases [
15,
16]. Current histopathological analysis based on frozen sections is unable to identify molecular characteristics, including IDH mutation, within the time frame of surgery [
17], thus creating an imminent need for new solutions.
We have previously shown that differential mobility spectrometry (DMS) is able to identify different brain tumors ex vivo [
18]. DMS characterizes substances based on the mobility differences of ionized particles in high-frequency electrical fields, resulting in a substance-specific dispersion spectrum, or “smell fingerprint” [
19]. The simplicity, quickness and cost-effectiveness of DMS makes it a compelling emerging technology for clinical applications [
18]. In this study, we demonstrate the rapid, preparation-free analysis of a tumor’s IDH mutation status with DMS.
2. Materials and Methods
We prospectively obtained tumor samples from 22 patients who had neurosurgical operations at Tampere University Hospital between the years 2018 and 2021, and at Helsinki University Hospital in 2020. Patient recruitment was continued until we had a sufficient number of IDH-mutated tumors, which are rarer. To make balanced classes, an equal number of IDH wild-type tumors were randomly selected for the experiment. Eventually, we had 11 IDH-mutated tumors and 11 IDH wild-type tumors. IDH-mutated tumors included 5 WHO gr. II–III astrocytomas, 3 gr. II–III oligodendrogliomas, and 3 gr. IV glioblastomas (GBM). IDH wild-type tumors included 1 gr. III astrocytoma and 10 GBMs. Diagnoses were made by an experienced neuropathologist and IDH mutation was identified with immunohistochemistry. The study was approved by the ethics review board of Pirkanmaa Hospital District, Finland. The patients gave their written consent for the study.
All samples were stored in a freezer at −70 °C. The samples were carefully cut into 88 (44 IDH-mutated and 44 IDH wild-type) smaller specimens of macroscopically equal sizes. Blood, if macroscopically visible, was carefully rinsed from the samples before the analysis. The samples were randomly placed in a plastic well plate with each well containing 0.18 mL of agar in the bottom. Each sample was incised with a custom-built, computer-controlled, 40 W, 10.6 μm CO
2 laser evaporator four times in a quadratic manner, with 1 mm gaps between the incisions. The total number of incisions was 352. The laser sampling was controlled by a graphical user interface. To provide a clean and controlled supply of carrier gas for the analyte gas, purified and humidified pressurized air was introduced to the sampling stage via a sampling nozzle. The sampling nozzle provided a protective stream of carrier gas around the sampling area and, after sample vaporization, transported the sample gas to the DMS inlet. The DMS used in the study was a commercial IonVision instrument (Olfactomics Oy, Finland). The measurement parameters for the DMS spectrum were: separation voltage (Usv), 200–1000 V with 20 increments; compensation voltage (Ucv), −2–10 V with 60 increments; separation field frequency, 1 MHz; and duty cycle of the field, 22%. With these parameters, the DMS measurement produced a total of 1200 data points and the duration of the measurement was approximately 13 s, during which 250 2 ms laser pulses were used to provide a sample stream of vaporized tissue to the DMS. A gross appearance of the setup (A–D) and examples of the dispersion spectra (G) are presented in
Figure 1.
We evaluated the accuracy of several machine learning algorithms for the detection of differences in dispersion spectra and the classification of the analyzed samples. Linear discriminant analysis (LDA) was found to be the best performing algorithm. The main idea of training an LDA algorithm is the projection of data points to a lower dimensional space so that the between-class distance of class centers is maximized, and the within-class distance of data points is minimized, defining a decision boundary between the classes that is used to classify new samples. The other algorithms tested were K-nearest neighbors (KNN), random forest (RF), decision tree (DT), support vector machines (SVM) and XGBoost (XGB).
3. Results
The data set revealed a temperature rise, which caused baseline drift during the measurement of one well plate, making the data biased. Thus, a necessary preprocessing method was to remove the dimension-wise linear trend which belonged the well plate from each part of the data set. This preprocessing step improved the classification results compared to the classification of the raw data. The data set contained 352 samples taken from 22 patients. Group cross-validation was utilised to estimate the classification performance. Group cross-validation is implemented so that, at every iteration, it leaves one group of samples only for testing. The other groups are used for training. In this case, the nested group cross-validation technique was used. This algorithm leaves one group for testing and the other groups are used for training and validating. For the next iteration, the second group is used for testing and the others for training and validating, and so on. This approach ensures that there are no data leakages into the training phase. With the nested group cross-validation training, the LDA algorithm reached a classification accuracy of 86%, with 86% sensitivity and 85% specificity (
Table 1). The workflow of the LDA algorithm is presented in
Figure 1F. Further details of the cross-validation and classification results reached with other algorithms are presented in the
Supplementary File.
In terms of the samples, out of the original 22 tumor samples (352 incisions), 8 samples had all their incisions correctly classified. In five samples, less than 10% of incisions were erroneous. In four samples, 10–20% were wrong. In five samples, 20–50% of the incisions were incorrectly classified. The tumors that had incorrectly clustered incisions included eight IDH wild-type tumors and six IDH-mutated tumors. The most difficult tumor type for the classifier was gr. IV GBM.
4. Discussion
Our results show that the smoke generated from the IDH-mutated and IDH wild-type gliomas had distinct DMS profiles, and the DMS could differentiate them with good sensitivity and specificity. The laser evaporator platform is compact enough to be placed in the operating room and used for intermittent analysis of the tumor samples during surgery. The duration of measurement was approximately 13 s, so the DMS operates in almost real time. The DMS is also simpler and more economical than conventional mass spectrometer-based solutions. Conventional frozen section analysis is unable to identify molecular alterations in tumors, such as IDH mutation. In the latest WHO tumor classification, these alterations have become ever more prominent. This creates an increasing need for novel tumor identification methods in neurosurgical departments worldwide.
Recently, Raman spectroscopy has also been used for genotyping unprocessed glioma samples [
20]. Raman spectroscopy is a modality that gives spectral tissue characteristics based on molecular signatures resulting from the inelastic scattering of incident light. Our results equal those achieved with Raman spectroscopy, and the workflow in DMS is at least as fast and straightforward.
Our tumor sample set included both IDH-mutated and IDH wild-type gr. IV GBMs and gr. III malignant astrocytomas. Out of the tumors with an unusual IDH mutation status given their histology, one GBM had 25% (9 out of 36) of the incisions erroneously classified, but all the other tumors (two IDH-mutated gr. IV GBMs and one IDH wild-type gr. III astrocytoma) had all their incisions correct classified, even though the opposite cluster had multiple histologically similar tumors. This indirectly indicates that the divisive features in the classification process were actually due to the cellular metabolic changes driven by an IDH mutation. The phospholipid content of tissue has previously been identified as a key distinguishing factor in DMS analysis [
18]. The metabolic changes associated with an IDH mutation include aberrations in phospholipid composition [
10], which constitutes a plausible theoretical basis for the detection of IDH mutation by DMS.
A potential source of error in DMS analysis is intratumoral heterogeneity. This is especially true in GBMs, which vary in terms of cellular density, nuclear pleomorphism, necrosis, histologic architecture, vasculature, mitoses, and multifaceted microenvironments [
21,
22]. This can cause variance in tissue impedance and disturb the classifier [
23]. An additional confounding factor in our study was 5-ALA, which was used only in the resections of tumors that radiologically appeared as malignant. However, all three IDH-mutated GBMs were resected with 5-ALA guidance, and still the classifier was able to classify them correctly.
Our study was limited by a relatively small number of samples that we multiplied into smaller specimens. In order to achieve a setup resembling actual intraoperative use, we only minimally prepared the tumor samples for the analysis. This inevitably caused spatial variance in the specimens that affected the DMS signal strength, thus creating an additional confounding factor to the classifier. This issue could be addressed in future studies by processing the samples into a more homogeneous cell suspension by a centrifuge before the analysis. The suspension could then be pipetted into the well plate to obtain precisely equal sample sizes. We also used frozen samples instead of fresh tumors. In our earlier unpublished experiments, freezing of the samples was not found to affect the classification results. However, this should be verified in peer-reviewed studies in the future.
5. Conclusions
Our results show that the DMS is able to differentiate IDH-mutated and IDH wild-type tumors with good accuracy in a setting suitable for intraoperative use. The role of molecular alterations in classifying brain tumors and evaluating their prognosis is increasing. Additionally, the degree of survival benefit achieved with a gross-total resection varies even in histologically similar tumors based on their IDH mutation status, which is impossible to identify with conventional frozen section analysis. This makes the DMS a promising novel tool for neurosurgical practice.
Author Contributions
Conceptualization, I.H., J.H., A.R. (Antti Roine), A.V. and N.O.; data curation, A.R. (Anton Rauhameri), A.K., M.K. and M.M.; formal analysis, A.R. (Anton Rauhameri), A.K., M.K. and M.M.; funding acquisition, J.H., K.N., A.R. (Antti Roine) and N.O.; investigation, I.H., A.K. and M.M.; methodology, I.H., J.H., A.K., M.K., A.R. (Antti Roine), N.O. and A.V.; project administration, J.H., K.N., A.R. (Antti Roine), A.V. and N.O.; resources, H.H., A.R. (Antti Roine), A.V., N.O., A.L. and P.K.-P.; software, A.K., M.K., A.R. (Antti Roine) and N.O.; supervision, J.H., K.N., A.R. (Antti Roine), A.V., N.O., A.L. and P.K.-P.; validation, A.V., A.R. (Antti Roine) and N.O.; visualization, I.H., A.R. (Anton Rauhameri), A.K. and M.M.; writing—original draft, I.H. and A.K.; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.
Funding
Ilkka Haapala declares funding from The Finnish Medical Foundation, Eka -grant, decision number 3535. Anton Kondratev declares funding from Academy of Finland, Programmable Scent Environments -project, decision number 323530. Anton Kontunen declares funding from the Doctoral School of Tampere University and Emil Aaltonen Foundation (Grant number 210073). Markus Karjalainen declares funding from The Finnish Cultural Foundation, Pirkanmaa Regional Fund. Kristiina Nordfors declares funding from Finnish Pediatric Research Foundation, Pediatric Cancer Research Foundation Väre and The Finnish Ministry of Social Affairs and Health. Hannu Haapasalo declares funding from Competitive State Research Financing of the Expert Responsibility area of Tampere University Hospital. Niku Oksala declares funding from Competitive State Research Financing of the Expert Responsibility area of Tampere University Hospital (Grant numbers 9s045, 9T044, 9U042, 150618, 9V044, 9X040, 9AA057, 9AB052, and MK301); from Competitive funding to strengthen university research profiles funded by Academy of Finland, decision number 292477; and from Tampere Tuberculosis Foundation. Joonas Haapasalo declares funding from Emil Aaltonen Foundation, Finnish Pediatric Research Foundation, Competitive State Research Financing of the Expert Responsibility area of Tampere University Hospital, Pediatric Cancer Research Foundation Väre and The Finnish Ministry of Social Affairs and Health. The study sponsors did not have any involvement in the study design; collection, analysis, and interpretation of data; the writing of the manuscript; or the decision to submit the manuscript for publication.
Institutional Review Board Statement
The study was approved by the Ethics Committee of Tampere University Hospital (ETL R10066) but was not listed in ClinicalTrials Database, because it was a non-randomized, non-interventional study.
Informed Consent Statement
Written informed consent has been obtained from the patients.
Data Availability Statement
The data presented in this study are available on request from the corresponding author.
Conflicts of Interest
A.R., N.O., A.Kontunen and M.K. are shareholders and employees of Olfactomics Ltd. M.M. is an employee of Olfactomics. Other authors do not declare any conflict of interest with respect to this work.
Correction Statement
This article has been republished with a minor correction to the Author Contributions. This change does not affect the scientific content of the article.
References
- Eckel-Passow, J.E.; Lachance, D.H.; Molinaro, A.M.; Walsh, K.M.; Decker, P.A.; Sicotte, H.; Pekmezci, M.; Rice, T.W.; Kosel, M.L.; Smirnov, I.V.; et al. Glioma Groups Based on 1p/19q, IDH, and TERT Promoter Mutations in Tumors. N. Engl. J. Med. 2015, 372, 2499–2508. [Google Scholar] [CrossRef] [PubMed]
- Houillier, C.; Wang, X.; Kaloshi, G.; Mokhtari, K.; Guillevin, R.; Laffaire, J.; Paris, S.; Boisselier, B.; Idbaih, A.; Laigle-Donadey, F.; et al. IDH1 or IDH2 mutations predict longer survival and response to temozolomide in low-grade gliomas. Neurology 2010, 75, 1560–1566. Available online: https://n.neurology.org/content/75/17/1560.short (accessed on 7 April 2021). [CrossRef] [PubMed]
- Metellus, P.; Coulibaly, B.; Colin, C.; de Paula, A.M.; Vasiljevic, A.; Taieb, D.; Barlier, A.; Boisselier, B.; Mokhtari, K.; Wang, X.W.; et al. Absence of IDH mutation identifies a novel radiologic and molecular subtype of WHO grade II gliomas with dismal prognosis. Acta Neuropathol. 2010, 120, 719–729. [Google Scholar] [CrossRef] [PubMed]
- Songtao, Q.; Lei, Y.; Si, G.; Yanqing, D.; Huixia, H.; Xuelin, Z.; Lanxiao, W.; Fei, Y. IDH mutations predict longer survival and response to temozolomide in secondary glioblastoma. Wiley Online Libr. 2012, 103, 269–273. [Google Scholar] [CrossRef] [PubMed]
- Juratli, T.A.; Kirsch, M.; Geiger, K.; Klink, B.; Leipnitz, E.; Pinzer, T.; Soucek, S.; Schrok, E.; Schackert, G.; Krex, D. The prognostic value of IDH mutations and MGMT promoter status in secondary high-grade gliomas. J. Neurooncol. 2012, 110, 325–333. [Google Scholar] [CrossRef] [PubMed]
- Rizzo, A.; Ricci, A.D.; Tober, N.; Nigro, M.C.; Mosca, M.; Palloni, A.; Abbati, F.; Frega, G.; De Lorenzo, S.; Tavolari, S.; et al. Second-line Treatment in Advanced Biliary Tract Cancer: Today and Tomorrow. Anticancer Res. 2020, 40, 3013–3030. [Google Scholar] [CrossRef] [PubMed]
- Rizzo, A.; Ricci, A.D.; Brandi, G. IDH inhibitors in advanced cholangiocarcinoma: Another arrow in the quiver? Cancer Treat. Res. Commun. 2021, 27, 100356. [Google Scholar] [CrossRef]
- Rizzo, A.; Brandi, G. First-line Chemotherapy in Advanced Biliary Tract Cancer Ten Years after the ABC-02 Trial: “And Yet It Moves!”. Cancer Treat. Res. Commun. 2021, 27, 100335. [Google Scholar] [CrossRef]
- Waitkus, M.S.; Diplas, B.H.; Yan, H. Isocitrate dehydrogenase mutations in gliomas. Neuro-Oncology 2015, 18, 16–26. Available online: https://academic.oup.com/neuro-oncology/article-abstract/18/1/16/2509155 (accessed on 7 April 2021). [CrossRef]
- Fack, F.; Tardito, S.; Hochart, G.; Oudin, A.; Zheng, L.; Fritah, S.; Golebiewska, A.; Nazarov, P.; Bernard, A.; Hau, A.; et al. Altered metabolic landscape in IDH-mutant gliomas affects phospholipid, energy, and oxidative stress pathways. EMBO Mol. Med. 2017, 9, 1681–1695. [Google Scholar] [CrossRef]
- Garrett, M.; Sperry, J.; Braas, D.; Yan, W.; Le, T.M.; Mottahedeh, J.; Ludwig, K.; Eskin, A.; Qin, Y.; Levy, R.; et al. Metabolic characterization of isocitrate dehydrogenase (IDH) mutant and IDH wildtype gliomaspheres uncovers cell type-specific vulnerabilities. Cancer Metab. 2018, 6, 4. [Google Scholar] [CrossRef] [PubMed]
- Beiko, J.; Suki, D.; Hess, K.R.; Fox, B.D.; Cheung, V.; Cabral, M.; Shonka, N.; Gilbert, M.R.; Sawaya, R.; Prabhu, S.S.; et al. IDH1 mutant malignant astrocytomas are more amenable to surgical resection and have a survival benefit associated with maximal surgical resection. Neuro-Oncology 2014, 16, 81–91. [Google Scholar] [CrossRef] [PubMed]
- Kawaguchi, T.; Sonoda, Y.; Shibahara, I.; Saito, R.; Kanamori, M.; Kumabe, T.; Tominaga, T. Impact of gross total resection in patients with WHO grade III glioma harboring the IDH 1/2 mutation without the 1p/19q co-deletion. J. Neurooncol. 2016, 129, 505–514. [Google Scholar] [CrossRef] [PubMed]
- Patel, T.; Bander, E.D.; Venn, R.A.; Powell, T.; Cederquist, G.Y.-M.; Schaefer, P.M.; Puchi, L.A.; Akhmerov, A.; Ogilvie, S.; Reiner, A.S.; et al. The role of extent of resection in IDH1 wild-type or mutant low-grade gliomas. Neurosurgery 2018, 82, 808–814. [Google Scholar] [CrossRef]
- Montemurro, N.; Fanelli, G.N.; Scatena, C.; Ortenzi, V.; Pasqualetti, F.; Mazzanti, C.M.; Morganti, R.; Paiar, F.; Naccarato, A.G.; Perrini, P. Surgical outcome and molecular pattern characterization of recurrent glioblastoma multiforme: A single-center retrospective series. Clin. Neurol Neurosurg. 2021, 207, 106735. [Google Scholar] [CrossRef]
- Van Tellingen, O.; Yetkin-Arik, B.; De Gooijer, M.C.; Wesseling, P.; Wurdinger, T.; De Vries, H.E. Overcoming the blood–brain tumor barrier for effective glioblastoma treatment. Drug Resist. Updat. 2015, 19, 1–12. [Google Scholar] [CrossRef]
- Longuespée, R.; Wefers, A.K.; De Vita, E.; Miller, A.K.; Reuss, D.E.; Wick, W.; Herold-Mende, C.; Kriegsmann, M.; Schirmacher, P.; Von Deimling, A.; et al. Rapid detection of 2-hydroxyglutarate in frozen sections of IDH mutant tumors by MALDI-TOF mass spectrometry. Acta Neuropathol. Commun. 2018, 6, 21. [Google Scholar] [CrossRef]
- Haapala, I.; Karjalainen, M.; Kontunen, A.S.; Vehkaoja, A.; Nordfors, K.; Haapasalo, H.; Haapasalo, J.; Oksala, N.; Roine, A. Identifying brain tumors by differential mobility spectrometry analysis of diathermy smoke. J. Neurosurg. 2020, 133, 100–106. [Google Scholar] [CrossRef]
- Schneider, B.B.; Nazarov, E.G.; Londry, F.; Vouros, P.; Covey, T.R. Differential mobility spectrometry/mass spectrometry history, theory, design optimization, simulations, and applications: Differential Mobility Spectrometry/Mass Spectrometry Non-radioactive Ion Source for Operation in Ambient Pressure View project Properties of Gas Phase Molecular Clusters View project Differential Mobility Spectrometry/Mass Spectrometry History, Theory, Design Optimization, Simulations, and Applications. Wiley Online Libr. 2016, 35, 687–737. [Google Scholar] [CrossRef]
- Sciortino, T.; Secoli, R.; D’Amico, E.; Moccia, S.; Nibali, M.C.; Gay, L.; Rossi, M.; Pecco, N.; Castellano, A.; De Momi, E.; et al. Raman Spectroscopy and Machine Learning for IDH Genotyping of Unprocessed Glioma Biopsies. Cancers 2021, 13, 4196. [Google Scholar] [CrossRef]
- Decordova, S.; Shastri, A.; Tsolaki, A.G.; Yasmin, H.; Klein, L.; Singh, S.K.; Kishore, U. Molecular Heterogeneity and Immunosuppressive Microenvironment in Glioblastoma. Front. Immunol. 2020, 11, 1402. [Google Scholar] [CrossRef] [PubMed]
- Campos, B.; Olsen, L.R.; Urup, T.; Poulsen, H.S. A comprehensive profile of recurrent glioblastoma. Oncogene 2016, 35, 5819–5825. [Google Scholar] [CrossRef] [PubMed]
- Zhang, Y.; Xu, S.; Min, W.; Shen, L.; Zhang, Y.; Yue, Z. SURG-25. A novel bio-impedance spectroscopy system real-time intraoperatively discriminates glioblastoma from brain tissue in mice. Neuro-Oncology 2017, 19 (Suppl. 6), vi240. [Google Scholar] [CrossRef]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |