{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,25]],"date-time":"2026-02-25T03:45:49Z","timestamp":1771991149853,"version":"3.50.1"},"reference-count":29,"publisher":"Springer Science and Business Media LLC","issue":"1","license":[{"start":{"date-parts":[[2016,8,26]],"date-time":"2016-08-26T00:00:00Z","timestamp":1472169600000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2016,8,26]],"date-time":"2016-08-26T00:00:00Z","timestamp":1472169600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"funder":[{"DOI":"10.13039\/100000001","name":"National Science Foundation","doi-asserted-by":"publisher","award":["DMS-1440037"],"award-info":[{"award-number":["DMS-1440037"]}],"id":[{"id":"10.13039\/100000001","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000001","name":"National Science Foundation","doi-asserted-by":"publisher","award":["DMS-1438957"],"award-info":[{"award-number":["DMS-1438957"]}],"id":[{"id":"10.13039\/100000001","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000001","name":"National Science Foundation","doi-asserted-by":"publisher","award":["DMS-1440038"],"award-info":[{"award-number":["DMS-1440038"]}],"id":[{"id":"10.13039\/100000001","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["1R01GM113242-01"],"award-info":[{"award-number":["1R01GM113242-01"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["1R01GM122080-01"],"award-info":[{"award-number":["1R01GM122080-01"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["BMC Bioinformatics"],"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Background<\/jats:title><jats:p>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.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>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.<\/jats:p><\/jats:sec><jats:sec><jats:title>Conclusions<\/jats:title><jats:p>The proposed method is powerful and efficient to detect biologically relevant DE genes and gene sets in time course RNA-Seq data.<\/jats:p><\/jats:sec>","DOI":"10.1186\/s12859-016-1180-9","type":"journal-article","created":{"date-parts":[[2016,8,26]],"date-time":"2016-08-26T13:07:02Z","timestamp":1472216822000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":24,"title":["Statistical inference for time course RNA-Seq data using a negative binomial mixed-effect model"],"prefix":"10.1186","volume":"17","author":[{"given":"Xiaoxiao","family":"Sun","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"David","family":"Dalpiaz","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Di","family":"Wu","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Jun","family":"S. Liu","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Wenxuan","family":"Zhong","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Ping","family":"Ma","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"297","published-online":{"date-parts":[[2016,8,26]]},"reference":[{"issue":"1","key":"1180_CR1","doi-asserted-by":"publisher","first-page":"62","DOI":"10.1214\/10-STS343","volume":"26","author":"J Salzman","year":"2011","unstructured":"Salzman J, Jiang H, Wong WH. Statistical modeling of RNA-seq data. Stat Sci. 2011; 26(1):62\u201383.","journal-title":"Stat Sci."},{"issue":"7","key":"1180_CR2","doi-asserted-by":"publisher","first-page":"621","DOI":"10.1038\/nmeth.1226","volume":"5","author":"A Mortazavi","year":"2008","unstructured":"Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-seq. Nature Methods. 2008; 5(7):621\u20138.","journal-title":"Nature Methods"},{"issue":"5881","key":"1180_CR3","doi-asserted-by":"publisher","first-page":"1344","DOI":"10.1126\/science.1158441","volume":"320","author":"U Nagalakshmi","year":"2008","unstructured":"Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M. The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008; 320(5881):1344\u20131349.","journal-title":"Science"},{"issue":"7199","key":"1180_CR4","doi-asserted-by":"publisher","first-page":"1239","DOI":"10.1038\/nature07002","volume":"453","author":"BT Wilhelm","year":"2008","unstructured":"Wilhelm BT, Marguerat S, Watt S, Schubert F, Wood V, Goodhead I, Penkett CJ, Rogers J, B\u00e4hler J. Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution. Nature. 2008; 453(7199):1239\u20131243.","journal-title":"Nature"},{"issue":"1","key":"1180_CR5","doi-asserted-by":"publisher","first-page":"139","DOI":"10.1093\/bioinformatics\/btp616","volume":"26","author":"MD Robinson","year":"2010","unstructured":"Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinforma. 2010; 26(1):139\u201340.","journal-title":"Bioinforma"},{"key":"1180_CR6","doi-asserted-by":"publisher","first-page":"106","DOI":"10.1186\/gb-2010-11-10-r106","volume":"11","author":"S Anders","year":"2010","unstructured":"Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biology. 2010; 11:106.","journal-title":"Genome Biology"},{"key":"1180_CR7","first-page":"1","volume":"2013","author":"S Oh","year":"2013","unstructured":"Oh S, Song S, Grabowski G, Zhao H, Noonan JP. Time series expression analyses using RNA-seq: a statistical approach. BioMed Res Int. 2013; 2013:1\u201316.","journal-title":"BioMed Res Int."},{"issue":"2","key":"1180_CR8","doi-asserted-by":"publisher","first-page":"144","DOI":"10.1007\/s12561-009-9014-1","volume":"1","author":"P Ma","year":"2009","unstructured":"Ma P, Zhong W, Liu JS. Identifying differentially expressed genes in time course microarray data. Stat Biosci. 2009; 1(2):144\u201359.","journal-title":"Stat Biosci"},{"issue":"2","key":"1180_CR9","doi-asserted-by":"publisher","first-page":"485","DOI":"10.1198\/106186005X47651","volume":"14","author":"C Gu","year":"2005","unstructured":"Gu C, Ma P. Generalized nonparametric mixed-effect models: computation and smoothing parameter selection. J Comput Graph Stat. 2005; 14(2):485\u2013504.","journal-title":"J Comput Graph Stat"},{"issue":"4","key":"1180_CR10","doi-asserted-by":"publisher","first-page":"377","DOI":"10.1007\/BF01404567","volume":"31","author":"P Craven","year":"1978","unstructured":"Craven P, Wahba G. Smoothing noisy data with spline functions. Numerische Mathematik. 1978; 31(4):377\u2013403.","journal-title":"Numerische Mathematik"},{"key":"1180_CR11","doi-asserted-by":"publisher","DOI":"10.1007\/978-1-4614-5369-7","volume-title":"Smoothing Spline ANOVA Models","author":"C Gu","year":"2013","unstructured":"Gu C. Smoothing Spline ANOVA Models, Second. New York: Springer; 2013."},{"issue":"3","key":"1180_CR12","doi-asserted-by":"publisher","first-page":"25","DOI":"10.1186\/gb-2010-11-3-r25","volume":"11","author":"M Robinson","year":"2010","unstructured":"Robinson M, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010; 11(3):25.","journal-title":"Genome Biol"},{"key":"1180_CR13","doi-asserted-by":"publisher","DOI":"10.1137\/1.9781611970128","volume-title":"Spline Models for Observational Data","author":"G Wahba","year":"1990","unstructured":"Wahba G. Spline Models for Observational Data. Philadelphia: Society for Industrial and Applied Mathematics; 1990."},{"key":"1180_CR14","doi-asserted-by":"publisher","first-page":"15","DOI":"10.1214\/ss\/1177011926","volume":"6","author":"GK Robinson","year":"1991","unstructured":"Robinson GK. That BLUP is a good thing: The estimation of random effects. Stat Sci. 1991; 6:15\u201332.","journal-title":"Stat Sci."},{"issue":"4","key":"1180_CR15","doi-asserted-by":"publisher","first-page":"347","DOI":"10.2307\/3316020","volume":"32","author":"C Gu","year":"2004","unstructured":"Gu C. Model diagnostics for smoothing spline ANOVA models. Canadian J Stat. 2004; 32(4):347\u201358.","journal-title":"Canadian J Stat"},{"issue":"43","key":"1180_CR16","doi-asserted-by":"publisher","first-page":"15545","DOI":"10.1073\/pnas.0506580102","volume":"102","author":"A Subramanian","year":"2005","unstructured":"Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005; 102(43):15545\u201315550.","journal-title":"Proc Natl Acad Sci"},{"issue":"1","key":"1180_CR17","doi-asserted-by":"publisher","first-page":"107","DOI":"10.1214\/07-AOAS101","volume":"1","author":"B Efron","year":"2007","unstructured":"Efron B, Tibshirani R. On testing the significance of sets of genes. Annals Appl Stat. 2007; 1(1):107\u201329.","journal-title":"Annals Appl Stat"},{"key":"1180_CR18","doi-asserted-by":"crossref","first-page":"289","DOI":"10.1111\/j.2517-6161.1995.tb02031.x","volume":"57","author":"Y Benjamini","year":"1995","unstructured":"Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Royal Stat Soc Ser B (Methodological), J. 1995; 57:289\u2013300.","journal-title":"Royal Stat Soc Ser B (Methodological), J"},{"issue":"18","key":"1180_CR19","doi-asserted-by":"publisher","first-page":"2598","DOI":"10.1093\/bioinformatics\/btu333","volume":"30","author":"MJ Nueda","year":"2014","unstructured":"Nueda MJ, Tarazona S, Conesa A. Next maSigPro: updating maSigPro bioconductor package for RNA-seq time series. Bioinformatics. 2014; 30(18):2598\u2013602.","journal-title":"Bioinformatics"},{"issue":"12","key":"1180_CR20","doi-asserted-by":"publisher","first-page":"113","DOI":"10.1093\/bioinformatics\/btu274","volume":"30","author":"T \u00c4ij\u00f6","year":"2014","unstructured":"\u00c4ij\u00f6 T, Butty V, Chen Z, Salo V, Tripathi S, Burge CB, Lahesmaa R, L\u00e4hdesm\u00e4ki H. Methods for time series analysis of RNA-seq data with application to human Th17 cell differentiation. Bioinformatics. 2014; 30(12):113\u201320.","journal-title":"Bioinformatics"},{"issue":"17","key":"1180_CR21","doi-asserted-by":"publisher","first-page":"2778","DOI":"10.1093\/bioinformatics\/btv272","volume":"31","author":"AC Frazee","year":"2015","unstructured":"Frazee AC, Jaffe AE, Langmead B, Leek JT. Polyester: simulating RNA-seq datasets with differential transcript expression. Bioinformatics. 2015; 31(17):2778\u2013784.","journal-title":"Bioinformatics"},{"issue":"3","key":"1180_CR22","doi-asserted-by":"publisher","first-page":"562","DOI":"10.1038\/nprot.2012.016","volume":"7","author":"C Trapnell","year":"2012","unstructured":"Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nature Protocols. 2012; 7(3):562\u201378.","journal-title":"Nature Protocols"},{"issue":"16","key":"1180_CR23","doi-asserted-by":"publisher","first-page":"2078","DOI":"10.1093\/bioinformatics\/btp352","volume":"25","author":"H Li","year":"2009","unstructured":"Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, et al. The sequence alignment\/map format and samtools. Bioinformatics. 2009; 25(16):2078\u2013079.","journal-title":"Bioinformatics"},{"issue":"10","key":"1180_CR24","doi-asserted-by":"publisher","first-page":"2008","DOI":"10.1101\/gr.133744.111","volume":"22","author":"S Anders","year":"2012","unstructured":"Anders S, Reyes A, Huber W. Detecting differential usage of exons from RNA-seq data. Genome Res. 2012; 22(10):2008\u2013017.","journal-title":"Genome Res"},{"issue":"7339","key":"1180_CR25","doi-asserted-by":"publisher","first-page":"473","DOI":"10.1038\/nature09715","volume":"471","author":"BR Graveley","year":"2011","unstructured":"Graveley BR, Brooks AN, Carlson JW, Duff MO, Landolin JM, Yang L, Artieri CG, van Baren MJ, Boley N, Booth BW, et al.The developmental transcriptome of Drosophila melanogaster. Nature. 2011; 471(7339):473\u20139.","journal-title":"Nature"},{"issue":"3","key":"1180_CR26","doi-asserted-by":"publisher","first-page":"25","DOI":"10.1186\/gb-2009-10-3-r25","volume":"10","author":"B Langmead","year":"2009","unstructured":"Langmead B, Trapnell C, Pop M, Salzberg S. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009; 10(3):25.","journal-title":"Genome Biol"},{"key":"1180_CR27","unstructured":"Campos-Ortega JA, Hartenstein V. The Embryonic Development of Drosophila Melanogaster. Berlin: Springer."},{"issue":"1","key":"1180_CR28","doi-asserted-by":"publisher","first-page":"25","DOI":"10.1038\/75556","volume":"25","author":"M Ashburner","year":"2000","unstructured":"Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. Nature Genetics. 2000; 25(1):25\u20139.","journal-title":"Nature Genetics"},{"issue":"1","key":"1180_CR29","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1093\/nar\/gkn923","volume":"37","author":"DW Huang","year":"2009","unstructured":"Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009; 37(1):1\u201313.","journal-title":"Nucleic Acids Res"}],"container-title":["BMC Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-016-1180-9.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1186\/s12859-016-1180-9\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"http:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-016-1180-9","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-016-1180-9.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,6,19]],"date-time":"2024-06-19T02:32:19Z","timestamp":1718764339000},"score":1,"resource":{"primary":{"URL":"https:\/\/bmcbioinformatics.biomedcentral.com\/articles\/10.1186\/s12859-016-1180-9"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2016,8,26]]},"references-count":29,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2016,12]]}},"alternative-id":["1180"],"URL":"https:\/\/doi.org\/10.1186\/s12859-016-1180-9","relation":{},"ISSN":["1471-2105"],"issn-type":[{"value":"1471-2105","type":"electronic"}],"subject":[],"published":{"date-parts":[[2016,8,26]]},"assertion":[{"value":"3 March 2016","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"11 August 2016","order":2,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"26 August 2016","order":3,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}}],"article-number":"324"}}