{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,8]],"date-time":"2026-02-08T08:28:02Z","timestamp":1770539282169,"version":"3.49.0"},"reference-count":36,"publisher":"Oxford University Press (OUP)","issue":"8","license":[{"start":{"date-parts":[[2020,1,6]],"date-time":"2020-01-06T00:00:00Z","timestamp":1578268800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2020,4,15]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Motivation<\/jats:title><jats:p>Longitudinal metagenomics data, including both 16S rRNA and whole-metagenome shotgun sequencing data, enhanced our abilities to understand the dynamic associations between the human microbiome and various diseases. However, analytic tools have not been fully developed to simultaneously address the main challenges of longitudinal metagenomics data, i.e. high-dimensionality, dependence among samples and zero-inflation of observed counts.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>We propose a fast zero-inflated negative binomial mixed modeling (FZINBMM) approach to analyze high-dimensional longitudinal metagenomic count data. The FZINBMM approach is based on zero-inflated negative binomial mixed models (ZINBMMs) for modeling longitudinal metagenomic count data and a fast EM-IWLS algorithm for fitting ZINBMMs. FZINBMM takes advantage of a commonly used procedure for fitting linear mixed models, which allows us to include various types of fixed and random effects and within-subject correlation structures and quickly analyze many taxa. We found that FZINBMM remarkably outperformed in computational efficiency and was statistically comparable with two R packages, GLMMadaptive and glmmTMB, that use numerical integration to fit ZINBMMs. Extensive simulations and real data applications showed that FZINBMM outperformed other previous methods, including linear mixed models, negative binomial mixed models and zero-inflated Gaussian mixed models.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and implementation<\/jats:title><jats:p>FZINBMM has been implemented in the R package NBZIMM, available in the public GitHub repository http:\/\/github.com\/\/nyiuab\/\/NBZIMM.<\/jats:p><\/jats:sec><jats:sec><jats:title>Supplementary information<\/jats:title><jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/btz973","type":"journal-article","created":{"date-parts":[[2020,1,1]],"date-time":"2020-01-01T12:09:17Z","timestamp":1577880557000},"page":"2345-2351","source":"Crossref","is-referenced-by-count":37,"title":["Fast zero-inflated negative binomial mixed modeling approach for analyzing longitudinal metagenomics data"],"prefix":"10.1093","volume":"36","author":[{"given":"Xinyan","family":"Zhang","sequence":"first","affiliation":[{"name":"Department of Biostatistics, Jiann-Ping Hsu College of Public Health, Georgia Southern University , Statesboro, GA 30458, USA"}]},{"given":"Nengjun","family":"Yi","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, University of Alabama at Birmingham , Birmingham, AL 35294, USA"}]}],"member":"286","published-online":{"date-parts":[[2020,1,6]]},"reference":[{"key":"2023013110170591900_btz973-B1","doi-asserted-by":"crossref","first-page":"18933","DOI":"10.1073\/pnas.1007028107","article-title":"Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors","volume":"107","author":"Benson","year":"2010","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023013110170591900_btz973-B2","doi-asserted-by":"crossref","first-page":"127","DOI":"10.1016\/j.tree.2008.10.008","article-title":"Generalized linear mixed models: a practical guide for ecology and evolution","volume":"24","author":"Bolker","year":"2009","journal-title":"Trends Ecol. Evol"},{"key":"2023013110170591900_btz973-B3","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. Am. Stat. Assoc"},{"key":"2023013110170591900_btz973-B4","doi-asserted-by":"crossref","first-page":"378","DOI":"10.32614\/RJ-2017-066","article-title":"glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling","volume":"9","author":"Brooks","year":"2017","journal-title":"R J"},{"key":"2023013110170591900_btz973-B5","doi-asserted-by":"crossref","first-page":"2611","DOI":"10.1093\/bioinformatics\/btw308","article-title":"A two-part mixed-effects model for analyzing longitudinal microbiome compositional data","volume":"32","author":"Chen","year":"2016","journal-title":"Bioinformatics"},{"key":"2023013110170591900_btz973-B6","doi-asserted-by":"crossref","first-page":"4131","DOI":"10.1016\/j.febslet.2014.02.037","article-title":"The dynamic microbiome","volume":"588","author":"Gerber","year":"2014","journal-title":"FEBS Lett"},{"key":"2023013110170591900_btz973-B7","doi-asserted-by":"crossref","first-page":"1323","DOI":"10.1080\/01621459.2016.1215989","article-title":"Joint selection in mixed models using regularized PQL","volume":"112","author":"Hui","year":"2017","journal-title":"J. Am. Stat. Assoc"},{"key":"2023013110170591900_btz973-B8","doi-asserted-by":"crossref","first-page":"276","DOI":"10.1016\/j.chom.2014.08.014","article-title":"The Integrative Human Microbiome Project: dynamic analysis of microbiome-host omics profiles during periods of human health and disease","volume":"16","year":"2014","journal-title":"Cell Host Microbe"},{"key":"2023013110170591900_btz973-B9","doi-asserted-by":"crossref","first-page":"459","DOI":"10.3389\/fmicb.2016.00459","article-title":"Characterization of the gut microbiome using 16S or shotgun metagenomics","volume":"7","author":"Jovel","year":"2016","journal-title":"Front. Microbiol"},{"key":"2023013110170591900_btz973-B10","doi-asserted-by":"crossref","first-page":"12522","DOI":"10.1073\/pnas.1409497111","article-title":"Patterned progression of bacterial populations in the premature infant gut","volume":"111","author":"La Rosa","year":"2014","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023013110170591900_btz973-B11","doi-asserted-by":"crossref","first-page":"552","DOI":"10.1186\/s13059-014-0552-6","article-title":"Host genetics and diet, but not immunoglobulin A expression, converge to shape compositional features of the gut microbiome in an advanced intercross population of mice","volume":"15","author":"Leamy","year":"2014","journal-title":"Genome Biol"},{"key":"2023013110170591900_btz973-B12","doi-asserted-by":"crossref","first-page":"4985","DOI":"10.1038\/srep04985","article-title":"Dynamic microbe and molecule networks in a mouse model of colitis-associated colorectal cancer","volume":"4","author":"Liang","year":"2014","journal-title":"Sci. Rep"},{"key":"2023013110170591900_btz973-B13","doi-asserted-by":"crossref","first-page":"186","DOI":"10.1016\/j.anaerobe.2011.03.016","article-title":"Interpreting the epidemiology and natural history of bacterial vaginosis: are we still confused?","volume":"17","author":"Marrazzo","year":"2011","journal-title":"Anaerobe"},{"key":"2023013110170591900_btz973-B14","volume-title":"Generalized, Linear, and Mixed Models","author":"McCulloch","year":"2001"},{"key":"2023013110170591900_btz973-B15","doi-asserted-by":"crossref","first-page":"88","DOI":"10.1111\/ppe.12106","article-title":"Early pregnancy changes in bacterial vaginosis-associated bacteria and preterm delivery","volume":"28","author":"Nelson","year":"2014","journal-title":"Paediatr. Perinat. Epidemiol"},{"key":"2023013110170591900_btz973-B16","doi-asserted-by":"crossref","first-page":"118","DOI":"10.1111\/j.2041-210X.2010.00021.x","article-title":"Do not log-transform count data","volume":"1","author":"O'Hara","year":"2010","journal-title":"Methods Ecol. Evol"},{"key":"2023013110170591900_btz973-B17","doi-asserted-by":"crossref","first-page":"75","DOI":"10.4137\/BBI.S12462","article-title":"Metagenomics: tools and insights for analyzing next-generation sequencing data derived from biodiversity studies","volume":"9","author":"Oulas","year":"2015","journal-title":"Bioinform. Biol. Insights"},{"key":"2023013110170591900_btz973-B18","doi-asserted-by":"crossref","first-page":"1023","DOI":"10.1038\/nmeth.4468","article-title":"Accessible, curated metagenomic data through ExperimentHub","volume":"14","author":"Pasolli","year":"2017","journal-title":"Nat. Methods"},{"key":"2023013110170591900_btz973-B19","doi-asserted-by":"crossref","DOI":"10.1007\/978-1-4419-0318-1","volume-title":"Mixed-Effects Models in S and S-PLUS","author":"Pinheiro","year":"2000"},{"key":"2023013110170591900_btz973-B20","doi-asserted-by":"crossref","first-page":"58","DOI":"10.1198\/106186006X96962","article-title":"Efficient Laplacian and adaptive Gaussian quadrature algorithms for multilevel generalized linear mixed models","volume":"15","author":"Pinheiro","year":"2006","journal-title":"J. Comput. Graph. Stat"},{"key":"2023013110170591900_btz973-B21","doi-asserted-by":"crossref","first-page":"507","DOI":"10.1111\/j.1467-9876.2005.00510.x","article-title":"Generalized additive models for location, scale and shape","volume":"54","author":"Rigby","year":"2005","journal-title":"J. R. Stat. Soc. C Appl"},{"key":"2023013110170591900_btz973-B22","author":"Rizopoulos","year":"2019"},{"key":"2023013110170591900_btz973-B23","doi-asserted-by":"crossref","first-page":"4","DOI":"10.1186\/2049-2618-2-4","article-title":"The composition and stability of the vaginal microbiota of normal pregnant women is different from that of non-pregnant women","volume":"2","author":"Romero","year":"2014","journal-title":"Microbiome"},{"key":"2023013110170591900_btz973-B24","doi-asserted-by":"crossref","first-page":"18","DOI":"10.1186\/2049-2618-2-18","article-title":"The vaginal microbiota of pregnant women who subsequently have spontaneous preterm labor and delivery and those with a normal delivery at term","volume":"2","author":"Romero","year":"2014","journal-title":"Microbiome"},{"key":"2023013110170591900_btz973-B25","doi-asserted-by":"crossref","first-page":"719","DOI":"10.1093\/biomet\/78.4.719","article-title":"Estimation in generalized linear models with random effects","volume":"78","author":"Schall","year":"1991","journal-title":"Biometrika"},{"key":"2023013110170591900_btz973-B26","doi-asserted-by":"crossref","first-page":"e00036","DOI":"10.1128\/mSystems.00036-17","article-title":"Development of inflammatory Bowel disease is linked to a longitudinal restructuring of the gut metagenome in mice","volume":"2","author":"Sharpton","year":"2017","journal-title":"mSystems"},{"key":"2023013110170591900_btz973-B27","doi-asserted-by":"crossref","first-page":"2462","DOI":"10.1038\/ncomms3462","article-title":"Genome-wide mapping of gene-microbiota interactions in susceptibility to autoimmune skin blistering","volume":"4","author":"Srinivas","year":"2013","journal-title":"Nat. Commun"},{"key":"2023013110170591900_btz973-B28","doi-asserted-by":"crossref","first-page":"75","DOI":"10.1186\/s40168-017-0295-1","article-title":"Longitudinal development of the gut microbiome and metabolome in preterm neonates with late onset sepsis and healthy controls","volume":"5","author":"Stewart","year":"2017","journal-title":"Microbiome"},{"key":"2023013110170591900_btz973-B29","doi-asserted-by":"crossref","first-page":"S38","DOI":"10.1111\/j.1753-4887.2012.00493.x","article-title":"Defining the human microbiome","volume":"70","author":"Ursell","year":"2012","journal-title":"Nutr. Rev"},{"key":"2023013110170591900_btz973-B30","doi-asserted-by":"crossref","first-page":"1551","DOI":"10.1016\/j.cell.2016.05.056","article-title":"Variation in microbiome LPS immunogenicity contributes to autoimmunity in humans","volume":"165","author":"Vatanen","year":"2016","journal-title":"Cell"},{"key":"2023013110170591900_btz973-B31","doi-asserted-by":"crossref","DOI":"10.1007\/978-0-387-21706-2","volume-title":"Modern Applied Statistics with S","author":"Venables","year":"2002"},{"key":"2023013110170591900_btz973-B32","doi-asserted-by":"crossref","first-page":"6440","DOI":"10.1038\/ncomms7440","article-title":"Analysis of intestinal microbiota in hybrid house mice reveals evolutionary divergence in a vertebrate hologenome","volume":"6","author":"Wang","year":"2015","journal-title":"Nat. Commun"},{"key":"2023013110170591900_btz973-B33","doi-asserted-by":"crossref","first-page":"1548","DOI":"10.1080\/01621459.2016.1180986","article-title":"Smoothing parameter and model selection for general smooth models","volume":"111","author":"Wood","year":"2016","journal-title":"J. Am. Stat. Assoc"},{"key":"2023013110170591900_btz973-B34","doi-asserted-by":"crossref","first-page":"441","DOI":"10.3389\/fnagi.2017.00441","article-title":"Longitudinal analysis of fecal microbiome and pathologic processes in a rotenone induced mice model of Parkinson's disease","volume":"9","author":"Yang","year":"2018","journal-title":"Front. Aging Neurosci"},{"key":"2023013110170591900_btz973-B35","doi-asserted-by":"crossref","first-page":"4","DOI":"10.1186\/s12859-016-1441-7","article-title":"Negative binomial mixed models for analyzing microbiome count data","volume":"18","author":"Zhang","year":"2017","journal-title":"BMC Bioinformatics"},{"key":"2023013110170591900_btz973-B36","doi-asserted-by":"crossref","first-page":"1683","DOI":"10.3389\/fmicb.2018.01683","article-title":"Negative binomial mixed models for analyzing longitudinal microbiome data","volume":"9","author":"Zhang","year":"2018","journal-title":"Front. Microbiol"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btz973\/31721162\/btz973.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/36\/8\/2345\/48984348\/bioinformatics_36_8_2345.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/36\/8\/2345\/48984348\/bioinformatics_36_8_2345.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,7,29]],"date-time":"2024-07-29T14:15:24Z","timestamp":1722262524000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/36\/8\/2345\/5697093"}},"subtitle":[],"editor":[{"given":"Alfonso","family":"Valencia","sequence":"additional","affiliation":[]}],"short-title":[],"issued":{"date-parts":[[2020,1,6]]},"references-count":36,"journal-issue":{"issue":"8","published-print":{"date-parts":[[2020,4,15]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btz973","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2020,4,15]]},"published":{"date-parts":[[2020,1,6]]}}}