Pipeline for amplifying and analyzing amplicons of the V1-V3 region of the 16S rRNA gene

被引：54

作者：

Allen H.K. ^{[1
]}

Bayles D.O. ^{[2
]}

Looft T. ^{[1
]}

Trachsel J. ^{[1
,3
]}

Bass B.E. ^{[1
,4
]}

Alt D.P. ^{[2
]}

Bearson S.M.D. ^{[1
]}

Nicholson T. ^{[5
]}

Casey T.A. ^{[1
]}

机构：

[1] Food Safety and Enteric Pathogens Research Unit, National Animal Disease Center, Agricultural Research Service, Ames, 50010, IA

[2] Infectious Bacterial Diseases Research Unit, National Animal Disease Center, Agricultural Research Service, Ames, 50010, IA

[3] Interdepartmental Microbiology Graduate Program, Iowa State University, Ames, 50010, IA

[4] Diamond v, Cedar Rapids, 52404, IA

[5] Virus and Prion Research Unit, National Animal Disease Center, Agricultural Research Service, USDA, Ames, 50010, IA

来源：

BMC Research Notes | / 9卷 / 1期

关键词：

16S rRNA gene; Microbial ecology; MiSeq; Mock community; V1-V3;

D O I：

10.1186/s13104-016-2172-6

中图分类号：

学科分类号：

摘要：

Background: Profiling of 16S rRNA gene sequences is an important tool for testing hypotheses in complex microbial communities, and analysis methods must be updated and validated as sequencing technologies advance. In host-associated bacterial communities, the V1-V3 region of the 16S rRNA gene is a valuable region to profile because it provides a useful level of taxonomic resolution; however, use of Illumina MiSeq data for experiments targeting this region needs validation. Results: Using a MiSeq machine and the version 3 (300 × 2) chemistry, we sequenced the V1-V3 region of the 16S rRNA gene within a mock community. Nineteen bacteria and one archaeon comprised the mock community, and 12 replicate amplifications of the community were performed and sequenced. Sequencing the large fragment (490 bp) that encompasses V1-V3 yielded a higher error rate (3.6 %) than has been reported when using smaller fragment sizes. This higher error rate was due to a large number of sequences that occurred only one or two times among all mock community samples. Removing sequences that occurred one time among all samples (singletons) reduced the error rate to 1.4 %. Diversity estimates of the mock community containing all sequences were inflated, whereas estimates following singleton removal more closely reflected the actual mock community membership. A higher percentage of the sequences could be taxonomically assigned after singleton and doubleton sequences were removed, and the assignments reflected the membership of the input DNA. Conclusions: Sequencing the V1-V3 region of the 16S rRNA gene on the MiSeq platform may require additional sequence curation in silico, and improved error rates and diversity estimates show that removing low-frequency sequences is reasonable. When datasets have a high number of singletons, these singletons can be removed from the analysis without losing statistical power while reducing error and improving microbiota assessment. © 2016 The Author(s).

引用

共 30 条

[1]

Evaluation of 16S rDNA-based community profiling for human microbiome research, PLoS One, 7, 6, (2012)

[2]

Kozich J.J., Westcott S.L., Baxter N.T., Highlander S.K., Schloss P.D., Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform, Appl Environ Microbiol, 79, 17, pp. 5112-5120, (2013)

[3]

Edgar R.C., Haas B.J., Clemente J.C., Quince C., Knight R., UCHIME improves sensitivity and speed of chimera detection, Bioinformatics, 27, 16, pp. 2194-2200, (2011)

[4]

Haas B.J., Gevers D., Earl A.M., Feldgarden M., Ward D.V., Giannoukos G., Ciulla D., Tabbaa D., Highlander S.K., Sodergren E., Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons, Genome Res, 21, 3, pp. 494-504, (2011)

[5]

Schloss P.D., Gevers D., Westcott S.L., Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies, PLoS One, 6, 12, (2011)

[6]

Schloss P.D., Westcott S.L., Assessing and improving methods used in operational taxonomic unit-based approaches for 16S rRNA gene sequence analysis, Appl Environ Microbiol, 77, 10, pp. 3219-3226, (2011)

[7]

May A., Abeln S., Crielaard W., Heringa J., Brandt B.W., Unraveling the outcome of 16S rDNA-based taxonomy analysis through mock data and simulations, Bioinformatics, 30, 11, pp. 1530-1538, (2014)

[8]

Shade A., Caporaso J.G., Handelsman J., Knight R., Fierer N., A meta-analysis of changes in bacterial and archaeal communities with time, ISME J, 7, 8, pp. 1493-1506, (2013)

[9]

Conlan S., Kong H.H., Segre J.A., Species-level analysis of DNA sequence data from the NIH human microbiome project, PLoS One, 7, 10, (2012)

[10]

Looft T., Allen H.K., Casey T.A., Alt D.P., Stanton T.B., Carbadox has both temporary and lasting effects on the swine gut microbiota, Front Microbiol, 5, (2014)

← 1 2 3 →