Statistical inference for time course RNA-Seq data using a negative binomial mixed-effect model
View/ Open
Author
Sun, Xiaoxiao
Dalpiaz, David
Wu, Di
S. Liu, Jun
Zhong, Wenxuan
Ma, Ping
Published Version
https://doi.org/10.1186/s12859-016-1180-9Metadata
Show full item recordCitation
Sun, Xiaoxiao, David Dalpiaz, Di Wu, Jun S. Liu, Wenxuan Zhong, and Ping Ma. 2016. “Statistical inference for time course RNA-Seq data using a negative binomial mixed-effect model.” BMC Bioinformatics 17 (1): 324. doi:10.1186/s12859-016-1180-9. http://dx.doi.org/10.1186/s12859-016-1180-9.Abstract
Background: Accurate identification of differentially expressed (DE) genes in time course RNA-Seq data is crucial for understanding the dynamics of transcriptional regulatory network. However, most of the available methods treat gene expressions at different time points as replicates and test the significance of the mean expression difference between treatments or conditions irrespective of time. They thus fail to identify many DE genes with different profiles across time. In this article, we propose a negative binomial mixed-effect model (NBMM) to identify DE genes in time course RNA-Seq data. In the NBMM, mean gene expression is characterized by a fixed effect, and time dependency is described by random effects. The NBMM is very flexible and can be fitted to both unreplicated and replicated time course RNA-Seq data via a penalized likelihood method. By comparing gene expression profiles over time, we further classify the DE genes into two subtypes to enhance the understanding of expression dynamics. A significance test for detecting DE genes is derived using a Kullback-Leibler distance ratio. Additionally, a significance test for gene sets is developed using a gene set score. Results: Simulation analysis shows that the NBMM outperforms currently available methods for detecting DE genes and gene sets. Moreover, our real data analysis of fruit fly developmental time course RNA-Seq data demonstrates the NBMM identifies biologically relevant genes which are well justified by gene ontology analysis. Conclusions: The proposed method is powerful and efficient to detect biologically relevant DE genes and gene sets in time course RNA-Seq data. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-1180-9) contains supplementary material, which is available to authorized users.Other Sources
http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5002174/pdf/Terms of Use
This article is made available under the terms and conditions applicable to Other Posted Material, as set forth at http://nrs.harvard.edu/urn-3:HUL.InstRepos:dash.current.terms-of-use#LAACitable link to this page
http://nrs.harvard.edu/urn-3:HUL.InstRepos:29407854
Collections
- FAS Scholarly Articles [18292]
Contact administrator regarding this item (to report mistakes or request changes)