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

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

Substitution by Inosine at the 3��-Ultimate and Penultimate Positions of 16S rRNA Gene Universal Primers

Universal 16S rRNA gene primers (8F and 518R) bearing inosine substitutions at either the 3??-ultimate or the 3??-ultimate and penultimate base positions were exploited for the first time to study the bacterial community associated with coral polymicrobial Black Band Disease (BBD). Inosine-modified universal primer pairs display some shifting in the composition of 16S rRNA gene libraries, as well as expanding the observed diversity of a BBD bacterial community at the family/class level. Possible explanations for the observed shifts are discussed. These results thus point to the need for adopting multiple approaches in designing 16S rRNA universal primers for PCR amplification and subsequent construction of 16S rRNA gene libraries or pyrosequencing in the exploration of complex microbial communities.     

Biased Diversity Metrics Revealed by Bacterial 16S Pyrotags Derived from Different Primer Sets

In recent years, PCR-based pyrosequencing of 16S rRNA genes has continuously increased our understanding of complex microbial communities in various environments of the Earth. However, there is always concern on the potential biases of diversity determination using different 16S rRNA gene primer sets and covered regions. Here, we first report how bacterial 16S rRNA gene pyrotags derived from a series of different primer sets resulted in the biased diversity metrics. In total, 14 types of pyrotags were obtained from two-end pyrosequencing of 7 amplicon pools generated by 7 primer sets paired by 1 of 4 forward primers (V1F, V3F, V5F, and V7F) and 1 of 4 reverse primers (V2R, V4R, V6R, and V9R), respectively. The results revealed that: i) the activated sludge exhibited a large bacterial diversity that represented a broad range of bacterial populations and served as a good sample in this methodology research; ii) diversity metrics highly depended on the selected primer sets and covered regions; iii) paired pyrotags obtained from two-end pyrosequencing of each short amplicon displayed different diversity metrics; iv) relative abundance analysis indicated the sequencing depth affected the determination of rare bacteria but not abundant bacteria; v) the primer set of V1F and V2R significantly underestimated the diversity of activated sludge; and vi) the primer set of V3F and V4R was highly recommended for future studies due to its advantages over other primer sets. All of these findings highlight the significance of this methodology research and offer a valuable reference for peer researchers working on microbial diversity determination.     

Investigation of bacterial and archaeal communities: novel protocols using modern sequencing by Illumina MiSeq and traditional DGGE-cloning

Different protocols based on Illumina high-throughput DNA sequencing and denaturing gradient gel electrophoresis (DGGE)-cloning were developed and applied for investigating hot spring related samples. The study was focused on three target genes: archaeal and bacterial 16S rRNA and mcrA of methanogenic microflora. Shorter read lengths of the currently most popular technology of sequencing by Illumina do not allow analysis of the complete 16S rRNA region, or of longer gene fragments, as was the case of Sanger sequencing. Here, we demonstrate that there is no need for special indexed or tailed primer sets dedicated to short variable regions of 16S rRNA since the presented approach allows the analysis of complete bacterial 16S rRNA amplicons (V1–V9) and longer archaeal 16S rRNA and mcrA sequences. Sample augmented with transposon is represented by a set of approximately 300 bp long fragments that can be easily sequenced by Illumina MiSeq. Furthermore, a low proportion of chimeric sequences was observed. DGGE-cloning based strategies were performed combining semi-nested PCR, DGGE and clone library construction. Comparing both investigation methods, a certain degree of complementarity was observed confirming that the DGGE-cloning approach is not obsolete. Novel protocols were created for several types of laboratories, utilizing the traditional DGGE technique or using the most modern Illumina sequencing.     

Coamplification of Eukaryotic DNA with 16S rRNA Gene-Based PCR Primers: Possible Consequences for Population Fingerprinting of Complex Microbial Communities

