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.

Advantages of phylogenetic distance based constrained ordination analyses for the examination of microbial communities

Bayesian Nonparametric Ordination for the Analysis of Microbial Communities

Application of multivariate statistical techniques in microbial ecology

A New Approach for Scalable Analysis of Microbial Communities

Deciphering Centimeter-Scale Microbial Biogeographic Patterns and Community Assembly Processes From A 3D Perspective

PhyloFunc: Phylogeny-informed Functional Distance as a New Ecological Metric for Metaproteomic Data Analysis

Phylogenetic association analysis with conditional rank correlation

Towards a Taxonomical Consensus: Diversity and Richness Inference from Large Scale rRNA gene Analysis

Use of Direct Gradient Analysis to Uncover Biological Hypotheses in 16S Survey Data and Beyond

High-resolution microbial community reconstruction by integrating short reads from multiple 16S rRNA regions

[Phylogenetic Diversity of Microbes and Its Perspectives in Conservation Biology].

Revealing the hierarchical structure of microbial communities

Transformation and differential abundance analysis of microbiome data incorporating phylogeny

Testing ecological theories with sequence similarity networks: marine ciliates exhibit similar geographic dispersal patterns as multicellular organisms

Principal Amalgamation Analysis for Microbiome Data

A Statistical Perspective on the Challenges in Molecular Microbial Biology

Molecular Techniques for Determining Microbial Diversity and Community Structure in Natural Environments

A rapid phylogeny-based method for accurate community profiling of large-scale metabarcoding datasets

Preserving Hidden Hierarchical Structure: Poincaré Distance for Enhanced Genomic Sequence Analysis

A geometric view of Biodiversity: scaling to metagenomics