Causal integration of multi‐omics data with prior knowledge to generate mechanistic hypotheses

Multi-omics datasets can provide molecular insights beyond the sum of individual omics. Diverse tools have been recently developed to integrate such datasets, but there are limited strategies to systematically extract mechanistic hypotheses from them. Here, we present COSMOS (Causal Oriented Search of Multi-Omics Space), a method that integrates phosphoproteomics, transcriptomics, and metabolics datasets. COSMOS combines extensive prior knowledge of signaling, metabolic, and gene regulatory networks with computational methods to estimate activities of transcription factors and kinases as well as network-level causal reasoning. COSMOS provides mechanistic hypotheses for experimental observations across multi-omics datasets. We applied COSMOS to a dataset comprising transcriptomics, phosphoproteomics, and metabolomics data from healthy and cancerous tissue from nine renal cell carcinoma patients. We used COSMOS to generate novel hypotheses such as the impact of Androgen Receptor on nucleoside metabolism and the influence of the JAK-STAT pathway on propionyl coenzyme A production. We expect that our freely available method will be broadly useful to extract mechanistic insights from multi-omics studies.

[1]  Nils Blüthgen,et al.  Classification of gene signatures for their information value and functional redundancy , 2017, npj Systems Biology and Applications.

[2]  E. Barillot,et al.  Classification of gene signatures for their information value and functional redundancy , 2017, bioRxiv.

[3]  Martin Vingron,et al.  Variance stabilization applied to microarray data calibration and to the quantification of differential expression , 2002, ISMB.

[4]  M. Ivan,et al.  HIFα Targeted for VHL-Mediated Destruction by Proline Hydroxylation: Implications for O2 Sensing , 2001, Science.

[5]  R. Rosell,et al.  AURKB as a target in non-small cell lung cancer with acquired resistance to anti-EGFR therapy , 2019, Nature Communications.

[6]  E. Gottlieb,et al.  Analysis of Cell Metabolism Using LC-MS and Isotope Tracers. , 2015, Methods in enzymology.

[7]  V. Frouin,et al.  Variable selection for generalized canonical correlation analysis. , 2014, Biostatistics.

[8]  Laura L. Elo,et al.  A systematic evaluation of normalization methods in quantitative label-free proteomics , 2016, Briefings Bioinform..

[9]  Daniel C. Zielinski,et al.  Recon3D enables a three-dimensional view of gene variation in human metabolism , 2018 .

[10]  C. Wykoff,et al.  The tumour suppressor protein VHL targets hypoxia-inducible factors for oxygen-dependent proteolysis , 1999, Nature.

[11]  C. Lindskog,et al.  A pathology atlas of the human cancer transcriptome , 2017, Science.

[12]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[13]  Chawnshang Chang,et al.  Androgen receptor (AR) promotes clear cell renal cell carcinoma (ccRCC) migration and invasion via altering the circHIAT1/miR-195-5p/29a-3p/29c-3p/CDC42 signals. , 2017, Cancer letters.

[14]  M. Falasca,et al.  Targeting PDK1 for Chemosensitization of Cancer Cells , 2017, Cancers.

[15]  Mariano J. Alvarez,et al.  Network-based inference of protein activity helps functionalize the genetic landscape of cancer , 2016, Nature Genetics.

[16]  Julio Saez-Rodriguez,et al.  OmniPath: guidelines and gateway for literature-curated signaling pathway resources , 2016, Nature Methods.

[17]  Kim-Anh Lê Cao,et al.  DIABLO: an integrative approach for identifying key molecular drivers from multi-omics assays , 2019, Bioinform..

[18]  C. Olson,et al.  Peer review of the biomedical literature. , 1990, The American journal of emergency medicine.

[19]  M. Sayed-Ahmed Role of carnitine in cancer chemotherapy-induced multiple organ toxicity. , 2010, Saudi pharmaceutical journal : SPJ : the official publication of the Saudi Pharmaceutical Society.

[20]  A. Viale,et al.  Epigenetic expansion of VHL-HIF signal output drives multi-organ metastasis in renal cancer , 2012, Nature Medicine.

[21]  S. Leung,et al.  STAT1 drives tumor progression in serous papillary endometrial cancer. , 2014, Cancer research.

[22]  S. Vanharanta,et al.  NF-κB-Dependent Lymphoid Enhancer Co-option Promotes Renal Carcinoma Metastasis. , 2018, Cancer discovery.

[23]  Emanuel J. V. Gonçalves,et al.  A Landscape of Pharmacogenomic Interactions in Cancer , 2016, Cell.

[24]  Angela N. Brooks,et al.  A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles , 2017, Cell.

[25]  Emanuel J. V. Gonçalves,et al.  Fumarate is an epigenetic modifier that elicits epithelial-to-mesenchymal transition , 2016, Nature.

