Efficient inference for sparse latent variable models of transcriptional regulation

Abstract Motivation Regulation of gene expression in prokaryotes involves complex co-regulatory mechanisms involving large numbers of transcriptional regulatory proteins and their target genes. Uncovering these genome-scale interactions constitutes a major bottleneck in systems biology. Sparse latent factor models, assuming activity of transcription factors (TFs) as unobserved, provide a biologically interpretable modelling framework, integrating gene expression and genome-wide binding data, but at the same time pose a hard computational inference problem. Existing probabilistic inference methods for such models rely on subjective filtering and suffer from scalability issues, thus are not well-suited for realistic genome-scale applications. Results We present a fast Bayesian sparse factor model, which takes input gene expression and binding sites data, either from ChIP-seq experiments or motif predictions, and outputs active TF-gene links as well as latent TF activities. Our method employs an efficient variational Bayes scheme for model inference enabling its application to large datasets which was not feasible with existing MCMC-based inference methods for such models. We validate our method on synthetic data against a similar model in the literature, employing MCMC for inference, and obtain comparable results with a small fraction of the computational time. We also apply our method to large-scale data from Mycobacterium tuberculosis involving ChIP-seq data on 113 TFs and matched gene expression data for 3863 putative target genes. We evaluate our predictions using an independent transcriptomics experiment involving over-expression of TFs. Availability and implementation An easy-to-use Jupyter notebook demo of our method with data is available at https://github.com/zhenwendai/SITAR. Supplementary information Supplementary data are available at Bioinformatics online.

[1]  Fabio Rinaldi,et al.  RegulonDB version 9.0: high-level integration of gene regulation, coexpression, motif clustering and beyond , 2015, Nucleic Acids Res..

[2]  Neil D. Lawrence,et al.  Probabilistic inference of transcription factor concentrations and gene-specific regulatory activities , 2006, Bioinform..

[3]  Kenta Nakai,et al.  DBTBS: a database of transcriptional regulation in Bacillus subtilis containing upstream intergenic conservation information , 2007, Nucleic Acids Res..

[4]  M. West,et al.  High-Dimensional Sparse Factor Modeling: Applications in Gene Expression Genomics , 2008, Journal of the American Statistical Association.

[5]  Mario L. Arrieta-Ortiz,et al.  An experimentally supported model of the Bacillus subtilis global transcriptional regulatory network , 2015, Molecular systems biology.

[6]  Guido Sanguinetti,et al.  Large-scale learning of combinatorial transcriptional dynamics from gene expression , 2011, Bioinform..

[7]  Chiara Sabatti,et al.  Bayesian sparse hidden components analysis for transcription regulation networks , 2005, Bioinform..

[8]  Neil D. Lawrence,et al.  Probabilistic Non-linear Principal Component Analysis with Gaussian Process Latent Variable Models , 2005, J. Mach. Learn. Res..

[9]  Oliver Stegle,et al.  Scalable latent-factor models applied to single-cell RNA-seq data separate biological drivers from confounding effects , 2016, bioRxiv.

[10]  Kyle J. Minch,et al.  A comprehensive map of genome-wide gene regulation in Mycobacterium tuberculosis , 2015, Scientific Data.

[11]  Martin Tompa,et al.  Rv3133c/dosR is a transcription factor that mediates the hypoxic response of Mycobacterium tuberculosis , 2003, Molecular microbiology.

[12]  Nathan D. Price,et al.  The DNA-binding network of Mycobacterium tuberculosis , 2015, Nature Communications.

[13]  Lorenz Wernisch,et al.  Factor analysis for gene regulatory networks and transcription factor activity profiles , 2007, BMC Bioinformatics.

[14]  Chiara Sabatti,et al.  Network component analysis: Reconstruction of regulatory signals in biological systems , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[15]  David J Studholme,et al.  Bmc Microbiology Bioinformatic Identification of Novel Regulatory Dna Sequence Motifs in Streptomyces Coelicolor , 2004 .

[16]  Diogo M. Camacho,et al.  Wisdom of crowds for robust gene network inference , 2012, Nature Methods.

[17]  Eric D Siggia,et al.  Identification of the binding sites of regulatory proteins in bacterial genomes , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Kyle J. Minch,et al.  Mapping and manipulating the Mycobacterium tuberculosis transcriptome using a transcription factor overexpression-derived regulatory network , 2014, Genome Biology.

[19]  Yves Van de Peer,et al.  The Mycobacterium tuberculosis regulatory network and hypoxia , 2013, Nature.

[20]  Neil D. Lawrence,et al.  Bayesian Gaussian Process Latent Variable Model , 2010, AISTATS.

[21]  David A. Hodgson,et al.  Extracting regulator activity profiles by integration of de novo motifs and expression data: characterizing key regulators of nutrient depletion responses in Streptomyces coelicolor , 2012, Nucleic acids research.

[22]  James E. Galagan,et al.  Genomic insights into tuberculosis , 2014, Nature Reviews Genetics.

[23]  Michalis K. Titsias,et al.  Variational Learning of Inducing Variables in Sparse Gaussian Processes , 2009, AISTATS.

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

[25]  Fabian J Theis,et al.  Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells , 2015, Nature Biotechnology.

[26]  Kay Nieselt,et al.  The dynamic architecture of the metabolic switch in Streptomyces coelicolor , 2010, BMC Genomics.