Manifold Learning in MR spectroscopy using nonlinear dimensionality reduction and unsupervised clustering

Purpose To investigate whether nonlinear dimensionality reduction improves unsupervised classification of 1H MRS brain tumor data compared with a linear method. Methods In vivo single-voxel 1H magnetic resonance spectroscopy (55 patients) and 1H magnetic resonance spectroscopy imaging (MRSI) (29 patients) data were acquired from histopathologically diagnosed gliomas. Data reduction using Laplacian eigenmaps (LE) or independent component analysis (ICA) was followed by k-means clustering or agglomerative hierarchical clustering (AHC) for unsupervised learning to assess tumor grade and for tissue type segmentation of MRSI data. Results An accuracy of 93% in classification of glioma grade II and grade IV, with 100% accuracy in distinguishing tumor and normal spectra, was obtained by LE with unsupervised clustering, but not with the combination of k-means and ICA. With 1H MRSI data, LE provided a more linear distribution of data for cluster analysis and better cluster stability than ICA. LE combined with k-means or AHC provided 91% accuracy for classifying tumor grade and 100% accuracy for identifying normal tissue voxels. Color-coded visualization of normal brain, tumor core, and infiltration regions was achieved with LE combined with AHC. Conclusion Purpose To investigate whether nonlinear dimensionality reduction improves unsupervised classification of 1H MRS brain tumor data compared with a linear method. Methods In vivo single-voxel 1H magnetic resonance spectroscopy (55 patients) and 1H magnetic resonance spectroscopy imaging (MRSI) (29 patients) data were acquired from histopathologically diagnosed gliomas. Data reduction using Laplacian eigenmaps (LE) or independent component analysis (ICA) was followed by k-means clustering or agglomerative hierarchical clustering (AHC) for unsupervised learning to assess tumor grade and for tissue type segmentation of MRSI data. Results An accuracy of 93% in classification of glioma grade II and grade IV, with 100% accuracy in distinguishing tumor and normal spectra, was obtained by LE with unsupervised clustering, but not with the combination of k-means and ICA. With 1H MRSI data, LE provided a more linear distribution of data for cluster analysis and better cluster stability than ICA. LE combined with k-means or AHC provided 91% accuracy for classifying tumor grade and 100% accuracy for identifying normal tissue voxels. Color-coded visualization of normal brain, tumor core, and infiltration regions was achieved with LE combined with AHC. Conclusion The LE method is promising for unsupervised clustering to separate brain and tumor tissue with automated color-coding for visualization of 1H MRSI data after cluster analysis

^{1}H MRS spectroscopy in brain tumors : in search of the highest efficacy

Many thousands of published papers, a lot of communications and daily practice show the outstanding role of proton spectroscopy (^{1}H MRS) in neurooncology. Paradoxically, in several clinical centres (neurosurgery ones) we can observe the lack of full acceptance of the reliability of this method. Presumably, the reason for this attitude lies in non-uniform methodology and quite superficial treatment of spectroscopic results. At the beginning of this paper, we would like to stress that fortunately, there is a certain number of publications which order criteria of tumor's malignancy, its tissue characterisation and the degree of malignancy. In this paper we intend to point out the proper ways of improving reliability of ^{1}H MRS in brain tumor diagnosis. Obviously, the quality of spectroscopic exams is influenced by methodology and manner of result interpretation (e.g. just simple visual inspection of MRS spectra is not sufficient for correct diagnosis). Currently, imaging diagnosis of brain tumors is based on many modalities (DTI, perfusion CT, MR, PET, SPECT) which are sufficiently effective in detection and differentiation of the lesions; in this context, we can say that ^{1}H MRS spectroscopy constitutes the element of multimodal diagnostic approach for brain tumors. Nevertheless, the current papers concerning MR spectroscopy stress its significant role in making therapeutic decision, particularly as an alternative for surgery treatment and solution of the recurrence of tumor or/and extent of the tumor's spreading. Over 15 years of experience with ^{1}H MRS in diagnosis of brain tumors convinced us that these methods have very high diagnostic and decision-making value provided that it is properly applied and interpreted correctly. It is also important to choose the statistic methods that will show the best discriminators (markers) for particular pathologies. Authors' intention is to recommend the best ways to form the right diagnosis considering technological requirements, methodology and result interpretation

