Influence maximization in time bounded network identifies transcription factors regulating perturbed pathways

Motivation: To understand the dynamic nature of the biological process, it is crucial to identify perturbed pathways in an altered environment and also to infer regulators that trigger the response. Current time-series analysis methods, however, are not powerful enough to identify perturbed pathways and regulators simultaneously. Widely used methods include methods to determine gene sets such as differentially expressed genes or gene clusters and these genes sets need to be further interpreted in terms of biological pathways using other tools. Most pathway analysis methods are not designed for time series data and they do not consider gene-gene influence on the time dimension. Results: In this article, we propose a novel time-series analysis method TimeTP for determining transcription factors (TFs) regulating pathway perturbation, which narrows the focus to perturbed sub-pathways and utilizes the gene regulatory network and protein–protein interaction network to locate TFs triggering the perturbation. TimeTP first identifies perturbed sub-pathways that propagate the expression changes along the time. Starting points of the perturbed sub-pathways are mapped into the network and the most influential TFs are determined by influence maximization technique. The analysis result is visually summarized in TF-Pathway map in time clock. TimeTP was applied to PIK3CA knock-in dataset and found significant sub-pathways and their regulators relevant to the PIP3 signaling pathway. Availability and Implementation: TimeTP is implemented in Python and available at http://biohealth.snu.ac.kr/software/TimeTP/. Supplementary information: Supplementary data are available at Bioinformatics online. Contact: sunkim.bioinfo@snu.ac.kr

[1]  Rafael A. Irizarry,et al.  Bioinformatics and Computational Biology Solutions using R and Bioconductor , 2005 .

[2]  R. Dickson,et al.  Calmodulin modulates Akt activity in human breast cancer cell lines , 2009, Breast Cancer Research and Treatment.

[3]  中尾 光輝,et al.  KEGG(Kyoto Encyclopedia of Genes and Genomes)〔和文〕 (特集 ゲノム医学の現在と未来--基礎と臨床) -- (データベース) , 2000 .

[4]  Winston Haynes,et al.  Differential Expression Analysis for Pathways , 2013, PLoS Comput. Biol..

[5]  Kathleen A. Kennedy,et al.  Systems biology approaches identify ATF3 as a negative regulator of Toll-like receptor 4 , 2006, Nature.

[6]  Jason H. Moore,et al.  Pathway analysis of genomic data: concepts, methods, and prospects for future development. , 2012, Trends in genetics : TIG.

[7]  Reddanna Pallu,et al.  Role of PI3K-AKT-mTOR and Wnt Signaling Pathways in Transition of G1-S Phase of Cell Cycle in Cancer Cells , 2013, Front. Oncol..

[8]  B. Eickholt,et al.  Disruption of epithelial architecture caused by loss of PTEN or by oncogenic mutant p110α/PIK3CA but not by HER2 or mutant AKT1 , 2013, Oncogene.

[9]  Daniel Spies,et al.  Dynamics in Transcriptomics: Advancements in RNA-seq Time Course and Downstream Analysis , 2015, Computational and structural biotechnology journal.

[10]  Sheng Ding,et al.  Cooperation between both Wnt/{beta}-catenin and PTEN/PI3K/Akt signaling promotes primitive hematopoietic stem cell self-renewal and expansion. , 2011, Genes & development.

[11]  Taesung Park,et al.  Statistical tests for identifying differentially expressed genes in time-course microarray experiments , 2003, Bioinform..

[12]  I. Simon,et al.  Reconstructing dynamic regulatory maps , 2007, Molecular systems biology.

[13]  Sandrine Dudoit,et al.  More power via graph-structured tests for differential expression of gene networks , 2012, 1206.6980.

[14]  Ana Conesa,et al.  Gene expression maSigPro : a method to identify significantly differential expression profiles in time-course microarray experiments , 2006 .

[15]  Éva Tardos,et al.  Maximizing the Spread of Influence through a Social Network , 2015, Theory Comput..

[16]  Nicolas Le Novère,et al.  Perturbations of PIP3 signalling trigger a global remodelling of mRNA landscape and reveal a transcriptional feedback loop , 2015, Nucleic acids research.

[17]  Cheng-Te Li,et al.  Labeled Influence Maximization in Social Networks for Target Marketing , 2011, 2011 IEEE Third Int'l Conference on Privacy, Security, Risk and Trust and 2011 IEEE Third Int'l Conference on Social Computing.

[18]  Steven J. M. Jones,et al.  Comprehensive molecular portraits of human breast tumours , 2013 .

[19]  Jason B. Ernst,et al.  Integrating multiple evidence sources to predict transcription factor binding in the human genome. , 2010, Genome research.

