Abstract:Microorganisms play critical roles in human health and disease. It is well known that microbes live in diverse communities in which they interact synergistically or antagonistically. Thus for estimating microbial associations with clinical covariates, multivariate statistical models are preferred. Multivariate models allow one to estimate and exploit complex interdependencies among multiple taxa, yielding more powerful tests of exposure or treatment effects than application of taxon-specific univariate analyses. In addition, the analysis of microbial count data requires special attention because data commonly exhibit zero inflation. To meet these needs, we developed a Bayesian variable selection model for multivariate count data with excess zeros that incorporates information on the covariance structure of the outcomes (counts for multiple taxa), while estimating associations with the mean levels of these outcomes. Although there has been a great deal of effort in zero-inflated models for longitudinal data, little attention has been given to high-dimensional multivariate zero-inflated data modeled via a general correlation structure. Through simulation, we compared performance of the proposed method to that of existing univariate approaches, for both the binary and count parts of the model. When outcomes were correlated the proposed variable selection method maintained type I error while boosting the ability to identify true associations in the binary component of the model. For the count part of the model, in some scenarios the the univariate method had higher power than the multivariate approach. This higher power was at a cost of a highly inflated false discovery rate not observed with the proposed multivariate method. We applied the approach to oral microbiome data from the Pediatric HIV/AIDS Cohort Oral Health Study and identified five species (of 44) associated with HIV infection.

Model-Based Microbiome Data Ordination: A Variational Approximation Approach.

A Zero-Inflated Logistic Normal Multinomial Model for Extracting Microbial Compositions

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

Latent Variable Modeling for the Microbiome

Variational inference for microbiome survey data with application to global ocean data

Bayesian Nonparametric Ordination for the Analysis of Microbial Communities

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

Random Graphical Model of Microbiome Interactions in Related Environments

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

Mbdenoise: Microbiome Data Denoising Using Zero-Inflated Probabilistic Principal Components Analysis

Zero-inflated Poisson Factor Model with Application to Microbiome Absolute Abundance Data

Modelling Microbiome Association with Host Phenotypes Using a Bayesian Dirichlet Process Model

Bayesian compositional regression with microbiome features via variational inference

A Bayesian model of microbiome data for simultaneous identification of covariate associations and prediction of phenotypic outcomes

Applications and Comparison of Dimensionality Reduction Methods for Microbiome Data

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

It's All Relative: New Regression Paradigm for Microbiome Compositional Data

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

Negative Binomial factor regression with application to microbiome data analysis

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

Bayesian Modeling of Microbiome Data for Differential Abundance Analysis