A total crapshoot? Evaluating bioinformatic decisions in animal diet metabarcoding analyses

Citations Over TimeTop 10% of 2020 papers

Abstract

Metabarcoding studies provide a powerful approach to estimate the diversity and abundance of organisms in mixed communities in nature. While strategies exist for optimizing sample and sequence library preparation, best practices for bioinformatic processing of amplicon sequence data are lacking in animal diet studies. Here we evaluate how decisions made in core bioinformatic processes, including sequence filtering, database design, and classification, can influence animal metabarcoding results. We show that denoising methods have lower error rates compared to traditional clustering methods, although these differences are largely mitigated by removing low-abundance sequence variants. We also found that available reference datasets from GenBank and BOLD for the animal marker gene cytochrome oxidase I (COI) can be complementary, and we discuss methods to improve existing databases to include versioned releases. Taxonomic classification methods can dramatically affect results. For example, the commonly used Barcode of Life Database (BOLD) Classification API assigned fewer names to samples from order through species levels using both a mock community and bat guano samples compared to all other classifiers (vsearch-SINTAX and q2-feature-classifier's BLAST + LCA, VSEARCH + LCA, and Naive Bayes classifiers). The lack of consensus on bioinformatics best practices limits comparisons among studies and may introduce biases. Our work suggests that biological mock communities offer a useful standard to evaluate the myriad computational decisions impacting animal metabarcoding accuracy. Further, these comparisons highlight the need for continual evaluations as new tools are adopted to ensure that the inferences drawn reflect meaningful biology instead of digital artifacts.

Related Papers

• → GenBank is a reliable resource for 21st century biodiversity research(2019)200 cited
• → Many Genbank entries for complete microbial genomes violate the Genbank standard(2001)9 cited
• → AnnotationBustR: An R package to extract subsequences from GenBank annotations(2017)1 cited
• → An automated workflow to assess completeness and curate GenBank for eDNA metabarcoding: the marine fish assemblage as case study(2022)3 cited
• → Exploration of genetic data from GenBank using web services(2011)