Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance

The abundance of different SSU rRNA (“16S”) gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly – from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.     

High-Density Universal 16S rRNA Microarray Analysis Reveals Broader Diversity than Typical Clone Library When Sampling the Environment

Molecular approaches aimed at detection of a broad-range of prokaryotes in the environment routinely rely on classifying heterogeneous 16S rRNA genes amplified by polymerase chain reaction (PCR) using primers with broad specificity. The general method of sampling and categorizing DNA has been to clone then sequence the PCR products. However, the number of clones required to adequately catalog the majority of taxa in a sample is unwieldy. Alternatively, hybridizing target sequences to a universal 16S rRNA gene microarray may provide a more rapid and comprehensive view of prokaryotic community composition. This study investigated the breadth and accuracy of a microarray in detecting diverse 16S rRNA gene sequence types compared to clone-and-sequencing using three environmental samples: urban aerosol, subsurface soil, and subsurface water. PCR products generated from universal 16S rRNA gene-targeted primers were classified by using either the clone-and-sequence method or by hybridization to a novel high-density microarray of 297,851 probes complementary to 842 prokaryotic subfamilies. The three clone libraries comprised 1391 high-quality sequences. Approximately 8% of the clones could not be placed into a known subfamily and were considered novel. The microarray results confirmed the majority of clone-detected subfamilies and additionally demonstrated greater amplicon diversity extending into phyla not observed by the cloning method. Sequences matching operational taxonomic units within the phyla Nitrospira, Planctomycetes, and TM7, which were uniquely detected by the array, were verified with specific primers and subsequent amplicon sequencing. Subfamily richness detected by the array corresponded well with nonparametric richness predictions extrapolated from clone libraries except in the water community where clone-based richness predictions were greatly exceeded. It was concluded that although the microarray is unreliable in identifying novel prokaryotic taxa, it reveals greater diversity in environmental samples than sequencing a typically sized clone library. Furthermore, the microarray allowed samples to be rapidly evaluated with replication, a significant advantage in studies of microbial ecology.     

Evaluation of PCR Primer Selectivity and Phylogenetic Specificity by Using Amplification of 16S rRNA Genes from Betaproteobacterial Ammonia-Oxidizing Bacteria in Environmental Samples

The effect of primer specificity for studying the diversity of ammonia-oxidizing betaproteobacteria (βAOB) was evaluated. βAOB represent a group of phylogenetically related organisms for which the 16S rRNA gene approach is especially suitable. We used experimental comparisons of primer performance with water samples, together with an in silico analysis of published sequences and a literature review of clone libraries made with four specific PCR primers for the βAOB 16S rRNA gene. With four aquatic samples, the primers NitA/NitB produced the highest frequency of ammonia-oxidizing-bacterium-like sequences compared to clone libraries with products amplified with the primer combinations βAMOf/βAMOr, βAMOf/Nso1255g, and NitA/Nso1225g. Both the experimental examination of ammonia-oxidizing-bacterium-specific 16S rRNA gene primers and the literature search showed that neither specificity nor sensitivity of primer combinations can be evaluated reliably only by sequence comparison. Apparently, the combination of sequence comparison and experimental data is the best approach to detect possible biases of PCR primers. Although this study focused on βAOB, the results presented here more generally exemplify the importance of primer selection and potential primer bias when analyzing microbial communities in environmental samples.     

Diversity and Distribution in Hypersaline Microbial Mats of Bacteria Related to Chloroflexus spp.

Filamentous bacteria containing bacteriochlorophylls c and a were enriched from hypersaline microbial mats. Based on phylogenetic analyses of 16S rRNA gene sequences, these organisms form a previously undescribed lineage distantly related to Chloroflexus spp. We developed and tested a set of PCR primers for the specific amplification of 16S rRNA genes from filamentous phototrophic bacteria within the kingdom of “green nonsulfur bacteria.” PCR products recovered from microbial mats in a saltern in Guerrero Negro, Mexico, were subjected to cloning or denaturing gradient gel electrophoresis and then sequenced. We found evidence of a high diversity of bacteria related to Chloroflexus which exhibit different distributions along a gradient of salinity from 5.5 to 16%.     