High-resolution magic angle spinning magnetic resonance spectroscopy detects glycine as a biomarker in brain tumors

The non-essential amino acid neurotransmitter glycine (Gly) may serve as a biomarker for brain tumors. Using 36 biopsies from patients with brain tumors [12 glioblastoma multiforme (GBM); 10 low-grade (LG), including 7 schwannoma and 3 pylocytic astrocytoma; 7 meningioma (MN); 7 brain metastases (MT), including 3 adenocarcinoma and 4 breast cancer] and 9 control biopsies from patients undergoing surgery for epilepsy, we tested the hypothesis that the presence of glycine may distinguish among these brain tumor types. Using high-resolution magic angle spinning (HRMAS) 1H magnetic resonance spectroscopy (MRS), we determined a theoretically optimum echo time (TE) of 50 ms for distinguishing Gly signals from overlapping myo-inositol (Myo) signals and tested our methodology in phantom and biopsy specimens. Quantitative analysis revealed higher levels of Gly in tumor biopsies (all combined) relative to controls; Gly levels were significantly elevated in LG, MT and GBM biopsies (P≤0.05). Residual Myo levels were elevated in LG and MT and reduced in MN and GBM (P<0.05 vs. control levels). We observed higher levels of Gly in GBM as compared to LG tumors (P=0.05). Meanwhile, although Gly levels in GBM and MT did not differ significantly from each other, the Gly:Myo ratio did distinguish GBM from MT (P<0.003) and from all other groups, a distinction that has not been adequately made previously. We conclude from these findings that Gly can serve as a biomarker for brain tumors and that the Gly:Myo ratio may be a useful index for brain tumor classification

Metabolic profiling of human brain metastases using in vivo proton MR spectroscopy at 3T

<p>Abstract</p> <p>Background</p> <p>Metastases to the central nervous system from different primary cancers are an oncologic challenge as the overall prognosis for these patients is generally poor. The incidence of brain metastases varies with type of primary cancer and is probably increasing due to improved therapies of extracranial metastases prolonging patient's overall survival and thereby time for brain metastases to develop. In addition, the greater access to improved neuroimaging techniques can provide earlier diagnosis. The aim of this study was to investigate the feasibility of using proton magnetic resonance spectroscopy (MRS) and multivariate analyses to characterize brain metastases originating from different primary cancers, to assess changes in spectra during radiation treatment and to correlate the spectra to clinical outcome after treatment.</p> <p>Methods</p> <p>Patients (n = 26) with brain metastases were examined using single voxel MRS at a 3T clinical MR system. Five patients were excluded due to poor spectral quality. The spectra were obtained before start (n = 21 patients), immediately after (n = 6 patients) and two months after end of treatment (n = 4 patients). Principal component analysis (PCA) and partial least square regression analysis (PLS) were applied in order to identify clustering of spectra due to origin of metastases and to relate clinical outcome (survival) of the patients to spectral data from the first MR examination.</p> <p>Results</p> <p>The PCA results indicated that brain metastases from primary lung and breast cancer were separated into two clusters, while the metastases from malignant melanomas showed no uniformity. The PLS analysis showed a significant correlation between MR spectral data and survival five months after MRS before start of treatment.</p> <p>Conclusion</p> <p>MRS determined metabolic profiles analysed by PCA and PLS might give valuable clinical information when planning and evaluating the treatment of brain metastases, and also when deciding to terminate further therapies.</p

Feature selection in proton magnetic resonance spectroscopy data of brain tumors

