Abstract:Abstract Motivation Microbiome data have proven extremely useful for understanding microbial communities and their impacts in health and disease. Although microbiome analysis methods and standards are evolving rapidly, obtaining meaningful and interpretable results from microbiome studies still requires careful statistical treatment. In particular, many existing and emerging methods for differential abundance (DA) analysis fail to account for the fact that microbiome data are high-dimensional and sparse, compositional, negatively and positively correlated and phylogenetically structured. To better describe microbiome data and improve the power of DA testing, there is still a great need for the continued development of appropriate statistical methodology. Results In this article, we propose a model-based approach for microbiome data transformation, and a phylogenetically informed procedure for DA testing based on the transformed data. First, we extend the Dirichlet-tree multinomial (DTM) to zero-inflated DTM for multivariate modeling of microbial counts, addressing data sparsity and correlation and phylogeny among bacterial taxa. Then, within this framework and using a Bayesian formulation, we introduce posterior mean transformation to convert raw counts into non-zero relative abundances that sum to one, accounting for the compositionality nature of microbiome data. Second, using the transformed data, we propose adaptive analysis of composition of microbiomes (adaANCOM) for DA testing by constructing log-ratios adaptively on the tree for each taxon, greatly reducing the computational complexity of ANCOM in high dimensions. Finally, we present extensive simulation studies, an analysis of HMP data across 18 body sites and 2 visits, and an application to a gut microbiome and malnutrition study, to investigate the performance of posterior mean transformation and adaANCOM. Comparisons with ANCOM and other DA testing procedures show that adaANCOM controls the false discovery rate well, allows for easy interpretation of the results, and is computationally efficient for high-dimensional problems. Availability and implementation The developed R package is available at https://github.com/ZRChao/adaANCOM. For replicability purposes, scripts for our simulations and data analysis are available at https://github.com/ZRChao/Papers_supplementary. Supplementary information Supplementary data are available at Bioinformatics online.

A Bayesian Nonparametric Approach for Identifying Differentially Abundant Taxa in Multigroup Microbiome Data with Covariates

Bayesian Modeling of Microbiome Data for Differential Abundance Analysis

A Bayesian Zero-Inflated Negative Binomial Regression Model for the Integrative Analysis of Microbiome Data

Testing for differential abundance in compositional counts data, with application to microbiome studies

An integrative Bayesian Dirichlet-multinomial regression model for the analysis of taxonomic abundances in microbiome data

Bayesian Mixed Effects Models for Zero-inflated Compositions in Microbiome Data Analysis

A Zero‐inflated Beta‐binomial Model for Microbiome Data Analysis

Zero is not absence: censoring-based differential abundance analysis for microbiome data

An Empirical Bayes Approach to Normalization and Differential Abundance Testing for Microbiome Data

A Survey of Statistical Methods for Microbiome Data Analysis

Association of Body Index with Fecal Microbiome in Children Cohorts with Ethnic-Geographic Factor Interaction: Accurately Using a Bayesian Zero-inflated Negative Binomial Regression Model

Powerful and robust non-parametric association testing for microbiome data via a zero-inflated quantile approach (ZINQ)

Bayesian Variable Selection for Multivariate Zero-Inflated Models: Application to Microbiome Count Data

Bayesian biclustering for microbial metagenomic sequencing data via multinomial matrix factorization

MZINBVA: Variational Approximation for Multilevel Zero-Inflated Negative-Binomial Models for Association Analysis in Microbiome Surveys

Negative Binomial factor regression with application to microbiome data analysis

Transformation and differential abundance analysis of microbiome data incorporating phylogeny

Conditional regression based on a multivariate zero-inflated logistic normal model for microbiome relative abundance data

A two-part mixed-effects model for analyzing longitudinal microbiome compositional data

A strategy for differential abundance analysis of sparse microbiome data with group-wise structured zeros

A nonparametric spatial test to identify factors that shape a microbiome