Abstract:Short-amplicon 16S rRNA gene sequencing is currently the method of choice for studies investigating microbiomes. However, comparative studies on differences in procedures are scarce. We sequenced human stool samples and mock communities with increasing complexity using a variety of commonly used protocols. Short amplicons targeting different variable regions (V-regions) or ranges thereof (V1-V2, V1-V3, V3-V4, V4, V4-V5, V6-V8, and V7-V9) were investigated for differences in the composition outcome due to primer choices. Next, the influence of clustering (operational taxonomic units [OTUs], zero-radius OTUs [zOTUs], and amplicon sequence variants [ASVs]), different databases (GreenGenes, the Ribosomal Database Project, Silva, the genomic-based 16S rRNA Database, and The All-Species Living Tree), and bioinformatic settings on taxonomic assignment were also investigated. We present a systematic comparison across all typically used V-regions using well-established primers. While it is known that the primer choice has a significant influence on the resulting microbial composition, we show that microbial profiles generated using different primer pairs need independent validation of performance. Further, comparing data sets across V-regions using different databases might be misleading due to differences in nomenclature (e.g., Enterorhabdus versus Adlercreutzia ) and varying precisions in classification down to genus level. Overall, specific but important taxa are not picked up by certain primer pairs (e.g., Bacteroidetes is missed using primers 515F-944R) or due to the database used (e.g., Acetatifactor in GreenGenes and the genomic-based 16S rRNA Database). We found that appropriate truncation of amplicons is essential and different truncated-length combinations should be tested for each study. Finally, specific mock communities of sufficient and adequate complexity are highly recommended. IMPORTANCE In 16S rRNA gene sequencing, certain bacterial genera were found to be underrepresented or even missing in taxonomic profiles when using unsuitable primer combinations, outdated reference databases, or inadequate pipeline settings. Concerning the last, quality thresholds as well as bioinformatic settings (i.e., clustering approach, analysis pipeline, and specific adjustments such as truncation) are responsible for a number of observed differences between studies. Conclusions drawn by comparing one data set to another (e.g., between publications) appear to be problematic and require independent cross-validation using matching V-regions and uniform data processing. Therefore, we highlight the importance of a thought-out study design including sufficiently complex mock standards and appropriate V-region choice for the sample of interest. The use of processing pipelines and parameters must be tested beforehand.

Epistemic misalignments in microbiome research

Identifying biases and their potential solutions in human microbiome studies

Putting science back into microbial ecology: a question of approach

Variability and bias in microbiome metagenomic sequencing: an interlaboratory study comparing experimental protocols

Reagent and laboratory contamination can critically impact sequence-based microbiome analyses

Insights into microbiome research 2: Experimental design, sample collection, and shipment

Microbiomes Through The Looking Glass

Microbiome definition re-visited: old concepts and new challenges

Assessment of statistical methods from single cell, bulk RNA-seq, and metagenomics applied to microbiome data

Emergent Ecological Patterns and Modelling of Gut Microbiomes in Health and in Disease

Emerging Priorities for Microbiome Research

A host of issues: pseudoreplication in host–microbiota studies

Evaluating the accuracy of amplicon-based microbiome computational pipelines on simulated human gut microbial communities

Characterizing Microbiomes via Sequencing of Marker Loci: Techniques To Improve Throughput, Account for Cross-Contamination, and Reduce Cost

Human microbiome myths and misconceptions

Primer, Pipelines, Parameters: Issues in 16S rRNA Gene Sequencing

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

Accounting for reciprocal host–microbiome interactions in experimental science

Integrating biological knowledge for mechanistic inference in the host-associated microbiome

Networks as tools for defining emergent properties of microbiomes and their stability