High-throughput ANI Analysis of 90K Prokaryotic Genomes Reveals Clear Species Boundaries

Abstract

A fundamental question in microbiology is whether there is a continuum of genetic diversity among genomes or clear species boundaries prevail instead. Answering this question requires robust measurement of whole-genome relatedness among thousands of genomes and from diverge phylogenetic lineages. Whole-genome similarity metrics such as Average Nucleotide Identity (ANI) can provide the resolution needed for this task, overcoming several limitations of traditional techniques used for the same purposes. Although the number of genomes currently available may be adequate, the associated bioinformatics tools for analysis are lagging behind these developments and cannot scale to large datasets. Here, we present a new method, FastANI, to compute ANI using alignmentfree approximate sequence mapping. Our analyses demonstrate that FastANI produces an accurate ANI estimate and is up to three orders of magnitude faster when compared to an alignment (e.g., BLAST)-based approach. We leverage FastANI to compute pairwise ANI values among all prokaryotic genomes available in the NCBI database. Our results reveal a clear genetic discontinuity among the database genomes, with 99.8% of the total 8 billion genome pairs analyzed showing either >95% intra-species ANI or <83% inter-species ANI values. We further show that this discontinuity is recovered with or without the most frequently represented species in the database and is robust to historic additions in the public genome databases. Therefore, 95% ANI represents an accurate threshold for demarcating almost all currently named prokaryotic species, and wide species boundaries may exist for prokaryotes. species concept | nucleotide identity | minhash | prokaryotic diversity

Host‐specific competitiveness to form nodules in Rhizobium leguminosarum symbiovar viciae

“…) and digital DNA–DNA hybridisation (dDDH; Fig. S3) comparisons, these isolates are mainly related to three genospecies of R. leguminosarum , gsE (27 strains), gsC (10 strains) and gsB (3 strains), that have been previously identified in European isolates (Kumar et al , ; Jain et al , ; Cavassim et al , ). Eight strains belonged to two other genospecies that were not observed in these previous studies.…”

Microbial composition of Kombucha determined using amplicon sequencing and shotgun metagenomics

“…Contig sequences were aligned to raw data using BWA (Li & Durbin, ) and placed into taxonomic bins with MetaBat2 (Kang, Froula, Egan, & Wang, ) using a minimum contig length of 2,000 bp. The comparison and dereplication of recovered genomes bins were completed using dRep (Olm, Brown, Brooks, & Banfield, ) with minimum primary Average Nucleotide Identity (ANI) 60% (Varghese et al., ) and minimum secondary ANI 95% (Jain, Rodriguez‐R, Phillippy, Konstantinidis, & Aluru, ). CheckM (Parks, Imelfort, Skennerton, Hugenholtz, & Tyson, ) was used to determine completeness and contamination of the best genomes selected by dRep analysis.…”

Switching on the light: using metagenomic shotgun sequencing to characterize the intestinal microbiome of Atlantic cod

“…The relatedness between the 15 reference genomes recruiting the highest portion of reads (Supporting Information Table S7) was then estimated based on Average Nucleotide Identity (ANI) and Mash genome distances using FastANI (ver. 1.1) (Jain et al ., ) and Mash (ver. 2.1) (Ondov et al ., ) (Supporting Information Fig.…”