Search

This article investigates methods for the accurate and robust differentiation of metastases from glioblastomas on the basis of single-voxel 1H MRS information. Single-voxel 1H MR spectra from a total of 109 patients (78 glioblastomas and 31 metastases) from the multicenter, international INTERPRET database, plus a test set of 40 patients (30 glioblastomas and 10 metastases) from three different centers in the Barcelona (Spain) metropolitan area, were analyzed using a robust method for feature (spectral frequency) selection coupled with a linear-in-the-parameters single-layer perceptron classifier. For the test set, a parsimonious selection of five frequencies yielded an area under the receiver operating characteristic curve of 0.86, and an area under the convex hull of the receiver operating characteristic curve of 0.91. Moreover, these accurate results for the discrimination between glioblastomas and metastases were obtained using a small number of frequencies that are amenable to metabolic interpretation, which should ease their use as diagnostic markers. Importantly, the prediction can be expressed as a simple formula based on a linear combination of these frequencies. As a result, new cases could be straightforwardly predicted by integrating this formula into a computer-based medical decision support system. This work also shows that the combination of spectra acquired at different TEs (short TE, 20-32 ms; long TE, 135-144 ms) is key to the successful discrimination between glioblastomas and metastases from single-voxel 1H MRS. Copyright © 2011 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

Abstract Background In-vivo single voxel proton magnetic resonance spectroscopy (SV 1H-MRS), coupled with supervised pattern recognition (PR) methods, has been widely used in clinical studies of discrimination of brain tumour types and follow-up of patients bearing abnormal brain masses. SV 1H-MRS provides useful biochemical information about the metabolic state of tumours and can be performed at short (< 45 ms) or long (> 45 ms) echo time (TE), each with particular advantages. Short-TE spectra are more adequate for detecting lipids, while the long-TE provides a much flatter signal baseline in between peaks but also negative signals for metabolites such as lactate. Both, lipids and lactate, are respectively indicative of specific metabolic processes taking place. Ideally, the information provided by both TE should be of use for clinical purposes. In this study, we characterise the performance of a range of Non-negative Matrix Factorisation (NMF) methods in two respects: first, to derive sources correlated with the mean spectra of known tissue types (tumours and normal tissue); second, taking the best performing NMF method for source separation, we compare its accuracy for class assignment when using the mixing matrix directly as a basis for classification, as against using the method for dimensionality reduction (DR). For this, we used SV 1H-MRS data with positive and negative peaks, from a widely tested SV 1H-MRS human brain tumour database. Results The results reported in this paper reveal the advantage of using a recently described variant of NMF, namely Convex-NMF, as an unsupervised method of source extraction from SV1H-MRS. Most of the sources extracted in our experiments closely correspond to the mean spectra of some of the analysed tumour types. This similarity allows accurate diagnostic predictions to be made both in fully unsupervised mode and using Convex-NMF as a DR step previous to standard supervised classification. The obtained results are comparable to, or more accurate than those obtained with supervised techniques. Conclusions The unsupervised properties of Convex-NMF place this approach one step ahead of classical label-requiring supervised methods for the discrimination of brain tumour types, as it accounts for their increasingly recognised molecular subtype heterogeneity. The application of Convex-NMF in computer assisted decision support systems is expected to facilitate further improvements in the uptake of MRS-derived information by clinicians.

Abstract Background Proton Magnetic Resonance (MR) Spectroscopy (MRS) is a widely available technique for those clinical centres equipped with MR scanners. Unlike the rest of MR-based techniques, MRS yields not images but spectra of metabolites in the tissues. In pathological situations, the MRS profile changes and this has been particularly described for brain tumours. However, radiologists are frequently not familiar to the interpretation of MRS data and for this reason, the usefulness of decision-support systems (DSS) in MRS data analysis has been explored. Results This work presents the INTERPRET DSS version 3.0, analysing the improvements made from its first release in 2002. Version 3.0 is aimed to be a program that 1st, can be easily used with any new case from any MR scanner manufacturer and 2nd, improves the initial analysis capabilities of the first version. The main improvements are an embedded database, user accounts, more diagnostic discrimination capabilities and the possibility to analyse data acquired under additional data acquisition conditions. Other improvements include a customisable graphical user interface (GUI). Most diagnostic problems included have been addressed through a pattern-recognition based approach, in which classifiers based on linear discriminant analysis (LDA) were trained and tested. Conclusions The INTERPRET DSS 3.0 allows radiologists, medical physicists, biochemists or, generally speaking, any person with a minimum knowledge of what an MR spectrum is, to enter their own SV raw data, acquired at 1.5 T, and to analyse them. The system is expected to help in the categorisation of MR Spectra from abnormal brain masses.

Abstract Background SpectraClassifier (SC) is a Java solution for designing and implementing Magnetic Resonance Spectroscopy (MRS)-based classifiers. The main goal of SC is to allow users with minimum background knowledge of multivariate statistics to perform a fully automated pattern recognition analysis. SC incorporates feature selection (greedy stepwise approach, either forward or backward), and feature extraction (PCA). Fisher Linear Discriminant Analysis is the method of choice for classification. Classifier evaluation is performed through various methods: display of the confusion matrix of the training and testing datasets; K-fold cross-validation, leave-one-out and bootstrapping as well as Receiver Operating Characteristic (ROC) curves. Results SC is composed of the following modules: Classifier design, Data exploration, Data visualisation, Classifier evaluation, Reports, and Classifier history. It is able to read low resolution in-vivo MRS (single-voxel and multi-voxel) and high resolution tissue MRS (HRMAS), processed with existing tools (jMRUI, INTERPRET, 3DiCSI or TopSpin). In addition, to facilitate exchanging data between applications, a standard format capable of storing all the information needed for a dataset was developed. Each functionality of SC has been specifically validated with real data with the purpose of bug-testing and methods validation. Data from the INTERPRET project was used. Conclusions SC is a user-friendly software designed to fulfil the needs of potential users in the MRS community. It accepts all kinds of pre-processed MRS data types and classifies them semi-automatically, allowing spectroscopists to concentrate on interpretation of results with the use of its visualisation tools.

