Repetitive Sequences in Plant Nuclear DNA:Types,Distribution, Evolution and Function

Repetitive DNA sequences are a major component of eukaryotic genomes and may account for up to 90% of the genome size. They can be divided into minisatellite, microsatellite and satellite sequences. Satellite DNA sequences are considered to be a fast-evolving component of eukaryotic genomes, comprising tandemly-arrayed, highly-repetitive and highly-conserved monomer sequences. The monomer unit of satellite DNA is 150–400 base pairs(bp) in length.Repetitive sequences may be species- or genus-specific, and may be centromeric or subtelomeric in nature. They exhibit cohesive and concerted evolution caused by molecular drive, leading to high sequence homogeneity. Repetitive sequences accumulate variations in sequence and copy number during evolution, hence they are important tools for taxonomic and phylogenetic studies, and are known as ‘‘tuning knobs' ' in the evolution. Therefore, knowledge of repetitive sequences assists our understanding of the organization, evolution and behavior of eukaryotic genomes. Repetitive sequences have cytoplasmic, cellular and developmental effects and play a role in chromosomal recombination. In the post-genomics era, with the introduction of next-generation sequencing technology, it is possible to evaluate complex genomes for analyzing repetitive sequences and deciphering the yet unknown functional potential of repetitive sequences.     

Genomic evolution of the class Acidithiobacillia: deep-branching Proteobacteria living in extreme acidic conditions

Members of the genus Acidithiobacillus, now ranked within the class Acidithiobacillia, are model bacteria for the study of chemolithotrophic energy conversion under extreme conditions. Knowledge of the genomic and taxonomic diversity of Acidithiobacillia is still limited. Here, we present a systematic analysis of nearly 100 genomes from the class sampled from a wide range of habitats. Some of these genomes are new and others have been reclassified on the basis of advanced genomic analysis, thus defining 19 Acidithiobacillia lineages ranking at different taxonomic levels. This work provides the most comprehensive classification and pangenomic analysis of this deep-branching class of Proteobacteria to date. The phylogenomic framework obtained illuminates not only the evolutionary past of this lineage, but also the molecular evolution of relevant aerobic respiratory proteins, namely the cytochrome bo₃ ubiquinol oxidases.Subject terms: Bacterial evolution, Bacterial genomics, Phylogenetics     

A Mitogenomic Phylogeny of Living Primates

Primates, the mammalian order including our own species, comprise 480 species in 78 genera. Thus, they represent the third largest of the 18 orders of eutherian mammals. Although recent phylogenetic studies on primates are increasingly built on molecular datasets, most of these studies have focused on taxonomic subgroups within the order. Complete mitochondrial (mt) genomes have proven to be extremely useful in deciphering within-order relationships even up to deep nodes. Using 454 sequencing, we sequenced 32 new complete mt genomes adding 20 previously not represented genera to the phylogenetic reconstruction of the primate tree. With 13 new sequences, the number of complete mt genomes within the parvorder Platyrrhini was widely extended, resulting in a largely resolved branching pattern among New World monkey families. We added 10 new Strepsirrhini mt genomes to the 15 previously available ones, thus almost doubling the number of mt genomes within this clade. Our data allow precise date estimates of all nodes and offer new insights into primate evolution. One major result is a relatively young date for the most recent common ancestor of all living primates which was estimated to 66-69 million years ago, suggesting that the divergence of extant primates started close to the K/T-boundary. Although some relationships remain unclear, the large number of mt genomes used allowed us to reconstruct a robust primate phylogeny which is largely in agreement with previous publications. Finally, we show that mt genomes are a useful tool for resolving primate phylogenetic relationships on various taxonomic levels.     

EvolMarkers: a database for mining exon and intron markers for evolution, ecology and conservation studies

