An integrative somatic mutation analysis to identify pathways linked with survival outcomes across 19 cancer types

MOTIVATION Identification of altered pathways that are clinically relevant across human cancers is a key challenge in cancer genomics. Precise identification and understanding of these altered pathways may provide novel insights into patient stratification, therapeutic strategies and the development of new drugs. However, a challenge remains in accurately identifying pathways altered by somatic mutations across human cancers, due to the diverse mutation spectrum. We developed an innovative approach to integrate somatic mutation data with gene networks and pathways, in order to identify pathways altered by somatic mutations across cancers. RESULTS We applied our approach to The Cancer Genome Atlas (TCGA) dataset of somatic mutations in 4790 cancer patients with 19 different types of tumors. Our analysis identified cancer-type-specific altered pathways enriched with known cancer-relevant genes and targets of currently available drugs. To investigate the clinical significance of these altered pathways, we performed consensus clustering for patient stratification using member genes in the altered pathways coupled with gene expression datasets from 4870 patients from TCGA, and multiple independent cohorts confirmed that the altered pathways could be used to stratify patients into subgroups with significantly different clinical outcomes. Of particular significance, certain patient subpopulations with poor prognosis were identified because they had specific altered pathways for which there are available targeted therapies. These findings could be used to tailor and intensify therapy in these patients, for whom current therapy is suboptimal. AVAILABILITY AND IMPLEMENTATION The code is available at: http://www.taehyunlab.org CONTACT jhcheong@yuhs.ac or taehyun.hwang@utsouthwestern.edu or taehyun.cs@gmail.com SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  D. Duan,et al.  Inhibition of transcription elongation by the VHL tumor suppressor protein , 1995, Science.

[2]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[3]  H. Sebastian Seung,et al.  Algorithms for Non-negative Matrix Factorization , 2000, NIPS.

[4]  L. Poellinger,et al.  Mechanism of regulation of the hypoxia‐inducible factor‐1α by the von Hippel‐Lindau tumor suppressor protein , 2000, The EMBO journal.

[5]  M. Ivan,et al.  Ubiquitination of hypoxia-inducible factor requires direct binding to the β-domain of the von Hippel–Lindau protein , 2000, Nature Cell Biology.

[6]  S. A. Strel’tsov,et al.  [Interaction of topotecan--a DNA topoisomerase I inhibitor--with dual-stranded polydeoxyribonucleotides. II. Formation of a complex containing several DNA molecules in the presence of topotecan]. , 2001, Molekuliarnaia biologiia.

[7]  Strel'tsov Sa,et al.  [Interaction of topotecan--a DNA topoisomerase I inhibitor--with dual-stranded polydeoxyribonucleotides. II. Formation of a complex containing several DNA molecules in the presence of topotecan]. , 2001 .

[8]  J. Fenton,et al.  Effect of glucosamine on interleukin-1-conditioned articular cartilage. , 2010, Equine veterinary journal. Supplement.

[9]  R. Körfer,et al.  Cardiac remodelling in end stage heart failure: upregulation of matrix metalloproteinase (MMP) irrespective of the underlying disease, and evidence for a direct inhibitory effect of ACE inhibitors on MMP , 2002, Heart.

[10]  S. Jimenez,et al.  Glucosamine sulfate modulates the levels of aggrecan and matrix metalloproteinase-3 synthesized by cultured human osteoarthritis articular chondrocytes. , 2003, Osteoarthritis and cartilage.

[11]  A. Harris,et al.  Combination antiangiogenesis therapy with marimastat, captopril and fragmin in patients with advanced cancer , 2004, British Journal of Cancer.

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

[13]  R. Williams,et al.  Inhibition of matrix metalloproteinase activity and growth of gastric adenocarcinoma cells by an angiotensin converting enzyme inhibitor in in vitro and murine models. , 2005, European journal of surgical oncology : the journal of the European Society of Surgical Oncology and the British Association of Surgical Oncology.

[14]  N. Rajapakse,et al.  Carboxy derivatized glucosamine is a potent inhibitor of matrix metalloproteinase-9 in HT1080 cells. , 2006, Bioorganic & medicinal chemistry letters.

[15]  John P. Overington,et al.  How many drug targets are there? , 2006, Nature Reviews Drug Discovery.

[16]  P. Imming,et al.  Drugs, their targets and the nature and number of drug targets , 2006, Nature Reviews Drug Discovery.

[17]  M. Chou,et al.  Glucosamine sulfate suppresses the expressions of urokinase plasminogen activator and inhibitor and gelatinases during the early stage of osteoarthritis. , 2006, Clinica chimica acta; international journal of clinical chemistry.

[18]  Sakae Tanaka,et al.  Tissue-type Plasminogen Activator Acts as a Cytokine That Triggers Intracellular Signal Transduction and Induces Matrix Metalloproteinase-9 Gene Expression* , 2006, Journal of Biological Chemistry.

[19]  N. Rajapakse,et al.  Sulfated glucosamine inhibits MMP-2 and MMP-9 expressions in human fibrosarcoma cells. , 2007, Bioorganic & medicinal chemistry.

[20]  P. Berna,et al.  Metalloelastase (MMP-12) induced inflammatory response in mice airways: effects of dexamethasone, rolipram and marimastat. , 2007, European journal of pharmacology.

[21]  M. Okada,et al.  Captopril attenuates matrix metalloproteinase-2 and -9 in monocrotaline-induced right ventricular hypertrophy in rats. , 2008, Journal of pharmacological sciences.

