A hierarchical Bayesian mixture model for inferring the expression state of genes in transcriptomes
Open Access
- 24 July 2020
- journal article
- research article
- Published by Proceedings of the National Academy of Sciences in Proceedings of the National Academy of Sciences of the United States of America
- Vol. 117 (32), 19339-19346
- https://doi.org/10.1073/pnas.1919748117
Abstract
Transcriptomes are key to understanding the relationship between genotype and phenotype. The ability to infer the expression state (active or inactive) of genes in the transcriptome offers unique benefits for addressing this issue. For example, qualitative changes in gene expression may underly the origin of novel phenotypes, and expression states are readily comparable between tissues and species. However, inferring the expression state of genes is a surprisingly difficult problem, owing to the complex biological and technical processes that give rise to observed transcriptomic datasets. Here, we develop a hierarchical Bayesian mixture model that describes this complex process and allows us to infer expression state of genes from replicate transcriptomic libraries. We explore the statistical behavior of this method with analyses of simulated datasets—where we demonstrate its ability to correctly infer true (known) expression states—and empirical-benchmark datasets, where we demonstrate that the expression states inferred from RNA-sequencing (RNA-seq) datasets using our method are consistent with those based on independent evidence. The power of our method to correctly infer expression states is generally high and remarkably, approaches the maximum possible power for this inference problem. We present an empirical analysis of primate-brain transcriptomes, which identifies genes that have a unique expression state in humans. Our method is implemented in the freely available R package zigzag.Keywords
Funding Information
- Foundation for the National Institutes of Health (5F32GM125107-02)
- Foundation for the National Institutes of Health (R35GM122592)
- National Science Foundation (DEB-0842181)
- National Science Foundation (DEB-0919529)
- National Science Foundation (DBI-1356737)
This publication has 47 references indexed in Scilit:
- Multiplatform single-sample estimates of transcriptional activationProceedings of the National Academy of Sciences of the United States of America, 2013
- Quantitative Analysis of Fission Yeast Transcriptomes and Proteomes in Proliferating and Quiescent CellsCell, 2012
- Landscape of transcription in human cellsNature, 2012
- Insights into the regulation of protein abundance from proteomic and transcriptomic analysesNature Reviews Genetics, 2012
- Global analysis of genome, transcriptome and proteome reveals the response to aneuploidy in human cellsMolecular Systems Biology, 2012
- Global quantification of mammalian gene expression controlNature, 2011
- Mapping and analysis of chromatin state dynamics in nine human cell typesNature, 2011
- Protein abundances are more conserved than mRNA abundances across diverse taxaProteomics, 2010
- Mapping and quantifying mammalian transcriptomes by RNA-SeqNature Methods, 2008
- Reversible jump Markov chain Monte Carlo computation and Bayesian model determinationBiometrika, 1995