Dynamic gene and transcriptional regulatory networks inferring with multi-Laplacian prior from time-course gene microarray data

This paper presents a dynamic gene and transcriptional regulatory network inferring method by using the time-varying autoregressive (TVAR) model. It employs the Li-based regularization terms with spatial sparsity, temporal continuity and proposed multi-Laplacian prior (MLP) for key transcriptional factors (TFs) and their interactions identification. The MLP regularization allows the connections of a gene to be better preserved as a group so that putative TFs can be identified in dynamic gene network. Furthermore, an ADMM-based method is proposed to solve the problem by using the augmented Lagrangian multiplier technique. The simulation using DREAM 4 datasets shows the proposed method performs better than other well-established algorithms for gene network inferring. This enables us to apply the proposed method to a yeast cell cycle microarray datasets containing 215 genes and 17 timepoints more effectively. We are able to identify key genes and gene interactions align well with the natural of yeast cell cycle and related literatures. These suggest that the proposed method can serve as a useful exploratory tool for putative TFs and dynamic gene/TFs networks identification using microarray data.

[1]  Anna Shevchenko,et al.  The budding yeast proteins Spc24p and Spc25p interact with Ndc80p and Nuf2p at the kinetochore and are important for kinetochore clustering and checkpoint control , 2001, The EMBO journal.

[2]  Volker Roth,et al.  The Bayesian group-Lasso for analyzing contingency tables , 2009, ICML '09.

[3]  Robert D. Leclerc Survival of the sparsest: robust gene networks are parsimonious , 2008, Molecular systems biology.

[4]  I R Gibbons,et al.  Saccharomyces cerevisiae kinesin- and dynein-related proteins required for anaphase chromosome segregation , 1995, The Journal of cell biology.

[5]  A. Yakovlev,et al.  A New Type of Stochastic Dependence Revealed in Gene Expression Data , 2006, Statistical applications in genetics and molecular biology.

[6]  Paul P. Wang,et al.  Advances to Bayesian network inference for generating causal networks from observational biological data , 2004, Bioinform..

[7]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[8]  E. Johansson,et al.  The Pol32 Subunit of DNA Polymerase δ Contains Separable Domains for Processive Replication and Proliferating Cell Nuclear Antigen (PCNA) Binding* , 2004, Journal of Biological Chemistry.

[9]  R. Chanet,et al.  Characterization of mutations that are synthetic lethal with pol3-13, a mutated allele of DNA polymerase delta in Saccharomyces cerevisiae , 2003, Current Genetics.

[10]  A. Barabasi,et al.  Network biology: understanding the cell's functional organization , 2004, Nature Reviews Genetics.

[11]  F. Turano,et al.  Expression of a Glutamate Decarboxylase Homologue Is Required for Normal Oxidative Stress Tolerance in Saccharomyces cerevisiae * , 2001, The Journal of Biological Chemistry.

[12]  Li Zhang,et al.  A Maximum A Posteriori Probability and Time-Varying Approach for Inferring Gene Regulatory Networks from Time Course Gene Microarray Data , 2015, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[13]  I R Gibbons,et al.  Cytoplasmic dynein is required for normal nuclear segregation in yeast. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[14]  P. Burgers,et al.  Eukaryotic DNA polymerases in DNA replication and DNA repair , 1998, Chromosoma.

[15]  Ronald W. Davis,et al.  A genome-wide transcriptional analysis of the mitotic cell cycle. , 1998, Molecular cell.

[16]  E. Stadtman,et al.  The regulation of glutamine synthesis in microorganisms. , 1970, Annual review of microbiology.

[17]  U. Hjorth Model Selection and Forward Validation , 1982 .

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

[19]  David A Agard,et al.  Spindle pole body duplication in fission yeast occurs at the G1/S boundary but maturation is blocked until exit from S by an event downstream of cdc10+. , 2004, Molecular biology of the cell.

[20]  D. Floreano,et al.  Revealing strengths and weaknesses of methods for gene network inference , 2010, Proceedings of the National Academy of Sciences.

[21]  Kaoru Irie,et al.  Clb/Cdc28 kinases promote nuclear export of the replication initiator proteins Mcm2–7 , 2000, Current Biology.

[22]  B. Magasanik,et al.  Role of NAD-linked glutamate dehydrogenase in nitrogen metabolism in Saccharomyces cerevisiae , 1990, Journal of bacteriology.

[23]  H. Robichon-Szulmajster,et al.  Methionine biosynthesis in Saccharomyces cerevisiae , 1975, Molecular and General Genetics MGG.

[24]  Carolyn A. Larabell,et al.  The Yeast Polo Kinase Cdc5 Regulates the Shape of the Mitotic Nucleus , 2014, Current Biology.