Estimation of bacterial species phylogeny through oligonucleotide frequency distances

Classification of bacteria is mainly based on sequence comparisons of certain homologous genes such as 16S rRNA. Recently there are challenges to classify bacteria using oligonucleotide frequency pattern of nonhomologous sequences. However, the evolutionary significance of oligonucleotides longer than tetra-nucleotide is not studied well. We performed phylogenetic analysis by using the Euclidean distances calculated from the di to deca-nucleotide frequencies in bacterial genomes, and compared these oligonucleotide frequency-based tree topologies with those for 16S rRNA gene and concatenated seven genes. When oligonucleotide frequency-based trees were constructed for bacterial species with similar GC content, their topologies at genus and family level were congruent with those based on homologous genes. Our results suggest that oligonucleotide frequency is useful not only for classification of bacteria, but also for estimation of their phylogenetic relationships for closely related species.     

Cyanoseq: A database of cyanobacterial 16S rRNA gene sequences with curated taxonomy

Cyanobacteria are photosynthetic bacteria that occupy various habitats across the globe, playing critical roles in many of Earth's biogeochemical cycles both in both aquatic and terrestrial systems. Despite their well-known significance, their taxonomy remains problematic and is the subject of much research. Taxonomic issues of Cyanobacteria have consequently led to inaccurate curation within known reference databases, ultimately leading to problematic taxonomic assignment during diversity studies. Recent advances in sequencing technologies have increased our ability to characterize and understand microbial communities, leading to the generation of thousands of sequences that require taxonomic assignment. We herein propose CyanoSeq ( https://zenodo.org/record/7569105 ), a database of cyanobacterial 16S rRNA gene sequences with curated taxonomy. The taxonomy of CyanoSeq is based on the current state of cyanobacterial taxonomy, with ranks from the domain to genus level. Files are provided for use with common naive Bayes taxonomic classifiers, such as those included in DADA2 or the QIIME2 platform. Additionally, FASTA files are provided for creation of de novo phylogenetic trees with (near) full-length 16S rRNA gene sequences to determine the phylogenetic relationship of cyanobacterial strains and/or ASV/OTUs. The database currently consists of 5410 cyanobacterial 16S rRNA gene sequences along with 123 Chloroplast, Bacterial, and Vampirovibrionia (formally Melainabacteria) sequences.     

The threonine-tRNA ligase gene region is applicable in classification,typing, and phylogenetic analysis of bifidobacteria

In the modern era, molecular genetic techniques are crucial in ecological studies, as well as in the classification, typing, and phylogenetic analysis of prokaryotes. These techniques are mainly aimed at whole genome comparisons and PCR-derived experiments, including amplifying the 16S rRNA and other various housekeeping genes used in taxonomy, as well as MLST (multilocus sequence typing) and MLSA (multilocus sequence analysis) of different taxonomic bacterial groups. The gene encoding threonine-tRNA ligase (thrS) is a gene potentially applicable as an identification and phylogenetic marker in bacteria. It is widely distributed in bacterial genomes and is subject to evolutionary selection pressure due to its important function in protein synthesis. In this study, specific primers were used to amplify a thrS gene fragment (~740 bp) in 36 type and 30 wild strains classified under family Bifidobacteriaceae. The full-length gene has not yet been considered as a possible identification, classification, and phylogenetic marker in bifidobacteria. The thrS sequences revealed higher sequence variability (82.7% of pairwise identities) among members of the family than that shown by 16S rRNA gene sequences (96.0%). Although discrepancies were found between the thrS-derived and previously reported whole genome phylogenetic analyses, the main phylogenetic groups of bifidobacteria were properly assigned. Most wild strains of bifidobacteria were better differentiated based on their thrS sequences than on their 16S rRNA gene identities. Phylogenetic confidence of the evaluated gene with respect to other alternative genetic markers widely used in taxonomy of bifidobacteria (fusA, GroELhsp60, pyrG, and rplB genes) was confirmed using the localized incongruence difference - Templeton analysis.     