Review and re-analysis of domain-specific 16S primers

The Polymerase Chain Reaction (PCR) has facilitated the detection of unculturable microorganisms in virtually any environmental source and has thus been used extensively in the assessment of environmental microbial diversity. This technique relies on the assumption that the gene sequences present in the environment are complementary to the "universal" primers used in their amplification. The recent discovery of new taxa with 16S rDNA sequences not complementary to standard universal primers suggests that current 16S rDNA libraries are not representative of true prokaryotic biodiversity. Here we re-assess the specificity of commonly used 16S rRNA gene primers and present these data in tabular form designed as a tool to aid simple analysis, selection and implementation. In addition, we present two new primer pairs specifically designed for effective "universal" Archaeal 16S rDNA sequence amplification. These primers are found to amplify sequences from Crenarchaeote and Euryarchaeote type strains and environmental DNA.     

引物应用策略：基于2对古菌16S rRNA基因通用引物对环境样品古菌多样性分析

【目的】找到适宜的16S rRNA基因通用引物应用策略,应对复杂环境微生物多样性调查,尤其目前高速发展的高通量测序技术带来的巨大挑战。【方法】用Oligocheck软件分别将两对应试的古菌16S rRNA基因通用引物与RDP(Ribosomal database project)数据库中古菌16S rRNA基因序列进行匹配比对。用两对应试引物分别构建海洋沉积物样品的古菌16S rRNA基因文库。【结果】软件匹配结果显示引物f109/r958与目的基因的匹配程度高于引物f21/r958。该结果与古菌16S rRNA基因文库RFLP分析、古菌多样性指数分析结果相吻合。数据还表明,2对引物的综合文库能更好满足该沉积物样品的古菌多样性分析。【结论】选用与数据库中目的基因匹配性高的通用引物和多个引物的联合使用,可以有效提高环境样品微生物多样性调查的分辨率。     

Diversity and geographical distribution of members of the Streptomyces violaceusniger 16S rRNA gene clade detected by clade-specific PCR primers

The Streptomyces violaceusniger 16S rRNA gene clade contains organisms that are of ecological interest and a rich source of novel bioactive metabolites. Improvements in the classification of members of the S. violaceusniger clade made it possible to design, evaluate and use an oligonucleotide primer set to gain an insight into the presence, distribution and taxonomic diversity of members of this taxon in environmental samples. In silico testing showed that the primers had a perfect match with representatives of the S. violaceusniger clade. The primers, designated S-S-Svio-66-a-S-20 and S-S-Svio-1274-a-A-20, amplified an approximately 1190-bp stretch of 16S rRNA gene from authenticated members of the S. violaceusniger clade, but not from representatives of other actinomycete taxa. Following amplification of DNA extracted from sediment and soil samples, the sequences of cloned PCR products confirmed the specific amplification of target sequences in 87% of the clones; the use of 16S rRNA gene fragment similarity correlations indicated that the clones represented new species. The primers can be used to facilitate the isolation of novel members of the S. violaceusniger 16S rRNA gene clade by allowing prescreening of environmental samples and the subsequent detection and retrieval of targetted strains through the use of selective isolation procedures.     

16S rRNA Gene-Based Oligonucleotide Microarray for Environmental Monitoring of the Betaproteobacterial Order “Rhodocyclales”