[20]  Antti Honkela,et al.  Model-based method for transcription factor target identification with limited data , 2010, Proceedings of the National Academy of Sciences.

[21]  Paolo G. V. Martini,et al.  timeClip: pathway analysis for time course data without replicates , 2014, BMC Bioinformatics.

[22]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[23]  Hendrik C Korswagen,et al.  Functional Interaction Between ß-Catenin and FOXO in Oxidative Stress Signaling , 2005, Science.

[24]  Adam D Hoppe,et al.  Coordination of Fc receptor signaling regulates cellular commitment to phagocytosis , 2010, Proceedings of the National Academy of Sciences.

[25]  Pooja Mittal,et al.  A novel signaling pathway impact analysis , 2009, Bioinform..

[26]  Seungjin Choi,et al.  TEMPI: probabilistic modeling time-evolving differential PPI networks with multiPle information , 2014, Bioinform..

[27]  Sangsoo Kim,et al.  GSA-SNP: a general approach for gene set analysis of polymorphisms , 2010, Nucleic Acids Res..

[28]  Davide Heller,et al.  STRING v10: protein–protein interaction networks, integrated over the tree of life , 2014, Nucleic Acids Res..

[29]  J. Uhm,et al.  The transcriptional network for mesenchymal transformation of brain tumours , 2010 .

[30]  Ney Lemke,et al.  HTRIdb: an open-access database for experimentally verified human transcriptional regulation interactions , 2012 .

[31]  Monica Chiogna,et al.  Along signal paths: an empirical gene set approach exploiting pathway topology , 2012, Nucleic acids research.

[32]  Mohammed J. Zaki,et al.  TRICLUSTER: an effective algorithm for mining coherent clusters in 3D microarray data , 2005, SIGMOD '05.

[33]  J. Ianniello,et al.  Time delay estimation via cross-correlation in the presence of large estimation errors , 1982 .

[34]  Kyuri Jo,et al.  Time-series RNA-seq analysis package (TRAP) and its application to the analysis of rice, Oryza sativa L. ssp. Japonica, upon drought stress. , 2014, Methods.

[35]  Yuan Li,et al.  EBSeq-HMM: a Bayesian approach for identifying gene-expression changes in ordered RNA-seq experiments , 2015, Bioinform..

[36]  Naimei Tang,et al.  Akt, FoxO and regulation of apoptosis. , 2011, Biochimica et biophysica acta.

[37]  Seungjin Choi,et al.  Inference of dynamic networks using time-course data , 2014, Briefings Bioinform..

[38]  Steven J. M. Jones,et al.  Comprehensive molecular portraits of human breast tumors , 2012, Nature.

[39]  Christopher B. Burge,et al.  Methods for time series analysis of RNA-seq data with application to human Th17 cell differentiation , 2014, Bioinform..

[40]  H. Lane,et al.  ERBB receptors and cancer: the complexity of targeted inhibitors , 2005, Nature Reviews Cancer.

[41]  Pascal Nsoh,et al.  Large-scale temporal gene expression mapping of central nervous system development , 2007 .

[42]  Gordon K. Smyth,et al.  limma: Linear Models for Microarray Data , 2005 .

[43]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[44]  Joaquín Dopazo,et al.  Gene set-based analysis of polymorphisms: finding pathways or biological processes associated to traits in genome-wide association studies , 2009, Nucleic Acids Res..

[45]  R. Dickson,et al.  Calmodulin-mediated Activation of Akt Regulates Survival of c-Myc-overexpressing Mouse Mammary Carcinoma Cells* , 2004, Journal of Biological Chemistry.

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

[47]  John D. Storey,et al.  Significance analysis of time course microarray experiments. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[48]  T. Jaakkola,et al.  Comparing the continuous representation of time-series expression profiles to identify differentially expressed genes , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[49]  Mario Mellado,et al.  Role of the Pi3k Regulatory Subunit in the Control of Actin Organization and Cell Migration , 2000, The Journal of cell biology.

[50]  Alexander Schliep,et al.  Using hidden Markov models to analyze gene expression time course data , 2003, ISMB.

[51]  Atul J. Butte,et al.  Ten Years of Pathway Analysis: Current Approaches and Outstanding Challenges , 2012, PLoS Comput. Biol..

[52]  M. Neville,et al.  Epithelial cells as phagocytes: apoptotic epithelial cells are engulfed by mammary alveolar epithelial cells and repress inflammatory mediator release , 2005, Cell Death and Differentiation.

[53]  A. Frigessi,et al.  Principles and methods of integrative genomic analyses in cancer , 2014, Nature Reviews Cancer.

[54]  Paola Sebastiani,et al.  Cluster analysis of gene expression dynamics , 2002, Proceedings of the National Academy of Sciences of the United States of America.