Inferring stable genetic networks from steady-state data

Gene regulatory networks capture the interactions between genes and other cell substances, resulting from the fundamental biological process of transcription and translation. In some applications, the topology of the regulatory network is not known, and has to be inferred from experimental data. The experimental data consist of expression levels of the genes, which are typically measured as mRNA concentrations in micro-array experiments. In a so-called genetic perturbation experiment, small perturbations are applied to equilibrium states and the resulting changes in expression activity are measured. This paper develops novel algorithms that identify a sparse and stable genetic network that explains data obtained from noisy genetic perturbation experiments. Our identification algorithm is based on convex relaxations of the sparsity and stability constraints and can also incorporate a variety of prior knowledge of the network structure. Such knowledge can be either qualitative, specifying positive, negative or no interactions between genes, or quantitative, specifying a range of interaction strengths. Our approach is applied to both synthetic and experimental data, obtained for the SOS pathway in Escherichia coli, and the results show that the stability specification not only ensures consistency with the steady-state assumptions, but also significantly increases the identification performance. Since the method is based on convex optimization, it can be efficiently applied to large scale networks.

[1]  P. McSharry,et al.  Reconstructing biochemical pathways from time course data , 2007, Proteomics.

[2]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

[3]  Diego di Bernardo,et al.  Robust Identification of Large Genetic Networks , 2003, Pacific Symposium on Biocomputing.

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

[5]  Stephen P. Boyd,et al.  Low-Authority Controller Design by Means of Convex Optimization , 1999 .

[6]  R. Woodgate,et al.  Identification of additional genes belonging to the LexA regulon in Escherichia coli , 2000, Molecular microbiology.

[7]  Antonis Papachristodoulou,et al.  Efficient, sparse biological network determination , 2009, BMC Systems Biology.

[8]  Joel A. Tropp,et al.  Just relax: convex programming methods for identifying sparse signals in noise , 2006, IEEE Transactions on Information Theory.

[9]  Timothy S Gardner,et al.  Reverse-engineering transcription control networks. , 2005, Physics of life reviews.

[10]  Isaac Dialsingh,et al.  Basic Statistics and Pharmaceutical Statistical Applications, 2nd edn , 2007 .

[11]  D. di Bernardo,et al.  How to infer gene networks from expression profiles , 2007, Molecular systems biology.

[12]  Michal Linial,et al.  Using Bayesian Networks to Analyze Expression Data , 2000, J. Comput. Biol..

[13]  A. Papachristodoulou,et al.  Determining Interconnections in Chemical Reaction Networks , 2007, 2007 American Control Conference.

[14]  Max Donath,et al.  American Control Conference , 1993 .

[15]  E. Davidson,et al.  The hardwiring of development: organization and function of genomic regulatory systems. , 1997, Development.

[16]  Kwang-Hyun Cho,et al.  Inferring biomolecular interaction networks based on convex optimization , 2007, Comput. Biol. Chem..

[17]  S. Schreiber,et al.  Target-oriented and diversity-oriented organic synthesis in drug discovery. , 2000, Science.

[18]  J. Hasty,et al.  Reverse engineering gene networks: Integrating genetic perturbations with dynamical modeling , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Ning Li,et al.  Basic Statistics and Pharmaceutical Statistical Applications , 2000 .

[20]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[21]  Hidde de Jong,et al.  Structural Identification of Piecewise-Linear Models of Genetic Regulatory Networks , 2008, J. Comput. Biol..

[22]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[23]  Diego di Bernardo,et al.  LMI-based Algorithm for the Reconstruction of Biological Networks , 2007, 2007 American Control Conference.

[24]  Eduardo D. Sontag,et al.  Inferring dynamic architecture of cellular networks using time series of gene expression, protein and metabolite data , 2004, Bioinform..

[25]  Carsten O. Daub,et al.  The mutual information: Detecting and evaluating dependencies between variables , 2002, ECCB.

[26]  D. Botstein,et al.  Cluster analysis and display of genome-wide expression patterns. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Stephen P. Boyd,et al.  Enhancing Sparsity by Reweighted ℓ1 Minimization , 2007, 0711.1612.

[28]  George J. Pappas,et al.  Genetic network identification using convex programming. , 2009, IET systems biology.

[29]  M. Govindaraju,et al.  The Linear System , 1998 .

[30]  G.J. Pappas,et al.  Identification of stable genetic networks using convex programming , 2008, 2008 American Control Conference.

[31]  Antonino Staiano,et al.  A multi-step approach to time series analysis and gene expression clustering , 2006, Bioinform..

[32]  Diego di Bernardo,et al.  Inference of gene regulatory networks and compound mode of action from time course gene expression profiles , 2006, Bioinform..

[33]  John Lygeros,et al.  Canalizing structure of genetic network dynamics: modelling and identification via mixed-integer programming , 2009, Proceedings of the 48h IEEE Conference on Decision and Control (CDC) held jointly with 2009 28th Chinese Control Conference.

[34]  Roy M. Howard,et al.  Linear System Theory , 1992 .

[35]  Araceli M. Huerta,et al.  From specific gene regulation to genomic networks: a global analysis of transcriptional regulation in Escherichia coli. , 1998, BioEssays : news and reviews in molecular, cellular and developmental biology.

[36]  Stephen P. Boyd,et al.  Low-authority controller design via convex optimization , 1998, Proceedings of the 37th IEEE Conference on Decision and Control (Cat. No.98CH36171).