Coverage evaluation of universal bacterial primers using the metagenomic datasets
© Mao et al.; licensee BioMed Central Ltd. 2012
Received: 28 October 2011
Accepted: 3 May 2012
Published: 3 May 2012
The coverage of universal primers for the bacterial 16S rRNA gene plays a crucial role in the correct understanding of microbial community structure. However, existing studies on primer coverage are limited by the lack of appropriate databases and are restricted to the domain level. Additionally, most studies do not account for the positional effect of single primer-template mismatches. In this study, we used 7 metagenomic datasets as well as the Ribosomal Database Project (RDP) to assess the coverage of 8 widely used bacterial primers.
The coverage rates for bacterial primers were found to be overestimated by previous studies that only investigated the RDP because of PCR amplification bias in the sequence composition of the dataset. In the RDP, the non-coverage rates for all primers except 27F were ≪6%, while in the metagenomic datasets, most were ≫10%. If one considers that a single mismatch near the 3′ end of the primer might greatly reduce PCR efficiency, then some phylum non-coverage rates would change by more than 20%. Primer binding-site sequence variants that could not pair with their corresponding primers are discussed.
Our study revealed the potential bias introduced by the use of universal bacterial primers in the assessment of microbial communities. With the development of high-throughput, next-generation sequencing techniques, it will become feasible to sequence more of the hypervariable regions of the bacterial 16S rRNA gene. This, in turn, will lead to the more frequent use of the primers discussed here.
In the field of microbial ecology, the polymerase chain reaction (PCR) has been widely used for the amplification, detection and quantification of DNA targets since its introduction [1, 2], resulting in increased knowledge of the microbial world [3, 4]. However, the efficiency and accuracy of PCR can be diminished by many factors including primer-template mismatches, reactant concentrations, the number of PCR cycles, annealing temperature, the complexity of the DNA template, and others. [5–7]. Primer-template mismatches are the most important because they can lead to selective amplification which prevents the correct assessment of microbial diversity [8, 9]. Target sequences that cannot match the primers precisely will be amplified to a lesser extent, possibly even below the detection limit. The relative content of the sequences achieved is therefore changed, resulting in a deviation from the true community composition. Hence a comprehensive evaluation of bacterial primer coverage is critical to the interpretation of PCR results in microbial ecology research.
Many related studies on primer coverage have been performed previously, but most are qualitative or semi-quantitative studies restricted to the domain level [10, 11]. Low coverage rates in some rare phyla might have been overlooked.
Although Wang et al.  investigated primer coverage rates at the phylum level, only sequences from the Ribosomal Database Project (RDP) were used. This sole reliance on the RDP is another common limitation of previous studies. The RDP is a professional database containing more than one million 16S rRNA gene sequences. It also provides a series of data analysis services [13, 14], including Probe Match, which is often used in primer studies. However, despite the RDP’s large collection of sequences and extensive application, most of its sequences were generated through PCR amplification. Sequences that fail to match the universal primers may become lost in the PCR results, and so are not included in the RDP. Consequently, primer coverage rates in the RDP appear to be higher than they actually are.
Fortunately, with the rapid development of sequencing techniques, many large-scale metagenomic datasets have become available. Metagenomic sequences are generated directly from sequencing environmental samples and are free of PCR bias; thus, the resulting datasets faithfully reflect microbial composition, especially in the case of rare biospheres. The Community Cyberinfrastructure for Advanced Microbial Ecology Research and Analysis (CAMERA) is not only a repository for rich and distinctive metagenomic data, but it also provides a set of bioinformatic tools for research.
Another shortcoming of previous primer-coverage studies has recently been illuminated through studies on the PCR mechanism. In the past, it was assumed that a single primer-template mismatch would not obstruct amplification under proper annealing temperature so long as the mismatch did not occur at the 3′ end of the primer. However, recent studies have shown that a single mismatch within the last 3–4 nucleotides of the 3′ end could also significantly reduce PCR amplification efficiency, even under optimal annealing temperature [16, 17]. This changed the criteria for judging whether a primer binding-site sequence could be amplified faithfully by PCR. In this study, we define sequences that “match with” the primers as having either no mismatch with the primer, or as having only one mismatch that is not located within the last 4 nucleotides of the 3′ end.
All of the primers in this study are frequently used in molecular microbial ecology research. The most common primer pairs are 27F and 1390R/1492R, which are mainly used for constructing clone libraries of the full-length 16S rDNA sequence . The primers such as 338F and 338R are frequently used in pyrosequencing [19–21]. The remaining primers are most commonly used for fingerprint analyses, but the development of next-generation sequencing techniques will likely broaden their roles in future studies [22, 23]. Pyrosequencing has extended the read length from 100bp to 800bp , and as a result, hypervariable regions in 16S rDNA other than V6 and V3 will be able to be sequenced. Those primers that can cover these hypervariable regions will become more frequently used.
The aim of this study was to assess the coverage rates of 8 common primers (27F, 338F, 338R, 519F, 519R, 907R, 1390R and 1492R), which target different regions of the bacterial 16S rRNA gene, using sequences from the RDP and 7 metagenomic datasets. We used the non-coverage rate, the percentage of sequences that could not match with the primer, as the major indicator in this study. Non-coverage rates were calculated at both the domain and phylum levels, and the influence of a single mismatched position on the non-coverage rate was analyzed. By comparing the RDP and the metagenomic datasets, we found that the non-coverage rates were seriously underestimated when only the RDP dataset was used.
Results and discussion
Influence of a single mismatch in the last 4 nucleotides
Since the beginning of the 1990s, it has been widely acknowledged that PCR amplification is significantly inhibited by a single mismatch occurring at the 3′ end of the primer [25–27]. Even when the last nucleotide was substituted with inosine, which is capable of binding to all four nucleotides, primers still failed to amplify all of the expected sequences in the microbial community . Recently, Bru et al.  and Wu et al.  demonstrated that the efficiency of PCR amplification was also inhibited if a single mismatch occurred within the last 3–4 nucleotides of the 3′ end of primer, even when the annealing temperature was decreased for optimal efficiency. These single mismatches have not been considered in previous primer coverage studies [12, 18, 29].
Influence of a single mismatch near the 3′ end in the RDP dataset
Non-coverage rate 4+ (%)
Non-coverage rate 4- (%)
Non-coverage rates of 8 primers at the domain level
Non-coverage rates for 8 primers at the phylum level
Because each dataset is a mixture of sequences from various microbes occurring in various proportions according to different phyla, low coverage of minor phyla could be easily masked by the higher coverage of the dominant phyla. Moreover, the compositions of microbial communities differ greatly with environments; Minor microbes found in common environments may in fact be major components in other ecological niches. It is therefore necessary to assess the non-coverage rates at the phylum level in the different metagenomic datasets.
338F and 338R
Non-coverage rates for the primers 338F and 338R varied among different phyla (Additional file 2: Figure S2.). In the RDP dataset, the non-coverage rates for 338F in 4 phyla (Aquificae, Planctomycetes, Verrucomicrobia and OD1) were ≫95%. Primer binding-site sequences that could not match with primer 338F are listed in Additional file 3: Table S2.
In the RDP dataset, the most frequent sequence variant retrieved (3,587 sequences) was 338F-3A12T (3A indicates that the 3rd base is the nucleotide A, and 12T that the 12th base is the nucleotide T). This sequence was the major variant in the Verrucomicrobia, accounting for 97.8% of the sequences in the RDP dataset and 85.7% in the GOS (Global Ocean Sampling Expedition) dataset; it also predominated in the phyla Chloroflexi, BRC1, OP10 and OP11. The second variant, 338F-16T, was the major variant in the Lentisphaerae but also appeared in many other phyla. The third variant, 338F-3A12T16T, was specific for Planctomycetes and OD1, and accounted for approximately 50% of Planctomycetes in both the RDP and GOS datasets. The variants 338F-4T11A and 338F-12G were distributed in various phyla, while 338F-3C12G was specific for Aquificae and 338F-3C4T11A12G for Cyanobacteria.
Also significant was the non-coverage rate for 338F in the Actinobacteria. In the RDP dataset, this rate was only 1.3%, but in the metagenomic datasets, the results were substantially different. The non-coverage rates in the GOS and HOT datasets, for example, were 60.4% and 66.7%, respectively. We observed that the absolute number of 338F-16T sequences from Actinobacteria in the RDP dataset was 631, which was much larger than the numbers in the GOS and HOT datasets. The implication is that the 338F-16T Actinobacteria sequences in the RDP most likely came from environments similar to those from which the GOS and HOT sequences were sampled.
For the primer 338R, the reverse complement of 338F, the homologous variants 338F-16T and 338F-16C had no effect on the non-coverage rate, while three other variants (338R-16G, 338R-18C and 338R-15A) warranted further attention (Additional file 3: Table S3). Although hundreds of sequences for each variant were found, they accounted for low percentages of the major phyla (Actinobacteria, Bacteroidetes, Firmicutes and Proteobacteria). Variants with more than one mismatch were similar to those of 338F.
The BisonMetagenome dataset was dominated by Aquificae and the non-coverage rates for both 338F and 338R in Aquificae were 100%. The sequence variant 338F-3C12G (338R-7C16G) was the major type. Thus, the primers 338F/338R might not be appropriate for the analysis of hot spring samples or the detection of Aquificae.
519F and 519R
The coverage of primer 519R was quite “universal” except for its high non-coverage rate in the phylum OD1 in the AntarcticaAquatic dataset, where the primer binding-site sequence variant 519R-14T-11T12C had a rate of 84.6%. Although non-coverage rates of approximately 20% were found scattered across other phyla, these rates resulted from variants with only one or two sequences, and no dominating variant was found. Overall, primer 519R could authentically amplify sequences from most phyla.
A substantial difference was found between the non-coverage rates of 519F and 519R. Five sequence variants were mainly responsible for the high non-coverage rate for 519F (Additional file 3: Table S4). Notably, the 3 most dominant variants had one trait in common – a single mismatch at the 16th nucleotide (the 3rd nucleotide from the 3′ end of 519F). This mismatch did not influence the non-coverage rate of 519R.
Further analysis showed that the high non-coverage rate of 519F was caused primarily by sequences from the phylum Nitrospirae. The AcidMine metagenome is dominated by Leptospirillum species of the Nitrospirae, and therefore forms an ideal dataset for Nitrospirae studies . Of the 519F-binding sequences in the dataset, 89% were from Nitrospirae, and none could match with 519F. The non-coverage rate in the RDP dataset was also high (68%) in Nitrospirae, whereas the total non-coverage rate for 519F in the RDP dataset was only 6%. Similar sample analyses should therefore be focused on the use of primer 519F.
Frank et al.  have studied the 27F and 1492R primer pair and have proposed 27F-YM + 3 as a modification of the common 27F primer. Our results support this modification as being necessary (Additional file 3: Table S1). The non-coverage rates for 1390R and 1492R were quite low, even at the phylum level. For primer 907R, only one sequence variant that could not match with the primer (907R-11C-15A16T) was observed. It resulted in the high non-coverage rate observed in phylum TM7 (Additional file 3: Table S5).
The 16S rRNA gene is an important genetic marker for the characterization of microbial community structure by 16S rRNA gene amplicon sequencing with conserved primers . Because of the increase in read length with the development of pyrosequencing (454 sequencing) technology, different multi-hypervariable regions can be selected for amplification. In this strategy, different pairs of “universal” primers are used for barcoded pyrosequencing . However, even with pyrosequencing, the bias caused by primer-template mismatch may misrepresent the real community composition of environmental samples. Therefore, the assessment of primer coverage to perfect the use of universal primers is urgently required.
In this study, we assessed the non-coverage rates for 8 common universal bacterial primers in the RDP dataset and 7 metagenomic datasets. Comparisons of non-coverage rates, with or without constraining the position of a single mismatch, emphasized the importance of further study of the mechanism of PCR. Metagenomic dataset analysis revealed that some sequence variants, which appeared to be minor in the public databases, were actually dominant in some ecological niches. These results are of great practical significance for studies on similar environmental samples, and new primer formulations could be designed using our results. One strategy is to increase coverage through the introduction of proper degenerate nucleotides.
Although the total number of sequences in a metagenomic dataset may be very large, the number of 16S rRNA gene sequences is limited, and may account for only approximately 0.2% of all sequence reads [33, 34]. In contrast, the metatranscriptomic analysis of environmental samples generates a large number of small subunit sequences . Although the short length (approximately 200bp) of the sequences currently deposited in metatranscriptomic datasets are not appropriate for assessing primer coverage, the further development of pyrosequencing will make such assessments possible in the near future.
Retrieval of 16S rRNA gene sequences from the RDP
A FASTA file for all bacterial 16S rRNA gene sequences was downloaded from the “RESOURCES” section of the RDP website (release 10.18; http://rdp.cme.msu.edu/) . With the help of the service “BROWSERS”, good quality, almost full-length (size ≥ 1200bp) sequences were obtained. These sequences were extracted from the FASTA file by Perl scripts. A final dataset with 462,719 bacterial 16S rRNA gene sequences was constructed (referred to as the “RDP dataset”).
Elimination of primer contamination in the RDP dataset
Most sequences deposited in the RDP dataset were generated by PCR. However, as described by Frank et al. , many of these sequences lack correct primer trimming. Only sequence fragments extending at least 3 nucleotides past the start (the 5′ end) of the longest version of each primer were considered uncontaminated by the PCR primers. Because the sequences selected from the RDP were all longer than 1200bp, only the primer-binding sites for 27F, 1390R and 1492R could be contaminated (Additional file 4: Figure S3). Thus, 15,045, 188,792 and 35,462 sequences were selected for the primers 27F, 1390R and 1492R, respectively, as containing authentic primer-binding sites.
Retrieval of 16S rDNA sequences from the metagenomic datasets
Selection of metagenomic datasets
Descriptions of the metagenomic datasets
Antarctica Aquatic Microbial Metagenome
(All Metagenomic Shotgun Reads)
Acid Mine Drainage Metagenome
(All Metagenomic Shotgun Reads)
Metagenome from Yellowstone Bison Hot Spring
(All Metagenomic Shotgun Reads)
Global Ocean Sampling Expedition
(All Metagenomic Sequence Reads)
Mediterranean Gutless Worm Metagenome
(All Metagenomic Sequence Reads)
Human Distal Gut Biome project
Microbial Community Genomics at the Hawaii Ocean Time-series (HOT) station ALOHA
(All Metagenomic Sequence Reads)
Retrieval of 16S rDNA homologs
The Basic Local Alignment Search Tool (BLAST) was used to acquire as many 16S rRNA gene homologs as possible for the low content of such sequences in the metagenomic datasets. A query set of 34 representative and almost full-length 16S rRNA gene sequences from 34 bacterial phyla was constructed. BLAST searches using the query set and each selected dataset were performed using the CAMERA interface (db alignments per query, 50000; e-value exponent (1Ex), -5; filter low-complexity seq, T; lower case filtering, False). For the GOS dataset, BLAST was performed using each query sequence separately because the subjects exceeded the threshold of “db alignments per query” when BLAST was performed using the complete query set. After removing reads containing the nucleotide “N”, sequence reads were merged into one file without duplication. Seven files were obtained, one from each of the 7 datasets.
Further filtration of 16S rDNA homologs
The software program Mothur (http://www.mothur.org) was used for further filtration . Sequences and their reverse complements were aligned separately via the command “align.seqs”. One reference file containing large subunit rRNA gene sequences was downloaded from Silva (http://www.arb-silva.de/) . The second reference file was a combination of Silva reference files of small subunit rRNA gene sequences downloaded from Mothur. According to the alignment scores, the origin and direction of the sequences were ascertained. Sequences whose scores were always ≪30 might represent non-rRNA genes and were therefore removed.
For the RDP dataset, the alignment with the reference file of small subunit rDNA sequences was run first, and sequences with alignment scores ≪30 were removed.
The 16S rRNA gene sequences from both the RDP dataset and the metagenomic datasets were assigned to different taxonomic groups by Mothur, with the confidence threshold set at 80%. Sequences classified as belonging to the domain Bacteria were listed and extracted.
Identification of primer-binding sites in 16S rDNA sequences
Because the alignment using the Silva template sequences did not include the entire length of the gene, thus missing the primer-binding sites for 27F and 1492R, alignment with another reference file (the “Core Set” of the Greengenes database) was used to identify the primer-binding sites . A full-length 16S rRNA gene sequence from Escherichia coli (GenBank ID: J01695) was added for base positioning.
Detailed information for the 8 primers evaluated
Sequence of primer
Position in Escherichia coli
27 F (8 F)
5′- AGA GTT TGA TYM TGG CTC AG-3′
5′-ACT CCT ACG GGA GGC AGC-3′
5′-GCT GCC TCC CGT AGG AGT-3′
5′-CAG CM G CCG CGG TAA TAC-3′
5′-GTA TTA CCG CGG CK G CTG-3′
5′-CCG TCA ATT CM T TTG AGT TT-3′
5′-ACG GGC GGT GTG TR C AA-3′
5′-TAC CTT GTT AY G ACT T-3′
Primer binding-site sequences with more than one mismatch, or with a single mismatch within the last 4 nucleotides of the 3′ end, were considered unmatched with the primer. Non-coverage rates were calculated as the percentage of such sequences. The non-coverage rates of phyla with sequence numbers of less than 50 in the RDP dataset or less than 10 in the metagenomic datasets were not shown in Figure 1 and Additional file 2: Figure S2.
Because different phyla vary considerably in the numbers of sequences reported, we attempted a normalization approach to calculate the non-coverage rates for each dataset. Phyla with less than 10 sequences or 1% of the total of each dataset were merged into a new “phylum”. The domain non-coverage rate was computed as the arithmetical average of the phylum non-coverage rates.
This work was supported by the National Key Technology R&D Program of China (2006BAI19B02) and the National High Technology Research and Development Program of China (2008AA062501-2).
- Olsen GJ, Lane DJ, Giovannoni SJ, Pace NR, Stahl DA: Microbial ecology and evolution: a ribosomal RNA approach. Annu Rev Microbiol. 1986, 40: 337-365. 10.1146/annurev.mi.40.100186.002005.PubMedView ArticleGoogle Scholar
- Schmidt TM, Delong EF, Pace NR: Analysis of a marine picoplankton community by 16S rRNA gene cloning and sequencing. J Bacteriol. 1991, 173: 4371-4378.PubMedPubMed CentralGoogle Scholar
- Sharkey FH, Banat IM, Marchant R: Detection and quantification of gene expression in environmental bacteriology. Appl Environ Microb. 2004, 70: 3795-3806. 10.1128/AEM.70.7.3795-3806.2004.View ArticleGoogle Scholar
- Steffan RJ, Atlas RM: Polymerase chain reaction: applications in environmental microbiology. Annu Rev Microbiol. 1991, 45: 137-161. 10.1146/annurev.mi.45.100191.001033.PubMedView ArticleGoogle Scholar
- Forney LJ, Zhou X, Brown CJ: Molecular microbial ecology: land of the one-eyed king. Curr Opin Microbiol. 2004, 7: 210-220. 10.1016/j.mib.2004.04.015.PubMedView ArticleGoogle Scholar
- Smith S, Vigilant L, Morin PA: The effects of sequence length and oligonucleotide mismatches on 5' exonuclease assay efficiency. Nucleic Acids Res. 2002, 30: e111-10.1093/nar/gnf110.PubMedPubMed CentralView ArticleGoogle Scholar
- von Wintzingerode F, Gobel UB, Stackebrandt E: Determination of microbial diversity in environmental samples: pitfalls of PCR-based rRNA analysis. FEMS Microbiol Rev. 1997, 21: 213-229. 10.1111/j.1574-6976.1997.tb00351.x.PubMedView ArticleGoogle Scholar
- Polz MF, Cavanaugh CM: Bias in template-to-product ratios in multitemplate PCR. Appl Environ Microb. 1998, 64: 3724-3730.Google Scholar
- Reysenbach AL, Giver LJ, Wickham GS, Pace NR: Differential amplification of rRNA genes by polymerase chain reaction. Appl Environ Microb. 1992, 58: 3417-3418.Google Scholar
- Baker GC, Smith JJ, Cowan DA: Review and re-analysis of domain-specific 16S primers. J Microbiol Meth. 2003, 55: 541-555. 10.1016/j.mimet.2003.08.009.View ArticleGoogle Scholar
- Huws SA, Edwards JE, Kim EJ, Scollan ND: Specificity and sensitivity of eubacterial primers utilized for molecular profiling of bacteria within complex microbial ecosystems. J Microbiol Meth. 2007, 70: 565-569. 10.1016/j.mimet.2007.06.013.View ArticleGoogle Scholar
- Wang Y, Qian PY: Conservative fragments in bacterial 16S rRNA genes and primer design for 16S ribosomal DNA amplicons in metagenomic studies. PLoS One. 2009, 4: e7401-10.1371/journal.pone.0007401.PubMedPubMed CentralView ArticleGoogle Scholar
- Cole JR, Chai B, Farris RJ, Wang Q, Kulam SA, McGarrell DM, Garrity GM, Tiedje JM: The Ribosomal Database Project (RDP-II): sequences and tools for high-throughput rRNA analysis. Nucleic Acids Res. 2005, 33: D294-D296.PubMedPubMed CentralView ArticleGoogle Scholar
- Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-Syed-Mohideen AS, McGarrell DM, Marsh T, Garrity GM, Tiedje JM: The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009, 37: D141-D145. 10.1093/nar/gkn879.PubMedPubMed CentralView ArticleGoogle Scholar
- Seshadri R, Kravitz SA, Smarr L, Gilna P, Frazier M: CAMERA: a community resource for metagenomics. PLoS Biol. 2007, 5: 394-397.View ArticleGoogle Scholar
- Bru D, Martin-Laurent F, Philippot L: Quantification of the detrimental effect of a single primer-template mismatch by real-time PCR using the 16S rRNA gene as an example. Appl Environ Microb. 2008, 74: 1660-1663. 10.1128/AEM.02403-07.View ArticleGoogle Scholar
- Wu JH, Hong PY, Liu WT: Quantitative effects of position and type of single mismatch on single base primer extension. J Microbiol Meth. 2009, 77: 267-275. 10.1016/j.mimet.2009.03.001.View ArticleGoogle Scholar
- Frank JA, Reich CI, Sharma S, Weisbaum JS, Wilson BA, Olsen GJ: Critical evaluation of two primers commonly used for amplification of bacterial 16S rRNA genes. Appl Environ Microb. 2008, 74: 2461-2470. 10.1128/AEM.02272-07.View ArticleGoogle Scholar
- Humblot C, Guyot J-P: Pyrosequencing of tagged 16S rRNA gene amplicons for rapid deciphering of the microbiomes of fermented foods such as pearl millet slurries. Appl Environ Microb. 2009, 75: 4354-4361. 10.1128/AEM.00451-09.View ArticleGoogle Scholar
- Forney LJ, Gajer P, Williams CJ, Schneider GM, Koenig SSK, McCulle SL, Karlebach S, Brotman RM, Davis CC, Ault K, Ravel J: Comparison of self-collected and physician-collected vaginal swabs for microbiome analysis. J Clin Microbiol. 2010, 48: 1741-1748. 10.1128/JCM.01710-09.PubMedPubMed CentralView ArticleGoogle Scholar
- Lauber CL, Hamady M, Knight R, Fierer N: Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl Environ Microb. 2009, 75: 5111-5120. 10.1128/AEM.00335-09.View ArticleGoogle Scholar
- Bai YH, Sun QH, Zhao C, Wen DH, Tang XY: Bioaugmentation treatment for coking wastewater containing pyridine and quinoline in a sequencing batch reactor. Appl Microbiol Biot. 2010, 87: 1943-1951. 10.1007/s00253-010-2670-8.View ArticleGoogle Scholar
- Tan YF, Ji GD: Bacterial community structure and dominant bacteria in activated sludge from a 70 degrees C ultrasound-enhanced anaerobic reactor for treating carbazole-containing wastewater. Bioresource Technol. 2010, 101: 174-180. 10.1016/j.biortech.2009.08.044.View ArticleGoogle Scholar
- Miller W, Hayes VM, Ratan A, Petersen DC, Wittekindt NE, Miller J, Walenz B, Knight J, Qi J, Zhao F, et al.: Genetic diversity and population structure of the endangered marsupialSarcophilus harrisii(Tasmanian devil). P Natl Acad Sci USA. 2011, 108: 12348-12353. 10.1073/pnas.1102838108.View ArticleGoogle Scholar
- Ayyadevara S, Thaden JJ, Reis RJS: Discrimination of primer 3'-nucleotide mismatch by Taq DNA polymerase during polymerase chain reaction. Anal Biochem. 2000, 284: 11-18. 10.1006/abio.2000.4635.PubMedView ArticleGoogle Scholar
- Huang MM, Arnheim N, Goodman MF: Extension of base mispairs by Taq DNA polymerase: implications for single nucleotide discrimination in PCR. Nucleic Acids Res. 1992, 20: 4567-4573. 10.1093/nar/20.17.4567.PubMedPubMed CentralView ArticleGoogle Scholar
- Kwok S, Kellogg DE, McKinney N, Spasic D, Goda L, Levenson C, Sninsky JJ: Effects of Primer Template Mismatches on the Polymerase Chain-Reaction - Human-Immunodeficiency-Virus Type-1 Model Studies. Nucleic Acids Res. 1990, 18: 999-1005. 10.1093/nar/18.4.999.PubMedPubMed CentralView ArticleGoogle Scholar
- Brands B, Vianna ME, Seyfarth I, Conrads G, Horz HP: Complementary retrieval of 16S rRNA gene sequences using broad-range primers with inosine at the 3'-terminus: implications for the study of microbial diversity. FEMS Microbiol Ecol. 2009, 71: 157-167.View ArticleGoogle Scholar
- Daims H, Bruhl A, Amann R, Schleifer KH, Wagner M: The domain-specific probe EUB338 is insufficient for the detection of allBacteria: development and evaluation of a more comprehensive probe set. Syst Appl Microbiol. 1999, 22: 434-444. 10.1016/S0723-2020(99)80053-8.PubMedView ArticleGoogle Scholar
- Tyson GW, Chapman J, Hugenholtz P, Allen EE, Ram RJ, Richardson PM, Solovyev VV, Rubin EM, Rokhsar DS, Banfield JF: Community structure and metabolism through reconstruction of microbial genomes from the environment. Nature. 2004, 428: 37-43. 10.1038/nature02340.PubMedView ArticleGoogle Scholar
- Schmalenberger A, Schwieger F, Tebbe CC: Effect of Primers Hybridizing to Different Evolutionarily Conserved Regions of the Small-Subunit rRNA Gene in PCR-Based Microbial Community Analyses and Genetic Profiling. Appl Environ Microb. 2001, 67: 3557-3563. 10.1128/AEM.67.8.3557-3563.2001.View ArticleGoogle Scholar
- Petrosino JF, Highlander S, Luna RA, Gibbs RA, Versalovic J: Metagenomic Pyrosequencing and Microbial Identification. Clin Chem. 2009, 55: 856-866. 10.1373/clinchem.2008.107565.PubMedPubMed CentralView ArticleGoogle Scholar
- Biers EJ, Sun SL, Howard EC: Prokaryotic genomes and diversity in surface ocean waters: interrogating the global ocean sampling metagenome. Appl Environ Microb. 2009, 75: 2221-2229. 10.1128/AEM.02118-08.View ArticleGoogle Scholar
- Mou XZ, Sun SL, Edwards RA, Hodson RE, Moran MA: Bacterial carbon processing by generalist species in the coastal ocean. Nature. 2008, 451: 708-711. 10.1038/nature06513.PubMedView ArticleGoogle Scholar
- Urich T, Lanzen A, Qi J, Huson DH, Schleper C, Schuster SC: Simultaneous assessment of soil microbial community structure and function through analysis of the meta-transcriptome. PLoS One. 2008, 3: e2527-10.1371/journal.pone.0002527.PubMedPubMed CentralView ArticleGoogle Scholar
- Lauro FM, DeMaere MZ, Yau S, Brown MV, Ng C, Wilkins D, Raftery MJ, Gibson JAE, Andrews-Pfannkoch C, Lewis M, et al.: An integrative study of a meromictic lake ecosystem in Antarctica. ISME J. 2011, 5: 879-895. 10.1038/ismej.2010.185.PubMedPubMed CentralView ArticleGoogle Scholar
- Swingley WD, Alsop EB, Falenski HD, Raymond J: The 470 megabase metagenome of the Bison Pool (Yellowstone National Park) Alkaline Hot Spring Outflow Channel. Ab Sci Con. 2010, 2010: 5525-Google Scholar
- Yutin N, Suzuki MT, Teeling H, Weber M, Venter JC, Rusch DB, Béjà O: Assessing diversity and biogeography of aerobic anoxygenic phototrophic bacteria in surface waters of the Atlantic and Pacific Oceans using the Global Ocean Sampling expedition metagenomes. Environ Microbiol. 2007, 9: 1464-1475. 10.1111/j.1462-2920.2007.01265.x.PubMedView ArticleGoogle Scholar
- Woyke T, Teeling H, Ivanova NN, Huntemann M, Richter M, Gloeckner FO, Boffelli D, Anderson IJ, Barry KW, Shapiro HJ, et al.: Symbiosis insights through metagenomic analysis of a microbial consortium. Nature. 2006, 443: 950-955. 10.1038/nature05192.PubMedView ArticleGoogle Scholar
- Gloux K, Berteau O, El oumami H, Béguet F, Leclerc M, Doré J: A metagenomic β-glucuronidase uncovers a core adaptive function of the human intestinal microbiome. P Natl Acad Sci USA. 2011, 108: 4539-4546. 10.1073/pnas.1000066107.View ArticleGoogle Scholar
- DeLong EF, Preston CM, Mincer T, Rich V, Hallam SJ, Frigaard N-U, Martinez A, Sullivan MB, Edwards R, Brito BR, et al.: Community genomics among stratified microbial assemblages in the ocean's interior. Science. 2006, 311: 496-503. 10.1126/science.1120250.PubMedView ArticleGoogle Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al.: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microb. 2009, 75: 7537-7541. 10.1128/AEM.01541-09.View ArticleGoogle Scholar
- Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glöckner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007, 35: 7188-7196. 10.1093/nar/gkm864.PubMedPubMed CentralView ArticleGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microb. 2006, 72: 5069-5072. 10.1128/AEM.03006-05.View ArticleGoogle Scholar
- Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar , Buchner A, Lai T, Steppi S, Jobb G, et al.: ARB: a software environment for sequence data. Nucleic Acids Res. 2004, 32: 1363-1371. 10.1093/nar/gkh293.PubMedPubMed CentralView ArticleGoogle Scholar
- Ludwig W, Mittenhuber G, Friedrich CG: Transfer ofThiosphaera pantotrophatoParacoccus denitrificans. Int J Syst Bacteriol. 1993, 43: 363-367. 10.1099/00207713-43-2-363.PubMedView ArticleGoogle Scholar
- Whiteley AS, Bailey MJ: Bacterial community structure and physiological state within an industrial phenol bioremediation system. Appl Environ Microb. 2000, 66: 2400-2407. 10.1128/AEM.66.6.2400-2407.2000.View ArticleGoogle Scholar
- Suzuki MT, Giovannoni SJ: Bias caused by template annealing in the amplification of mixtures of 16S rRNA genes by PCR. Appl Environ Microb. 1996, 62: 625-630.Google Scholar
- Burggraf S, Huber H, Stetter KO: Reclassification of the crenarchaeal orders and families in accordance with 16S rRNA sequence data. Int J Syst Bacteriol. 1997, 47: 657-660. 10.1099/00207713-47-3-657.PubMedView ArticleGoogle Scholar
- Ruffroberts AL, Kuenen JG, Ward DM: Distribution of cultivated and uncultivatedcyanobacteriaandChloroflexus-like bacteria in hot spring microbial mats. Appl Environ Microb. 1994, 60: 697-704.Google Scholar
- Muyzer G, Teske A, Wirsen CO, Jannasch HW: Phylogenetic relationships ofThiomicrospiraspecies and their identification in deep-sea hydrothermal vent samples by denaturing gradient gel electrophoresis of 16S rDNA fragments. Arch Microbiol. 1995, 164: 165-172. 10.1007/BF02529967.PubMedView ArticleGoogle Scholar
- Brunk CF, Eis N: Quantitative measure of small-subunit rRNA gene sequences of the kingdomKorarchaeota. Appl Environ Microb. 1998, 64: 5064-5066.Google Scholar
- Wilson KH, Blitchington RB, Greene RC: Amplification of bacterial 16S ribosomal DNA with polymerase chain reaction. J Clin Microbiol. 1990, 28: 1942-1946.PubMedPubMed CentralGoogle Scholar
- Oguntoyinbo FA: Monitoring of marineBacillusdiversity among the bacteria community of sea water. Afr J Biotechnol. 2007, 6: 163-166.Google Scholar