Recent innovations in next-generation sequencing have lowered the cost of genome projects. Nevertheless, sequencing entire genomes for all representatives in a study remains expensive and unnecessary for most studies in ecology, evolution and conservation. It is still more cost-effective and efficient to target and sequence single-copy nuclear gene markers for such studies. Many tools have been developed for identifying nuclear markers, but most of these have focused on particular taxonomic groups. We have built a searchable database, EvolMarkers, for developing single-copy coding sequence (CDS) and exon-primed-intron-crossing (EPIC) markers that is designed to work across a broad range of phylogenetic divergences. The database is made up of single-copy CDS derived from BLAST searches of a variety of metazoan genomes. Users can search the database for different types of markers (CDS or EPIC) that are common to different sets of input species with different divergence characteristics. EvolMarkers can be applied to any taxonomic group for which genome data are available for two or more species. We included 82 genomes in the first version of EvolMarkers and have found the methods to be effective across Placozoa, Cnidaria, Arthropod, Nematoda, Annelida, Mollusca, Echinodermata, Hemichordata, Chordata and plants. We demonstrate the effectiveness of searching for CDS markers within annelids and show how to find potentially useful intronic markers within the lizard Anolis.     

Using the taxon-specific genes for the taxonomic classification of bacterial genomes

Background

The correct taxonomic assignment of bacterial genomes is a primary and challenging task. With the availability of whole genome sequences, the gene content based approaches appear promising in inferring the bacterial taxonomy. The complete genome sequencing of a bacterial genome often reveals a substantial number of unique genes present only in that genome which can be used for its taxonomic classification.

Results

In this study, we have proposed a comprehensive method which uses the taxon-specific genes for the correct taxonomic assignment of existing and new bacterial genomes. The taxon-specific genes identified at each taxonomic rank have been successfully used for the taxonomic classification of 2,342 genomes present in the NCBI genomes, 36 newly sequenced genomes, and 17 genomes for which the complete taxonomy is not yet known. This approach has been implemented for the development of a tool ‘Microtaxi’ which can be used for the taxonomic assignment of complete bacterial genomes.

Conclusion

The taxon-specific gene based approach provides an alternate valuable methodology to carry out the taxonomic classification of newly sequenced or existing bacterial genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1542-0) contains supplementary material, which is available to authorized users.     

被子植物叶绿体基因组的结构变异研究进展

叶绿体是绿色植物特有的细胞器,其基因组信息被广泛应用于植物系统发育和比较基因组学研究。目前,越来越多的物种有了叶绿体全基因组序列,人们对叶绿体基因组的结构及其变异规律有了更深入的了解。该文对近年来国内外有关被子植物叶绿体基因组插入/缺失、短片段倒位与重复、基因组结构重排以及基因丢失等结构变异式样的研究进展进行综述,并分析了叶绿体基因组结构研究中仍存在的问题以及该领域未来的发展趋势。     

Comparative analysis of the complete mitochondrial genomes of three rockfishes (Scorpaeniformes,Sebastiscus) and insights into the phylogenetic relationships of Sebastidae

Mitochondrial genome is a powerful molecule marker to provide information for phylogenetic relationships and revealing molecular evolution in ichthyological studies. Sebastiscus species, a marine rockfish, are of essential economic value. However, the taxonomic status and phylogenetic relationships of Sebastidae have been controversial so far. Here, the mitochondrial genomes (mitogenomes) of three species, S. tertius, S. albofasciatus, and S. marmoratus, were systemically investigated. The lengths of the mitogenomes’ sequences of S. tertius, S. albofasciatus, and S. marmoratus were 16910, 17056, and 17580 bp, respectively. It contained 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNA (tRNA) genes, and one identical control region (D-loop) among the three species. The genetic distance and K_a/K_s ratio analyses indicated 13 PCGs were suffering purifying selection and the selection pressures were different from certain deep-sea fishes, which were most likely due to the difference in their living environment. The phylogenetic tree was constructed by Bayesian Inference (BI) and Maximum Likelihood (ML). Most interestingly, the results indicated that Sebastidae and Scorpaenidae were grouped into a separate branch, so the taxonomic status of Sebastidae should be classified into subfamily Sebastinae. Our results may lead to a taxonomic revision of Scorpaenoidei.     