A novel method for determining microflora composition using dynamic phylogenetic analysis of 16S ribosomal RNA deep sequencing data

Deep sequencing of the 16S rRNA gene provides a comprehensive view of bacterial communities in a particular environment and has expanded our ability to study the impact of the microflora on human health and disease. Current analysis methods rely on comparisons of the sequences generated with an expanding but limited set of annotated 16S rRNA sequences or phylogenic clustering of sequences based on arbitrary similarity cutoffs. We describe a novel approach to characterize bacterial composition using deep sequencing of 16S rRNA gene. Our method defines operational taxonomic units based on phylogenetic tree reconstruction and dynamic clustering of sequences using solely sequencing data. These OTUs can be used to identify differences in bacteria abundance between environments. This approach can perform better than previous phylogenetic methods and will significantly improve our understanding of the microfloral role on human diseases by providing a comprehensive analysis of the microbial composition from various bacterial communities.     

Benchmarking of 16S rRNA gene databases using known strain sequences

16S rRNA gene analysis is the most convenient and robust method for microbiome studies. Inaccurate taxonomic assignment of bacterial strains could have deleterious effects as all downstream analyses rely heavily on the accurate assessment of microbial taxonomy. The use of mock communities to check the reliability of the results has been suggested. However, often the mock communities used in most of the studies represent only a small fraction of taxa and are used mostly as validation of sequencing run to estimate sequencing artifacts. Moreover, a large number of databases and tools available for classification and taxonomic assignment of the 16S rRNA gene make it challenging to select the best-suited method for a particular dataset. In the present study, we used authentic and validly published 16S rRNA gene type strain sequences (full length, V3-V4 region) and analyzed them using a widely used QIIME pipeline along with different parameters of OTU clustering and QIIME compatible databases. Data Analysis Measures (DAM) revealed a high discrepancy in ratifying the taxonomy at different taxonomic hierarchies. Beta diversity analysis showed clear segregation of different DAMs. Limited differences were observed in reference data set analysis using partial (V3-V4) and full-length 16S rRNA gene sequences, which signify the reliability of partial 16S rRNA gene sequences in microbiome studies. Our analysis also highlights common discrepancies observed at various taxonomic levels using various methods and databases.     

Phylogenetic clusters of rhizobia revealed by genome structures

Rhizobia, bacteria that fix atmospheric nitrogen, are important agricultural resources. In order to establish the evolutionary relationships among rhizobia isolated from different geographic regions and different plant hosts for systematic studies, we evaluated the use of physical structure of the rhizobial genomes as a phylogenetic marker to categorize these bacteria. In this work, we analyzed the features of genome structures of 64 rhizobial strains. These rhizobial strains were divided into 21 phylogenetic clusters according to the features of genome structures evaluated by the endonuclease I-Ceul. These clusters were supported by 16S rRNA comparisons and genomic sequences of four rhizobial strains, but they are largely different from those based on the current taxonomic scheme (except 16S rRNA).     

VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

Background

The 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.

Results

We developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.

Conclusions

The improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

     

Taxonomic Classification of Bacterial 16S rRNA Genes Using Short Sequencing Reads: Evaluation of Effective Study Designs

Massively parallel high throughput sequencing technologies allow us to interrogate the microbial composition of biological samples at unprecedented resolution. The typical approach is to perform high-throughout sequencing of 16S rRNA genes, which are then taxonomically classified based on similarity to known sequences in existing databases. Current technologies cause a predicament though, because although they enable deep coverage of samples, they are limited in the length of sequence they can produce. As a result, high-throughout studies of microbial communities often do not sequence the entire 16S rRNA gene. The challenge is to obtain reliable representation of bacterial communities through taxonomic classification of short 16S rRNA gene sequences. In this study we explored properties of different study designs and developed specific recommendations for effective use of short-read sequencing technologies for the purpose of interrogating bacterial communities, with a focus on classification using naïve Bayesian classifiers. To assess precision and coverage of each design, we used a collection of ∼8,500 manually curated 16S rRNA gene sequences from cultured bacteria and a set of over one million bacterial 16S rRNA gene sequences retrieved from environmental samples, respectively. We also tested different configurations of taxonomic classification approaches using short read sequencing data, and provide recommendations for optimal choice of the relevant parameters. We conclude that with a judicious selection of the sequenced region and the corresponding choice of a suitable training set for taxonomic classification, it is possible to explore bacterial communities at great depth using current technologies, with only a minimal loss of taxonomic resolution.     

