Imaging-genomics reveals driving pathways of MRI derived volumetric tumor phenotype features in Glioblastoma

BackgroundGlioblastoma (GBM) tumors exhibit strong phenotypic differences that can be quantified using magnetic resonance imaging (MRI), but the underlying biological drivers of these imaging phenotypes remain largely unknown. An Imaging-Genomics analysis was performed to reveal the mechanistic associations between MRI derived quantitative volumetric tumor phenotype features and molecular pathways.MethodsOne hundred fourty one patients with presurgery MRI and survival data were included in our analysis. Volumetric features were defined, including the necrotic core (NE), contrast-enhancement (CE), abnormal tumor volume assessed by post-contrast T1w (tumor bulk or TB), tumor-associated edema based on T2-FLAIR (ED), and total tumor volume (TV), as well as ratios of these tumor components. Based on gene expression where available (n = 91), pathway associations were assessed using a preranked gene set enrichment analysis. These results were put into context of molecular subtypes in GBM and prognostication.ResultsVolumetric features were significantly associated with diverse sets of biological processes (FDR < 0.05). While NE and TB were enriched for immune response pathways and apoptosis, CE was associated with signal transduction and protein folding processes. ED was mainly enriched for homeostasis and cell cycling pathways. ED was also the strongest predictor of molecular GBM subtypes (AUC = 0.61). CE was the strongest predictor of overall survival (C-index = 0.6; Noether test, p = 4x10−4).ConclusionGBM volumetric features extracted from MRI are significantly enriched for information about the biological state of a tumor that impacts patient outcomes. Clinical decision-support systems could exploit this information to develop personalized treatment strategies on the basis of noninvasive imaging.

[1]  Keith L. Ligon,et al.  Somatic mutations associated with MRI-derived volumetric features in glioblastoma , 2015, Neuroradiology.

[2]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Shawn M. Gillespie,et al.  Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma , 2014, Science.

[4]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[5]  Aedín C. Culhane,et al.  survcomp: an R/Bioconductor package for performance assessment and comparison of survival models , 2011, Bioinform..

[6]  C. Kitanaka,et al.  Apoptosis in cancer. , 1996, Human cell.

[7]  L. Deangelis,et al.  Glioblastoma and other malignant gliomas: a clinical review. , 2013, JAMA.

[8]  V. P. Collins,et al.  Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics , 2013, Proceedings of the National Academy of Sciences.

[9]  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.

[10]  M. Gilbert Faculty Opinions recommendation of Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. , 2014 .

[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]  S. Plevritis,et al.  Glioblastoma multiforme: exploratory radiogenomic analysis by using quantitative image features. , 2014, Radiology.

[13]  Jill S. Barnholtz-Sloan,et al.  Epidemiologic and Molecular Prognostic Review of Glioblastoma , 2014, Cancer Epidemiology, Biomarkers & Prevention.

[14]  H. Schreiber,et al.  Innate and adaptive immune cells in the tumor microenvironment , 2013, Nature Immunology.

[15]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[16]  Olivier Gevaert,et al.  Addition of MR imaging features and genetic biomarkers strengthens glioblastoma survival prediction in TCGA patients. , 2015, Journal of neuroradiology. Journal de neuroradiologie.

[17]  Lei Wang,et al.  Identifying the Association of Contrast Enhancement with Vascular Endothelia Growth Factor Expression in Anaplastic Gliomas: A Volumetric Magnetic Resonance Imaging Analysis , 2015, PloS one.

[18]  Tom Fawcett,et al.  An introduction to ROC analysis , 2006, Pattern Recognit. Lett..

[19]  Helga Thorvaldsdóttir,et al.  Molecular signatures database (MSigDB) 3.0 , 2011, Bioinform..

[20]  Ferenc A. Jolesz,et al.  Radiogenomic Mapping of Edema/Cellular Invasion MRI-Phenotypes in Glioblastoma Multiforme , 2011, PloS one.

[21]  E. Wong,et al.  Personalized medicine for glioblastoma: current challenges and future opportunities. , 2013, Current molecular medicine.

[22]  William D. Dunn,et al.  MR imaging predictors of molecular profile and survival: multi-institutional study of the TCGA glioblastoma data set. , 2013, Radiology.

[23]  Benjamin E. Gross,et al.  Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal , 2013, Science Signaling.

[24]  Ritesh Kumar,et al.  Evaluation of outcome and prognostic factors in patients of glioblastoma multiforme: A single institution experience , 2013, Journal of neurosciences in rural practice.

[25]  Sabrina S Wilson Radiology , 1938, Glasgow Medical Journal.

[26]  Paul Kinahan,et al.  Radiomics: Images Are More than Pictures, They Are Data , 2015, Radiology.

[27]  R. Mirimanoff,et al.  MGMT gene silencing and benefit from temozolomide in glioblastoma. , 2005, The New England journal of medicine.

[28]  Jill S Barnholtz-Sloan,et al.  Whole-genome and multisector exome sequencing of primary and post-treatment glioblastoma reveals patterns of tumor evolution , 2015, Genome research.

[29]  K. Aldape,et al.  Identification of noninvasive imaging surrogates for brain tumor gene-expression modules , 2008, Proceedings of the National Academy of Sciences.

[30]  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.

[31]  Neema Jamshidi,et al.  Illuminating radiogenomic characteristics of glioblastoma multiforme through integration of MR imaging, messenger RNA expression, and DNA copy number variation. , 2013, Radiology.

[32]  Paul S Mischel,et al.  Relationship between gene expression and enhancement in glioblastoma multiforme: exploratory DNA microarray analysis. , 2008, Radiology.

[33]  A. Levitzki,et al.  Signal transduction therapy of cancer. , 2010, Molecular aspects of medicine.

[34]  Sylvia Drabycz,et al.  An analysis of image texture, tumor location, and MGMT promoter methylation in glioblastoma using magnetic resonance imaging , 2010, NeuroImage.

[35]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[36]  Z L Gokaslan,et al.  A multivariate analysis of 416 patients with glioblastoma multiforme: prognosis, extent of resection, and survival. , 2001, Journal of neurosurgery.

[37]  David Gutman,et al.  Genomic mapping and survival prediction in glioblastoma: molecular subclassification strengthened by hemodynamic imaging biomarkers. , 2013, Radiology.

[38]  F. Harrell,et al.  Evaluating the yield of medical tests. , 1982, JAMA.

[39]  D. Kaufman,et al.  Cell cycle control, DNA repair and initiation of carcinogenesis , 1993, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[40]  Benjamin M Ellingson,et al.  Identifying the mesenchymal molecular subtype of glioblastoma using quantitative volumetric analysis of anatomic magnetic resonance images. , 2013, Neuro-oncology.