Identification of potential biomarkers related to glioma survival by gene expression profile analysis

BackgroundRecent studies have proposed several gene signatures as biomarkers for different grades of gliomas from various perspectives. However, most of these genes can only be used appropriately for patients with specific grades of gliomas.MethodsIn this study, we aimed to identify survival-relevant genes shared between glioblastoma multiforme (GBM) and lower-grade glioma (LGG), which could be used as potential biomarkers to classify patients into different risk groups. Cox proportional hazard regression model (Cox model) was used to extract relative genes, and effectiveness of genes was estimated against random forest regression. Finally, risk models were constructed with logistic regression.ResultsWe identified 104 key genes that were shared between GBM and LGG, which could be significantly correlated with patients’ survival based on next-generation sequencing data obtained from The Cancer Genome Atlas for gene expression analysis. The effectiveness of these genes in the survival prediction of GBM and LGG was evaluated, and the average receiver operating characteristic curve (ROC) area under the curve values ranged from 0.7 to 0.8. Gene set enrichment analysis revealed that these genes were involved in eight significant pathways and 23 molecular functions. Moreover, the expressions of ten (CTSZ, EFEMP2, ITGA5, KDELR2, MDK, MICALL2, MAP 2 K3, PLAUR, SERPINE1, and SOCS3) of these genes were significantly higher in GBM than in LGG, and comparing their expression levels to those of the proposed control genes (TBP, IPO8, and SDHA) could have the potential capability to classify patients into high- and low- risk groups, which differ significantly in the overall survival. Signatures of candidate genes were validated, by multiple microarray datasets from Gene Expression Omnibus, to increase the robustness of using these potential prognostic factors. In both the GBM and LGG cohort study, most of the patients in the high-risk group had the IDH1 wild-type gene, and those in the low-risk group had IDH1 mutations. Moreover, most of the high-risk patients with LGG possessed a 1p/19q-noncodeletion.ConclusionIn this study, we identified survival relevant genes which were shared between GBM and LGG, and those enabled to classify patients into high- and low-risk groups based on expression level analysis. Both the risk groups could be correlated with the well-known genetic variants, thus suggesting their potential prognostic value in clinical application.

[1]  J. Mesirov,et al.  The Molecular Signatures Database (MSigDB) hallmark gene set collection. , 2015, Cell systems.

[2]  A. Luini,et al.  A Golgi-based KDELR-dependent signalling pathway controls extracellular matrix degradation , 2014, Oncotarget.

[3]  N. Hashimoto,et al.  Gene Expression-Based Molecular Diagnostic System for Malignant Gliomas Is Superior to Histological Diagnosis , 2007, Clinical Cancer Research.

[4]  Baohui Liu,et al.  EFEMP2 is upregulated in gliomas and promotes glioma cell proliferation and invasion. , 2015, International journal of clinical and experimental pathology.

[5]  A. Bittner,et al.  Comparison of RNA-Seq and Microarray in Transcriptome Profiling of Activated T Cells , 2014, PloS one.

[6]  Yu Yao,et al.  Isocitrate Dehydrogenase (IDH)1/2 Mutations as Prognostic Markers in Patients With Glioblastomas , 2016, Medicine.

[7]  R. Tibshirani,et al.  Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[8]  P. Kleihues,et al.  Epidemiology and etiology of gliomas , 2005, Acta Neuropathologica.

[9]  Pedro Martínez,et al.  Identification of survival‐related genes of the phosphatidylinositol 3′‐kinase signaling pathway in glioblastoma multiforme , 2008, Cancer.

[10]  D. Cox.  Regression Models and Life-Tables , 1972 .

[11]  D. Louis.  WHO classification of tumours of the central nervous system , 2007 .

[12]  Jun Du,et al.  MICAL1 facilitates breast cancer cell proliferation via ROS‐sensitive ERK/cyclin D pathway , 2018, Journal of cellular and molecular medicine.

[13]  Hyunsoo Kim,et al.  Predicting survial by cancer pathway gene expression profiles in the TCGA , 2012, 2012 IEEE International Conference on Bioinformatics and Biomedicine Workshops.

[14]  Pei Yang,et al.  Classifying lower grade glioma cases according to whole genome gene expression , 2016, Oncotarget.