The standard treatment in glioblastoma includes maximal safe resection followed by concomitant radiotherapy plus chemotherapy and adjuvant temozolomide. The first follow-up study to evaluate treatment response is performed 1 month after concomitant treatment, when contrast-enhancing regions may appear that can correspond to true progression or pseudoprogression. We retrospectively evaluated 31 consecutive patients at the first follow-up after concomitant treatment to check whether the metabolic pattern assessed with multivoxel MRS was predictive of treatment response 2 months later. We extracted the underlying metabolic patterns of the contrast-enhancing regions with a blind-source separation method and mapped them over the reference images. Pattern heterogeneity was calculated using entropy, and association between patterns and outcomes was measured with Cramér's V. We identified three distinct metabolic patterns-proliferative, necrotic, and responsive, which were associated with status 2 months later. Individually, 70% of the patients showed metabolically heterogeneous patterns in the contrast-enhancing regions. Metabolic heterogeneity was not related to the regions' size and only stable patients were less heterogeneous than the rest. Contrast-enhancing regions are also metabolically heterogeneous 1 month after concomitant treatment. This could explain the reported difficulty in finding robust pseudoprogression biomarkers., (© 2024 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.)

Machine Learning is entering a phase of maturity, but its medical applications still lag behind in terms of practical use. The field of oncological radiology (and neuro-oncology in particular) is at the forefront of these developments, now boosted by the success of Deep-Learning methods for the analysis of medical images. This paper reviews in detail some of the most recent advances in the use of Deep Learning in this field, from the broader topic of the development of Machine-Learning-based analytical pipelines to specific instantiations of the use of Deep Learning in neuro-oncology; the latter including its use in the groundbreaking field of ultra-low field magnetic resonance imaging.

Magnetic resonance spectroscopy (MRS) is an MR technique that provides information about the biochemistry of tissues in a noninvasive way. MRS has been widely used for the study of brain tumors, both preoperatively and during follow-up. In this study, we investigated the performance of a range of variants of unsupervised matrix factorization methods of the non-negative matrix underapproximation (NMU) family, namely, sparse NMU, global NMU, and recursive NMU, and compared them with convex non-negative matrix factorization (C-NMF), which has previously shown a good performance on brain tumor diagnostic support problems using MRS data. The purpose of the investigation was 2-fold: first, to ascertain the differences among the sources extracted by these methods; and second, to compare the influence of each method in the diagnostic accuracy of the classification of brain tumors, using them as feature extractors. We discovered that, first, NMU variants found meaningful sources in terms of biological interpretability, but representing parts of the spectrum, in contrast to C-NMF; and second, that NMU methods achieved better classification accuracy than C-NMF for the classification tasks when one class was not meningioma., (© 2023 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.)

In vivo magnetic resonance spectroscopy (MRS) has two modalities, single-voxel (SV) and multivoxel (MV), in which one or more contiguous grids of SVs are acquired., Purpose: To test whether MV grids can be classified with models trained with SV., Methods: Retrospective study. Training dataset: Multicenter multiformat SV INTERPRET, 1.5T. Testing dataset: MV eTumour, 3T. Two classification tasks were completed: 3-class (meningioma vs. aggressive vs. normal) and 4-class (meningioma vs. low-grade glioma vs. aggressive vs. normal). Five different methods were tested for feature selection. The classification was implemented using linear discriminant analysis (LDA), random forest, and support vector machines. The evaluation was completed with balanced error rate (BER) and area under the curve (AUC) on both sets. The accuracy in class prediction was calculated by developing a solid tumor index (STI) and segmentation accuracy with the Dice score., Results: The best method was sequential forward feature selection combined with LDA, with AUCs = 0.95 (meningioma), 0.89 (aggressive), 0.82 (low-grade glioma), and 0.82 (normal). STI was 66% (4-class task) and 71% (3-class task) because two cases failed completely and two more had suboptimal STI as defined by us., Discussion: The reasons for failure in the classification of the MV test set were related to the presence of artifacts.

Despite the success of automated pattern recognition methods in problems of human brain tumor diagnostic classification, limited attention has been paid to the issue of automated data quality assessment in the field of MRS for neuro-oncology. Beyond some early attempts to address this issue, the current standard in practice is MRS quality control through human (expert-based) assessment. One aspect of automatic quality control is the problem of detecting artefacts in MRS data. Artefacts, whose variety has already been reviewed in some detail and some of which may even escape human quality control, have a negative influence in pattern recognition methods attempting to assist tumor characterization. The automatic detection of MRS artefacts should be beneficial for radiology as it guarantees more reliable tumor characterizations, as well as the development of more robust pattern recognition-based tumor classifiers and more trustable MRS data processing and analysis pipelines. Feature extraction methods have previously been used to help distinguishing between good and bad quality spectra to apply subsequent supervised pattern recognition techniques. In this study, we apply feature extraction differently and use a variant of a method for blind source separation, namely Convex Non-Negative Matrix Factorization, to unveil MRS signal sources in a completely unsupervised way. We hypothesize that, while most sources will correspond to the different tumor patterns, some of them will reflect signal artefacts. The experimental work reported in this paper, analyzing a combined short and long echo time 1 H-MRS database of more than 2000 spectra acquired at 1.5T and corresponding to different tumor types and other anomalous masses, provides a first proof of concept that points to the possible validity of this approach., (© 2019 John Wiley & Sons, Ltd.)

