Comparing diversity levels in environmental samples: DNA sequence capture and metabarcoding approaches using 18S and COI genes

Environmental DNA studies targeting multiple taxa using metabarcoding provide remarkable insights into levels of species diversity in any habitat. The main drawbacks are the presence of primer bias and difficulty in identifying rare species. We tested a DNA sequence‐capture method in parallel with the metabarcoding approach to reveal possible advantages of one method over the other. Both approaches were performed using the same eDNA samples and the same 18S and COI regions, followed by high throughput sequencing. Metabarcoded eDNA libraries were PCR amplified with one primer pair from 18S and COI genes. DNA sequence‐capture libraries were enriched with 3,639 baits targeting the same gene regions. We tested amplicon sequence variants (ASVs) and operational taxonomic units (OTUs) in silico approaches for both markers and methods, using for this purpose the metabarcoding data set. ASVs methods uncovered more species for the COI gene, whereas the opposite occurred for the 18S gene, suggesting that clustering reads into OTUs could bias diversity richness especially using 18S with relaxed thresholds. Additionally, metabarcoding and DNA sequence‐capture recovered 80%–90% of the control sample species. DNA sequence‐capture was 8x more expensive, nonetheless it identified 1.5x more species for COI and 13x more genera for 18S than metabarcoding. Both approaches offer reliable results, sharing ca. 40% species and 72% families and retrieve more taxa when nuclear and mitochondrial markers are combined. eDNA metabarcoding is quite well established and low‐cost, whereas DNA‐sequence capture for biodiversity assessment is still in its infancy, is more time‐consuming but provides more taxonomic assignments.     

eDNA metabarcoding survey reveals fine‐scale coral reef community variation across a remote,tropical island ecosystem

Environmental DNA (eDNA) metabarcoding, a technique for retrieving multispecies DNA from environmental samples, can detect a diverse array of marine species from filtered seawater samples. There is a growing potential to integrate eDNA alongside existing monitoring methods in order to establish or improve the assessment of species diversity. Remote island reefs are increasingly vulnerable to climate‐related threats and as such there is a pressing need for cost‐effective whole‐ecosystem surveying to baseline biodiversity, study assemblage changes and ultimately develop sustainable management plans. We investigated the utility of eDNA metabarcoding as a high‐resolution, multitrophic biomonitoring tool at the Cocos (Keeling) Islands, Australia (CKI)—a remote tropical coral reef atoll situated within the eastern Indian Ocean. Metabarcoding assays targeting the mitochondrial 16S rRNA and CO1 genes, as well as the 18S rRNA nuclear gene, were applied to 252 surface seawater samples collected from 42 sites within a 140 km² area. Our assays successfully detected a wide range of bony fish and elasmobranchs (244 taxa), crustaceans (88), molluscs (37) and echinoderms (7). Assemblage composition varied significantly between sites, reflecting habitat partitioning across the island ecosystem and demonstrating the localisation of eDNA signals, despite extensive tidal and oceanic movements. In addition, we document putative new occurrence records for 46 taxa and compare the efficiency of our eDNA approach to visual survey techniques at CKI. Our study demonstrates the utility of a multimarker metabarcoding approach in capturing multitrophic biodiversity across an entire coral reef atoll and sets an important baseline for ongoing monitoring and management.     

eDNA metabarcoding as a new surveillance approach for coastal Arctic biodiversity

Because significant global changes are currently underway in the Arctic, creating a large‐scale standardized database for Arctic marine biodiversity is particularly pressing. This study evaluates the potential of aquatic environmental DNA (eDNA) metabarcoding to detect Arctic coastal biodiversity changes and characterizes the local spatio‐temporal distribution of eDNA in two locations. We extracted and amplified eDNA using two COI primer pairs from ~80 water samples that were collected across two Canadian Arctic ports, Churchill and Iqaluit, based on optimized sampling and preservation methods for remote regions surveys. Results demonstrate that aquatic eDNA surveys have the potential to document large‐scale Arctic biodiversity change by providing a rapid overview of coastal metazoan biodiversity, detecting nonindigenous species, and allowing sampling in both open water and under the ice cover by local northern‐based communities. We show that DNA sequences of ~50% of known Canadian Arctic species and potential invaders are currently present in public databases. A similar proportion of operational taxonomic units was identified at the species level with eDNA metabarcoding, for a total of 181 species identified at both sites. Despite the cold and well‐mixed coastal environment, species composition was vertically heterogeneous, in part due to river inflow in the estuarine ecosystem, and differed between the water column and tide pools. Thus, COI‐based eDNA metabarcoding may quickly improve large‐scale Arctic biomonitoring using eDNA, but we caution that aquatic eDNA sampling needs to be standardized over space and time to accurately evaluate community structure changes.     

