Short- and long-read metabarcoding of the eukaryotic rRNA operon: evaluation of primers and comparison to shotgun metagenomics sequencing

High-throughput sequencing for analysis of environmental microbial diversity has evolved vastly over the last decade. Currently the go-to method for microbial eukaryotes is short-read metabarcoding of variable regions of the 18S rRNA gene with <500 bp amplicons. However, there is a growing interest in long-read sequencing of amplicons covering the rRNA operon for improving taxonomic resolution. For both methods, the choice of primers is crucial. It determines if community members are covered, if they can be identified at a satisfactory taxonomic level, and if the obtained community profile is representative. Here, we designed new primers targeting 18S and 28S rRNA based on 177,934 and 21,072 database sequences, respectively. The primers were evaluated in silico along with published primers on reference sequence databases and marine metagenomics datasets. We further evaluated a subset of the primers for short- and long-read sequencing on environmental samples in vitro and compared the obtained community profile with primer-unbiased metagenomic sequencing. Of the short-read pairs, a new V6-V8 pair and the V4_Balzano pair used with a simplified PCR protocol provided good results in silico and in vitro. Fewer differences were observed between the long-read primer pairs. The long-read amplicons and ITS1 alone provided higher taxonomic resolution than V4. Together, our results represent a reference and guide for selection of robust primers for research on and environmental monitoring of microbial eukaryotes.

Complete Genome Sequence of a Novel Thermoacidophilic and Facultatively Anaerobic Crenarchaeon, Stygiolobus sp. Strain KN-1

The complete genome sequence of the thermoacidophilic crenarchaeon Stygiolobus sp. strain KN-1 was determined and annotated. The genome was 2,958,410 bp in size, with a GC content of 40.1%. It contained 2,973 coding sequences, a single copy of the 16S-23S rRNA operon, 47 tRNA genes, and 9 CRISPR repeat sequences.

Establishment and assessment of an amplicon sequencing method targeting the 16S-ITS-23S rRNA operon for analysis of the equine gut microbiome

Microbial communities are commonly studied by using amplicon sequencing of part of the 16S rRNA gene. Sequencing of the full-length 16S rRNA gene can provide higher taxonomic resolution and accuracy. To obtain even higher taxonomic resolution, with as few false-positives as possible, we assessed a method using long amplicon sequencing targeting the rRNA operon combined with a CCMetagen pipeline. Taxonomic assignment had > 90% accuracy at the species level in a mock sample and at the family level in equine fecal samples, generating similar taxonomic composition as shotgun sequencing. The rRNA operon amplicon sequencing of equine fecal samples underestimated compositional percentages of bacterial strains containing unlinked rRNA genes by a fourth to a third, but unlinked rRNA genes had a limited effect on the overall results. The rRNA operon amplicon sequencing with the A519F + U2428R primer set was able to detect some kind of archaeal genomes such as Methanobacteriales and Methanomicrobiales, whereas full-length 16S rRNA with 27F + 1492R could not. Therefore, we conclude that amplicon sequencing targeting the rRNA operon captures more detailed variations of equine microbiota.

Establishment and Assessment of An Amplicon Sequencing Method Targeting The 16S-ITS-23S rRNA Operon For Analysis of The Equine Gut Microbiome

Microbial communities are commonly studied by using amplicon sequencing of part of the 16S rRNA gene. Sequencing of the full-length 16S rRNA gene can provide higher taxonomic resolution and accuracy. To obtain even higher taxonomic resolution, with as few false-positives as possible, we assessed a method using long amplicon sequencing targeting the rRNA operon combined with a CCMetagen pipeline. Taxonomic assignment had >90% accuracy at the species level in a mock sample and at the family level in equine fecal samples, generating similar taxonomic composition as shotgun sequencing. The rRNA operon amplicon sequencing of equine fecal samples underestimated compositional percentages of bacterial strains containing unlinked rRNA genes by a third to almost a half, but unlinked rRNA genes had a limited effect on the overall results. The rRNA operon amplicon sequencing with the A519F + U2428R primer set was able to reflect archaeal genomes, whereas full-length 16S rRNA with 27F + 1492R could not. Therefore, we conclude that amplicon sequencing targeting the rRNA operon captures more detailed variations of bacterial and archaeal microbiota.

Draft Genome Sequence of Pseudomonas sp. Strain T2.31D-1, Isolated from a Drilling Core Sample Obtained 414 Meters below Surface in the Iberian Pyrite Belt

