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.     

Quantification of mesocosm fish and amphibian species diversity via environmental DNA metabarcoding

Freshwater fauna are particularly sensitive to environmental change and disturbance. Management agencies frequently use fish and amphibian biodiversity as indicators of ecosystem health and a way to prioritize and assess management strategies. Traditional aquatic bioassessment that relies on capture of organisms via nets, traps and electrofishing gear typically has low detection probabilities for rare species and can injure individuals of protected species. Our objective was to determine whether environmental DNA (eDNA) sampling and metabarcoding analysis can be used to accurately measure species diversity in aquatic assemblages with differing structures. We manipulated the density and relative abundance of eight fish and one amphibian species in replicated 206‐L mesocosms. Environmental DNA was filtered from water samples, and six mitochondrial gene fragments were Illumina‐sequenced to measure species diversity in each mesocosm. Metabarcoding detected all nine species in all treatment replicates. Additionally, we found a modest, but positive relationship between species abundance and sequencing read abundance. Our results illustrate the potential for eDNA sampling and metabarcoding approaches to improve quantification of aquatic species diversity in natural environments and point the way towards using eDNA metabarcoding as an index of macrofaunal species abundance.     

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.     

Targeted and passive environmental DNA approaches outperform established methods for detection of quagga mussels,Dreissena rostriformis bugensis in flowing water

1. The early detection of invasive non‐native species (INNS) is important for informing management actions. Established monitoring methods require the collection or observation of specimens, which is unlikely at the beginning of an invasion when densities are likely to be low. Environmental DNA (eDNA) analysis is a highly promising technique for the detection of INNS—particularly during the early stages of an invasion.
2. Here, we compared the use of traditional kick‐net sampling with two eDNA approaches (targeted detection using both conventional and quantitative PCR and passive detection via metabarcoding with conserved primers) for detection of quagga mussel, Dreissena rostriformis bugensis, a high priority INNS, along a density gradient on the River Wraysbury, UK.
3. All three molecular tools outperformed traditional sampling in terms of detection. Conventional PCR and qPCR both had 100% detection rate in all samples and outperformed metabarcoding when the target species was at low densities. Additionally, quagga mussel DNA copy number (qPCR) and relative read count (metabarcoding) were significantly influenced by both mussel density and distance from source population, with distance being the most significant predictor.
4. Synthesis and application. All three molecular approaches were more sensitive than traditional kick‐net sampling for the detection of the quagga mussel in flowing water, and both qPCR and metabarcoding enabled estimates of relative abundance. Targeted approaches were more sensitive than metabarcoding, but metabarcoding has the advantage of providing information on the wider community and consequently the impacts of INNS.

     

Environmental DNA metabarcoding primers for freshwater fish detection and quantification: In silico and in tanks

Environmental DNA (eDNA) techniques refer to utilizing the organisms’ DNA extracted from environment samples to genetically identify target species without capturing actual organisms. eDNA metabarcoding via high‐throughput sequencing can simultaneously detect multiple fish species from a single water sample, which is a powerful tool for the qualitative detection and quantitative estimates of multiple fish species. However, sequence counts obtained from eDNA metabarcoding may be influenced by many factors, of which primer bias is one of the foremost causes of methodological error. The performance of 18 primer pairs for COI, cytb, 12S rRNA, and 16S rRNA mitochondrial genes, which are all frequently used in fish eDNA metabarcoding, were evaluated in the current study. The ribosomal gene markers performed better than the protein‐coding gene markers during in silico screening, resulting in higher taxonomic coverage and appropriate barcode lengths. Four primer pairs—AcMDB07, MiFish‐U, Ve16S1, and Ve16S3—designed for various regions of the 12S and 16S rRNA genes were screened for tank metabarcoding in a case study targeting six freshwater fish species. The four primer pairs were able to accurately detect all six species in different tanks, while only MiFish‐U, Ve16S1, and Ve16S3 revealed a significant positive relationship between species biomass and read count for the pooled tank data. The positive relationship could not be found in all species within the tanks. Additionally, primer efficiency differed depending on the species while primer preferential species varied in different fish assemblages. This case study supports the potential for eDNA metabarcoding to assess species diversity in natural ecosystems and provides an alternative strategy to evaluate the performance of candidate primers before application of eDNA metabarcoding in natural ecosystems.     

