Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin

Abstract

BackgroundTaxonomic classification of marker-gene sequences is an important step in microbiome analysis.ResultsWe present q2-feature-classifier (https://github.com/qiime2/q2-feature-classifier), a QIIME 2 plugin containing several novel machine-learning and alignment-based methods for taxonomy classification. We evaluated and optimized several commonly used classification methods implemented in QIIME 1 (RDP, BLAST, UCLUST, and SortMeRNA) and several new methods implemented in QIIME 2 (a scikit-learn naive Bayes machine-learning classifier, and alignment-based taxonomy consensus methods based on VSEARCH, and BLAST+) for classification of bacterial 16S rRNA and fungal ITS marker-gene amplicon sequence data. The naive-Bayes, BLAST+-based, and VSEARCH-based classifiers implemented in QIIME 2 meet or exceed the species-level accuracy of other commonly used methods designed for classification of marker gene sequences that were evaluated in this work. These evaluations, based on 19 mock communities and error-free sequence simulations, including classification of simulated “novel” marker-gene sequences, are available in our extensible benchmarking framework, tax-credit (https://github.com/caporaso-lab/tax-credit-data).ConclusionsOur results illustrate the importance of parameter tuning for optimizing classifier performance, and we make recommendations regarding parameter choices for these classifiers under a range of standard operating conditions. q2-feature-classifier and tax-credit are both free, open-source, BSD-licensed packages available on GitHub.

Anacapa Toolkit: An environmental DNA toolkit for processing multilocus metabarcode datasets

“…However, the lack of accurate and easily customizable bioinformatic pipelines is limiting the broader application of eDNA approaches. The Anacapa Toolkit provides enhanced functionality for eDNA projects and can be used for other common applications such as gut content analysis (Leray et al, ), autonomous reef monitoring structures (Ransome et al, ), and microbiomes (Bokulich et al, ). Importantly, the Anacapa Toolkit is modular and its parameters are easily modifiable, making it easily adapted to user specific needs in several important ways.…”

“…Importantly, the Anacapa Toolkit is modular and its parameters are easily modifiable, making it easily adapted to user specific needs in several important ways. First, CRUX reference databases are compatible with alternative classifiers (Bokulich et al, ), and users can append their own reference sequences to CRUX databases as needed. Second, the Q uality Control and ASV Parsing module is designed to process pooled metabarcoding libraries and automatically sort them by barcode and sample.…”

Gut microbiota dysbiosis in a cohort of patients with psoriasis

“…The resulting sequences were stored in a fastq file generated for each sample and processed using QIIME2 v2018·11 software, importing data through a manifest file. Then, single‐end raw reads were filtered, quality retained for quality score > 20 and chimeras were removed using DADA2 v2018·4 .…”

Rhizosphere modelling reveals spatiotemporal distribution of daidzein shaping soybean rhizosphere bacterial community

“…The α‐ and β‐diversity metrics were estimated from a subsampled ASV dataset, with 45,000 sequences per sample (Table S2). Each ASV was identified using a Naïve Bayes classifier trained on 16S rRNA gene sequences from the Greengenes release 13_8 dataset (Bokulich et al, ; Table S3). Principal coordinate analysis and other statistical analyses were performed using custom R and Python scripts.…”