Phylogenetic clusters of rhizobia revealed by genome structures

Rhizobiaareagriculturallyandenvironmentallyimportantbacteria.Theirsymbiosiswithleguminousplantsisresponsibleformostoftheatmosphericni-trogenfixedonland.Classificationofthesebacteriabasedontheirnaturalrelationshipswillpromotetheirapplication.Thisresearchemploysanewphylogeneticmethod,i.e.,revelationandcomparisonofgenomestructure,tocategorizerhizobia.Phylogenyisthestudyoftheevolutionaryrelationshipsamongorgan-isms[1].Currentlyphylogeneticrelationshipsamongrhizobiaaremostlyinferredfromcomparisons…     

Targeted recovery of novel phylogenetic diversity from next-generation sequence data

Next-generation sequencing technologies have led to recognition of a so-called ‘rare biosphere''. These microbial operational taxonomic units (OTUs) are defined by low relative abundance and may be specifically adapted to maintaining low population sizes. We hypothesized that mining of low-abundance next-generation 16S ribosomal RNA (rRNA) gene data would lead to the discovery of novel phylogenetic diversity, reflecting microorganisms not yet discovered by previous sampling efforts. Here, we test this hypothesis by combining molecular and bioinformatic approaches for targeted retrieval of phylogenetic novelty within rare biosphere OTUs. We combined BLASTN network analysis, phylogenetics and targeted primer design to amplify 16S rRNA gene sequences from unique potential bacterial lineages, comprising part of the rare biosphere from a multi-million sequence data set from an Arctic tundra soil sample. Demonstrating the feasibility of the protocol developed here, three of seven recovered phylogenetic lineages represented extremely divergent taxonomic entities. These divergent target sequences correspond to (a) a previously unknown lineage within the BRC1 candidate phylum, (b) a sister group to the early diverging and currently recognized monospecific Cyanobacteria Gloeobacter, a genus containing multiple plesiomorphic traits and (c) a highly divergent lineage phylogenetically resolved within mitochondria. A comparison to twelve next-generation data sets from additional soils suggested persistent low-abundance distributions of these novel 16S rRNA genes. The results demonstrate this sequence analysis and retrieval pipeline as applicable for exploring underrepresented phylogenetic novelty and recovering taxa that may represent significant steps in bacterial evolution.     

Multiple copies of 16s rRNA gene affect the restriction patterns and DGGE profile as revealed by analysis of genome database

The use of 16S rRNA gene has been a "golden" method to determine the diversity of microbial communities in environmental samples, phylogenetic relationships of prokaryotes and taxonomic position of newly isolated organisms. However due to the presence of multiple heterogeneous 16S rRNA gene copies in many strains, the interpretation of microbial ecology via 16S rRNA sequences is complicated. Purpose of present paper is to demonstrate the extent to which the multiple heterogeneous 16S rRNA gene copies affect RFLP patterns and DGG E profiles by using the genome database. In present genome database, there are 782 bacterial strains in total whose genomes have been completely sequenced and annotated. Among the total strains, 639 strains (82%) possess multiple 16S rRNA gene copies, 415 strains (53%) whose multiple copies are heterogeneous in sequences as revealed by alignment, 236 strains (30%) whose multiple copies show different restrict patterns by CSP61 + Hinfl, MspI + Rsal or HhaI as analyzed in silico. Polymorphisms of the multiple copies in certain strains were further characterized by G + C% and phylogentic distances based on the sequences of V3 region, which are linked to DGGE patters. Polymorphisms of a few strains were shown as examples. Using artificial communities, it is demonstrated that the presence of multiple heterogeneous 16S rRNA gene copies potentially leads to over-estimation of the diversity of a community. It is suggested that care must be taken when interpreting 16S rRNA-based RFLP and DGGE data and profiling an environmental community.     