Environmental DNA (eDNA) metabarcoding reveals strong discrimination among diverse marine habitats connected by water movement

While in recent years environmental DNA (eDNA) metabarcoding surveys have shown great promise as an alternative monitoring method, the integration into existing marine monitoring programs may be confounded by the dispersal of the eDNA signal. Currents and tidal influences could transport eDNA over great distances, inducing false‐positive species detection, leading to inaccurate biodiversity assessments and, ultimately, mismanagement of marine environments. In this study, we determined the ability of eDNA metabarcoding surveys to distinguish localized signals obtained from four marine habitats within a small spatial scale (<5 km) subject to significant tidal and along‐shore water flow. Our eDNA metabarcoding survey detected 86 genera, within 77 families and across 11 phyla using three established metabarcoding assays targeting fish (16S rRNA gene), crustacean (16S rRNA gene) and eukaryotic (cytochrome oxidase subunit 1) diversity. Ordination and cluster analyses for both taxonomic and OTU data sets show distinct eDNA signals between the sampled habitats, suggesting dispersal of eDNA among habitats was limited. Individual taxa with strong habitat preferences displayed localized eDNA signals in accordance with their respective habitat, whereas taxa known to be less habitat‐specific generated more ubiquitous signals. Our data add to evidence that eDNA metabarcoding surveys in marine environments detect a broad range of taxa that are spatially discrete. Our work also highlights that refinement of assay choice is essential to realize the full potential of eDNA metabarcoding surveys in marine biodiversity monitoring programs.     

Net overboard: Comparing marine eDNA sampling methodologies at sea to unravel marine biodiversity

Environmental DNA (eDNA) analyses are powerful for describing marine biodiversity but must be optimized for their effective use in routine monitoring. To maximize eDNA detection probabilities of sparsely distributed populations, water samples are usually concentrated from larger volumes and filtered using fine-pore membranes, often a significant cost–time bottleneck in the workflow. This study aimed to streamline eDNA sampling by investigating plankton net versus bucket sampling, direct versus sequential filtration including self-preserving filters. Biodiversity was assessed using metabarcoding of the small ribosomal subunit (18S rRNA) and mitochondrial cytochrome c oxidase I (COI) genes. Multispecies detection probabilities were estimated for each workflow using a probabilistic occupancy modelling approach. Significant workflow-related differences in biodiversity metrics were reported. Highest amplicon sequence variant (ASV) richness was attained by the bucket sampling combined with self-preserving filters, comprising a large portion of microplankton. Less diversity but more metazoan taxa were captured in the net samples combined with 5 μm pore size filters. Prefiltered 1.2 μm samples yielded few or no unique ASVs. The highest average (~32%) metazoan detection probabilities in the 5 μm pore size net samples confirmed the effectiveness of preconcentration plankton for biodiversity screening. These results contribute to streamlining eDNA sampling protocols for uptake and implementation in marine biodiversity research and surveillance.     

基于环境DNA宏条形码技术的辽东湾典型围海养殖池塘内水母多样性研究

研究使用环境DNA宏条形码技术(eDNA metabarcoding)检测辽东湾东北部河口区围海养殖池塘水母种类多样性,探索适用于水母种类物种鉴定和监测的新方法。利用环境DNA宏条形码技术,分别基于18S rDNA和COI宏条形码检测了辽东湾东北部河口区围海养殖池塘水母种类多样性,通过水样采集、过滤、eDNA提取、遗传标记扩增、测序与生物信息分析的环境DNA宏条形码标准化分析流程,从围海养殖池塘7个采样点中获得可检测的采样点数据。结果显示,基于18S rDNA宏条形码检测出8种水母种类,其中钵水母纲大型水母2种、水螅水母总纲小型水母6种;基于COI宏条形码技术共检测出19种水母种类,其中钵水母纲大型水母5种、水螅水母总纲小型水母14种;两种DNA条形码标记都显示养殖种类海蜇(Rhopilema esculentum)为优势种。研究结果表明,环境DNA宏条形码技术作为一种新兴的生物多样性监测手段可用于快速检测水母种类多样性,在水母类物种鉴定、监测及早期预警中有较大的应用潜能。     

