Design and Experimental Application of a Novel Non-Degenerate Universal Primer Set that Amplifies Prokaryotic 16S rRNA Genes with a Low Possibility to Amplify Eukaryotic rRNA Genes

The deep sequencing of 16S rRNA genes amplified by universal primers has revolutionized our understanding of microbial communities by allowing the characterization of the diversity of the uncultured majority. However, some universal primers also amplify eukaryotic rRNA genes, leading to a decrease in the efficiency of sequencing of prokaryotic 16S rRNA genes with possible mischaracterization of the diversity in the microbial community. In this study, we compared 16S rRNA gene sequences from genome-sequenced strains and identified candidates for non-degenerate universal primers that could be used for the amplification of prokaryotic 16S rRNA genes. The 50 identified candidates were investigated to calculate their coverage for prokaryotic and eukaryotic rRNA genes, including those from uncultured taxa and eukaryotic organelles, and a novel universal primer set, 342F-806R, covering many prokaryotic, but not eukaryotic, rRNA genes was identified. This primer set was validated by the amplification of 16S rRNA genes from a soil metagenomic sample and subsequent pyrosequencing using the Roche 454 platform. The same sample was also used for pyrosequencing of the amplicons by employing a commonly used primer set, 338F-533R, and for shotgun metagenomic sequencing using the Illumina platform. Our comparison of the taxonomic compositions inferred by the three sequencing experiments indicated that the non-degenerate 342F-806R primer set can characterize the taxonomic composition of the microbial community without substantial bias, and is highly expected to be applicable to the analysis of a wide variety of microbial communities.     

OTU Analysis Using Metagenomic Shotgun Sequencing Data

Because of technological limitations, the primer and amplification biases in targeted sequencing of 16S rRNA genes have veiled the true microbial diversity underlying environmental samples. However, the protocol of metagenomic shotgun sequencing provides 16S rRNA gene fragment data with natural immunity against the biases raised during priming and thus the potential of uncovering the true structure of microbial community by giving more accurate predictions of operational taxonomic units (OTUs). Nonetheless, the lack of statistically rigorous comparison between 16S rRNA gene fragments and other data types makes it difficult to interpret previously reported results using 16S rRNA gene fragments. Therefore, in the present work, we established a standard analysis pipeline that would help confirm if the differences in the data are true or are just due to potential technical bias. This pipeline is built by using simulated data to find optimal mapping and OTU prediction methods. The comparison between simulated datasets revealed a relationship between 16S rRNA gene fragments and full-length 16S rRNA sequences that a 16S rRNA gene fragment having a length >150 bp provides the same accuracy as a full-length 16S rRNA sequence using our proposed pipeline, which could serve as a good starting point for experimental design and making the comparison between 16S rRNA gene fragment-based and targeted 16S rRNA sequencing-based surveys possible.     

A nested array of rRNA targeted probes for the detection and identification of enterococci by reverse hybridization

Complete 23S and almost complete 16S rRNA gene sequences were determined for the type strains of the validly described Enterococcus species, Melissococcus pluton and Tetragenococcus halophilus. A comprehensive set of rRNA targeted specific oligonucleotide hybridization probes was designed according to the multiple probe concept. In silico probe design and evaluation was performed using the respective tools of the ARB program package in combination with the ARB databases comprising the currently available 16S as well as 23S rRNA primary structures. The probes were optimized with respect to their application for reverse hybridization in microplate format. The target comprising 16S and 23S rDNA was amplified and labeled by PCR (polymerase chain reaction) using general primers targeting a wide spectrum of bacteria. Alternatively, amplification of two adjacent rDNA fragments of enterococci was performed by using specific primers. In vitro evaluation of the probe set was done including all Enterococcus type strains, and a selection of other representatives of the gram-positive bacteria with a low genomic DNA G+C content. The optimized probe set was used to analyze enriched drinking water samples as well as original samples from waste water treatment plants.     

Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such "heterogeneity hot spots" occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.     

Reconstruction of ribosomal RNA genes from metagenomic data

Direct sequencing of environmental DNA (metagenomics) has a great potential for describing the 16S rRNA gene diversity of microbial communities. However current approaches using this 16S rRNA gene information to describe community diversity suffer from low taxonomic resolution or chimera problems. Here we describe a new strategy that involves stringent assembly and data filtering to reconstruct full-length 16S rRNA genes from metagenomicpyrosequencing data. Simulations showed that reconstructed 16S rRNA genes provided a true picture of the community diversity, had minimal rates of chimera formation and gave taxonomic resolution down to genus level. The strategy was furthermore compared to PCR-based methods to determine the microbial diversity in two marine sponges. This showed that about 30% of the abundant phylotypes reconstructed from metagenomic data failed to be amplified by PCR. Our approach is readily applicable to existing metagenomic datasets and is expected to lead to the discovery of new microbial phylotypes.     

