Abstract:Targeted amplicon sequencing is widely used in microbial ecology studies. However, sequencing artifacts and amplification biases are of great concern. To identify sources of these artifacts, a systematic analysis was performed using mock communities comprised of 16S rRNA genes from 33 bacterial strains. Our results indicated that while sequencing errors were generally isolated to low-abundance operational taxonomic units, chimeric sequences were a major source of artifacts. Singleton and doubleton sequences were primarily chimeras. Formation of chimeric sequences was significantly correlated with the GC content of the targeted sequences. Low-GC-content mock community members exhibited lower rates of chimeric sequence formation. GC content also had a large impact on sequence recovery. The quantitative capacity was notably limited, with substantial recovery variations and weak correlation between anticipated and observed strain abundances. The mock community strains with higher GC content had higher recovery rates than strains with lower GC content. Amplification bias was also observed due to the differences in primer affinity. A two-step PCR strategy reduced the number of chimeric sequences by half. In addition, comparative analyses based on the mock communities showed that several widely used sequence processing pipelines/methods, including DADA2, Deblur, UCLUST, UNOISE, and UPARSE, had different advantages and disadvantages in artifact removal and rare species detection. These results are important for improving sequencing quality and reliability and developing new algorithms to process targeted amplicon sequences.IMPORTANCEAmplicon sequencing of targeted genes is the predominant approach to estimate the membership and structure of microbial communities. However, accurate reconstruction of community composition is difficult due to sequencing errors, and other methodological biases and effective approaches to overcome these challenges are essential. Using a mock community of 33 phylogenetically diverse strains, this study evaluated the effect of GC content on sequencing results and tested different approaches to improve overall sequencing accuracy while characterizing the pros and cons of popular amplicon sequence data processing approaches. The sequencing results from this study can serve as a benchmarking data set for future algorithmic improvements. Furthermore, the new insights on sequencing error, chimera formation, and GC bias from this study will help enhance the quality of amplicon sequencing studies and support the development of new data analysis approaches. Amplicon sequencing of targeted genes is the predominant approach to estimate the membership and structure of microbial communities. However, accurate reconstruction of community composition is difficult due to sequencing errors, and other methodological biases and effective approaches to overcome these challenges are essential. Using a mock community of 33 phylogenetically diverse strains, this study evaluated the effect of GC content on sequencing results and tested different approaches to improve overall sequencing accuracy while characterizing the pros and cons of popular amplicon sequence data processing approaches. The sequencing results from this study can serve as a benchmarking data set for future algorithmic improvements. Furthermore, the new insights on sequencing error, chimera formation, and GC bias from this study will help enhance the quality of amplicon sequencing studies and support the development of new data analysis approaches.

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

Application of Next Generation Sequencing Technology on Contamination Monitoring in Microbiology Laboratory

Diverse and widespread contamination evident in the unmapped depths of high throughput sequencing data

Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data

Identification and Mitigation of microbial contaminant reads from long-read metagenomic samples

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

Recognizing the reagent microbiome

Identification and Removal of Contaminant Sequences From Ribosomal Gene Databases: Lessons From the Census of Deep Life

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

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

Toward Standards in Clinical Microbiota Studies: Comparison of Three DNA Extraction Methods and Two Bioinformatic Pipelines

Clearing the plate: a strategic approach to mitigate well-to-well contamination in large-scale microbiome studies

The impact of sequencing depth on the inferred taxonomic composition and AMR gene content of metagenomic samples

Comparison of commercial DNA extraction kits for whole metagenome sequencing of human oral, vaginal, and rectal microbiome samples

The impact of kit, environment, and sampling contamination on the observed microbiome of bovine milk

Unveiling errors in soil microbial community sequencing: a case for reference soils and improved diagnostics for nanopore sequencing

Variation in the metagenomic analysis of fecal microbiome composition calls for a standardized operating approach

Sequencing Introduced False Positive Rare Taxa Lead to Biased Microbial Community Diversity, Assembly, and Interaction Interpretation in Amplicon Studies

Small amounts of misassembly can have disproportionate effects on pangenome-based metagenomic analyses

Effects of Error, Chimera, Bias, and GC Content on the Accuracy of Amplicon Sequencing.

In silico analyses identify sequence contamination thresholds for Nanopore-generated SARS-CoV-2 sequences