ONCOhabitats: A system for glioblastoma heterogeneity assessment through MRI

BACKGROUND Neuroimaging analysis is currently crucial for an early assessment of glioblastoma, to help improving treatment and tumor follow-up. To this end, multiple functional and morphological MRI sequences are usually employed, requiring the development of automated tools capable to extract the relevant information from these sources. In this work we present ONCOhabitats (https://www.oncohabitats.upv.es): an online open access system for glioblastoma analysis based on MRI data. METHODS ONCOhabitats provides two main services for untreated glioblastomas: (1) malignant tissue segmentation, and (2) vascular heterogeneity assessment of the tumor. The segmentation service implements a deep patch-wise 3D Convolutional Neural Network with residual connections. The vascular heterogeneity assessment service implements the Hemodynamic Tissue Signature (HTS) method patented in P201431289, which aims to identify habitats within the tumor with early prognostic capabilities. RESULTS The segmentation service was validated against the BRATS 2017 reference dataset, showing comparable results with current state-of-the-art methods (whole tumor Dice segmentation: 0.89). The vascular heterogeneity assessment service was validated in a retrospective cohort of 50 patients, in a study focused on predicting patient overall survival based on the HTS habitats. Cox proportional hazard regression analysis and Kaplan-Meier survival study showed significant positive correlations (p-value <.05) between the HTS habitats and patient overall survival. ONCOhabitats system also generates radiological reports for each service, including volumetries and perfusion measurements of the different regions of the lesion. CONCLUSION ONCOhabitats system provides open-access services for glioblastoma heterogeneity assessment, implementing consolidated state-of-the-art techniques for medical image analysis. Additionally, we also give access to the scientific community to our computational resources, offering a computational capacity of about 300 cases per day.

[1]  J. Seoane,et al.  Glioblastoma Multiforme: A Look Inside Its Heterogeneous Nature , 2014, Cancers.

[2]  Yang Liu,et al.  Relationship between Glioblastoma Heterogeneity and Survival Time: An MR Imaging Texture Analysis , 2017, American Journal of Neuroradiology.

[3]  G Johnson,et al.  Glial neoplasms: dynamic contrast-enhanced T2*-weighted MR imaging. , 1999, Radiology.

[4]  Guang Yang,et al.  Automatic Brain Tumor Detection and Segmentation Using U-Net Based Fully Convolutional Networks , 2017, MIUA.

[5]  Simona Marzi,et al.  Early perfusion changes in patients with recurrent high-grade brain tumor treated with Bevacizumab: preliminary results by a quantitative evaluation , 2012, Journal of experimental & clinical cancer research : CR.

[6]  Thomas Brox,et al.  U-Net: Convolutional Networks for Biomedical Image Segmentation , 2015, MICCAI.

[7]  Jianhui Zhong,et al.  Changes in relative cerebral blood volume 1 month after radiation-temozolomide therapy can help predict overall survival in patients with glioblastoma. , 2010, Radiology.

[8]  Brian B. Avants,et al.  N4ITK: Improved N3 Bias Correction , 2010, IEEE Transactions on Medical Imaging.

[9]  Luis Martí-Bonmatí,et al.  Glioblastoma: Vascular Habitats Detected at Preoperative Dynamic Susceptibility-weighted Contrast-enhanced Perfusion MR Imaging Predict Survival. , 2018, Radiology.

[10]  Dafna Ben Bashat,et al.  Differentiation between vasogenic-edema versus tumor-infiltrative area in patients with glioblastoma during bevacizumab therapy: a longitudinal MRI study. , 2014, European journal of radiology.

[11]  Scott N. Hwang,et al.  Outcome prediction in patients with glioblastoma by using imaging, clinical, and genomic biomarkers: focus on the nonenhancing component of the tumor. , 2014, Radiology.

[12]  Carlos Sáez,et al.  Stability metrics for multi-source biomedical data based on simplicial projections from probability distribution distances , 2017, Statistical methods in medical research.

[13]  R. Mirimanoff,et al.  Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. , 2005, The New England journal of medicine.

[14]  Guang Yang,et al.  MRI Brain Tumor Segmentation and Patient Survival Prediction Using Random Forests and Fully Convolutional Networks , 2017, BrainLes@MICCAI.

[15]  Glyn Johnson,et al.  Relative cerebral blood volume measurements in intracranial mass lesions: interobserver and intraobserver reproducibility study. , 2002, Radiology.

[16]  Guang Yang,et al.  Supervised Learning based Multimodal MRI Brain Tumour Segmentation using Texture Features from Supervoxels# , 2018, Comput. Methods Programs Biomed..

[17]  Brian B. Avants,et al.  Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain , 2008, Medical Image Anal..

[18]  Susan M. Chang,et al.  Dynamic susceptibility-weighted perfusion imaging of high-grade gliomas: characterization of spatial heterogeneity. , 2005, AJNR. American journal of neuroradiology.

[19]  S. Bauer,et al.  A survey of MRI-based medical image analysis for brain tumor studies , 2013, Physics in medicine and biology.

[20]  Brian B. Avants,et al.  The Multimodal Brain Tumor Image Segmentation Benchmark (BRATS) , 2015, IEEE Transactions on Medical Imaging.

[21]  Nikolas P. Galatsanos,et al.  A Class-Adaptive Spatially Variant Mixture Model for Image Segmentation , 2007, IEEE Transactions on Image Processing.

[22]  Guang Yang,et al.  Tissue-type mapping of gliomas , 2018, NeuroImage: Clinical.

[23]  A. Bjørnerud,et al.  Glioma grading by using histogram analysis of blood volume heterogeneity from MR-derived cerebral blood volume maps. , 2008, Radiology.

[24]  Wanyong Shin,et al.  Quantitative cerebral MR perfusion imaging: Preliminary results in stroke , 2010, Journal of magnetic resonance imaging : JMRI.

[25]  J. Gee,et al.  The Insight ToolKit image registration framework , 2014, Front. Neuroinform..

[26]  S. Gabriel,et al.  Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. , 2010, Cancer cell.

[27]  Linda Knutsson,et al.  Absolute quantification of perfusion using dynamic susceptibility contrast MRI: pitfalls and possibilities , 2010, Magnetic Resonance Materials in Physics, Biology and Medicine.

[28]  B. Rosen,et al.  Tracer arrival timing‐insensitive technique for estimating flow in MR perfusion‐weighted imaging using singular value decomposition with a block‐circulant deconvolution matrix , 2003, Magnetic resonance in medicine.

[29]  Pierrick Coupé,et al.  Author manuscript, published in "Journal of Magnetic Resonance Imaging 2010;31(1):192-203" DOI: 10.1002/jmri.22003 Adaptive Non-Local Means Denoising of MR Images with Spatially Varying Noise Levels , 2010 .

[30]  R M Weisskoff,et al.  Relative cerebral blood volume maps corrected for contrast agent extravasation significantly correlate with glioma tumor grade, whereas uncorrected maps do not. , 2006, AJNR. American journal of neuroradiology.

[31]  Luke Macyszyn,et al.  Pattern analysis of dynamic susceptibility contrast-enhanced MR imaging demonstrates peritumoral tissue heterogeneity. , 2014, Radiology.

[32]  J. Raizer,et al.  Glioblastoma: a method for predicting response to antiangiogenic chemotherapy by using MR perfusion imaging--pilot study. , 2010, Radiology.

[33]  D. Busam,et al.  An Integrated Genomic Analysis of Human Glioblastoma Multiforme , 2008, Science.

[34]  Max Wintermark,et al.  Multicenter imaging outcomes study of The Cancer Genome Atlas glioblastoma patient cohort: imaging predictors of overall and progression-free survival. , 2015, Neuro-oncology.

[35]  Brian B. Avants,et al.  An Open Source Multivariate Framework for n-Tissue Segmentation with Evaluation on Public Data , 2011, Neuroinformatics.

[36]  Arno Klein,et al.  Large-scale evaluation of ANTs and FreeSurfer cortical thickness measurements , 2014, NeuroImage.

[37]  D. Zagzag,et al.  Mechanisms of glioma-associated neovascularization. , 2012, The American journal of pathology.

[38]  Ninon Burgos,et al.  New advances in the Clinica software platform for clinical neuroimaging studies , 2019 .

[39]  J. Barnholtz-Sloan,et al.  CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2007-2011. , 2012, Neuro-oncology.

[40]  J. Lemée,et al.  Intratumoral heterogeneity in glioblastoma: don't forget the peritumoral brain zone. , 2015, Neuro-oncology.

[41]  Jian Sun,et al.  Deep Residual Learning for Image Recognition , 2015, 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR).

