Should names reflect the evolution of bacterial species?

A Genomic redefinition of Pseudomonas avellanae species

The circumscription of bacterial species is a complex task. So far, DNA-DNA hybridization (DDH), 16S rRNA gene sequencing, and multiocus sequence typing analysis (MLSA) are currently the preferred techniques for their genetic determination. However, the average nucleotide identity (ANI) analysis of conserved and shared genes between two bacterial strains based on the pair-wise genome comparisons, with support of the tetranucleotide frequency correlation coefficients (TETRA) value, has recently been proposed as a reliable substitute for DDH. The species demarcation boundary has been set to a value of 95-96% of the ANI identity, with further confirmation through the assessment of the corresponding TETRA value. In this study, we performed a genome-wide MLSA of 14 phytopathogenic pseudomonads genomes, and assessed the ANI and TETRA values of 27 genomes, representing seven out of the nine genomospecies of Pseudomonas spp. sensu Gardan et alii, and their phylogenetic relationships using maximum likelihood and Bayesian approaches. The results demonstrate the existence of a well demarcated genomic cluster that includes strains classified as P. avellanae, P. syringae pv. theae, P. s. pv. actinidiae and one P. s. pv. morsprunorum strain all belonging to the single species P. avellanae. In addition, when compared with P. avellanae, five strains of P. s. pv. tomato, including the model strain DC3000, and one P. s. pv. lachrymans strain, appear as very closely related to P. avellanae, with ANI values of nearly 96% as confirmed by the TETRA analysis. Conversely, one representative strain, previously classified as P. avellanae and isolated in central Italy, is a genuine member of the P. syringae species complex and can be defined as P. s. pv. avellanae. Currently. The core and pan genomes of P. avellanae species consist of 3,995 and 5,410 putative protein-coding genes, respectively.

Postgenomic taxonomy of human ureaplasmas -- a case study based on multiple gene sequences

In 2000, the full genome sequence of Ureaplasma parvum (previously known as Ureaplasma urealyticum) serovar 3 was released. In 2002, after prolonged debate, it was agreed that the former U. urealyticum should be divided into two species -- U. parvum and U. urealyticum. To provide additional support for this decision and improve our understanding of the relationship between these two species, the authors studied four 'core' genes or gene clusters in ATCC reference strains of all 14 serovars of U. parvum and U. urealyticum. These 'core' regions were the rRNA gene clusters, the EF-Tu genes (tuf), urease gene clusters and multiple-banded antigen genes (mba). The known U. parvum genome sequences (GenBank accession no. NC_002162) were used as reference. DNA insertions and deletions (indels) were found in all of the gene regions studied, except tuf, but they were found only between, not within, the two species. An incidental finding was that there was inter-copy heterogeneity for rRNA gene cluster sequences. Sequence analysis (sequence heterogeneity and especially indels) of all four selected targets consistently supported the separation of human ureaplasmas into two species. Except for multiple-banded antigen, there was less heterogeneity in amino acid sequences of proteins, between species, than in the nucleic acid sequences of the corresponding genes. The degrees of heterogeneity at the 5' end of the species-specific regions of multiple-banded antigen were almost identical for both amino acid and nucleotide sequences. Analysis of the authors' results provided an interesting case study to help resolve some common problems in the use of sequence data to infer phylogenetic relationships and support taxonomic changes. It is recommended that, to avoid confusion, the new nomenclature be used for human ureaplasmas in future publications.