RiboFR-Seq: a novel approach to linking 16S rRNA amplicon profiles to metagenomes

16S rRNA amplicon analysis and shotgun metagenome sequencing are two main culture-independent strategies to explore the genetic landscape of various microbial communities. Recently, numerous studies have employed these two approaches together, but downstream data analyses were performed separately, which always generated incongruent or conflict signals on both taxonomic and functional classifications. Here we propose a novel approach, RiboFR-Seq (Ribosomal RNA gene flanking region sequencing), for capturing both ribosomal RNA variable regions and their flanking protein-coding genes simultaneously. Through extensive testing on clonal bacterial strain, salivary microbiome and bacterial epibionts of marine kelp, we demonstrated that RiboFR-Seq could detect the vast majority of bacteria not only in well-studied microbiomes but also in novel communities with limited reference genomes. Combined with classical amplicon sequencing and shotgun metagenome sequencing, RiboFR-Seq can link the annotations of 16S rRNA and metagenomic contigs to make a consensus classification. By recognizing almost all 16S rRNA copies, the RiboFR-seq approach can effectively reduce the taxonomic abundance bias resulted from 16S rRNA copy number variation. We believe that RiboFR-Seq, which provides an integrated view of 16S rRNA profiles and metagenomes, will help us better understand diverse microbial communities.     

Evaluating Bias of Illumina-Based Bacterial 16S rRNA Gene Profiles

Massively parallel sequencing of 16S rRNA genes enables the comparison of terrestrial, aquatic, and host-associated microbial communities with sufficient sequencing depth for robust assessments of both alpha and beta diversity. Establishing standardized protocols for the analysis of microbial communities is dependent on increasing the reproducibility of PCR-based molecular surveys by minimizing sources of methodological bias. In this study, we tested the effects of template concentration, pooling of PCR amplicons, and sample preparation/interlane sequencing on the reproducibility associated with paired-end Illumina sequencing of bacterial 16S rRNA genes. Using DNA extracts from soil and fecal samples as templates, we sequenced pooled amplicons and individual reactions for both high (5- to 10-ng) and low (0.1-ng) template concentrations. In addition, all experimental manipulations were repeated on two separate days and sequenced on two different Illumina MiSeq lanes. Although within-sample sequence profiles were highly consistent, template concentration had a significant impact on sample profile variability for most samples. Pooling of multiple PCR amplicons, sample preparation, and interlane variability did not influence sample sequence data significantly. This systematic analysis underlines the importance of optimizing template concentration in order to minimize variability in microbial-community surveys and indicates that the practice of pooling multiple PCR amplicons prior to sequencing contributes proportionally less to reducing bias in 16S rRNA gene surveys with next-generation sequencing.     

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.     

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.     

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.     

Analysis of microbial communities using culture-dependent and culture-independent approaches in an anaerobic/aerobic SBR reactor

Comparative analysis of microbial communities in a sequencing batch reactor which performed enhanced biological phosphorus removal (EBPR) was carried out using a cultivation-based technique and 16S rRNA gene clone libraries. A standard PCR protocol and a modified PCR protocol with low PCR cycle was applied to the two clone libraries of the 16S rRNA gene sequences obtained from EBPR sludge, respectively, and the resulting 424 clones were analyzed using restriction fragment length polymorphisms (RFLPs) on 16S rRNA gene inserts. Comparison of two clone libraries showed that the modified PCR protocol decreased the incidence of distinct fragment patterns from about 63% (137 of 217) in the standard PCR method to about 34% (70 of 207) under the modified protocol, suggesting that just a low level of PCR cycling (5 cycles after 15 cycles) can significantly reduce the formation of chimeric DNA in the final PCR products. Phylogenetic analysis of 81 groups with distinct RFLP patterns that were obtained using the modified PCR method revealed that the clones were affiliated with at least 11 phyla or classes of the domain Bacteria. However, the analyses of 327 colonies, which were grouped into just 41 distinct types by RFLP analysis, showed that they could be classified into five major bacterial lineages: alpha, beta, gamma- Proteobacteria, Actinobacteria, and the phylum Bacteroidetes, which indicated that the microbial community yielded from the cultivation-based method was still much simpler than that yielded from the PCR-based molecular method. In this study, the discrepancy observed between the communities obtained from PCR-based and cultivation-based methods seems to result from low culturabilities of bacteria or PCR bias even though modified culture and PCR methods were used. Therefore, continuous development of PCR protocol and cultivation techniques is needed to reduce this discrepancy.     

