Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time-series and global field samples

Microbial community analysis via high-throughput sequencing of amplified 16S rRNA genes is an essential microbiology tool. We found the popular primer pair 515F (515F-C) and 806R greatly underestimated (e.g. SAR11) or overestimated (e.g. Gammaproteobacteria) common marine taxa. We evaluated marine samples and mock communities (containing 11 or 27 marine 16S clones), showing alternative primers 515F-Y (5’- GTGYCAGCMGCCGCGGTAA) and 926R (5’- CCGYCAATTYMTTTRAGTTT) yield more accurate estimates of mock community abundances, produce longer amplicons that can differentiate taxa unresolvable with 515F-C/806R, and amplify eukaryotic 18S rRNA. Mock communities amplified with 515F-Y/926R yielded closer observed community composition vs. expected (r2=0.95) compared to 515F-Y/806R (r2 ~0.5). Unexpectedly, biases with 515F-Y/806R against SAR11 in field samples (~4-10-fold) were stronger than in mock communities (~2-fold). Correcting a mismatch to Thaumarchaea in the 515F-C increased their apparent abundance in field samples, but not as much as using 926R rather than 806R. With plankton samples rich in eukaryotic DNA (>1?m size fraction), 18S sequences averaged ~17% of all sequences. A single mismatch can strongly bias amplification, but even perfectly matched primers can exhibit preferential amplification. We show that beyond in silico predictions, testing with mock communities and field samples is important in primer selection.

Identifier
Source https://data.blue-cloud.org/search-details?step=~012BABD9A832479CBDD56999B6F819263F2189F16EF
Metadata Access https://data.blue-cloud.org/api/collections/BABD9A832479CBDD56999B6F819263F2189F16EF
Provenance
Instrument Illumina MiSeq; ILLUMINA
Publisher Blue-Cloud Data Discovery & Access service; ELIXIR-ENA
Publication Year 2024
OpenAccess true
Contact blue-cloud-support(at)maris.nl
Representation
Discipline Marine Science
Spatial Coverage (-157.960W, -63.000S, -50.070E, 33.550N)
Temporal Coverage Begin 1989-09-25T00:00:00Z
Temporal Coverage End 2015-09-01T00:00:00Z