{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,7,18]],"date-time":"2026-07-18T16:43:10Z","timestamp":1784392990925,"version":"3.55.0"},"reference-count":29,"publisher":"Oxford University Press (OUP)","issue":"8","license":[{"start":{"date-parts":[[2023,8,1]],"date-time":"2023-08-01T00:00:00Z","timestamp":1690848000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"National Institutes of Health gr","award":["5P30AG072959-02"],"award-info":[{"award-number":["5P30AG072959-02"]}]},{"name":"National Institutes of Health gr","award":["3R01DK042191-30S1"],"award-info":[{"award-number":["3R01DK042191-30S1"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,8,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Summary<\/jats:title>\n                  <jats:p>Microbiome research is now moving beyond the compositional analysis of microbial taxa in a sample. Increasing evidence from large human microbiome studies suggests that functional consequences of changes in the intestinal microbiome may provide more power for studying their impact on inflammation and immune responses. Although 16S rRNA analysis is one of the most popular and a cost-effective method to profile the microbial compositions, marker-gene sequencing cannot provide direct information about the functional genes that are present in the genomes of community members. Bioinformatic tools have been developed to predict microbiome function with 16S rRNA gene data. Among them, PICRUSt2 (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) has become one of the most popular functional profile prediction tools, which generates community-wide pathway abundances. However, no state-of-art inference tools are available to test the differences in pathway abundances between comparison groups. We have developed ggpicrust2, an R package, for analyzing functional profiles derived from 16S rRNA sequencing. This powerful tool enables researchers to conduct extensive differential abundance analyses and generate visually appealing visualizations that effectively highlight functional signals. With ggpicrust2, users can obtain publishable results and gain deeper insights into the functional composition of their microbial communities.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availability and implementation<\/jats:title>\n                  <jats:p>The package is open-source under the MIT and file license and is available at CRAN and https:\/\/github.com\/cafferychen777\/ggpicrust2. Its shiny web is available at https:\/\/a95dps-caffery-chen.shinyapps.io\/ggpicrust2_shiny\/.<\/jats:p>\n               <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btad470","type":"journal-article","created":{"date-parts":[[2023,8,1]],"date-time":"2023-08-01T16:23:46Z","timestamp":1690907026000},"source":"Crossref","is-referenced-by-count":264,"title":["ggpicrust2: an R package for PICRUSt2 predicted functional profile analysis and visualization"],"prefix":"10.1093","volume":"39","author":[{"ORCID":"https:\/\/orcid.org\/0009-0009-9151-2901","authenticated-orcid":false,"given":"Chen","family":"Yang","sequence":"first","affiliation":[{"name":"Department of Biostatistics, Southern Medical University , Guangzhou 510515, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jiahao","family":"Mai","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, Southern Medical University , Guangzhou 510515, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xuan","family":"Cao","sequence":"additional","affiliation":[{"name":"Department of Mathematical Sciences, University of Cincinnati , Cincinnati, OH 45221, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Aaron","family":"Burberry","sequence":"additional","affiliation":[{"name":"Department of Pathology, School of Medicine, Case Western Reserve University , Cleveland, OH 44106, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-1571-1548","authenticated-orcid":false,"given":"Fabio","family":"Cominelli","sequence":"additional","affiliation":[{"name":"Department of Pathology, School of Medicine, Case Western Reserve University , Cleveland, OH 44106, United States"},{"name":"Case Digestive Health Research Institute, Case Western Reserve University , Cleveland, OH 44016, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-5986-2260","authenticated-orcid":false,"given":"Liangliang","family":"Zhang","sequence":"additional","affiliation":[{"name":"Department of Population and Quantitative Health Sciences, Case Western Reserve University , Cleveland, OH 44106, United States"},{"name":"Case Comprehensive Cancer Center , Cleveland, OH 44106, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2023,8,1]]},"reference":[{"key":"2023081420342609000_btad470-B1","doi-asserted-by":"crossref","first-page":"89","DOI":"10.1038\/s41586-020-2288-7","article-title":"C9orf72 suppresses systemic and neural inflammation induced by gut bacteria","volume":"582","author":"Burberry","year":"2020","journal-title":"Nature"},{"key":"2023081420342609000_btad470-B2","doi-asserted-by":"crossref","first-page":"191","DOI":"10.1186\/s13059-020-02104-1","article-title":"Assessment of statistical methods from single cell, bulk RNA-seq, and metagenomics applied to microbiome data","volume":"21","author":"Calgaro","year":"2020","journal-title":"Genome Biol"},{"key":"2023081420342609000_btad470-B3","doi-asserted-by":"crossref","first-page":"4027","DOI":"10.1093\/bioinformatics\/btac438","article-title":"microbiomeMarker: an R\/bioconductor package for microbiome marker identification and visualization","volume":"38","author":"Cao","year":"2022","journal-title":"Bioinformatics"},{"key":"2023081420342609000_btad470-B4","doi-asserted-by":"crossref","first-page":"799","DOI":"10.1038\/s41596-019-0264-1","article-title":"Using MicrobiomeAnalyst for comprehensive statistical, functional, and meta-analysis of microbiome data","volume":"15","author":"Chong","year":"2020","journal-title":"Nat Protoc"},{"key":"2023081420342609000_btad470-B5","doi-asserted-by":"crossref","first-page":"685","DOI":"10.1038\/s41587-020-0548-6","article-title":"PICRUSt2 for prediction of metagenome functions","volume":"38","author":"Douglas","year":"2020","journal-title":"Nat Biotechnol"},{"key":"2023081420342609000_btad470-B6","doi-asserted-by":"crossref","first-page":"e67019","DOI":"10.1371\/journal.pone.0067019","article-title":"ANOVA-Like differential expression (ALDEx) analysis for mixed population RNA-Seq","volume":"8","author":"Fernandes","year":"2013","journal-title":"PLoS One"},{"key":"2023081420342609000_btad470-B7","doi-asserted-by":"crossref","first-page":"66","DOI":"10.1007\/978-1-4612-4380-9_6","volume-title":"Breakthroughs in Statistics","author":"Fisher","year":"1992"},{"key":"2023081420342609000_btad470-B8","doi-asserted-by":"crossref","first-page":"D587","DOI":"10.1093\/nar\/gkac963","article-title":"KEGG for taxonomy-based analysis of pathways and genomes","volume":"51","author":"Kanehisa","year":"2023","journal-title":"Nucleic Acids Res"},{"key":"2023081420342609000_btad470-B9","doi-asserted-by":"crossref","first-page":"2114","DOI":"10.3389\/fmicb.2017.02114","article-title":"Analysis of microbiome data in the presence of excess zeros","volume":"8","author":"Kaul","year":"2017","journal-title":"Front Microbiol"},{"key":"2023081420342609000_btad470-B10","doi-asserted-by":"crossref","first-page":"583","DOI":"10.1080\/01621459.1952.10483441","article-title":"Use of ranks in one-criterion variance analysis","volume":"47","author":"Kruskal","year":"1952","journal-title":"J Amer Stat Assoc"},{"key":"2023081420342609000_btad470-B11","doi-asserted-by":"crossref","first-page":"814","DOI":"10.1038\/nbt.2676","article-title":"Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences","volume":"31","author":"Langille","year":"2013","journal-title":"Nat Biotechnol"},{"key":"2023081420342609000_btad470-B12","doi-asserted-by":"crossref","first-page":"3514","DOI":"10.1038\/s41467-020-17041-7","article-title":"Analysis of compositions of microbiomes with bias correction","volume":"11","author":"Lin","year":"2020","journal-title":"Nat Commun"},{"key":"2023081420342609000_btad470-B13","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":"2023081420342609000_btad470-B14","doi-asserted-by":"crossref","first-page":"e1009442","DOI":"10.1371\/journal.pcbi.1009442","article-title":"Multivariable association discovery in population-scale meta-omics studies","volume":"17","author":"Mallick","year":"2021","journal-title":"PLoS Comput Biol"},{"issue":"1","key":"2023081420342609000_btad470-B15","article-title":"Analysis of composition of microbiomes: a novel method for studying microbial composition","volume":"26","author":"Mandal","year":"2015","journal-title":"Microb Ecol Health Dis"},{"key":"2023081420342609000_btad470-B16","doi-asserted-by":"crossref","first-page":"3635","DOI":"10.1016\/j.ygeno.2021.08.016","article-title":"MicFunPred: a conserved approach to predict functional profiles from 16S rRNA gene sequence data","volume":"113","author":"Mongad","year":"2021","journal-title":"Genomics"},{"key":"2023081420342609000_btad470-B17","doi-asserted-by":"crossref","first-page":"342","DOI":"10.1038\/s41467-022-28034-z","article-title":"Microbiome differential abundance methods produce different results across 38 datasets","volume":"13","author":"Nearing","year":"2022","journal-title":"Nat Commun"},{"key":"2023081420342609000_btad470-B18","doi-asserted-by":"crossref","first-page":"3123","DOI":"10.1093\/bioinformatics\/btu494","article-title":"STAMP: statistical analysis of taxonomic and functional profiles","volume":"30","author":"Parks","year":"2014","journal-title":"Bioinformatics"},{"key":"2023081420342609000_btad470-B19","doi-asserted-by":"crossref","first-page":"1200","DOI":"10.1038\/nmeth.2658","article-title":"Differential abundance analysis for microbial marker-gene surveys","volume":"10","author":"Paulson","year":"2013","journal-title":"Nat Methods"},{"key":"2023081420342609000_btad470-B20","doi-asserted-by":"crossref","first-page":"1317","DOI":"10.1093\/bioinformatics\/btaa838","article-title":"MicrobiomeExplorer\u202f: an R package for the analysis and visualization of microbial communities","volume":"37","author":"Reeder","year":"2021","journal-title":"Bioinformatics"},{"key":"2023081420342609000_btad470-B21","doi-asserted-by":"crossref","first-page":"e47","DOI":"10.1093\/nar\/gkv007","article-title":"Limma powers differential expression analyses for RNA-sequencing and microarray studies","volume":"43","author":"Ritchie","year":"2015","journal-title":"Nucleic Acids Res"},{"key":"2023081420342609000_btad470-B22","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":"2023081420342609000_btad470-B23","doi-asserted-by":"crossref","first-page":"R60","DOI":"10.1186\/gb-2011-12-6-r60","article-title":"Metagenomic biomarker discovery and explanation","volume":"12","author":"Segata","year":"2011","journal-title":"Genome Biol"},{"key":"2023081420342609000_btad470-B24","doi-asserted-by":"crossref","first-page":"1","DOI":"10.2307\/2331554","article-title":"The probable error of a mean","volume":"6","author":"Student","year":"1908","journal-title":"Biometrika"},{"key":"2023081420342609000_btad470-B25","doi-asserted-by":"crossref","first-page":"350","DOI":"10.1093\/biomet\/29.3-4.350","article-title":"The significance of the difference between two means when the population variances are unequal","volume":"29","author":"Welch","year":"1938","journal-title":"Biometrika"},{"key":"2023081420342609000_btad470-B26","doi-asserted-by":"crossref","first-page":"11","DOI":"10.1186\/s40793-020-00358-7","article-title":"Tax4Fun2: prediction of habitat-specific functional profiles and functional redundancy based on 16S rRNA gene sequences","volume":"15","author":"Wemheuer","year":"2020","journal-title":"Environ Microbiome"},{"key":"2023081420342609000_btad470-B27","doi-asserted-by":"crossref","first-page":"817","DOI":"10.2307\/1912934","article-title":"A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity","volume":"48","author":"White","year":"1980","journal-title":"Econometrica"},{"key":"2023081420342609000_btad470-B28","doi-asserted-by":"crossref","first-page":"130","DOI":"10.1186\/s40168-022-01320-0","article-title":"A comprehensive evaluation of microbial differential abundance analysis methods: current status and potential solutions","volume":"10","author":"Yang","year":"2022","journal-title":"Microbiome"},{"key":"2023081420342609000_btad470-B29","doi-asserted-by":"crossref","first-page":"95","DOI":"10.1186\/s13059-022-02655-5","article-title":"LinDA: linear models for differential abundance analysis of microbiome compositional data","volume":"23","author":"Zhou","year":"2022","journal-title":"Genome Biol"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btad470\/51016576\/btad470.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/8\/btad470\/51111779\/btad470.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/8\/btad470\/51111779\/btad470.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,8,14]],"date-time":"2023-08-14T20:59:38Z","timestamp":1692046778000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btad470\/7234609"}},"subtitle":[],"editor":[{"given":"Arne","family":"Elofsson","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2023,8,1]]},"references-count":29,"journal-issue":{"issue":"8","published-print":{"date-parts":[[2023,8,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btad470","relation":{},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,8,1]]},"published":{"date-parts":[[2023,8,1]]},"article-number":"btad470"}}