Glioblastoma is the most frequent aggressive primary brain tumor amongst human adults. Its standard treatment involves chemotherapy, for which the drug temozolomide is a common choice. These are heterogeneous and variable tumors which might benefit from personalized, data-based therapy strategies, and for which there is room for improvement in therapy response follow-up, investigated with preclinical models. This study addresses a preclinical question that involves distinguishing between treated and control (untreated) mice bearing glioblastoma, using machine learning techniques, from magnetic resonance-based data in two modalities: MRI and MRSI. It aims to go beyond the comparison of methods for such discrimination to provide an analytical pipeline that could be used in subsequent human studies. This analytical pipeline is meant to be a usable and interpretable tool for the radiology expert in the hope that such interpretation helps revealing new insights about the problem itself. For that, we propose coupling source extraction-based and radiomics-based data transformations with feature selection. Special attention is paid to the generation of radiologist-friendly visual nosological representations of the analyzed tumors.

The aim of the paper is two-fold. First, we show that structure finding with the PC algorithm can be inherently unstable and requires further operational constraints in order to consistently obtain models that are faithful to the data. We propose a methodology to stabilise the structure finding process, minimising both false positive and false negative error rates. This is demonstrated with synthetic data. Second, to apply the proposed structure finding methodology to a data set comprising single-voxel Magnetic Resonance Spectra of normal brain and three classes of brain tumours, to elucidate the associations between brain tumour types and a range of observed metabolites that are known to be relevant for their characterisation. The data set is bootstrapped in order to maximise the robustness of feature selection for nominated target variables. Specifically, Conditional Independence maps (CI-maps) built from the data and their derived Bayesian networks have been used. A Directed Acyclic Graph (DAG) is built from CI-maps, being a major challenge the minimization of errors in the graph structure. This work presents empirical evidence on how to reduce false positive errors via the False Discovery Rate, and how to identify appropriate parameter settings to improve the False Negative Reduction. In addition, several node ordering policies are investigated that transform the graph into a DAG. The obtained results show that ordering nodes by strength of mutual information can recover a representative DAG in a reasonable time, although a more accurate graph can be recovered using a random order of samples at the expense of increasing the computation time., Competing Interests: The authors have declared that no competing interests exist.

Glioblastomas (GB) are brain tumours with poor prognosis even after aggressive therapy. Improvements in both therapeutic and follow-up strategies are urgently needed. In previous work we described an oscillatory pattern of response to Temozolomide (TMZ) using a standard administration protocol, detected through MRSI-based machine learning approaches. In the present work, we have introduced the Immune-Enhancing Metronomic Schedule (IMS) with an every 6-d TMZ administration at 60 mg/kg and investigated the consistence of such oscillatory behaviour. A total of n = 17 GL261 GB tumour-bearing C57BL/6j mice were studied with MRI/MRSI every 2 d, and the oscillatory behaviour (6.2 ± 1.5 d period from the TMZ administration day) was confirmed during response. Furthermore, IMS-TMZ produced significant improvement in mice survival (22.5 ± 3.0 d for controls vs 135.8 ± 78.2 for TMZ-treated), outperforming standard TMZ treatment. Histopathological correlation was investigated in selected tumour samples (n = 6) analyzing control and responding fields. Significant differences were found for CD3+ cells (lymphocytes, 3.3 ± 2.5 vs 4.8 ± 2.9, respectively) and Iba-1 immunostained area (microglia/macrophages, 16.8% ± 9.7% and 21.9% ± 11.4%, respectively). Unexpectedly, during IMS-TMZ treatment, tumours from some mice (n = 6) fully regressed and remained undetectable without further treatment for 1 mo. These animals were considered "cured" and a GL261 re-challenge experiment performed, with no tumour reappearance in five out of six cases. Heterogeneous therapy response outcomes were detected in tumour-bearing mice, and a selected group was investigated (n = 3 non-responders, n = 6 relapsing tumours, n = 3 controls). PD-L1 content was found ca. 3-fold increased in the relapsing group when comparing with control and non-responding groups, suggesting that increased lymphocyte inhibition could be associated to IMS-TMZ failure. Overall, data suggest that host immune response has a relevant role in therapy response/escape in GL261 tumours under IMS-TMZ therapy. This is associated to changes in the metabolomics pattern, oscillating every 6 d, in agreement with immune cycle length, which is being sampled by MRSI-derived nosological images., (© 2020 John Wiley & Sons, Ltd.)

The contribution of MRS(I) to the in vivo evaluation of cancer-metabolism-derived metrics, mostly since 2016, is reviewed here. Increased carbon consumption by tumour cells, which are highly glycolytic, is now being sampled by 13 C magnetic resonance spectroscopic imaging (MRSI) following the injection of hyperpolarized [1- 13 C] pyruvate (Pyr). Hot-spots of, mostly, increased lactate dehydrogenase activity or flow between Pyr and lactate (Lac) have been seen with cancer progression in prostate (preclinical and in humans), brain and pancreas (both preclinical) tumours. Therapy response is usually signalled by decreased Lac/Pyr 13 C-labelled ratio with respect to untreated or non-responding tumour. For therapeutic agents inducing tumour hypoxia, the 13 C-labelled Lac/bicarbonate ratio may be a better metric than the Lac/Pyr ratio. 31 P MRSI may sample intracellular pH changes from brain tumours (acidification upon antiangiogenic treatment, basification at fast proliferation and relapse). The steady state tumour metabolome pattern is still in use for cancer evaluation. Metrics used for this range from quantification of single oncometabolites (such as 2-hydroxyglutarate in mutant IDH1 glial brain tumours) to selected metabolite ratios (such as total choline to N-acetylaspartate (plain ratio or CNI index)) or the whole 1 H MRSI(I) pattern through pattern recognition analysis. These approaches have been applied to address different questions such as tumour subtype definition, following/predicting the response to therapy or defining better resection or radiosurgery limits., (© 2019 John Wiley & Sons, Ltd.)

