Learning gene network structure from time laps cell imaging in RNAi Knock downs

MOTIVATION As RNA interference is becoming a standard method for targeted gene perturbation, computational approaches to reverse engineer parts of biological networks based on measurable effects of RNAi become increasingly relevant. The vast majority of these methods use gene expression data, but little attention has been paid so far to other data types. RESULTS Here we present a method, which can infer gene networks from high-dimensional phenotypic perturbation effects on single cells recorded by time-lapse microscopy. We use data from the Mitocheck project to extract multiple shape, intensity and texture features at each frame. Features from different cells and movies are then aligned along the cell cycle time. Subsequently we use Dynamic Nested Effects Models (dynoNEMs) to estimate parts of the network structure between perturbed genes via a Markov Chain Monte Carlo approach. Our simulation results indicate a high reconstruction quality of this method. A reconstruction based on 22 gene knock downs yielded a network, where all edges could be explained via the biological literature. AVAILABILITY The implementation of dynoNEMs is part of the Bioconductor R-package nem.

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

[2]  Peter Bühlmann,et al.  Predicting causal effects in large-scale systems from observational data , 2010, Nature Methods.

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

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

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

[6]  M V Boland,et al.  Automated recognition of patterns characteristic of subcellular structures in fluorescence microscopy images. , 1998, Cytometry.

[7]  Holger Fröhlich,et al.  MC EMiNEM Maps the Interaction Landscape of the Mediator , 2012, PLoS Comput. Biol..

[8]  Richard Gran,et al.  On the Convergence of Random Search Algorithms In Continuous Time with Applications to Adaptive Control , 1970, IEEE Trans. Syst. Man Cybern..

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

[10]  Anne E Carpenter,et al.  Improved structure, function and compatibility for CellProfiler: modular high-throughput image analysis software , 2011, Bioinform..

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

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

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

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

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

[16]  C. Bakal,et al.  Quantitative Morphological Signatures Define Local Signaling Networks Regulating Cell Morphology , 2007, Science.

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

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

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

[20]  M. Maathuis,et al.  Estimating high-dimensional intervention effects from observational data , 2008, 0810.4214.

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

[22]  Xiaobo Zhou,et al.  Automated segmentation, classification, and tracking of cancer cell nuclei in time-lapse microscopy , 2006, IEEE Transactions on Biomedical Engineering.

[23]  Joshua M. Stuart,et al.  Inferring Disease-Related Pathways Using A Probabilistic Epistasis Model , 2008, Pacific Symposium on Biocomputing.

[24]  Robert J. Flassig,et al.  TRANSWESD: inferring cellular networks with transitive reduction , 2010, Bioinform..

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

[26]  S. Reed,et al.  Nuclear Accumulation of p21Cip1 at the Onset of Mitosis: a Role at the G2/M-Phase Transition , 1998, Molecular and Cellular Biology.

[27]  Holger Fröhlich,et al.  Fast and efficient dynamic nested effects models , 2011, Bioinform..

[28]  R. Durbin,et al.  Phenotypic profiling of the human genome by time-lapse microscopy reveals cell division genes , 2010, Nature.

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

[30]  Lars Kaderali,et al.  Reconstructing signaling pathways from RNAi data using probabilistic Boolean threshold networks , 2009, Bioinform..

[31]  N. Otsu.  A threshold selection method from gray level histograms , 1979 .

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

[33]  J J Vaquero,et al.  Applying watershed algorithms to the segmentation of clustered nuclei. , 1998, Cytometry.

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

[35]  J. R. Smith,et al.  Inhibition of E2F activity by the cyclin-dependent protein kinase inhibitor p21 in cells expressing or lacking a functional retinoblastoma protein , 1996, Molecular and cellular biology.