Genome Taxonomy Database

The Genome Taxonomy Database (GTDB) is an online database that maintains information on a proposed nomenclature of prokaryotes, following a phylogenomic approach based on a set of conserved single-copy proteins. In addition to resolving paraphyletic groups, this method also reassigns taxonomic ranks algorithmically, updating names in both cases. ^[1] Information for archaea was added in 2020, ^[2] along with a species classification based on average nucleotide identity. ^[3] Each update incorporates new genomes as well as automated and manual curation of the taxonomy. ^[4]

Contents

• Methodology
• Contents of the database
• Effects on the accepted taxonomy
• See also
• References
• Further reading
• External links

An open-source tool called GTDB-Tk is available to classify draft genomes into the GTDB hierarchy. ^[5] The GTDB system, via GTDB-Tk, has been used to catalogue not-yet-named bacteria in the human gut microbiome and other metagenomic sources. ^{[6] [7]}

The GTDB is incorporated into the Bergey's Manual of Systematics of Archaea and Bacteria in 2019 as its phylogenomic resource. ^[8]

Methodology

The genomes used to construct the phylogeny are obtained from NCBI (RefSeq and Genbank), and GTDB releases are indexed to RefSeq releases, starting with release 76. Importantly and increasingly, this dataset includes draft genomes of uncultured microorganisms obtained from metagenomes and single cells, ensuring improved genomic representation of the microbial world. All genomes are independently quality controlled using CheckM before inclusion in GTDB. ^[9]

Genomes first undergo gene calling to extract genes. The taxonomy is based on trees inferred with FastTree from an aligned concatenated set of 120 single copy marker proteins for Bacteria under a WAG model, and with IQ-TREE from a concatenated set of 53 (since RS207; 122 before) marker proteins for Archaea under the PMSF model. Additional marker sets are also used to cross-validate tree topologies including concatenated ribosomal proteins and ribosomal RNA genes. ^[9] The relative evolutionary divergence (RED) metric, which determines the taxonomic ranks used, is derived from the two main trees by the PhyloRank program. ^[1]

Species are deliminated using average nucleotide identity and alignment fraction, both calculated by skani. For species existing in a previous release, GTDB compares the quality and position of two genomes and may decide to switch to a new species representative genome. ^[9]

Taxomony comes from the following sources:

• A previous release, if available for the neighborhood of genomic similarity.
• National Center for Biotechnology Information (NCBI) taxonomy was initially used to decorate the genome tree via tax2tree. ^[1]
• The 16S rRNA-based Greengenes taxonomy is used to supplement the taxonomy particularly in regions of the tree with no cultured representatives. ^[1]
• List of Prokaryotic names with Standing in Nomenclature (LPSN) is used as the primary taxonomic authority for establishing naming priorities. ^[1]

GTDB personnel curates the taxonomy from the aforementioned sources by checking them against the results of PhyloRank and the tree.

• The tree node corresponding to a taxon name may have a RED inappropriate for its rank. The name may either be moved onto another node or (by changing the Latin suffix) into a different rank. ^[1]
  • Splitting may happen on the level of species or genera if the divergence turns out too high. Doing so creates new taxa. ^[3]
• The taxon may turn out to be polyphyletic. The curator first restricts the taxon to the clade containing its type material. A new taxon is created for each of the other clades. ^[1]

For the each new taxon, the curators try to find a proposed name in literature for it. If there is no name proposed, the taxon is given a placeholder name by adding a suffix to the original name, e.g. Lactobacillus gasseri_A. After "Z" comes "AA". ^[1]

Contents of the database

Each release contains: ^[10]