Purpose: To investigate and compare human judgment and machine learning tools for quality assessment of clinical MR spectra of brain tumors., Methods: A very large set of 2574 single voxel spectra with short and long echo time from the eTUMOUR and INTERPRET databases were used for this analysis. Original human quality ratings from these studies as well as new human guidelines were used to train different machine learning algorithms for automatic quality control (AQC) based on various feature extraction methods and classification tools. The performance was compared with variance in human judgment., Results: AQC built using the RUSBoost classifier that combats imbalanced training data performed best. When furnished with a large range of spectral and derived features where the most crucial ones had been selected by the TreeBagger algorithm it showed better specificity (98%) in judging spectra from an independent test-set than previously published methods. Optimal performance was reached with a virtual three-class ranking system., Conclusion: Our results suggest that feature space should be relatively large for the case of MR tumor spectra and that three-class labels may be beneficial for AQC. The best AQC algorithm showed a performance in rejecting spectra that was comparable to that of a panel of human expert spectroscopists. Magn Reson Med 79:2500-2510, 2018. © 2017 International Society for Magnetic Resonance in Medicine., (© 2017 International Society for Magnetic Resonance in Medicine.)

Glioblastoma (GBM) causes poor survival in patients even when applying aggressive treatment. Temozolomide (TMZ) is the standard chemotherapeutic choice for GBM treatment, but resistance always ensues. In previous years, efforts have focused on new therapeutic regimens with conventional drugs to activate immune responses that may enhance tumor regression and prevent regrowth, for example the "metronomic" approaches. In metronomic scheduling studies, cyclophosphamide (CPA) in GL261 GBM growing subcutaneously in C57BL/6 mice was shown not only to activate antitumor CD8 + T-cell response, but also to induce long-term specific T-cell tumor memory. Accordingly, we have evaluated whether metronomic CPA or TMZ administration could increase survival in orthotopic GL261 in C57BL/6 mice, an immunocompetent model. Longitudinal in vivo studies with CPA (140 mg/kg) or TMZ (range 140-240 mg/kg) metronomic administration (every 6 days) were performed in tumor-bearing mice. Tumor evolution was monitored at 7 T with MRI (T 2 -weighted, diffusion-weighted imaging) and MRSI-based nosological images of response to therapy. Obtained results demonstrated that both treatments resulted in increased survival (38.6 ± 21.0 days, n = 30) compared with control (19.4 ± 2.4 days, n = 18). Best results were obtained with 140 mg/kg TMZ (treated, 44.9 ± 29.0 days, n = 12, versus control, 19.3 ± 2.3 days, n = 12), achieving a longer survival rate than previous group work using three cycles of TMZ therapy at 60 mg/kg (33.9 ± 11.7 days, n = 38). Additional interesting findings were, first, clear edema appearance during chemotherapeutic treatment, second, the ability to apply the semi-supervised source analysis previously developed in our group for non-invasive TMZ therapy response monitoring to detect CPA-induced response, and third, the necropsy findings in mice cured from GBM after high TMZ cumulative dosage (980-1400 mg/kg), which demonstrated lymphoma incidence. In summary, every 6 day administration schedule of TMZ or CPA improves survival in orthotopic GL261 GBM with respect to controls or non-metronomic therapy, in partial agreement with previous work on subcutaneous GL261., (Copyright © 2017 John Wiley & Sons, Ltd.)

Glioblastoma (GBM) is the most common aggressive primary brain tumor in adults, with a short survival time even after aggressive therapy. Non-invasive surrogate biomarkers of therapy response may be relevant for improving patient survival. Previous work produced such biomarkers in preclinical GBM using semi-supervised source extraction and single-slice Magnetic Resonance Spectroscopic Imaging (MRSI). Nevertheless, GBMs are heterogeneous and single-slice studies could prevent obtaining relevant information. The purpose of this work was to evaluate whether a multi-slice MRSI approach, acquiring consecutive grids across the tumor, is feasible for preclinical models and may produce additional insight into therapy response. Nosological images were analyzed pixel-by-pixel and a relative responding volume, the Tumor Responding Index ( TRI ), was defined to quantify response. Heterogeneous response levels were observed and treated animals were ascribed to three arbitrary predefined groups: high response (HR, n = 2), TRI = 68.2 ± 2.8%, intermediate response (IR, n = 6), TRI = 41.1 ± 4.2% and low response (LR, n = 2), TRI = 13.4 ± 14.3%, producing therapy response categorization which had not been fully registered in single-slice studies. Results agreed with the multi-slice approach being feasible and producing an inverse correlation between TRI and Ki67 immunostaining. Additionally, ca. 7-day oscillations of TRI were observed, suggesting that host immune system activation in response to treatment could contribute to the responding patterns detected., Competing Interests: The authors declare no conflict of interest.