A novel method for determining microflora composition using dynamic phylogenetic analysis of 16S ribosomal RNA deep sequencing data

Deep sequencing of the 16S rRNA gene provides a comprehensive view of bacterial communities in a particular environment and has expanded our ability to study the impact of the microflora on human health and disease. Current analysis methods rely on comparisons of the sequences generated with an expanding but limited set of annotated 16S rRNA sequences or phylogenic clustering of sequences based on arbitrary similarity cutoffs. We describe a novel approach to characterize bacterial composition using deep sequencing of 16S rRNA gene. Our method defines operational taxonomic units based on phylogenetic tree reconstruction and dynamic clustering of sequences using solely sequencing data. These OTUs can be used to identify differences in bacteria abundance between environments. This approach can perform better than previous phylogenetic methods and will significantly improve our understanding of the microfloral role on human diseases by providing a comprehensive analysis of the microbial composition from various bacterial communities.     

Comparison of PCR primer-based strategies for characterization of ammonia oxidizer communities in environmental samples

PCR-based techniques are commonly used to characterize microbial communities, but are subject to bias that is difficult to assess. This study aimed to evaluate bias of several PCR primer-based strategies used to study diversity of autotrophic ammonia oxidizers. 16S rRNA genes from soil- or sediment-DNA were amplified using primers considered either selective or specific for betaproteobacterial ammonia oxidizers. Five approaches were assessed: (a) amplification with primers betaAMO143f-betaAMO1315r; (b) amplification with primers CTO189f-CTO654r; (c) nested amplification with betaAMO143f-betaAMO1315r followed by CTO189f-CTO654r primers; (d) nested amplification with betaAMO143f-betaAMO1315r and CTO189f-Pf1053r primers; (e) nested amplification with 27f-1492r and CTO189f-CTO654r primers. Amplification products were characterized by denaturing gradient gel electrophoresis (DGGE) analysis after further amplification with 357f-GC-518r primers. DGGE profiles of soil communities were heterogeneous and depended on the approach followed. Ammonia oxidizer diversity was higher using approaches (b), (c) and (e) than using (a) and (d), where sequences of the most prominent bands showed similarities to nonammonia oxidizers. Profiles from marine sediments were more consistent, regardless of the approach adopted, and sequence analysis of excised bands indicated that these consisted of ammonia oxidizers only. The study demonstrates the importance of choice of primer, of screening for sequences of nontarget organisms and use of several approaches when characterizing microbial communities in natural environments.     

3 次连续重复提取DNA 能较好反映土壤微生物丰度

【目的】研究同一个土壤需要反复提取几次才能在最大程度上反映土壤微生物的丰度,探讨风干土壤代替新鲜土壤用于微生物丰度研究的可行性。【方法】针对两种理化性质具有较大差异的旱地和稻田新鲜土壤及其风干土壤,分别对土壤微生物进行5次连续裂解提取DNA。通过实时荧光定量PCR技术分析连续反复提取对土壤古菌和细菌16S rRNA gene数量、氨氧化古菌和细菌功能基因amoA数量的影响。【结果】3次连续提取DNA占5次提取DNA总量的76%以上,氨氧化古菌、氨氧化细菌、古菌和细菌4类微生物的3次连续提取最低回收率为77.5%;与新鲜土壤相比,风干处理导致氨氧化古菌、氨氧化细菌、古菌、细菌的数量分别降低84.3%、81.2%、12.5%和90.3%,然而,2种土壤风干过程中主要微生物类群的数量变化规律基本一致,表明土壤微生物对风干处理的响应可能受土壤类型的影响较小。【结论】土壤微生物连续3次裂解能较好反映微生物丰度。与新鲜土壤相比,风干过程显著降低了土壤微生物丰度,然而,通过风干土壤中微生物丰度的变化趋势反映新鲜土壤中微生物数量变化规律具有一定的可行性。     

