Target Region Selection Is a Critical Determinant of Community Fingerprints Generated by 16S Pyrosequencing

Abstract

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, Campylobacer 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 providesd 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.

About the Authors:

Purnima S. Kumar

* E-mail: [email protected]

Affiliation: Division of Periodontology, College of Dentistry, The Ohio State University, Columbus, Ohio, United States of America

Michael R. Brooker

Affiliation: Division of Periodontology, College of Dentistry, The Ohio State University, Columbus, Ohio, United States of America

Scot E. Dowd

Affiliation: Research Testing Laboratory, Lubbock, Texas, United States of America

Terry Camerlengo

Affiliation: Comprehensive Cancer Center, The Ohio State University, Columbus, Ohio, United States of America

Introduction

Molecular approaches have revealed the presence of large numbers of as-yet-uncultivated organisms in the subgingival microbiome; creating a paradigm shift in our understanding of periodontal health and disease [1], [2], [3], [4]. In recent years, sequencing of 16S rRNA genes by the Sanger method (16S cloning and sequencing) has been widely used to examine subgingival microbial profiles in periodontal health and disease, as well as to characterize compositional shifts in these communities [5], [6], [7], [8], [9]. However, recent studies suggest that next-generation sequencing methodologies provide an economical and significantly higher-throughput alternative to Sanger sequencing for comparative genomics [10], [11].

Pyrosequencing of PCR-amplified 16S rDNA (‘16S pyrotags’) is a next-generation sequencing methodology that is capable of generating thousands of sequences from several samples simultaneously. The unprecedented sampling depth provided by this deep-sequencing approach allows the identification of several numerically minor or rare species within a community and has revealed a significantly greater level of microbial diversity than was previously apparent with Sanger sequencing [12], [13].

Unlike Sanger sequencing, which is capable of sequencing the entire gene, pyrosequencing is currently limited to generating sequences that are usually 350–500 bp in length. In order to improve community coverage, various investigations have employed primers that target different regions of the gene [12], [13], [14]. It has previously been shown, using Sanger sequencing, that the region of the 16S gene that is targeted for sequencing as well as the degeneracy of the sequencing primers introduce a level of bias into the community profile [2], [15]. Since pyrosequencing provides an enormously increased depth-of-coverage, it is important to understand the extent and severity of bias introduced by primer selection on the profile of any given community.

Previous studies have examined this bias using simulated datasets obtained by truncating full-length sequences, in silico testing of primer sequences for community coverage rates or by analyzing artificial bacterial communities created by mixing bacterial isolates [15], [16], [17], [18], [19]. However, it is logical to expect that fragment length (∼1.5 kb with Sanger sequencing and 150–500 bp with pyrosequencing) and as well as sequencing chemistry will affect amplification efficiency; therefore, profiles derived from artificially generated sequences may not accurately represent the coverage obtained from naturally occurring microbial communities. In fact, a recent investigation comparing 454 and Illumina sequencing has found significant divergence between in silico predictions and experimental results, emphasizing the need for experimental validation of primer pairs [20]. Hence, it is important to investigate the extent of this bias using sequences derived from clinical samples.

The purpose of this investigation, therefore, was to examine the bias introduced by target region selection and as well as by primer degeneracy on coverage of subgingival microbial communities using pyrosequencing.

Methods

Subject selection

Approval for this study was obtained from the Office of Responsible Research Practices at The Ohio State University. 10 current smokers with generalized moderate to severe chronic periodontitis were identified following clinical and radiographic examination and written informed consent was obtained. Exclusion criteria included diabetes, HIV infection, use of immunosuppressant medications, bisphosphonates or steroids, antibiotic therapy or oral prophylactic procedures within the last three months and less than 20 teeth in the dentition.

Sample collection and DNA isolation

Subgingival plaque samples were collected and pooled from four non-adjacent proximal sites demonstrating at least 6 mm of attachment loss and 5 mm of probe depths. Samples were collected by inserting 4 sterile endodontic paper points (Caulk-Dentsply) into each of the 4 sites for 10 seconds, following isolation and supragingival plaque removal. Samples were placed in 1.5 ml microcentrifuge tubes and frozen until further analysis. Bacteria were separated from the paper points by adding 200 µl of phosphate buffered saline (PBS) to the tubes and vortexing. The points were then removed, and DNA was isolated with a Qiagen DNA MiniAmp kit (Qiagen, Valencia, CA) using the tissue protocol according to the manufacturer's instructions.