Assessing insect biodiversity with automatic light traps in Brazil: Pearls and pitfalls of metabarcoding samples in preservative ethanol

Automated species identification based on data produced with metabarcoding offers an alternative for assessing biodiversity of bulk insect samples obtained with traps. We used a standard two‐step PCR approach to amplify a 313 bp fragment of the barcoding region of the mitochondrial COI gene. The PCR products were sequenced on an Illumina MiSeq platform, and the OTUs production and taxonomic identifications were performed with a customized pipeline and database. The DNA used in the PCR procedures was extracted directly from the preservative ethanol of bulk insect samples obtained with automatic light traps in 12 sampling areas located in different biomes of Brazil, during wet and dry seasons. Agricultural field and forest edge habitats were collected for all sampling areas. A total of 119 insect OTUs and nine additional OTUs assigned to other arthropod taxa were obtained at a ≥97% sequence similarity level. The alpha and beta diversity analyses comparing biomes, habitats, and seasons were mostly inconclusive, except for a significant difference in beta diversity between biomes. In this study, we were able to metabarcode and HTS adult insects from their preservative medium. Notwithstanding, our results underrepresent the true magnitude of insect diversity expected from samples obtained with automatic light traps in Brazil. Although biological and technical factors might have impacted our results, measures to optimize and standardize eDNA HTS should be in place to improve taxonomic coverage of samples of unknown diversity and stored in suboptimal conditions, which is the case of most eDNA samples.     

Environmental DNA metabarcoding studies are critically affected by substrate selection

Effective biomonitoring is critical for driving management outcomes that ensure long‐term sustainability of the marine environment. In recent years, environmental DNA (eDNA), coupled with metabarcoding methodologies, has emerged as a promising tool for generating biotic surveys of marine ecosystems, including those under anthropogenic pressure. However, more empirical data are needed on how to best implement eDNA field sampling approaches to maximize their utility for each specific application. The effect of the substrate chosen for eDNA sampling on the diversity of marine taxa detected by DNA metabarcoding has not yet been systematically analysed, despite aquatic systems being those most commonly targeted for eDNA studies. We investigated the effect of four commonly used eDNA substrates to explore taxonomic diversity: (a) surface water, (b) marine sediment, (c) settlement plates and (d) planktonic tows. With a focus on coastal ports, 332 eDNA samples from Australia (Indian and Southern oceans) and Kazakhstan (Caspian Sea) were collected and analysed by multi‐assay DNA metabarcoding. Across study locations, between 30% and 52% of eukaryotic families detected were unique to a particular substrate and <6% of families were found in all four substrates. Taxonomic composition varied significantly depending on the substrate sampled implying that the suitability (and bias) of an eDNA substrate will depend on the focal taxa. These findings demonstrate that single substrate eDNA metabarcoding likely underestimates the total eukaryotic diversity. Future eDNA experimental design should consider incorporating multiple substrates or select substrate(s) best suited to the specific detection of target taxa.     

Metabarcoding of fungal communities associated with bark beetles

Many species of fungi are closely allied with bark beetles, including many tree pathogens, but their species richness and patterns of distribution remain largely unknown. We established a protocol for metabarcoding of fungal communities directly from total genomic DNA extracted from individual beetles, showing that the ITS3/4 primer pair selectively amplifies the fungal ITS. Using three specimens of bark beetle from different species, we assess the fungal diversity associated with these specimens and the repeatability of these estimates in PCRs conducted with different primer tags. The combined replicates produced 727 fungal Operational Taxonomic Units (OTUs) for the specimen of Hylastes ater, 435 OTUs for Tomicus piniperda, and 294 OTUs for Trypodendron lineatum, while individual PCR reactions produced on average only 229, 54, and 31 OTUs for the three specimens, respectively. Yet, communities from PCR replicates were very similar in pairwise comparisons, in particular when considering species abundance, but differed greatly among the three beetle specimens. Different primer tags or the inclusion of amplicons in separate libraries did not impact the species composition. The ITS2 sequences were identified with the Lowest Common Ancestor approach and correspond to diverse lineages of fungi, including Ophiostomaceae and Leotiomycetes widely found to be tree pathogens. We conclude that Illumina MiSeq metabarcoding reliably captures fungal diversity associated with bark beetles, although numerous PCR replicates are recommended for an exhaustive sample. Direct PCR from beetle DNA extractions provides a rapid method for future surveys of fungal species diversity and their associations with bark beetles and environmental variables.     