Evidence from GC-TRFLP that bacterial communities in soil are lognormally distributed

The Species Abundance Distribution (SAD) is a fundamental property of ecological communities and the form and formation of SADs have been examined for a wide range of communities including those of microorganisms. Progress in understanding microbial SADs, however, has been limited by the remarkable diversity and vast size of microbial communities. As a result, few microbial systems have been sampled with sufficient depth to generate reliable estimates of the community SAD. We have used a novel approach to characterize the SAD of bacterial communities by coupling genomic DNA fractionation with analysis of terminal restriction fragment length polymorphisms (GC-TRFLP). Examination of a soil microbial community through GC-TRFLP revealed 731 bacterial operational taxonomic units (OTUs) that followed a lognormal distribution. To recover the same 731 OTUs through analysis of DNA sequence data is estimated to require analysis of 86,264 16S rRNA sequences. The approach is examined and validated through construction and analysis of simulated microbial communities in silico. Additional simulations performed to assess the potential effects of PCR bias show that biased amplification can cause a community whose distribution follows a power-law function to appear lognormally distributed. We also show that TRFLP analysis, in contrast to GC-TRFLP, is not able to effectively distinguish between competing SAD models. Our analysis supports use of the lognormal as the null distribution for studying the SAD of bacterial communities as for plant and animal communities.     

Species Identification and Profiling of Complex Microbial Communities Using Shotgun Illumina Sequencing of 16S rRNA Amplicon Sequences

The high throughput and cost-effectiveness afforded by short-read sequencing technologies, in principle, enable researchers to perform 16S rRNA profiling of complex microbial communities at unprecedented depth and resolution. Existing Illumina sequencing protocols are, however, limited by the fraction of the 16S rRNA gene that is interrogated and therefore limit the resolution and quality of the profiling. To address this, we present the design of a novel protocol for shotgun Illumina sequencing of the bacterial 16S rRNA gene, optimized to amplify more than 90% of sequences in the Greengenes database and with the ability to distinguish nearly twice as many species-level OTUs compared to existing protocols. Using several in silico and experimental datasets, we demonstrate that despite the presence of multiple variable and conserved regions, the resulting shotgun sequences can be used to accurately quantify the constituents of complex microbial communities. The reconstruction of a significant fraction of the 16S rRNA gene also enabled high precision (>90%) in species-level identification thereby opening up potential application of this approach for clinical microbial characterization.     

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 gamma-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.     

Quantification of mcrA by fluorescent PCR in methanogenic and methanotrophic microbial communities

A quantitative fluorogenic PCR method for detecting methanogenic and methanotrophic orders was established using a refined primer set for the methyl coenzyme M reductase subunit A gene (mcrA). The method developed was applied to several microbial communities in which diversity and abundance of methanogens or anaerobic methanotrophs (ANMEs) was identified by 16S rRNA gene clone analysis, and strong correlations between the copy numbers of mcrA with those of archaeal 16S rRNA genes in the communities were observed. The assay can be applied to detecting and assessing the abundance of methanogens and/or ANMEs in anoxic environments that could not be detected by 16S rRNA gene sequence analyses.     

Novel approach to quantitative detection of specific rRNA in a microbial community, using catalytic DNA