• Taxonomy tables containing the assignment of all included genome assemblies to the phylum-to-species taxonomy. (One per domain.)
• Files containing the metadata given to each genome assembly, including original taxonomy from NCBI, original strain identifier, GTDB taxonomy, quality estimates, and presence of important genes (tRNA and rRNA). (One per domain.)
• Species tree Newick files containing the species-representative genomes (1 per species), built as described in the previous section. (One per domain.)
• For species-representative genomes:
  • alignments of marker genes identified from these genomes
  • file containing one 16S rRNA sequence from each species
  • tarballs containing amino acid and nucleotide versions of all predicted genes in these genomes
  • tarball containing the full contents of all these genomes
• For all genomes that pass quality check:
  • alignments of marker genes identified from these genomes
  • file containing all 16S rRNA sequences identified from these genomes
• Auxiliary files; see the full FILE_DESCRIPTIONS.txt。

The web interface displays a tree based on the taxonomy (not the entire Newick file), down to the genome assembly level. Each genome assembly has a page detailing its metadata and a history of how it's classified in each GTDB release. There is a search functionality.

Effects on the accepted taxonomy

GTDB "has now become an important resource for prokaryotic taxonomy". Both its species tree and elements of its methodology are used by taxonomists to improve the current, accepted taxonomy under the Prokaryotic Code. For example, a taxonomist may make references to the GTDB tree on top of their own phylogenetic tree to further support a taxonomic proposal. ^[11]

There has been even more ambitious proposals to import large parts of the database into the accepted taxonomy. A 2022 article in the IJSEM, written by third-party authors, proposes to assign names based on meaningless Latin syllable to over 65 thousand GTDB taxa, ^[12] though none of these names have made their way into the LPSN. A 2023 article by the GTDB team proposes to import 223 higher-order taxa into the Prokaryotic Code system and 49 under the SeqCode system. ^[13] Many of the names published under the under the Prokaryotic Code have already been validated. ^[14] (The SeqCode requires registration of the names for valid publication, which has also been done.)

See also

• PhyloCode
• National Center for Biotechnology Information
• SILVA ribosomal RNA database
• List of Prokaryotic names with Standing in Nomenclature

Related Research Articles

The Geobacteraceae are a family within the Thermodesulfobacteriota.

Desulfovibrionales are a taxonomic order of bacteria belonging to the phylum Thermodesulfobacteriota, with four families. They are Gram-negative. The majority are sulfate-reducing, with the exception of Lawsonia and Bilophila. All members of this order are obligately anaerobic. Most species are mesophilic, but some are moderate thermophiles.

The Syntrophobacterales are an order of Thermodesulfobacteriota. All genera are strictly anaerobic. Many of the family Syntrophobacteraceae are sulfate-reducing. Some species are motile by using one polar flagellum.

The Syntrophobacteraceae are a family of Thermodesulfobacteriota.

Desulfovibrionaceae is a family of bacteria belonging to the phylum Thermodesulfobacteriota.

Desulfohalobiaceae is a family of bacteria belonging to the phylum Thermodesulfobacteriota.

The Myxococcota are a phylum of bacteria known as the fruiting gliding bacteria. All species of this group are Gram-negative. They are predominantly aerobic genera that release myxospores in unfavorable environments.

The "Aigarchaeota" are a proposed archaeal phylum of which the main representative is Caldiarchaeum subterraneum. It is not yet clear if this represents a new phylum or a Nitrososphaerota order, since the genome of Caldiarchaeum subterraneum encodes several Nitrososphaerota-like features. The name "Aigarchaeota" comes from the Greek αυγή, avgí, meaning "dawn" or "aurora", for the intermediate features of hyperthermophilic and mesophilic life during the evolution of its lineage.

Parvarchaeota is a phylum of archaea belonging to the DPANN archaea. They have been discovered in acid mine drainage waters and later in marine sediments. The cells of these organisms are extremely small consistent with small genomes. Metagenomic techniques allow obtaining genomic sequences from non-cultured organisms, which were applied to determine this phylum.