Microbial systematics in the post-genomics era

Microbial systematics and phylogeny should form the foundation and guiding light for a comprehensive understanding of different aspects of microbiology. However, there are many critical issues in microbial systematics that are currently not resolved. Some of these include: how to define and delimit a prokaryotic species; development of rationale criteria for the assignment of higher taxonomic ranks; understanding what unique properties distinguish species from different groups; and understanding the branching order and interrelationship among higher prokaryotic clades. The sequencing of genomes from large numbers of cultured as well as uncultured microbes covering prokaryotic diversity provides unique means to achieve these important objectives. Prokaryotic genomes are found to be very diverse and dynamic and horizontal gene transfers (HGTs) are indicated to have played important role in species/genome evolution. Although HGT adds a layer of complexity in terms of understanding the genomes and species evolution, it is contended that vast majority of genes and genetic characteristics that are distinctive characteristics of higher prokaryotic taxa are vertically inherited and based on them a solid foundation for microbial systematics can be developed. We describe two kinds of molecular markers consisting of conserved indels in protein sequences and whole proteins that are specific for different groups that are proving particularly valuable in defining different prokaryotic groups in clear molecular terms and in understanding their interrelationships. The genetic and biochemical studies on these taxa-specific molecular markers also open the way to discover novel biochemical and physiological characteristics that are unique properties of these groups.     

Metagenomic insights into the metabolism and evolution of a new Thermoplasmata order (Candidatus Gimiplasmatales)

Thermoplasmata is a widely distributed and ecologically important archaeal class in the phylum Euryarchaeota. Because few cultures and genomes are available, uncharacterized Thermoplasmata metabolisms remain unexplored. In this study, we obtained four medium- to high-quality archaeal metagenome-assembled genomes (MAGs) from the filamentous fragments of black-odorous aquatic sediments (Foshan, Guangdong, China). Based on their 16S rRNA gene and ribosomal protein phylogenies, the four MAGs belong to the previously unnamed Thermoplasmata UBA10834 clade. We propose that this clade (five reference genomes from the Genome Taxonomy Database (GTDB) and four MAGs from this study) be considered a new order, Candidatus Gimiplasmatales. Metabolic pathway reconstructions indicated that the Ca. Gimiplasmatales MAGs can biosynthesize isoprenoids and nucleotides de novo. Additionally, some taxa have genes for formaldehyde and acetate assimilation, and the Wood–Ljungdahl CO₂-fixation pathway, indicating a mixotrophic lifestyle. Sulfur reduction, hydrogen metabolism, and arsenic detoxification pathways were predicted, indicating sulfur-, hydrogen-, and arsenic-transformation potentials. Comparative genomics indicated that the H₄F Wood–Ljungdahl pathway of both Ca. Gimiplasmatales and Methanomassiliicoccales was likely obtained by the interdomain lateral gene transfer from the Firmicutes. Collectively, this study elucidates the taxonomic and potential metabolic diversity of the new order Ca. Gimiplasmatales and the evolution of this subgroup and its sister lineage Methanomassiliicoccales.     

Intertypic modular exchanges of genomic segments by homologous recombination at universally conserved segments in human adenovirus species D

Human adenovirus species D (HAdV-D), which is composed of clinically and epidemiologically important pathogens worldwide, contains more taxonomic “types” than any other species of the genus Mastadenovirus, although the mechanisms accounting for the high level of diversity remain to be disclosed. Recent studies of known and new types of HAdV-D have indicated that intertypic recombination between distant types contributes to the increasing diversity of the species. However, such findings raise the question as to how homologous recombination events occur between diversified types since homologous recombination is suppressed as nucleotide sequences diverge. In order to address this question, we investigated the distribution of the recombination boundaries in comparison with the landscape of intergenomic sequence conservation assessed according to the synonymous substitution rate (d_S). The results revealed that specific genomic segments are conserved between even the most distantly related genomes; we call these segments “universally conserved segments” (UCSs). These findings suggest that UCSs facilitate homologous recombination, resulting in intergenomic segmental exchanges of UCS-flanking genomic regions as recombination modules. With the aid of such a mechanism, the haploid genomes of HAdV-Ds may have been reshuffled, resulting in chimeric genomes out of diversified repertoires in the HAdV-D population analogous to the MHC region reshuffled via crossing over in vertebrates. In addition, some HAdVs with chimeric genomes may have had the opportunity to avoid host immune responses thereby causing epidemics.     