A comparison of droplet digital polymerase chain reaction (PCR), quantitative PCR and metabarcoding for species‐specific detection in environmental DNA

Targeted species‐specific and community‐wide molecular diagnostics tools are being used with increasing frequency to detect invasive or rare species. Few studies have compared the sensitivity and specificity of these approaches. In the present study environmental DNA from 90 filtered seawater and 120 biofouling samples was analyzed with quantitative PCR (qPCR), droplet digital PCR (ddPCR) and metabarcoding targeting the cytochrome c oxidase I (COI) and 18S rRNA genes for the Mediterranean fanworm Sabella spallanzanii. The qPCR analyses detected S. spallanzanii in 53% of water and 85% of biofouling samples. Using ddPCR S. spallanzanii was detected in 61% of water of water and 95% of biofouling samples. There were strong relationships between COI copy numbers determined via qPCR and ddPCR (water R² = 0.81, p < .001, biofouling R² = 0.68, p < .001); however, qPCR copy numbers were on average 125‐fold lower than those measured using ddPCR. Using metabarcoding there was higher detection in water samples when targeting the COI (40%) compared to 18S rRNA (5.4%). The difference was less pronounced in biofouling samples (25% COI, 29% 18S rRNA). Occupancy modelling showed that although the occupancy estimate was higher for biofouling samples (ψ = 1.0), higher probabilities of detection were derived for water samples. Detection probabilities of ddPCR (1.0) and qPCR (0.93) were nearly double metabarcoding (0.57 to 0.27 marker dependent). Studies that aim to detect specific invasive or rare species in environmental samples should consider using targeted approaches until a detailed understanding of how community and matrix complexity, and primer biases affect metabarcoding data.     

DNA宏条形码技术在海洋浮游动物多样性和生态学研究中的应用

浮游动物是海洋生态系统的关键类群,其覆盖门类广泛,多样性高。传统形态鉴定技术需要检测人员具备专业的形态鉴定知识,且费时费力。宏条形码技术无需分离生物个体,而是提取拖网采集到的浮游动物混合样本的总DNA,或者水体中的环境DNA （eDNA）,依托高通量测序平台测序,能够实现对大规模样本快速、准确、经济的分析,在海洋浮游动物生态学研究中得到越来越广泛的应用。分析了DNA宏条形码技术常用的核糖体和线粒体分子标记,在浮游动物多样性和数量研究中的可靠性和不足,并给出在海洋浮游动物群落监测,食物关系分析及生物入侵早期预警等研究中的应用。未来,开发多基因片段组合条形码,发展完备的参考数据库及实现准确的量化研究是DNA宏条形码技术发展的重要方向。     

Optimization and validation of a cost‐effective protocol for biosurveillance of invasive alien species

Environmental DNA (eDNA) metabarcoding has revolutionized biodiversity monitoring and invasive pest biosurveillance programs. The introduction of insect pests considered invasive alien species (IAS) into a non‐native range poses a threat to native plant health. The early detection of IAS can allow for prompt actions by regulating authorities, thereby mitigating their impacts. In the present study, we optimized and validated a fast and cost‐effective eDNA metabarcoding protocol for biosurveillance of IAS and characterization of insect and microorganism diversity. Forty‐eight traps were placed, following the CFIA''s annual forest insect trapping survey, at four locations in southern Ontario that are high risk for forest IAS. We collected insects and eDNA samples using Lindgren funnel traps that contained a saturated salt (NaCl) solution in the collection jar. Using cytochrome c oxidase I (COI) as a molecular marker, a modified Illumina protocol effectively identified 2,535 Barcode Index Numbers (BINs). BINs were distributed among 57 Orders and 304 Families, with the vast majority being arthropods. Two IAS (Agrilus planipennis and Lymantria dispar) are regulated by the Canadian Food Inspection Agency (CFIA) as plant health pests, are known to occur in the study area, and were identified through eDNA in collected traps. Similarly, using 16S ribosomal RNA and nuclear ribosomal internal transcribed spacer (ITS), five bacterial and three fungal genera, which contain species of regulatory concern across several Canadian jurisdictions, were recovered from all sampling locations. Our study results reaffirm the effectiveness and importance of integrating eDNA metabarcoding as part of identification protocols in biosurveillance programs.     

