{"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":1771991149864,"version":"3.50.1"},"reference-count":31,"publisher":"Oxford University Press (OUP)","issue":"4","license":[{"start":{"date-parts":[[2020,11,5]],"date-time":"2020-11-05T00:00:00Z","timestamp":1604534400000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by-nc\/4.0\/"}],"funder":[{"name":"Duchenne Parent Project Netherlands foundation"},{"DOI":"10.13039\/100013465","name":"Association Fran\u00e7aise contre les Myopathies","doi-asserted-by":"publisher","award":["19118"],"award-info":[{"award-number":["19118"]}],"id":[{"id":"10.13039\/100013465","id-type":"DOI","asserted-by":"publisher"}]},{"name":"Duchenne Centrum Nederland and European Commission"},{"name":"Duchenne Parent Project Netherlands Foundation"},{"name":"Duchenne Centrum Nederland"},{"name":"European Commission FP7-HEALTH through the Project Neuromics","award":["305121"],"award-info":[{"award-number":["305121"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2021,7,20]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:p>Time-course RNAseq experiments, where tissues are repeatedly collected from the same subjects, e.g. humans or animals over time or under several different experimental conditions, are becoming more popular due to the reducing sequencing costs. Such designs offer the great potential to identify genes that change over time or progress differently in time across experimental groups. Modelling of the longitudinal gene expression in such time-course RNAseq data is complicated by the serial correlations, missing values due to subject dropout or sequencing errors, long follow up with potentially non-linear progression in time and low number of subjects. Negative Binomial mixed models can address all these issues. However, such models under the maximum likelihood (ML) approach are less popular for RNAseq data due to convergence issues (see, e.g. [1]). We argue in this paper that it is the use of an inaccurate numerical integration method in combination with the typically small sample sizes which causes such mixed models to fail for a great portion of tested genes. We show that when we use the accurate adaptive Gaussian quadrature approach to approximate the integrals over the random-effects terms, we can successfully estimate the model parameters with the maximum likelihood method. Moreover, we show that the boostrap method can be used to preserve the type I error rate in small sample settings. We evaluate empirically the small sample properties of the test statistics and compare with state-of-the-art approaches. The method is applied on a longitudinal mice experiment to study the dynamics in Duchenne Muscular Dystrophy. Contact:s.tsonaka@lumc.nl<\/jats:p><jats:p>Roula Tsonaka is an assistant professor at the Medical Statistics, Department of Biomedical Data Sciences, Leiden University Medical Center. Her research focuses on statistical methods for longitudinal omics data. Pietro Spitali is an assistant professor at the Department of Human Genetics, Leiden University Medical Center. His research focuses on the identification of biomarkers for neuromuscular disorders.<\/jats:p>","DOI":"10.1093\/bib\/bbaa264","type":"journal-article","created":{"date-parts":[[2020,9,18]],"date-time":"2020-09-18T11:23:27Z","timestamp":1600428207000},"source":"Crossref","is-referenced-by-count":11,"title":["Negative Binomial mixed models estimated with the maximum likelihood method can be used for longitudinal RNAseq data"],"prefix":"10.1093","volume":"22","author":[{"given":"Roula","family":"Tsonaka","sequence":"first","affiliation":[{"name":"Medical Statistics section, Department of Biomedical Data Sciences, Leiden University Medical Center"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Pietro","family":"Spitali","sequence":"additional","affiliation":[{"name":"Department of Human Genetics, Leiden University Medical Center"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2020,11,5]]},"reference":[{"key":"2021072112310943900_ref1","doi-asserted-by":"crossref","first-page":"87","DOI":"10.1515\/sagmb-2015-0011","article-title":"What if we ignore the random effects when analyzing RNA-seq data in a multifactor experiment","volume":"15","author":"Cui","year":"2016","journal-title":"Stat Appl Genet Mol Biol"},{"key":"2021072112310943900_ref2","volume-title":"Applied Longitudinal Analysis","author":"Fitzmaurice","year":"2004"},{"key":"2021072112310943900_ref3","doi-asserted-by":"crossref","first-page":"139","DOI":"10.1093\/bioinformatics\/btp616","article-title":"edgeR: a bioconductor package for differential expression analysis of digital gene expression data","volume":"26","author":"Robinson","year":"2010","journal-title":"Bioinformatics"},{"key":"2021072112310943900_ref4","doi-asserted-by":"crossref","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","article-title":"Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2","volume":"15","author":"Love","year":"2014","journal-title":"Genome Biol"},{"key":"2021072112310943900_ref5","doi-asserted-by":"crossref","first-page":"R29","DOI":"10.1186\/gb-2014-15-2-r29","article-title":"Voom: precision weights unlock linear model analysis tools for RNA-seq read counts","volume":"15","author":"Law","year":"2014","journal-title":"Genome Biol"},{"key":"2021072112310943900_ref6","doi-asserted-by":"crossref","first-page":"139","DOI":"10.1146\/annurev-biodatasci-072018-021255","article-title":"RNA sequencing data: Hitchhiker\u2019s guide to expression analysis","volume":"2","author":"van den Berge","year":"2019","journal-title":"Annu Rev Biomed Data Sci"},{"key":"2021072112310943900_ref7","doi-asserted-by":"crossref","first-page":"2598","DOI":"10.1093\/bioinformatics\/btu333","article-title":"Next maSigPro: updating maSigPro bioconductor package for RNA-seq time series","volume":"30","author":"Nueda","year":"2014","journal-title":"Bioinformatics"},{"key":"2021072112310943900_ref8","doi-asserted-by":"crossref","first-page":"324","DOI":"10.1186\/s12859-016-1180-9","article-title":"Statistical inference for time course RNA-Seq data using a negative binomial mixed-effect model","volume":"17","author":"Sun","year":"2016","journal-title":"BMC Bioinformatics"},{"key":"2021072112310943900_ref9","doi-asserted-by":"crossref","first-page":"116","DOI":"10.1186\/1471-2105-15-116","article-title":"ShrinkBayes: a versatile R-package for analysis of count-based sequencing data in complex study designs","volume":"15","author":"van de Wiel","year":"2014","journal-title":"BMC Bioinformatics"},{"key":"2021072112310943900_ref10","doi-asserted-by":"crossref","first-page":"2067","DOI":"10.1093\/bioinformatics\/bti270","article-title":"The use of within-array replicate spots for assessing differential expression in microarray experiments","volume":"21","author":"Smyth","year":"2005","journal-title":"Bioinformatics"},{"key":"2021072112310943900_ref11","volume-title":"edgeR: Differential Analysis of Sequence Read Count Data User\u2019s Guide","author":"Chen","year":"2020"},{"key":"2021072112310943900_ref12","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1145\/3195833","article-title":"Triclustering algorithms for three-dimensional data analysis: a comprehensive survey","volume":"51","author":"Henriques","year":"2018","journal-title":"ACM Comput Surv"},{"key":"2021072112310943900_ref13","doi-asserted-by":"crossref","first-page":"265","DOI":"10.1111\/1467-9868.00176","article-title":"Maximizing generalized linear mixed model likelihoods with an automated Monte Carlo EM algorithm","volume":"61","author":"Booth","year":"1999","journal-title":"J R Stat Soc Series B Stat Methodology"},{"key":"2021072112310943900_ref14","doi-asserted-by":"crossref","first-page":"12","DOI":"10.1080\/10618600.1995.10474663","article-title":"Approximations to the log-likelihood function in the nonlinear mixed-effects model","volume":"4","author":"Pinheiro","year":"1995","journal-title":"J Comput Graph Stat"},{"key":"2021072112310943900_ref15","doi-asserted-by":"crossref","first-page":"82","DOI":"10.1080\/01621459.1986.10478240","article-title":"Accurate approximations for posterior moments and marginal densities","volume":"81","author":"Tierny","year":"1986","journal-title":"J Am Stat Assoc"},{"key":"2021072112310943900_ref16","doi-asserted-by":"crossref","first-page":"1","DOI":"10.18637\/jss.v067.i01","article-title":"Fitting linear mixed-effects models using lme4","volume":"67","author":"Bates","year":"2015","journal-title":"J Stat Softw"},{"key":"2021072112310943900_ref17","volume-title":"Generalized, Linear, and Mixed Models","author":"McCulloch","year":"2008"},{"key":"2021072112310943900_ref18","volume-title":"GLMMadaptive: Generalized Linear Mixed Models Using Adaptive Gaussian Quadrature. R Package Version 0.6-5","author":"Rizopoulos","year":"2019"},{"key":"2021072112310943900_ref19","article-title":"Generalized additive models","volume-title":"Statistical Models in S","author":"Hastie","year":"1992"},{"key":"2021072112310943900_ref20","doi-asserted-by":"crossref","first-page":"R106","DOI":"10.1186\/gb-2010-11-10-r106","article-title":"Differential expression analysis for sequence count data","volume":"11","author":"Anders","year":"2010","journal-title":"Genome Biol"},{"key":"2021072112310943900_ref21","doi-asserted-by":"crossref","first-page":"R25","DOI":"10.1186\/gb-2010-11-3-r25","article-title":"A scaling normalization method for differential expression analysis of RNA-seq data","volume":"11","author":"Robinson","year":"2010","journal-title":"Genome Biol"},{"key":"2021072112310943900_ref22","doi-asserted-by":"crossref","DOI":"10.1201\/9781315370279","volume-title":"Generalized Additive Models: An Introduction with R","author":"Wood","year":"2017"},{"key":"2021072112310943900_ref23","doi-asserted-by":"crossref","first-page":"9","DOI":"10.1080\/01621459.1993.10594284","article-title":"Approximate inference in generalized linear mixed models","volume":"88","author":"Breslow","year":"1993","journal-title":"J Amer Statist Soc"},{"key":"2021072112310943900_ref24","doi-asserted-by":"crossref","first-page":"719","DOI":"10.1093\/biomet\/78.4.719","article-title":"Estimation in generalised linear models with random effects","volume":"78","author":"Schall","year":"1991","journal-title":"Biometrika"},{"key":"2021072112310943900_ref25","doi-asserted-by":"crossref","first-page":"961","DOI":"10.2307\/2531147","article-title":"Random-effects model for serial observations with binary response","volume":"40","author":"Stiratelli","year":"1984","journal-title":"Biometrics"},{"key":"2021072112310943900_ref26","doi-asserted-by":"crossref","first-page":"45","DOI":"10.1093\/biomet\/78.1.45","article-title":"Nonlinear multilevel models, with an application to discrete response data","volume":"78","author":"Goldstein","year":"1991","journal-title":"Biometrika"},{"key":"2021072112310943900_ref27","doi-asserted-by":"crossref","first-page":"1","DOI":"10.18637\/jss.v070.i05","article-title":"TMB: automatic differentiation and Laplace approximation","volume":"70","author":"Kristensen","year":"2016","journal-title":"J Stat Softw"},{"key":"2021072112310943900_ref28","doi-asserted-by":"crossref","DOI":"10.1007\/978-1-4899-4541-9","volume-title":"An Introduction to the Bootstrap","author":"Efron","year":"1993"},{"key":"2021072112310943900_ref29","article-title":"Springer Series in Statistics","volume-title":"Linear Mixed Models for Longitudinal Data","author":"Verbeke","year":"2000"},{"key":"2021072112310943900_ref30","volume-title":"Theoretical Statistics","author":"Cox","year":"1990"},{"key":"2021072112310943900_ref31","doi-asserted-by":"crossref","first-page":"289","DOI":"10.1111\/j.2517-6161.1995.tb02031.x","article-title":"Controlling the false discovery rate: a practical and powerful approach to multiple testing","volume":"57","author":"Benjamini","year":"1995","journal-title":"J. Roy. Statist. Soc. Ser. B"}],"container-title":["Briefings in Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/academic.oup.com\/bib\/article-pdf\/22\/4\/bbaa264\/39135990\/bbaa264.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"http:\/\/academic.oup.com\/bib\/article-pdf\/22\/4\/bbaa264\/39135990\/bbaa264.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,8,14]],"date-time":"2024-08-14T03:36:00Z","timestamp":1723606560000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bib\/article\/doi\/10.1093\/bib\/bbaa264\/5956897"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2020,11,5]]},"references-count":31,"journal-issue":{"issue":"4","published-print":{"date-parts":[[2021,7,20]]}},"URL":"https:\/\/doi.org\/10.1093\/bib\/bbaa264","relation":{},"ISSN":["1467-5463","1477-4054"],"issn-type":[{"value":"1467-5463","type":"print"},{"value":"1477-4054","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2021,7]]},"published":{"date-parts":[[2020,11,5]]},"article-number":"bbaa264"}}