石松类和蕨类植物质体基因组结构演化研究进展

近年来, 随着测序技术的发展, 石松类和蕨类植物的核基因组、质体基因组以及线粒体基因组研究发展迅速, 质体基因组研究工作更是呈爆发式增长。截至2019年3月1日, GenBank公布的石松类和蕨类植物的175个质体基因组中, 约3/4为最近两年新增。研究内容从早期对个别质体基因组结构和序列特征的简单报道, 逐渐发展到综合性的比较基因组学和系统发育基因组学研究。目前已发表的质体基因组覆盖了石松类和蕨类植物的所有目和大部分科, 这两大类群的质体基因组结构变异和系统发育的基本框架已逐渐清晰。这些研究为我们理解维管植物的早期演化提供了重要参考。本文对石松类和蕨类植物的质体基因组结构特征进行了系统梳理, 发现其结构变异主要包括大片段倒位、IR区边界变动、基因或内含子丢失等, 其中一些结构变异可作为较高分类阶元的共衍征。RNA编辑和长片段非编码序列插入普遍存在于石松类和蕨类植物的质体基因组中, 但其起源、演化机制和功能等仍不清楚。我们对质体基因组的应用、系统发育研究中质体和核基因组的优劣性, 以及系统发育基因组学的前景进行了评述。     

Simple structural differences between coding and noncoding DNA

Background

The study of large-scale genome structure has revealed patterns suggesting the influence of evolutionary constraints on genome evolution. However, the results of these studies can be difficult to interpret due to the conceptual complexity of the analyses. This makes it difficult to understand how observed statistical patterns relate to the physical distribution of genomic elements. We use a simpler and more intuitive approach to evaluate patterns of genome structure.

Methodology/Principal Findings

We used randomization tests based on Morisita''s Index of aggregation to examine average differences in the distribution of purines and pyrimidines among coding and noncoding regions of 261 chromosomes from 223 microbial genomes representing 21 phylum level groups. Purines and pyrimidines were aggregated in the noncoding DNA of 86% of genomes, but were only aggregated in the coding regions of 52% of genomes. Coding and noncoding DNA differed in aggregation in 94% of genomes. Noncoding regions were more aggregated than coding regions in 91% of these genomes. Genome length appears to limit aggregation, but chromosome length does not. Chromosomes from the same species are similarly aggregated despite substantial differences in length. Aggregation differed among taxonomic groups, revealing support for a previously reported pattern relating genome structure to environmental conditions.

Conclusions/Significance

Our approach revealed several patterns of genome structure among different types of DNA, different chromosomes of the same genome, and among different taxonomic groups. Similarity in aggregation among chromosomes of varying length from the same genome suggests that individual chromosome structure has not evolved independently of the general constraints on genome structure as a whole. These patterns were detected using simple and readily interpretable methods commonly used in other areas of biology.     

Intragenomic DNA sequence homologies in the chicken and other members of the classAves: DNA re-association under reduced stringency conditions

