Improved Bacterial 16S rRNA Gene (V4 and V4-5) and Fungal Internal Transcribed Spacer Marker Gene Primers for Microbial Community Surveys

William Walters, Embriette R Hyde, Donna Berg-Lyons, Gail Ackermann, Greg Humphrey, Alma Parada, Jack A Gilbert, Janet K Jansson, J Gregory Caporaso, Jed A Fuhrman, Amy Apprill, Rob Knight, William Walters, Embriette R Hyde, Donna Berg-Lyons, Gail Ackermann, Greg Humphrey, Alma Parada, Jack A Gilbert, Janet K Jansson, J Gregory Caporaso, Jed A Fuhrman, Amy Apprill, Rob Knight

Abstract

Designing primers for PCR-based taxonomic surveys that amplify a broad range of phylotypes in varied community samples is a difficult challenge, and the comparability of data sets amplified with varied primers requires attention. Here, we examined the performance of modified 16S rRNA gene and internal transcribed spacer (ITS) primers for archaea/bacteria and fungi, respectively, with nonaquatic samples. We moved primer bar codes to the 5' end, allowing for a range of different 3' primer pairings, such as the 515f/926r primer pair, which amplifies variable regions 4 and 5 of the 16S rRNA gene. We additionally demonstrated that modifications to the 515f/806r (variable region 4) 16S primer pair, which improves detection of Thaumarchaeota and clade SAR11 in marine samples, do not degrade performance on taxa already amplified effectively by the original primer set. Alterations to the fungal ITS primers did result in differential but overall improved performance compared to the original primers. In both cases, the improved primers should be widely adopted for amplicon studies. IMPORTANCE We continue to uncover a wealth of information connecting microbes in important ways to human and environmental ecology. As our scientific knowledge and technical abilities improve, the tools used for microbiome surveys can be modified to improve the accuracy of our techniques, ensuring that we can continue to identify groundbreaking connections between microbes and the ecosystems they populate, from ice caps to the human body. It is important to confirm that modifications to these tools do not cause new, detrimental biases that would inhibit the field rather than continue to move it forward. We therefore demonstrated that two recently modified primer pairs that target taxonomically discriminatory regions of bacterial and fungal genomic DNA do not introduce new biases when used on a variety of sample types, from soil to human skin. This confirms the utility of these primers for maintaining currently recommended microbiome research techniques as the state of the art.

Keywords: 16S; ITS; marker genes; microbial ecology; primers.

Figures

FIG 1
FIG 1
Comparison of the original 515f/806r primer pair and the new, modified 515f/806rB primer pair. (A) Procrustes plot of original and modified 515f/806r constructs, with unweighted UniFrac metric, M2 = 0.111. (B) Original and modified 515f/806r constructs, weighted UniFrac metric M2 = 0.196. (C) Pie charts illustrating the mean relative abundance of phyla present (all studies combined) in samples amplified with the old 515f/806r construct or with the modified 515f/806rB construct. (D) Taxa scatterplots for the original and modified 515f/806r primers. Phyla plots are shown for American Gut fecal, American Gut skin, agricultural soils, EMP Rice Rhizome, Body Farm 1, Body Farm 2, mouse decomposition, and Sloan built environment samples. Outlier samples have been removed from the results shown.
FIG 2
FIG 2
Comparison of the modified 515f/806rB V4 primer pair and the 515f (modified)/926r primer pair. (A) Procrustes plot of modified 515f/806r and 515f/926r constructs, Bray-Curtis dissimilarity M2 = 0.058. (B) Pie charts illustrating the mean relative abundance of phyla present (all studies combined) in samples amplified with the modified 515f/806rB construct or with the 515f/926r construct. (C) Taxa scatterplots for the modified 515f/806rB construct and the 515f/926r construct. Phyla plots are shown for American Gut fecal, American Gut skin, agricultural soils, EMP Rice Rhizome, Body Farm 1, Body Farm 2, mouse decomposition, and Sloan built environment samples. Outlier samples have not been removed from the data shown.
FIG 3
FIG 3
Comparison of the original ITS primer pair and the new, modified ITS primer pair. (A) Procrustes plot of original and modified ITS1 constructs, Bray-Curtis dissimilarity M2 = 0.363. (B) Pie charts illustrating the mean relative abundance of phyla present (all studies combined) in samples amplified with the old ITS construct and the modified ITS construct. (C) Taxa scatterplots for the original and modified ITS primers. Phyla plots are shown for American Gut fecal, American Gut skin, agricultural soils, EMP Rice Rhizome, Body Farm 1, Body Farm 2, mouse decomposition, and Sloan built environment samples. Outlier samples have been removed from the data shown.