For simultaneous identification of members of the betaproteobacterial order “Rhodocyclales” in environmental samples, a 16S rRNA gene-targeted oligonucleotide microarray (RHC-PhyloChip) consisting of 79 probes was developed. Probe design was based on phylogenetic analysis of available 16S rRNA sequences from all cultured and as yet uncultured members of the “Rhodocyclales.” The multiple nested probe set was evaluated for microarray hybridization with 16S rRNA gene PCR amplicons from 29 reference organisms. Subsequently, the RHC-PhyloChip was successfully used for cultivation-independent “Rhodocyclales” diversity analysis in activated sludge from an industrial wastewater treatment plant. The implementation of a newly designed “Rhodocyclales”-selective PCR amplification system prior to microarray hybridization greatly enhanced the sensitivity of the RHC-PhyloChip and thus enabled the detection of “Rhodocyclales” populations with relative abundances of less than 1% of all bacteria (as determined by fluorescence in situ hybridization) in the activated sludge. The presence of as yet uncultured Zoogloea-, Ferribacterium/Dechloromonas-, and Sterolibacterium-related bacteria in the industrial activated sludge, as indicated by the RHC-PhyloChip analysis, was confirmed by retrieval of their 16S rRNA gene sequences and subsequent phylogenetic analysis, demonstrating the suitability of the RHC-PhyloChip as a novel monitoring tool for environmental microbiology.     

Developing a Novel Approach of rpoB Gene as a Powerful Biomarker for the Environmental Microbial Diversity

Use of the rpoB gene (the encoding the β-subunit of RNA Polymerase gene), a potential alternative biomarker to the 16S rRNA gene, has been limited to environmental microbial investigation for a long time because of the lack of effective primers. Here we developed a novel rpoB gene-based approach using a newly designed primer pair and tested it in three different environmental water samples with the traditional library method. The results showed that our novel approach presented different microbial diversity patterns from the different environmental samples. Compared to previous rpoB gene based reports, the first retrieved groups in our approach included mainly α-Proteobacteria, δ-Proteobacteria, Planctomycetes, Verrucomicrobia, Firmicutes, Chlorofexi and Actinobacteria, which greatly expanded the potential ability of the rpoB gene approach for environmental surveys. Most importantly, the use of the traditional clone library approach with the novel rpoB gene primers greatly supplement the microbial diversity based on the 16S rRNA gene approach with the universal primer pair (27f and 1492r), at all phylum, class, order, family and genus levels, indicating a powerful complementary method and a potential alternative biomarker of the current popular NGS (next-generation sequences) technologies for the environmental microbial investigation.     

Scratching the surface of the rare biosphere with ribosomal sequence tag primers

Increasingly large datasets of 16S rRNA gene sequences reveal new information about the extent of microbial diversity and the surprising extent of the rare biosphere. Currently, many of the largest datasets are represented by short and variable ribosomal sequence tags (RSTs) that are limited in their ability to accurately assign sequences to broad-scale phylogenetic trees. In this study, we selected 30 rare RSTs from existing sequence datasets and designed primers to amplify c. 1400 bases of the 16S rRNA gene to determine whether these sequences were represented by existing databases or if they might reveal new lineages within the Bacteria. Approximately one-third of the RST primers successfully amplified longer portions of these low-abundance 16S rRNA genes in a specific manner. Subsequent phylogenetic analysis demonstrated that most of these sequences were (1) distantly related to existing cultivated microorganisms and (2) closely related to uncultivated clone sequences that were recently deposited in GenBank. The presence of so many recently collected 16S rRNA gene reference sequences in existing databases suggests that progress is being made quickly towards a microbial census, one which has begun scratching the surface of the 'rare biosphere'.     

Intragenomic Heterogeneity of 16S rRNA Genes Causes Overestimation of Prokaryotic Diversity

Ever since Carl Woese introduced the use of 16S rRNA genes for determining the phylogenetic relationships of prokaryotes, this method has been regarded as the “gold standard” in both microbial phylogeny and ecology studies. However, intragenomic heterogeneity within 16S rRNA genes has been reported in many investigations and is believed to bias the estimation of prokaryotic diversity. In the current study, 2,013 completely sequenced genomes of bacteria and archaea were analyzed and intragenomic heterogeneity was found in 952 genomes (585 species), with 87.5% of the divergence detected being below the 1% level. In particular, some extremophiles (thermophiles and halophiles) were found to harbor highly divergent 16S rRNA genes. Overestimation caused by 16S rRNA gene intragenomic heterogeneity was evaluated at different levels using the full-length and partial 16S rRNA genes usually chosen as targets for pyrosequencing. The result indicates that, at the unique level, full-length 16S rRNA genes can produce an overestimation of as much as 123.7%, while at the 3% level, an overestimation of 12.9% for the V6 region may be introduced. Further analysis showed that intragenomic heterogeneity tends to concentrate in specific positions, with the V1 and V6 regions suffering the most intragenomic heterogeneity and the V4 and V5 regions suffering the least intragenomic heterogeneity in bacteria. This is the most up-to-date overview of the diversity of 16S rRNA genes within prokaryotic genomes. It not only provides general guidance on how much overestimation can be introduced when applying 16S rRNA gene-based methods, due to its intragenomic heterogeneity, but also recommends that, for bacteria, this overestimation be minimized using primers targeting the V4 and V5 regions.     