The main aim of this study was to evaluate the specificity of three commonly used 16S rRNA gene-based polymerase chain reaction (PCR) primer sets for bacterial community analysis of samples contaminated with eukaryotic DNA. The specificity of primer sets targeting the V₃, V₃-V₅, and V₆-V₈ hypervariable regions of the 16S rRNA gene was investigated in silico and by community fingerprinting of human and fish intestinal samples. Both in silico and PCR-based analysis revealed cross-reactivity of the V₃ and V₃-V₅ primers with the 18S rRNA gene of human and sturgeon. The consequences of this primer anomaly were illustrated by denaturing gradient gel electrophoresis (DGGE) profiling of microbial communities in human feces and mixed gut of Siberian sturgeon. DGGE profiling indicated that the cross-reactivity of 16S rRNA gene primers with nontarget eukaryotic DNA might lead to an overestimation of bacterial biodiversity. This study has confirmed previous sporadic indications in literature indicating that several commonly applied 16S rRNA gene primer sets lack specificity toward bacteria in the presence of eukaryotic DNA. The phenomenon of cross-reactivity is a potential source of systematic error in all biodiversity studies where no subsequent analysis of individual community amplicons by cloning and sequencing is performed.     

Designing specific bacterial 16S primers to sequence and quantitate plant endo-bacteriome

Plant endophytic bacteria colonize the internal tissues of plants and interact with plants closely. The past two decades have witnessed the increasing application of next-generation 16S rRNA gene sequencing in the investigation of bacterial communities. However, deciphering plant endo-bacterial communities by this method is difficult because of the co-amplification of massive plant organellar DNAs with bacterial 16S. Here, we designed polymerase chain reaction (PCR) primer sets, including 799F/1107R, 322F/796R, and 322F-Dr/796Rs (primer pair 322F/796R with a penultimate-base substitution in 322F), that can specifically amplify bacterial 16S from plant total DNAs. We computationally and experimentally evaluated the specificity, coverage, and accuracy of the newly designed primer sets. Both 799F/1107R and 322F-Dr/796Rs produced plant DNA-free 16S amplicon libraries or reduced plant DNA contamination to lower than 5% for the plant materials with extremely-low-abundance bacterial communities. The primer set 322F-A/796R was used through absolute quantitative PCR to quantitate the population size of rice leaf or root endo-bacteriome, which revealed 10⁶–10⁷ and 10⁹–10¹⁰ bacteria per gram fresh weight, respectively. These 16S primer sets and amplification methods enable the simple and inexpensive next-generation sequencing and quantification of plant endo-bacteriome, which will significantly advance studies on the plant-related microbiome.

     

Target region selection is a critical determinant of community fingerprints generated by 16S pyrosequencing

Pyrosequencing of 16S rRNA genes allows for in-depth characterization of complex microbial communities. Although it is known that primer selection can influence the profile of a community generated by sequencing, the extent and severity of this bias on deep-sequencing methodologies is not well elucidated. We tested the hypothesis that the hypervariable region targeted for sequencing and primer degeneracy play important roles in influencing the composition of 16S pyrotag communities. Subgingival plaque from deep sites of current smokers with chronic periodontitis was analyzed using Sanger sequencing and pyrosequencing using 4 primer pairs. Greater numbers of species were detected by pyrosequencing than by Sanger sequencing. Rare taxa constituted nearly 6% of each pyrotag community and less than 1% of the Sanger sequencing community. However, the different target regions selected for pyrosequencing did not demonstrate a significant difference in the number of rare and abundant taxa detected. The genera Prevotella, Fusobacterium, Streptococcus, Granulicatella, Bacteroides, Porphyromonas and Treponema were abundant when the V1-V3 region was targeted, while Streptococcus, Treponema, Prevotella, Eubacterium, Porphyromonas, Campylobacter and Enterococcus predominated in the community generated by V4-V6 primers, and the most numerous genera in the V7-V9 community were Veillonella, Streptococcus, Eubacterium, Enterococcus, Treponema, Catonella and Selenomonas. Targeting the V4-V6 region failed to detect the genus Fusobacterium, while the taxa Selenomonas, TM7 and Mycoplasma were not detected by the V7-V9 primer pairs. The communities generated by degenerate and non-degenerate primers did not demonstrate significant differences. Averaging the community fingerprints generated by V1-V3 and V7-V9 primers provided results similar to Sanger sequencing, while allowing a significantly greater depth of coverage than is possible with Sanger sequencing. It is therefore important to use primers targeted to these two regions of the 16S rRNA gene in all deep-sequencing efforts to obtain representational characterization of complex microbial communities.     

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.     

