Bayesian analysis of signaling networks governing embryonic stem cell fate decisions

MOTIVATION Signaling events that direct mouse embryonic stem (ES) cell self-renewal and differentiation are complex and accordingly difficult to understand in an integrated manner. We address this problem by adapting a Bayesian network learning algorithm to model proteomic signaling data for ES cell fate responses to external cues. Using this model we were able to characterize the signaling pathway influences as quantitative, logic-circuit type interactions. Our experimental dataset includes measurements for 28 signaling protein phosphorylation states across 16 different factorial combinations of cytokine and matrix stimuli as reported previously. RESULTS The Bayesian network modeling approach allows us to uncover previously reported signaling activities related to mouse ES cell self-renewal, such as the roles of LIF and STAT3 in maintaining undifferentiated ES cell populations. Furthermore, the network predicts novel influences such as between ERK phosphorylation and differentiation, or RAF phosphorylation and differentiated cell proliferation. Visualization of the influences detected by the Bayesian network provides intuition about the underlying physiology of the signaling pathways. We demonstrate that the Bayesian networks can capture the linear, nonlinear and multistate logic interactions that connect extracellular cues, intracellular signals and consequent cell functional responses.

[1]  Wendy A. Prudhomme,et al.  Cell population dynamics model for deconvolution of murine embryonic stem cell self‐renewal and differentiation responses to cytokines and extracellular matrix , 2004, Biotechnology and bioengineering.

[2]  Douglas A Lauffenburger,et al.  Multivariate proteomic analysis of murine embryonic stem cell self-renewal versus differentiation signaling. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[3]  J. Frost,et al.  MEKK1 phosphorylates MEK1 and MEK2 but does not cause activation of mitogen-activated protein kinase. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[4]  Tommi S. Jaakkola,et al.  Combining Location and Expression Data for Principled Discovery of Genetic Regulatory Network Models , 2001, Pacific Symposium on Biocomputing.

[5]  John G. Albeck,et al.  Cue-Signal-Response Analysis of TNF-Induced Apoptosis by Partial Least Squares Regression of Dynamic Multivariate Data , 2004, J. Comput. Biol..

[6]  P. D’Eustachio,et al.  Essential role of STAT3 for embryonic stem cell pluripotency. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Nir Friedman,et al.  Data Analysis with Bayesian Networks: A Bootstrap Approach , 1999, UAI.

[8]  Daphne Koller,et al.  Active Learning for Structure in Bayesian Networks , 2001, IJCAI.

[9]  Donald Metcalf,et al.  Myeloid leukaemia inhibitory factor maintains the developmental potential of embryonic stem cells , 1988, Nature.

[10]  Anastasios Bezerianos,et al.  Gene expression data analysis with a dynamically extended self-organized map that exploits class information , 2002, Bioinform..

[11]  R. Eckert,et al.  Green Tea Polyphenol Stimulates a Ras, MEKK1, MEK3, and p38 Cascade to Increase Activator Protein 1 Factor-dependent Involucrin Gene Expression in Normal Human Keratinocytes* , 2002, The Journal of Biological Chemistry.

[12]  E. Krebs,et al.  Phosphorylation of the human leukemia inhibitory factor (LIF) receptor by mitogen-activated protein kinase and the regulation of LIF receptor function by heterologous receptor activation. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[13]  J. Collins,et al.  Inferring Genetic Networks and Identifying Compound Mode of Action via Expression Profiling , 2003, Science.

[14]  Gregory F. Cooper,et al.  Causal Discovery from a Mixture of Experimental and Observational Data , 1999, UAI.

[15]  M. White,et al.  ERK5 and ERK2 cooperate to regulate NF-kappaB and cell transformation. , 2001, The Journal of biological chemistry.

[16]  T. Jaakkola,et al.  Bayesian Network Approach to Cell Signaling Pathway Modeling , 2002, Science's STKE.

[17]  G. Martin,et al.  The Role of Basement Membranes in Vascular Development , 1990, Annals of the New York Academy of Sciences.

[18]  P. Lazarovici,et al.  Signaling Pathways for PC12 Cell Differentiation: Making the Right Connections , 2002, Science.

[19]  K. Kawauchi,et al.  Constitutive Activation of Extracellular Signal-Regulated Kinase and p38 Mitogen-Activated Protein Kinase in B-Cell Lymphoproliferative Disorders , 2003, International journal of hematology.

[20]  Michael I. Jordan Learning in Graphical Models , 1999, NATO ASI Series.

[21]  David Maxwell Chickering,et al.  Learning Bayesian Networks is NP-Complete , 2016, AISTATS.

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

[23]  Channing J Der,et al.  The dark side of Ras: regulation of apoptosis , 2003, Oncogene.