Use of 16S rRNA and rpoB Genes as Molecular Markers for Microbial Ecology Studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such “heterogeneity hot spots” occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.     

Confident phylogenetic identification of uncultured prokaryotes through long read amplicon sequencing of the 16S-ITS-23S rRNA operon

Amplicon sequencing of the 16S rRNA gene is the predominant method to quantify microbial compositions and to discover novel lineages. However, traditional short amplicons often do not contain enough information to confidently resolve their phylogeny. Here we present a cost-effective protocol that amplifies a large part of the rRNA operon and sequences the amplicons with PacBio technology. We tested our method on a mock community and developed a read-curation pipeline that reduces the overall read error rate to 0.18%. Applying our method on four environmental samples, we captured near full-length rRNA operon amplicons from a large diversity of prokaryotes. The method operated at moderately high-throughput (22286–37,850 raw ccs reads) and generated a large amount of putative novel archaeal 23S rRNA gene sequences compared to the archaeal SILVA database. These long amplicons allowed for higher resolution during taxonomic classification by means of long (∼1000 bp) 16S rRNA gene fragments and for substantially more confident phylogenies by means of combined near full-length 16S and 23S rRNA gene sequences, compared to shorter traditional amplicons (250 bp of the 16S rRNA gene). We recommend our method to those who wish to cost-effectively and confidently estimate the phylogenetic diversity of prokaryotes in environmental samples at high throughput.     

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.     

Shotgun metagenomics and microscopy indicate diverse cyanophytes,other bacteria,and microeukaryotes in the epimicrobiota of a northern Chilean wetland Nostoc (Cyanobacteria)

Prokaryotic Nostoc, one of the world's most conspicuous and widespread algal genera (similar to eukaryotic algae, plants, and animals) is known to support a microbiome that influences host ecological roles. Past taxonomic characterizations of surface microbiota (epimicrobiota) of free‐living Nostoc sampled from freshwater systems employed 16S rRNA genes, typically amplicons. We compared taxa identified from 16S, 18S, 23S, and 28S rRNA gene sequences filtered from shotgun metagenomic sequence and used microscopy to illuminate epimicrobiota diversity for Nostoc sampled from a wetland in the northern Chilean Altiplano. Phylogenetic analysis and rRNA gene sequence abundance estimates indicated that the host was related to Nostoc punctiforme PCC 73102. Epimicrobiota were inferred to include 18 epicyanobacterial genera or uncultured taxa, six epieukaryotic algal genera, and 66 anoxygenic bacterial genera, all having average genomic coverage ≥90X. The epicyanobacteria Geitlerinemia, Oscillatoria, Phormidium, and an uncultured taxon were detected only by 16S rRNA gene; Gloeobacter and Pseudanabaena were detected using 16S and 23S; and Phormididesmis, Neosynechococcus, Symphothece, Aphanizomenon, Nodularia, Spirulina, Nodosilinea, Synechococcus, Cyanobium, and Anabaena (the latter corroborated by microscopy), plus two uncultured cyanobacterial taxa (JSC12, O77) were detected only by 23S rRNA gene sequences. Three chlamydomonad and two heterotrophic stramenopiles genera were inferred from 18S; the streptophyte green alga Chaetosphaeridium globosum was detected by microscopy and 28S rRNA genes, but not 18S rRNA genes. Overall, >60% of epimicrobial taxa were detected by markers other than 16S rRNA genes. Some algal taxa observed microscopically were not detected from sequence data. Results indicate that multiple taxonomic markers derived from metagenomic sequence data and microscopy increase epimicrobiota detection.     

人工神经原网络处理PCR数据在细菌鉴定中的应用

目的16SrRNA和16S-23SrRNA间区片段是常用细菌分类鉴定靶点,本研究探讨人工神经原网络（ANN）对上述位点PCR扩增产物数据分析在细菌快速鉴定方面的价值。方法2对15SrRNA基因荧光引物和1对16S-23SrRNA区间基因引物用于扩增血液标本中分离出的317株细菌。相关毛细管电泳（CE）限制性片段长度多态性（RFLP）和单链构象多态性（SSCP）数据进行人工神经原网络分析。结果16S-23SrRNA基因的RFLP数据对未知菌鉴定的准确率高于16SrRNA基因的SSCP数据,分别为98．0％和79．6％。结论实验证明了人工神经原网络作为一种模式识别方法对于简化细菌鉴定十分有价值。     

