A Mathematical Methodology for Determining the Temporal Order of Pathway Alterations Arising during Gliomagenesis

Human cancer is caused by the accumulation of genetic alterations in cells. Of special importance are changes that occur early during malignant transformation because they may result in oncogene addiction and thus represent promising targets for therapeutic intervention. We have previously described a computational approach, called Retracing the Evolutionary Steps in Cancer (RESIC), to determine the temporal sequence of genetic alterations during tumorigenesis from cross-sectional genomic data of tumors at their fully transformed stage. Since alterations within a set of genes belonging to a particular signaling pathway may have similar or equivalent effects, we applied a pathway-based systems biology approach to the RESIC methodology. This method was used to determine whether alterations of specific pathways develop early or late during malignant transformation. When applied to primary glioblastoma (GBM) copy number data from The Cancer Genome Atlas (TCGA) project, RESIC identified a temporal order of pathway alterations consistent with the order of events in secondary GBMs. We then further subdivided the samples into the four main GBM subtypes and determined the relative contributions of each subtype to the overall results: we found that the overall ordering applied for the proneural subtype but differed for mesenchymal samples. The temporal sequence of events could not be identified for neural and classical subtypes, possibly due to a limited number of samples. Moreover, for samples of the proneural subtype, we detected two distinct temporal sequences of events: (i) RAS pathway activation was followed by TP53 inactivation and finally PI3K2 activation, and (ii) RAS activation preceded only AKT activation. This extension of the RESIC methodology provides an evolutionary mathematical approach to identify the temporal sequence of pathway changes driving tumorigenesis and may be useful in guiding the understanding of signaling rearrangements in cancer development.

[1]  E. Domany,et al.  Stem cell-related "self-renewal" signature and high epidermal growth factor receptor expression associated with resistance to concomitant chemoradiotherapy in glioblastoma. , 2008, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[2]  Jens Lagergren,et al.  New Probabilistic Network Models and Algorithms for Oncogenesis , 2006, J. Comput. Biol..

[3]  Maode Lai,et al.  -8p12-23 and +20q are predictors of subtypes and metastatic pathways in colorectal cancer: construction of tree models using comparative genomic hybridization data. , 2011, Omics : a journal of integrative biology.

[4]  C. Sander,et al.  Automated Network Analysis Identifies Core Pathways in Glioblastoma , 2010, PloS one.

[5]  M. Berger,et al.  Akt pathway activation converts anaplastic astrocytoma to glioblastoma multiforme in a human astrocyte model of glioma. , 2001, Cancer research.

[6]  W Arap,et al.  Loss of P16INK4 expression is frequent in high grade gliomas. , 1995, Cancer research.

[7]  K. Lim,et al.  Reduction in the requirement of oncogenic Ras signaling to activation of PI3K/AKT pathway during tumor maintenance. , 2005, Cancer cell.

[8]  Joshua M. Korn,et al.  Comprehensive genomic characterization defines human glioblastoma genes and core pathways , 2008, Nature.

[9]  B. Iacopetta,et al.  p53 and disease: when the guardian angel fails , 2006, Cell Death and Differentiation.

[10]  F Vogel,et al.  Spontaneous mutation in man. , 1975, Advances in human genetics.

[11]  Alejandro A Schäffer,et al.  Genetic differences detected by comparative genomic hybridization in head and neck squamous cell carcinomas from different tumor sites: construction of oncogenetic trees for tumor progression , 2002, Genes, chromosomes & cancer.

[12]  G. Parmigiani,et al.  The Consensus Coding Sequences of Human Breast and Colorectal Cancers , 2006, Science.

[13]  Ji Luo,et al.  Principles of Cancer Therapy: Oncogene and Non-oncogene Addiction , 2009, Cell.

[14]  Gary D. Bader,et al.  Pathway Commons, a web resource for biological pathway data , 2010, Nucleic Acids Res..

[15]  A. Jackson,et al.  The mutation rate and cancer. , 1998, Genetics.