Advantage of Using Inosine at the 3′ Termini of 16S rRNA Gene Universal Primers for the Study of Microbial Diversity

To overcome the shortcomings of universal 16S rRNA gene primers 8F and 907R when studying the diversity of complex microbial communities, the 3′ termini of both primers were replaced with inosine. A comparison of the clone libraries derived using both primer sets showed seven bacterial phyla amplified by the altered primer set (8F-I/907R-I) whereas the original set amplified sequences belonging almost exclusively to Proteobacteria (95.8%). Sequences belonging to Firmicutes (42.6%) and Thermotogae (9.3%) were more abundant in a library obtained by using 8F-I/907R-I at a PCR annealing temperature of 54°C, while Proteobacteria sequences were more frequent (62.7%) in a library obtained at 50°C, somewhat resembling the result obtained using the original primer set. The increased diversity revealed by using primers 8F-I/907R-I confirms the usefulness of primers with inosine at the 3′ termini in studying the microbial diversity of environmental samples.     

Simultaneous analysis of seven 16S rRNA hypervariable gene regions increases efficiency in marine bacterial diversity detection

Environmental DNA sequencing is the gold standard to reveal microbial community structures. In most applications, a one-fragment PCR approach is applied to amplify a taxonomic marker gene, usually a hypervariable region of the 16S rRNA gene. We used a new reverse complement (RC)-PCR-based assay that amplifies seven out of the nine hypervariable regions of the 16S rRNA gene, to interrogate bacterial communities in sediment samples collected from different coastal marine sites with an impact gradient. In parallel, we employed a traditional one-fragment analysis of the hypervariable V3–V4 region to investigate whether the RC-PCR reveals more of the ‘unseen’ diversity obtained by the one-fragment approach. As a benchmark for the full deck of diversity, we subjected the samples to PCR-free metagenomic sequencing. None of the two PCR-based approaches recorded the full taxonomic repertoire obtained from the metagenomics datasets. However, the RC-PCR approach detected 2.8 times more bacterial genera compared to the near-saturation sequenced V3–V4 samples. RC-PCR is an ideal compromise between the standard one-fragment approach and metagenomics sequencing and may guide future environmental sequencing studies, in which bacterial diversity is a central subject.     

Frequency of formation of chimeric molecules as a consequence of PCR coamplification of 16S rRNA genes from mixed bacterial genomes.

PCR is routinely used in amplification and cloning of rRNA genes from environmental DNA samples for studies of microbial community structure and identification of novel organisms. There have been concerns about generation of chimeric sequences as a consequence of PCR coamplification of highly conserved genes, because such sequences may lead to reports of nonexistent organisms. To quantify the frequency of chimeric molecule formation, mixed genomic DNAs from eight actinomycete species whose 16S rRNA sequences had been determined were used for PCR coamplification of 16S rRNA genes. A large number of cloned 16S ribosomal DNAs were examined by sequence analysis, and chimeric molecules were identified by multiple-sequence alignment with reference species. Here, we report that the level of occurrence of chimeric sequences after 30 cycles of PCR amplification was 32%. We also show that PCR-induced chimeras were formed between different rRNA gene copies from the same organism. Because of the wide use of PCR for direct isolation of 16S rRNA sequences from environmental DNA to assess microbial diversity, the extent of chimeric molecule formation deserves serious attention.     