ABSTRACT We report the draft genome of Pseudomonas sp. strain T2.31D-1, which was isolated from a drilling core sample obtained 414 m below surface in the Iberian Pyrite Belt. The genome consists of a 4.7-Mb chromosome with 4,428 coding sequences, 1 rRNA operon, 59 tRNA genes, and a 31.8-kb plasmid.

Draft Genome Sequence of a Filamentous Anoxygenic Phototrophic Bacterium, “Candidatus Roseilinea sp. Strain NK_OTU-006,” Recovered from Metagenomic Data of a Hot Spring Microbial Mat

ABSTRACT We report here the metagenome-assembled draft genome of an uncultured filamentous anoxygenic phototroph of the phylum Chloroflexota (formerly Chloroflexi), “Candidatus Roseilinea sp. strain NK_OTU-006,” recovered from hot spring-associated microbial mats. The 3.6-Mb genome is estimated to be 94% complete and comprises 117 contigs encoding 3,203 predicted genes, including a full-length rRNA operon.

Enhanced Symbiotic Characteristics in Bacterial Genomes with the Disruption of rRNA Operon

Ribosomal RNA is an indispensable molecule in living organisms that plays an essential role in protein synthesis. Especially in bacteria, 16S, 23S, and 5S rRNAs are usually co-transcribed as operons. Despite the positive effects of rRNA co-transcription on growth and reproduction rate, a recent study revealed that bacteria with unlinked rRNA operons are more widespread than expected. However, it is still unclear why the rRNA operon is broken. Here, we explored rRNA operon linkage status in 15,898 bacterial genomes and investigated whether they have common features or lifestyles; 574 genomes were found to have unlinked rRNA operons and tended to be phylogenetically conserved. Most of them were symbionts and showed enhanced symbiotic genomic features such as reduced genome size and high adenine–thymine (AT) content. In an eggNOG-mapper analysis, they were also found to have significantly fewer genes than rRNA operon-linked bacteria in the “transcription” and “energy production and conversion in metabolism” categories. These genomes also tend to decrease RNases related to the synthesis of ribosomes and tRNA processing. Based on these results, the disruption of the rRNA operon seems to be one of the tendencies associated with the characteristics of bacteria requiring a low dynamic range.

Draft Genome Sequence of Shewanella sp. Strain T2.3D-1.1, Isolated from 121.8 Meters Deep in the Subsurface of the Iberian Pyrite Belt

ABSTRACT Shewanella sp. strain T2.3D-1.1 was isolated from the deep subsurface of the Iberian Pyrite Belt. We report its draft genome sequence, consisting of 49 scaffolds, with a chromosome of ≈4.6 Mb and a 23.8-kb plasmid. The chromosome annotation identified 4,068 coding DNA sequences, 1 rRNA operon, and 108 tRNA genes.

Next-Generation Genome Sequencing of Sedum plumbizincicola Sheds Light on the Structural Evolution of Plastid rRNA Operon and Phylogenetic Implications within Saxifragales

The genus Sedum, with about 470 recognized species, is classified in the family Crassulaceae of the order Saxifragales. Phylogenetic relationships within the Saxifragales are still unresolved and controversial. In this study, the plastome of S. plumbizincicola was firstly presented, with a focus on the structural analysis of rrn operon and phylogenetic implications within the order Saxifragaceae. The assembled complete plastome of S. plumbizincicola is 149,397 bp in size, with a typical circular, double-stranded, and quadripartite structure of angiosperms. It contains 133 genes, including 85 protein-coding genes (PCGs), 36 tRNA genes, 8 rRNA genes, and four pseudogenes (one ycf1, one rps19, and two ycf15). The predicted secondary structure of S. plumbizincicola 16S rRNA includes three main domains organized in 74 helices. Further, our results confirm that 4.5S rRNA of higher plants is associated with fragmentation of 23S rRNA progenitor. Notably, we also found the sequence of putative rrn5 promoter has some evolutionary implications within the order Saxifragales. Moreover, our phylogenetic analyses suggested that S. plumbizincicola had a closer relationship with S. sarmentosum than S. oryzifolium, and supported the taxonomic revision of Phedimus. Our findings of the present study will be useful for further investigation of the evolution of plastid rRNA operon and phylogenetic relationships within Saxifragales.