[16]  Carlotta Costa,et al.  Taming the PI3K team to hold inflammation and cancer at bay. , 2008, Pharmacology & therapeutics.

[17]  A. Schäffer,et al.  Chromosome abnormalities in ovarian adenocarcinoma: III. Using breakpoint data to infer and test mathematical models for oncogenesis , 2000, Genes, chromosomes & cancer.

[18]  E. Lander,et al.  Assessing the significance of chromosomal aberrations in cancer: Methodology and application to glioma , 2007, Proceedings of the National Academy of Sciences.

[19]  Y. Yonekawa,et al.  Overexpression of the EGF Receptor and p53 Mutations are Mutually Exclusive in the Evolution of Primary and Secondary Glioblastomas , 1996 .

[20]  D. Louis,et al.  Molecular Genetic Correlates of p16, cdk4, and pRb Immunohistochemistry in Glioblastomas , 1998, Journal of neuropathology and experimental neurology.

[21]  C. Brennan,et al.  Molecular subclassification of diffuse gliomas: Seeing order in the chaos , 2011, Glia.

[22]  Thomas Lengauer,et al.  Learning multiple evolutionary pathways from cross-sectional data , 2004, J. Comput. Biol..

[23]  J. Mesirov,et al.  GenePattern 2.0 , 2006, Nature Genetics.

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

[25]  Axel Benner,et al.  Oncogenetic tree modeling of human hepatocarcinogenesis , 2012, International journal of cancer.

[26]  Yasuhiro Yonekawa,et al.  Alterations of cell cycle regulatory genes in primary (de novo) and secondary glioblastomas , 1997, Acta Neuropathologica.

[27]  Niko Beerenwinkel,et al.  Quantifying cancer progression with conjunctive Bayesian networks , 2009, Bioinform..

[28]  Niko Beerenwinkel,et al.  Construction of oncogenetic tree models reveals multiple pathways of oral cancer progression , 2009, International journal of cancer.

[29]  Kei-Hoi Cheung,et al.  BioPAX – A community standard for pathway data sharing , 2010, Nature Biotechnology.

[30]  P. Rousseeuw.  Silhouettes: a graphical aid to the interpretation and validation of cluster analysis , 1987 .

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

[32]  Camille Stephan-Otto Attolini,et al.  A mathematical framework to determine the temporal sequence of somatic genetic events in cancer , 2010, Proceedings of the National Academy of Sciences.

[33]  D. Louis,et al.  CDKN2/p16 or RB alterations occur in the majority of glioblastomas and are inversely correlated. , 1996, Cancer research.

[34]  W Arap,et al.  Cyclin-dependent kinase 6 (CDK6) amplification in human gliomas identified using two-dimensional separation of genomic DNA. , 1997, Cancer research.

[35]  J. Lupski,et al.  Genomic rearrangements and sporadic disease , 2007, Nature Genetics.

[36]  Roland Penzel,et al.  Oncogenetic tree models based on cytogenetic data: new insights into the development of epithelial tumors of the thymus. , 2005, Cancer genetics and cytogenetics.

[37]  Y Sawamura,et al.  Roles of the functional loss of p53 and other genes in astrocytoma tumorigenesis and progression. , 1999, Neuro-oncology.

[38]  I. Weinstein.  Addiction to Oncogenes--the Achilles Heal of Cancer , 2002, Science.

[39]  A. Liebman,et al.  Systematic Reviews: Synthesis of Best Evidence for Health Care Decisions , 1998, Annals of Internal Medicine.

[40]  A. W. F. Edwards,et al.  The statistical processes of evolutionary theory , 1963 .

[41]  Santosh Kesari,et al.  Molecular pathogenesis of adult brain tumors and the role of stem cells. , 2007, Neurologic clinics.

[42]  Wael M Abdel-Rahman,et al.  Genomic Instability and Carcinogenesis: An Update , 2008, Current genomics.

[43]  Jayant P. Menon,et al.  Neuronal and glioma-derived stem cell factor induces angiogenesis within the brain. , 2006, Cancer cell.