[15]  Brad T. Sherman,et al.  Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists , 2008, Nucleic acids research.

[16]  Weixing Hu,et al.  miR-330-5p suppresses glioblastoma cell proliferation and invasiveness through targeting ITGA5 , 2017, Bioscience reports.

[17]  Doron Lancet,et al.  GeneAnalytics: An Integrative Gene Set Analysis Tool for Next Generation Sequencing, RNAseq and Microarray Data , 2016, Omics : a journal of integrative biology.

[18]  W. Alexander,et al.  Suppressor of cytokine signaling 3 (SOCS3) limits damage-induced crypt hyper-proliferation and inflammation-associated tumorigenesis in the colon , 2007, Oncogene.

[19]  N. Hu,et al.  Role of chromosomal 1p/19q co-deletion on the prognosis of oligodendrogliomas: A systematic review and meta-analysis , 2016 .

[20]  R. Verhaak,et al.  Prognostically useful gene-expression profiles in acute myeloid leukemia. , 2004, The New England journal of medicine.

[21]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[22]  M. McClelland,et al.  The Accuracy of Survival Time Prediction for Patients with Glioma Is Improved by Measuring Mitotic Spindle Checkpoint Gene Expression , 2011, PloS one.

[23]  J. Mesirov,et al.  The Molecular Signatures Database Hallmark Gene Set Collection , 2015 .

[24]  Y. Nagamine,et al.  Regulation by p38 mitogen-activated protein kinase of adenylate- and uridylate-rich element-mediated urokinase-type plasminogen activator (uPA) messenger RNA stability and uPA-dependent in vitro cell invasion. , 1999, Cancer research.

[25]  Kim-Anh Do,et al.  Bayesian ensemble methods for survival prediction in gene expression data , 2011, Bioinform..

[26]  M. Cowperthwaite,et al.  Molecular Predictors of Long-Term Survival in Glioblastoma Multiforme Patients , 2016, PloS one.

[27]  H. Kretzschmar,et al.  Identification of valid endogenous control genes for determining gene expression in human glioma. , 2010, Neuro-oncology.

[28]  Hemant Ishwaran,et al.  Random Survival Forests , 2008, Wiley StatsRef: Statistics Reference Online.

[29]  Ben-Ari FuchsShani,et al.  GeneAnalytics: An Integrative Gene Set Analysis Tool for Next Generation Sequencing, RNAseq and Microarray Data , 2016 .

[30]  R. Mirimanoff,et al.  Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-year analysis of the EORTC-NCIC trial. , 2009, The Lancet. Oncology.

[31]  Y. Fujii,et al.  Gene expression signature‐based prognostic risk score in patients with glioblastoma , 2013, Cancer science.

[32]  Paola Pisani,et al.  Genetic Pathways to Glioblastoma , 2004, Cancer Research.

[33]  M. West,et al.  Gene expression profiling and genetic markers in glioblastoma survival. , 2005, Cancer research.

[34]  Steven D Chang,et al.  Magnetic resonance perfusion image features uncover an angiogenic subgroup of glioblastoma patients with poor survival and better response to antiangiogenic treatment , 2016, Neuro-oncology.

[35]  R. Verhaak,et al.  Multigene signature for predicting prognosis of patients with 1p19q co-deletion diffuse glioma , 2017, Neuro-oncology.

[36]  Udaya B. Kogalur,et al.  Random Survival Forests for R , 2007 .

[37]  Shalom Madar,et al.  Mutant p53 gain of function is interwoven into the hallmarks of cancer , 2011, The Journal of pathology.

[38]  Tom C. Freeman,et al.  Improved grading and survival prediction of human astrocytic brain tumors by artificial neural network analysis of gene expression microarray data , 2008, Molecular Cancer Therapeutics.

[39]  T. Golub,et al.  Gene expression-based classification of malignant gliomas correlates better with survival than histological classification. , 2003, Cancer research.

[40]  Yu Wang,et al.  CGCG clinical practice guidelines for the management of adult diffuse gliomas. , 2016, Cancer letters.