基于环境DNA宏条码和底拖网的珠江河口鱼类多样性

文章采用环境DNA宏条码和底拖网对珠江河口鱼类多样性进行了研究, 并对两种方法进行了比较。利用环境DNA宏条码检测到了175种鱼类, 而利用底拖网采集到了47种鱼类, 结合两种方法共检测出179种鱼类, 隶属于15 目63科128属。其中两种方法共同识别了鱼类43种, 占总检测物种的24.02%, 基于底拖网的调查未能收集到基于环境DNA宏条码检测到的大多数物种。根据Shannon指数和Simpson指数显示, DNA宏条码所检测珠江河口鱼类群落α多样性显著高于底拖网方法(P<0.05)。两种方法的PCoA结果均显示珠江河口鱼类群落存在空间结构, 基于环境DNA宏条码的分析显示空间重叠更多。两种方法基于冗余分析均显示溶解氧和盐度是影响鱼类群落结构的主要环境因子。研究表明, 环境DNA 宏条形码是一种环保且可靠的评估方法, 将其搭载到现有调查可以更好地了解河口鱼类多样性。     

基于环境DNA-宏条形码技术的底栖动物监测及水质评价研究进展

底栖动物是淡水生态系统中物种多样性最高的类群,也是应用最广泛的水质监测指示生物之一。传统的底栖动物监测以形态学为基础,耗时费力,无法满足流域尺度大规模监测的需求。环境DNA-宏条形码技术是一种新兴的生物监测方法,其与传统方法相比优势在于采样方法简单、低成本、高灵敏度,不受生物样本和环境状况的影响,不依赖分类专家和鉴定资料,能够快速准确地对多个类群进行大规模、高通量的物种鉴定。然而,在实际应用中该方法的效果受诸多因素的影响,不同的方法、流程往往会产生差异较大的结果。鉴于此,着重分析总结了应用环境DNA-宏条形码技术监测底栖动物的关键影响因素,包括样品采集与处理流程、分子标记选择、引物设计、PCR偏好性、参考数据库的完整性及相应的优化。并基于此探讨了提高环境DNA-宏条形码技术在底栖动物监测效率和准确率的途径,以期为底栖动物环境DNA-宏条形码监测方案的制定提供可靠的参考。最后对该技术在底栖动物监测和水质评价中的最新发展方向进行了展望。     

A DNA mini‐barcode for marine macrophytes

Studies focusing on marine macrophyte metabarcoding from environmental samples are scarce, due to the lack of a universal barcode for these taxa, and to their poor representation in DNA databases. Here, we searched for a short barcode able to identify marine macrophytes from tissue samples; then, we created a DNA reference library which was used to identify macrophytes in eDNA from coastal sediments. Barcoding of seagrasses, mangroves and marine macroalgae (Chlorophyta, Rhodophyta and Phaeophyceae) was tested using 18 primer pairs from six barcoding genes: the plant barcodes rbcL, matK and trnL, plus the genes ITS2, COI and 18S. The 18S gene showed the highest universality among marine macrophytes, amplifying 95%–100% of samples; amplification performance of the other barcodes was limited. Taxonomy was assigned using a phylogeny‐based approach to create an 18S DNA reference library. Macrophyte tissue sequences were accurately identified within their phyla (88%), order (76%), genus (71%) and species (23%). Nevertheless, out of 86 macrophytes tested, only 48% and 15% had a reference sequence at genus and at species level, respectively. Identification at these levels can be improved by more inclusive reference libraries. Using the 18S mini‐barcode and the reference library, we recovered eDNA from 21 marine macrophytes in sediments, demonstrating the barcode's ability to trace primary producers that contribute to blue carbon. We expect this barcode to also be useful for other ecological questions, such as tracing macro primary producers in marine food webs.     

Assessing vertebrate biodiversity in a kelp forest ecosystem using environmental DNA

