skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Comprehensive single‐PCR 16S and 18S rRNA community analysis validated with mock communities, and estimation of sequencing bias against 18S
Summary Universal primers for SSU rRNA genes allow profiling of natural communities by simultaneously amplifying templates from Bacteria, Archaea, and Eukaryota in a single PCR reaction. Despite the potential to show relative abundance for all rRNA genes, universal primers are rarely used, due to various concerns including amplicon length variation and its effect on bioinformatic pipelines. We thus developed 16S and 18S rRNA mock communities and a bioinformatic pipeline to validate this approach. Using these mocks, we show that universal primers (515Y/926R) outperformed eukaryote‐specific V4 primers in observed versus expected abundance correlations (slope = 0.88 vs. 0.67–0.79), and mock community members with single mismatches to the primer were strongly underestimated (threefold to eightfold). Using field samples, both primers yielded similar 18S beta‐diversity patterns (Mantel test,p < 0.001) but differences in relative proportions of many rarer taxa. To test for length biases, we mixed mock communities (16S + 18S) before PCR and found a twofold underestimation of 18S sequences due to sequencing bias. Correcting for the twofold underestimation, we estimate that, in Southern California field samples (1.2–80 μm), there were averages of 35% 18S, 28% chloroplast 16S, and 37% prokaryote 16S rRNA genes. These data demonstrate the potential for universal primers to generate comprehensive microbiome profiles.  more » « less
Award ID(s):
1737409
PAR ID:
10367512
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Environmental Microbiology
Volume:
23
Issue:
6
ISSN:
1462-2912
Page Range / eLocation ID:
p. 3240-3250
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Gilbert, Jack A. (Ed.)
    ABSTRACT Small subunit rRNA (SSU rRNA) amplicon sequencing can quantitatively and comprehensively profile natural microbiomes, representing a critically important tool for studying diverse global ecosystems. However, results will only be accurate if PCR primers perfectly match the rRNA of all organisms present. To evaluate how well marine microorganisms across all 3 domains are detected by this method, we compared commonly used primers with >300 million rRNA gene sequences retrieved from globally distributed marine metagenomes. The best-performing primers compared to 16S rRNA of bacteria and archaea were 515Y/926R and 515Y/806RB, which perfectly matched over 96% of all sequences. Considering cyanobacterial and chloroplast 16S rRNA, 515Y/926R had the highest coverage (99%), making this set ideal for quantifying marine primary producers. For eukaryotic 18S rRNA sequences, 515Y/926R also performed best (88%), followed by V4R/V4RB (18S rRNA specific; 82%)—demonstrating that the 515Y/926R combination performs best overall for all 3 domains. Using Atlantic and Pacific Ocean samples, we demonstrate high correspondence between 515Y/926R amplicon abundances (generated for this study) and metagenomic 16S rRNA (median R 2 = 0.98, n  = 272), indicating amplicons can produce equally accurate community composition data compared with shotgun metagenomics. Our analysis also revealed that expected performance of all primer sets could be improved with minor modifications, pointing toward a nearly completely universal primer set that could accurately quantify biogeochemically important taxa in ecosystems ranging from the deep sea to the surface. In addition, our reproducible bioinformatic workflow can guide microbiome researchers studying different ecosystems or human health to similarly improve existing primers and generate more accurate quantitative amplicon data. IMPORTANCE PCR amplification and sequencing of marker genes is a low-cost technique for monitoring prokaryotic and eukaryotic microbial communities across space and time but will work optimally only if environmental organisms match PCR primer sequences exactly. In this study, we evaluated how well primers match globally distributed short-read oceanic metagenomes. Our results demonstrate that primer sets vary widely in performance, and that at least for marine systems, rRNA amplicon data from some primers lack significant biases compared to metagenomes. We also show that it is theoretically possible to create a nearly universal primer set for diverse saline environments by defining a specific mixture of a few dozen oligonucleotides, and present a software pipeline that can guide rational design of primers for any environment with available meta’omic data. 
    more » « less
  2. The microbiomes of tropical corals are actively studied using 16S rRNA gene amplicons to understand microbial roles in coral health, metabolism, and disease resistance. However, primers targeting bacterial and archaeal 16S rRNA genes may also amplify organelle rRNA genes from the coral, associated microbial eukaryotes, and encrusting organisms. In this manuscript, we demonstrate that standard workflows for annotating microbial taxonomy severely under-annotate mitochondrial sequences in 1272 coral microbiomes from the Earth Microbiome Project. This issue prevents annotation of >95% of reads in some samples and persists when using either Greengenes or SILVA taxonomies. Worse, mitochondrial under-annotation varies between species and across anatomy, biasing comparisons of α- and β-diversity. By supplementing existing taxonomies with diverse mitochondrial rRNA sequences, we resolve ~97% of unique unclassified sequences as mitochondrial, without increasing misannotation in mock communities. We recommend using these extended taxonomies for coral microbiome analysis and encourage vigilance regarding similar issues in other hosts. 
    more » « less
  3. The microbiomes of tropical corals are actively studied using 16S rRNA gene amplicons to understand microbial roles in coral health, metabolism, and disease resistance. However, due to the prokaryotic origins of mitochondria, primers targeting bacterial and archaeal 16S rRNA genes may also amplify homologous 12S mitochondrial rRNA genes from the host coral, associated microbial eukaryotes, and encrusting organisms. Standard microbial bioinformatics pipelines attempt to identify and remove these sequences by comparing them to reference taxonomies. However, commonly used tools have severely under-annotated mitochondrial sequences in 1440 coral microbiomes from the Global Coral Microbiome Project, preventing annotation of over 95% of reads in some samples. This issue persists when using Greengenes or SILVA prokaryotic reference taxonomies, and in other hosts, including 16S studies of vertebrates, and of marine sponges. Worse, mitochondrial under-annotation varies between coral families and across coral compartments, biasing comparisons of  - and  -diversity. By supplementing existing reference taxonomies with over 3000 animal mitochondrial rRNA gene sequences, we resolved roughly 97% of unique unclassified sequences as mitochondrial. These additional sequences did not cause a false elevation in mitochondrial annotations in mock communities with known compositions. We recommend using these extended taxonomies for coral microbiome analysis and whenever eukaryotic contamination may be a concern. 
    more » « less
  4. Biodiversity monitoring based on DNA metabarcoding depends on primer performance. Here, we develop a new metabarcoding primer pair that targets a ~ 318 bp fragment of the 28S rRNA gene. We validate the primer pair in assessing sponges, a notoriously challenging group for coral reef metabarcoding studies, by using mock and natural complex reef communities to examine its performance in species detection, amplification efficiency, and quantitative potential. Mock community experiments revealed a high number of sponge species (n = 94) spanning a broad taxonomic scope (15 orders), limited taxon-specific primer biases (only a single species exceeded a two-fold deviation from the expected number of reads), and its suitability for quantitative metabarcoding – there was a significant relationship between read abundance and visual percent coverage of sponge taxa (R = 0.76). In the natural complex coral reef community experiments, commonly used COI metabarcoding primers detected only 30.9% of sponge species, while the new 28S primer increased detection to 79.4%. These new 28S primers detect a broader taxonomic array of species across phyla and classes within the complex cryptobiome of coral reef communities than the Leray-Geller COI primers. As biodiversity assessments using metabarcoding tools are increasingly being leveraged for environmental monitoring and guide policymaking, these new 28S rRNA primers can improve biodiversity assessments for complex ecological coral reef communities. 
    more » « less
  5. The study of the thyroid is an emerging topic, particularly in postmortem microbiome studies, due to the organ’s ability to affect the endocrine system. Also, the submandibular gland is a promising, emerging gland of study due to its position relative to the oral cavity. Previous thanatomicrobiome studies have demonstrated that bacteria belonging to the phyla Firmicutes, Proteobacteria, Bacteroides, and Pseudomonadota predominate internal organs and have been considered an important biomarker for postmortem interval. Further, Clostridium species that dominate in internal organs are linked to the hypoxic change that occurs after death, which leads to the switch of bacteria to become obligate anaerobes. Therefore, obligate anaerobes dominate the body after death due to their ability to thrive off fermentation products. 16S rRNA gene sequencing has been critical in thanatomicrobiome studies, which refers to the human microbiome (microorganisms within the body) after death. Currently, it has not been elucidated regarding the microorganisms that are associated with the decay of submandibular and thyroid glands. We hypothesized that through sequencing of the 16S rRNA gene of the submandibular and thyroid glands, the presence of Firmicutes and Proteobacteria will indicate potential biomarkers for postmortem interval. The present study revealed the postmortem microbial signatures of the submandibular and thyroid glands using the 16S rRNA gene, specifically the V3-V4 hypervariable regions, using universal primers 341F and 805R. We investigated a total of 37 cadavers obtained from ongoing criminal casework, 17 submandibular samples and 20 thyroid samples, and found that there is a correlation between microbial abundance in these postmortem glands. The predominating phyla of interest found in both glands were Firmicutes and Proteobacteria. The predominating genera were Paeniclostridium and Streptococcus in both glands, respectively. Further experimentation of the submandibular and thyroid glands will help to link oral thanatomicrobiome communities to “microbial clock” determinations, thus enhancing postmortem interval estimation. 
    more » « less