Background: Magnetic resonance spectroscopy (MRS) provides non-invasive information about the metabolic pattern of the brain parenchyma in vivo. The SpectraClassifier software performs MRS pattern-recognition by determining the spectral features (metabolites) which can be used objectively to classify spectra. Our aim was to develop an Infarct Evolution Classifier and a Brain Regions Classifier in a rat model of focal ischemic stroke using SpectraClassifier., Results: A total of 164 single-voxel proton spectra obtained with a 7 Tesla magnet at an echo time of 12 ms from non-infarcted parenchyma, subventricular zones and infarcted parenchyma were analyzed with SpectraClassifier ( http://gabrmn.uab.es/?q=sc ). The spectra corresponded to Sprague-Dawley rats (healthy rats, n = 7) and stroke rats at day 1 post-stroke (acute phase, n = 6 rats) and at days 7 ± 1 post-stroke (subacute phase, n = 14). In the Infarct Evolution Classifier, spectral features contributed by lactate + mobile lipids (1.33 ppm), total creatine (3.05 ppm) and mobile lipids (0.85 ppm) distinguished among non-infarcted parenchyma (100% sensitivity and 100% specificity), acute phase of infarct (100% sensitivity and 95% specificity) and subacute phase of infarct (78% sensitivity and 100% specificity). In the Brain Regions Classifier, spectral features contributed by myoinositol (3.62 ppm) and total creatine (3.04/3.05 ppm) distinguished among infarcted parenchyma (100% sensitivity and 98% specificity), non-infarcted parenchyma (84% sensitivity and 84% specificity) and subventricular zones (76% sensitivity and 93% specificity)., Conclusion: SpectraClassifier identified candidate biomarkers for infarct evolution (mobile lipids accumulation) and different brain regions (myoinositol content).

Background: Magnetic resonance spectroscopy provides metabolic information about living tissues in a non-invasive way. However, there are only few multi-centre clinical studies, mostly performed on a single scanner model or data format, as there is no flexible way of documenting and exchanging processed magnetic resonance spectroscopy data in digital format. This is because the DICOM standard for spectroscopy deals with unprocessed data. This paper proposes a plugin tool developed for jMRUI, namely jMRUI2XML, to tackle the latter limitation. jMRUI is a software tool for magnetic resonance spectroscopy data processing that is widely used in the magnetic resonance spectroscopy community and has evolved into a plugin platform allowing for implementation of novel features., Results: jMRUI2XML is a Java solution that facilitates common preprocessing of magnetic resonance spectroscopy data across multiple scanners. Its main characteristics are: 1) it automates magnetic resonance spectroscopy preprocessing, and 2) it can be a platform for outputting exchangeable magnetic resonance spectroscopy data. The plugin works with any kind of data that can be opened by jMRUI and outputs in extensible markup language format. Data processing templates can be generated and saved for later use. The output format opens the way for easy data sharing- due to the documentation of the preprocessing parameters and the intrinsic anonymization--for example for performing pattern recognition analysis on multicentre/multi-manufacturer magnetic resonance spectroscopy data., Conclusions: jMRUI2XML provides a self-contained and self-descriptive format accounting for the most relevant information needed for exchanging magnetic resonance spectroscopy data in digital form, as well as for automating its processing. This allows for tracking the procedures the data has undergone, which makes the proposed tool especially useful when performing pattern recognition analysis. Moreover, this work constitutes a first proposal for a minimum amount of information that should accompany any magnetic resonance processed spectrum, towards the goal of achieving better transferability of magnetic resonance spectroscopy studies.

Glioblastoma (Gb) is one of the most deadly tumors. Its molecular subtypes are yet to be fully characterized while the attendant efforts for personalized medicine need to be intensified in relation to glioblastoma diagnosis, treatment, and prognosis. Several molecular signatures based on gene expression microarrays were reported, but the use of microarrays for routine clinical practice is challenged by attendant economic costs. Several authors have proposed discriminant equations based on RT-PCR. Still, the discriminant threshold is often incompletely described, which makes proper validation difficult. In a previous work, we have reported two Gb subtypes based on the expression levels of four genes: CHI3L1, LDHA, LGALS1, and IGFBP3. One Gb subtype presented with low expression of the four genes mentioned, and of MGMT in a large portion of the patients (with anticipated high methylation of its promoter), and mutated IDH1. Here, we evaluate the robustness of the equations fitted with these genes using RT-PCR values in a set of 64 cases and importantly, define an unequivocal discriminant threshold with a view to prognostic implications. We developed two approaches to generate the discriminant equations: 1) using the expression level of the four genes mentioned above, and 2) using those genes displaying the highest correlation with survival among the aforementioned four ones, plus MGMT, as an attempt to further reduce the number of genes. The ease of equations' applicability, reduction in cost for raw data, and robustness in terms of resampling-based classification accuracy warrant further evaluation of these equations to discern Gb tumor biopsy heterogeneity at molecular level, diagnose potential malignancy, and prognosis of individual patients with glioblastomas.

Non-invasive monitoring of response to treatment of glioblastoma (GB) is nowadays carried out using MRI. MRS and MR spectroscopic imaging (MRSI) constitute promising tools for this undertaking. A temozolomide (TMZ) protocol was optimized for GL261 GB. Sixty-three mice were studied by MRI/MRS/MRSI. The spectroscopic information was used for the classification of control brain and untreated and responding GB, and validated against post-mortem immunostainings in selected animals. A classification system was developed, based on the MRSI-sampled metabolome of normal brain parenchyma, untreated and responding GB, with a 93% accuracy. Classification of an independent test set yielded a balanced error rate of 6% or less. Classifications correlated well both with tumor volume changes detected by MRI after two TMZ cycles and with the histopathological data: a significant decrease (p < 0.05) in the proliferation and mitotic rates and a 4.6-fold increase in the apoptotic rate. A surrogate response biomarker based on the linear combination of 12 spectral features has been found in the MRS/MRSI pattern of treated tumors, allowing the non-invasive classification of growing and responding GL261 GB. The methodology described can be applied to preclinical treatment efficacy studies to test new antitumoral drugs, and begets translational potential for early response detection in clinical studies., (Copyright © 2014 John Wiley & Sons, Ltd.)

