Abstract:Studying the human microbiome has gained substantial interest in recent years, and a common task in the analysis of these data is to cluster microbiome compositions into subtypes. This subdivision of samples into subgroups serves as an intermediary step in achieving personalized diagnosis and treatment. In applying existing clustering methods to modern microbiome studies including the American Gut Project (AGP) data, we found that this seemingly standard task, however, is very challenging in the microbiome composition context due to several key features of such data. Standard distance-based clustering algorithms generally do not produce reliable results as they do not take into account the heterogeneity of the cross-sample variability among the bacterial taxa, while existing model-based approaches do not allow sufficient flexibility for the identification of complex within-cluster variation from cross-cluster variation. Direct applications of such methods generally lead to overly dispersed clusters in the AGP data and such phenomenon is common for other microbiome data. To overcome these challenges, we introduce Dirichlet-tree multinomial mixtures (DTMM) as a Bayesian generative model for clustering amplicon sequencing data in microbiome studies. DTMM models the microbiome population with a mixture of Dirichlet-tree kernels that utilizes the phylogenetic tree to offer a more flexible covariance structure in characterizing within-cluster variation, and it provides a means for identifying a subset of signature taxa that distinguish the clusters. We perform extensive simulation studies to evaluate the performance of DTMM and compare it to state-of-the-art model-based and distance-based clustering methods in the microbiome context, and carry out a validation study on a publicly available longitudinal data set to confirm the biological relevance of the clusters. Finally, we report a case study on the fecal data from the AGP to identify compositional clusters among individuals with inflammatory bowel disease and diabetes. Among our most interesting findings is that enterotypes (i.e., gut microbiome clusters) are not always defined by the most dominant species as previous analyses had assumed, but can involve a number of less abundant OTUs, which cannot be identified with existing distance-based and method-based approaches.

Poisson hurdle model-based method for clustering microbiome features

Dirichlet-tree multinomial mixtures for clustering microbiome compositions

Bayesian biclustering for microbial metagenomic sequencing data via multinomial matrix factorization

Sparse tree-based clustering of microbiome data to characterize microbiome heterogeneity in pancreatic cancer

Bayesian Modeling of Microbiome Data for Differential Abundance Analysis

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Clustering on Human Microbiome Sequencing Data: A Distance-Based Unsupervised Learning Model

Cluster analysis of microbiome data by using mixtures of Dirichlet–multinomial regression models

MSClust: A Multi-Seeds based Clustering algorithm for microbiome profiling using 16S rRNA sequence

A Multivariate Poisson-Log Normal Mixture Model for Clustering Transcriptome Sequencing Data

Modeling Sparse Data Using MLE with Applications to Microbiome Data

The Poisson distribution model fits UMI-based single-cell RNA-sequencing data

A Bi-Poisson Model for Clustering Gene Expression Profiles by RNA-seq.

A Flexible Zero-Inflated Poisson-Gamma model with application to microbiome read counts

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

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

DIMM-SC: a Dirichlet mixture model for clustering droplet-based single cell transcriptomic data

A New Unsupervised Binning Approachfor Metagenomic Sequences Based onN-grams and Automatic Feature Weighting.

A semi-parametric multiple imputation method for high-sparse, high-dimensional, compositional data

Accurate Reconstruction of Microbial Strains from Metagenomic Sequencing Using Representative Reference Genomes

Variable selection for sparse Dirichlet-multinomial regression with an application to microbiome data analysis