Towards a prokaryotic genomic taxonomy

One of the most interesting developments in the field of modern-day microbiology is the ever increasing number of whole-genome sequences that is publicly available. There is an increasing interest in the use of these genome sequences to assess evolutionary relationships among microbial taxa, as it is anticipated that much additional taxonomic information can be extracted from these sequences. In a first part of the present review, mechanisms that are responsible for the evolution of genomes will be discussed. Subsequently, we will give an overview of approaches that are presently available to assess the taxonomic relationships between prokaryotic species based on complete genome sequences, followed by a brief discussion of the potential implications of these novel approaches for bacterial taxonomy in general and our thinking about the bacterial species concept in particular.     

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.     

Outer membrane efflux protein (OMEP) is a suitable molecular marker for resolving the phylogeny and taxonomic status of closely related cyanobacteria

Taxonomy of Cyanobacteria, the oldest phototrophic prokaryotes, is problematic for many years due to their simple morphology, high variability and adaptability to diverse ecological niches. After introduction of the polyphasic approach, which is based on the combination of several criteria (molecular sequencing, morphological and ecological), the whole classification system of these organisms is subject to reorganization. The aim of this study was to evaluate whether the outer membrane efflux protein (OMEP) sequences can be used as a molecular marker for resolving the phylogeny and taxonomic status of closely related cyanobacteria. We have performed phylogenetic analyses based on the amino acid sequences of the OMEP and the DNA sequences of the 16S rRNA gene from 86 cyanobacterial species/strains with completely sequenced genomes. Phylogenetic trees based on the OMEP showed that most of the cyanobacterial species/strains belonging to different genera are clustered in separate clades supported by high bootstrap values. Comparing the OMEP trees with the 16S rDNA tree clearly showed that the OMEP is more suitable marker in resolving phylogenetic relationships within Cyanobacteria at generic and species level.     

RibAlign: a software tool and database for eubacterial phylogeny based on concatenated ribosomal protein subunits

Background  

Until today, analysis of 16S ribosomal RNA (rRNA) sequences has been the de-facto gold standard for the assessment of phylogenetic relationships among prokaryotes. However, the branching order of the individual phlya is not well-resolved in 16S rRNA-based trees. In search of an improvement, new phylogenetic methods have been developed alongside with the growing availability of complete genome sequences. Unfortunately, only a few genes in prokaryotic genomes qualify as universal phylogenetic markers and almost all of them have a lower information content than the 16S rRNA gene. Therefore, emphasis has been placed on methods that are based on multiple genes or even entire genomes. The concatenation of ribosomal protein sequences is one method which has been ascribed an improved resolution. Since there is neither a comprehensive database for ribosomal protein sequences nor a tool that assists in sequence retrieval and generation of respective input files for phylogenetic reconstruction programs, RibAlign has been developed to fill this gap.     

Phylogenetic analysis of nitrite, nitric oxide, and nitrous oxide respiratory enzymes reveal a complex evolutionary history for denitrification