Succession of Microbial Communities during Hot Composting as Detected by PCR–Single-Strand-Conformation Polymorphism-Based Genetic Profiles of Small-Subunit rRNA Genes

A cultivation-independent technique for genetic profiling of PCR-amplified small-subunit rRNA genes (SSU rDNA) was chosen to characterize the diversity and succession of microbial communities during composting of an organic agricultural substrate. PCR amplifications were performed with DNA directly extracted from compost samples and with primers targeting either (i) the V4–V5 region of eubacterial 16S rRNA genes, (ii) the V3 region in the 16S rRNA genes of actinomycetes, or (iii) the V8–V9 region of fungal 18S rRNA genes. Homologous PCR products were converted to single-stranded DNA molecules by exonuclease digestion and were subsequently electrophoretically separated by their single-strand-conformation polymorphism (SSCP). Genetic profiles obtained by this technique showed a succession and increasing diversity of microbial populations with all primers. A total of 19 single products were isolated from the profiles by PCR reamplification and cloning. DNA sequencing of these molecular isolates showed similarities in the range of 92.3 to 100% to known gram-positive bacteria with a low or high G+C DNA content and to the SSU rDNA of γ-Proteobacteria. The amplified 18S rRNA gene sequences were related to the respective gene regions of Candida krusei and Candida tropicalis. Specific molecular isolates could be attributed to different composting stages. The diversity of cultivated bacteria isolated from samples taken at the end of the composting process was low. A total of 290 isolates were related to only 6 different species. Two or three of these species were also detectable in the SSCP community profiles. Our study indicates that community SSCP profiles can be highly useful for the monitoring of bacterial diversity and community successions in a biotechnologically relevant process.     

Discovery of the Novel Candidate Phylum “Poribacteria” in Marine Sponges

Marine sponges (Porifera) harbor large amounts of commensal microbial communities within the sponge mesohyl. We employed 16S rRNA gene library construction using specific PCR primers to provide insights into the phylogenetic identity of an abundant sponge-associated bacterium that is morphologically characterized by the presence of a membrane-bound nucleoid. In this study, we report the presence of a previously unrecognized evolutionary lineage branching deeply in the domain Bacteria that is moderately related to the Planctomycetes, Verrucomicrobia, and Chlamydia lines of decent. Because members of this lineage showed <75% 16S rRNA gene sequence similarity to known bacterial phyla, we suggest the status of a new candidate phylum, named “Poribacteria”, to acknowledge the affiliation of the new bacterium with sponges. The affiliation of the morphologically conspicuous sponge bacterium with the novel phylogenetic lineage was confirmed by fluorescence in situ hybridization with newly designed probes targeting different sites of the poribacterial 16S rRNA. Consistent with electron microscopic observations of cell compartmentalization, the fluorescence signals appeared in a ring-shaped manner. PCR screening with “Poribacteria”-specific primers gave positive results for several other sponge species, while samples taken from the environment (seawater, sediments, and a filter-feeding tunicate) were PCR negative. In addition to a report for Planctomycetes, this is the second report of cell compartmentalization, a feature that was considered exclusive to the eukaryotic domain, in prokaryotes.     

Detection of Actinobacteria cultivated from environmental samples reveals bias in universal primers

AIMS: The aims of this study were to develop media to cultivate actinomycetes, screen the resulting isolates with Actinobacteria-specific primers, and examine the efficacy of detection of the actinobacterial isolates with universal primers. METHODS AND RESULTS: Soil-extract medium was developed for a terrestrial bluff environment. Recovered isolates were subjected to polymerase chain reaction (PCR) with taxon-specific primers to identify Actinobacteria. Universal bacterial primers 24f and 1492r (modified and original versions) were used to amplify the 16S rRNA gene from the putative Actinobacteria. While both reverse primers failed to provide amplification products from 20% to 50% of the isolates, the 1492r primer detected Actinobacteria more effectively than 1492r-mod. The region of the gene containing the annealing site for the 1492r primers from 15 isolates that failed to amplify showed no differences in nucleotide sequence to the original 1492r primer. CONCLUSIONS: Universal 16S rRNA gene primers are not capable of amplifying this gene from all bacteria within an environmental sample. Some Actinobacteria may share 100% sequence similarity to universal primers but remain undetected. SIGNIFICANCE AND IMPACT OF THE STUDY: These findings are important for studies of particular taxa in environmental samples where reactions utilizing universal primers may not reveal the extent of their presence and diversity.     

