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.

Analysis of Barcode sequence features to find anomalies due to amplification Bias

Too many needles in this haystack: algorithms for the analysis of next generation sequence data

Improving the Accuracy of Bulk Fitness Assays by Correcting Barcode Processing Biases

Highly Scalable Algorithms for Robust String Barcoding

Single Nucleotide Recognition and Mutation Site Sequencing Based on a Barcode Assay and Rolling Circle Amplification

Pacybara: accurate long-read sequencing for barcoded mutagenized allelic libraries

DNA Barcodes using a Double Nanopore System

The efficacy and further functional advantages of random-base molecular barcodes for absolute and digital quantification of nucleic acid molecules

Counting unique molecular identifiers in sequencing using a multitype branching process with immigration

Sequencing-induced artefacts in NGS STR data

Improving Detection Of Copy-Number Variation By Simultaneous Bias Correction And Read-Depth Segmentation

Optimal sequencing depth for measuring the concentrations of molecular barcodes

Quantifying molecular bias in DNA data storage

A novel method addressing NGS-based mappability bias for sensitive detection of DNA alterations

Digital Rna Sequencing Minimizes Sequence-Dependent Bias And Amplification Noise With Optimized Single-Molecule Barcodes

Barcode demultiplexing of nanopore sequencing raw signals by unsupervised machine learning

Analysis and correction of compositional bias in sparse sequencing count data

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

Mixed Sequence Reader: A Program for Analyzing DNA Sequences with Heterozygous Base Calling

An Online Copy Number Variant Detection Method for Short Sequencing Reads

Fundamental Bounds and Approaches to Sequence Reconstruction from Nanopore Sequencers