References

    1. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R. 2011. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A 108(Suppl):4516–4522. doi:10.1073/pnas.1000080107.
    1. Hugerth LW, Wefer HA, Lundin S, Jakobsson HE, Lindberg M, Rodin S, Engstrand L, Andersson AF. 2014. DegePrime, a program for degenerate primer design for broad-taxonomic-range PCR in microbial ecology studies. Appl Environ Microbiol 80:5116–5123. doi:10.1128/AEM.01403-14.
    1. Morris RM, Rappé MS, Connon SA, Vergin KL, Siebold WA, Carlson CA, Giovannoni SJ. 2002. SAR11 clade dominates ocean surface bacterioplankton communities. Nature 420:806–810. doi:10.1038/nature01240.
    1. Parada AE, Needham DM, Fuhrman JA. 14 August 2015. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol doi:10.1111/1462-2920.13023.
    1. Apprill A, McNally S, Parsons R, Weber L. 2015. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat Microb Ecol 75:129–137. doi:10.3354/ame01753.
    1. Gardes M, Bruns TD. 1993. ITS primers with enhanced specificity for basidiomycetes: application to the identification of mycorrhizae and rusts. Mol Ecol 2:113–118. doi:10.1111/j.1365-294X.1993.tb00005.x.
    1. White TJ, Bruns S, Lee S, Taylor J. 1990. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics, p. 315–322. In Innis M, Gelfand D, Sninsky K, White T (ed.), PCR protocols: a guide to methods and applications. Academic Press, New York, NY.
    1. Smith DP, Peay KG. 2014. Sequence depth, not PCR replication, improves ecological inference from next generation DNA sequencing. PLoS One 9:e90234. doi:10.1371/journal.pone.0090234.
    1. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. 2010. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336. doi:10.1038/nmeth.f.303.
    1. McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, Andersen GL, Knight R, Hugenholtz P. 2012. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J 6:610–618. doi:10.1038/ismej.2011.139.
    1. Lauber CL, Ramirez KS, Aanderud Z, Lennon J, Fierer N. 2013. Temporal variability in soil microbial communities across land-use types. ISME J 7:1641–1650. doi:10.1038/ismej.2013.50.
    1. Metcalf JL, Xu ZZ, Weiss S, Lax S, Van Treuren W, Hyde ER, Song SJ, Amir A, Larsen P, Sangwan N, Haarmann D, Humphrey GC, Ackermann G, Thompson LR, Lauber C, Bibat A, Nicholas C, Gebert MJ, Petrosino JF, Reed SC, Gilbert JA, Lynne AM, Bucheli SR, Carter DO, Knight R. 10 December 2015. Microbial community assembly and metabolic function during mammalian corpse decomposition. Science doi:10.1126/science.aad2646.
    1. Metcalf JL, Parfrey LW, Gonzalez A, Lauber CL, Knights D, Ackermann G, Humphrey GC, Gebert MJ, Van Treuren W, Berg-Lyons D, Keepers K, Guo Y, Bullard J, Fierer N, Carter DO, Knight R. 2013. A microbial clock provides an accurate estimate of the postmortem interval in a mouse model system. eLife 2:e01104. doi:10.7554/eLife.01104.
    1. Hamady M, Walker JJ, Harris JK, Gold NJ, Knight R. 2008. Error-correcting barcoded primers for pyrosequencing hundreds of samples in multiplex. Nat Methods 5:235–237. doi:10.1038/nmeth.1184.
    1. Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. 2013. 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:5112–5120. doi:10.1128/AEM.01043-13.
    1. Fadrosh DW, Ma B, Gajer P, Sengamalay N, Ott S, Brotman RM, Ravel J. 2014. An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. Microbiome 2:6. doi:10.1186/2049-2618-2-6.
    1. Walters WA, Caporaso JG, Lauber CL, Berg-Lyons D, Fierer N, Knight R. 2011. PrimerProspector: de novo design and taxonomic analysis of barcoded polymerase chain reaction primers. Bioinformatics 27:1159–1161. doi:10.1093/bioinformatics/btr087.
    1. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO. 2013. The SILVA ribosomal RNA gene database project: improved data processing and Web-based tools. Nucleic Acids Res 41:D590–D596. doi:10.1093/nar/gks1219.

Source: PubMed

3
Prenumerera