Preserving biodiversity is a global challenge requiring data on species’ distribution and abundance over large geographic and temporal scales. However, traditional methods to survey mobile species’ distribution and abundance in marine environments are often inefficient, environmentally destructive, or resource‐intensive. Metabarcoding of environmental DNA (eDNA) offers a new means to assess biodiversity and on much larger scales, but adoption of this approach for surveying whole animal communities in large, dynamic aquatic systems has been slowed by significant unknowns surrounding error rates of detection and relevant spatial resolution of eDNA surveys. Here, we report the results of a 2.5 km eDNA transect surveying the vertebrate fauna present along a gradation of diverse marine habitats associated with a kelp forest ecosystem. Using PCR primers that target the mitochondrial 12S rRNA gene of marine fishes and mammals, we generated eDNA sequence data and compared it to simultaneous visual dive surveys. We find spatial concordance between individual species’ eDNA and visual survey trends, and that eDNA is able to distinguish vertebrate community assemblages from habitats separated by as little as ~60 m. eDNA reliably detected vertebrates with low false‐negative error rates (1/12 taxa) when compared to the surveys, and revealed cryptic species known to occupy the habitats but overlooked by visual methods. This study also presents an explicit accounting of false negatives and positives in metabarcoding data, which illustrate the influence of gene marker selection, replication, contamination, biases impacting eDNA count data and ecology of target species on eDNA detection rates in an open ecosystem.     

DNA metabarcoding and the cytochrome c oxidase subunit I marker: not a perfect match

DNA metabarcoding enables efficient characterization of species composition in environmental DNA or bulk biodiversity samples, and this approach is making significant and unique contributions in the field of ecology. In metabarcoding of animals, the cytochrome c oxidase subunit I (COI) gene is frequently used as the marker of choice because no other genetic region can be found in taxonomically verified databases with sequences covering so many taxa. However, the accuracy of metabarcoding datasets is dependent on recovery of the targeted taxa using conserved amplification primers. We argue that COI does not contain suitably conserved regions for most amplicon-based metabarcoding applications. Marker selection deserves increased scrutiny and available marker choices should be broadened in order to maximize potential in this exciting field of research.     

Effect of marker choice and thermal cycling protocol on zooplankton DNA metabarcoding studies

DNA metabarcoding is a promising approach for rapidly surveying biodiversity and is likely to become an important tool for measuring ecosystem responses to environmental change. Metabarcoding markers need sufficient taxonomic coverage to detect groups of interest, sufficient sequence divergence to resolve species, and will ideally indicate relative abundance of taxa present. We characterized zooplankton assemblages with three different metabarcoding markers (nuclear 18S rDNA, mitochondrial COI, and mitochondrial 16S rDNA) to compare their performance in terms of taxonomic coverage, taxonomic resolution, and correspondence between morphology‐ and DNA‐based identification. COI amplicons sequenced on separate runs showed that operational taxonomic units representing >0.1% of reads per sample were highly reproducible, although slightly more taxa were detected using a lower annealing temperature. Mitochondrial COI and nuclear 18S showed similar taxonomic coverage across zooplankton phyla. However, mitochondrial COI resolved up to threefold more taxa to species compared to 18S. All markers revealed similar patterns of beta‐diversity, although different taxa were identified as the greatest contributors to these patterns for 18S. For calanoid copepod families, all markers displayed a positive relationship between biomass and sequence reads, although the relationship was typically strongest for 18S. The use of COI for metabarcoding has been questioned due to lack of conserved primer‐binding sites. However, our results show the taxonomic coverage and resolution provided by degenerate COI primers, combined with a comparatively well‐developed reference sequence database, make them valuable metabarcoding markers for biodiversity assessment.     

Estimating species richness using environmental DNA

The foundation for any ecological study and for the effective management of biodiversity in natural systems requires knowing what species are present in an ecosystem. We assessed fish communities in a stream using two methods, depletion‐based electrofishing and environmental DNA metabarcoding (eDNA) from water samples, to test the hypothesis that eDNA provides an alternative means of determining species richness and species identities for a natural ecosystem. In a northern Indiana stream, electrofishing yielded a direct estimate of 12 species and a mean estimated richness (Chao II estimator) of 16.6 species with a 95% confidence interval from 12.8 to 42.2. eDNA sampling detected an additional four species, congruent with the mean Chao II estimate from electrofishing. This increased detection rate for fish species between methods suggests that eDNA sampling can enhance estimation of fish fauna in flowing waters while having minimal sampling impacts on fish and their habitat. Modern genetic approaches therefore have the potential to transform our ability to build a more complete list of species for ecological investigations and inform management of aquatic ecosystems.     

How many replicates to accurately estimate fish biodiversity using environmental DNA on coral reefs?