[24]  John K. Heath,et al.  Inhibition of pluripotential embryonic stem cell differentiation by purified polypeptides , 1988, Nature.

[25]  Zoubin Ghahramani,et al.  An Introduction to Hidden Markov Models and Bayesian Networks , 2001, Int. J. Pattern Recognit. Artif. Intell..

[26]  Judea Pearl,et al.  Probabilistic reasoning in intelligent systems - networks of plausible inference , 1991, Morgan Kaufmann series in representation and reasoning.

[27]  G. Merlino,et al.  Hepatocyte growth factor/scatter factor activates proliferation in melanoma cells through p38 MAPK, ATF-2 and cyclin D1 , 2002, Oncogene.

[28]  Nir Friedman,et al.  Probabilistic models for identifying regulation networks , 2003, ECCB.

[29]  M. White,et al.  ERK5 and ERK2 Cooperate to Regulate NF-κB and Cell Transformation* , 2001, The Journal of Biological Chemistry.

[30]  K L Morand,et al.  Advances in high-throughput mass spectrometry. , 2001, Current opinion in drug discovery & development.

[31]  John R. Koza,et al.  Hidden Order: How Adaptation Builds Complexity. , 1995, Artificial Life.

[32]  Bruce D'Ambrosio,et al.  Real-Time Inference with Large-Scale Temporal Bayes Nets , 2002, UAI.

[33]  S. Birken,et al.  Preparation and analysis of the common urinary forms of human chorionic gonadotropin. , 2000, Methods.

[34]  Gregory F. Cooper,et al.  A Bayesian method for the induction of probabilistic networks from data , 1992, Machine Learning.

[35]  Jason A. Papin,et al.  Genome-scale microbial in silico models: the constraints-based approach. , 2003, Trends in biotechnology.

[36]  K. Kishi,et al.  Overexpression of wild-type Akt1 promoted insulin-stimulated p70S6 kinase (p70S6K) activity and affected GSK3 beta regulation, but did not promote insulin-stimulated GLUT4 translocation or glucose transport in L6 myotubes. , 2000, The journal of medical investigation : JMI.

[37]  D. Lauffenburger,et al.  Leukemia inhibitory factor (LIF) concentration modulates embryonic stem cell self-renewal and differentiation independently of proliferation. , 2000, Biotechnology and bioengineering.

[38]  P. Uetz,et al.  High-throughput screening for protein-protein interactions using two-hybrid assay. , 2000, Methods in enzymology.

[39]  M. Katsuki,et al.  STAT3 activation is sufficient to maintain an undifferentiated state of mouse embryonic stem cells , 1999, The EMBO journal.

[40]  James C Maliakal,et al.  Quantitative high throughput endothelial cell migration and invasion assay system. , 2002, Methods in enzymology.

[41]  Michael A. Beer,et al.  Predicting Gene Expression from Sequence , 2004, Cell.

[42]  Doug Fisher,et al.  Learning from Data: Artificial Intelligence and Statistics V , 1996 .

[43]  G. Johnson,et al.  Mitogen-Activated Protein Kinase Pathways Mediated by ERK, JNK, and p38 Protein Kinases , 2002, Science.

[44]  Richard E. Neapolitan,et al.  Learning Bayesian networks , 2007, KDD '07.

[45]  Pavel Kovarik,et al.  Serine phosphorylation of STATs , 2000, Oncogene.

[46]  L. Peters,et al.  The mouse adducin gene family: alternative splicing and chromosomal localization , 2000, Mammalian Genome.

[47]  Y. Fujio,et al.  Activation of Phosphatidylinositol 3-Kinase through Glycoprotein 130 Induces Protein Kinase B and p70 S6 Kinase Phosphorylation in Cardiac Myocytes* , 1998, The Journal of Biological Chemistry.

[48]  R. Albert,et al.  The large-scale organization of metabolic networks , 2000, Nature.

[49]  A. Smith,et al.  Self-renewal of pluripotent embryonic stem cells is mediated via activation of STAT3. , 1998, Genes & development.

[50]  A. Dunn,et al.  gp130-mediated Signal Transduction in Embryonic Stem Cells Involves Activation of Jak and Ras/Mitogen-activated Protein Kinase Pathways* , 1996, The Journal of Biological Chemistry.

[51]  R. de Waal Malefyt,et al.  Flow cytometric analysis of immunoprecipitates: high-throughput analysis of protein phosphorylation and protein-protein interactions. , 2000, Cytometry.

[52]  Daphne Koller,et al.  Nonuniform Dynamic Discretization in Hybrid Networks , 1997, UAI.

[53]  Brian N. Johnson,et al.  An integrated nanoliter DNA analysis device. , 1998, Science.

[54]  E. Krebs,et al.  p38 MAPK is required for CD40-induced gene expression and proliferation in B lymphocytes. , 1998, Journal of immunology.