Taking eDNA underground: Factors affecting eDNA detection of subterranean fauna in groundwater

Stygofauna are aquatic fauna that have evolved to live underground. The impacts of anthropogenic climate change, extraction and pollution on groundwater pose major threats to groundwater health, prompting the need for efficient and reliable means to detect and monitor stygofaunal communities. Conventional survey techniques for these species rely on morphological identification and can be biased, labour-intensive and often indeterminate to lower taxonomic levels. By contrast, environmental DNA (eDNA)-based methods have the potential to dramatically improve on existing stygofaunal survey methods in a large range of habitats and for all life stages, reducing the need for the destructive manual collection of often critically endangered species or for specialized taxonomic expertise. We compared eDNA and haul-net samples collected in 2020 and 2021 from 19 groundwater bores and a cave on Barrow Island, northwest Western Australia, and assessed how sampling factors influenced the quality of eDNA detection of stygofauna. The two detection methods were complementary; eDNA metabarcoding was able to detect soft-bodied taxa and fish often missed by nets, but only detected seven of the nine stygofaunal crustacean orders identified from haul-net specimens. Our results also indicated that eDNA metabarcoding could detect 54%–100% of stygofauna from shallow-water samples and 82%–90% from sediment samples. However, there was significant variation in stygofaunal diversity between sample years and sampling types. The findings of this study demonstrate that haul-net sampling has a tendency to underestimate stygofaunal diversity and that eDNA metabarcoding of groundwater can substantially improve the efficiency of stygofaunal surveys.     

Use of invertebrate traits for the biomonitoring of European large rivers: the effects of sampling effort on genus richness and functional diversity

1. Studies on biodiversity and ecosystem function require considering metrics for accurately describing the functional diversity of communities. The number of taxa (richness) is commonly used to characterise biological diversity. The disadvantage of richness as a measure of biological diversity is that all taxa are taken into account on an equal basis regardless of their abundance, their biological characteristics or their function in the ecosystem. 2. To circumvent this problem, we applied a recently described measure of biological diversity that incorporates dissimilarities among taxa. Dissimilarities were defined from biological traits (e.g. life history, morphology, physiology and behaviour) of stream invertebrate taxa and the resulting biological diversity index was considered as a surrogate for functional diversity. 3. As sampling effort is known to affect the number of taxa collected within a reach, we investigated how change in functional diversity is affected by sampling effort. We used stream invertebrate community data from three large European rivers to model accumulation curves and to assess the number of samples required to estimate (i.e. closeness to the maximal value) functional diversity and genera richness. We further evaluated the precision of estimates (i.e. similarity of temporal or spatial replicates) of the total functional diversity. 4. As expected, richness estimates were strongly dependent on sampling effort, and 10 replicate samples were found to underestimate actual richness. Moreover, richness estimates showed much variation with season and location. In contrast, functional diversity had greater accuracy with less sampling effort and the precision of the estimates was higher than richness both across sampling occasions and sampling reaches. These results are further arguments towards conducting research on the design of a biomonitoring tool based on biological traits.     

Toward accurate species‐level metabarcoding of arthropod communities from the tropical forest canopy

相似文献   

Assessment of fish communities using environmental DNA: Effect of spatial sampling design in lentic systems of different sizes

