UFCG: database of universal fungal core genes and pipeline for genome-wide phylogenetic analysis of fungi.

In phylogenomics the evolutionary relationship of organisms is studied by their genomic information. A common approach to phylogenomics is to extract related genes from each organism, build a multiple sequence alignment and then reconstruct evolution relations through a phylogenetic tree. Often a set of highly conserved genes occurring in single-copy, called core genes, are used for this analysis, as they allow efficient automation within a taxonomic clade. Here we introduce the Universal Fungal Core Genes (UFCG) database and pipeline for genome-wide phylogenetic analysis of fungi. The UFCG database consists of 61 curated fungal marker genes, including a novel set of 41 computationally derived core genes and 20 canonical genes derived from literature, as well as marker gene sequences extracted from publicly available fungal genomes. Furthermore, we provide an easy-to-use, fully automated and open-source pipeline for marker gene extraction, training and phylogenetic tree reconstruction. The UFCG pipeline can identify marker genes from genomic, proteomic and transcriptomic data, while producing phylogenies consistent with those previously reported, and is publicly available together with the UFCG database at https://ufcg.steineggerlab.com.

EzBioCloud: a genome-driven database and platform for microbiome identification and discovery.

With the continued evolution of DNA sequencing technologies, the role of genome sequence data has become more integral in the classification and identification of Bacteria and Archaea. Six years after introducing EzBioCloud, an integrated platform representing the taxonomic hierarchy of Bacteria and Archaea through quality-controlled 16S rRNA gene and genome sequences, we present an updated version, that further refines and expands its capabilities. The current update recognizes the growing need for accurate taxonomic information as defining a species increasingly relies on genome sequence comparisons. We also incorporated an advanced strategy for addressing underrepresented or less studied lineages, bolstering the comprehensiveness and accuracy of our database. Our rigorous quality control protocols remain, where whole-genome assemblies from the NCBI Assembly Database undergo stringent screening to remove low-quality sequence data. These are then passed through our enhanced identification bioinformatics pipeline which initiates a 16S rRNA gene similarity search and then calculates the average nucleotide identity (ANI). For genome sequences lacking a 16S rRNA sequence and without a closely related genomic representative for ANI calculation, we apply a different ANI approach using bacterial core genes for improved taxonomic placement (core gene ANI, cgANI). Because of the increase in genome sequences available in NCBI and our newly introduced cgANI method, EzBioCloud now encompasses a total of 109 835 species, of which 21 964 have validly published names. 47 896 are candidate species identified either through 16S rRNA sequence similarity (phylotypes) or through whole genome ANI (genomospecies), and the remaining 39 975 were positioned in the taxonomic tree by cgANI (species clusters). Our EzBioCloud database is accessible at www.ezbiocloud.net/db.

AMDB: a database of animal gut microbial communities with manually curated metadata.

Variations in gut microbiota can be explained by animal host characteristics, including host phylogeny and diet. However, there are currently no databases that allow for easy exploration of the relationship between gut microbiota and diverse animal hosts. The Animal Microbiome Database (AMDB) is the first database to provide taxonomic profiles of the gut microbiota in various animal species. AMDB contains 2530 amplicon data from 34 projects with manually curated metadata. The total data represent 467 animal species and contain 10 478 bacterial taxa. This novel database provides information regarding gut microbiota structures and the distribution of gut bacteria in animals, with an easy-to-use interface. Interactive visualizations are also available, enabling effective investigation of the relationship between the gut microbiota and animal hosts. AMDB will contribute to a better understanding of the gut microbiota of animals. AMDB is publicly available without login requirements at http://leb.snu.ac.kr/amdb.

Genomic and evolutionary study from SARS-CoV-2 virus isolates from Bangladesh during the early stage of pandemic strongly correlate with European origin and not with China.

The goal of this study was to identify the genomic variants and determine molecular epidemiology of SARS-CoV-2 virus during the early pandemic stage in Bangladesh. Viral RNA was extracted, converted to cDNA, and amplified using Ion AmpliSeq™ SARS-CoV-2 Research Panel. 413 unique mutants from 151 viral isolates were identified. 80% of cases belongs to 8 mutants: 241C toT, 1163A toT, 3037C toT, 14408C toT, 23403A toG, 28881G toA, 28,882 G toA, and 28883G toC. Observed dominance of GR clade variants that have strong presence in Europe, suggesting European channel a possible entry route. Among 37 genomic mutants significantly associated with clinical symptoms, 3916CtoT (associated with sore-throat), 14408C to T (associated with cough-protection), 28881G to A, 28882G to A, and 28883G to C (associated with chest pain) were notable. These findings may inform future research platforms for disease management and epidemiological study.

