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.

Correction for Regueira-Iglesias et al., "Impact of 16S rRNA Gene Redundancy and Primer Pair Selection on the Quantification and Classification of Oral Microbiota in Next-Generation Sequencing"

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

Correction for Morón-López et al., "Comparison of Reverse Transcription (RT)-Quantitative PCR and RT-Droplet Digital PCR for Detection of Genomic and Subgenomic SARS-CoV-2 RNA"

16S rRNA gene sequencing of mock microbial populations- impact of DNA extraction method, primer choice and sequencing platform

Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies

In Silico and Experimental Evaluation of Primer Sets for Species-Level Resolution of the Vaginal Microbiota Using 16S Ribosomal RNA Gene Sequencing

Optimized bacterial community characterization through full-length 16S rRNA gene sequencing utilizing MinION nanopore technology

A 16S rDNA PCR-based theoretical to actual delta approach on culturable mock communities revealed severe losses of diversity information

Improving Species Level‐taxonomic Assignment from 16S rRNA Sequencing Technologies

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

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

Performance of a Shotgun Prediction Model for Colorectal Cancer When Using 16S rRNA Sequencing Data

Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys

Comparison between 16S rRNA and shotgun sequencing in colorectal cancer, advanced colorectal lesions, and healthy human gut microbiota

Accurate quantitation of 16S gene copies in low biomass samples post-antibiotic treatment through deep sequencing with a balanced nucleotide synthetic spike-in approach

Targeting the 16S rRNA Gene for Bacterial Identification in Complex Mixed Samples: Comparative Evaluation of Second (Illumina) and Third (Oxford Nanopore Technologies) Generation Sequencing Technologies

The Influences of Bioinformatics Tools and Reference Databases in Analyzing the Human Oral Microbial Community

An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform

The effect of taxonomic classification by full-length 16S rRNA sequencing with a synthetic long-read technology

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

Equivolumetric Protocol Generates Library Sizes Proportional to Total Microbial Load in 16S Amplicon Sequencing