Network-based comparison of temporal gene expression patterns

MOTIVATION In the pursuits of mechanistic understanding of cell differentiation, it is often necessary to compare multiple differentiation processes triggered by different external stimuli and internal perturbations. Available methods for comparing temporal gene expression patterns are limited to a gene-by-gene approach, which ignores co-expression information and thus is sensitive to measurement noise. METHODS We present a method for co-expression network based comparison of temporal expression patterns (NACEP). NACEP compares the temporal patterns of a gene between two experimental conditions, taking into consideration all of the possible co-expression modules that this gene may participate in. The NACEP program is available at http://biocomp.bioen.uiuc.edu/nacep. RESULTS We applied NACEP to analyze retinoid acid (RA)-induced differentiation of embryonic stem (ES) cells. The analysis suggests that RA may facilitate neural differentiation by inducing the shh and insulin receptor pathways. NACEP was also applied to compare the temporal responses of seven RNA inhibition (RNAi) experiments. As proof of concept, we demonstrate that the difference in the temporal responses to RNAi treatments can be used to derive interaction relationships of transcription factors (TFs), and therefore infer regulatory modules within a transcription network. In particular, the analysis suggested a novel regulatory relationship between two pluripotency regulators, Esrrb and Tbx3, which was supported by in vivo binding of Esrrb to the promoter of Tbx3. AVAILABILITY The NACEP program and the supplementary documents are available at http://biocomp.bioen.uiuc.edu/nacep. CONTACT szhong@illinois.edu SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  M. Tsai,et al.  Chicken Ovalbumin Upstream Promoter-Transcription Factor II (COUP-TFII) regulates growth and patterning of the postnatal mouse cerebellum. , 2009, Developmental biology.

[2]  Thuy T. Truong,et al.  Activin A-Induced Differentiation of Embryonic Stem Cells into Endoderm and Pancreatic Progenitors—The Influence of Differentiation Factors and Culture Conditions , 2009, Stem Cell Reviews and Reports.

[3]  John Quackenbush Microarrays--Guilt by Association , 2003, Science.

[4]  Tao Wang,et al.  Esrrb Activates Oct4 Transcription and Sustains Self-renewal and Pluripotency in Embryonic Stem Cells* , 2008, Journal of Biological Chemistry.

[5]  Zhaohui S. Qin,et al.  Clustering microarray gene expression data using weighted Chinese restaurant process , 2006, Bioinform..

[6]  K. Mikoshiba,et al.  Molecular Properties of Zic Proteins as Transcriptional Regulators and Their Relationship to GLI Proteins* , 2001, The Journal of Biological Chemistry.

[7]  R. Shamir,et al.  Regulatory networks define phenotypic classes of human stem cell lines , 2008, Nature.

[8]  M. Rudnicki,et al.  Retinoic acid-induced neural differentiation of embryonal carcinoma cells , 1983, Molecular and cellular biology.

[9]  John Quackenbush Genomics. Microarrays--guilt by association. , 2003, Science.

[10]  Andrew P. McMahon,et al.  Genomic characterization of Gli-activator targets in sonic hedgehog-mediated neural patterning , 2007, Development.

[11]  S. MacEachern,et al.  A semiparametric Bayesian model for randomised block designs , 1996 .

[12]  Ian Chambers,et al.  Estrogen-Related Receptor Beta Interacts with Oct4 To Positively Regulate Nanog Gene Expression , 2008, Molecular and Cellular Biology.

[13]  T. Ferguson A Bayesian Analysis of Some Nonparametric Problems , 1973 .

[14]  Sheng Zhong,et al.  Dissecting Early Differentially Expressed Genes in a Mixture of Differentiating Embryonic Stem Cells , 2009, PLoS Comput. Biol..

[15]  Megan F. Cole,et al.  Core Transcriptional Regulatory Circuitry in Human Embryonic Stem Cells , 2005, Cell.

[16]  Huda Y Zoghbi,et al.  Deletion of Atoh1 Disrupts Sonic Hedgehog Signaling in the Developing Cerebellum and Prevents Medulloblastoma , 2009, Science.

[17]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[18]  Anbupalam Thalamuthu,et al.  Gene expression Evaluation and comparison of gene clustering methods in microarray analysis , 2006 .

[19]  Radford M. Neal Markov Chain Sampling Methods for Dirichlet Process Mixture Models , 2000 .

[20]  S. MacEachern,et al.  Estimating mixture of dirichlet process models , 1998 .

[21]  Oliver Brüstle,et al.  Retinoic acid induction of ES-cell-derived neurons: the radial glia connection , 2005, Trends in Neurosciences.

[22]  F. Markowetz,et al.  Systems-level dynamic analyses of fate change in murine embryonic stem cells , 2009, Nature.

[23]  C. Antoniak Mixtures of Dirichlet Processes with Applications to Bayesian Nonparametric Problems , 1974 .

[24]  Donatello Telesca,et al.  Differential Expression and Network Inferences through Functional Data Modeling , 2009, Biometrics.

[25]  Jamie M. Verheyden,et al.  FGF-regulated Etv genes are essential for repressing Shh expression in mouse limb buds. , 2009, Developmental cell.

[26]  M. Pirity,et al.  Tbx3: another important piece fitted into the pluripotent stem cell puzzle , 2010, Stem Cell Research & Therapy.

[27]  Ivan Dikic,et al.  PC12 cells overexpressing the insulin receptor undergo insulin-dependent neuronal differentiation , 1994, Current Biology.

[28]  Alexei A. Sharov,et al.  Dissecting Oct3/4-Regulated Gene Networks in Embryonic Stem Cells by Expression Profiling , 2006, PloS one.

[29]  Yuriy L. Orlov,et al.  Tbx3 improves the germ-line competency of induced pluripotent stem cells , 2010, Nature.

[30]  Raivo Kolde,et al.  The FunGenES Database: A Genomics Resource for Mouse Embryonic Stem Cell Differentiation , 2009, PloS one.

[31]  Yuedong Wang Mixed effects smoothing spline analysis of variance , 1998 .

[32]  Ya Zhang,et al.  Clustering of Time-Course Gene Expression Data , 2004 .

[33]  Bruce A. Hay,et al.  Inactivation of Both foxo and reaper Promotes Long-Term Adult Neurogenesis in Drosophila , 2010, Current Biology.

[34]  Radu Dobrin,et al.  Dissecting self-renewal in stem cells with RNA interference , 2006, Nature.

[35]  Hitoshi Niwa,et al.  A parallel circuit of LIF signalling pathways maintains pluripotency of mouse ES cells , 2009, Nature.

[36]  W. Wong,et al.  A gene regulatory network in mouse embryonic stem cells , 2007, Proceedings of the National Academy of Sciences.

[37]  N. D. Clarke,et al.  Integration of External Signaling Pathways with the Core Transcriptional Network in Embryonic Stem Cells , 2008, Cell.

[38]  Adrian E. Raftery,et al.  MCLUST Version 3 for R: Normal Mixture Modeling and Model-Based Clustering † , 2007 .

[39]  Wenxuan Zhong,et al.  A data-driven clustering method for time course gene expression data , 2006, Nucleic acids research.

[40]  John D. Storey,et al.  Significance analysis of time course microarray experiments. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[41]  Kevin Struhl,et al.  Quantitative sequential chromatin immunoprecipitation, a method for analyzing co-occupancy of proteins at genomic regions in vivo. , 2004, Nucleic acids research.

[42]  Fengzhu Sun,et al.  MARD: a new method to detect differential gene expression in treatment-control time courses , 2006, Bioinform..