skip to main content

Title: DNA metabarcoding marker choice skews perception of marine eukaryotic biodiversity

DNA metabarcoding is an increasingly popular technique to investigate biodiversity; however, many methodological unknowns remain, especially concerning the biases resulting from marker choice. Regions of the cytochromecoxidase subunit I (COI) and 18S rDNA (18S) genes are commonly employed “universal” markers for eukaryotes, but the extent of taxonomic biases introduced by these markers and how such biases may impact metabarcoding performance is not well quantified. Here, focusing on macroeukaryotes, we use standardized sampling from autonomous reef monitoring structures (ARMS) deployed in the world's most biodiverse marine ecosystem, the Coral Triangle, to compare the performance of COI and 18S markers. We then compared metabarcoding data to image‐based annotations of ARMS plates. Although both markers provided similar estimates of taxonomic richness and total sequence reads, marker choice skewed estimates of eukaryotic diversity. The COI marker recovered relative abundances of the dominant sessile phyla consistent with image annotations. Both COI and the image annotations provided higher relative abundance estimates of Bryozoa and Porifera and lower estimates of Chordata as compared to 18S, but 18S recovered 25% more phyla than COI. Thus, while COI more reliably reflects the occurrence of dominant sessile phyla, 18S provides a more holistic representation of overall taxonomic diversity. Ideal marker choice is, therefore, contingent on study system and research question, especially in relation to desired taxonomic resolution, and a multimarker approach provides the greatest application across a broad range of research objectives. As metabarcoding becomes an essential tool to monitor biodiversity in our changing world, it is critical to evaluate biases associated with marker choice.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Environmental DNA
Page Range / eLocation ID:
p. 1229-1246
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Biodiversity is changing at an accelerating rate at both local and regional scales. Beta diversity, which quantifies species turnover between these two scales, is emerging as a key driver of ecosystem function that can inform spatial conservation. Yet measuring biodiversity remains a major challenge, especially in aquatic ecosystems. Decoding environmental DNA (eDNA) left behind by organisms offers the possibility of detecting species sans direct observation, a Rosetta Stone for biodiversity. While eDNA has proven useful to illuminate diversity in aquatic ecosystems, its utility for measuring beta diversity over spatial scales small enough to be relevant to conservation purposes is poorly known. Here we tested how eDNA performs relative to underwater visual census (UVC) to evaluate beta diversity of marine communities. We paired UVC with 12S eDNA metabarcoding and used a spatially structured hierarchical sampling design to assess key spatial metrics of fish communities on temperate rocky reefs in southern California. eDNA provided a more-detailed picture of the main sources of spatial variation in both taxonomic richness and community turnover, which primarily arose due to strong species filtering within and among rocky reefs. As expected, eDNA detected more taxa at the regional scale (69 vs. 38) which accumulated quickly with space and plateaued at only ~ 11 samples. Conversely, the discovery rate of new taxa was slower with no sign of saturation for UVC. Based on historical records in the region (2000–2018) we found that 6.9 times more UVC samples would be required to detect 50 taxa compared to eDNA. Our results show that eDNA metabarcoding can outperform diver counts to capture the spatial patterns in biodiversity at fine scales with less field effort and more power than traditional methods, supporting the notion that eDNA is a critical scientific tool for detecting biodiversity changes in aquatic ecosystems.

    more » « less
  2. Fields, David (Ed.)
    Abstract Metabarcoding of environmental DNA has provided striking insights into marine microbial diversity. With this approach, we assessed the diversity of metazoan assemblages and their temporal variations at the Long-Term Ecological Research site MareChiara (LTER-MC) in the Gulf of Naples (Mediterranean Sea). The Illumina sequencing of the V4-18S rRNA fragment from 48 surface samples collected from 2011 to 2013 produced a total of 5 011 047 marine metazoan reads. The normalized dataset was generally dominated by copepods (60.3%), followed by annelids (34.7%) mostly represented by the invasive benthic polychaete Hydroides elegans. Non-copepod holoplankton was mainly represented by siphonophores, rotifers, and appendicularians, with occasional mass occurrences of jellyfish. The rest of meroplankton (mainly molluscs, annelids, and anthozoans) showed a high diversity, covering all 11 zooplankton phyla, from Porifera to Craniata. A high number of copepod genera were identified, with seasonal recurrence matching patterns observed in 30 years of studies in the Gulf of Naples. Despite limitations related to the molecular marker resolution and reference dataset, the study provided valuable insights into diversity and seasonal patterns of the whole metazoan assemblage, expanding the knowledge on rare or hardly identifiable taxa and confirming DNA metabarcoding as a powerful approach to be integrated at LTER sites. 
    more » « less
  3. Abstract

    High-throughput sequencing of a targeted genetic marker is being widely used to analyze biodiversity across taxa and environments. Amid a multitude of exciting findings, scientists have also identified and addressed technical and biological limitations. Improved study designs and alternative sampling, lab and bioinformatic procedures have progressively enhanced data quality, but some problems persist. This article provides a framework to recognize and bypass the main types of errors that can affect metabarcoding data: false negatives, false positives, artifactual variants, disproportions and incomplete or incorrect taxonomic identifications. It is crucial to discern potential error impacts on different ecological parameters (e.g. taxon distribution, community structure, alpha and beta-diversity), as error management implies compromises and is thus directed by the research question. Synthesis of multiple plankton metabarcoding evaluations (mock sample sequencing or microscope comparisons) shows that high-quality data for qualitative and some semiquantitative goals can be achieved by implementing three checkpoints: first, rigorous protocol optimization; second, error minimization; and third, downstream analysis that considers potentially remaining biases. Conclusions inform us about the reliability of metabarcoding for plankton studies and, because plankton provides unique chances to compare genotypes and phenotypes, the robustness of this method in general.

    more » « less
  4. Sergio Stefanni (Ed.)

    Zooplankton diversity in the deep “midnight zone” (>1000 m), where sunlight does not reach, remains largely unknown. Uncovering such diversity has been challenging because of the major difficulties in sampling deep pelagic fauna and identifying many (unknown) species that belong to these complex swimmer assemblages. In this study, we evaluated zooplankton diversity using two taxonomic marker genes: mitochondrial cytochrome oxidase subunit 1 (COI) and nuclear 18S ribosomal RNA (18S). We collected samples from plankton net tows, ranging from the surface to a depth of 5000 m above the Atacama Trench in the Southeast Pacific. Our study aimed to assess the zooplankton diversity among layers from the upper 1000 m to the ultra-deep abyssopelagic zone to test the hypothesis of decreasing diversity with depth resulting from limited carbon sources. The results showed unique, highly vertically structured communities within the five depth strata sampled, with maximal species richness observed in the upper bathypelagic layer (1000–2000 m). The high species richness of zooplankton (>750 OTUS) at these depths was higher than that found in the upper 1000 m. The vertical diversity trend exhibited a pattern similar to the well-known vertical pattern described for the benthic system. However, a large part of this diversity was either unknown (>50%) or could not be assigned to any known species in current genetic diversity databases. DNA analysis showed that the Calanoid copepods, mostly represented bySubeucalanus monachus, the Euphausiacea,Euphausia mucronata, and the halocypridade,Paraconchoecia dasyophthalma, dominated the community. Water column temperature, dissolved oxygen, particulate carbon, and nitrogen appeared to be related to the observed vertical diversity pattern. Our findings revealed rich and little-known zooplankton diversity in the deep sea, emphasizing the importance of further exploration of this ecosystem to conserve and protect its unique biota.

    more » « less
  5. ProkaryoticNostoc, one of the world's most conspicuous and widespread algal genera (similar to eukaryotic algae, plants, and animals) is known to support a microbiome that influences host ecological roles. Past taxonomic characterizations of surface microbiota (epimicrobiota) of free‐livingNostocsampled from freshwater systems employed 16S rRNA genes, typically amplicons. We compared taxa identified from 16S, 18S, 23S, and 28S rRNA gene sequences filtered from shotgun metagenomic sequence and used microscopy to illuminate epimicrobiota diversity forNostocsampled from a wetland in the northern Chilean Altiplano. Phylogenetic analysis and rRNA gene sequence abundance estimates indicated that the host was related toNostoc punctiformePCC 73102. Epimicrobiota were inferred to include 18 epicyanobacterial genera or uncultured taxa, six epieukaryotic algal genera, and 66 anoxygenic bacterial genera, all having average genomic coverage ≥90X. The epicyanobacteriaGeitlerinemia,Oscillatoria,Phormidium, and an uncultured taxon were detected only by 16S rRNA gene;GloeobacterandPseudanabaenawere detected using 16S and 23S; andPhormididesmis,Neosynechococcus,Symphothece,Aphanizomenon,Nodularia,Spirulina,Nodosilinea,Synechococcus,Cyanobium, andAnabaena(the latter corroborated by microscopy), plus two uncultured cyanobacterial taxa (JSC12, O77) were detected only by 23S rRNA gene sequences. Three chlamydomonad and two heterotrophic stramenopiles genera were inferred from 18S; the streptophyte green algaChaetosphaeridium globosumwas detected by microscopy and 28S rRNA genes, but not 18S rRNA genes. Overall, >60% of epimicrobial taxa were detected by markers other than 16S rRNA genes. Some algal taxa observed microscopically were not detected from sequence data. Results indicate that multiple taxonomic markers derived from metagenomic sequence data and microscopy increase epimicrobiota detection.

    more » « less