In a previous study, we have shown the added value of (1) H MRS for the neuroradiological characterisation of adult human brain tumours. In that study, several methods of MRS analysis were used, and a software program, the International Network for Pattern Recognition of Tumours Using Magnetic Resonance Decision Support System 1.0 (INTERPRET DSS 1.0), with a short-TE classifier, provided the best results. Since then, the DSS evolved into a version 2.0 that contains an additional long-TE classifier. This study has two objectives. First, to determine whether clinicians with no experience of spectroscopy are comparable with spectroscopists in the use of the system, when only minimum training in the use of the system was given. Second, to assess whether or not a version with another TE is better than the initial version. We undertook a second study with the same cases and nine evaluators to assess whether the diagnostic accuracy of DSS 2.0 was comparable with the values obtained with DSS 1.0. In the second study, the analysis protocol was flexible in comparison with the first one to mimic a clinical environment. In the present study, on average, each case required 5.4 min by neuroradiologists and 9 min by spectroscopists for evaluation. Most classes and superclasses of tumours gave the same results as with DSS 1.0, except for astrocytomas of World Health Organization (WHO) grade III, in which performance measured as the area under the curve (AUC) decreased: AUC = 0.87 (0.72-1.02) with DSS 1.0 and AUC = 0.62 (0.55-0.70) with DSS 2.0. When analysing the performance of radiologists and spectroscopists with respect to DSS 1.0, the results were the same for most classes. Having data with two TEs instead of one did not affect the results of the evaluation., (Copyright © 2014 John Wiley & Sons, Ltd.)

Background: The clinical investigation of human brain tumors often starts with a non-invasive imaging study, providing information about the tumor extent and location, but little insight into the biochemistry of the analyzed tissue. Magnetic Resonance Spectroscopy can complement imaging by supplying a metabolic fingerprint of the tissue. This study analyzes single-voxel magnetic resonance spectra, which represent signal information in the frequency domain. Given that a single voxel may contain a heterogeneous mix of tissues, signal source identification is a relevant challenge for the problem of tumor type classification from the spectroscopic signal., Methodology/principal Findings: Non-negative matrix factorization techniques have recently shown their potential for the identification of meaningful sources from brain tissue spectroscopy data. In this study, we use a convex variant of these methods that is capable of handling negatively-valued data and generating sources that can be interpreted as tumor class prototypes. A novel approach to convex non-negative matrix factorization is proposed, in which prior knowledge about class information is utilized in model optimization. Class-specific information is integrated into this semi-supervised process by setting the metric of a latent variable space where the matrix factorization is carried out. The reported experimental study comprises 196 cases from different tumor types drawn from two international, multi-center databases. The results indicate that the proposed approach outperforms a purely unsupervised process by achieving near perfect correlation of the extracted sources with the mean spectra of the tumor types. It also improves tissue type classification., Conclusions/significance: We show that source extraction by unsupervised matrix factorization benefits from the integration of the available class information, so operating in a semi-supervised learning manner, for discriminative source identification and brain tumor labeling from single-voxel spectroscopy data. We are confident that the proposed methodology has wider applicability for biomedical signal processing.

This article investigates methods for the accurate and robust differentiation of metastases from glioblastomas on the basis of single-voxel (1)H MRS information. Single-voxel (1)H MR spectra from a total of 109 patients (78 glioblastomas and 31 metastases) from the multicenter, international INTERPRET database, plus a test set of 40 patients (30 glioblastomas and 10 metastases) from three different centers in the Barcelona (Spain) metropolitan area, were analyzed using a robust method for feature (spectral frequency) selection coupled with a linear-in-the-parameters single-layer perceptron classifier. For the test set, a parsimonious selection of five frequencies yielded an area under the receiver operating characteristic curve of 0.86, and an area under the convex hull of the receiver operating characteristic curve of 0.91. Moreover, these accurate results for the discrimination between glioblastomas and metastases were obtained using a small number of frequencies that are amenable to metabolic interpretation, which should ease their use as diagnostic markers. Importantly, the prediction can be expressed as a simple formula based on a linear combination of these frequencies. As a result, new cases could be straightforwardly predicted by integrating this formula into a computer-based medical decision support system. This work also shows that the combination of spectra acquired at different TEs (short TE, 20-32 ms; long TE, 135-144 ms) is key to the successful discrimination between glioblastomas and metastases from single-voxel (1)H MRS., (Copyright © 2011 John Wiley & Sons, Ltd.)

The purpose of this study was to evaluate whether single-voxel (1)H MRS could add useful information to conventional MRI in the preoperative characterisation of the type and grade of brain tumours. MRI and MRS examinations from a prospective cohort of 40 consecutive patients were analysed double blind by radiologists and spectroscopists before the histological diagnosis was known. The spectroscopists had only the MR spectra, whereas the radiologists had both the MR images and basic clinical details (age, sex and presenting symptoms). Then, the radiologists and spectroscopists exchanged their predictions and re-evaluated their initial opinions, taking into account the new evidence. Spectroscopists used four different systems of analysis for (1)H MRS data, and the efficacy of each of these methods was also evaluated. Information extracted from (1)H MRS significantly improved the radiologists' MRI-based characterisation of grade IV tumours (glioblastomas, metastases, medulloblastomas and lymphomas) in the cohort [area under the curve (AUC) in the MRI re-evaluation 0.93 versus AUC in the MRI evaluation 0.85], and also of the less malignant glial tumours (AUC in the MRI re-evaluation 0.93 versus AUC in the MRI evaluation 0.81). One of the MRS analysis systems used, the INTERPRET (International Network for Pattern Recognition of Tumours Using Magnetic Resonance) decision support system, outperformed the others, as well as being better than the MRI evaluation for the characterisation of grade III astrocytomas. Thus, preoperative MRS data improve the radiologists' performance in diagnosing grade IV tumours and, for those of grade II-III, MRS data help them to recognise the glial lineage. Even in cases in which their diagnoses were not improved, the provision of MRS data to the radiologists had no negative influence on their predictions., (Copyright © 2011 John Wiley & Sons, Ltd.)