Summary We have investigated the intragenomic DNA sequence homologies of twelve species of birds representing five orders, and emphasizing Galliformes. This study differs in two important ways from the classical approaches taken in constructing and evaluating phylogenies based on DNA sequence similarities. Comparisons are made on the basis of sequence homologieswithin genomes of related birds, rather than between genomes. DNA is reassociated at 50°C in 0.5M phosphate buffer; these conditions allow formation and detection of duplexes containing more mismatch than would normally be permitted using more stringent conditions, affording an opportunity to observe more ancient sequence homologies. Thermal stability profiles of DNA duplexes formed under these conditions are the basis of comparison; three general patterns were observed. This approach emphasizes differences in sequence composition between genomes while the more traditional method of intergenomic tracer DNA hybridization at higher stringency emphasizes sequence similarities.No correlation was found between taxonomic position and intragenomic sequence composition, either within or between lineages. The thermal stability profiles of DNA duplexes formed within avian genomes did not reflect the biological similarities inferred from morphology, karyotype, and studies of interspecific hybridization. While all of the differences observed could have occurred over geological time, it was surprising that the genomes of the domestic chicken and the Red Jungle Fowl (Gallus gallus) differ in their sequence compositions. It appears that amplification/reduction events and/or positional changes occur rather often during evolution of a lineage.Abbreviations SDS sodium dodecyl sulphate - PB equimolar sodium phosphate buffer pH 6.8 - Cot concentration of DNA in moles of nucleotide per liter times the incubation time in seconds - Equiv. or Equivalent Cot Cot corrected for the monovalent cation concentration effect on re-association rate - HAP hydroxylapatite - Te_1/2 temperature at which one-half the DNA has eluted from HAP - SSC 0.15M sodium chloride-0.015M sodium citrate     

Identification of novel phospholipase A2 group IX members in metazoans

Sequence homologues of the bacterium Streptomyces violaceoruber and sea anemone Nematostella vectensis PLA₂ pfam09056 members were identified in several bacteria, fungi and metazoans illustrating the evolution of this PLA₂ sub-family. Comparison of their molecular structures revealed that bacteria and fungi members are part of the GXIV of PLA₂s while metazoan representatives are similar with GIX PLA₂ of the marine snail Conus magus. Members of GXIV and GIX PLA₂s show modest overall sequence similarity (21–35%) but considerable motif conservation within the putative Ca²⁺-binding, catalytic sites and cysteine residue positions which are essential for enzyme function. GXIV PLA₂s of bacteria and fungi typically contain four cysteine residues composing two intramolecular disulphide bonds. GIX PLA₂ homologues were identified in cnidarians and molluscs and in a single tunicate but appear to be absent from other metazoan genomes. The mature GIX PLA₂ deduced peptides contain up to ten cysteine residues capable of forming five putative disulphide bonds. Three disulphide bonds were identified in GIX PLA₂s, two of which correspond to those localized in GXIV PLA₂s. Phylogenetic analysis demonstrates that metazoan GIX PLA₂s cluster separate from the bacterial and fungal GXIV PLA₂s and both pfam09056 members form a group separate from the prokaryote and eukaryote GXIIA PLA₂ pfam06951. Duplicate PLA₂ pfam09056 genes were identified in the genomes of sea anemone N. vectensis and oyster Crassostrea gigas suggest that members of this family evolved via species-specific duplication events. These observations indicate that the newly identified metazoan pfam09056 members may be classified as GIX PLA₂s and support the idea of the common evolutionary origin of GXIV and GIX PLA₂ pfam09056 members, which emerged early in bacteria and were maintained in the genomes of fungi and selected extant metazoan taxa.     

Complete Chloroplast Genome of Sedum sarmentosum and Chloroplast Genome Evolution in Saxifragales

Comparative chloroplast genome analyses are mostly carried out at lower taxonomic levels, such as the family and genus levels. At higher taxonomic levels, chloroplast genomes are generally used to reconstruct phylogenies. However, little attention has been paid to chloroplast genome evolution within orders. Here, we present the chloroplast genome of Sedum sarmentosum and take advantage of several available (or elucidated) chloroplast genomes to examine the evolution of chloroplast genomes in Saxifragales. The chloroplast genome of S. sarmentosum is 150,448 bp long and includes 82,212 bp of a large single-copy (LSC) region, 16.670 bp of a small single-copy (SSC) region, and a pair of 25,783 bp sequences of inverted repeats (IRs).The genome contains 131 unique genes, 18 of which are duplicated within the IRs. Based on a comparative analysis of chloroplast genomes from four representative Saxifragales families, we observed two gene losses and two pseudogenes in Paeonia obovata, and the loss of an intron was detected in the rps16 gene of Penthorum chinense. Comparisons among the 72 common protein-coding genes confirmed that the chloroplast genomes of S. sarmentosum and Paeonia obovata exhibit accelerated sequence evolution. Furthermore, a strong correlation was observed between the rates of genome evolution and genome size. The detected genome size variations are predominantly caused by the length of intergenic spacers, rather than losses of genes and introns, gene pseudogenization or IR expansion or contraction. The genome sizes of these species are negatively correlated with nucleotide substitution rates. Species with shorter duration of the life cycle tend to exhibit shorter chloroplast genomes than those with longer life cycles.     