Quantifying fish species diversity in rich tropical marine environments remains challenging. Environmental DNA (eDNA) metabarcoding is a promising tool to face this challenge through the filtering, amplification, and sequencing of DNA traces from water samples. However, because eDNA concentration is low in marine environments, the reliability of eDNA to detect species diversity can be limited. Using an eDNA metabarcoding approach to identify fish Molecular Taxonomic Units (MOTUs) with a single 12S marker, we aimed to assess how the number of sampling replicates and filtered water volume affect biodiversity estimates. We used a paired sampling design of 30 L per replicate on 68 reef transects from 8 sites in 3 tropical regions. We quantified local and regional sampling variability by comparing MOTU richness, compositional turnover, and compositional nestedness. We found strong turnover of MOTUs between replicated pairs of samples undertaken in the same location, time, and conditions. Paired samples contained non‐overlapping assemblages rather than subsets of one another. As a result, non‐saturated localized diversity accumulation curves suggest that even 6 replicates (180 L) in the same location can underestimate local diversity (for an area <1 km). However, sampling regional diversity using ~25 replicates in variable locations (often covering 10 s of km) often saturated biodiversity accumulation curves. Our results demonstrate variability of diversity estimates possibly arising from heterogeneous distribution of eDNA in seawater, highly skewed frequencies of eDNA traces per MOTU, in addition to variability in eDNA processing. This high compositional variability has consequences for using eDNA to monitor temporal and spatial biodiversity changes in local assemblages. Avoiding false‐negative detections in future biomonitoring efforts requires increasing replicates or sampled water volume to better inform management of marine biodiversity using eDNA.     

Environmental DNA metabarcoding of wild flowers reveals diverse communities of terrestrial arthropods

Terrestrial arthropods comprise the most species‐rich communities on Earth, and grassland flowers provide resources for hundreds of thousands of arthropod species. Diverse grassland ecosystems worldwide are threatened by various types of environmental change, which has led to decline in arthropod diversity. At the same time, monitoring grassland arthropod diversity is time‐consuming and strictly dependent on declining taxonomic expertise. Environmental DNA (eDNA) metabarcoding of complex samples has demonstrated that information on species compositions can be efficiently and non‐invasively obtained. Here, we test the potential of wild flowers as a novel source of arthropod eDNA. We performed eDNA metabarcoding of flowers from several different plant species using two sets of generic primers, targeting the mitochondrial genes 16S rRNA and COI. Our results show that terrestrial arthropod species leave traces of DNA on the flowers that they interact with. We obtained eDNA from at least 135 arthropod species in 67 families and 14 orders, together representing diverse ecological groups including pollinators, parasitoids, gall inducers, predators, and phytophagous species. Arthropod communities clustered together according to plant species. Our data also indicate that this experiment was not exhaustive, and that an even higher arthropod richness could be obtained using this eDNA approach. Overall, our results demonstrate that it is possible to obtain information on diverse communities of insects and other terrestrial arthropods from eDNA metabarcoding of wild flowers. This novel source of eDNA represents a vast potential for addressing fundamental research questions in ecology, obtaining data on cryptic and unknown species of plant‐associated arthropods, as well as applied research on pest management or conservation of endangered species such as wild pollinators.     

Environmental DNA metabarcoding: Transforming how we survey animal and plant communities

The genomic revolution has fundamentally changed how we survey biodiversity on earth. High‐throughput sequencing (“HTS”) platforms now enable the rapid sequencing of DNA from diverse kinds of environmental samples (termed “environmental DNA” or “eDNA”). Coupling HTS with our ability to associate sequences from eDNA with a taxonomic name is called “eDNA metabarcoding” and offers a powerful molecular tool capable of noninvasively surveying species richness from many ecosystems. Here, we review the use of eDNA metabarcoding for surveying animal and plant richness, and the challenges in using eDNA approaches to estimate relative abundance. We highlight eDNA applications in freshwater, marine and terrestrial environments, and in this broad context, we distill what is known about the ability of different eDNA sample types to approximate richness in space and across time. We provide guiding questions for study design and discuss the eDNA metabarcoding workflow with a focus on primers and library preparation methods. We additionally discuss important criteria for consideration of bioinformatic filtering of data sets, with recommendations for increasing transparency. Finally, looking to the future, we discuss emerging applications of eDNA metabarcoding in ecology, conservation, invasion biology, biomonitoring, and how eDNA metabarcoding can empower citizen science and biodiversity education.