Classifiers based on statistical pattern recognition analysis of MRSI data are becoming important tools for the non-invasive diagnosis of human brain tumors. Here we investigate the potential interest of perturbation-enhanced MRSI (PE-MRSI), in this case acute hyperglycemia, for improving the discrimination between mouse brain MRS patterns of glioblastoma multiforme (GBM), oligodendroglioma (ODG), and non-tumor brain parenchyma (NT). Six GBM-bearing mice and three ODG-bearing mice were scanned at 7 Tesla by PRESS-MRSI with 12 and 136 ms echo-time, during euglycemia (Eug) and also during induced acute hyperglycemia (Hyp), generating altogether four datasets per animal (echo time + glycemic condition): 12Eug, 136Eug, 12Hyp, and 136Hyp. For classifier development all spectral vectors (spv) selected from the MRSI matrix were unit length normalized (UL2) and used either as a training set (76 GBM spv, four mice; 70 ODG spv, two mice; 54 NT spv) or as an independent testing set (61 GBM spv, two mice; 31 ODG, one mouse; 23 NT spv). All Fisher's LDA classifiers obtained were evaluated as far as their descriptive performance-correctly classified cases of the training set (bootstrapping)-and predictive accuracy-balanced error rate of independent testing set classification. MRSI-based classifiers at 12Hyp were consistently more efficient in separating GBM, ODG, and NT regions, with overall accuracies always >80% and up to 95-96%; remaining classifiers were within the 48-85% range. This was also confirmed by user-independent selection of training and testing sets, using leave-one-out (LOO). This highlights the potential interest of perturbation-enhanced MRSI protocols for improving the non-invasive characterization of preclinical brain tumors., (This journal is © The Royal Society of Chemistry 2012)

Background: Pattern Recognition techniques can provide invaluable insights in the field of neuro-oncology. This is because the clinical analysis of brain tumors requires the use of non-invasive methods that generate complex data in electronic format. Magnetic Resonance (MR), in the modalities of spectroscopy (MRS) and spectroscopic imaging (MRSI), has been widely applied to this purpose. The heterogeneity of the tissue in the brain volumes analyzed by MR remains a challenge in terms of pathological area delimitation., Methodology/principal Findings: A pre-clinical study was carried out using seven brain tumor-bearing mice. Imaging and spectroscopy information was acquired from the brain tissue. A methodology is proposed to extract tissue type-specific sources from these signals by applying Convex Non-negative Matrix Factorization (Convex-NMF). Its suitability for the delimitation of pathological brain area from MRSI is experimentally confirmed by comparing the images obtained with its application to selected target regions, and to the gold standard of registered histopathology data. The former showed good accuracy for the solid tumor region (proliferation index (PI)>30%). The latter yielded (i) high sensitivity and specificity in most cases, (ii) acquisition conditions for safe thresholds in tumor and non-tumor regions (PI>30% for solid tumoral region; ≤5% for non-tumor), and (iii) fairly good results when borderline pixels were considered., Conclusions/significance: The unsupervised nature of Convex-NMF, which does not use prior information regarding the tumor area for its delimitation, places this approach one step ahead of classical label-requiring supervised methods for discrimination between tissue types, minimizing the negative effect of using mislabeled voxels. Convex-NMF also relaxes the non-negativity constraints on the observed data, which allows for a natural representation of the MRSI signal. This should help radiologists to accurately tackle one of the main sources of uncertainty in the clinical management of brain tumors, which is the difficulty of appropriately delimiting the pathological area.

