Abstract:BackgroundMicrobiome studies incorporate next-generation sequencing to obtain profiles of microbial communities. Data generated from these experiments are high-dimensional with a rich correlation structure but modest sample sizes. A statistical model that utilizes these microbiome profiles to explain a clinical or biological endpoint needs to tackle high-dimensionality resulting from the very large space of variable configurations. Ensemble models are a class of approaches that can address high-dimensionality by aggregating information across large model spaces. Although such models are popular in fields as diverse as economics and genetics, their performance on microbiome data has been largely unexplored.ResultsWe developed a simulation framework that accurately captures the constraints of experimental microbiome data. Using this setup, we systematically evaluated a selection of both frequentist and Bayesian regression modeling ensembles. These are represented by variants of stability selection in conjunction with elastic net and spike-and-slab Bayesian model averaging (BMA), respectively. BMA ensembles that explore a larger space of models relative to stability selection variants performed better and had lower variability across simulations. However, stability selection ensembles were able to match the performance of BMA in scenarios of low sparsity where several variables had large regression coefficients.ConclusionsGiven a microbiome dataset of interest, we present a methodology to generate simulated data that closely mimics its characteristics in a manner that enables meaningful evaluation of analytical strategies. Our evaluation demonstrates that the largest ensembles yield the strongest performance on microbiome data with modest sample sizes and high-dimensional measurements. We also demonstrate the ability of these ensembles to identify microbiome signatures that are associated with opportunistic Candida albicans colonization during antibiotic exposure. As the focus of microbiome research evolves from pilot to translational studies, we anticipate that our strategy will aid investigators in making evaluation-based decisions for selecting appropriate analytical methods.

Latent Variable Modeling for the Microbiome

Robust and Scalable Models of Microbiome Dynamics

Latent Space Temporal Model of Microbial Abundance to Predict Domination and Bacteremia

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

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

Inference of Dynamic Regimes in the Microbiome

Statistical challenges in longitudinal microbiome data analysis

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

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

Negative Binomial factor regression with application to microbiome data analysis

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

Bayesian Modeling of Microbiome Data for Differential Abundance Analysis

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

Random Graphical Model of Microbiome Interactions in Related Environments

A systematic evaluation of high-dimensional, ensemble-based regression for exploring large model spaces in microbiome analyses

Bayesian Nonparametric Ordination for the Analysis of Microbial Communities

Bayesian modeling of co-occurrence microbial interaction networks

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

Microbiome subcommunity learning with logistic-tree normal latent Dirichlet allocation

mbtransfer: Microbiome intervention analysis using transfer functions and mirror statistics

Dirichlet distribution parameter estimation with applications in microbiome analyses