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.     

Higher-level snake phylogeny inferred from mitochondrial DNA sequences of 12S rRNA and 16S rRNA genes

Portions of two mitochondrial genes (12S and 16S ribosomal RNA) were sequenced to determine the phylogenetic relationships among the major clades of snakes. Thirty-six species, representing nearly all extant families, were examined and compared with sequences of a tuatara and three families of lizards. Snakes were found to constitute a monophyletic group (confidence probability [CP] = 96%), with the scolecophidians (blind snakes) as the most basal lineages (CP = 99%). This finding supports the hypothesis that snakes underwent a subterranean period early in their evolution. Caenophidians (advanced snakes), excluding Acrochordus, were found to be monophyletic (CP = 99%). Among the caenophidians, viperids were monophyletic (CP = 98%) and formed the sister group to the elapids plus colubrids (CP = 94%). Within the viperids, two monophyletic groups were identified: true vipers (CP = 98%) and pit vipers plus Azemiops (CP = 99%). The elapids plus Atractaspis formed a monophyletic clade (CP = 99%). Within the paraphyletic Colubridae, the largely Holarctic Colubrinae was found to be a monophyletic assemblage (CP = 98%), and the Xenodontinae was found to be polyphyletic (CP = 91%). Monophyly of the henophidians (primitive snakes) was neither supported nor rejected because of the weak resolution of relationships among those taxa, except for the clustering of Calabaria with a uropeltid, Rhinophis (CP = 94%).      

Novel 16S rRNA gene sequences retrieved from highly saline brine sediments of Kebrit Deep, Red Sea

In this study, we report on first 16S rRNA gene sequences from highly saline brine sediments taken at a depth of 1,515 m in the Kebrit Deep, northern Red Sea. Microbial DNA extracted directly from the sediments was subjected to PCR amplification with primers specific for bacterial and archaeal 16S rRNA gene sequences. The PCR products were cloned, and a total of 11 (6 bacterial and 5 archaeal) clone types were determined by restriction endonuclease digestion. Phylogenetic analysis revealed that most of the cloned sequences were unique, showing no close association with sequences of cultivated organisms or sequences derived from environmental samples. The bacterial clone sequences form a novel phylogenetic lineage (KB1 group) that branches between the Aquificales and the Thermotogales. The archaeal clone sequences group within the Euryarchaeota. Some of the sequences cluster with the group II and group III uncultivated archaea sequence clones, while two clone groups form separate branches. Our results suggest that hitherto unknown archaea and bacteria may thrive in highly saline brines of the Red Sea under extreme environmental conditions. Received: 5 February 1999 / Accepted: 14 July 1999     

The PCR amplification of non-tuberculous mycobacterial 16S rRNA sequences from soil

Non-tuberculous mycobacteria are free living saprophytic organisms commonly found in soil and water. Some are major causes of opportunistic infection, particularly in immuno-compromised patients, and may influence the efficacy of bacille Calmette-Guérin vaccinations. Many of these organisms are not amenable to culture, so information about their distribution is limited. PCR primers designed to amplify part of the mycobacterial 16S rRNA gene were applied to DNA extracted from cultured organisms and soil. The PCR products from soil contained sequences with similarity to slow growing mycobacteria similar to Mycobacterium lentiflavum, and to fast growing mycobacteria such as the xenobiotic degraders PYR-I and RJGII.     

Novel haloarchaeal 16S rRNA gene sequences from Alpine Permo-Triassic rock salt

Prokaryotic diversity in Alpine salt sediments was investigated by polymerase chain reaction (PCR) amplification of 16S rRNA genes, sequencing of cloned products, and comparisons with culturable strains. DNA was extracted from the residue following filtration of dissolved Permo-Triassic rock salt. Fifty-four haloarchaeal sequences were obtained, which could be grouped into at least five distinct clusters. Similarity values of three clusters to known 16S rRNA genes were less than 90%-95%, suggesting the presence of uncultured novel taxa; two clusters were 98% and 99% similar to isolates from Permo-Triassic or Miocene salt from England and Poland, and to Halobacterium salinarum, respectively. Some rock salt samples, including drilling cores, yielded no amplifiable DNA and no cells or only a few culturable cells. This result suggested a variable distribution of haloarchaea within different strata, probably consistent with the known geologic heterogeneity of Alpine salt deposits. We recently reported identical culturable Halococcus salifodinae strains in Permo-Triassic salt sediments from England, Germany, and Austria; together with the data presented here, those results suggest one plausible scenario to be an ancient continuous hypersaline ocean (Zechstein sea) populated by haloarchaea, whose descendants are found today in the salt sediments. The novelty of the sequences also suggested avoidance of haloarchaeal contaminants during our isolation of strains, preparation of DNA, and PCR reactions.     

A Staphylococcus aureus-specific oligonucleotide probe derived from 16S rRNA gene sequences

