Unsupervised nosologic imaging for glioma diagnosis

Establishing Imaging Biomarkers of Host Immune System Efficacy during Glioblastoma Therapy Response: Challenges, Obstacles and Future Perspectives

This hypothesis proposal addresses three major questions: (1) Why do we need imaging biomarkers for assessing the efficacy of immune system participation in glioblastoma therapy response? (2) Why are they not available yet? and (3) How can we produce them? We summarize the literature data supporting the claim that the immune system is behind the efficacy of most successful glioblastoma therapies but, unfortunately, there are no current short-term imaging biomarkers of its activity. We also discuss how using an immunocompetent murine model of glioblastoma, allowing the cure of mice and the generation of immune memory, provides a suitable framework for glioblastoma therapy response biomarker studies. Both magnetic resonance imaging and magnetic resonance-based metabolomic data (i.e., magnetic resonance spectroscopic imaging) can provide non-invasive assessments of such a system. A predictor based in nosological images, generated from magnetic resonance spectroscopic imaging analyses and their oscillatory patterns, should be translational to clinics. We also review hurdles that may explain why such an oscillatory biomarker was not reported in previous imaging glioblastoma work. Single shot explorations that neglect short-term oscillatory behavior derived from immune system attack on tumors may mislead actual response extent detection. Finally, we consider improvements required to properly predict immune system-mediated early response (1–2 weeks) to therapy. The sensible use of improved biomarkers may enable translatable evidence-based therapeutic protocols, with the possibility of extending preclinical results to human patients.

Recovering Hidden Diagonal Structures via Non-Negative Matrix Factorization with Multiple Constraints

Revealing data with intrinsically diagonal block structures is particularly useful for analyzing groups of highly correlated variables. Earlier researches based on non-negative matrix factorization (NMF) have been shown to be effective in representing such data by decomposing the observed data into two factors, where one factor is considered to be the feature and the other the expansion loading from a linear algebra perspective. If the data are sampled from multiple independent subspaces, the loading factor would possess a diagonal structure under an ideal matrix decomposition. However, the standard NMF method and its variants have not been reported to exploit this type of data via direct estimation. To address this issue, a non-negative matrix factorization with multiple constraints model is proposed in this paper. The constraints include an sparsity norm on the feature matrix and a total variational norm on each column of the loading matrix. The proposed model is shown to be capable of efficiently recovering diagonal block structures hidden in observed samples. An efficient numerical algorithm using the alternating direction method of multipliers model is proposed for optimizing the new model. Compared with several benchmark models, the proposed method performs robustly and effectively for simulated and real biological data.

Nonnegative Canonical Polyadic Decomposition for Tissue-Type Differentiation in Gliomas

Magnetic resonance spectroscopic imaging (MRSI) reveals chemical information that characterizes different tissue types in brain tumors. Blind source separation techniques are used to extract the tissue-specific profiles and their corresponding distribution from the MRSI data. We focus on automatic detection of the tumor, necrotic and normal brain tissue types by constructing a 3D MRSI tensor from in vivo 2D-MRSI data of individual glioma patients. Nonnegative canonical polyadic decomposition (NCPD) is applied to the MRSI tensor to differentiate various tissue types. An in vivo study shows that NCPD has better performance in identifying tumor and necrotic tissue type in glioma patients compared to previous matrix-based decompositions, such as nonnegative matrix factorization and hierarchical nonnegative matrix factorization.

Hierarchical non-negative matrix factorization applied to three-dimensional 3 T MRSI data for automatic tissue characterization of the prostate

In this study non-negative matrix factorization (NMF) was hierarchically applied to simulated and in vivo three-dimensional 3 T MRSI data of the prostate to extract patterns for tumour and benign tissue and to visualize their spatial distribution. Our studies show that the hierarchical scheme provides more reliable tissue patterns than those obtained by performing only one NMF level. We compared the performance of three different NMF implementations in terms of pattern detection accuracy and efficiency when embedded into the same kind of hierarchical scheme. The simulation and in vivo results show that the three implementations perform similarly, although one of them is more robust and better pinpoints the most aggressive tumour voxel(s) in the dataset. Furthermore, they are able to detect tumour and benign tissue patterns even in spectra with lipid artefacts. Copyright © 2016 John Wiley & Sons, Ltd.

Tensor based tumor tissue type differentiation using magnetic resonance spectroscopic imaging

