DNA-DNA hybridization (DDH) has been used for nearly 50 years as the gold standard for prokaryotic species circumscriptions at the genomic level. It has been the only taxonomic method that offered a numerical and relatively stable species boundary, and its use has had a paramount influence on how the current classification has been constructed. However, now, in the era of genomics, DDH appears to be an outdated method for classification that needs to be substituted. The average nucleotide identity (ANI) between two genomes seems the most promising method since it mirrors DDH closely. Here we examine the work package JSpecies as a user-friendly, biologist-oriented interface to calculate ANI and the correlation of the tetranucleotide signatures between pairwise genomic comparisons. The results agreed with the use of ANI to substitute DDH, with a narrowed boundary that could be set at ≈95–96%. In addition, the JSpecies package implemented the tetranucleotide signature correlation index, an alignment-free parameter that generally correlates with ANI and that can be of help in deciding when a given pair of organisms should be classified in the same species. Moreover, for taxonomic purposes, the analyses can be produced by simply randomly sequencing at least 20% of the genome of the query strains rather than obtaining their full sequence.

https://www.pnas.org/content/pnas/106/45/19126.full.pdf

Shifting the genomic gold standard for the prokaryotic species definition

For the last 25 years species delimitation in prokaryotes (Archaea and Bacteria) was to a large extent based on DNA-DNA hybridization (DDH), a tedious lab procedure designed in the early 1970s that served its purpose astonishingly well in the absence of deciphered genome sequences. With the rapid progress in genome sequencing time has come to directly use the now available and easy to generate genome sequences for delimitation of species. GBDP (Genome Blast Distance Phylogeny) infers genome-to-genome distances between pairs of entirely or partially sequenced genomes, a digital, highly reliable estimator for the relatedness of genomes. Its application as an in-silico replacement for DDH was recently introduced. The main challenge in the implementation of such an application is to produce digital DDH values that must mimic the wet-lab DDH values as close as possible to ensure consistency in the Prokaryotic species concept. Correlation and regression analyses were used to determine the best-performing methods and the most influential parameters. GBDP was further enriched with a set of new features such as confidence intervals for intergenomic distances obtained via resampling or via the statistical models for DDH prediction and an additional family of distance functions. As in previous analyses, GBDP obtained the highest agreement with wet-lab DDH among all tested methods, but improved models led to a further increase in the accuracy of DDH prediction. Confidence intervals yielded stable results when inferred from the statistical models, whereas those obtained via resampling showed marked differences between the underlying distance functions. Despite the high accuracy of GBDP-based DDH prediction, inferences from limited empirical data are always associated with a certain degree of uncertainty. It is thus crucial to enrich in-silico DDH replacements with confidence-interval estimation, enabling the user to statistically evaluate the outcomes. Such methodological advancements, easily accessible through the web service at http://ggdc.dsmz.de
 , are crucial steps towards a consistent and truly genome sequence-based classification of microorganisms.

/pdf/genome-sequence-based-species-delimitation-with-confidence-31cdx9rh52.pdf

Genome sequence-based species delimitation with confidence intervals and improved distance functions

DNA-DNA hybridization (DDH) values have been used by bacterial taxonomists since the 1960s to determine relatedness between strains and are still the most important criterion in the delineation of bacterial species. Since the extent of hybridization between a pair of strains is ultimately governed by their respective genomic sequences, we examined the quantitative relationship between DDH values and genome sequence-derived parameters, such as the average nucleotide identity (ANI) of common genes and the percentage of conserved DNA. A total of 124 DDH values were determined for 28 strains for which genome sequences were available. The strains belong to six important and diverse groups of bacteria for which the intra-group 16S rRNA gene sequence identity was greater than 94 %. The results revealed a close relationship between DDH values and ANI and between DNA-DNA hybridization and the percentage of conserved DNA for each pair of strains. The recommended cut-off point of 70 % DDH for species delineation corresponded to 95 % ANI and 69 % conserved DNA. When the analysis was restricted to the protein-coding portion of the genome, 70 % DDH corresponded to 85 % conserved genes for a pair of strains. These results reveal extensive gene diversity within the current concept of "species". Examination of reciprocal values indicated that the level of experimental error associated with the DDH method is too high to reveal the subtle differences in genome size among the strains sampled. It is concluded that ANI can accurately replace DDH values for strains for which genome sequences are available.

DNA–DNA hybridization values and their relationship to whole-genome sequence similarities

Among available genome relatedness indices, average nucleotide identity (ANI) is one of the most robust measurements of genomic relatedness between strains, and has great potential in the taxonomy of bacteria and archaea as a substitute for the labour-intensive DNA–DNA hybridization (DDH) technique. An ANI threshold range (95–96 %) for species demarcation had previously been suggested based on comparative investigation between DDH and ANI values, albeit with rather limited datasets. Furthermore, its generality was not tested on all lineages of prokaryotes. Here, we investigated the overall distribution of ANI values generated by pairwise comparison of 6787 genomes of prokaryotes belonging to 22 phyla to see whether the suggested range can be applied to all species. There was an apparent distinction in the overall ANI distribution between intra- and interspecies relationships at around 95–96 % ANI. We went on to determine which level of 16S rRNA gene sequence similarity corresponds to the currently accepted ANI threshold for species demarcation using over one million comparisons. A twofold cross-validation statistical test revealed that 98.65 % 16S rRNA gene sequence similarity can be used as the threshold for differentiating two species, which is consistent with previous suggestions (98.2–99.0 %) derived from comparative studies between DDH and 16S rRNA gene sequence similarity. Our findings should be useful in accelerating the use of genomic sequence data in the taxonomy of bacteria and archaea.

Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes

A fundamental question in microbiology is whether there is continuum of genetic diversity among genomes, or clear species boundaries prevail instead. Whole-genome similarity metrics such as Average Nucleotide Identity (ANI) help address this question by facilitating high resolution taxonomic analysis of thousands of genomes from diverse phylogenetic lineages. To scale to available genomes and beyond, we present FastANI, a new method to estimate ANI using alignment-free approximate sequence mapping. FastANI is accurate for both finished and draft genomes, and is up to three orders of magnitude faster compared to alignment-based approaches. We leverage FastANI to compute pairwise ANI values among all prokaryotic genomes available in the NCBI database. Our results reveal clear genetic discontinuity, with 99.8% of the total 8 billion genome pairs analyzed conforming to >95% intra-species and <83% inter-species ANI values. This discontinuity is manifested with or without the most frequently sequenced species, and is robust to historic additions in the genome databases. Average Nucleotide Identity (ANI) is a robust and useful measure to gauge genetic relatedness between two genomes. Here, the authors develop FastANI, a method to compute ANI using alignment-free approximate sequence mapping, and show 95% ANI is an accurate threshold for demarcating prokaryotic species by analyzing about 90,000 prokaryotic genomes.

/pdf/high-throughput-ani-analysis-of-90k-prokaryotic-genomes-120jzpm1w0.pdf

High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries.

A strain identified as Comamonas testosteroni I2 was isolated from activated sludge and found to be able to mineralize 3-chloroaniline (3-CA). During the mineralization, a yellow intermediate accumulated temporarily, due to the distal meta-cleavage of chlorocatechol. This strain was tested for its ability to clean wastewater containing 3-CA upon inoculation into activated sludge. To monitor its survival, the strain was chromosomally marked with the gfp gene and designated I2gfp. After inoculation into a lab-scale semicontinuous activated-sludge (SCAS) system, the inoculated strain maintained itself in the sludge for at least 45 days and was present in the sludge flocs. After an initial adaptation period of 6 days, complete degradation of 3-CA was obtained during 2 weeks, while no degradation at all occurred in the noninoculated control reactor. Upon further operation of the SCAS system, only 50% 3-CA removal was observed. Denaturing gradient gel electrophoresis (DGGE) of 16S rRNA genes revealed a dynamic change in the microbial community structure of the activated sludge. The DGGE patterns of the noninoculated and the inoculated reactors evolved after 7 days to different clusters, which suggests an effect of strain inoculation on the microbial community structure. The results indicate that bioaugmentation, even with a strain originating from that ecosystem and able to effectively grow on a selective substrate, is not permanent and will probably require regular resupplementation.

Bioaugmentation of Activated Sludge by an Indigenous 3-Chloroaniline-Degrading Comamonas testosteroni Strain, I2gfp

Burkholderia tuberum sp. nov. and Burkholderia phymatum sp. nov., Nodulate the Roots of Tropical Legumes

Previously reported hybridization data, obtained with the initial renaturation method (De Ley et al. 1970), were compared with data obtained with a microplate DNA - DNA hybridization method (Ezaki ...

Evaluation of a microplate DNA : DNA hybridization method compared with the initial renaturation method

The bacterial community composition of a linuron-degrading enrichment culture and the role of the individual strains in linuron degradation have been determined by a combination of methods, such as denaturing gradient gel electrophoresis of the total 16S rRNA gene pool, isolation and identification of strains, and biodegradation assays. Three strains, Variovorax sp. strain WDL1, Delftia acidovorans WDL34, and Pseudomonas sp. strain WDL5, were isolated directly from the linuron-degrading culture. In addition, subculture of this enrichment culture on potential intermediates in the degradation pathway of linuron (i.e., N,O-dimethylhydroxylamine and 3-chloroaniline) resulted in the isolation of, respectively, Hyphomicrobium sulfonivorans WDL6 and Comamonas testosteroni WDL7. Of these five strains, only Variovorax sp. strain WDL1 was able to use linuron as the sole source of C, N, and energy. WDL1 first converted linuron to 3,4-dichloroaniline (3,4-DCA), which transiently accumulated in the medium but was subsequently degraded. To the best of our knowledge, this is the first report of a strain that degrades linuron further than the aromatic intermediates. Interestingly, the rate of linuron degradation by strain WDL1 was lower than that for the consortium, but was clearly increased when WDL1 was coinoculated with each of the other four strains. D. acidovorans WDL34 and C. testosteroni WDL7 were found to be responsible for degradation of the intermediate 3,4-DCA, and H. sulfonivorans WDL6 was the only strain able to degrade N,O-dimethylhydroxylamine. The role of Pseudomonas sp. strain WDL5 needs to be further elucidated. The degradation of linuron can thus be performed by a single isolate, Variovorax sp. strain WDL1, but is stimulated by a synergistic interaction with the other strains isolated from the same linuron-degrading culture.

https://biblio.ugent.be/publication/211356/file/869564.pdf

Johan Goris

Papers

DNA–DNA hybridization values and their relationship to whole-genome sequence similarities

Bioaugmentation of Activated Sludge by an Indigenous 3-Chloroaniline-Degrading Comamonas testosteroni Strain, I2gfp

Burkholderia tuberum sp. nov. and Burkholderia phymatum sp. nov., Nodulate the Roots of Tropical Legumes

Evaluation of a microplate DNA : DNA hybridization method compared with the initial renaturation method

Synergistic degradation of linuron by a bacterial consortium and isolation of a single linuron-degrading variovorax strain.