Freshwater fish biodiversity is quickly decreasing and requires effective monitoring and conservation. Environmental DNA (eDNA)‐based methods have been shown to be highly sensitive and cost‐efficient for aquatic biodiversity surveys, but few studies have systematically investigated how spatial sampling design affects eDNA‐detected fish communities across lentic systems of different sizes. We compared the spatial patterns of fish diversity determined using eDNA in three lakes of small (SL; 3 ha), medium (ML; 122 ha) and large (LL; 4,343 ha) size using a spatially explicit grid sampling method. A total of 100 water samples (including nine, 17 and 18 shoreline samples and six, 14 and 36 interior samples from SL, ML and LL, respectively) were collected, and fish communities were analysed using eDNA metabarcoding of the mitochondrial 12S region. Together, 30, 35 and 41 fish taxa were detected in samples from SL, ML, and LL, respectively. We observed that eDNA from shoreline samples effectively captured the majority of the fish diversity of entire waterbodies, and pooled samples recovered fewer species than individually processed samples. Significant spatial autocorrelations between fish communities within 250 m and 2 km of each other were detected in ML and LL, respectively. Additionally, the relative sequence abundances of many fish species exhibited spatial distribution patterns that correlated with their typical habitat occupation. Overall, our results support the validity of a shoreline sampling strategy for eDNA‐based fish community surveys in lentic systems but also suggest that a spatially comprehensive sampling design can reveal finer distribution patterns of individual species.     

The importance of including imperfect detection models in eDNA experimental design

Environmental DNA (eDNA) is DNA that has been isolated from field samples, and it is increasingly used to infer the presence or absence of particular species in an ecosystem. However, the combination of sampling procedures and subsequent molecular amplification of eDNA can lead to spurious results. As such, it is imperative that eDNA studies include a statistical framework for interpreting eDNA presence/absence data. We reviewed published literature for studies that utilized eDNA where the species density was known and compared the probability of detecting the focal species to the sampling and analysis protocols. Although biomass of the target species and the volume per sample did not impact detectability, the number of field replicates and number of samples from each replicate were positively related to detection. Additionally, increased number of PCR replicates and increased primer specificity significantly increased detectability. Accordingly, we advocate for increased use of occupancy modelling as a method to incorporate effects of sampling effort and PCR sensitivity in eDNA study design. Based on simulation results and the hierarchical nature of occupancy models, we suggest that field replicates, as opposed to molecular replicates, result in better detection probabilities of target species.     

Needle in a haystack? A comparison of eDNA metabarcoding and targeted qPCR for detection of the great crested newt (Triturus cristatus)

Environmental DNA (eDNA) analysis is a rapid, cost‐effective, non‐invasive biodiversity monitoring tool which utilises DNA left behind in the environment by organisms for species detection. The method is used as a species‐specific survey tool for rare or invasive species across a broad range of ecosystems. Recently, eDNA and “metabarcoding” have been combined to describe whole communities rather than focusing on single target species. However, whether metabarcoding is as sensitive as targeted approaches for rare species detection remains to be evaluated. The great crested newt Triturus cristatus is a flagship pond species of international conservation concern and the first UK species to be routinely monitored using eDNA. We evaluate whether eDNA metabarcoding has comparable sensitivity to targeted real‐time quantitative PCR (qPCR) for T. cristatus detection. Extracted eDNA samples (N = 532) were screened for T. cristatus by qPCR and analysed for all vertebrate species using high‐throughput sequencing technology. With qPCR and a detection threshold of 1 of 12 positive qPCR replicates, newts were detected in 50% of ponds. Detection decreased to 32% when the threshold was increased to 4 of 12 positive qPCR replicates. With metabarcoding, newts were detected in 34% of ponds without a detection threshold, and in 28% of ponds when a threshold (0.028%) was applied. Therefore, qPCR provided greater detection than metabarcoding but metabarcoding detection with no threshold was equivalent to qPCR with a stringent detection threshold. The proportion of T. cristatus sequences in each sample was positively associated with the number of positive qPCR replicates (qPCR score) suggesting eDNA metabarcoding may be indicative of eDNA concentration. eDNA metabarcoding holds enormous potential for holistic biodiversity assessment and routine freshwater monitoring. We advocate this community approach to freshwater monitoring to guide management and conservation, whereby entire communities can be initially surveyed to best inform use of funding and time for species‐specific surveys.     

Metabarcoding Is Powerful yet Still Blind: A Comparative Analysis of Morphological and Molecular Surveys of Seagrass Communities

