Database size positively correlates with the loss of species-level taxonomic resolution for the 16S rRNA and other prokaryotic marker genes

Abstract

For decades, the 16S rRNA gene has been used to taxonomically classify prokaryotic species and to taxonomically profile microbial communities. The 16S rRNA gene has been criticized for being too conserved to differentiate between distinct species. We argue that the inability to differentiate between species is not a unique feature of the 16S rRNA gene. Rather, we observe the gradual loss of species-level resolution for other marker genes as the number of gene sequences increases in reference databases. We demonstrate this effect through the analysis of three commonly used databases of nearly-universal prokaryotic marker genes: the SILVA 16S rRNA gene database, the Genome Taxonomy Database (GTDB), and a set of 40 taxonomically-informative single-copy genes. Our results reflect a more fundamental property of the taxonomies themselves and have broad implications for bioinformatic analyses beyond taxonomic classification. Effective solutions for fine-level taxonomic classification require a more precise, and operationally-relevant, definition of the taxonomic labels being sought, and the use of combinations of genomic markers in the classification process. ImportanceThe use of reference databases for assigning taxonomic labels to genomic and metagenomic sequences is a fundamental bioinformatic task in the characterization of microbial communities. The increasing accessibility of high throughput sequencing has led to a rapid increase in the size and number of sequences in databases. This has been beneficial for improving our understanding of the global microbial genetic diversity. However, there is evidence that as the microbial diversity is more densely sampled, increasingly longer genomic segments are needed to differentiate between distinct species. The scientific community needs to be aware of this issue and needs to develop methods that better account for it when assigning taxonomic labels to metagenomic sequences from microbial communities.