Competitive PCR-DGGE analysis of bacterial mixtures: an internal standard and an appraisal of template enumeration accuracy

Analysis of polymerase chain reaction (PCR) amplified 16S rDNA fragments from environmental samples by denaturing gradients of chemicals or heat [denaturing gradient gel electrophoresis (DGGE) and thermal gradient gel electrophoresis (TGGE)] within polyacrylamide gels is a popular tool in microbial ecology. Difficulties in acceptance of the technique and interpretation of the results remain, due to its qualitative nature. In this study we have addressed this problem by the construction and evaluation of a quantitative standard for incorporation into test DNA samples. The standard was based on a naturally occurring 16S rRNA gene carried by the X-endosymbiont of the psyllid Anomoneura mori, a gamma-proteobacterium. This sequence is the most AT-rich 16S rDNA gene recovered from any cultured organism or environmental sample described to date, and a specifically amplified rDNA fragment denatured under exceptionally low stringency denaturing conditions. The native sequence was modified to incorporate perfect matches to the PCR primers used. The efficiency of amplification of this standard in comparison to a range of 16S rDNA sequences and the errors involved in enumerating template molecules under a range of PCR conditions are demonstrated and quantified. Tests indicated that highly accurate counts of released target molecules from a range of bacterial cells could be achieved in both laboratory mixtures and compost.     

Next-generation Sequencing of 16S Ribosomal RNA Gene Amplicons

One of the major questions in microbial ecology is “who is there?” This question can be answered using various tools, but one of the long-lasting gold standards is to sequence 16S ribosomal RNA (rRNA) gene amplicons generated by domain-level PCR reactions amplifying from genomic DNA. Traditionally, this was performed by cloning and Sanger (capillary electrophoresis) sequencing of PCR amplicons. The advent of next-generation sequencing has tremendously simplified and increased the sequencing depth for 16S rRNA gene sequencing. The introduction of benchtop sequencers now allows small labs to perform their 16S rRNA sequencing in-house in a matter of days. Here, an approach for 16S rRNA gene amplicon sequencing using a benchtop next-generation sequencer is detailed. The environmental DNA is first amplified by PCR using primers that contain sequencing adapters and barcodes. They are then coupled to spherical particles via emulsion PCR. The particles are loaded on a disposable chip and the chip is inserted in the sequencing machine after which the sequencing is performed. The sequences are retrieved in fastq format, filtered and the barcodes are used to establish the sample membership of the reads. The filtered and binned reads are then further analyzed using publically available tools. An example analysis where the reads were classified with a taxonomy-finding algorithm within the software package Mothur is given. The method outlined here is simple, inexpensive and straightforward and should help smaller labs to take advantage from the ongoing genomic revolution.     

Molecular Evidence for the Broad Distribution of Anaerobic Ammonium-Oxidizing Bacteria in Freshwater and Marine Sediments

Previously available primer sets for detecting anaerobic ammonium-oxidizing (anammox) bacteria are inefficient, resulting in a very limited database of such sequences, which limits knowledge of their ecology. To overcome this limitation, we designed a new primer set that was 100% specific in the recovery of ~700-bp 16S rRNA gene sequences with >96% homology to the “Candidatus Scalindua” group of anammox bacteria, and we detected this group at all sites studied, including a variety of freshwater and marine sediments and permafrost soil. A second primer set was designed that exhibited greater efficiency than previous primers in recovering full-length (1,380-bp) sequences related to “Ca. Scalindua,” “Candidatus Brocadia,” and “Candidatus Kuenenia.” This study provides evidence for the widespread distribution of anammox bacteria in that it detected closely related anammox 16S rRNA gene sequences in 11 geographically and biogeochemically diverse freshwater and marine sediments.