Denitrification is a facultative respiratory pathway in which nitrite (NO2(-)), nitric oxide (NO), and nitrous oxide (N2O) are successively reduced to nitrogen gas (N(2)), effectively closing the nitrogen cycle. The ability to denitrify is widely dispersed among prokaryotes, and this polyphyletic distribution has raised the possibility of horizontal gene transfer (HGT) having a substantial role in the evolution of denitrification. Comparisons of 16S rRNA and denitrification gene phylogenies in recent studies support this possibility; however, these results remain speculative as they are based on visual comparisons of phylogenies from partial sequences. We reanalyzed publicly available nirS, nirK, norB, and nosZ partial sequences using Bayesian and maximum likelihood phylogenetic inference. Concomitant analysis of denitrification genes with 16S rRNA sequences from the same organisms showed substantial differences between the trees, which were supported by examining the posterior probability of monophyletic constraints at different taxonomic levels. Although these differences suggest HGT of denitrification genes, the presence of structural variants for nirK, norB, and nosZ makes it difficult to determine HGT from other evolutionary events. Additional analysis using phylogenetic networks and likelihood ratio tests of phylogenies based on full-length sequences retrieved from genomes also revealed significant differences in tree topologies among denitrification and 16S rRNA gene phylogenies, with the exception of the nosZ gene phylogeny within the data set of the nirK-harboring genomes. However, inspection of codon usage and G + C content plots from complete genomes gave no evidence for recent HGT. Instead, the close proximity of denitrification gene copies in the genomes of several denitrifying bacteria suggests duplication. Although HGT cannot be ruled out as a factor in the evolution of denitrification genes, our analysis suggests that other phenomena, such gene duplication/divergence and lineage sorting, may have differently influenced the evolution of each denitrification gene.     

基因组时代的放线菌系统分类学及其研究进展

由于放线菌在生物技术、医药卫生和环境治理等领域中的实际重要性,已经成为人们研究得最为深入的微生物类群之一。现代放线菌分类学是建立在对16S rRNA基因保守分子序列的系统发育学分析基础上、综合利用多种微生物信息的多相分类学(polyphasic taxonomy)。随着大规模基因测序技术的飞速发展,已经完成了超过100株放线菌的全基因组序列。近年来发展起来的"基因组系统发育学"有助让人们更加系统地理解放线菌类群的发生、演化,以及不同生境类群之间的相互关系。随着越来越多的以系统发育学为指导的基因组测序研究项目,例如"细菌和古菌基因组百科全书(Genomic Encyclopedia of Bacteria andArchaea,GEBA)计划"的出现,标志着放线菌系统学进入了基因组学研究时代。本文对基因组时代的放线菌系统学方法,以及国内外的相关研究成果进行综述。     

Divergence and redundancy of 16S rRNA sequences in genomes with multiple rrn operons

The level of sequence heterogeneity among rrn operons within genomes determines the accuracy of diversity estimation by 16S rRNA-based methods. Furthermore, the occurrence of widespread horizontal gene transfer (HGT) between distantly related rrn operons casts doubt on reconstructions of phylogenetic relationships. For this study, patterns of distribution of rrn copy numbers, interoperonic divergence, and redundancy of 16S rRNA sequences were evaluated. Bacterial genomes display up to 15 operons and operon numbers up to 7 are commonly found, but ~40% of the organisms analyzed have either one or two operons. Among the Archaea, a single operon appears to dominate and the highest number of operons is five. About 40% of sequences among 380 operons in 76 bacterial genomes with multiple operons were identical to at least one other 16S rRNA sequence in the same genome, and in 38% of the genomes all 16S rRNAs were invariant. For Archaea, the number of identical operons was only 25%, but only five genomes with 21 operons are currently available. These considerations suggest an upper bound of roughly threefold overestimation of bacterial diversity resulting from cloning and sequencing of 16S rRNA genes from the environment; however, the inclusion of genomes with a single rrn operon may lower this correction factor to ~2.5. Divergence among operons appears to be small overall for both Bacteria and Archaea, with the vast majority of 16S rRNA sequences showing <1% nucleotide differences. Only five genomes with operons with a higher level of nucleotide divergence were detected, and Thermoanaerobacter tengcongensis exhibited the highest level of divergence (11.6%) noted to date. Overall, four of the five extreme cases of operon differences occurred among thermophilic bacteria, suggesting a much higher incidence of HGT in these bacteria than in other groups.