[22]  S. Takai,et al.  Inhibitory profiles of captopril on matrix metalloproteinase-9 activity. , 2008, European journal of pharmacology.

[23]  J. Netterville,et al.  Phase II trial of irinotecan plus cisplatin in patients with recurrent or metastatic squamous carcinoma of the head and neck , 2008, Cancer.

[24]  F. Waldman,et al.  Improved Identification of von Hippel-Lindau Gene Alterations in Clear Cell Renal Tumors , 2008, Clinical Cancer Research.

[25]  Heesang Song,et al.  Low-density lipoprotein receptor-related protein 1 promotes cancer cell migration and invasion by inducing the expression of matrix metalloproteinases 2 and 9. , 2009, Cancer research.

[26]  Sandhya Rani,et al.  Human Protein Reference Database—2009 update , 2008, Nucleic Acids Res..

[27]  Joel Dudley,et al.  Network-Based Elucidation of Human Disease Similarities Reveals Common Functional Modules Enriched for Pluripotent Drug Targets , 2010, PLoS Comput. Biol..

[28]  David Haussler,et al.  Inference of patient-specific pathway activities from multi-dimensional cancer genomics data using PARADIGM , 2010, Bioinform..

[29]  B. Langlois,et al.  LRP-1 Promotes Cancer Cell Invasion by Supporting ERK and Inhibiting JNK Signaling Pathways , 2010, PloS one.

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

[31]  C. Gialeli,et al.  Roles of matrix metalloproteinases in cancer progression and their pharmacological targeting , 2011, The FEBS journal.

[32]  M. Daly,et al.  Proteins Encoded in Genomic Regions Associated with Immune-Mediated Disease Physically Interact and Suggest Underlying Biology , 2011, PLoS genetics.

[33]  M. Rask-Andersen,et al.  Trends in the exploitation of novel drug targets , 2011, Nature Reviews Drug Discovery.

[34]  J. Baselga Targeting the phosphoinositide-3 (PI3) kinase pathway in breast cancer. , 2011, The oncologist.

[35]  Eli Upfal,et al.  Algorithms for Detecting Significantly Mutated Pathways in Cancer , 2010, RECOMB.

[36]  Juan Liu,et al.  A novel computational framework for simultaneous integration of multiple types of genomic data to identify microRNA-gene regulatory modules , 2011, Bioinform..

[37]  Georgios B. Giannakis,et al.  Sparse robust matrix tri-factorization with application to cancer genomics , 2012, 2012 3rd International Workshop on Cognitive Information Processing (CIP).

[38]  Matthew B. Callaway,et al.  MuSiC: Identifying mutational significance in cancer genomes , 2012, Genome research.

[39]  Vipin Kumar,et al.  Large-scale integrative network-based analysis identifies common pathways disrupted by copy number alterations across cancers , 2013, BMC Genomics.

[40]  N. Grishin,et al.  BAP1 loss defines a new class of renal cell carcinoma , 2012, Nature Genetics.

[41]  Vipin Kumar,et al.  Co-clustering phenome–genome for phenotype classification and disease gene discovery , 2012, Nucleic acids research.

[42]  Benjamin E. Gross,et al.  The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. , 2012, Cancer discovery.

[43]  Tamara G. Kolda,et al.  On Tensors, Sparsity, and Nonnegative Factorizations , 2011, SIAM J. Matrix Anal. Appl..

[44]  Steven J. M. Jones,et al.  Comprehensive molecular characterization of clear cell renal cell carcinoma , 2013, Nature.

[45]  Andrew M. Gross,et al.  Network-based stratification of tumor mutations , 2013, Nature Methods.

[46]  S. Vandenberg,et al.  Myeloid cell receptor LRP1/CD91 regulates monocyte recruitment and angiogenesis in tumors. , 2013, Cancer research.

[47]  Baolin Wu,et al.  Network-based Survival Analysis Reveals Subnetwork Signatures for Predicting Outcomes of Ovarian Cancer Treatment , 2013, PLoS Comput. Biol..

[48]  A. Sato,et al.  Inhibition of MMP-9 using a pyrrole-imidazole polyamide reduces cell invasion in renal cell carcinoma. , 2013, International journal of oncology.

[49]  Benjamin J. Raphael,et al.  Mutational landscape and significance across 12 major cancer types , 2013, Nature.

[50]  E. Lander,et al.  Lessons from the Cancer Genome , 2013, Cell.

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

[52]  Steven A. Roberts,et al.  Mutational heterogeneity in cancer and the search for new cancer-associated genes , 2013 .

[53]  H. Aburatani,et al.  Integrated molecular analysis of clear-cell renal cell carcinoma , 2013, Nature Genetics.

[54]  Gary D Bader,et al.  Comprehensive identification of mutational cancer driver genes across 12 tumor types , 2013, Scientific Reports.

[55]  J. Brugarolas,et al.  Simultaneous isolation of high-quality DNA, RNA, miRNA and proteins from tissues for genomic applications , 2013, Nature Protocols.

[56]  J. Valcárcel,et al.  Synonymous Mutations Frequently Act as Driver Mutations in Human Cancers , 2014, Cell.

[57]  S. Gabriel,et al.  Discovery and saturation analysis of cancer genes across 21 tumor types , 2014, Nature.

[58]  C. Leslie,et al.  Linking signaling pathways to transcriptional programs in breast cancer , 2014, Genome research.

[59]  Mingming Jia,et al.  COSMIC: exploring the world's knowledge of somatic mutations in human cancer , 2014, Nucleic Acids Res..