From defaults to databases: parameter and database choice dramatically impact the performance of metagenomic taxonomic classification tools
Top Cited Papers
Open Access
- 3 March 2023
- journal article
- research article
- Published by Microbiology Society in Microbial Genomics
- Vol. 9 (3) , 000949
Abstract
In metagenomic analyses of microbiomes, one of the first steps is usually the taxonomic classification of reads by comparison to a database of previously taxonomically classified genomes. While different studies comparing metagenomic taxonomic classification methods have determined that different tools are ‘best’, there are two tools that have been used the most to-date: Kraken (k-mer-based classification against a user-constructed database) and MetaPhlAn (classification by alignment to clade-specific marker genes), the latest versions of which are Kraken2 and MetaPhlAn 3, respectively. We found large discrepancies in both the proportion of reads that were classified as well as the number of species that were identified when we used both Kraken2 and MetaPhlAn 3 to classify reads within metagenomes from human-associated or environmental datasets. We then investigated which of these tools would give classifications closest to the real composition of metagenomic samples using a range of simulated and mock samples and examined the combined impact of tool–parameter–database choice on the taxonomic classifications given. This revealed that there may not be a one-size-fits-all ‘best’ choice. While Kraken2 can achieve better overall performance, with higher precision, recall and F1 scores, as well as alpha- and beta-diversity measures closer to the known composition than MetaPhlAn 3, the computational resources required for this may be prohibitive for many researchers, and the default database and parameters should not be used. We therefore conclude that the best tool–parameter–database choice for a particular application depends on the scientific question of interest, which performance metric is most important for this question and the limit of available computational resources.Keywords
All Related Versions
Funding Information
- Natural Sciences and Engineering Research Council of Canada (2016-05039)
- Dalhousie Medical Research Foundation
This publication has 59 references indexed in Scilit:
- The SILVA ribosomal RNA gene database project: improved data processing and web-based toolsNucleic Acids Research, 2012
- A framework for human microbiome researchNature, 2012
- Metagenomic microbial community profiling using unique clade-specific marker genesNature Methods, 2012
- Fast gapped-read alignment with Bowtie 2Nature Methods, 2012
- Taxonomic classification of metagenomic shotgun sequences with CARMA3Nucleic Acids Research, 2011
- RAIphy: Phylogenetic classification of metagenomics samples using iterative refinement of relative abundance index profilesBMC Bioinformatics, 2011
- Use of simulated data sets to evaluate the fidelity of metagenomic processing methodsNature Methods, 2007
- UniFrac: a New Phylogenetic Method for Comparing Microbial CommunitiesApplied and Environmental Microbiology, 2005
- Conservation evaluation and phylogenetic diversityBiological Conservation, 1992
- An Ordination of the Upland Forest Communities of Southern WisconsinEcological Monographs, 1957