Effects of Experimental Choices and Analysis Noise on Surveys of the “Rare Biosphere”

When planning a survey of 16S rRNA genes from a complex environment, investigators face many choices including which primers to use and how to taxonomically classify sequences. In this study, we explored how these choices affected a survey of microbial diversity in a sample taken from the aerobic basin of the activated sludge of a North Carolina wastewater treatment plant. We performed pyrosequencing reactions on PCR products generated from primers targeting the V1-V2, V6, and V6-V7 variable regions of the 16S rRNA gene. We compared these sequences to 16S rRNA gene sequences found in a whole-genome shotgun pyrosequencing run performed on the same sample. We found that sequences generated from primers targeting the V1-V2 variable region had the best match to the whole-genome shotgun reaction across a range of taxonomic classifications from phylum to family. Pronounced differences between primer sets, however, occurred in the “rare biosphere” involving taxa that we observed in fewer than 11 sequences. We also examined the results of analysis strategies comparing a classification scheme using a nearest-neighbor approach to directly classifying sequences with a naïve Bayesian algorithm. Again, we observed pronounced differences between these analysis schemes in infrequently observed taxa. We conclude that if a study is meant to probe the rare biosphere, both the experimental conditions and analysis choices will have a profound impact on the observed results.For nearly 3 decades, investigations of the distribution of microbes in complex environments have focused on the use of rRNA genes (1, 2, 4, 11, 16, 18, 19, 22, 24). Because the full-length 16S rRNA sequence can be obtained with paired-end reads via traditional Sanger sequencing, until recently most studies of the 16S rRNA gene captured most or nearly most of the 16S sequence length. New pyrosequencing technologies, however, have recently been introduced that greatly reduce the per base cost of sequencing but with shorter read lengths than traditional Sanger sequencing (17). This new approach has proven powerful, yielding a previously unobtainable view of rare taxa (7, 12-14, 25).The shorter reads produced by pyrosequencing require the choice of a particular region of the 16S rRNA gene to target for pyrosequencing as well as the choice of an algorithm to classify the taxonomy of the shorter reads. In their initial surveys of microbial diversity with pyrosequencing (12, 14, 25), Sogin and colleagues targeted the V6 variable region, in part because it is was small enough to be captured with the 100-bp reads of the pyrosequencing technology available at the time. Recently, the read length of 454 pyrosequencing machines has been increased to an average of ∼250 bp. This allows for more flexibility in primer design and opens up the possibility of targeting regions of the 16S rRNA gene other than V6. In recent work, Huse et al. took advantage of this new capability to compare the classifications made for the human gut microbiome with the V6 and longer V3 regions (13). Plotting the taxonomic abundance of these two sequence sets against each other yielded an excellent correlation (r² = 0.99), suggesting that the choice of which variable region to target makes little difference. In this report, we introduce a data set examining the performance of sets of primers targeting the V1-V2, V6, and V6-V7 regions. By using a sample for which we have also generated a whole-genome shotgun sequencing run with 250 bp reads, we were able to compare the observed 16S rRNA genes in samples with and without an initial PCR step targeting the 16S rRNA gene. Our results demonstrate that experimental choices such as which region of the 16S rRNA gene to sequence and which algorithm to use to classify taxa are much more likely to affect observations of the “rare biosphere” than more commonly observed taxa.     