Prokaryotic taxonomy in the sequencing era--the polyphasic approach revisited

The ultimate goal of taxonomy is to establish a system that mirrors the 'order in nature'. In prokaryote microbiology, almost all taxonomic concepts try to mirror the whole evolutionary order back to the origin of life with the cell as basic unit. The introduction of the 16S rRNA gene as molecular marker allowed for the first time the creation of a hierarchical taxonomic system based on one practical molecular marker. With the development of new and rapid sequencing technologies a wealth of new data can and will be used for critical evaluation of the taxonomic system. Comprehensive analyses of other molecular markers as well as total or partial genome comparisons confirmed the 16S rRNA based hierarchical system as 'backbone of prokaryote taxonomy' at least at the genus level and above. A tendency is visible to classify novel taxa more and more based on the genotype, i.e. comparative analyses of 16S rRNA and/or other gene sequence data (in multilocus sequence analysis, MLSA) at the genus and the species level, sometimes contrary to the indications of other (often phenotypic) data. The understanding of all the information behind these data is lagging far behind their accumulation. Genes and genomes do not function on its own and can only display their potential within the cell as the basic unit of evolution (and hence taxonomy). It is the phenotype and the natural selection that 'drive' evolution in a given environment. In this context, the 'polyphasic taxonomic approach' should be revisited again, taking into account the novel insights into genomes and other 'omic' sciences in a more strict and detailed context with the phenotype. This approach allows a more holistic view and provides a sound basis for describing the diversity of prokaryotes and has the potential to become the foundation of a more stable, in-depth taxonomy of the prokaryotes.     

Recombination Rate and AT-content Show Opposite Correlations in Mammalian and Other Animal Genomes

Genomes are often characterized by their genome size, recombination rate and nucleotide content. Negative correlations of recombination rate and AT-content have been reported for selected species. Here we investigate 26 completely sequenced animal genomes with respect to their genomic characteristics. After correcting for phylogenetic correlations we find significant correlations between these characters, showing that these do not evolve independently. The genomes of mammals differ when compared to other animal genomes. Correlations between AT-content and recombination rate are negative, but only for mammals. All other animals show a positive correlation for these characters. These results show that studies on the relationship between genome characteristics using single species might lead to the right conclusion, but cannot be applied as a general rule to all taxonomic groups. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

New features of Asian Crassostrea oyster mitochondrial genomes: A novel alloacceptor tRNA gene recruitment and two novel ORFs