In cancer diagnosis, classification of the different tumor types is of great importance. An accurate prediction of different tumor types provides better treatment and may minimize the negative impact of incorrectly targeted toxic or aggressive treatments. Moreover, the correct prediction of cancer types using non-invasive information –e.g. 1H-MRS data– could avoid patients to suffer collateral problems derived from exploration techniques that require surgery. A Feature Selection Algorithm specially designed to be use in 1H-MRS Proton Magnetic Resonance Spectroscopy data of brain tumors is presented. It takes advantage of a highly distinctive aspect in this data: some metabolite levels are notoriously different between types of tumors. Experimental read- ings on an international dataset show highly competitive models in terms of accuracy, complexity and medical interpretability.Postprint (author’s final draft

Improving the Clinical Use of Magnetic Resonance Spectroscopy for the Analysis of Brain Tumours using Machine Learning and Novel Post-Processing Methods

Magnetic Resonance Spectroscopy (MRS) provides unique and clinically relevant information for the assessment of several diseases. However, using the currently available tools, MRS processing and analysis is time-consuming and requires profound expert knowledge. For these two reasons, MRS did not gain general acceptance as a mainstream diagnostic technique yet, and the currently available clinical tools have seen little progress during the past years. MRS provides localized chemical information non-invasively, making it a valuable technique for the assessment of various diseases and conditions, namely brain, prostate and breast cancer, and metabolic diseases affecting the brain. In brain cancer, MRS is normally used for: (1.) differentiation between tumors and non-cancerous lesions, (2.) tumor typing and grading, (3.) differentiation between tumor-progression and radiation necrosis, and (4.) identification of tumor infiltration. Despite the value of MRS for these tasks, susceptibility differences associated with tissue-bone and tissue-air interfaces, as well as with the presence of post-operative paramagnetic particles, affect the quality of brain MR spectra and consequently reduce their clinical value. Therefore, the proper quality management of MRS acquisition and processing is essential to achieve unambiguous and reproducible results. In this thesis, special emphasis was placed on this topic. This thesis addresses some of the major problems that limit the use of MRS in brain tumors and focuses on the use of machine learning for the automation of the MRS processing pipeline and for assisting the interpretation of MRS data. Three main topics were investigated: (1.) automatic quality control of MRS data, (2.) identification of spectroscopic patterns characteristic of different tissue-types in brain tumors, and (3.) development of a new approach for the detection of tumor-related changes in GBM using MRSI data. The first topic tackles the problem of MR spectra being frequently affected by signal artifacts that obscure their clinical information content. Manual identification of these artifacts is subjective and is only practically feasible for single-voxel acquisitions and in case the user has an extensive experience with MRS. Therefore, the automatic distinction between data of good or bad quality is an essential step for the automation of MRS processing and routine reporting. The second topic addresses the difficulties that arise while interpreting MRS results: the interpretation requires expert knowledge, which is not available at every site. Consequently, the development of methods that enable the easy comparison of new spectra with known spectroscopic patterns is of utmost importance for clinical applications of MRS. The third and last topic focuses on the use of MRSI information for the detection of tumor-related effects in the periphery of brain tumors. Several research groups have shown that MRSI information enables the detection of tumor infiltration in regions where structural MRI appears normal. However, many of the approaches described in the literature make use of only a very limited amount of the total information contained in each MR spectrum. Thus, a better way to exploit MRSI information should enable an improvement in the detection of tumor borders, and consequently improve the treatment of brain tumor patients. The development of the methods described was made possible by a novel software tool for the combined processing of MRS and MRI: SpectrIm. This tool, which is currently distributed as part of the jMRUI software suite (www.jmrui.eu), is ubiquitous to all of the different methods presented and was one of the main outputs of the doctoral work. Overall, this thesis presents different methods that, when combined, enable the full automation of MRS processing and assist the analysis of MRS data in brain tumors. By allowing clinical users to obtain more information from MRS with less effort, this thesis contributes to the transformation of MRS into an important clinical tool that may be available whenever its information is of relevance for patient management

The INTERPRET Decision-Support System version 3.0 for evaluation of Magnetic Resonance Spectroscopy data from human brain tumours and other abnormal brain masses.

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

A Novel Semi-Supervised Methodology for Extracting Tumor Type-Specific MRS Sources in Human Brain Data

BackgroundThe 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 FindingsNon-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/SignificanceWe 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