Detection of bacterial and mycoplasma contamination in cell cultures by polymerase chain reaction

A fast and simple method to detect bacterial and especially mycoplasma contamination in tissue culture by means of polymerase chain reaction (PCR) amplification is described. In a first step the universal primer pairs P1/P2 (190-bp fragment) and P3/P4 (120-bp fragment) directed to different conserved parts of the prokaryotic 16S rRNA gene are used. A positive signal after amplification on cell culture DNA with these primers provides an indication of bacterial infection. Using the internal primers IP1, IP3 and IP'3 complementary to a part of the V4 and V8 variable regions of the 16S rRNA gene, in combination with a universal primer, cultures contaminated with mycoplasma could be identified. Six mycoplasma species, typical contaminants in tissue cultures, were investigated: Mycoplasma orale, M. fermentans, M. arginini, M. hyorhinis, M. hominis and Aeromonas laidlawii. This mycoplasma test is an easy, specific and sensitive assay which should be extremely useful in any tissue culture setting.     

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.     

Analysis,Optimization and Verification of Illumina-Generated 16S rRNA Gene Amplicon Surveys

The exploration of microbial communities by sequencing 16S rRNA genes has expanded with low-cost, high-throughput sequencing instruments. Illumina-based 16S rRNA gene sequencing has recently gained popularity over 454 pyrosequencing due to its lower costs, higher accuracy and greater throughput. Although recent reports suggest that Illumina and 454 pyrosequencing provide similar beta diversity measures, it remains to be demonstrated that pre-existing 454 pyrosequencing workflows can transfer directly from 454 to Illumina MiSeq sequencing by simply changing the sequencing adapters of the primers. In this study, we modified 454 pyrosequencing primers targeting the V4-V5 hyper-variable regions of the 16S rRNA gene to be compatible with Illumina sequencers. Microbial communities from cows, humans, leeches, mice, sewage, and termites and a mock community were analyzed by 454 and MiSeq sequencing of the V4-V5 region and MiSeq sequencing of the V4 region. Our analysis revealed that reference-based OTU clustering alone introduced biases compared to de novo clustering, preventing certain taxa from being observed in some samples. Based on this we devised and recommend an analysis pipeline that includes read merging, contaminant filtering, and reference-based clustering followed by de novo OTU clustering, which produces diversity measures consistent with de novo OTU clustering analysis. Low levels of dataset contamination with Illumina sequencing were discovered that could affect analyses that require highly sensitive approaches. While moving to Illumina-based sequencing platforms promises to provide deeper insights into the breadth and function of microbial diversity, our results show that care must be taken to ensure that sequencing and processing artifacts do not obscure true microbial diversity.     

Selection of primers for optimal taxonomic classification of environmental 16S rRNA gene sequences

Microbial community profiling using 16S rRNA gene sequences requires accurate taxonomy assignments. ‘Universal'' primers target conserved sequences and amplify sequences from many taxa, but they provide variable coverage of different environments, and regions of the rRNA gene differ in taxonomic informativeness—especially when high-throughput short-read sequencing technologies (for example, 454 and Illumina) are used. We introduce a new evaluation procedure that provides an improved measure of expected taxonomic precision when classifying environmental sequence reads from a given primer. Applying this measure to thousands of combinations of primers and read lengths, simulating single-ended and paired-end sequencing, reveals that these choices greatly affect taxonomic informativeness. The most informative sequence region may differ by environment, partly due to variable coverage of different environments in reference databases. Using our Rtax method of classifying paired-end reads, we found that paired-end sequencing provides substantial benefit in some environments including human gut, but not in others. Optimal primer choice for short reads totaling 96 nt provides 82–100% of the confident genus classifications available from longer reads.     

Critical evaluation of two primers commonly used for amplification of bacterial 16S rRNA genes