Genome insights into the pharmaceutical and plant growth promoting features of the novel species Nocardia alni sp. nov.

BACKGROUND: Recent studies highlighted the biosynthetic potential of nocardiae to produce diverse novel natural products comparable to that of Streptomyces, thereby making them an attractive source of new drug leads. Many of the 119 Nocardia validly named species were isolated from natural habitats but little is known about the diversity and the potential of the endophytic nocardiae of root nodule of actinorhizal plants. RESULTS: The taxonomic status of an actinobacterium strain, designated ncl2T, was established in a genome-based polyphasic study. The strain was Gram-stain-positive, produced substrate and aerial hyphae that fragmented into coccoid and rod-like elements and showed chemotaxonomic properties that were also typical of the genus Nocardia. It formed a distinct branch in the Nocardia 16S rRNA gene tree and was most closely related to the type strains of Nocardia nova (98.6%), Nocardia jiangxiensis (98.4%), Nocardia miyuensis (97.8%) and Nocardia vaccinii (97.7%). A comparison of the draft genome sequence generated for the isolate with the whole genome sequences of its closest phylogenetic neighbours showed that it was most closely related to the N. jiangxiensis, N. miyuensis and N. vaccinii strains, a result underpinned by average nucleotide identity and digital DNA-DNA hybridization data. Corresponding taxogenomic data, including those from a pan-genome sequence analysis showed that strain ncl2T was most closely related to N. vaccinii DSM 43285T. A combination of genomic, genotypic and phenotypic data distinguished these strains from one another. Consequently, it is proposed that strain ncl2T (= DSM 110931T = CECT 30122T) represents a new species within the genus Nocardia, namely Nocardia alni sp. nov. The genomes of the N. alni and N. vaccinii strains contained 36 and 29 natural product-biosynthetic gene clusters, respectively, many of which were predicted to encode for a broad range of novel specialised products, notably antibiotics. Genome mining of the N. alni strain and the type strains of its closest phylogenetic neighbours revealed the presence of genes associated with direct and indirect mechanisms that promote plant growth. The core genomes of these strains mainly consisted of genes involved in amino acid transport and metabolism, energy production and conversion and transcription. CONCLUSIONS: Our genome-based taxonomic study showed that isolate ncl2T formed a new centre of evolutionary variation within the genus Nocardia. This novel endophytic strain contained natural product biosynthetic gene clusters predicted to synthesize novel specialised products, notably antibiotics and genes associated with the expression of plant growth promoting compounds.

Taxonomic composition and variation in the gut microbiota of laboratory mice.

The gut microbiota can affect host health, including humans. Mouse models have been used extensively to study the relationships between the host and the gut microbiota. With the development of cost-effective high-throughput DNA sequencing, several methods have been used to identify members of the gut microbiota of laboratory mice. In recent years, the amount of research and knowledge about the mouse gut microbiota has exploded, leading to significant breakthroughs in understanding of the taxonomic composition of and variation in this community. In addition, the rapidly increasing volume of data has allowed the development of public resources for exploring the mouse gut microbiota. In this review, we describe the concepts and pros and cons of basic methodologies that can be used to determine the gut bacterial profile in laboratory mice. We also present the key bacterial components of the mouse gut microbiota from the phylum to the species level and then compare them with those identified in other references. Additionally, we discuss variations in the mouse gut microbiota and their association with experiments using mice. Finally, we summarize the properties and functions of currently available public resources for exploring the mouse gut microbiota.

Epidemiologic Linkage of COVID-19 Outbreaks at Two University-affiliated Hospitals in the Seoul Metropolitan Area in March 2020.