[42]  Helena Borges,et al.  Glioblastoma cells: a heterogeneous and fatal tumor interacting with the parenchyma. , 2011, Life sciences.

[43]  M R Segal,et al.  Distinguishing Recurrent Intra-Axial Metastatic Tumor from Radiation Necrosis Following Gamma Knife Radiosurgery Using Dynamic Susceptibility-Weighted Contrast-Enhanced Perfusion MR Imaging , 2008, American Journal of Neuroradiology.

[44]  Glyn Johnson,et al.  Glioma grading: sensitivity, specificity, and predictive values of perfusion MR imaging and proton MR spectroscopic imaging compared with conventional MR imaging. , 2003, AJNR. American journal of neuroradiology.

[45]  J E Heiserman,et al.  Relative Cerebral Blood Volume Values to Differentiate High-Grade Glioma Recurrence from Posttreatment Radiation Effect: Direct Correlation between Image-Guided Tissue Histopathology and Localized Dynamic Susceptibility-Weighted Contrast-Enhanced Perfusion MR Imaging Measurements , 2009, American Journal of Neuroradiology.

[46]  Deric M. Park,et al.  The Evidence of Glioblastoma Heterogeneity , 2015, Scientific Reports.

[47]  Stephen L. Brown,et al.  Acute Temporal Changes of MRI-Tracked Tumor Vascular Parameters after Combined Anti-angiogenic and Radiation Treatments in a Rat Glioma Model: Identifying Signatures of Synergism , 2016, Radiation Research.

[48]  Guang Yang,et al.  Automated brain tumour detection and segmentation using superpixel-based extremely randomized trees in FLAIR MRI , 2016, International Journal of Computer Assisted Radiology and Surgery.

[49]  A. Kyritsis,et al.  Mechanisms of angiogenesis in gliomas , 2006, Journal of Neuro-Oncology.