rRNA-based studies, which have become the most common method for assessing microbial communities, rely upon faithful amplification of the corresponding genes from the original DNA sample. We report here an analysis and reevaluation of commonly used primers for amplifying the DNA between positions 27 and 1492 of bacterial 16S rRNA genes (numbered according to the Escherichia coli rRNA). We propose a formulation for a forward primer (27f) that includes three sequences not usually present. We compare our proposed formulation to two common alternatives by using linear amplification-providing an assessment that is independent of a reverse primer-and in combination with the 1492 reverse primer (1492r) under the PCR conditions appropriate for making community rRNA gene clone libraries. For analyses of DNA from human vaginal samples, our formulation was better at maintaining the original rRNA gene ratio of Lactobacillus spp. to Gardnerella spp., particularly under stringent amplification conditions. Because our 27f formulation remains relatively simple, having seven distinct primer sequences, there is minimal loss of overall amplification efficiency and specificity.     

基于不同通用引物比较分析肠道微生物群落变化

肠道微生物对于人体健康的重要作用已经得到广泛证实,目前,对肠道微生物的研究大多采用基于扩增细菌16S rRNA基因V3-V4区的高通量测序分析,对古菌的关注较少。本研究选择了一对可以同时扩增细菌和古菌16S rRNA基因的引物,通过比较人为干扰肠道微生物前后的群落变化,说明这对引物适宜分析人类肠道细菌和古菌群落变化并具有一定优越性。采集志愿者粪便样品,同时用仅能扩增细菌引物 (B引物) 和细菌古菌通用引物 (AB引物) 进行扩增和高通量测序;使用几个常用的rRNA数据库判断引物对细菌的覆盖度和对古菌的扩增能力。结果表明,AB引物在可以展示B引物扩增出的细菌群落的基础上,可以得到肠道中常见的产甲烷古菌的序列,同时也展示出人为干扰肠道微生物前后的群落结构变化。AB引物可以仅通过一次扩增和测序同时分析肠道中的细菌和古菌群落,更加全面展示肠道微生物群落结构,适用于肠道微生物相关研究。     

Systematic Design of 18S rRNA Gene Primers for Determining Eukaryotic Diversity in Microbial Consortia

High-throughput sequencing of ribosomal RNA gene (rDNA) amplicons has opened up the door to large-scale comparative studies of microbial community structures. The short reads currently produced by massively parallel sequencing technologies make the choice of sequencing region crucial for accurate phylogenetic assignments. While for 16S rDNA, relevant regions have been well described, no truly systematic design of 18S rDNA primers aimed at resolving eukaryotic diversity has yet been reported. Here we used 31,862 18S rDNA sequences to design a set of broad-taxonomic range degenerate PCR primers. We simulated the phylogenetic information that each candidate primer pair would retrieve using paired- or single-end reads of various lengths, representing different sequencing technologies. Primer pairs targeting the V4 region performed best, allowing discrimination with paired-end reads as short as 150 bp (with 75% accuracy at genus level). The conditions for PCR amplification were optimised for one of these primer pairs and this was used to amplify 18S rDNA sequences from isolates as well as from a range of environmental samples which were then Illumina sequenced and analysed, revealing good concordance between expected and observed results. In summary, the reported primer sets will allow minimally biased assessment of eukaryotic diversity in different microbial ecosystems.     

Soil Bacterial Diversity Screening Using Single 16S rRNA Gene V Regions Coupled with Multi-Million Read Generating Sequencing Technologies

