Characterization of glioblastoma (GB) response to treatment is a key factor for improving patients' survival and prognosis. MRI and magnetic resonance spectroscopic imaging (MRSI) provide morphologic and metabolic profiles of GB but usually fail to produce unequivocal biomarkers of response. The purpose of this work is to provide proof of concept of the ability of a semi‐supervised signal source extraction methodology to produce images with robust recognition of response to temozolomide (TMZ) in a preclinical GB model. A total of 38 female C57BL/6 mice were used in this study. The semi‐supervised methodology extracted the required sources from a training set consisting of MRSI grids from eight GL261 GBs treated with TMZ, and six control untreated GBs. Three different sources (normal brain parenchyma, actively proliferating GB and GB responding to treatment) were extracted and used for calculating nosologic maps representing the spatial response to treatment. These results were validated with an independent test set (7 control and 17 treated cases) and correlated with histopathology. Major differences between the responder and non‐responder sources were mainly related to the resonances of mobile lipids (MLs) and polyunsaturated fatty acids in MLs (0.9, 1.3 and 2.8 ppm). Responding tumors showed significantly lower mitotic (3.3 ± 2.9 versus 14.1 ± 4.2 mitoses/field) and proliferation rates (29.8 ± 10.3 versus 57.8 ± 5.4%) than control untreated cases. The methodology described in this work is able to produce nosological images of response to TMZ in GL261 preclinical GBs and suitably correlates with the histopathological analysis of tumors. A similar strategy could be devised for monitoring response to treatment in patients. Copyright © 2016 John Wiley & Sons, Ltd.

BackgroundPattern 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 FindingsA 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/SignificanceThe 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.