The candidate division SR1 and gracilibacteria code is used in two groups of uncultivated bacteria found in marine and fresh-water environments and in the intestines and oral cavities of mammals among others. The difference to the standard and the bacterial code is that UGA represents an additional glycine codon and does not code for termination. A survey of many genomes with the codon assignment software Codetta, analyzed through the GTDB taxonomy system shows that this genetic code is limited to the Patescibacteria order BD1-5, not what are now termed Gracilibacteria, and that the SR1 genome assembly GCA_000350285.1 for which the table 25 code was originally defined is actually using the Absconditibacterales genetic code and has the associated three special recoding tRNAs. Thus this code may now be better named the "BD1-5 code".

Syntrophus is a Gram negative bacterial genus from the family of Syntrophaceae.

The candidate phyla radiation is a large evolutionary radiation of bacterial lineages whose members are mostly uncultivated and only known from metagenomics and single cell sequencing. They have been described as nanobacteria or ultra-small bacteria due to their reduced size (nanometric) compared to other bacteria.

NC10 is a bacterial phylum with candidate status, meaning its members remain uncultured to date. The difficulty in producing lab cultures may be linked to low growth rates and other limiting growth factors.

There are several models of the branching order of bacterial phyla, one of these is the Genome Taxonomy Database (GTDB).

Bdellovibrionota is a phylum of bacteria.

The Enterosoma genetic code translates AGG to methionine, as determined by the codon assignment software Codetta; it was further shown that this recoding is associated with a special tRNA with the appropriate anticodon and tRNA identity elements. The code is found in a small clade of species within the Enterosoma genus, according to the GTDB taxonomy system release 220. Codetta called the Enterosoma code for the following genome assemblies: GCA_002431755.1, GCA_002439645.1, GCA_002436825.1, GCA_002451385.1, GCA_002297105.1, GCA_002297045.1, GCA_002404995.1, and GCA_900549915.1.

The Anaerococcus and Onthovivens genetic code translates CGG to tryptophan, as determined by the codon assignment software Codetta; it was further shown that this recoding is associated with a special tRNA with the appropriate anticodon and tRNA identity elements appropriate for such decoding. As currently known, this code is limited to two distinct clades, the genus Anaerococcus in the class Clostridia and the genus Onthovivens in the class Bacilli, as defined by the GTDB taxonomy system release 220. Codetta originally called the Anaerococcus and Onthovivens code for the following genome assemblies: GCA_000024105.1, GCA_900445285.1, GCA_902500265.1, GCA_900258475.1, GCA_002399785.1, GCA_004558005.1, GCA_900540365.1, GCA_900540395.1, GCA_900545015.1.

The Absconditabacterales genetic code translates UGA to glycine, and CGG and GCA to tryptophan, as determined by the codon assignment software Codetta; it was further shown that these recodings are associated with three special tRNAs with appropriate anticodons and tRNA identity elements. Codetta called the Absconditibacterales code for the following genome assemblies: GCA_002792495.1, GCA_001007975.1, GCA_003488625.1, GCA_003260355.1, GCA_003242865.1, GCA_000350285.1, GCA_002746475.1, GCA_007116275.1, GCA_007115995.1, GCA_002361595.1, GCA_000503875.1, GCA_003543185.1, GCA_002441085.1, and GCA_002791215.1. Review of the GTDB taxonomy system for the order Absconditabacterales left two questionable genome assemblies ; spot-checking these two genomes shows that they both have all three special tRNAs, suggesting that the code is universal across the order.

References

Further reading

• Zhu, Q; Mai, U; Pfeiffer, W; et al. (2 December 2019). "Phylogenomics of 10,575 genomes reveals evolutionary proximity between domains Bacteria and Archaea". Nature Communications. 10 (1): 5477. Bibcode:2019NatCo..10.5477Z. doi: 10.1038/s41467-019-13443-4 . PMC   6889312 . PMID   31792218.– a similar independent effort that also evaluates GTDB taxonomic quality.

External links

• AnnoTree - third-party tool for visualization of genome annotations using the GTDB (R95 or R214) species tree.