Genes encoding 16S rRNA were sequenced from 16 species of Staphylococcus. Sequence analysis highlighted a potential Staph. aureus -specific region and a complementary oligonucleotide probe was synthesized and its specificity tested. Northern blotting indicated molecular specificity, and dot blots to RNA from Staph. aureus, Staph. capitis, Staph. caprae, Staph. carnosus, Staph. caseolyticus, Staph. cohnii, Staph. epidermidis, Staph. gallinarum, Staph. haemolyticus, Staph. hominis, Staph. hyicus, Staph. saprophyticus, Staph. sciuri, Staph. simulans, Staph. warneri and Staph. xylosus indicated species-specificity.     

Identification of Vibrio anguillarum in fish by using partial 16S rRNA sequences and a specific 16S rRNA oligonucleotide probe

16S rRNA from seven different Vibrio anguillarum strains was partially sequenced and compared. From this sequence information we could design a 25-base-long oligonucleotide and use it as a specific probe for identification of V. anguillarum. This was determined by RNA-DNA colony hybridization and slot-blot hybridization. Strong, specific hybridization to the probe was observed for all V. anguillarum strains tested. Furthermore, no cross-hybridization could be seen against five other bacterial species. The detection limit was 5 x 10(3) bacteria per ml. It was even possible to detect V. anguillarum, by slot-blot hybridization, directly in a homogenized kidney from a fish that had died of vibriosis. The partial sequence information revealed small but significant differences between strains of the same species. These sequence differences are sufficiently significant to allow serotyping on the RNA level. Comparing strains of different serotypes revealed a 10-base and an 11-base difference in V. anguillarum serotypes O8 and O9, respectively, in a 122-base partial sequence.     

Identification of Vibrio anguillarum in fish by using partial 16S rRNA sequences and a specific 16S rRNA oligonucleotide probe.

16S rRNA from seven different Vibrio anguillarum strains was partially sequenced and compared. From this sequence information we could design a 25-base-long oligonucleotide and use it as a specific probe for identification of V. anguillarum. This was determined by RNA-DNA colony hybridization and slot-blot hybridization. Strong, specific hybridization to the probe was observed for all V. anguillarum strains tested. Furthermore, no cross-hybridization could be seen against five other bacterial species. The detection limit was 5 x 10(3) bacteria per ml. It was even possible to detect V. anguillarum, by slot-blot hybridization, directly in a homogenized kidney from a fish that had died of vibriosis. The partial sequence information revealed small but significant differences between strains of the same species. These sequence differences are sufficiently significant to allow serotyping on the RNA level. Comparing strains of different serotypes revealed a 10-base and an 11-base difference in V. anguillarum serotypes O8 and O9, respectively, in a 122-base partial sequence.     

Retrieval of nearly complete 16S rRNA gene sequences from environmental DNA following 16S rRNA-based community fingerprinting

16S rRNA-based fingerprinting techniques allow rapid analyses of overall bacterial community structure but suffer from a lack of phylogenetic information hitherto retrievable from the short 16S rRNA gene sequences obtained from excised bands. An approach is presented that allows nearly complete 16S rRNA gene sequences to be retrieved for abundant components of the bacterial community as obtained by the community fingerprint, i.e. those reflected by major fingerprint bands. This was achieved by designing a pair of highly specific primers derived from the sequence of an excised band. Combined with universal 16S rRNA primers, these specific primers were applied directly to environmental DNA serving as template. This procedure allowed the generation of a nearly complete 16S rRNA gene sequence of the target taxon by specific polymerase chain reaction (PCR) followed by cycle sequencing down to a relative abundance of at least 1.5% of the environmental DNA. The procedure was exemplified for an epsilonproteobacterium related to Thiomicrospira denitrificans occurring in the central Baltic Sea. This approach is based only on PCR without any cloning step involved. It allows focussing on specific target taxa and is thus rather efficient. This approach should be applicable in general to 16S rRNA or 16S rRNA gene-based fingerprinting techniques and their respective environmental DNA.     

Comparison of partial 16S rRNA gene sequences obtained from activated sludge bacteria

The cultivated and uncultivated bacterial communities of an activated sludge plant were studied. Two samples were taken and a total of 516 bacterial isolates were classified into groups using their whole-cell protein patterns. The distribution of bacteria into protein-pattern groups differed significantly between the two samples, suggesting variation in culturable bacterial flora. Partial 16S rRNA gene sequences were determined for representatives of the commonest protein-pattern groups. Most of the sequences obtained were previously unknown, but relatively closely related to known sequences of organisms belonging to the α, β or γ subclasses of the proteobacteria, the first two subclasses being predominant. This classification of bacteria isolated on a diluted nutrient-rich medium differed from recent culture-dependent studies using nutrient-rich media. The uncultivated bacterial community was studied by analyzing ten partial 16S rRNA gene sequences cloned directly from activated sludge. None of the cloned sequences was identical to those determined for culturable organisms; or to those in the GenBank database. They were, however, related to the α or β subclasses of the proteobacteria, or to the gram-positive bacteria with a high G+C DNA content. Received: 4 November 1996 / Received revision: 24 February 1997 / Accepted: 28 February 1997     

