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.

Metagenomic coverage bias at transcription start sites is correlated with gene expression

Codon usages in different gene classes of the Escherichia coli genome

Multiple Factors Drive Replicating Strand Composition Bias in Bacterial Genomes

Codon bias as a factor in regulating expression via translation rate in the human genome

Relationships among Stop Codon Usage Bias, Its Context, Isochores, and Gene Expression Level in Various Eukaryotes

Microbial Diversity Biased Estimation Caused by Intragenomic Heterogeneity and Interspecific Conservation of 16S Rrna Genes.

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

Identifying biases and their potential solutions in human microbiome studies

Systematic dissection of biases in whole-exome and whole-genome sequencing reveals major determinants of coding sequence coverage

The impact of RNA secondary structure on read start locations on the Illumina sequencing platform

Genomic GC-Content Affects the Accuracy of 16S rRNA Gene Sequencing Based Microbial Profiling due to PCR Bias

Codon usage of highly expressed genes affects proteome-wide translation efficiency

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression

Intragenic spatial patterns of codon usage bias in prokaryotic and eukaryotic genomes

Gene length and detection bias in single cell RNA sequencing protocols

Transcriptome-wide meta-analysis of codon usage in Escherichia coli

Selection on codon bias in yeast: a transcriptional hypothesis

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

Growth-dependent gene expression variation influences the strength of codon usage biases

Gene representation bias in spatial transcriptomics