[41]  Caterina Giannini,et al.  Panel Review of Anaplastic Oligodendroglioma From European Organization for Research and Treatment of Cancer Trial 26951: Assessment of Consensus in Diagnosis, Influence of 1p/19q Loss, and Correlations With Outcome , 2007, Journal of neuropathology and experimental neurology.

[42]  Tej D. Azad,et al.  Magnetic resonance image features identify glioblastoma phenotypic subtypes with distinct molecular pathway activities , 2015, Science Translational Medicine.

[43]  Tong Zhou,et al.  Ion channel gene expression predicts survival in glioma patients , 2015, Scientific Reports.

[44]  L. Stewart,et al.  Chemotherapy in adult high-grade glioma: a systematic review and meta-analysis of individual patient data from 12 randomised trials , 2002, The Lancet.

[45]  F. Langmark,et al.  Progress in long-term survival in adult patients with supratentorial low-grade gliomas: a population-based study of 993 patients in whom tumors were diagnosed between 1970 and 1993. , 2003, Journal of neurosurgery.

[46]  M. Vinciguerra,et al.  ARNTL2 and SERPINE1: potential biomarkers for tumor aggressiveness in colorectal cancer , 2012, Journal of Cancer Research and Clinical Oncology.

[47]  T. Muramatsu.  Midkine, a heparin-binding cytokine with multiple roles in development, repair and diseases , 2010, Proceedings of the Japan Academy. Series B, Physical and biological sciences.

[48]  B. Fisher,et al.  Pretreatment factors predict overall survival for patients with low-grade glioma: a recursive partitioning analysis. , 1999, International journal of radiation oncology, biology, physics.

[49]  Howard Colman,et al.  IDH1 and IDH2 Mutations in Gliomas , 2013, Current Neurology and Neuroscience Reports.

[50]  G. Prager,et al.  Integrins in the Spotlight of Cancer , 2016, International journal of molecular sciences.

[51]  G. Dittmar,et al.  RNA sequencing and transcriptome arrays analyses show opposing results for alternative splicing in patient derived samples , 2017, BMC Genomics.

[52]  O. Chinot,et al.  High expression of cathepsin B and plasminogen activator inhibitor type-1 are strong predictors of survival in glioblastomas , 2009, Acta Neuropathologica.

[53]  B. Scheithauer,et al.  The 2007 WHO classification of tumours of the central nervous system , 2007, Acta Neuropathologica.

[54]  S. Horvath,et al.  Gene Expression Profiling of Gliomas Strongly Predicts Survival , 2004, Cancer Research.

[55]  Terence P. Speed,et al.  A comparison of normalization methods for high density oligonucleotide array data based on variance and bias , 2003, Bioinform..

[56]  Steven J. M. Jones,et al.  Molecular Profiling Reveals Biologically Discrete Subsets and Pathways of Progression in Diffuse Glioma , 2016, Cell.

[57]  C. Gondi,et al.  Synergistic down-regulation of urokinase plasminogen activator receptor and matrix metalloproteinase-9 in SNB19 glioblastoma cells efficiently inhibits glioma cell invasion, angiogenesis, and tumor growth. , 2003, Cancer research.

[58]  L. Girard,et al.  SHOX2 is a Potent Independent Biomarker to Predict Survival of WHO Grade II–III Diffuse Gliomas , 2016, EBioMedicine.

[59]  S. Ishii,et al.  Using gene expression profiling to identify a prognostic molecular spectrum in gliomas , 2008, Cancer science.

[60]  Y. Fujii,et al.  Identification and validation of a gene expression signature that predicts outcome in malignant glioma patients. , 2011, International journal of oncology.

[61]  Cole Trapnell,et al.  Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. , 2010, Nature biotechnology.

[62]  P. Blackshear,et al.  The RNA-binding zinc-finger protein tristetraprolin regulates AU-rich mRNAs involved in breast cancer-related processes , 2010, Oncogene.

[63]  K. Khabar.  Hallmarks of cancer and AU‐rich elements , 2016, Wiley interdisciplinary reviews. RNA.

[64]  Thomas D. Wu,et al.  Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. , 2006, Cancer cell.

[65]  Timothy C Ryken,et al.  Trends in brain cancer incidence and survival in the United States: Surveillance, Epidemiology, and End Results Program, 1973 to 2001. , 2006, Neurosurgical focus.