Dataset Information

AmpliCI: A High-resolution Model-Based Approach for Denoising Illumina Amplicon Data.

ABSTRACT:

Motivation

Next-generation amplicon sequencing is a powerful tool for investigating microbial communities. A main challenge is to distinguish true biological variants from errors caused by amplification and sequencing. In traditional analyses, such errors are eliminated by clustering reads within a sequence similarity threshold, usually 97%, and constructing operational taxonomic units, but the arbitrary threshold leads to low resolution and high false positive rates. Recently developed "denoising" methods have proven able to resolve single-nucleotide amplicon variants, but they still miss low frequency sequences, especially those near more frequent sequences, because they ignore the sequencing quality information.

Results

We introduce AmpliCI, a reference-free, model-based method for rapidly resolving the number, abundance and identity of error-free sequences in massive Illumina amplicon datasets. AmpliCI takes into account quality information and allows the data, not an arbitrary threshold or an external database, to drive conclusions. AmpliCI estimates a finite mixture model, using a greedy strategy to gradually select error-free sequences and approximately maximize the likelihood. AmpliCI has better performance than three popular denoising methods, with acceptable computation time and memory usage.

Availability

Source code is available at https://github.com/DormanLab/AmpliCI.

Supplementary information

Supplementary material are available at Bioinformatics online.

SUBMITTER: Peng X

PROVIDER: S-EPMC7850112 | biostudies-literature | 2020 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

AmpliCI: a high-resolution model-based approach for denoising Illumina amplicon data.

Peng Xiyu X Dorman Karin S KS

Bioinformatics (Oxford, England) 20210101 21

<h4>Motivation</h4>Next-generation amplicon sequencing is a powerful tool for investigating microbial communities. A main challenge is to distinguish true biological variants from errors caused by amplification and sequencing. In traditional analyses, such errors are eliminated by clustering reads within a sequence similarity threshold, usually 97%, and constructing operational taxonomic units, but the arbitrary threshold leads to low resolution and high false-positive rates. Recently developed ...[more]

PMID: 32697845

Similar Datasets

Project description:BACKGROUND:Sequencing of the 16S rRNA gene has been the standard for studying the composition of microbial communities. While it allows identification of bacteria at the level of species, this method does not usually provide sufficient information to resolve communities at the sub-species level. Species-level resolution is not adequate for studies of transmission or stability or for exploring subspecies variation in disease association. Strain level analysis using whole metagenome shotgun sequencing has significant limitations that can make it unsuitable for large-scale studies. Achieving sufficient depth of sequencing can be cost-prohibitive, and even with adequate coverage, deconvoluting complex communities such as the oral microbiota is computationally very challenging. Thus, there is a need for high-resolution, yet cost-effective, high-throughput methods for characterizing microbial communities. RESULTS:Significant improvement in resolution for amplicon-based bacterial community analysis was achieved by combining amplicon sequencing of a high-diversity marker gene, the ribosomal 16-23S intergenic spacer region (ISR), with a probabilistic error modeling based denoising algorithm, DADA2. The resolving power of this new approach was compared to that of both standard and high-resolution 16S-based approaches using a set of longitudinal subgingival plaque samples. The ISR strategy resulted in a 5.2-fold increase in community resolution compared to reference-based 16S rRNA gene analysis and showed 100% accuracy in predicting the correct source of a clinical sample. Individuals' microbial communities were highly personalized, and although they exhibited some drift in membership and levels over time, that difference was always smaller than the differences between any two subjects, even after 1 year. The construction of an ISR database from publicly available genomic sequences allowed us to explore genomic variation within species, resulting in the identification of multiple variants of the ISR for most species. CONCLUSIONS:The ISR approach resulted in significantly improved resolution of communities and revealed a highly personalized human oral microbiota that was stable over 1 year. Multiple ISR types were observed for all species examined, demonstrating a high level of subspecies variation in the oral microbiota. The approach is high-throughput, high-resolution yet cost-effective, allowing subspecies-level community fingerprinting at a cost comparable to that of 16S rRNA gene amplicon sequencing. It will be useful for a range of applications that require high-resolution identification of organisms, including microbial tracking, community fingerprinting, and potentially for identification of virulence-associated strains.