A feasible way to perform evolutionary analyses is to compare characters divergent enough to observe significant differences, but sufficiently similar to exclude saturation of the differences that occurred. Thus, comparisons of invertebrate mitochondrial (mt) genomes at low taxonomic levels can be extremely helpful in investigating patterns of variation and evolutionary dynamics of genomes, as intermediate stages of the process may be identified. Fortunately, in this study, we newly sequenced the mt genome of the eighth member of Asian Crassostrea oysters which can provide necessary intermediate characters for us to believe that the variation of Crassostrea mt genomes is considerably greater than previously acknowledged. Several new features of Asian Crassostrea oyster mitochondrial genomes were revealed, and our results are particularly significant as they 1) suggest a novel model of alloacceptor tRNA gene recruitment, namely "vertical" tRNA gene recruitment, which can be successfully used to explain the origination of the unusually additional trnK and trnQ genes (annotated as trnK(2) and trnQ(2) respectively) in the mt genomes of the five Asian oysters, and we speculate that this recruitment progress may be a common phenomenon in the evolution of the tRNA multigene family; 2) reveal the existence of two additional, lineage-specific, mtDNA-encoded genes that may originate from duplication of nad2 followed by rapid evolutionary change. Each of these two genes encodes a unique amino terminal signal peptide, thus each might possess an unknown function; and 3) identify for the first time the atp8 gene in oysters. The present study thus gives further credence to the comparison of congeneric bivalves as a meaningful strategy to investigate mt genomic evolutionary trends in genome organization, tRNA multigene family, and gene loss and/or duplication that are difficult to undertake at higher taxonomic levels. In particular, our study provides new evidence for the identification and characterization of ORFs in the "non-coding region" of animal mt genomes.     

Strategies to Avoid Wrongly Labelled Genomes Using as Example the Detected Wrong Taxonomic Affiliation for Aeromonas Genomes in the GenBank Database

Around 27,000 prokaryote genomes are presently deposited in the Genome database of GenBank at the National Center for Biotechnology Information (NCBI) and this number is exponentially growing. However, it is not known how many of these genomes correspond correctly to their designated taxon. The taxonomic affiliation of 44 Aeromonas genomes (only five of these are type strains) deposited at the NCBI was determined by a multilocus phylogenetic analysis (MLPA) and by pairwise average nucleotide identity (ANI). Discordant results in relation to taxa assignation were found for 14 (35.9%) of the 39 non-type strain genomes on the basis of both the MLPA and ANI results. Data presented in this study also demonstrated that if the genome of the type strain is not available, a genome of the same species correctly identified can be used as a reference for ANI calculations. Of the three ANI calculating tools compared (ANI calculator, EzGenome and JSpecies), EzGenome and JSpecies provided very similar results. However, the ANI calculator provided higher intra- and inter-species values than the other two tools (differences within the ranges 0.06–0.82% and 0.92–3.38%, respectively). Nevertheless each of these tools produced the same species classification for the studied Aeromonas genomes. To avoid possible misinterpretations with the ANI calculator, particularly when values are at the borderline of the 95% cutoff, one of the other calculation tools (EzGenome or JSpecies) should be used in combination. It is recommended that once a genome sequence is obtained the correct taxonomic affiliation is verified using ANI or a MLPA before it is submitted to the NCBI and that researchers should amend the existing taxonomic errors present in databases.     

Genome-informed Bradyrhizobium taxonomy: where to from here?

Bradyrhizobium is thought to be the largest and most diverse rhizobial genus, but this is not reflected in the number of described species. Although it was one of the first rhizobial genera recognised, its taxonomy remains complex. Various contemporary studies are showing that genome sequence information may simplify taxonomic decisions. Therefore, the growing availability of genomes for Bradyrhizobium will likely aid in the delineation and characterization of new species. In this study, we addressed two aims: first, we reviewed the availability and quality of available genomic resources for Bradyrhizobium. This was achieved by comparing genome sequences in terms of sequencing technologies used and estimated level of completeness for inclusion in genome-based phylogenetic analyses. Secondly, we utilized these genomes to investigate the taxonomic standing of Bradyrhizobium in light of its diverse lifestyles. Although genome sequences differed in terms of their quality and completeness, our data indicate that the use of these genome sequences is adequate for taxonomic purposes. By using these resources, we inferred a fully resolved, well-supported phylogeny. It separated Bradyrhizobium into seven lineages, three of which corresponded to the so-called supergroups known for the genus. Wide distribution of key lifestyle traits such as nodulation, nitrogen fixation and photosynthesis revealed that these traits have complicated evolutionary histories. We present the first robust Bradyrhizobium species phylogeny based on genome sequence information for investigating the evolution of this important assemblage of bacteria. Furthermore, this study provides the basis for using genome sequence information as a resource to make important taxonomic decisions, particularly at the species and genus levels.