Magnetic resonance spectroscopic imaging (MRSI) has the potential to characterise different tissue types in brain tumors. Blind source separation techniques are used to extract the specific tissue profiles and their corresponding distribution from the MRSI data. A 3-dimensional MRSI tensor is constructed from in vivo 2D-MRSI data of individual tumor patients. Non-negative canonical polyadic decomposition (NCPD) with common factor in mode-1 and mode-2 and l(1) regularization on mode-3 is applied on the MRSI tensor to differentiate various tissue types. Initial in vivo study shows that NCPD has better performance in identifying tumor and necrotic tissue type in high grade glioma patients compared to previous matrix-based decompositions, such as non-negative matrix factorization and hierarchical non-negative matrix factorization.

Hierarchical non-negative matrix factorization to characterize brain tumor heterogeneity using multi-parametric MRI

Tissue characterization in brain tumors and, in particular, in high-grade gliomas is challenging as a result of the co-existence of several intra-tumoral tissue types within the same region and the high spatial heterogeneity. This study presents a method for the detection of the relevant tumor substructures (i.e. viable tumor, necrosis and edema), which could be of added value for the diagnosis, treatment planning and follow-up of individual patients. Twenty-four patients with glioma [10 low-grade gliomas (LGGs), 14 high-grade gliomas (HGGs)] underwent a multi-parametric MRI (MP-MRI) scheme, including conventional MRI (cMRI), perfusion-weighted imaging (PWI), diffusion kurtosis imaging (DKI) and short-TE (1)H MRSI. MP-MRI parameters were derived: T2, T1 + contrast, fluid-attenuated inversion recovery (FLAIR), relative cerebral blood volume (rCBV), mean diffusivity (MD), fractional anisotropy (FA), mean kurtosis (MK) and the principal metabolites lipids (Lip), lactate (Lac), N-acetyl-aspartate (NAA), total choline (Cho), etc. Hierarchical non-negative matrix factorization (hNMF) was applied to the MP-MRI parameters, providing tissue characterization on a patient-by-patient and voxel-by-voxel basis. Tissue-specific patterns were obtained and the spatial distribution of each tissue type was visualized by means of abundance maps. Dice scores were calculated by comparing tissue segmentation derived from hNMF with the manual segmentation by a radiologist. Correlation coefficients were calculated between each pathologic tissue source and the average feature vector within the corresponding tissue region. For the patients with HGG, mean Dice scores of 78%, 85% and 83% were obtained for viable tumor, the tumor core and the complete tumor region. The mean correlation coefficients were 0.91 for tumor, 0.97 for necrosis and 0.96 for edema. For the patients with LGG, a mean Dice score of 85% and mean correlation coefficient of 0.95 were found for the tumor region. hNMF was also applied to reduced MRI datasets, showing the added value of individual MRI modalities.

Multimodal imaging of gliomas in the context of evolving cellular and molecular therapies

The vast majority of malignant gliomas relapse after surgery and standard radio-chemotherapy. Novel molecular and cellular therapies are thus being developed, targeting specific aspects of tumor growth. While histopathology remains the gold standard for tumor classification, neuroimaging has over the years taken a central role in the diagnosis and treatment follow up of brain tumors. It is used to detect and localize lesions, define the target area for biopsies, plan surgical and radiation interventions and assess tumor progression and treatment outcome. In recent years the application of novel drugs including anti-angiogenic agents that affect the tumor vasculature, has drastically modulated the outcome of brain tumor imaging. To properly evaluate the effects of emerging experimental therapies and successfully support treatment decisions, neuroimaging will have to evolve. Multi-modal imaging systems with existing and new contrast agents, molecular tracers, technological advances and advanced data analysis can all contribute to the establishment of disease relevant biomarkers that will improve disease management and patient care. In this review, we address the challenges of glioma imaging in the context of novel molecular and cellular therapies, and take a prospective look at emerging experimental and pre-clinical imaging techniques that bear the promise of meeting these challenges.

Microvascular MRI and unsupervised clustering yields histology-resembling images in two rat models of glioma

