Abstract:One of the major methods to identify microbial community composition, to unravel microbial population dynamics, and to explore microbial diversity in environmental samples is high-throughput DNA- or RNA-based 16S rRNA (gene) amplicon sequencing in combination with bioinformatics analyses. However, focusing on environmental samples from contrasting habitats, it was not systematically evaluated (i) which analysis methods provide results that reflect reality most accurately, (ii) how the interpretations of microbial community studies are biased by different analysis methods and (iii) if the most optimal analysis workflow can be implemented in an easy-to-use pipeline. Here, we compared the performance of 16S rRNA (gene) amplicon sequencing analysis tools (i.e., Mothur, QIIME1, QIIME2, and MEGAN) using three mock datasets with known microbial community composition that differed in sequencing quality, species number and abundance distribution (i.e., even or uneven), and phylogenetic diversity (i.e., closely related or well-separated amplicon sequences). Our results showed that QIIME2 outcompeted all other investigated tools in sequence recovery (>10 times fewer false positives), taxonomic assignments (>22% better F-score) and diversity estimates (>5% better assessment), suggesting that this approach is able to reflect the in situ microbial community most accurately. Further analysis of 24 environmental datasets obtained from four contrasting terrestrial and freshwater sites revealed dramatic differences in the resulting microbial community composition for all pipelines at genus level. For instance, at the investigated river water sites Sphaerotilus was only reported when using QIIME1 (8% abundance) and Agitococcus with QIIME1 or QIIME2 (2 or 3% abundance, respectively), but both genera remained undetected when analyzed with Mothur or MEGAN. Since these abundant taxa probably have implications for important biogeochemical cycles (e.g., nitrate and sulfate reduction) at these sites, their detection and semi-quantitative enumeration is crucial for valid interpretations. A high-performance computing conformant workflow was constructed to allow FAIR (Findable, Accessible, Interoperable, and Re-usable) 16S rRNA (gene) amplicon sequence analysis starting from raw sequence files, using the most optimal methods identified in our study. Our presented workflow should be considered for future studies, thereby facilitating the analysis of high-throughput 16S rRNA (gene) sequencing data substantially, while maximizing reliability and confidence in microbial community data analysis.

Identificación de Microorganismos en Muestras Ambientales: Análisis Bioinformático del Gen 16S RRNA Mediante QIIME2

A step-by-step procedure for analysing the 16S rRNA-based microbiome diversity using QIIME 2 and comprehensive PICRUSt2 illustration for functional prediction

Interpretations of Environmental Microbial Community Studies Are Biased by the Selected 16S rRNA (Gene) Amplicon Sequencing Pipeline

Benchmarking taxonomic assignments based on 16S rRNA gene profiling of the microbiota from commonly sampled environments

q2-metnet: QIIME2 package to analyze 16S rRNA data via high-quality metabolic reconstructions of the human gut microbiota

Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2

Comparison of Mothur and QIIME for the Analysis of Rumen Microbiota Composition Based on 16S rRNA Amplicon Sequences

Microbiome data analysis with applications to pre-clinical studies using QIIME2: Statistical considerations

Comparing bioinformatic pipelines for microbial 16S rRNA amplicon sequencing

Multi-amplicon microbiome data analysis pipelines for mixed orientation sequences using QIIME2: Assessing reference database, variable region and pre-processing bias in classification of mock bacterial community samples

Automation of QIIME2 Metagenomic Analysis Platform

A multi-amplicon 16S rRNA sequencing and analysis method for improved taxonomic profiling of bacterial communities

From reads to operational taxonomic units: an ensemble processing pipeline for MiSeq amplicon sequencing data

Using Mothur to Determine Bacterial Community Composition and Structure in 16S Ribosomal RNA Datasets

Snaq: A Dynamic Snakemake Pipeline for Microbiome Data Analysis With QIIME2

A Framework for Assessing 16S Rrna Marker-Gene Survey Data Analysis Methods Using Mixtures.

High-throughput identification and quantification of bacterial cells in the microbiota based on 16S rRNA sequencing with single-base accuracy using BarBIQ

Metagenomic profiling pipelines improve taxonomic classification for 16S amplicon sequencing data

MicrobiomePhylo: A New Tool for Metabarcoding Data Downstream Analysis - A Real-World Data Analysis Demonstration

Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis

RiboSnake - a user-friendly, robust, reproducible, multipurpose and documentation-extensive pipeline for 16S rRNA gene microbiome analysis