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.     

PhytoREF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy

Photosynthetic eukaryotes have a critical role as the main producers in most ecosystems of the biosphere. The ongoing environmental metabarcoding revolution opens the perspective for holistic ecosystems biological studies of these organisms, in particular the unicellular microalgae that often lack distinctive morphological characters and have complex life cycles. To interpret environmental sequences, metabarcoding necessarily relies on taxonomically curated databases containing reference sequences of the targeted gene (or barcode) from identified organisms. To date, no such reference framework exists for photosynthetic eukaryotes. In this study, we built the PhytoREF database that contains 6490 plastidial 16S rDNA reference sequences that originate from a large diversity of eukaryotes representing all known major photosynthetic lineages. We compiled 3333 amplicon sequences available from public databases and 879 sequences extracted from plastidial genomes, and generated 411 novel sequences from cultured marine microalgal strains belonging to different eukaryotic lineages. A total of 1867 environmental Sanger 16S rDNA sequences were also included in the database. Stringent quality filtering and a phylogeny‐based taxonomic classification were applied for each 16S rDNA sequence. The database mainly focuses on marine microalgae, but sequences from land plants (representing half of the PhytoREF sequences) and freshwater taxa were also included to broaden the applicability of PhytoREF to different aquatic and terrestrial habitats. PhytoREF, accessible via a web interface ( http://phytoref.fr ), is a new resource in molecular ecology to foster the discovery, assessment and monitoring of the diversity of photosynthetic eukaryotes using high‐throughput sequencing.     

16S rRNA gene and 16S-23S rRNA gene internal transcribed spacer sequences analysis of the genus Myxococcus.

Phylogenetic relationships of the species belonging to the genus Myxococcus were elucidated based on the sequences of 16S rRNA genes and 16S-23S rRNA gene internal transcribed spacer (ITS) regions. The Myxococcus species were consequently classified into four distinct groups. The type strain of Myxococcus coralloides occupied an independent position (Group 1); it has been recently reclassified as Corallococcus coralloides. Group 2 comprised the type strains of both Myxococcus virescens and Myxococcus xanthus, and some strains assigned to Myxococcus flavescens. The type strain of M. flavescens was contained in Group 3 along with the strains of Myxococcus fulvus. Group 4 included the strains belonging to C. coralloides, M. fulvus, and M. stipitatus. The type strain of M. fulvus that was allocated outside Group 4 in the 16S rRNA gene tree belonged to Group 3 in the ITS tree. These results strongly suggest that the morphological characteristics of Myxococcus species are not consistent with the phylogenetic relationships. The Myxococcus species must therefore be redefined according to the phylogenetic relationships revealed in this study.     

Release LTP_12_2020, featuring a new ARB alignment and improved 16S rRNA tree for prokaryotic type strains

The new release of the All-Species Living Tree Project (LTP) represents an important step forward in the reconstruction of 16S rRNA gene phylogenies, since we not only provide an updated set of type strain sequences until December 2020, but also a series of improvements that increase the quality of the database. An improved universal alignment has been introduced that is implemented in the ARB format. In addition, all low-quality sequences present in the previous releases have been substituted by new entries with higher quality, many of them as a result of whole genome sequencing. Altogether, the improvements in the dataset and 16S rRNA sequence alignment allowed us to reconstruct robust phylogenies. The trees made available through this current LTP release feature the best topologies currently achievable. The given nomenclature and taxonomic hierarchy reflect all the changes available up to December 2020. The aim is to regularly update the validly published nomenclatural classification changes and new taxa proposals. The new release can be found at the following URL: https://imedea.uib-csic.es/mmg/ltp/.     

Genome size of Eperythrozoon suis and hybridization with 16S rRNA gene

The genome size of Eperythrozoon suis, an unculturable haemotropic mycoplasma, was estimated using pulsed-field gel electrophoresis (PFGE). Gamma irradiation was used to introduce one (on the average) double-strand break in the E. suis Illinois chromosome. Restriction enzymes that cut infrequently were also used to analyze genome size. The size estimate for the full-length genome was 745 kilobases (kb), whereas the size estimates based on the summation of restriction fragments ranged from 730 to 770 kb. The 16S rRNA gene was located on the 120-kb MluI fragment, 128-kb NruI fragment, 25-kb SacII fragment, and 217-kb SalI fragment by Southern blotting.     

Improved group‐specific primers based on the full SILVA 16S rRNA gene reference database

Quantitative PCR (qPCR) and community fingerprinting methods, such as the Terminal Restriction Fragment Length Polymorphism (T‐RFLP) analysis, are well‐suited techniques for the examination of microbial community structures. The use of phylum‐ and class‐specific primers can provide enhanced sensitivity and phylogenetic resolution as compared with domain‐specific primers. To date, several phylum‐ and class‐specific primers targeting the 16S ribosomal RNA gene have been published. However, many of these primers exhibit low discriminatory power against non‐target bacteria in PCR. In this study, we evaluated the precision of certain published primers in silico and via specific PCR. We designed new qPCR and T‐RFLP primer pairs (for the classes Alphaproteobacteria and Betaproteobacteria, and the phyla Bacteroidetes, Firmicutes and Actinobacteria) by combining the sequence information from a public dataset (SILVA SSU Ref 102 NR) with manual primer design. We evaluated the primer pairs via PCR using isolates of the above‐mentioned groups and via screening of clone libraries from environmental soil samples and human faecal samples. As observed through theoretical and practical evaluation, the primers developed in this study showed a higher level of precision than previously published primers, thus allowing a deeper insight into microbial community dynamics.     

A renaissance for the pioneering 16S rRNA gene

Culture-independent molecular surveys using the 16S rRNA gene have become a mainstay for characterizing microbial community structure over the past quarter century. More recently this approach has been overshadowed by metagenomics, which provides a global overview of a community's functional potential rather than just an inventory of its inhabitants. However, the pioneering 16S rRNA gene is making a comeback in its own right thanks to a number of methodological advancements including higher resolution (more sequences), analysis of multiple related samples (e.g. spatial and temporal series) and improved metadata, and use of metadata. The standard conclusion that microbial ecosystems are remarkably complex and diverse is now being replaced by detailed insights into microbial ecology and evolution based only on this one historically important marker gene.     

Intrageneric structure of the genus Gluconobacter analyzed by the 16S rRNA gene and 16S-23S rRNA gene internal transcribed spacer sequences

Forty-nine strains belonging to the genus Gluconobacter were re-examined with respect to their species identification based on the sequences of the 16S rDNA and 16S-23S rDNA internal transcribed spacer regions (ITS). A phylogenetic tree constructed from the 16S rDNA sequences indicated the presence of five clusters corresponding, respectively, to the major five species of the genus Gluconobacter, namely G. albidus, G. cerinus, G. frateurii, G. oxydans (type species), and G. thailandicus. The type strain of G. asaii, NBRC 3276T (T=type strain) was included in the G. cerinus cluster, which is consistent with the report that G. asaii is a junior subjective synonym of G. cerinus. Existence of the G. albidus, G. cerinus, G. frateurii, G. oxydans, and G. thailandicus clusters was also recognized by the ITS sequence analysis. Both sequence analyses revealed that the G. cerinus and G. frateurii clusters were heterogeneous. The G. cerinus cluster comprised three strains of G. cerinus and one strain of G. frateurii, while the G. frateurii cluster included ten strains of G. frateurii, three of G. cerinus, and eleven of G. oxydans. These results suggest that phenotypic differences among Gluconobacter species are ambiguous and the species definition must be re-evaluated. The 16S rDNA and ITS sequences determined in this study are valuable for the identification and phylogenetic analysis of Gluconobacter species.     

一株耐热嗜盐菌的分离及16SrRNA基因序列分析

目的利用盐固体分离培养基,从西藏自治区澜沧江边康宁镇一个47℃的盐井样品中分离纯化到一株耐热嗜盐菌菌株YJ0232。方法通过形态观察、生理生化特性和16srRNA基因序列分析,鉴定嗜盐菌菌株YJ0232分类学地位。结果菌株YJ0232初步鉴定为中度嗜盐菌,属于盐单胞菌属（Halomonassp．）菌株。其16SrRNA基因序列已被GenBank数据库收录,序列号为EU029645。结论本研究对澜沧江高盐环境微生物资源进行了初步探索研究。可为今后研究同类极端环境中新的物种资源以及微生物多样性提供参考。     

Identification of Comamonas species using 16S rRNA gene sequence

A bacterial strain Bz02 was isolated from a water sample collected from river Gomti at the Indian city of Lucknow. We characterized the strain using 16S rRNA sequence. Phylogenetic analysis showed that the strain formed a monophyletic clade with members of the genus Comamonas. The closest phylogenetic relative was Comamonas testosteroni with 95% 16S rRNA gene sequence similarity. It is proposed that the identified strain Bz02 be assigned as the type strain of a species of the genus Comamonas (Comamonas sp Bz02) based on 16S rRNA gene sequence search in Ribosomal Database Project, small subunit rRNA and large subunit rRNA databases together with the phylogenetic tree analysis. The sequence is deposted in GenBank with the accession number {"type":"entrez-nucleotide","attrs":{"text":"FJ211417","term_id":"209164434"}}FJ211417.     

Multiple copies of 16S rRNA gene affect the restriction patterns and DGGE profile 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 DGGE 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 CSP6I+HinfI, MspI+RsaI or HhaI as analyzed in silico. Polymorphisms of the multiple copies in certain strains were further characterized by G+C% and phy-logentic 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.     

C16S - a Hidden Markov Model based algorithm for taxonomic classification of 16S rRNA gene sequences

Recent advances in high throughput sequencing technologies and concurrent refinements in 16S rDNA isolation techniques have facilitated the rapid extraction and sequencing of 16S rDNA content of microbial communities. The taxonomic affiliation of these 16S rDNA fragments is subsequently obtained using either BLAST-based or word frequency based approaches. However, the classification accuracy of such methods is observed to be limited in typical metagenomic scenarios, wherein a majority of organisms are hitherto unknown. In this study, we present a 16S rDNA classification algorithm, called C16S, that uses genus-specific Hidden Markov Models for taxonomic classification of 16S rDNA sequences. Results obtained using C16S have been compared with the widely used RDP classifier. The performance of C16S algorithm was observed to be consistently higher than the RDP classifier. In some scenarios, this increase in accuracy is as high as 34%. A web-server for the C16S algorithm is available at http://metagenomics.atc.tcs.com/C16S/.     

Bacterial,archaeal and eukaryotic diversity in Arctic sediment as revealed by 16S rRNA and 18S rRNA gene clone libraries analysis

We studied the microbial diversity in the sediment from the Kongsfjorden, Svalbard, Arctic, in the summer of 2005 based on the analysis of 16S rRNA and 18S rRNA gene clone libraries. The sequences of the cloned 16S rRNA and 18S rRNA gene inserts were used to determine the species identity or closest relatives by comparison with sequences of known species. Compared to the other samples acquired in Arctic and Antarctic, which are different from that of ours, the microbial diversity in our sediment is much higher. The bacterial sequences were grouped into 11 major lineages of the domain Bacteria: Proteobacteria (include α-, β-, γ-, δ-, and ε-Proteobacteria); Bacteroidetes; Fusobacteria; Firmicutes; Chloroflexi; Chlamydiae; Acidobacteria; Actinobacteria; Planctomycetes; Verrucomicrobiae and Lentisphaerae. Crenarchaeota were dominant in the archaeal clones containing inserts. In addition, six groups from eukaryotes including Cercozoa, Fungi, Telonema, Stramenopiles, Alveolata, and Metazoa were identified. Remarkably, the novel group Lentisphaerae was reported in Arctic sediment at the first time. Our study suggested that Arctic sediment as a unique habitat may contain substantial microbial diversity and novel species will be discovered.     

Nucleotide sequence of a Euglena gracilis chloroplast gene coding for the 16S rRNA: homologies to E. coli and Zea mays chloroplast 16S rRNA

The nucleotide sequence of 16S rDNA from Euglena gracilis chloroplasts has been determined representing the first complete sequence of an algal chloroplast rRNA gene. The structural part of the 16S rRNA gene has 1491 nucleotides according to a comparative analysis of our sequencing results with the published 5'- and 3'-terminal "T1-oligonucleotides" from 16S rRNA from E. gracilis. Alignment with 16S rDNA from Zea mays chloroplasts and E. coli reveals 80 to 72% sequence homology, respectively. Two deletions of 9 and 23 nucleotides are found which are identical in size and position with deletions observed in 16S rDNA of maize and tobacco chloroplasts and which seem to be characteristic for all chloroplast rRNA species. We also find insertions and deletions in E. gracilis not seen in 16S rDNA of higher plant chloroplasts. The 16S rRNA sequence of E. gracilis chloroplasts can be folded by base pairing according to the general 16S rRNA secondary structure model.     

Phylogenetic studies of the rRNA group II pseudomonads based on 16S rRNA gene sequences

Nearly complete sequences of 16S rRNA genes were determined for eight bacterial strains representing five species of the rRNA homology group II pseudomonads that are members of the beta subclass of the class Proteobacteria. Comparative analysis with published sequence data indicated that Pseudomonas andropogonis, Ps. caryophylli, Ps. gladioli pv. gladioli and Ps. cepacia aggregate in one coherent cluster at 94·2% sequence similarity; Ps. solanacearum and Ps. pickettii shared 95·3% and 92·8% similarity with Alcaligenes eutrophus in another cluster. Both clusters joined at 87·8% similarity, which is similar to that for genera in this subclass of Proteobacteria. Based on this study and on comparison with other works we suggest that these species are separated from authentic pseudomonads and constitute a new genus or possibly two related genera accommodating Ps. andropogonis, Ps. caryophylli, Ps. gladioli, Ps. cepacia, and Ps. solanacearum, Ps. pickettii and A. eutrophus, respectively. Four strains of Ps. solanacearum representing Biovars 1, 2, 3 and 4 were subdivided into two clusters at 99·1% sequence similarity, in agreement with other published phenotypic and genotypic studies. The two clusters may be potentially regarded as subspecies.     

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.