BACKGROUND: Coronavirus disease 2019 (COVID-19) outbreaks emerged at two university-affiliated hospitals in Seoul (hospital A) and Uijeongbu City (hospital S) in the metropolitan Seoul area in March 2020. The aim of this study was to investigate epidemiological links between the outbreaks using whole genome sequencing (WGS) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). METHODS: Fifteen patients were enrolled in the study, including four non-outbreak (A1-A4) and three outbreak cases (A5-A7) in hospital A and eight cases (S1-S8) in hospital S. Patients' hospital stays, COVID-19 symptoms, and transfer history were reviewed. RNA samples were submitted for WGS and genome-wide single nucleotide variants and phylogenetic relationships were analyzed. RESULTS: The index patient (A5) in hospital A was transferred from hospital S on 26 March. Patients A6 and A7 were the family caregiver and sister, respectively, of the patient who shared a room with A5 for 4 days. Prior to transfer, A5 was at the next bed to S8 in the emergency room on 25 March. Patient S6, a professional caregiver, took care of the patient in the room next to S8's room for 5 days until 22 March and then S5 for another 3 days. WGS revealed that SARS-CoV-2 in A2, A3, and A4 belong to clades V/B.2, S/A, and G/B.1, respectively, whereas that of A5-A7 and S1-S5 are of the V/B.2.1 clade and closely clustered. In particular, SARS-CoV-2 in patients A5 and S5 showed perfect identity. CONCLUSION: WGS is a useful tool to understand epidemiology of SARS-CoV-2. It is the first study to elucidate the role of patient transfer and caregivers as links of nosocomial outbreaks of COVID-19 in multiple hospitals.

Genome-based reclassification of Paenibacillus jamilae Aguilera et al. 2001 as a later heterotypic synonym of Paenibacillus polymyxa (Prazmowski 1880) Ash et al. 1994.

Paenibacillus is one of the genera that has high species diversity and Paenibacillus polymyxa, the type species of the genus, is mainly isolated from plant-associated environments. Among the plant-associated species, Paenibacillus jamilae B.3T (=CECT 5266T=DSM 13815T=KACC 10925T=KCTC 13919T) was proposed to be a novel species according to 16S rRNA gene similarity and DNA-DNA relatedness with related species, including Paenibacillus polymyxa. Nevertheless, in the description of Paenibacillus jamilae the used strain of Paenibacillus polymyxa was not the type strain of this species. In this work we found that the type strains of both species showed 16S rRNA gene similarity of 99.6 %. Therefore, in this study, we sequenced the genome of Paenibacillus jamilae KACC 10925T and compared it with those of the type strain of Paenibacillus polymyxa ATCC 842T and other phylogenetically related species. Genome relatedness value calculated by DNA-DNA hybridization between type strains of Paenibacillus polymyxa and Paenibacillus jamilae was 73.5 %, which is higher than the threshold value (70 %). For more objective and repeatable results of genome relatedness, we analysed an average nucleotide identity (ANI) between two strains. Our results showed that ANI value between the type strains of Paenibacillus jamilae and Paenibacillus polymyxa is 98.5 %, a phylogenetic distance also higher than the threshold values (95~96 %). These values were proposed by Yoon et al. (2017). In addition, their phylogenetic distance based on 92 bacterial core genes is highly close compared to other species. These mean that Paenibacillus jamilae and Paenibacillus polymyxa should be reclassified as a single species. Based on the results from genomic level comparison as well as reexamination results of physiological and chemotaxonomic features, we propose reclassification of Paenibacillus jamilae as a later heterotypic synonym of Paenibacillus polymyxa.

Replication of Vibrio cholerae classical CTX phage.

The toxigenic classical and El Tor biotype Vibrio cholerae serogroup O1 strains are generated by lysogenization of host-type-specific cholera toxin phages (CTX phages). Experimental evidence of the replication and transmission of an El Tor biotype-specific CTX phage, CTX-1, has explained the evolution of V. cholerae El Tor biotype strains. The generation of classical biotype strains has not been demonstrated in the laboratory, and the classical biotype-specific CTX phage, CTX-cla, is considered to be defective with regard to replication. However, the identification of atypical El Tor strains that contain CTX-cla-like phage, CTX-2, indicates that CTX-cla and CTX-2 replicate and can be transmitted to V. cholerae strains. The replication of CTX-cla and CTX-2 phages and the transduction of El Tor biotype strains by various CTX phages under laboratory conditions are demonstrated in this report. We have established a plasmid-based CTX phage replication system that supports the replication of CTX-1, CTX-cla, CTX-2, and CTX-O139. The replication of CTX-2 from the tandem repeat of lysogenic CTX-2 in Wave 2 El Tor strains is also presented. El Tor biotype strains can be transduced by CTX phages in vitro by introducing a point mutation in toxT, the transcriptional activator of the tcp (toxin coregulated pilus) gene cluster and the cholera toxin gene. This mutation also increases the expression of cholera toxin in El Tor strains in a sample single-phase culture. Our results thus constitute experimental evidence of the genetic mechanism of the evolution of V. cholerae.