Selection and optimization of primers

Four sets of primers were used to amplify each sample (A17 and 519R, 27F and 515R, 519F and 1114R, 1114F and 317). The primer sequences are listed in Table 1. Primer pairs were selected to generate 400–500 bp products from contiguous regions of the 16S rRNA gene. Previous sequencing-based investigations were examined and the primers most commonly used in these studies were selected [2], [6], [7], [8], [9], [13], [21], [22]. The universality of the primer pairs was assessed by comparing them to our locally hosted, curated database of 1800 nearly full-length 16S sequences derived from GenBank. MacVector was used for alignment and determining melting temperatures and GC ratios of the resulting amplicons. Complementary sequences were generated from the published sequences of primers 519 and 1114. Degeneracies were added to primer 515R following comparison to the oral bacterial database to maximize matches of primer against bacterial sequences.

[Figure omitted. See PDF.]

Table 1. Sequences of primers used in study.

https://doi.org/10.1371/journal.pone.0020956.t001

Pyrosequencing

Multiplexed bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP) was performed using the Titanium platform (Roche Applied Science, Indianapolis, IN) as previously described [22] in a commercial facility (Research and Testing Laboratories, Lubbock, TX). Briefly, a single step PCR with broad-range universal primers and 22 cycles of amplification was used to amplify the 16S rRNA genes as well as to introduce adaptor sequences and sample-specific 10-mer oligonucleotide tags into the DNA. The same bar codes were utilized for each primer set. Three regions of the 16S gene were sequenced from each sample (V1–V3, V4–V6, V7–V9). Adaptor sequences were trimmed from raw data with 98% or more of bases demonstrating a quality control of 30 and sequences binned into individual sample collections based on bar-code sequence tags, which were then trimmed. The resulting files were denoised with Pyronoise [23] and depleted of chimeras using B2C2 (http://www.researchandtesting.com/B2C2.html). Sequences less than <300 bp in length were deleted and the rest were clustered into species-level operational taxonomic units (s-OTUs) at 97% sequence similarity and assigned a taxonomic identity by alignment to locally hosted version of the Greengenes database [24] using the Blastn algorithm. Phylogenetic trees were generated and visualized using FastTree [25]. All analyses were conducted within the virtual environment provided by the QIIME pipeline [26].

Statistical analysis

Species-level OTUs (s-OTUs) were used to compute the Shannon Diversity and Equitability indices for each sample. EstimateS ((Version 7.5, R. K. Colwell, http://purl.oclc.org/estimates) was used to compute the indices and statistical analyses were carried out with JMP (SAS Institute Inc., Cary, NC). The indices were compared between groups using ANOVA. A variance stabilizing transformation was used to create normal distribution of the data as previously described [27], [28]. Two sample t-tests were used to compare the transformed values of species and genus-level OTUs between groups. Fisher's exact test was used to test for presence or absence of genera.

Results

The pyrotag sequences were compared to previously published data obtained by Sanger sequencing using the primer pairs A17 and 317 on the same samples [27]. A subset of the pyrosequencing data was created using a random number generator to select 100 pyrotag sequences from each primer set. This subset was compared to an equivalent number of Sanger sequences. A total of 1054 nearly full-length sequences (1300–1460 bp) were identified by Sanger sequencing, and 167,210 sequences by pyrosequencing, representing a 167-fold increase in depth-of-coverage with pyrosequencing.

Figure 1 shows the Shannon Diversity and Equitability indices for all primer sets. The Diversity Index was not different between groups; however, the community generated by Sanger sequencing demonstrated significantly greater equitability than all the pyrotag communities (p<0.01, ANOVA). Pyrotag communities generated by the 4 primer pairs demonstrated similar diversity.

[Figure omitted. See PDF.]

Figure 1. Shannon diversity and equitability indices of pyrotag and Sanger communities.

No differences were detected between any of the pyrotag communities; however, the Sanger community demonstrated significantly greater equitable than all the pyrotag communities (** p<0.01, ANOVA).

https://doi.org/10.1371/journal.pone.0020956.g001

Figures 2A and 2B show the distribution of rare and abundant taxa by primer pair and sequencing methodology. 1.9% of sequences could not be classified into any taxon below the level of domain. Taxa with less than 20 overall sequences were designated as rare. Sanger sequences demonstrated significantly lower coverage of rare as well as abundant species than pyrosequencing (p<0.001, ANOVA). However, there were no differences in the number of rare and abundant taxa in any of the pyrotag communities.

[Figure omitted. See PDF.]

Figure 2. Distribution of sequences by taxa.

Rare taxa are shown in Figure 2A and abundant taxa in Figure 2B. The Sanger community demonstrated significantly fewer species-level taxa than pyrosequencing (*** p<0.001, ANOVA). There were no differences between the pyrotag sequences.

https://doi.org/10.1371/journal.pone.0020956.g002

Figure 3 shows the distribution by genus of sequences generated by degenerate and non-degenerate primer pairs targeted to the V1–V3 region. There were no differences between the two groups (p>0.05, 2-sample t-test on transformed variable).

[Figure omitted. See PDF.]

Figure 3. Distribution of sequences generated by degenerate and non-degenerate primers by genus.

Percent mean abundances and standard deviations are shown. Genera are arranged in a gradient such that those predominant in the degenerate community are arranged on the left. There were no differences between the two communities in the relative abundance of any genus (p>0.05, 2-sample t-test on transformed variable).

https://doi.org/10.1371/journal.pone.0020956.g003

Table 2 shows the relative abundance of genera in sequences obtained by pyrosequencing different target regions. Genera accounting for 0.1% of total pyrosequences are shown. Overall, greater numbers of differences were detected in the levels of genera between the V1–V3 and V7–V9 regions (p<0.05, 2-sample t-test on transformed variable). The regions targeted significantly influenced community profiles generated by pyrosequencing. The genera Prevotella, Fusobacterium, Streptococcus, Granulicatella, Bacteroides, Porphyromonas and Treponema formed 65% of the community when the V1–V3 region was targeted, while Streptococcus, Treponema, Prevotella, Eubacterium, Porphyromonas, Campylobacter and Enterococcus accounted for the same abundance in the community generated by V4–V6 primers, and 65% of the V7–V9 community was formed by Veillonella, Streptococcus, Eubacterium, Enterococcus, Treponema, Catonella and Selenomonas. Among the predominant genera, Fusobacteria were not detected in any of the samples by the V7–V9 primers, while the V4–V6 primers did not detect the Selenomonads, Mycoplasma, or TM7 phylum in any sample (p<0.05, Fisher's exact test).

[Figure omitted. See PDF.]

Table 2. Relative abundances of genera in pyrotag sequences.

https://doi.org/10.1371/journal.pone.0020956.t002

Table 3 shows the relative abundance of genera obtained by concatenating data from pairs of target regions or by combining all three regions to provide near-full-length coverage of the 16S gene. Relative abundances of the same genera in near-full-length Sanger sequences are also shown for comparison. To arrive at these results, the subset pyrotag dataset was compared to an equivalent number of Sanger sequences from each sample. Concatenating data from V1–V3 and V7–V9 regions demonstrated the greatest similarity to Sanger data as well as to the averages of all 3 regions.

[Figure omitted. See PDF.]

Table 3. Relative abundances of genera in Sanger and concatenated pyrotag datasets.

https://doi.org/10.1371/journal.pone.0020956.t003

Discussion

It has been shown that sequences of 500–700 bp are required for phylogenetic discrimination at the species levels [9], [29]. However, previous reports have been equivocal on the level of community coverage achieved using the different hypervariable regions. While several investigations support using the V1, V2 and V3 regions for deep sequencing [17], others suggest that these regions overestimate species richness and promote the V4–V6 region as the most appropriate [19]. Yet others have demonstrated that V7–V8 fragments achieve representational characterization of a community [30]. Our previous investigations with Sanger sequencing have revealed that the subgingival microflora associated with periodontitis in smokers is extremely diverse, with several rare species/phylotypes [27]. Hence, plaque samples were collected and pooled from deep sites of current smokers with moderate to severe periodontitis to examine the extent to which primer design affects the community fingerprint of a highly complex and taxonomically heterogeneous microbial population. Using an adequately powered clinical study design to enable statistical analyses allowed an in-depth comparison of the community profiles generated by the different primer sets.

The Shannon Diversity index incorporates both the number of species (species richness) as well as the proportion of each species (species evenness) into a single value [31]. Thus, while a value of zero necessarily represents a mono-species community, a higher value may result either from the presence of several species at varying levels or from equitable distribution of a few species. Hence, the Equitability index is used to elucidate the relative contributions of species richness and evenness to the Diversity index. Pyrotag communities demonstrated similar diversity to the Sanger community, however, were significantly less equitable (Figure 1), suggesting that greater species richness contributed to the diversity. The increased species richness was apparent in both rare and abundant taxa (Figure 2). This is in contrast to previous investigations; which have suggested that pyrosequencing overestimates community diversity by overestimating the number of rare taxa [10], [32]. A single-step PCR with low cycle numbers and a high fidelity, proofreading polymerase were utilized in this study; and it is possible that this minimized over-representation of rare taxa in the present investigation. No differences were apparent in the number of rare and abundant taxa between the different hypervariable regions; suggesting that targeting a specific region for pyrosequencing does not affect species richness. Taken together, it appears that selecting a specific region for pyrosequencing is not a source of bias in the diversity of the resulting community or in the number of taxa detected.

Out of the four primer pairs selected, two pairs targeted the same region (V1–V3), one pair containing degenerate sequences and the other non-degenerate. Fragments encompassing the V1–V3 region have been the most common targets for both Sanger sequencing and pyrosequencing; and both non-degenerate and degenerate primers have been used to amplify this region [2], [6], [7], [8], [9], [33]. It has previously been suggested that inclusion of degenerate sequences improves the “universality” of primers (reviewed by Baker et al [15]), however, our data does not support a role for primer degeneracy in improving community coverage. This is in concordance with previous investigations that have reported no effect of primer degeneracy on profiles of naturally occurring microbial communities [34]. Although degenerate primers, by virtue of their lowered specificity, may amplify larger number of taxa within a community, it has been shown that this effect is magnified when large PCR cycle numbers are used [35]. The present investigation used 22 cycles to amplification to ensure representational amplification of the community template, and it is possible that the low cycle numbers precluded a possible influence by degenerate primers.

Our data suggest that the hypervariable region targeted for sequencing plays a critical role in influencing the composition of pyrotag communities. Previous investigations have reported that amplicon size and PCR kinetics may be a source of sequencing bias [36], [37]. To overcome this in the present study, sequencing primers were carefully selected to generate similar amplicon sizes (∼500 bp for V1–V3 amplicons, ∼550 bp for V4–V6 amplicons and ∼470 bp for V7–V9 amplicons). Identical PCR cycling conditions were also utilized for all primer sets, thereby reducing the possibility of bias from this source. Using a single pyrosequencing run to generate all sequences further reduced bias due to PCR and sequencing kinetics. Thus, the observed differences could not be attributed to these variables. It is especially striking that even though these samples were derived from sites with severe disease, the V7–V9 communities were dominated by Veillonella and the V4–V6 communities by Streptococci (Table 2), genera that have been previously associated with periodontal health [33]. Similarly, Treponema, a disease-associated genus; was found in high numbers in the V4–V6 and V7–V9 communities; while other disease-associated genera, for example, Prevotella, Porphyromonas, and Bacteroides were predominant in V1–V3 communities derived from the same run. Fusobacteria were undetected by the V7–V9 primers while forming nearly 19% of the V1–V3 community. Similarly, the Selenomonads were not detectable by the V4–V6 primers, while forming 6% of the V7–V9 community. Concatenated data from V1–V3 and V7–V9 regions resulted in community profiles that did not significantly differ from Sanger sequences or full-length pyrosequences for the predominant genera, while averages of the other two regions did not yield similar results (Table 3). It is also noteworthy that the greatest differences were observed in the community fingerprints generated by these two primer sets. The mechanism causing this difference is not clear and warrants further investigation. It could be hypothesized that presence and nature of secondary structures within the target regions as well as the GC ratios of the resultant fragments may have contributed to the differences. It is known that the V1, V4 and V7 regions exhibit differences in the number of stems as well as in nucleotide variations within these stems [38], [39], and while is possible that differential amplification efficiencies contributed to the compositional differences, it is not within the scope of this study to test this hypothesis. It has been shown that higher GC ratios result in higher amplification efficiencies [35], thereby altering PCR kinetics, with over-amplification of rare members and under-representation of dominant species [40]. In the present investigation, however, the GC ratios of the different amplicons were very similar; therefore, the observed discrepancies could not be attributable to this variable.

In summary, the hypervariable region targeted by the primer plays a critical role in determining the profile of a largely uncultivated, complex microbial community generated by pyrosequencing. This effect is significant, with the presence of certain dominant community members being masked and others being under-represented with different primer sets; thereby providing a critical source of error in microbial ecological studies. However, averaging the community fingerprints generated by V1–V3 and V7–V9 primers provides 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 characterize heterogeneous microbial communities.

Author Contributions

Conceived and designed the experiments: PSK. Performed the experiments: MRB SED. Analyzed the data: TC PSK MRB. Wrote the paper: PSK MRB.

Citation: Kumar PS, Brooker MR, Dowd SE, Camerlengo T (2011) Target Region Selection Is a Critical Determinant of Community Fingerprints Generated by 16S Pyrosequencing. PLoS ONE 6(6): e20956. https://doi.org/10.1371/journal.pone.0020956

References

1. Diaz PI, Chalmers NI, Rickard AH, Kong C, Milburn CL, et al. (2006) Molecular characterization of subject-specific oral microflora during initial colonization of enamel. Appl Environ Microbiol 72: 2837–2848.PI DiazNI ChalmersAH RickardC. KongCL Milburn2006Molecular characterization of subject-specific oral microflora during initial colonization of enamel.Appl Environ Microbiol7228372848

2. de Lillo A, Ashley FP, Palmer RM, Munson MA, Kyriacou L, et al. (2006) Novel subgingival bacterial phylotypes detected using multiple universal polymerase chain reaction primer sets. Oral Microbiol Immunol 21: 61–68.A. de LilloFP AshleyRM PalmerMA MunsonL. Kyriacou2006Novel subgingival bacterial phylotypes detected using multiple universal polymerase chain reaction primer sets.Oral Microbiol Immunol216168

3. Delima SL, McBride RK, Preshaw PM, Heasman PA, Kumar PS (2010) Response of subgingival bacteria to smoking cessation. J Clin Microbiol 48: 2344–2349.SL DelimaRK McBridePM PreshawPA HeasmanPS Kumar2010Response of subgingival bacteria to smoking cessation.J Clin Microbiol4823442349

4. Gomes SC, Piccinin FB, Oppermann RV, Susin C, Nonnenmacher CI, et al. (2006) Periodontal status in smokers and never-smokers: clinical findings and real-time polymerase chain reaction quantification of putative periodontal pathogens. J Periodontol 77: 1483–1490.SC GomesFB PiccininRV OppermannC. SusinCI Nonnenmacher2006Periodontal status in smokers and never-smokers: clinical findings and real-time polymerase chain reaction quantification of putative periodontal pathogens.J Periodontol7714831490

5. Aas JA, Paster BJ, Stokes LN, Olsen I, Dewhirst FE (2005) Defining the normal bacterial flora of the oral cavity. J Clin Microbiol 43: 5721–5732.JA AasBJ PasterLN StokesI. OlsenFE Dewhirst2005Defining the normal bacterial flora of the oral cavity.J Clin Microbiol4357215732

6. Hutter G, Schlagenhauf U, Valenza G, Horn M, Burgemeister S, et al. (2003) Molecular analysis of bacteria in periodontitis: evaluation of clone libraries, novel phylotypes and putative pathogens. Microbiology 149: 67–75.G. HutterU. SchlagenhaufG. ValenzaM. HornS. Burgemeister2003Molecular analysis of bacteria in periodontitis: evaluation of clone libraries, novel phylotypes and putative pathogens.Microbiology1496775

7. Kroes I, Lepp PW, Relman DA (1999) Bacterial diversity within the human subgingival crevice. Proc Natl Acad Sci U S A 96: 14547–14552.I. KroesPW LeppDA Relman1999Bacterial diversity within the human subgingival crevice.Proc Natl Acad Sci U S A961454714552

8. Kumar PS, Griffen AL, Moeschberger ML, Leys EJ (2005) Identification of candidate periodontal pathogens and beneficial species by quantitative 16S clonal analysis. J Clin Microbiol 43: 3944–3955.PS KumarAL GriffenML MoeschbergerEJ Leys2005Identification of candidate periodontal pathogens and beneficial species by quantitative 16S clonal analysis.J Clin Microbiol4339443955

9. Paster BJ, Boches SK, Galvin JL, Ericson RE, Lau CN, et al. (2001) Bacterial diversity in human subgingival plaque. J Bacteriol 183: 3770–3783.BJ PasterSK BochesJL GalvinRE EricsonCN Lau2001Bacterial diversity in human subgingival plaque.J Bacteriol18337703783

10. Kunin V, Engelbrektson A, Ochman H, Hugenholtz P. Wrinkles in the rare biosphere: pyrosequencing errors can lead to artificial inflation of diversity estimates. Environ Microbiol 12: 118–123.V. KuninA. EngelbrektsonH. OchmanP. HugenholtzWrinkles in the rare biosphere: pyrosequencing errors can lead to artificial inflation of diversity estimates.Environ Microbiol12118123

11. Zaura E, Keijser BJ, Huse SM, Crielaard W (2009) Defining the healthy “core microbiome” of oral microbial communities. BMC Microbiol 9: 259.E. ZauraBJ KeijserSM HuseW. Crielaard2009Defining the healthy “core microbiome” of oral microbial communities.BMC Microbiol9259

12. Keijser BJ, Zaura E, Huse SM, van der Vossen JM, Schuren FH, et al. (2008) Pyrosequencing analysis of the oral microflora of healthy adults. J Dent Res 87: 1016–1020.BJ KeijserE. ZauraSM HuseJM van der VossenFH Schuren2008Pyrosequencing analysis of the oral microflora of healthy adults.J Dent Res8710161020

13. Li L, Hsiao WW, Nandakumar R, Barbuto SM, Mongodin EF, et al. Analyzing Endodontic Infections by Deep Coverage Pyrosequencing. J Dent Res. L. LiWW HsiaoR. NandakumarSM BarbutoEF MongodinAnalyzing Endodontic Infections by Deep Coverage Pyrosequencing.J Dent Res

14. Dominguez-Bello MG, Costello EK, Contreras M, Magris M, Hidalgo G, et al. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci U S A 107: 11971–11975.MG Dominguez-BelloEK CostelloM. ContrerasM. MagrisG. HidalgoDelivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns.Proc Natl Acad Sci U S A1071197111975

15. Baker GC, Smith JJ, Cowan DA (2003) Review and re-analysis of domain-specific 16S primers. J Microbiol Methods 55: 541–555.GC BakerJJ SmithDA Cowan2003Review and re-analysis of domain-specific 16S primers.J Microbiol Methods55541555

16. Wang Y, Qian PY (2009) Conservative fragments in bacterial 16S rRNA genes and primer design for 16S ribosomal DNA amplicons in metagenomic studies. PLoS One 4: e7401.Y. WangPY Qian2009Conservative fragments in bacterial 16S rRNA genes and primer design for 16S ribosomal DNA amplicons in metagenomic studies.PLoS One4e7401

17. Chakravorty S, Helb D, Burday M, Connell N, Alland D (2007) A detailed analysis of 16S ribosomal RNA gene segments for the diagnosis of pathogenic bacteria. J Microbiol Methods 69: 330–339.S. ChakravortyD. HelbM. BurdayN. ConnellD. Alland2007A detailed analysis of 16S ribosomal RNA gene segments for the diagnosis of pathogenic bacteria.J Microbiol Methods69330339

18. Nossa CW, Oberdorf WE, Yang L, Aas JA, Paster BJ, et al. rRNA gene primers for 454 pyrosequencing of the human foregut microbiome. World J Gastroenterol 16: 4135–4144.CW NossaWE OberdorfL. YangJA AasBJ PasterrRNA gene primers for 454 pyrosequencing of the human foregut microbiome.World J Gastroenterol1641354144

19. Youssef N, Sheik CS, Krumholz LR, Najar FZ, Roe BA, et al. (2009) Comparison of species richness estimates obtained using nearly complete fragments and simulated pyrosequencing-generated fragments in 16S rRNA gene-based environmental surveys. Appl Environ Microbiol 75: 5227–5236.N. YoussefCS SheikLR KrumholzFZ NajarBA Roe2009Comparison of species richness estimates obtained using nearly complete fragments and simulated pyrosequencing-generated fragments in 16S rRNA gene-based environmental surveys.Appl Environ Microbiol7552275236

20. Claesson MJ, Wang Q, O'Sullivan O, Greene-Diniz R, Cole JR, et al. Comparison of two next-generation sequencing technologies for resolving highly complex microbiota composition using tandem variable 16S rRNA gene regions. Nucleic Acids Res 38: e200.MJ ClaessonQ. WangO. O'SullivanR. Greene-DinizJR ColeComparison of two next-generation sequencing technologies for resolving highly complex microbiota composition using tandem variable 16S rRNA gene regions.Nucleic Acids Res38e200

21. Riggio MP, Lennon A, Rolph HJ, Hodge PJ, Donaldson A, et al. (2008) Molecular identification of bacteria on the tongue dorsum of subjects with and without halitosis. Oral Dis 14: 251–258.MP RiggioA. LennonHJ RolphPJ HodgeA. Donaldson2008Molecular identification of bacteria on the tongue dorsum of subjects with and without halitosis.Oral Dis14251258

22. Dowd SE, Sun Y, Wolcott RD, Domingo A, Carroll JA (2008) Bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP) for microbiome studies: bacterial diversity in the ileum of newly weaned Salmonella-infected pigs. Foodborne Pathog Dis 5: 459–472.SE DowdY. SunRD WolcottA. DomingoJA Carroll2008Bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP) for microbiome studies: bacterial diversity in the ileum of newly weaned Salmonella-infected pigs.Foodborne Pathog Dis5459472

23. Quince C, Lanzen A, Davenport RJ, Turnbaugh PJ.Removing noise from pyrosequenced amplicons. BMC Bioinformatics 12: 38.C. QuinceA. LanzenRJ DavenportPJ. TurnbaughRemoving noise from pyrosequenced amplicons.BMC Bioinformatics1238

24. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, et al. (2006) Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol 72: 5069–5072.TZ DeSantisP. HugenholtzN. LarsenM. RojasEL Brodie2006Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB.Appl Environ Microbiol7250695072

25. Price MN, Dehal PS, Arkin AP (2009) FastTree: Computing Large Minimum Evolution Trees with Profiles instead of a Distance Matrix. Molecular Biology and Evolution 26: 1641–1650.MN PricePS DehalAP Arkin2009FastTree: Computing Large Minimum Evolution Trees with Profiles instead of a Distance Matrix.Molecular Biology and Evolution2616411650%R 1610.1093/molbev/msp1077. %R 1610.1093/molbev/msp1077.

26. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7: 335–336.JG CaporasoJ. KuczynskiJ. StombaughK. BittingerFD BushmanQIIME allows analysis of high-throughput community sequencing data.Nat Methods7335336

27. Shchipkova AY, Nagaraja HN, Kumar PS (2010) Subgingival Microbial Profiles of Smokers with Periodontitis. J Dent Res. AY ShchipkovaHN NagarajaPS Kumar2010Subgingival Microbial Profiles of Smokers with Periodontitis.J Dent Res

28. Osborne JW (2002) Notes on the use of data transformations. Practical Assessment, Research and Evaluation. JW Osborne2002Notes on the use of data transformations.Practical Assessment, Research and Evaluation

29. Clarridge JE 3rd (2004) Impact of 16S rRNA gene sequence analysis for identification of bacteria on clinical microbiology and infectious diseases. Clin Microbiol Rev 17: 840–862.JE Clarridge 3rd2004Impact of 16S rRNA gene sequence analysis for identification of bacteria on clinical microbiology and infectious diseases.Clin Microbiol Rev17840862table of contents. table of contents.

30. Trotha R, Hanck T, Konig W, Konig B (2001) Rapid ribosequencing–an effective diagnostic tool for detecting microbial infection. Infection 29: 12–16.R. TrothaT. HanckW. KonigB. Konig2001Rapid ribosequencing–an effective diagnostic tool for detecting microbial infection.Infection291216

31. Shannon CE (1997) The mathematical theory of communication. 1963. MD Comput 14: 306–317.CE Shannon1997The mathematical theory of communication. 1963.MD Comput14306317

32. Huse SM, Huber JA, Morrison HG, Sogin ML, Welch DM (2007) Accuracy and quality of massively parallel DNA pyrosequencing. Genome Biol 8: R143.SM HuseJA HuberHG MorrisonML SoginDM Welch2007Accuracy and quality of massively parallel DNA pyrosequencing.Genome Biol8R143

33. Kumar PS, Leys EJ, Bryk JM, Martinez FJ, Moeschberger ML, et al. (2006) Changes in periodontal health status are associated with bacterial community shifts as assessed by quantitative 16S cloning and sequencing. J Clin Microbiol 44: 3665–3673.PS KumarEJ LeysJM BrykFJ MartinezML Moeschberger2006Changes in periodontal health status are associated with bacterial community shifts as assessed by quantitative 16S cloning and sequencing.J Clin Microbiol4436653673

34. Huber JA, Morrison HG, Huse SM, Neal PR, Sogin ML, et al. (2009) Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure. Environ Microbiol 11: 1292–1302.JA HuberHG MorrisonSM HusePR NealML Sogin2009Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure.Environ Microbiol1112921302

35. Polz MF, Cavanaugh CM (1998) Bias in template-to-product ratios in multitemplate PCR. Appl Environ Microbiol 64: 3724–3730.MF PolzCM Cavanaugh1998Bias in template-to-product ratios in multitemplate PCR.Appl Environ Microbiol6437243730

36. Kleter B, van Doorn LJ, ter Schegget J, Schrauwen L, van Krimpen K, et al. (1998) Novel short-fragment PCR assay for highly sensitive broad-spectrum detection of anogenital human papillomaviruses. Am J Pathol 153: 1731–1739.B. KleterLJ van DoornJ. ter ScheggetL. SchrauwenK. van Krimpen1998Novel short-fragment PCR assay for highly sensitive broad-spectrum detection of anogenital human papillomaviruses.Am J Pathol15317311739

37. Suzuki MT, Giovannoni SJ (1996) Bias caused by template annealing in the amplification of mixtures of 16S rRNA genes by PCR. Appl Environ Microbiol 62: 625–630.MT SuzukiSJ Giovannoni1996Bias caused by template annealing in the amplification of mixtures of 16S rRNA genes by PCR.Appl Environ Microbiol62625630

38. Anton AI, Martinez-Murcia AJ, Rodriguez-Valera F, Dalet F (2001) Sequence microdiversity at the ribosomal RNA operons of Escherichia coli pyelonephritogenic strains. Clin Microbiol Infect 7: 345–351.AI AntonAJ Martinez-MurciaF. Rodriguez-ValeraF. Dalet2001Sequence microdiversity at the ribosomal RNA operons of Escherichia coli pyelonephritogenic strains.Clin Microbiol Infect7345351

39. Martinez-Murcia AJ, Anton AI, Rodriguez-Valera F (1999) Patterns of sequence variation in two regions of the 16S rRNA multigene family of Escherichia coli. Int J Syst Bacteriol 49 Pt 2: 601–610.AJ Martinez-MurciaAI AntonF. Rodriguez-Valera1999Patterns of sequence variation in two regions of the 16S rRNA multigene family of Escherichia coli.Int J Syst Bacteriol49 Pt 2601610

40. Suzuki M, Rappe MS, Giovannoni SJ (1998) Kinetic bias in estimates of coastal picoplankton community structure obtained by measurements of small-subunit rRNA gene PCR amplicon length heterogeneity. Appl Environ Microbiol 64: 4522–4529.M. SuzukiMS RappeSJ Giovannoni1998Kinetic bias in estimates of coastal picoplankton community structure obtained by measurements of small-subunit rRNA gene PCR amplicon length heterogeneity.Appl Environ Microbiol6445224529