[26]  J. Marioni,et al.  Multi‐Omics Factor Analysis—a framework for unsupervised integration of multi‐omics data sets , 2018, Molecular systems biology.

[27]  Michael I. Wilson,et al.  Targeting of HIF-α to the von Hippel-Lindau Ubiquitylation Complex by O2-Regulated Prolyl Hydroxylation , 2001, Science.

[28]  Nuno A. Fonseca,et al.  Transcription Factor Activities Enhance Markers of Drug Sensitivity in Cancer. , 2018, Cancer research.

[29]  David Haussler,et al.  Discovering causal pathways linking genomic events to transcriptional states using Tied Diffusion Through Interacting Events (TieDIE) , 2013, Bioinform..

[30]  Qiu-ming He,et al.  Low expression of PDK1 inhibits renal cell carcinoma cell proliferation, migration, invasion and epithelial mesenchymal transition through inhibition of the PI3K-PDK1-Akt pathway. , 2019, Cellular signalling.

[31]  Julio Saez-Rodriguez,et al.  Footprint-based functional analysis of multiomic data , 2019, Current opinion in systems biology.

[32]  Mikhail Chernov,et al.  Regulation of c‐myc expression by IFN‐γ through Stat1‐dependent and ‐independent pathways , 2000 .

[33]  P. Vinod,et al.  Network-based metabolic characterization of renal cell carcinoma , 2020, Scientific Reports.

[34]  Hyungwon Choi,et al.  Exploiting Interdata Relationships in Next-generation Proteomics Analysis* , 2019, Molecular & Cellular Proteomics.

[35]  J. Hengstler,et al.  Protein kinase C eta is associated with progression of renal cell carcinoma (RCC). , 2003, Anticancer research.

[36]  Mehmet Koyutürk,et al.  The KSEA App: a web‐based tool for kinase activity inference from quantitative phosphoproteomics , 2017, Bioinform..

[37]  D. Fabbro,et al.  Protein kinase C in human renal cell carcinomas: role in invasion and differential isoenzyme expression , 2000, British Journal of Cancer.

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

[39]  J. Blay,et al.  Adenosine upregulates CXCR4 and enhances the proliferative and migratory responses of human carcinoma cells to CXCL12/SDF‐1α , 2006, International journal of cancer.

[40]  P. O’Farrell Faculty Opinions recommendation of Targeting of HIF-alpha to the von Hippel-Lindau ubiquitylation complex by O2-regulated prolyl hydroxylation. , 2001 .

[41]  Damian Szklarczyk,et al.  STITCH 5: augmenting protein–chemical interaction networks with tissue and affinity data , 2015, Nucleic Acids Res..

[42]  Lu Zhang,et al.  Adenosine Stimulates the Migration of Human Endothelial Progenitor Cells. Role ofCXCR4 and MicroRNA-150 , 2013, PloS one.

[43]  Yutaka Suzuki,et al.  Trans-omic Analysis Reveals Selective Responses to Induced and Basal Insulin across Signaling, Transcriptional, and Metabolic Networks , 2018, iScience.

[44]  G. Ryffel,et al.  HNF4α reduces proliferation of kidney cells and affects genes deregulated in renal cell carcinoma , 2005, Oncogene.

[45]  Samuel H Payne,et al.  Extracting Pathway-level Signatures from Proteogenomic Data in Breast Cancer Using Independent Component Analysis* , 2019, Molecular & Cellular Proteomics.

[46]  J. Leppert,et al.  A Protective Role for Androgen Receptor in Clear Cell Renal Cell Carcinoma Based on Mining TCGA Data , 2016, PloS one.

[47]  Christian H. Holland,et al.  Robustness and applicability of transcription factor and pathway analysis tools on single-cell RNA-seq data , 2020, Genome Biology.

[48]  Joanna Polanska,et al.  Ranking metrics in gene set enrichment analysis: do they matter? , 2017, BMC Bioinformatics.

[49]  Evan O. Paull,et al.  Phosphoproteome Integration Reveals Patient-Specific Networks in Prostate Cancer , 2016, Cell.

[50]  Bangbei Wan,et al.  AURKB: a promising biomarker in clear cell renal cell carcinoma , 2019, PeerJ.

[51]  Panuwat Trairatphisan,et al.  From expression footprints to causal pathways: contextualizing large signaling networks with CARNIVAL , 2019, npj Systems Biology and Applications.

[52]  Junnian Zheng,et al.  Aurora kinases: novel therapy targets in cancers , 2017, Oncotarget.

[53]  S. Yeh,et al.  The expression and evaluation of androgen receptor in human renal cell carcinoma. , 2014, Urology.

[54]  T. Giordano,et al.  Follicular cell thyroid neoplasia: insights from genomics and The Cancer Genome Atlas research network , 2016, Current opinion in oncology.

[55]  Bertram Klinger,et al.  Perturbation-response genes reveal signaling footprints in cancer gene expression , 2016 .