We developed a novel method for the quantitative detection of the 16S rRNA of a specific bacterial species in the microbial community by using deoxyribozyme (DNAzyme), which possesses the catalytic function to cleave RNA in a sequence-specific manner. A mixture of heterogeneous 16S rRNA containing the target 16S rRNA was incubated with a species-specific DNAzyme. The cleaved target 16S rRNA was separated from the intact 16S rRNA by electrophoresis, and then their amounts were compared for the quantitative detection of target 16S rRNA. This method was used to determine the abundance of the 16S rRNA of a filamentous bacterium, Sphaerotilus natans, in activated sludge, which is a microbial mixture used in wastewater treatment systems. The result indicated that this DNAzyme-based approach would be applicable to actual microbial communities.     

Novel Approach to Quantitative Detection of Specific rRNA in a Microbial Community, Using Catalytic DNA

We developed a novel method for the quantitative detection of the 16S rRNA of a specific bacterial species in the microbial community by using deoxyribozyme (DNAzyme), which possesses the catalytic function to cleave RNA in a sequence-specific manner. A mixture of heterogeneous 16S rRNA containing the target 16S rRNA was incubated with a species-specific DNAzyme. The cleaved target 16S rRNA was separated from the intact 16S rRNA by electrophoresis, and then their amounts were compared for the quantitative detection of target 16S rRNA. This method was used to determine the abundance of the 16S rRNA of a filamentous bacterium, Sphaerotilus natans, in activated sludge, which is a microbial mixture used in wastewater treatment systems. The result indicated that this DNAzyme-based approach would be applicable to actual microbial communities.     

Selective Recovery of 16S rRNA Sequences from Natural Microbial Communities in the Form of cDNA

Cloning of cDNA obtained from 16S rRNA (16S rcDNA) selectively retrieves species-specific sequence information useful for analyzing the composition and structure of natural microbial communities. With this technique we obtained recombinant 16S rcDNA libraries from Escherichia coli and from a model hot-spring cyanobacterial-mat community. The recombinant plasmids contained exclusively 16S rRNA-derived inserts. This selective approach is independent of biasing culture techniques and eliminates the laborious screening required to locate 16S rRNA gene-bearing recombinants in genomic DNA libraries obtained from natural communities.     

Mathematical modeling of 16S ribosomal DNA amplification reveals optimal conditions for the interrogation of complex microbial communities with phylogenetic microarrays

MOTIVATION: Many current studies of complex microbial communities rely on the isolation of community genomic DNA, amplification of 16S ribosomal RNA genes (rDNA) and subsequent examination of community structure through interrogation of the amplified 16S rDNA pool by high-throughput sequencing, phylogenetic microarrays or quantitative PCR. RESULTS: Here we describe the development of a mathematical model aimed to simulate multitemplate amplification of 16S ribosomal DNA sample and subsequent detection of these amplified 16S rDNA species by phylogenetic microarray. Using parameters estimated from the experimental results obtained in the analysis of intestinal microbial communities with Microbiota Array, we show that both species detection and the accuracy of species abundance estimates depended heavily on the number of PCR cycles used to amplify 16S rDNA. Both parameters initially improved with each additional PCR cycle and reached optimum between 15 and 20 cycles of amplification. The use of more than 20 cycles of PCR amplification and/or more than 50 ng of starting genomic DNA template was, however, detrimental to both the fraction of detected community members and the accuracy of abundance estimates. Overall, the outcomes of the model simulations matched well available experimental data. Our simulations also showed that species detection and the accuracy of abundance measurements correlated positively with the higher sample-wide PCR amplification rate, lower template-to-template PCR bias and lower number of species in the interrogated community. The developed model can be easily modified to simulate other multitemplate DNA mixtures as well as other microarray designs and PCR amplification protocols.     

Comprehensive single-PCR 16S and 18S rRNA community analysis validated with mock communities,and estimation of sequencing bias against 18S