[44]  Alberto Mantovani,et al.  Cancer-related inflammation, the seventh hallmark of cancer: links to genetic instability. , 2009, Carcinogenesis.

[45]  Y. Yonekawa,et al.  Overexpression of the EGF receptor and p53 mutations are mutually exclusive in the evolution of primary and secondary glioblastomas. , 1996, Brain pathology.

[46]  Sanjay Ranka,et al.  Gene expression Distance-based clustering of CGH data , 2006 .

[47]  A. Jeffreys,et al.  Processes of copy-number change in human DNA: the dynamics of {alpha}-globin gene deletion. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[48]  Feng Jiang,et al.  Distance-Based Reconstruction of Tree Models for Oncogenesis , 2000, J. Comput. Biol..

[49]  T. Kunkel.  DNA Replication Fidelity* , 2004, Journal of Biological Chemistry.

[50]  T. Halazonetis,et al.  Genomic instability — an evolving hallmark of cancer , 2010, Nature Reviews Molecular Cell Biology.

[51]  Feng Jiang,et al.  Inferring Tree Models for Oncogenesis from Comparative Genome Hybridization Data , 1999, J. Comput. Biol..

[52]  Nader Sanai,et al.  Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma. , 2010, World neurosurgery.

[53]  J. Crow,et al.  Spontaneous mutation in man. , 1999, Mutation research.

[54]  Gunter Haroske,et al.  Comparative genomic hybridization-based oncogenetic tree model for genetic classification of breast cancer. , 2009, Analytical and quantitative cytology and histology.

[55]  Martin A. Nowak,et al.  Genetic Progression and the Waiting Time to Cancer , 2007, PLoS Comput. Biol..

[56]  K. Aldape,et al.  Formation of intracranial tumors by genetically modified human astrocytes defines four pathways critical in the development of human anaplastic astrocytoma. , 2001, Cancer research.

[57]  M Aickin,et al.  Chromosome abnormalities in ovarian adenocarcinoma: I. nonrandom chromosome abnormalities from 244 cases , 1999, Genes, chromosomes & cancer.

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

[59]  B. Vogelstein,et al.  A genetic model for colorectal tumorigenesis , 1990, Cell.

[60]  D. Hanahan,et al.  The Hallmarks of Cancer , 2000, Cell.

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

[62]  C. Brennan,et al.  Glioblastoma Subclasses Can Be Defined by Activity among Signal Transduction Pathways and Associated Genomic Alterations , 2009, PloS one.

[63]  K. Kinzler,et al.  Genetic instabilities in human cancers , 1998, Nature.

[64]  Robert A Beckman,et al.  Genetic instability is not a requirement for tumor development. , 2008, Cancer research.

[65]  D. Hanahan,et al.  Hallmarks of Cancer: The Next Generation , 2011, Cell.

[66]  Mattias Höglund,et al.  Statistical behavior of complex cancer karyotypes , 2005, Genes, chromosomes & cancer.

[67]  Eric C. Holland,et al.  Combined activation of Ras and Akt in neural progenitors induces glioblastoma formation in mice , 2000, Nature Genetics.

[68]  M. Stratton,et al.  The COSMIC (Catalogue of Somatic Mutations in Cancer) database and website , 2004, British Journal of Cancer.

[69]  L. Chin,et al.  Malignant astrocytic glioma: genetics, biology, and paths to treatment. , 2007, Genes & development.

[70]  F. Haller,et al.  An oncogenetic tree model in gastrointestinal stromal tumours (GISTs) identifies different pathways of cytogenetic evolution with prognostic implications , 2007, The Journal of pathology.

[71]  C. Sander,et al.  Functional Copy-Number Alterations in Cancer , 2008, PloS one.

[72]  Jill P. Mesirov,et al.  Consensus Clustering: A Resampling-Based Method for Class Discovery and Visualization of Gene Expression Microarray Data , 2003, Machine Learning.