Using Answer Set Programming to Integrate RNA Expression with Signalling Pathway Information to Infer How Mutations Affect Ageing

A challenge of systems biology is to integrate incomplete knowledge on pathways with existing experimental data sets and relate these to measured phenotypes. Research on ageing often generates such incomplete data, creating difficulties in integrating RNA expression with information about biological processes and the phenotypes of ageing, including longevity. Here, we develop a logic-based method that employs Answer Set Programming, and use it to infer signalling effects of genetic perturbations, based on a model of the insulin signalling pathway. We apply our method to RNA expression data from Drosophila mutants in the insulin pathway that alter lifespan, in a foxo dependent fashion. We use this information to deduce how the pathway influences lifespan in the mutant animals. We also develop a method for inferring the largest common sub-paths within each of our signalling predictions. Our comparisons reveal consistent homeostatic mechanisms across both long- and short-lived mutants. The transcriptional changes observed in each mutation usually provide negative feedback to signalling predicted for that mutation. We also identify an S6K-mediated feedback in two long-lived mutants that suggests a crosstalk between these pathways in mutants of the insulin pathway, in vivo. By formulating the problem as a logic-based theory in a qualitative fashion, we are able to use the efficient search facilities of Answer Set Programming, allowing us to explore larger pathways, combine molecular changes with pathways and phenotype and infer effects on signalling in in vivo, whole-organism, mutants, where direct signalling stimulation assays are difficult to perform. Our methods are available in the web-service NetEffects: http://www.ebi.ac.uk/thornton-srv/software/NetEffects.

[1]  Marius Thomas Lindauer,et al.  Potassco: The Potsdam Answer Set Solving Collection , 2011, AI Commun..

[2]  Luis Serrano,et al.  Correlation of mRNA and protein in complex biological samples , 2009, FEBS letters.

[3]  Martin Kircher,et al.  Deep proteome and transcriptome mapping of a human cancer cell line , 2011, Molecular systems biology.

[4]  Cynthia Kenyon,et al.  The first long-lived mutants: discovery of the insulin/IGF-1 pathway for ageing , 2011, Philosophical Transactions of the Royal Society B: Biological Sciences.

[5]  Martin Gebser,et al.  The BioASP Library: ASP Solutions for Systems Biology , 2010, 2010 22nd IEEE International Conference on Tools with Artificial Intelligence.

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

[7]  Irene Papatheodorou,et al.  Genome-wide dFOXO targets and topology of the transcriptomic response to stress and insulin signalling , 2011, Molecular systems biology.

[8]  Susumu Goto,et al.  KEGG for integration and interpretation of large-scale molecular data sets , 2011, Nucleic Acids Res..

[9]  Janet M Thornton,et al.  Evolutionary conservation of regulated longevity assurance mechanisms , 2007, Genome Biology.

[10]  R. Tjian,et al.  Control of cell number by Drosophila FOXO: downstream and feedback regulation of the insulin receptor pathway. , 2003, Genes & development.

[11]  Aurélien Naldi,et al.  Logical modelling of regulatory networks with GINsim 2.3 , 2009, Biosyst..

[12]  Peter K. Sorger,et al.  Logic-Based Models for the Analysis of Cell Signaling Networks† , 2010, Biochemistry.

[13]  L. Partridge,et al.  Lifespan extension by increased expression of the Drosophila homologue of the IGFBP7 tumour suppressor , 2011, Aging cell.

[14]  Janet M. Thornton,et al.  Ribosomal Protein S6 Kinase 1 Signaling Regulates Mammalian Life Span , 2009, Science.

[15]  N. Perrimon,et al.  Dynamic Switch of Negative Feedback Regulation in Drosophila Akt–TOR Signaling , 2010, PLoS genetics.

[16]  G. Thomas,et al.  Nutrient overload, insulin resistance, and ribosomal protein S6 kinase 1, S6K1. , 2006, Cell metabolism.

[17]  Daryl P. Shanley,et al.  Modelling the Response of FOXO Transcription Factors to Multiple Post-Translational Modifications Made by Ageing-Related Signalling Pathways , 2010, PloS one.

[18]  Janet M. Thornton,et al.  Regulation of Lifespan, Metabolism, and Stress Responses by the Drosophila SH2B Protein, Lnk , 2010, PLoS genetics.

[19]  E. Hafen,et al.  The Drosophila SH2B Family Adaptor Lnk Acts in Parallel to Chico in the Insulin Signaling Pathway , 2009, PLoS genetics.