Discordant 16S and 23S rRNA gene phylogenies for the genus Helicobacter: implications for phylogenetic inference and systematics

Analysis of 16S rRNA gene sequences has become the primary method for determining prokaryotic phylogeny. Phylogeny is currently the basis for prokaryotic systematics. Therefore, the validity of 16S rRNA gene-based phylogenetic analyses is of fundamental importance for prokaryotic systematics. Discrepancies between 16S rRNA gene analyses and DNA-DNA hybridization and phenotypic analyses have been noted in the genus Helicobacter. To clarify these discrepancies, we sequenced the 23S rRNA genes for 55 helicobacter strains representing 41 taxa (>2,700 bases per sequence). Phylogenetic-tree construction using neighbor-joining, parsimony, and maximum likelihood methods for 23S rRNA gene sequence data yielded stable trees which were consistent with other phenotypic and genotypic methods. The 16S rRNA gene sequence-derived trees were discordant with the 23S rRNA gene trees and other data. Discrepant 16S rRNA gene sequence data for the helicobacters are consistent with the horizontal transfer of 16S rRNA gene fragments and the creation of mosaic molecules with loss of phylogenetic information. These results suggest that taxonomic decisions must be supported by other phylogenetically informative macromolecules, such as the 23S rRNA gene, when 16S rRNA gene-derived phylogeny is discordant with other credible phenotypic and genotypic methods. This study found Wolinella succinogenes to branch with the unsheathed-flagellum cluster of helicobacters by 23S rRNA gene analyses and whole-genome comparisons. This study also found intervening sequences (IVSs) in the 23S rRNA genes of strains of 12 Helicobacter species. IVSs were found in helices 10, 25, and 45, as well as between helices 31' and 27'. Simultaneous insertion of IVSs at three sites was found in H. mesocricetorum.     

Accurate taxonomy assignments from 16S rRNA sequences produced by highly parallel pyrosequencers

The recent introduction of massively parallel pyrosequencers allows rapid, inexpensive analysis of microbial community composition using 16S ribosomal RNA (rRNA) sequences. However, a major challenge is to design a workflow so that taxonomic information can be accurately and rapidly assigned to each read, so that the composition of each community can be linked back to likely ecological roles played by members of each species, genus, family or phylum. Here, we use three large 16S rRNA datasets to test whether taxonomic information based on the full-length sequences can be recaptured by short reads that simulate the pyrosequencer outputs. We find that different taxonomic assignment methods vary radically in their ability to recapture the taxonomic information in full-length 16S rRNA sequences: most methods are sensitive to the region of the 16S rRNA gene that is targeted for sequencing, but many combinations of methods and rRNA regions produce consistent and accurate results. To process large datasets of partial 16S rRNA sequences obtained from surveys of various microbial communities, including those from human body habitats, we recommend the use of Greengenes or RDP classifier with fragments of at least 250 bases, starting from one of the primers R357, R534, R798, F343 or F517.     

Discrepancies in the widely applied GAM42a fluorescence in situ hybridisation probe for Gammaproteobacteria

A bacterial culture collection of 104 strains was obtained from an activated sludge wastewater treatment plant to pursue studies into microbial flocculation. Characterisation of the culture collection using a polyphasic approach indicated seven isolates, phylogenetically affiliated with the deep-branching Xanthomonas group of the class Gammaproteobacteria, were unable to hybridise the GAM42a fluorescence in situ hybridisation (FISH) probe for Gammaproteobacteria. The sequence of the GAM42a probe target region in the 23S rRNA gene of these isolates was determined to have mismatches to GAM42a. Probes perfectly targeting the mismatches (GAM42a_T1038_G1031, and GAM42a_T1038 and GAM42a_A1041_A1040) were synthesised, and used in conjunction with GAM42a in FISH to study the Gammaproteobacteria community structure in one full-scale activated sludge plant. Several bacteria in the activated sludge biomass bound the modified probes demonstrating their presence and the fact that these Gammaproteobacteria have been overlooked in community structure analyses of activated sludge.     

Evaluation of 23S rRNA PCR Primers for Use in Phylogenetic Studies of Bacterial Diversity

The availability of a diverse set of 23S rRNA gene sequences enabled evaluation of the specificity of 39 previously published and 4 newly designed primers specific for bacteria. An extensive clone library constructed using an optimized primer pair resulted in similar gene richness but slightly differing coverage of some phylogenetic groups, compared to a 16S rRNA gene library from the same environmental sample.     

Evaluation of 23S rRNA PCR primers for use in phylogenetic studies of bacterial diversity

