Deterministic Effects Propagation Networks for reconstructing protein signaling networks from multiple interventions

BackgroundModern gene perturbation techniques, like RNA interference (RNAi), enable us to study effects of targeted interventions in cells efficiently. In combination with mRNA or protein expression data this allows to gain insights into the behavior of complex biological systems.ResultsIn this paper, we propose Deterministic Effects Propagation Networks (DEPNs) as a special Bayesian Network approach to reverse engineer signaling networks from a combination of protein expression and perturbation data. DEPNs allow to reconstruct protein networks based on combinatorial intervention effects, which are monitored via changes of the protein expression or activation over one or a few time points. Our implementation of DEPNs allows for latent network nodes (i.e. proteins without measurements) and has a built in mechanism to impute missing data. The robustness of our approach was tested on simulated data. We applied DEPNs to reconstruct the ERBB signaling network in de novo trastuzumab resistant human breast cancer cells, where protein expression was monitored on Reverse Phase Protein Arrays (RPPAs) after knockdown of network proteins using RNAi.ConclusionDEPNs offer a robust, efficient and simple approach to infer protein signaling networks from multiple interventions. The method as well as the data have been made part of the latest version of the R package "nem" available as a supplement to this paper and via the Bioconductor repository.

[1]  Ezgi O. Booth,et al.  Epistasis analysis with global transcriptional phenotypes , 2005, Nature Genetics.

[2]  Y. Yarden,et al.  The ErbB signaling network in embryogenesis and oncogenesis: signal diversification through combinatorial ligand‐receptor interactions , 1997, FEBS letters.

[3]  Nir Friedman,et al.  Inferring subnetworks from perturbed expression profiles , 2001, ISMB.

[4]  Bartek Wilczynski,et al.  Applying dynamic Bayesian networks to perturbed gene expression data , 2006, BMC Bioinformatics.

[5]  G. Mills,et al.  Reverse phase protein array: validation of a novel proteomic technology and utility for analysis of primary leukemia specimens and hematopoietic stem cells , 2006, Molecular Cancer Therapeutics.

[6]  S. Muthuswamy,et al.  Controlled activation of ErbB1/ErbB2 heterodimers promote invasion of three-dimensional organized epithelia in an ErbB1-dependent manner: implications for progression of ErbB2-overexpressing tumors. , 2006, Cancer research.

[7]  Achim Tresch,et al.  Discrimination of Direct and Indirect Interactions in a Network of Regulatory Effects , 2007, J. Comput. Biol..

[8]  Holger Fröhlich,et al.  Large scale statistical inference of signaling pathways from RNAi and microarray data , 2007, BMC Bioinformatics.

[9]  Yoshihiro Yamanishi,et al.  KEGG for linking genomes to life and the environment , 2007, Nucleic Acids Res..

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

[11]  Stefan Wiemann,et al.  Combinatorial RNAi for quantitative protein network analysis , 2007, Proceedings of the National Academy of Sciences.

[12]  David J. Spiegelhalter,et al.  Probabilistic Networks and Expert Systems - Exact Computational Methods for Bayesian Networks , 1999, Information Science and Statistics.

[13]  Nir Friedman,et al.  The Bayesian Structural EM Algorithm , 1998, UAI.

[14]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[15]  Holger Fröhlich,et al.  Estimating Large Scale Scale Signaling Networks through Nested Effects Models from Intervention Effects in Microarray Data , 2007, German Conference on Bioinformatics.

[16]  Holger Fröhlich,et al.  Analyzing gene perturbation screens with nested effects models in R and bioconductor , 2008, Bioinform..

[17]  Holger Fröhlich,et al.  Estimating large-scale signaling networks through nested effect models with intervention effects from microarray data , 2008, Bioinform..

[18]  Anthony C. Davison,et al.  Bootstrap Methods and Their Application , 1998 .

[19]  J. Pearl.  Causality: Models, Reasoning and Inference , 2000 .

[20]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[21]  Jaak Vilo,et al.  Building and analysing genome-wide gene disruption networks , 2002, ECCB.

[22]  Holger Fröhlich,et al.  Modeling ERBB receptor-regulated G1/S transition to find novel targets for de novo trastuzumab resistance , 2009, BMC Systems Biology.

[23]  K. Sachs,et al.  Causal Protein-Signaling Networks Derived from Multiparameter Single-Cell Data , 2005, Science.

[24]  C. Sander,et al.  Models from experiments: combinatorial drug perturbations of cancer cells , 2008, Molecular systems biology.

[25]  Achim Tresch,et al.  Nested effects models for learning signaling networks from perturbation data , 2009, Biometrical journal. Biometrische Zeitschrift.

[26]  A. Fire,et al.  Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans , 1998, Nature.

[27]  Alfred V. Aho,et al.  The Transitive Reduction of a Directed Graph , 1972, SIAM J. Comput..

[28]  Holger Fröhlich,et al.  A Bayesian Network View on Nested Effects Models , 2008, EURASIP J. Bioinform. Syst. Biol..

[29]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

[30]  C. Albanese,et al.  Downregulation of cyclin D1 alters cdk 4- and cdk 2-specific phosphorylation of retinoblastoma protein. , 2000, Molecular cell biology research communications : MCBRC.

[31]  Sibylle Mittnacht,et al.  Differential Phosphorylation of the Retinoblastoma Protein by G1/S Cyclin-dependent Kinases* , 1997, The Journal of Biological Chemistry.

[32]  Achim Tresch,et al.  Structure Learning in Nested Effects Models , 2007, Statistical applications in genetics and molecular biology.

[33]  K Souček,et al.  Lineage specific composition of cyclin D–CDK4/CDK6–p27 complexes reveals distinct functions of CDK4, CDK6 and individual D‐type cyclins in differentiating cells of embryonic origin , 2008, Cell proliferation.

[34]  Rainer Spang,et al.  Non-transcriptional pathway features reconstructed from secondary effects of RNA interference , 2005, Bioinform..

[35]  Olga G. Troyanskaya,et al.  Nested effects models for high-dimensional phenotyping screens , 2007, ISMB/ECCB.

[36]  Roger E Bumgarner,et al.  Integrated genomic and proteomic analyses of a systematically perturbed metabolic network. , 2001, Science.

[37]  Achim Tresch,et al.  Modeling the temporal interplay of molecular signaling and gene expression by using dynamic nested effects models , 2009, Proceedings of the National Academy of Sciences.

[38]  J. Testa,et al.  AKT activation up-regulates insulin-like growth factor I receptor expression and promotes invasiveness of human pancreatic cancer cells. , 2001, Cancer research.

[39]  David J. Spiegelhalter,et al.  Probabilistic Networks and Expert Systems , 1999, Information Science and Statistics.

[40]  Lyndsay N Harris,et al.  Efficacy and safety of trastuzumab as a single agent in first-line treatment of HER2-overexpressing metastatic breast cancer. , 2002, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[41]  D. Pe’er,et al.  Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data , 2003, Nature Genetics.

[42]  Ron Shamir,et al.  A Probabilistic Methodology for Integrating Knowledge and Experiments on Biological Networks , 2006, J. Comput. Biol..

[43]  Simon Rogers,et al.  A Bayesian regression approach to the inference of regulatory networks from gene expression data , 2005, Bioinform..

[44]  Andreas Wagner,et al.  How to reconstruct a large genetic network from n gene perturbations in fewer than n2 easy steps , 2001, Bioinform..

[45]  Joshua M. Stuart,et al.  A Factor Graph Nested Effects Model To Identify Networks from Genetic Perturbations , 2009, PLoS Comput. Biol..