Comparison of 16S rRNA sequences from the family Pasteurellaceae: phylogenetic relatedness by cluster analysis

The taxonomy of the family Pasteurellaceae has remained controversial despite investigations of biochemistry, serology, and nucleic acid relatedness. In an attempt to resolve some of this confusion, we have partially sequenced the 16S rRNAs of seven members of the family, representing all three genera. The sequences were aligned, similarity scores calculated, and single, average and complete linkage cluster analysis of the resulting distance matrix performed. In this way, an evolutionary branching pattern of these closely related species was reconstructed, and the approximate phylogenetic position of the family determined. Actinobacillus (Haemophilus) actinomycetemcomitans clustered with Haemophilus instead of Actinobacillus, supporting transfer of this species to the genus Haemophilus. Thus cluster analysis of phylogenetic relatedness was found to be particularly useful for studying closely related organisms, and could be performed using a microcomputer.     

Cryptic anuran biodiversity in Bangladesh revealed by mitochondrial 16S rRNA gene sequences

To survey the diversity of anuran species in Bangladesh, we compared mitochondrial 16S rRNA gene sequences (approximately 1.4 kbp) from 107 Bangladesh frog specimens. The results of genetic divergence and phylogenetic analyses incorporating data from related species revealed the occurrence of at least eight cryptic species. Hoplobatrachus tigerinus from two districts diverged considerably, indicating the involvement of a cryptic species. Two Fejervarya sp. (large and medium types) and Hylarana cf. taipehensis formed lineages distinct from related species and are probably new species. Microhyla cf. ornata differed from M. ornata with respect to type locality area and involved two distinct species. In addition, we found that Hylarana sp. and Microhyla sp. did not match congeners examined to date in either morphology or 16S rRNA sequence. The occurrence of M. fissipes was tentatively suggested. Consequently, at least, 19 species were found from Bangladesh in this study. These findings revealed a rich anuran biodiversity in Bangladesh, which is unexpected considering the rather simple topographic features of the country.     

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.     

Identification of Escherichia coli through analysis of 16S rRNA and 16S-23S rRNA internal transcribed spacer region sequences

A bacterial strain, designated BzDS03 was isolated from water sample, collected from Dal Lake Srinagar. The strain was characterized by using 16S ribosomal RNA gene and 16S-23S rRNA internal transcribed spacer region sequences. Phylogenetic analysis showed that 16S rRNA sequence of the isolate formed a monophyletic clade with genera Escherichia. The closest phylogenetic relative was Escherichia coli with 99% 16S rRNA gene sequence similarity. The result of Ribosomal database project's classifier tool revealed that the strain BzDS03 belongs to genera Escherichia.16S rRNA sequence of isolate was deposited in GenBank with accession number FJ961336. Further analysis of 16S-23S rRNA sequence of isolate confirms that the identified strain BzDS03 be assigned as the type strain of Escherichia coli with 98% 16S-23S rRNA sequence similarity. The GenBank accession number allotted for 16S-23S rRNA intergenic spacer sequence of isolate is FJ961337.     

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.     

Characterization of the genus Bifidobacterium by automated ribotyping and 16S rRNA gene sequences

In order to characterize the genus Bifidobacterium, ribopatterns and approximately 500 bp (Escherichia coli positions 27 to 520) of 16S rRNA gene sequences of 28 type strains and 64 reference strains of the genus Bifidobacterium were determined. Ribopatterns obtained from Bifidobacterium strains were divided into nine clusters (clusters I-IX) with a similarity of 60%. Cluster V, containing 17 species, was further subdivided into 22 subclusters with a similarity of 90%. In the genus Bifidobacterium, four groups were shown according to Miyake et al.: (i) the Bifidobacterium longum infantis-longum-suis type group, (ii) the B. catenulatum-pseudocatenulatum group, (iii) the B. gallinarum-saeculare-pullorum group, and (iv) the B. coryneforme-indicum group, which showed higher than 97% similarity of the 16S rRNA gene sequences in each group. Using ribotyping analysis, unique ribopatterns were obtained from these species, and they could be separated by cluster analysis. Ribopatterns of six B. adolescentis strains were separated into different clusters, and also showed diversity in 16S rRNA gene sequences. B. adolescentis consisted of heterogeneous strains. The nine strains of B. pseudolongum subsp. pseudolongum were divided into five subclusters. Each type strain of B. pseudolongum subsp. pseudolongum and B. pseudolongum subsp. globosum and two intermediate groups, which were suggested by Yaeshima et al., consisted of individual clusters. B. animalis subsp. animalis and B. animalis subsp. lactis could not be separated by ribotyping using Eco RI. We conclude that ribotyping is able to provide another characteristic of Bifidobacterium strains in addition to 16S rRNA gene sequence phylogenetic analysis, and this information suggests that ribotyping analysis is a useful tool for the characterization of Bifidobacterium species in combination with other techniques for taxonomic characterization.     

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.