[20]  Chitta Baral,et al.  Knowledge Representation, Reasoning and Declarative Problem Solving , 2003 .

[21]  A. Teleman Molecular mechanisms of metabolic regulation by insulin in Drosophila. , 2009, The Biochemical journal.

[22]  Matthias Blüher,et al.  Extended Longevity in Mice Lacking the Insulin Receptor in Adipose Tissue , 2003, Science.

[23]  Robert Tjian,et al.  Transcriptional feedback control of insulin receptor by dFOXO/FOXO1. , 2005, Genes & development.

[24]  D. Lauffenburger,et al.  Discrete logic modelling as a means to link protein signalling networks with functional analysis of mammalian signal transduction , 2009, Molecular systems biology.

[25]  Torsten Schaub,et al.  Modeling Biological Networks by Action Languages via Answer Set Programming , 2008, Constraints.

[26]  Chris Cornelis,et al.  Modelling gene and protein regulatory networks with Answer Set Programming , 2011, Int. J. Data Min. Bioinform..

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

[28]  Chris Sander,et al.  Discovering modulators of gene expression , 2010, Nucleic acids research.

[29]  J. Thornton,et al.  Computational biology for ageing , 2011, Philosophical Transactions of the Royal Society B: Biological Sciences.

[30]  A. Newton,et al.  PHLPP: a phosphatase that directly dephosphorylates Akt, promotes apoptosis, and suppresses tumor growth. , 2005, Molecular cell.

[31]  M. Tatar,et al.  A Mutant Drosophila Insulin Receptor Homolog That Extends Life-Span and Impairs Neuroendocrine Function , 2001, Science.

[32]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[33]  E. Hafen,et al.  Susi, a negative regulator of Drosophila PI3-kinase. , 2005, Developmental cell.

[34]  T. Speed,et al.  Summaries of Affymetrix GeneChip probe level data. , 2003, Nucleic acids research.

[35]  Linda Partridge,et al.  Evidence for lifespan extension and delayed age–related biomarkers in insulin receptor substrate 1 null mice , 2008, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[36]  M. Nielsen,et al.  14‐3‐3ɛ antagonizes FoxO to control growth, apoptosis and longevity in Drosophila , 2008, Aging cell.

[37]  J. Dupont,et al.  IGF Type 1 Receptor: A Cell Cycle Progression Factor That Regulates Aging , 2003, Cell cycle.

[38]  E. Hafen,et al.  The hypoxia-induced paralogs Scylla and Charybdis inhibit growth by down-regulating S6K activity upstream of TSC in Drosophila. , 2004, Genes & development.

[39]  E. Hafen,et al.  Extension of Life-Span by Loss of CHICO, a Drosophila Insulin Receptor Substrate Protein , 2001, Science.

[40]  Stefan Woltran,et al.  Special issue on answer set programming , 2011, AI Commun..

[41]  H. Jasper,et al.  JNK Extends Life Span and Limits Growth by Antagonizing Cellular and Organism-Wide Responses to Insulin Signaling , 2005, Cell.

[42]  Mariano J. Alvarez,et al.  Genome-wide Identification of Post-translational Modulators of Transcription Factor Activity in Human B-Cells , 2009, Nature Biotechnology.

[43]  M. Hoch,et al.  The cytohesin Steppke is essential for insulin signalling in Drosophila , 2006, Nature.

[44]  F. Moreews,et al.  BioQuali Cytoscape plugin: analysing the global consistency of regulatory networks , 2009, BMC Genomics.

[45]  Julio Saez-Rodriguez,et al.  Revisiting the Training of Logic Models of Protein Signaling Networks with a Formal Approach based on Answer Set Programming , 2012, ArXiv.

[46]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[47]  M. Selbach,et al.  Global quantification of mammalian gene expression control , 2011, Nature.

[48]  Chris Cornelis,et al.  Modeling Protein Interaction Networks with Answer Set Programming , 2009, 2009 IEEE International Conference on Bioinformatics and Biomedicine.

[49]  Denis Thieffry,et al.  Mathematical Modelling of Cell-Fate Decision in Response to Death Receptor Engagement , 2010, PLoS Comput. Biol..

[50]  L. Partridge,et al.  A longer and healthier life with TOR down-regulation: genetics and drugs. , 2011, Biochemical Society transactions.

[51]  C. Kenyon,et al.  A C. elegans mutant that lives twice as long as wild type , 1993, Nature.

[52]  A. Sonntag,et al.  A Dynamic Network Model of mTOR Signaling Reveals TSC-Independent mTORC2 Regulation , 2012, Science Signaling.