MRI and MRS are established methodologies for evaluating intracranial lesions. One MR spectral feature suggested for in vivo grading of astrocytic tumours is the apparent myo-lnositol (ml) intensity (ca 3.55 ppm) at short echo times, although glycine (gly) may also contribute in vivo to this resonance. The purpose of this study was to quantitatively evaluate the ml + gly contribution to the recorded spectral pattern in vivo and correlate it with in vitro data obtained from perchloric acid extraction of tumour biopsies. Patient spectra (n = 95) at 1.5T at short (20-31 ms) and long (135-136 ms) echo times were obtained from the INTERPRET MRS database (http://gabrmn.uab.eslinterpretvalidateddbl). Phantom spectra were acquired with a comparable protocol. Spectra were automatically processed and the ratios of the (ml + gly) to Cr peak heights ((ml + gly)/Cr) calculated. Perchloric acid extracts of brain tumour biopsies were analysed by high-resolution NMR at 9.4T. The ratio (ml + gly)/Cr decreased significantly with astrocytic grade in vivo between low-grade astrocytoma (A2) and glioblastoma multiforme (GBM). In vitro results displayed a somewhat different tendency, with anaplastic astrocytomas having significantly higher (ml + gly)/Cr than A2 and GBM. The discrepancy between in vivo and in vitro data suggests that the NMR visibility of glycine in glial brain tumours is restricted in vivo.

In the last decade, machine learning (ML) techniques have been used for developing classifiers for automatic brain tumour diagnosis. However, the development of these ML models rely on a unique training set and learning stops once this set has been processed. Training these classifiers requires a representative amount of data, but the gathering, preprocess, and validation of samples is expensive and time-consuming. Therefore, for a classical, non-incremental approach to ML, it is necessary to wait long enough to collect all the required data. In contrast, an incremental learning approach may allow us to build an initial classifier with a smaller number of samples and update it incrementally when new data are collected. In this study, an incremental learning algorithm for Gaussian Discriminant Analysis (iGDA) based on the Graybill and Deal weighted combination of estimators is introduced. Each time a new set of data becomes available, a new estimation is carried out and a combination with a previous estimation is performed. iGDA does not require access to the previously used data and is able to include new classes that were not in the original analysis, thus allowing the customization of the models to the distribution of data at a particular clinical center. An evaluation using five benchmark databases has been used to evaluate the behaviour of the iGDA algorithm in terms of stability-plasticity, class inclusion and order effect. Finally, the iGDA algorithm has been applied to automatic brain tumour classification with magnetic resonance spectroscopy, and compared with two state-of-the-art incremental algorithms. The empirical results obtained show the ability of the algorithm to learn in an incremental fashion, improving the performance of the models when new information is available, and converging in the course of time. Furthermore, the algorithm shows a negligible instance and concept order effect, avoiding the bias that such effects could introduce., (Copyright © 2011 Elsevier Inc. All rights reserved.)

In order to evaluate the relevance of magnetic resonance (MR) features selected by automatic feature selection techniques to build classifiers for differential diagnosis and tissue segmentation two data sets containing MR spectroscopy data from patients with brain tumours were investigated. The automatically selected features were evaluated using literature and clinical experience. It was observed that a significant part of the automatically selected features correspond to what is known from the literature and clinical experience. We conclude that automatic feature selection is a useful tool to obtain relevant and possibly interesting features, but evaluation of the obtained features remains necessary., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

Object: This study demonstrates that 3T SV-MRS data can be used with the currently available automatic brain tumour diagnostic classifiers which were trained on databases of 1.5T spectra. This will allow the existing large databases of 1.5T MRS data to be used for diagnostic classification of 3T spectra, and perhaps also the combination of 1.5T and 3T databases., Materials and Methods: Brain tumour classifiers trained with 154 1.5T spectra to discriminate among high grade malignant tumours and common grade II glial tumours were evaluated with a subsequently-acquired set of 155 1.5T and 37 3T spectra. A similarity study between spectra and main brain tumour metabolite ratios for both field strengths (1.5T and 3T) was also performed., Results: Our results showed that classifiers trained with 1.5T samples had similar accuracy for both test datasets (0.87 ± 0.03 for 1.5T and 0.88 ± 0.03 for 3.0T). Moreover, non-significant differences were observed with most metabolite ratios and spectral patterns., Conclusion: These results encourage the use of existing classifiers based on 1.5T datasets for diagnosis with 3T (1)H SV-MRS. The large 1.5T databases compiled throughout many years and the prediction models based on 1.5T acquisitions can therefore continue to be used with data from the new 3T instruments.

Background and Purpose: There is a large range of survival times in patients with HGA that can only be partially explained by histologic grade and clinical aspects. This study aims to retrospectively assess the predictive value of single-voxel (1)H-MRS regarding survival in HGA., Materials and Methods: Pretreatment (1)H-MRS in 187 patients with HGA produced 180 spectra at STE (30 ms) and 182 at LTE (136 ms). Patients were dichotomized into 2 groups according to survival better or worse than the median. The spectra of the 2 groups were compared using the Mann-Whitney U test. The points on the spectrum with the most significant differences were selected for discriminating patients with good and poor prognosis. Thresholds were defined with ROC curves, and survival was analyzed by using the Kaplan-Meier method and the Cox proportional hazards model., Results: Four points on the spectrum showed the most significant differences: 0.98 and 3.67 ppm at STE; and 0.98 and 1.25 ppm at LTE (P between <.001 and .011). These points were useful for stratifying 2 prognostic groups (P between <.001 and .003, Kaplan-Meier). The Cox forward stepwise model selected 3 spectroscopic variables: the intensity values of the points 3.67 ppm at STE (hazard ratio, 2.132; 95% CI, 1.504-3.023), 0.98 ppm at LTE (hazard ratio, 0.499; 95% CI, 0.339-0.736), and 1.25 ppm at LTE (hazard ratio, 0.574; 95% CI, 0.368-0.897)., Conclusions: (1)H-MRS is of value in predicting the length of survival in patients with HGA and could be used to stratify prognostic groups.

Background: Proton Magnetic Resonance (MR) Spectroscopy (MRS) is a widely available technique for those clinical centres equipped with MR scanners. Unlike the rest of MR-based techniques, MRS yields not images but spectra of metabolites in the tissues. In pathological situations, the MRS profile changes and this has been particularly described for brain tumours. However, radiologists are frequently not familiar to the interpretation of MRS data and for this reason, the usefulness of decision-support systems (DSS) in MRS data analysis has been explored., Results: This work presents the INTERPRET DSS version 3.0, analysing the improvements made from its first release in 2002. Version 3.0 is aimed to be a program that 1st, can be easily used with any new case from any MR scanner manufacturer and 2nd, improves the initial analysis capabilities of the first version. The main improvements are an embedded database, user accounts, more diagnostic discrimination capabilities and the possibility to analyse data acquired under additional data acquisition conditions. Other improvements include a customisable graphical user interface (GUI). Most diagnostic problems included have been addressed through a pattern-recognition based approach, in which classifiers based on linear discriminant analysis (LDA) were trained and tested., Conclusions: The INTERPRET DSS 3.0 allows radiologists, medical physicists, biochemists or, generally speaking, any person with a minimum knowledge of what an MR spectrum is, to enter their own SV raw data, acquired at 1.5 T, and to analyse them. The system is expected to help in the categorisation of MR Spectra from abnormal brain masses.