The availability of a diverse set of 23S rRNA gene sequences enabled evaluation of the specificity of 39 previously published and 4 newly designed primers specific for bacteria. An extensive clone library constructed using an optimized primer pair resulted in similar gene richness but slightly differing coverage of some phylogenetic groups, compared to a 16S rRNA gene library from the same environmental sample.     

Comparative sequence analysis and oligonucleotide probe design based on 23S rRNA genes of Alphaproteobacteria from North Sea bacterioplankton

Almost complete 23S rRNA gene sequences were obtained from 11 Alphaproteobacteria isolated from marine surface water of the German Bight. Five of the strains belong to the "marine alpha" group, a phylogenetic cluster which encompasses members of the genus Roseobacter and closely related bacteria. Phylogenetic sequence analysis based on 52 published as well as unpublished complete 23S rDNA sequences from Alphaproteobacteria including the newly obtained was in general consistent with the 16S rRNA gene sequence-derived phylogeny. 16S and 23S rRNA based phylogenies both showed a distinct cluster for strains associated with the "marine alpha" group. The suitability of both markers for the design of oligonucleotide probes targeting selected groups of Alphaproteobacteria was systematically evaluated and compared in silico. Six clusters of sequences covering different phylogenetic levels as well as two strains were selected in a case study. To compensate for the quantitative difference in the two data sets, the 16S rRNA dataset was truncated to sequences with an equivalent in the 23S rRNA data set. Our results show, that the overall number of phylogenetically redundant probes available could be more than doubled by extending probe design to the 23S rRNA. For small clusters of high sequence similarity and single strains, up to 8 times more discriminating binding sites were provided by the 23S rRNA.     

Alterations in the number of rRNA operons within the Bacillus subtilis genome

Deletions and additions of rRNA gene sets in Bacillus subtilis were observed by Southern hybridizations using cloned radiolabeled rDNA sequences. Of the ten rRNA gene sets found in B. subtilis 168M or NCTC3610, one was deleted in strains possessing the leuB1, ilvC1, argA2 and pheA1 mutations. Among EcoRI restriction fragments of genomic DNA products, a 2.9-kb 23S rRNA homolog was missing. In HindIII digest, both 5.5- and 5.1-kb hybrid bands were lost with 16S and 23S probes, respectively. Similarly, genomic DNAs digested with SmaI showed the absence of both 2.1- and 2.0-kb fragments that hybridized to 16S and 5S sequences, respectively, in wild-type genomes. In contrast, B. subtilis strain 166 and its derivatives displayed a gain of a 3.3-kb HindIII fragment homologous to 16S rRNA. Transforming the ilvC1 and leuB1 mutations into new genetic backgrounds revealed in some clones the concomitant introduction of the ribosomal defect. Transformations with the slightly heterologous donor DNA from strain W23 yielded some Leu+ and Arg+ transformants with altered hybridization patterns when probed with cloned sequences. We propose that the deletion of the rRNA operon occurred in the ilv-leu gene cluster of the B. subtilis genome as a result of unequal recombination between redundant sequences.     

Oligonucleotide primers,probes and molecular methods for the environmental monitoring of methanogenic archaea

For the identification and quantification of methanogenic archaea (methanogens) in environmental samples, various oligonucleotide probes/primers targeting phylogenetic markers of methanogens, such as 16S rRNA, 16S rRNA gene and the gene for the α‐subunit of methyl coenzyme M reductase (mcrA), have been extensively developed and characterized experimentally. These oligonucleotides were designed to resolve different groups of methanogens at different taxonomic levels, and have been widely used as hybridization probes or polymerase chain reaction primers for membrane hybridization, fluorescence in situ hybridization, rRNA cleavage method, gene cloning, DNA microarray and quantitative polymerase chain reaction for studies in environmental and determinative microbiology. In this review, we present a comprehensive list of such oligonucleotide probes/primers, which enable us to determine methanogen populations in an environment quantitatively and hierarchically, with examples of the practical applications of the probes and primers.     

Detection and identification of gastrointestinal Lactobacillus species by using denaturing gradient gel electrophoresis and species-specific PCR primers

Denaturing gradient gel electrophoresis (DGGE) of DNA fragments obtained by PCR amplification of the V2-V3 region of the 16S rRNA gene was used to detect the presence of Lactobacillus species in the stomach contents of mice. Lactobacillus isolates cultured from human and porcine gastrointestinal samples were identified to the species level by using a combination of DGGE and species-specific PCR primers that targeted 16S-23S rRNA intergenic spacer region or 16S rRNA gene sequences. The identifications obtained by this approach were confirmed by sequencing the V2-V3 region of the 16S rRNA gene and by a BLAST search of the GenBank database.