Genome-based reclassification of Weissella jogaejeotgali as a later heterotypic synonym of Weissella thailandensis.

Members of the genus Weissella are mostly found in fermented plant material. Among the Weissella species, two species, Weissella thailandensis and Weissella jogaejeotgali, were isolated from foods fermented from marine animals. The two species showed a high level of 16S rRNA gene similarity (99.39 %), whereas they exhibited a moderate level of DNA-DNA hybridization relatedness (63.9 %) in an earlier study. In this study, we determined the whole genome sequence of W. thailandensis KCTC 3751T and compared it to those of W. jogaejeotgali FOL01T and other related species. The average nucleotide identity value between the type strains of W. thailandensis and W. jogaejeotgali was 96.4 %, which is clearly higher than the cut-off proposed for bacterial species. We, therefore, propose to reclassify W. jogaejeotgali as a later heterotypic synonym of W. thailandensis.

Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes.

Advancement of DNA sequencing technology allows the routine use of genome sequences in the various fields of microbiology. The information held in genome sequences proved to provide objective and reliable means in the taxonomy of prokaryotes. Here, we describe the minimal standards for the quality of genome sequences and how they can be applied for taxonomic purposes.

Deep-sea hydrothermal vent bacteria related to human pathogenic Vibrio species.

Vibrio species are both ubiquitous and abundant in marine coastal waters, estuaries, ocean sediment, and aquaculture settings worldwide. We report here the isolation, characterization, and genome sequence of a novel Vibrio species, Vibrio antiquarius, isolated from a mesophilic bacterial community associated with hydrothermal vents located along the East Pacific Rise, near the southwest coast of Mexico. Genomic and phenotypic analysis revealed V. antiquarius is closely related to pathogenic Vibrio species, namely Vibrio alginolyticus, Vibrio parahaemolyticus, Vibrio harveyi, and Vibrio vulnificus, but sufficiently divergent to warrant a separate species status. The V. antiquarius genome encodes genes and operons with ecological functions relevant to the environment conditions of the deep sea and also harbors factors known to be involved in human disease caused by freshwater, coastal, and brackish water vibrios. The presence of virulence factors in this deep-sea Vibrio species suggests a far more fundamental role of these factors for their bacterial host. Comparative genomics revealed a variety of genomic events that may have provided an important driving force in V. antiquarius evolution, facilitating response to environmental conditions of the deep sea.

Genomic encyclopedia of bacteria and archaea: sequencing a myriad of type strains.