Universal primers for SSU rRNA genes allow profiling of natural communities by simultaneously amplifying templates from Bacteria, Archaea, and Eukaryota in a single PCR reaction. Despite the potential to show relative abundance for all rRNA genes, universal primers are rarely used, due to various concerns including amplicon length variation and its effect on bioinformatic pipelines. We thus developed 16S and 18S rRNA mock communities and a bioinformatic pipeline to validate this approach. Using these mocks, we show that universal primers (515Y/926R) outperformed eukaryote-specific V4 primers in observed versus expected abundance correlations (slope = 0.88 vs. 0.67–0.79), and mock community members with single mismatches to the primer were strongly underestimated (threefold to eightfold). Using field samples, both primers yielded similar 18S beta-diversity patterns (Mantel test, p < 0.001) but differences in relative proportions of many rarer taxa. To test for length biases, we mixed mock communities (16S + 18S) before PCR and found a twofold underestimation of 18S sequences due to sequencing bias. Correcting for the twofold underestimation, we estimate that, in Southern California field samples (1.2–80 μm), there were averages of 35% 18S, 28% chloroplast 16S, and 37% prokaryote 16S rRNA genes. These data demonstrate the potential for universal primers to generate comprehensive microbiome profiles.     

Two-step denaturing gradient gel electrophoresis (2S-DGGE), a gel-based strategy to capture full-length 16S rRNA gene sequences

Obtaining full-length 16S rRNA gene sequences is important for generating accurate taxonomy assignments of bacteria, which normally is realized via clone library construction. However, the application of clone library has been hindered due to its limitations in sample throughput and in capturing minor populations (<1?% of total microorganisms). To overcome these limitations, a new strategy, two-step denaturing gradient gel electrophoresis (2S-DGGE), is developed to obtain full-length 16S rRNA gene sequences. 2S-DGGE can compare microbial communities based on its first-round DGGE profiles and generate partial 16S rRNA gene sequences (8-534?bp, Escherichia coli numbering). Then, strain-specific primers can be designed based on sequence information of bacteria of interest to PCR amplify their remaining 16S rRNA gene sequences (515-1541?bps, E. coli numbering). The second-round DGGE can confirm DNA sequence purity of these PCR products. Finally, the full-length 16S rRNA gene sequences can be obtained through combining the two partial DNA sequences. By employing 2S-DGGE, taxonomies of a group of dehalogenating bacteria have been assigned based on their full-length 16S rRNA gene sequences, several of which existed in dehalogenating microcosms as minor populations. In all, 2S-DGGE can be utilized as a medium throughput method for taxonomic identification of interested/minor populations from single or multiple microbial consortia.     

Microbial communities associated with geological horizons in coastal subseafloor sediments from the sea of okhotsk

Microbial communities from a subseafloor sediment core from the southwestern Sea of Okhotsk were evaluated by performing both cultivation-dependent and cultivation-independent (molecular) analyses. The core, which extended 58.1 m below the seafloor, was composed of pelagic clays with several volcanic ash layers containing fine pumice grains. Direct cell counting and quantitative PCR analysis of archaeal and bacterial 16S rRNA gene fragments indicated that the bacterial populations in the ash layers were approximately 2 to 10 times larger than those in the clays. Partial sequences of 1,210 rRNA gene clones revealed that there were qualitative differences in the microbial communities from the two different types of layers. Two phylogenetically distinct archaeal assemblages in the Crenarchaeota, the miscellaneous crenarchaeotic group and the deep-sea archaeal group, were the most predominant archaeal 16S rRNA gene components in the ash layers and the pelagic clays, respectively. Clones of 16S rRNA gene sequences from members of the gamma subclass of the class Proteobacteria dominated the ash layers, whereas sequences from members of the candidate division OP9 and the green nonsulfur bacteria dominated the pelagic clay environments. Molecular (16S rRNA gene sequence) analysis of 181 isolated colonies revealed that there was regional proliferation of viable heterotrophic mesophiles in the volcanic ash layers, along with some gram-positive bacteria and actinobacteria. The porous ash layers, which ranged in age from tens of thousands of years to hundreds of thousands of years, thus appear to be discrete microbial habitats within the coastal subseafloor clay sediment, which are capable of harboring microbial communities that are very distinct from the communities in the more abundant pelagic clays.