The novel multi-million read generating sequencing technologies are very promising for resolving the immense soil 16S rRNA gene bacterial diversity. Yet they have a limited maximum sequence length screening ability, restricting studies in screening DNA stretches of single 16S rRNA gene hypervariable (V) regions. The aim of the present study was to assess the effects of properties of four consecutive V regions (V3-6) on commonly applied analytical methodologies in bacterial ecology studies. Using an in silico approach, the performance of each V region was compared with the complete 16S rRNA gene stretch. We assessed related properties of the soil derived bacterial sequence collection of the Ribosomal Database Project (RDP) database and concomitantly performed simulations based on published datasets. Results indicate that overall the most prominent V region for soil bacterial diversity studies was V3, even though it was outperformed in some of the tests. Despite its high performance during most tests, V4 was less conserved along flanking sites, thus reducing its ability for bacterial diversity coverage. V5 performed well in the non-redundant RDP database based analysis. However V5 did not resemble the full-length 16S rRNA gene sequence results as well as V3 and V4 did when the natural sequence frequency and occurrence approximation was considered in the virtual experiment. Although, the highly conserved flanking sequence regions of V6 provide the ability to amplify partial 16S rRNA gene sequences from very diverse owners, it was demonstrated that V6 was the least informative compared to the rest examined V regions. Our results indicate that environment specific database exploration and theoretical assessment of the experimental approach are strongly suggested in 16S rRNA gene based bacterial diversity studies.     

Application of new primer-enzyme combinations to terminal restriction fragment length polymorphism profiling of bacterial populations in human feces

New primer-enzyme combinations for terminal restriction fragment length polymorphism (T-RFLP) targeting of the 16S rRNA gene were constructed by using the T-RFLP analysis program (designated TAP T-RFLP) located at the Ribosomal Database Project website, and their performance was examined empirically. By using the fluorescently labeled 516f primer (Escherichia coli positions 516 to 532) and 1510r primer (positions 1510 to 1492), the 16S rRNA gene was amplified from human fecal DNA. The resulting amplified product was digested with RsaI plus BfaI or with BslI. When the T-RFLP was carried out with fecal DNAs from eight individuals, eight predominant operational taxonomic units (OTUs) were detected with RsaI and BfaI digestion and 14 predominant OTUs were detected with BslI digestion. The distribution of the OTUs was consistent with the results of the computer simulations with TAP T-RFLP. The T-RFLP analyses of the fecal DNAs from individuals gave characteristic profiles, while the variability of the T-RFLP profiles between duplicate DNA preparations from the same samples were minimal. This new T-RFLP method made it easy to predict what kind of intestinal bacterial group corresponded to each OTU on the basis of the terminal restriction fragment length compared with the conventional T-RFLP and, moreover, made it possible to identify the bacterial species that an OTU represents by cloning and sequencing.     

Analysis of 16S rRNA amplicon sequencing options on the Roche/454 next-generation titanium sequencing platform

Background

16S rRNA gene pyrosequencing approach has revolutionized studies in microbial ecology. While primer selection and short read length can affect the resulting microbial community profile, little is known about the influence of pyrosequencing methods on the sequencing throughput and the outcome of microbial community analyses. The aim of this study is to compare differences in output, ease, and cost among three different amplicon pyrosequencing methods for the Roche/454 Titanium platform

Methodology/Principal Findings

The following three pyrosequencing methods for 16S rRNA genes were selected in this study: Method-1 (standard method) is the recommended method for bi-directional sequencing using the LIB-A kit; Method-2 is a new option designed in this study for unidirectional sequencing with the LIB-A kit; and Method-3 uses the LIB-L kit for unidirectional sequencing. In our comparison among these three methods using 10 different environmental samples, Method-2 and Method-3 produced 1.5–1.6 times more useable reads than the standard method (Method-1), after quality-based trimming, and did not compromise the outcome of microbial community analyses. Specifically, Method-3 is the most cost-effective unidirectional amplicon sequencing method as it provided the most reads and required the least effort in consumables management.

Conclusions

Our findings clearly demonstrated that alternative pyrosequencing methods for 16S rRNA genes could drastically affect sequencing output (e.g. number of reads before and after trimming) but have little effect on the outcomes of microbial community analysis. This finding is important for both researchers and sequencing facilities utilizing 16S rRNA gene pyrosequencing for microbial ecological studies.     

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.