In the context of the sixth wave of extinction, reliable surveys of biodiversity are increasingly needed to infer the cause and consequences of species and community declines, identify early warning indicators of tipping points, and provide reliable impact assessments before engaging in activities with potential environmental hazards. DNA metabarcoding has emerged as having potential to provide speedy assessment of community structure from environmental samples. Here we tested the reliability of metabarcoding by comparing morphological and molecular inventories of invertebrate communities associated with seagrasses through estimates of alpha and beta diversity, as well as the identification of the most abundant taxa. Sediment samples were collected from six Zostera marina seagrass meadows across Brittany, France. Metabarcoding surveys were performed using both mitochondrial (Cytochrome Oxidase I) and nuclear (small subunit 18S ribosomal RNA) markers, and compared to morphological inventories compiled by a long-term benthic monitoring network. A sampling strategy was defined to enhance performance and accuracy of results by preventing the dominance of larger animals, boosting statistical support through replicates, and using two genes to compensate for taxonomic biases. Molecular barcodes proved powerful by revealing a remarkable level of diversity that vastly exceeded the morphological survey, while both surveys identified congruent differentiation of the meadows. However, despite the addition of individual barcodes of common species into taxonomic reference databases, the retrieval of only 36% of these species suggest that the remaining were either not present in the molecular samples or not detected by the molecular screening. This finding exemplifies the necessity of comprehensive and well-curated taxonomic reference libraries and multi-gene surveys. Overall, results offer methodological guidelines and support for metabarcoding as a powerful and repeatable method of characterizing communities, while also presenting suggestions for improvement, including implementation of pilot studies prior to performing full “blind” metabarcoding assessments to optimize sampling and amplification protocols.     

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.     

Replication levels,false presences and the estimation of the presence/absence from eDNA metabarcoding data

Environmental DNA (eDNA) metabarcoding is increasingly used to study the present and past biodiversity. eDNA analyses often rely on amplification of very small quantities or degraded DNA. To avoid missing detection of taxa that are actually present (false negatives), multiple extractions and amplifications of the same samples are often performed. However, the level of replication needed for reliable estimates of the presence/absence patterns remains an unaddressed topic. Furthermore, degraded DNA and PCR/sequencing errors might produce false positives. We used simulations and empirical data to evaluate the level of replication required for accurate detection of targeted taxa in different contexts and to assess the performance of methods used to reduce the risk of false detections. Furthermore, we evaluated whether statistical approaches developed to estimate occupancy in the presence of observational errors can successfully estimate true prevalence, detection probability and false‐positive rates. Replications reduced the rate of false negatives; the optimal level of replication was strongly dependent on the detection probability of taxa. Occupancy models successfully estimated true prevalence, detection probability and false‐positive rates, but their performance increased with the number of replicates. At least eight PCR replicates should be performed if detection probability is not high, such as in ancient DNA studies. Multiple DNA extractions from the same sample yielded consistent results; in some cases, collecting multiple samples from the same locality allowed detecting more species. The optimal level of replication for accurate species detection strongly varies among studies and could be explicitly estimated to improve the reliability of results.     

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

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

Identifying error and accurately interpreting environmental DNA metabarcoding results: A case study to detect vertebrates at arid zone waterholes

Environmental DNA (eDNA) metabarcoding surveys enable rapid, noninvasive identification of taxa from trace samples with wide‐ranging applications from characterizing local biodiversity to identifying food‐web interactions. However, the technique is prone to error from two major sources: (a) contamination through foreign DNA entering the workflow, and (b) misidentification of DNA within the workflow. Both types of error have the potential to obscure true taxon presence or to increase taxonomic richness by incorrectly identifying taxa as present at sample sites, but multiple error sources can remain unaccounted for in metabarcoding studies. Here, we use data from an eDNA metabarcoding study designed to detect vertebrate species at waterholes in Australia's arid zone to illustrate where and how in the workflow errors can arise, and how to mitigate those errors. We detected the DNA of 36 taxa spanning 34 families, 19 orders and five vertebrate classes in water samples from waterholes, demonstrating the potential for eDNA metabarcoding surveys to provide rapid, noninvasive detection in remote locations, and to widely sample taxonomic diversity from aquatic through to terrestrial taxa. However, we initially identified 152 taxa in the samples, meaning there were many false positive detections. We identified the sources of these errors, allowing us to design a stepwise process to detect and remove error, and provide a template to minimize similar errors that are likely to arise in other metabarcoding studies. Our findings suggest eDNA metabarcoding surveys need to be carefully conducted and screened for errors to ensure their accuracy.