Microbes hold the key to life. They hold the secrets to our past (as the descendants of the earliest forms of life) and the prospects for our future (as we mine their genes for solutions to some of the planet's most pressing problems, from global warming to antibiotic resistance). However, the piecemeal approach that has defined efforts to study microbial genetic diversity for over 20 years and in over 30,000 genome projects risks squandering that promise. These efforts have covered less than 20% of the diversity of the cultured archaeal and bacterial species, which represent just 15% of the overall known prokaryotic diversity. Here we call for the funding of a systematic effort to produce a comprehensive genomic catalog of all cultured Bacteria and Archaea by sequencing, where available, the type strain of each species with a validly published name (currently∼11,000). This effort will provide an unprecedented level of coverage of our planet's genetic diversity, allow for the large-scale discovery of novel genes and functions, and lead to an improved understanding of microbial evolution and function in the environment.

Rejection of reclassification of Lactobacillus kimchii and Lactobacillus bobalius as later subjective synonyms of Lactobacillus paralimentarius using comparative genomics.

Lactobacillus bobalius, Lactobacillus kimchii and Lactobacillus paralimentarius belong to the genus Lactobacillus and show close phylogenetic relationships. In a previous study, L. bobalius and L. kimchii were proposed to be reclassified as later heterotypic synonyms of L. paralimentarius using high 16S rRNA gene sequence similarities (≥99.5 %) and DNA-DNA hybridization values (≥82 %). We determined high quality whole genome assemblies of the type strains of L. bobalius and L. kimchii, which were then compared with that of L. paralimentarius. Average nucleotide identity values among three genomes ranged from 91.4 to 92.3 % which are clearly below 95~96 %, the generally recognized cutoff value for bacterial species boundaries. On the basis of comparative genomic evidence, L. bobalius, L. kimchii, and L. paralimentarius should stand as separate species in the genus Lactobacillus. We therefore suggest rejecting the previous proposal to combine these three species into a single species.

ContEst16S: an algorithm that identifies contaminated prokaryotic genomes using 16S RNA gene sequences.

Thanks to the recent advancement of DNA sequencing technology, the cost and time of prokaryotic genome sequencing have been dramatically decreased. It has repeatedly been reported that genome sequencing using high-throughput next-generation sequencing is prone to contaminations due to its high depth of sequencing coverage. Although a few bioinformatics tools are available to detect potential contaminations, these have inherited limitations as they only use protein-coding genes. Here we introduce a new algorithm, called ContEst16S, to detect potential contaminations using 16S rRNA genes from genome assemblies. We screened 69 745 prokaryotic genomes from the NCBI Assembly Database using ContEst16S and found that 594 were contaminated by bacteria, human and plants. Of the predicted contaminated genomes, 8 % were not predicted by the existing protein-coding gene-based tool, implying that both methods can be complementary in the detection of contaminations. A web-based service of the algorithm is available at www.ezbiocloud.net/tools/contest16s.

Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies.

The recent advent of DNA sequencing technologies facilitates the use of genome sequencing data that provide means for more informative and precise classification and identification of members of the Bacteria and Archaea. Because the current species definition is based on the comparison of genome sequences between type and other strains in a given species, building a genome database with correct taxonomic information is of paramount need to enhance our efforts in exploring prokaryotic diversity and discovering novel species as well as for routine identifications. Here we introduce an integrated database, called EzBioCloud, that holds the taxonomic hierarchy of the Bacteria and Archaea, which is represented by quality-controlled 16S rRNA gene and genome sequences. Whole-genome assemblies in the NCBI Assembly Database were screened for low quality and subjected to a composite identification bioinformatics pipeline that employs gene-based searches followed by the calculation of average nucleotide identity. As a result, the database is made of 61 700 species/phylotypes, including 13 132 with validly published names, and 62 362 whole-genome assemblies that were identified taxonomically at the genus, species and subspecies levels. Genomic properties, such as genome size and DNA G+C content, and the occurrence in human microbiome data were calculated for each genus or higher taxa. This united database of taxonomy, 16S rRNA gene and genome sequences, with accompanying bioinformatics tools, should accelerate genome-based classification and identification of members of the Bacteria and Archaea. The database and related search tools are available at www.ezbiocloud.net/.

Dehalogenimonas formicexedens sp. nov., a chlorinated alkane-respiring bacterium isolated from contaminated groundwater.

A strictly anaerobic, Gram-stain-negative, non-spore-forming bacterium designated NSZ-14T, isolated from contaminated groundwater in Louisiana (USA), was characterized using a polyphasic approach. Strain NSZ-14T reductively dehalogenated a variety of polychlorinated aliphatic alkanes, producing ethene from 1,2-dichloroethane, propene from 1,2-dichloropropane, a mixture of cis- and trans-1,2-dichloroethene from 1,1,2,2-tetrachloroethane, vinyl chloride from 1,1,2-trichloroethane and allyl chloride (3-chloro-1-propene) from 1,2,3-trichloropropane. Formate or hydrogen could both serve as electron donors. Dechlorination occurred between pH 5.5 and 7.5 and over a temperature range of 20-37 °C. Major cellular fatty acids included C18 : 1ω9c, C14 : 0 and C16 : 0. 16S rRNA gene sequence-based phylogenetic analysis indicated that the strain clusters within the class Dehalococcoidia of the phylum Chloroflexi, most closely related to but distinct from type strains of the species Dehalogenimonas alkenigignens (97.63 % similarity) and Dehalogenimonas lykanthroporepellens (95.05 %). A complete genome sequence determined for strain NSZ-14T revealed a DNA G+C content of 53.96 mol%, which was corroborated by HPLC (54.1±0.2 mol% G+C). Genome-wide comparisons based on average nucleotide identity by orthology and estimated DNA-DNA hybridization values combined with phenotypic and chemotaxonomic traits and phylogenetic analysis indicate that strain NSZ-14T represents a novel species within the genus Dehalogenimonas, for which the name Dehalogenimonas formicexedens sp. nov. is proposed. The type strain is NSZ-14T (=HAMBI 3672T=JCM 19277T=VKM B-3058T). An emended description of Dehalogenimonas alkenigignens is also provided.

Evidence for several waves of global transmission in the seventh cholera pandemic.

Vibrio cholerae is a globally important pathogen that is endemic in many areas of the world and causes 3-5 million reported cases of cholera every year. Historically, there have been seven acknowledged cholera pandemics; recent outbreaks in Zimbabwe and Haiti are included in the seventh and ongoing pandemic. Only isolates in serogroup O1 (consisting of two biotypes known as 'classical' and 'El Tor') and the derivative O139 can cause epidemic cholera. It is believed that the first six cholera pandemics were caused by the classical biotype, but El Tor has subsequently spread globally and replaced the classical biotype in the current pandemic. Detailed molecular epidemiological mapping of cholera has been compromised by a reliance on sub-genomic regions such as mobile elements to infer relationships, making El Tor isolates associated with the seventh pandemic seem superficially diverse. To understand the underlying phylogeny of the lineage responsible for the current pandemic, we identified high-resolution markers (single nucleotide polymorphisms; SNPs) in 154 whole-genome sequences of globally and temporally representative V. cholerae isolates. Using this phylogeny, we show here that the seventh pandemic has spread from the Bay of Bengal in at least three independent but overlapping waves with a common ancestor in the 1950s, and identify several transcontinental transmission events. Additionally, we show how the acquisition of the SXT family of antibiotic resistance elements has shaped pandemic spread, and show that this family was first acquired at least ten years before its discovery in V. cholerae.

A large-scale evaluation of algorithms to calculate average nucleotide identity.

Average nucleotide identity (ANI) is a category of computational analysis that can be used to define species boundaries of Archaea and Bacteria. Calculating ANI usually involves the fragmentation of genome sequences, followed by nucleotide sequence search, alignment, and identity calculation. The original algorithm to calculate ANI used the BLAST program as its search engine. An improved ANI algorithm, called OrthoANI, was developed to accommodate the concept of orthology. Here, we compared four algorithms to compute ANI, namely ANIb (ANI algorithm using BLAST), ANIm (ANI using MUMmer), OrthoANIb (OrthoANI using BLAST) and OrthoANIu (OrthoANI using USEARCH) using >100,000 pairs of genomes with various genome sizes. By comparing values to the ANIb that is considered a standard, OrthoANIb and OrthoANIu exhibited good correlation in the whole range of ANI values. ANIm showed poor correlation for ANI of <90%. ANIm and OrthoANIu runs faster than ANIb by an order of magnitude. When genomes that are larger than 7 Mbp were analysed, the run-times of ANIm and OrthoANIu were shorter than that of ANIb by 53- and 22-fold, respectively. In conclusion, ANI calculation can be greatly sped up by the OrthoANIu method without losing accuracy. A web-service that can be used to calculate OrthoANIu between a pair of genome sequences is available at http://www.ezbiocloud.net/tools/ani . For large-scale calculation and integration in bioinformatics pipelines, a standalone JAVA program is available for download at http://www.ezbiocloud.net/tools/orthoaniu .

OrthoANI: An improved algorithm and software for calculating average nucleotide identity.

Species demarcation in Bacteria and Archaea is mainly based on overall genome relatedness, which serves a framework for modern microbiology. Current practice for obtaining these measures between two strains is shifting from experimentally determined similarity obtained by DNA-DNA hybridization (DDH) to genome-sequence-based similarity. Average nucleotide identity (ANI) is a simple algorithm that mimics DDH. Like DDH, ANI values between two genome sequences may be different from each other when reciprocal calculations are compared. We compared 63 690 pairs of genome sequences and found that the differences in reciprocal ANI values are significantly high, exceeding 1 % in some cases. To resolve this problem of not being symmetrical, a new algorithm, named OrthoANI, was developed to accommodate the concept of orthology for which both genome sequences were fragmented and only orthologous fragment pairs taken into consideration for calculating nucleotide identities. OrthoANI is highly correlated with ANI (using BLASTn) and the former showed approximately 0.1 % higher values than the latter. In conclusion, OrthoANI provides a more robust and faster means of calculating average nucleotide identity for taxonomic purposes. The standalone software tools are freely available at http://www.ezbiocloud.net/sw/oat.