Imaging heterogeneous cancer lesions is a real challenge. For diagnosis, histology often remains the reference, but it is widely acknowledged that biopsies are not reliable. There is thus a strong interest in establishing a link between clinical in vivo imaging and the biologic properties of tissues. In this study, we propose to construct histology-resembling images based on tissue microvascularization, a magnetic resonance imaging (MRI) accessible source of contrast. To integrate the large amount of information collected with microvascular MRI, we combined a manual delineation of a spatial region of interest with an unsupervised, model-based cluster analysis (Mclust). This approach was applied to two rat models of glioma (C6 and F98). Six MRI parameters were mapped: apparent diffusion coefficient, vessel wall permeability, cerebral blood volume fraction, cerebral blood flow, tissular oxygen saturation, and cerebral metabolic rate of oxygen. Five clusters, defined by their MRI features, were found to correspond to specific histologic features, and revealed intratumoral spatial structures. These results suggest that the presence of a cluster within a tumor can be used to assess the presence of a tissue type. In addition, the cluster composition, i.e., a signature of the intratumoral structure, could be used to characterize tumor models as histology does.

Fast spectroscopic multiple analysis (FASMA) for brain tumor classification: a clinical decision support system utilizing multi-parametric 3T MR data

INTRODUCTION

A clinical decision support system (CDSS) for brain tumor classification can be used to assist in the diagnosis and grading of brain tumors. A Fast Spectroscopic Multiple Analysis (FASMA) system that uses combinations of multiparametric MRI data sets was developed as a CDSS for brain tumor classification.

METHODS

MRI metabolic ratios and spectra, from long and short TE, respectively, as well as diffusion and perfusion data were acquired from the intratumoral and peritumoral area of 126 patients with untreated intracranial tumors. These data were categorized based on the pathology, and different machine learning methods were evaluated regarding their classification performance for glioma grading and differentiation of infiltrating versus non-infiltrating lesions. Additional databases were embedded to the system, including updated literature values of the related MR parameters and typical tumor characteristics (imaging and histological), for further comparisons. Custom Graphical User Interface (GUI) layouts were developed to facilitate classification of the unknown cases based on the user's available MR data.

RESULTS

The highest classification performance was achieved with a support vector machine (SVM) using the combination of all MR features. FASMA correctly classified 89 and 79% in the intratumoral and peritumoral area, respectively, for cases from an independent test set. FASMA produced the correct diagnosis, even in the misclassified cases, since discrimination between infiltrative versus non-infiltrative cases was possible.

CONCLUSIONS

FASMA is a prototype CDSS, which integrates complex quantitative MR data for brain tumor characterization. FASMA was developed as a diagnostic assistant that provides fast analysis, representation and classification for a set of MR parameters. This software may serve as a teaching tool on advanced MRI techniques, as it incorporates additional information regarding typical tumor characteristics derived from the literature.

Clinical decision support systems for brain tumor characterization using advanced magnetic resonance imaging techniques

In recent years, advanced magnetic resonance imaging (MRI) techniques, such as magnetic resonance spectroscopy, diffusion weighted imaging, diffusion tensor imaging and perfusion weighted imaging have been used in order to resolve demanding diagnostic problems such as brain tumor characterization and grading, as these techniques offer a more detailed and non-invasive evaluation of the area under study. In the last decade a great effort has been made to import and utilize intelligent systems in the so-called clinical decision support systems (CDSS) for automatic processing, classification, evaluation and representation of MRI data in order for advanced MRI techniques to become a part of the clinical routine, since the amount of data from the aforementioned techniques has gradually increased. Hence, the purpose of the current review article is two-fold. The first is to review and evaluate the progress that has been made towards the utilization of CDSS based on data from advanced MRI techniques. The second is to analyze and propose the future work that has to be done, based on the existing problems and challenges, especially taking into account the new imaging techniques and parameters that can be introduced into intelligent systems to significantly improve their diagnostic specificity and clinical application.

Data analysis and tissue type assignment for glioblastoma multiforme

Glioblastoma multiforme (GBM) is characterized by high infiltration. The interpretation of MRSI data, especially for GBMs, is still challenging. Unsupervised methods based on NMF by Li et al. (2013, NMR in Biomedicine) and Li et al. (2013, IEEE Transactions on Biomedical Engineering) have been proposed for glioma recognition, but the tissue types is still not well interpreted. As an extension of the previous work, a tissue type assignment method is proposed for GBMs based on the analysis of MRSI data and tissue distribution information. The tissue type assignment method uses the values from the distribution maps of all three tissue types to interpret all the information in one new map and color encodes each voxel to indicate the tissue type. Experiments carried out on in vivo MRSI data show the feasibility of the proposed method. This method provides an efficient way for GBM tissue type assignment and helps to display information of MRSI data in a way that is easy to interpret.