GenBank mining reveals novel insights into Rhizobium phylogeny: Identical 16S rRNA sequences are mainly uncoupled from species designation, host plant, and geographic origin: How this search suggested the definition of a direct ‘microbial h-index’
- 1Institute for the Animal Production System in Mediterranean Environment-National Research Council (ISPAAM CNR), Traversa La Crucca 3, Sassari, Italy
- 2Department of Agricultural Sciences, University of Sassari, Viale Italia, Sassari, Italy
- 3Interdepartmental Center for the Conservation and Valorization of Plant Biodiversity, University of Sassari, Alghero, Sassari, Italy
Abstract
16S rDNA is the historical gold standard for bacterial identification, particularly in metabarcoding approaches reliant on sequence similarity thresholds. We analyzed 6,660 Rhizobium 16S rRNA gene sequences from GenBank to examine the relationship between sequence identity and three metadata: species name, host plant, and geographic origin. Using an iterative BLAST-based pipeline, we detected 116,069 pairwise matches and assessed concordance among sequences (average length 1,328 bp) sharing 100% identity. For those in which the organism name, host plant and country of isolation were present in the record, surprisingly, 66.59% of identical sequence pairs showed full discordance across all three metadata, while only 1.40% shared the same name, host, and country. The most widespread sequence, detected 371 times, was associated with over 56 different host plants across 25 countries and bore multiple species name designations. These results highlight a striking mismatch between the 16S barcode and the taxonomic, ecological, and phenotypic variability it is assumed to reflect, likely arising from the slow evolution of rRNA genes contrasted with the mobility of ecologically relevant genes via horizontal transfer on plasmids, transposons, and phages. Our findings further challenge the limitations of relying on 16S rRNA alone for fine-scale taxonomic and metadata-based inference in capturing the true functional and ecological diversity of bacteria, endorsing the critical importance of polyphasic taxonomic approaches that integrate genomic, phenotypic, and ecological data. An interesting byproduct of the analysis was to realize the possibility of treating these data as if they were ‘citations.’ The more one finds the same query sequence, the more that sequence can be considered biologically ‘cited’, i.e., re-proposed elsewhere in the world. Thus, one can also analyze the h-index of such a ranking. In our Rhizobium dataset, we calculated an h-index = 201, meaning the sequence ranked 201st had 202 identical homologues in GenBank. Although the research effort on given species is directly connected with it, this number provides a quantitative indicator of a taxon’s sequence recurrence and distribution within public databases, independent of nomenclatural inconsistencies, offering a novel framework for assessing bacterial representation across global datasets.
Related
- Cyanogenic clovers invest more in rhizobia but do not affect rhizobium evolution
- A broad-host-range Rhizobium rhizogenes strain enables transient expression across diverse crop species
- How agroecological transition can jeopardize compliance with Protected Designation of Origin specifications: A living lab case study in Fourme de Montbrison cheese
- Field-Zone-Associated Transcriptomic and Metabolic Divergence Despite Nearly Identical Observed Chlorophyll in an Agrivoltaic Tea Setting
- Phylogeny of Termites Inferred from Mitogenomic Data with a Large Taxon Sampling
BibTeX
@article{muresu2026genbank,
title = {{GenBank mining reveals novel insights into Rhizobium phylogeny: Identical 16S rRNA sequences are mainly uncoupled from species designation, host plant, and geographic origin: How this search suggested the definition of a direct ‘microbial h-index’}},
author = {Rosella Muresu and Monica Rodriguez and Andrea Squartini},
journal = {PLOS One},
year = {2026},
month = sep,
doi = {10.1371/journal.pone.0357973},
url = {https://doi.org/10.1371/journal.pone.0357973},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - GenBank mining reveals novel insights into Rhizobium phylogeny: Identical 16S rRNA sequences are mainly uncoupled from species designation, host plant, and geographic origin: How this search suggested the definition of a direct ‘microbial h-index’
AU - Rosella Muresu
AU - Monica Rodriguez
AU - Andrea Squartini
PY - 2026
DA - 2026/09/11
JO - PLOS One
DO - 10.1371/journal.pone.0357973
UR - https://doi.org/10.1371/journal.pone.0357973
AB - 16S rDNA is the historical gold standard for bacterial identification, particularly in metabarcoding approaches reliant on sequence similarity thresholds. We analyzed 6,660 Rhizobium 16S rRNA gene sequences from GenBank to examine the relationship between sequence identity and three metadata: species name, host plant, and geographic origin. Using an iterative BLAST-based pipeline, we detected 116,069 pairwise matches and assessed concordance among sequences (average length 1,328 bp) sharing 100% identity. For those in which the organism name, host plant and country of isolation were present in the record, surprisingly, 66.59% of identical sequence pairs showed full discordance across all three metadata, while only 1.40% shared the same name, host, and country. The most widespread sequence, detected 371 times, was associated with over 56 different host plants across 25 countries and bore multiple species name designations. These results highlight a striking mismatch between the 16S barcode and the taxonomic, ecological, and phenotypic variability it is assumed to reflect, likely arising from the slow evolution of rRNA genes contrasted with the mobility of ecologically relevant genes via horizontal transfer on plasmids, transposons, and phages. Our findings further challenge the limitations of relying on 16S rRNA alone for fine-scale taxonomic and metadata-based inference in capturing the true functional and ecological diversity of bacteria, endorsing the critical importance of polyphasic taxonomic approaches that integrate genomic, phenotypic, and ecological data. An interesting byproduct of the analysis was to realize the possibility of treating these data as if they were ‘citations.’ The more one finds the same query sequence, the more that sequence can be considered biologically ‘cited’, i.e., re-proposed elsewhere in the world. Thus, one can also analyze the h-index of such a ranking. In our Rhizobium dataset, we calculated an h-index = 201, meaning the sequence ranked 201st had 202 identical homologues in GenBank. Although the research effort on given species is directly connected with it, this number provides a quantitative indicator of a taxon’s sequence recurrence and distribution within public databases, independent of nomenclatural inconsistencies, offering a novel framework for assessing bacterial representation across global datasets.
SN - 1932-6203
ER -