The effects of alignment quality, distance calculation method, sequence filtering, and region on the analysis of 16S rRNA gene-based studies

• Published in: PLoS computational biology Vol. 6; no. 7; p. e1000844
• Main Author: Schloss, Patrick D
• Format: Journal Article
• Language: English
• Published: United States Public Library of Science 08.07.2010; Public Library of Science (PLoS)
• Subjects: Algorithms; Community; Computational Biology - methods; Computational Biology/Comparative Sequence Analysis; Computational Biology/Genomics; Databases, Genetic; DNA sequencing; Escherichia coli - genetics; Genes, Bacterial; Genetic diversity; Genetic Variation; Methods; Microbiology; Microbiology/Environmental Microbiology; Nucleotide sequencing; Phylogenetics; Phylogeny; Physiological aspects; Polymerase chain reaction; Regression Analysis; Ribosomal RNA; RNA, Ribosomal, 16S - genetics; Sequence Alignment; Studies; Trees; United States; Genetic Variation; Genes, Bacterial; Regression Analysis; Sequence Alignment; Escherichia coli; Computational Biology; Databases, Genetic; Phylogeny; RNA, Ribosomal, 16S
• ISSN: 1553-7358; 1553-734X; 1553-7358
• DOI: 10.1371/journal.pcbi.1000844

Pyrosequencing of PCR-amplified fragments that target variable regions within the 16S rRNA gene has quickly become a powerful method for analyzing the membership and structure of microbial communities. This approach has revealed and introduced questions that were not fully appreciated by those carrying out traditional Sanger sequencing-based methods. These include the effects of alignment quality, the best method of calculating pairwise genetic distances for 16S rRNA genes, whether it is appropriate to filter variable regions, and how the choice of variable region relates to the genetic diversity observed in full-length sequences. I used a diverse collection of 13,501 high-quality full-length sequences to assess each of these questions. First, alignment quality had a significant impact on distance values and downstream analyses. Specifically, the greengenes alignment, which does a poor job of aligning variable regions, predicted higher genetic diversity, richness, and phylogenetic diversity than the SILVA and RDP-based alignments. Second, the effect of different gap treatments in determining pairwise genetic distances was strongly affected by the variation in sequence length for a region; however, the effect of different calculation methods was subtle when determining the sample's richness or phylogenetic diversity for a region. Third, applying a sequence mask to remove variable positions had a profound impact on genetic distances by muting the observed richness and phylogenetic diversity. Finally, the genetic distances calculated for each of the variable regions did a poor job of correlating with the full-length gene. Thus, while it is tempting to apply traditional cutoff levels derived for full-length sequences to these shorter sequences, it is not advisable. Analysis of beta-diversity metrics showed that each of these factors can have a significant impact on the comparison of community membership and structure. Taken together, these results urge caution in the design and interpretation of analyses using pyrosequencing data.