Next-generation freshwater bioassessment: eDNA metabarcoding with a conserved metazoan primer reveals species-rich and reservoir-specific communities

Development of environmental DNA metabarcoding primers for marine mollusks and comparison with published primers

Monitoring mollusk biodiversity is a great challenge due to their large diversity and broad distribution. Environmental DNA (eDNA) technology is increasingly applied for biodiversity monitoring, but relevant studies on marine mollusks are still limited. Although previous studies have developed several pairs of primers for mollusk eDNA analyses, most of them targeted only a small group of mollusks. In this study, seven primers were designed for the mollusk community and validated and compared with eight pairs of published primers to select the best candidates. After in silico test, MollCOI154 and MollCOI255 primers showed non-specific amplification, and same results were also obtained in published primers (COI204, Sepi, and veneroida). Moll12S100, Moll12S195 and Moll16S primers failed to amplify across all genomic DNA from selected mollusk. Except Moll16S, all developed and two published (unionoida and veneroida) primers were successfully amplified on four eDNA samples from Yangtze River estuary. After annotation of the amplified sequences, MollCOI253 showed higher annotation of the amplification results than the other primers. In conclusion, MollCOI253 had better performance in terms of amplification success and specificity, and can provide technical support for eDNA-based research, which will be beneficial for molluscan biodiversity investigation and conservation.

Sequence capture identifies fastidious chytrid fungi directly from host tissue

The chytrid fungus Batrachochytrium dendrobatidis (Bd) was discovered in 1998 as the cause of chytridiomycosis, an emerging infectious disease causing mass declines in amphibian populations worldwide. The rapid population declines of the 1970s-1990s were likely caused by the spread of a highly virulent lineage belonging to the Bd-GPL clade that was introduced to naïve susceptible populations. Multiple genetically distinct and regional lineages of Bd have since been isolated and sequenced, greatly expanding the known biological diversity within this fungal pathogen. To date, most Bd research has been restricted to the limited number of samples that could be isolated using culturing techniques, potentially causing a selection bias for strains that can grow on media and missing other unculturable or fastidious strains that are also present on amphibians. We thus attempted to characterize potentially non-culturable genetic lineages of Bd from distinct amphibian taxa using sequence capture technology on DNA extracted from host tissue and swabs. We focused our efforts on host taxa from two different regions that likely harbored distinct Bd clades: (1) wild-caught leopard frogs (Rana) from North America, and (2) a Japanese Giant Salamander (Andrias japonicus) at the Smithsonian Institution's National Zoological Park that exhibited signs of disease and tested positive for Bd using qPCR, but multiple attempts failed to isolate and culture the strain for physiological and genetic characterization. We successfully enriched for and sequenced thousands of fungal genes from both host clades, and Bd load was positively associated with number of recovered Bd sequences. Phylogenetic reconstruction placed all the Rana-derived strains in the Bd-GPL clade. In contrast, the A. japonicus strain fell within the Bd-Asia3 clade, expanding the range of this clade and generating additional genomic data to confirm its placement. The retrieved ITS locus matched public barcoding data from wild A. japonicus and Bd infections found on other amphibians in India and China, suggesting that this uncultured clade is widespread across Asia. Our study underscores the importance of recognizing and characterizing the hidden diversity of fastidious strains in order to reconstruct the spatiotemporal and evolutionary history of Bd. The success of the sequence capture approach highlights the utility of directly sequencing pathogen DNA from host tissue to characterize cryptic diversity that is missed by culture-reliant approaches.

Freshwater eDNA reveals dramatic biological shifts linked to deforestation of New Zealand

Deforestation is considered a major threat to biodiversity across many parts of the globe, but the biological impacts of this dramatic ecosystem disturbance often remain incompletely understood. In New Zealand - the world's last major landmass to be colonised by humans - widespread deforestation over recent centuries has left a highly fragmented suite of relict forest stands, ideal for assessing anthropogenic biological change. We hypothesise that this widespread environmental disturbance has underpinned repeated and predictable ecological shifts across distinct rivers and regions. Here we use freshwater environmental DNA (eDNA) data (113 samples across 38 locations; 89 insect taxa) to test for concordant biological shifts linked to this deforestation. eDNA analyses highlight consistent compositional and functional differentiation between forested versus deforested assemblages, including turnover of 'cryptic' congeneric taxa that are morphologically similar yet ecologically and genetically distinct. These dramatic biological shifts are evident even over fine spatial scales within streams, emphasising the widespread emergence of a novel 'deforested' assemblage. Our results illustrate that environmental change can drive predictable biological shifts across broad geographic regions, and highlight the power of eDNA for assessing anthropogenic ecosystem change over large geographic scales.

Network analysis with either Illumina or MinION reveals that detecting vertebrate species requires metabarcoding of iDNA from a diverse fly community

DNA obtained from invertebrates (iDNA) can be metabarcoded in order to survey vertebrate communities. However, little attention has been paid to the interaction between the invertebrate and vertebrate species. Here, we tested for specialization by sampling the dung and carrion fly community of a swamp forest remnant along a disturbance gradient (10 sites: 80-310 m from a road). Approximately, 60% of the baited 407 flies yielded 294 vertebrate identifications based on two COI fragments and 16S. A bipartite network analysis found no statistically significant specialization in the interactions between fly and vertebrate species, but uncommon fly species can carry the signal for vertebrate species that are otherwise difficult to detect with iDNA. A spatial analysis revealed that most of the 20 vertebrate species reported in this study could be detected within 150 m of the road (18 spp.) and that the fly community sourced for iDNA was unexpectedly rich (24 species, 3 families). They carried DNA for rare and common species inhabiting different layers of the forest (ground-dwelling: wild boar, Sunda pangolin, skinks, rats; arboreal: long-tailed macaque, Raffles' banded langur; flying: pin-striped tit-babbler, olive-winged bulbul). All our results were obtained with a new, greatly simplified iDNA protocol that eliminates DNA extraction by obtaining template directly through dissolving fly faeces and regurgitates with water. Lastly, we show that MinION- and Illumina-based metabarcoding yield similar results. We conclude by urging more studies that use different baits and involve experiments that are capable of revealing the dispersal capabilities of the flies carrying the iDNA.

Creating, curating and evaluating a mitogenomic reference database to improve regional species identification using environmental DNA

Species detection using eDNA is revolutionizing global capacity to monitor biodiversity. However, the lack of regional, vouchered, genomic sequence information-especially sequence information that includes intraspecific variation-creates a bottleneck for management agencies wanting to harness the complete power of eDNA to monitor taxa and implement eDNA analyses. eDNA studies depend upon regional databases of mitogenomic sequence information to evaluate the effectiveness of such data to detect and identify taxa. We created the Oregon Biodiversity Genome Project to create a database of complete, nearly error-free mitogenomic sequences for all of Oregon's fishes. We have successfully assembled the complete mitogenomes of 313 specimens of freshwater, anadromous and estuarine fishes representing 24 families, 55 genera and 129 species and lineages. Comparative analyses of these sequences illustrate that many regions of the mitogenome are taxonomically informative, that the short (~150 bp) mitochondrial 'barcode' regions typically used for eDNA assays do not consistently diagnose for species and that complete single or multiple genes of the mitogenome are preferable for identifying Oregon's fishes. This project provides a blueprint for other researchers to follow as they build regional databases, illustrates the taxonomic value and limits of complete mitogenomic sequences and offers clues as to how current eDNA assays and environmental genomics methods of the future can best leverage this information.

Environmental DNA surveillance of biocontamination in a drinking water treatment plant

A clean and adequate supply of drinking water is essential to life and good health. However, despite the risk of biologically derived contamination of drinking water, monitoring of invertebrate outbreaks has relied primarily on naked-eye inspections that are prone to errors. In this study, we applied environmental DNA (eDNA) metabarcoding as a biomonitoring tool at seven different stages of drinking water treatment, from prefiltration to release from household faucets. While the composition of invertebrate eDNA communities reflected the communities of the source water in earlier stages of the treatment, several predominant invertebrate taxa (e.g., rotifer) were shown to be introduced during purification, but most were eliminated in later treatment stages. In addition, the limit of detection/quantification of PCR assay and read capacity of high-throughput sequencing was assessed with further microcosm experiments to estimate the applicability eDNA metabarcoding to the biocontamination surveillance in drinking water treatment plants (DWTPs). Here we propose a novel eDNA-based approach for sensitive and efficient surveillance of invertebrate outbreaks in DWTPs.

Aquatic environmental DNA: A review of the macro-organismal biomonitoring revolution

Environmental DNA (eDNA) is the fastest growing biomonitoring tool fuelled by two key features: time efficiency and sensitivity. Technological advancements allow rapid biodiversity detection at both species and community levels with increasing accuracy. Concurrently, there has been a global demand to standardise eDNA methods, but this is only possible with an in-depth overview of the technological advancements and a discussion of the pros and cons of available methods. We therefore conducted a systematic literature review of 407 peer-reviewed papers on aquatic eDNA published between 2012 and 2021. We observed a gradual increase in the annual number of publications from four (2012) to 28 (2018), followed by a rapid growth to 124 publications in 2021. This was mirrored by a tremendous diversification of methods in all aspects of the eDNA workflow. For example, in 2012 only freezing was applied to preserve filter samples, whereas we recorded 12 different preservation methods in the 2021 literature. Despite an ongoing standardisation debate in the eDNA community, the field is seemingly moving fast in the opposite direction and we discuss the reasons and implications. Moreover, by compiling the largest PCR-primer database to date, we provide information on 522 and 141 published species-specific and metabarcoding primers targeting a wide range of aquatic organisms. This works as a user-friendly 'distillation' of primer information that was hitherto scattered across hundreds of papers, but the list also reflects which taxa are commonly studied with eDNA technology in aquatic environments such as fish and amphibians, and reveals that groups such as corals, plankton and algae are under-studied. Efforts to improve sampling and extraction methods, primer specificity and reference databases are crucial to capture these ecologically important taxa in future eDNA biomonitoring surveys. In a rapidly diversifying field, this review synthetises aquatic eDNA procedures and can guide eDNA users towards best practice.

Quantitative environmental DNA metabarcoding shows high potential as a novel approach to quantitatively assess fish community

The simultaneous conservation of species richness and evenness is important to effectively reduce biodiversity loss and keep ecosystem health. Environmental DNA (eDNA) metabarcoding has been used as a powerful tool for identifying community composition, but it does not necessarily provide quantitative information due to several methodological limitations. Thus, the quantification of eDNA through metabarcoding is an important frontier of eDNA-based biomonitoring. Particularly, the qMiSeq approach has recently been developed as a quantitative metabarcoding method and has attracted much attention due to its usefulness. The aim here was to evaluate the performance of the qMiSeq approach as a quantitative monitoring tool for fish communities by comparing the quantified eDNA concentrations with the results of fish capture surveys. The eDNA water sampling and the capture surveys using the electrical shocker were conducted at a total of 21 sites in four rivers in Japan. As a result, we found significant positive relationships between the eDNA concentrations of each species quantified by qMiSeq and both the abundance and biomass of each captured taxon at each site. Furthermore, for seven out of eleven taxa, a significant positive relationship was observed between quantified DNA concentrations by sample and the abundance and/or biomass. In total, our results demonstrated that eDNA metabarcoding with the qMiSeq approach is a suitable and useful tool for quantitative monitoring of fish communities. Due to the simplicity of the eDNA analysis, the eDNA metabarcoding with qMiSeq approach would promote further growth of quantitative monitoring of biodiversity.

The Effects of Water Level Fluctuation on Zooplankton Communities in Shahu Lake Based on DNA Metabarcoding and Morphological Methods

Background: The water level of Poyang Lake (China) fluctuates seasonally. Shahu Lake, a smaller body of water connected to Poyang Lake during the wet season, is separated in the dry season. Due to a special fishing method termed ‘lake enclosed in autumn’, the water level is lowered and reaches its lowest point in January, which is <0.5 m deep in the middle of the lake. Our research investigated the effect of water level changes on the zooplankton community composition in Shahu Lake. Methods: We used both DNA metabarcoding method (MBC) (18S rRNA gene V4 region) and morphological method (MOI) to track the zooplankton community structure over four seasons in Shahu Lake (China). Results: Totals of 90 and 98 species of zooplankton were detected by MOI and MBC, respectively, with rotifers being the main zooplankton component. The α-diversity index of both methods increased from spring to summer and decreased from summer to autumn, reaching the lowest value in winter. NMDS and a cluster analysis showed that all zooplankton communities detected by MOI and MBC were significantly separated by season. The zooplankton community in winter was separated from that of the other three seasons, but the summer and autumn communities were more similar. Conclusions: Changes in the water level had significant effects on the zooplankton community composition. We found that MBC was more able to detect the differences in the zooplankton composition than MOI. MBC also had more advantages in copepod recognition. In our study, 37 species of copepods were detected by MBC, but only 11 species were detected by MOI. We concluded that MBC should be used to research the seasonal variations of zooplankton.

Towards Improved Detection and Identification of Rust Fungal Pathogens in Environmental Samples Using a Metabarcoding Approach

The dispersion of fungal inocula such as the airborne spores of rust fungi (Pucciniales) can be monitored through metabarcoding of the internal transcribed spacer 2 (ITS2) of the rRNA gene in environmental DNAs. This method is largely dependent on a high-quality reference database (refDB) and primers with proper taxonomic coverage and specificity. For this study, a curated ITS2 reference database (named CR-ITS2-refDB) comprising representatives of the major cereal rust fungi and phylogenetically related species was compiled. Interspecific and intraspecific variation analyses suggested that the ITS2 region had reasonable discriminating power for the majority of the Puccinia species or species complexes in the database. In silico evaluation of nine forward and seven reverse ITS2 primers, including three newly designed, revealed marked variation in DNA amplification efficiency for the rusts. We validated the theoretical assessment of rust-enhanced (Rust2inv/ITS4var_H) and universal fungal (ITS9F/ITS4) ITS2 primer pairs by profiling the airborne rust fungal communities from environmental samples via a metabarcoding approach. Species- or subspecies-level identification of the rusts was improved by use of CR-ITS2-refDB and the Automated Oligonucleotide Design Pipeline (AODP), which identified all mutations distinguishing highly conserved DNA markers between close relatives. A generic bioinformatics pipeline was developed, including all steps used in this study from in silico evaluation of primers to accurate identification of short metabarcodes at the level of interest for defining phytopathogens. The results highlight the importance of primer selection, refDBs that are resolved to reflect phylogenetic relationships, and the use of AODP for improving the reliability of metabarcoding in phytopathogen biosurveillance.

Molecular Diet Analysis of Adélie Penguins (Pygoscelis adeliae) in the Ross Sea Using Fecal DNA

Simple Summary

The diet of Adélie penguins, Pygoscelis adeliae, in the Ross Sea was studied applying quantitative polymerase chain reaction (qPCR) protocols to their feces. Two krill species (Euphausia superba and Euphausia crystallorophias) and notothenioid fish (mainly Pleuragramma antarctica and Pagothenia borchgrevinki) were among the most abundant components of the diet. The composition of the two krill species and notothenioid fish was found to be strongly related to the geographic characteristics of the Ross Sea.

Abstract

The diet of Adélie penguins, Pygoscelis adeliae, is a useful indicator in understanding the ecological conditions of their habitats. The diets of Adélie penguins were studied using metabarcoding and quantitative PCR (qPCR) analyses of fecal DNA from seven habitats along the Ross Sea region. Using metabarcoding analysis with dual universal primers (18Sv9 and miniFish), the overall diet composition and detailed information about piscine prey were clearly elucidated. It was found that two krill species (Euphausia superba and Euphausia crystallorophias) and notothenioid fish were the most abundant in the diets of Adélie penguins. Among the notothenioid prey, Pleuragramma antarctica (56.50%) and Pagothenia borchgrevinki (18.21%) were the two most abundant species. qPCR analysis showed a significant geographic difference in the composition of main prey. Penguins inhabiting outbound parts of the Ross Sea (Capes Adare (CA) and Duke of York Island (DY)) mainly preyed on E. superba, without any significant changes in prey composition. By contrast, those inhabiting the inbound parts of the Ross Sea (Edmonson Point (EP) and Inexpressible Island (II)) preyed on E. crystallorophias and notothenioid fish rather than E. superba. Compared with the outbound habitats, prey compositions for penguins inhabiting the inbound regions were significantly different year to year, which was presumably due to the food availability based on the annual environmental and meteorological conditions of the coastal region along with the inbound parts of the Ross Sea.

Estimating the genetic diversity of Pacific salmon and trout using multigene eDNA metabarcoding

Genetic diversity underpins species conservation and management goals, and ultimately determines a species' ability to adapt. Using freshwater environmental DNA (eDNA) samples, we examined mitochondrial genetic diversity using multigene metabarcode sequence data from four Oncorhynchus species across 16 sites in Oregon and northern California. Our multigene metabarcode panel included targets commonly used in population genetic NADH dehydrogenase 2 (ND2), phylogenetic cytochrome c oxidase subunit 1 (COI) and eDNA (12S ribosomal DNA) screening. The ND2 locus showed the greatest within-species haplotype diversity for all species, followed by COI and then 12S rDNA for all species except Oncorhynchus kisutch. Sequences recovered for O. clarkii clarkii were either identical to, or one mutation different from, previously characterized haplotypes (95.3% and 4.5% of reads, respectively). The greatest diversity in O. c. clarkii was among coastal watersheds, and subsets of this diversity were shared with fish in inland watersheds. However, coastal streams and the Umpqua River watershed appear to harbour unique haplotypes. Sequences from O. mykiss revealed a disjunction between the Willamette watershed and southern watersheds suggesting divergent histories. We also identified similarities between populations in the northern Deschutes and southern Klamath watersheds, consistent with previously hypothesized connections between the two via inland basins. Oncorhynchus kisutch was only identified in coastal streams and the Klamath River watershed, with most diversity concentrated in the coastal Coquille watershed. Oncorhynchus tshawytscha was only observed at one site, but contained multiple haplotypes at each locus. The characterization of genetic diversity at multiple loci expands the knowledge gained from eDNA sampling and provides crucial information for conservation actions and genetic management.

ONTbarcoder and MinION barcodes aid biodiversity discovery and identification by everyone, for everyone

BACKGROUND

DNA barcodes are a useful tool for discovering, understanding, and monitoring biodiversity which are critical tasks at a time of rapid biodiversity loss. However, widespread adoption of barcodes requires cost-effective and simple barcoding methods. We here present a workflow that satisfies these conditions. It was developed via "innovation through subtraction" and thus requires minimal lab equipment, can be learned within days, reduces the barcode sequencing cost to < 10 cents, and allows fast turnaround from specimen to sequence by using the portable MinION sequencer.

RESULTS

We describe how tagged amplicons can be obtained and sequenced with the real-time MinION sequencer in many settings (field stations, biodiversity labs, citizen science labs, schools). We also provide amplicon coverage recommendations that are based on several runs of the latest generation of MinION flow cells ("R10.3") which suggest that each run can generate barcodes for > 10,000 specimens. Next, we present a novel software, ONTbarcoder, which overcomes the bioinformatics challenges posed by MinION reads. The software is compatible with Windows 10, Macintosh, and Linux, has a graphical user interface (GUI), and can generate thousands of barcodes on a standard laptop within hours based on only two input files (FASTQ, demultiplexing file). We document that MinION barcodes are virtually identical to Sanger and Illumina barcodes for the same specimens (> 99.99%) and provide evidence that MinION flow cells and reads have improved rapidly since 2018.

CONCLUSIONS

We propose that barcoding with MinION is the way forward for government agencies, universities, museums, and schools because it combines low consumable and capital cost with scalability. Small projects can use the flow cell dongle ("Flongle") while large projects can rely on MinION flow cells that can be stopped and re-used after collecting sufficient data for a given project.

Seeing through sedimented waters: environmental DNA reduces the phantom diversity of sharks and rays in turbid marine habitats

BACKGROUND

Sharks and rays are some of the most threatened marine taxa due to the high levels of bycatch and significant demand for meat and fin-related products in many Asian communities. At least 25% of shark and ray species are considered to be threatened with extinction. In particular, the density of reef sharks in the Pacific has declined to 3-10% of pre-human levels. Elasmobranchs are thought to be sparse in highly urbanised and turbid environments. Low visibility coupled with the highly elusive behaviour of sharks and rays pose a challenge to diversity estimation and biomonitoring efforts as sightings are limited to chance encounters or from carcasses ensnared in nets. Here we utilised an eDNA metabarcoding approach to enhance the precision of elasmobranch diversity estimates in urbanised marine environments.

RESULTS

We applied eDNA metabarcoding on seawater samples to detect elasmobranch species in the hyper-urbanised waters off Singapore. Two genes-vertebrate 12S and elasmobranch COI-were targeted and amplicons subjected to Illumina high-throughput sequencing. With a total of 84 water samples collected from nine localities, we found 47 shark and ray molecular operational taxonomic units, of which 16 had species-level identities. When data were compared against historical collections and contemporary sightings, eDNA detected 14 locally known species as well as two potential new records.

CONCLUSIONS

Local elasmobranch richness uncovered by eDNA is greater than the seven species sighted over the last two decades, thereby reducing phantom diversity. Our findings demonstrate that eDNA metabarcoding is effective in detecting shark and ray species despite the challenges posed by the physical environment, granting a more consistent approach to monitor these highly elusive and threatened species.

Pond Acoustic Sampling Scheme: A draft protocol for rapid acoustic data collection in small waterbodies

Freshwater conservation is vital to the maintenance of global biodiversity. Ponds are a critical, yet often under-recognized, part of this, contributing to overall ecosystem functioning and diversity. They provide habitats for a range of aquatic, terrestrial, and amphibious life, often including rare and declining species.Effective, rapid, and accessible survey methods are needed to enable evidence-based conservation action, but freshwater taxa are often viewed as "difficult"-and few specialist surveyors are available. Datasets on ponds are therefore limited in their spatiotemporal coverage.With the advent of new recording technologies, acoustic survey methods are becoming increasingly available to researchers, citizen scientists, and conservation practitioners. They can be an effective and noninvasive approach for gathering data on target species, assemblages, and environmental variables. However, freshwater applications are lagging behind those in terrestrial and marine spheres, and as an emergent method, research studies have employed a multitude of different sampling protocols.We propose the Pond Acoustic Sampling Scheme (PASS), a simple protocol to allow a standardized minimal sample to be collected rapidly from small waterbodies, alongside environmental and methodological metadata. This sampling scheme can be incorporated into a variety of survey designs and is intended to allow access to a wide range of participants, without requiring complicated or prohibitively expensive equipment.Adoption of this sampling protocol would enable consistent sound recordings to be gathered by researchers and conservation organizations, and allow the development of landscape-scale surveys, data sharing, and collaboration within an expanding freshwater ecoacoustic community-rather than individual approaches that produce incompatible datasets. The compilation of standardized data would improve the prospects for effective research into the soundscapes of small waterbodies and aid freshwater conservation efforts.

Validating eDNA measurements of the richness and abundance of anurans at a large scale

In some situations, environmental DNA (eDNA) metabarcoding is a more accurate tool for measuring species richness of a taxon than conventional methods. Whether this tool can reliably estimate the abundance of a taxon remains unresolved. We examined the reliability of metabarcoding for measuring anuran diversity compared to a commonly used traditional line transect method (TLTM) through the replicate sampling of three visits across 71 waterbodies (ponds or reservoirs) in Liuheng, China. We also investigated the relative contributions of species-specific characteristics and the physiochemical properties of a waterbody on the relative read count across species and waterbodies. We found that eDNA metabarcoding had a higher detection probability for each of seven anuran species found in the sampling region than TLTM. Furthermore, the relative read count estimated by metabarcoding was positively correlated with the density or relative density of individuals identified with the TLTM across waterbodies for every species. Species-specific characteristics of anurans, such as density, relative density, body mass, biomass and relative biomass, accounted for substantial variations in the read count across species and waterbodies, while physiochemical factors, including pH, temperature, water volume, vegetation and elevation, had little effect on the read count. Our results based on robust sampling suggest that metabarcoding enables more reliable and efficient measurements of anuran occurrence at a large scale during a short-term survey (within 15 days) than that obtained by the TLTM, and offers an alternative tool for quantifications of anuran abundance. Density or biomass is better and more reliable indicator of anuran abundance associated with read count than relative density or relative biomass.

DNA metabarcoding of zooplankton communities: species diversity and seasonal variation revealed by 18S rRNA and COI

Background

Zooplankton is an important component of aquatic organisms and has important biological and economical significance in freshwater ecosystems. However, traditional methods that rely on morphology to classify zooplankton require expert taxonomic skills. Moreover, traditional classification methods are time-consuming and labor-intensive, which is not practical for the design of conservation measures and ecological management tools based on zooplankton diversity assessment.

Methods

We used DNA metabarcoding technology with two different markers: the nuclear small subunit ribosomal RNA (18S rRNA) and mitochondrial cytochrome c oxidase (COI), to analyze 72 zooplankton samples collected in 4 seasons and 9 locations from the Sanmenxia Reservoir. We investigated seasonal changes in the zooplankton community and their relationship with water environmental factors.

Results

A total of 190 species of zooplankton were found, belonging to 12 phyla, 24 classes, 61 orders, 111 families, and 174 genera. Protozoa, especially ciliates, were the most diverse taxa. Richness and relative abundance of zooplankton showed significant seasonal changes. Both alpha and beta diversity showed seasonal trends: the diversity in summer and autumn was higher than that in winter and spring. The zooplankton diversity was most similar in winter and spring. By correlating metabarcoding data and water environmental factors, we proved that water temperature, chemical oxygen demand, total nitrogen and ammoniacal nitrogen were the main environmental factors driving the seasonal changes in zooplankton in the Sanmenxia Reservoir. Water temperature, followed by total nitrogen, were the most influential factors. This study highlights the advantages and some limitations of zooplankton molecular biodiversity assessment using two molecular markers.

Invasive species trait-based risk assessment for non-native freshwater fishes in a tropical city basin in Southeast Asia

Biological invasions have created detrimental impacts in freshwater ecosystems. As non-native freshwater species include economically beneficial, but also harmful, species, trait-based risk assessments can be used to identify and prevent the import of potentially invasive species. Freshwater fishes are one of the most evaluated freshwater taxa to date. However, such assessments have mostly been done in sub-temperate to temperate regions, with a general lack of such research in the tropics. In view of this knowledge gap, this study aims to determine if a different set of traits are associated with successful establishment of non-native fishes within the tropics. In tropical Southeast Asia, Singapore represents a suitable model site to perform an invasive species trait-based risk assessment for the tropical region given its susceptibility to the introduction and establishment of non-native freshwater fishes and lack of stringent fish import regulation. A quantitative trait-based risk assessment was performed using random forest to determine the relative importance of species attributes associated with the successful establishment of introduced freshwater fishes in Singapore. Species having a match in climate, prior invasion success, lower absolute fecundity, higher trophic level, and involvement in the aquarium trade were found to have higher establishment likelihood (as opposed to native distributional range and maximum size being among the commonly identified predictors in subtropical/temperate trait-based risk assessments). To minimize invasive risk, incoming freshwater fishes could be screened in future for such traits, allowing lists of prohibited or regulated species to be updated. The findings could also potentially benefit the development of invasive species action plans and inform management decisions in the Southeast Asian region. Considering a geographical bias in terms of having relatively less documentation of biological invasions in the tropics, particularly Asia, this study highlights the need to perform more of such risk assessments in other parts of the tropics.

Integration of DNA-Based Approaches in Aquatic Ecological Assessment Using Benthic Macroinvertebrates

Benthic macroinvertebrates are among the most used biological quality elements for assessing the condition of all types of aquatic ecosystems worldwide (i.e., fresh water, transitional, and marine). Current morphology-based assessments have several limitations that may be circumvented by using DNA-based approaches. Here, we present a comprehensive review of 90 publications on the use of DNA metabarcoding of benthic macroinvertebrates in aquatic ecosystems bioassessments. Metabarcoding of bulk macrozoobenthos has been preferentially used in fresh waters, whereas in marine waters, environmental DNA (eDNA) from sediment and bulk communities from deployed artificial structures has been favored. DNA extraction has been done predominantly through commercial kits, and cytochrome c oxidase subunit I (COI) has been, by far, the most used marker, occasionally combined with others, namely, the 18S rRNA gene. Current limitations include the lack of standardized protocols and broad-coverage primers, the incompleteness of reference libraries, and the inability to reliably extrapolate abundance data. In addition, morphology versus DNA benchmarking of ecological status and biotic indexes are required to allow general worldwide implementation and higher end-user confidence. The increased sensitivity, high throughput, and faster execution of DNA metabarcoding can provide much higher spatial and temporal data resolution on aquatic ecological status, thereby being more responsive to immediate management needs.

Relationship of taxonomic error to frequency of observation

Biological nomenclature is the entry point to a wealth of information related to or associated with living entities. When applied accurately and consistently, communication between and among researchers and investigators is enhanced, leading to advancements in understanding and progress in research programs. Based on freshwater benthic macroinvertebrate taxonomic identifications, inter-laboratory comparisons of >900 samples taken from rivers, streams, and lakes across the U.S., including the Great Lakes, provided data on taxon-specific error rates. Using the error rates in combination with frequency of observation (FREQ; as a surrogate for rarity), six uncertainty/frequency classes (UFC) are proposed for approximately 1,000 taxa. The UFC, error rates, FREQ each are potentially useful for additional analyses related to interpreting biological assessment results and/or stressor response relationships, as weighting factors for various aspects of ecological condition or biodiversity analyses and helping set direction for taxonomic research and refining identification tools.

A powerful long metabarcoding method for the determination of complex diets from faecal analysis of the European pond turtle (Emys orbicularis, L. 1758)

High‐throughput sequencing has become an accurate method for the identification of species present in soil, water, faeces, gut or stomach contents. However, information at the species level is limited due to the choice of short barcodes and based on the idea that DNA is too degraded to allow longer sequences to be amplified. We have therefore developed a long DNA metabarcoding method based on the sequencing of short reads followed by de novo assembly, which can precisely identify the taxonomic groups of organisms associated with complex diets, such as omnivorous individuals. The procedure includes 11 different primer pairs targeting the COI gene, the large subunit of the ribulose‐1,5‐bisphosphate carboxylase gene, the maturase K gene, the 28S rRNA and the trnL‐trnF chloroplastic region. We validated this approach using 32 faeces samples from an omnivorous reptile, the European pond turtle (Emys orbicularis, L. 1758). This metabarcoding approach was assessed using controlled experiments including mock communities and faecal samples from captive feeding trials. The method allowed us to accurately identify prey DNA present in the diet of the European pond turtles to the species level in most of the cases (82.4%), based on the amplicon lengths of multiple markers (168–1,379 bp, average 546 bp), and produced by de novo assembly. The proposed approach can be adapted to analyse various diets, in numerous conservation and ecological applications. It is consequently appropriate for detecting fine dietary variations among individuals, populations and species as well as for the identification of rare food items.

Multispecies models reveal that eDNA metabarcoding is more sensitive than backpack electrofishing for conducting fish surveys in freshwater streams

Environmental DNA (eDNA) sampling can provide accurate, cost-effective, landscape-level data on species distributions. Previous studies have compared the sensitivity of eDNA sampling to traditional sampling methods for single species, but similar comparative studies on multispecies eDNA metabarcoding are rare. Using hierarchical site occupancy detection models, we examined whether key choices associated with eDNA metabarcoding (primer selection, low-abundance read filtering and the number of positive water samples used to classify a species as present at a site) affect the sensitivity of metabarcoding, relative to backpack electrofishing for fish in freshwater streams. Under all scenarios (teleostei and vertebrate primers; 0%, 0.1% and 1% read filtering thresholds; one or two positive samples required to classify species as present), we found that eDNA metabarcoding is, on average, more sensitive than electrofishing. Combining vertebrate and teleostei markers resulted in higher detection probabilities relative to the use of either marker in isolation. Increasing the threshold used to filter low-abundance reads decreased species detection probabilities but did not change our overall finding that eDNA metabarcoding was more sensitive than electrofishing. Using a threshold of two positive water samples (out of five) to classify a species as present typically had negligible effects on detection probabilities compared to using one positive water sample. Our findings demonstrate that eDNA metabarcoding is generally more sensitive than electrofishing for conducting fish surveys in freshwater streams, and that this outcome is not sensitive to methodological decisions associated with metabarcoding.

Thermal stratification and fish thermal preference explain vertical eDNA distributions in lakes

Significant advances have been made towards surveying animal and plant communities using DNA isolated from environmental samples. Despite rapid progress, we lack a comprehensive understanding of the "ecology" of environmental DNA (eDNA), particularly its temporal and spatial distribution and how this is shaped by abiotic and biotic processes. Here, we tested how seasonal variation in thermal stratification and animal habitat preferences influences the distribution of eDNA in lakes. We sampled eDNA depth profiles of five dimictic lakes during both summer stratification and autumn turnover, each containing warm- and cool-water fishes as well as the cold-water stenotherm, lake trout (Salvelinus namaycush). Habitat use by S. namaycush was validated by acoustic telemetry and was significantly related to eDNA distribution during stratification. Fish eDNA became "stratified" into layers during summer months, reflecting lake stratification and the thermal niches of the species. During summer months, S. namaycush, which rarely ventured into shallow waters, could only be detected at the deepest layers of the lakes, whereas the eDNA of warm-water fishes was much more abundant above the thermocline. By contrast, during autumn lake turnover, the fish species assemblage as detected by eDNA was homogenous throughout the water column. These findings contribute to our overall understanding of the "ecology" of eDNA within lake ecosystems, illustrating how the strong interaction between seasonal thermal structure in lakes and thermal niches of species on very localized spatial scales influences our ability to detect species.

Longer is Not Always Better: Optimizing Barcode Length for Large-Scale Species Discovery and Identification

New techniques for the species-level sorting of millions of specimens are needed in order to accelerate species discovery, determine how many species live on earth, and develop efficient biomonitoring techniques. These sorting methods should be reliable, scalable, and cost-effective, as well as being largely insensitive to low-quality genomic DNA, given that this is usually all that can be obtained from museum specimens. Mini-barcodes seem to satisfy these criteria, but it is unclear how well they perform for species-level sorting when compared with full-length barcodes. This is here tested based on 20 empirical data sets covering ca. 30,000 specimens (5500 species) and six clade-specific data sets from GenBank covering ca. 98,000 specimens ($>$20,000 species). All specimens in these data sets had full-length barcodes and had been sorted to species-level based on morphology. Mini-barcodes of different lengths and positions were obtained in silico from full-length barcodes using a sliding window approach (three windows: 100 bp, 200 bp, and 300 bp) and by excising nine mini-barcodes with established primers (length: 94-407 bp). We then tested whether barcode length and/or position reduces species-level congruence between morphospecies and molecular operational taxonomic units (mOTUs) that were obtained using three different species delimitation techniques (Poisson Tree Process, Automatic Barcode Gap Discovery, and Objective Clustering). Surprisingly, we find no significant differences in performance for both species- or specimen-level identification between full-length and mini-barcodes as long as they are of moderate length ($>$200 bp). Only very short mini-barcodes (<200 bp) perform poorly, especially when they are located near the 5$^\prime$ end of the Folmer region. The mean congruence between morphospecies and mOTUs was ca. 75% for barcodes $>$200 bp and the congruent mOTUs contain ca. 75% of all specimens. Most conflict is caused by ca. 10% of the specimens that can be identified and should be targeted for re-examination in order to efficiently resolve conflict. Our study suggests that large-scale species discovery, identification, and metabarcoding can utilize mini-barcodes without any demonstrable loss of information compared to full-length barcodes. [DNA barcoding; metabarcoding; mini-barcodes; species discovery.].

Description of Neoperla mindoroensis sp. nov., the first record of a stonefly from Mindoro, Philippines (Plecoptera, Perlidae), and identification of its life stages using COI barcodes

The new stonefly species, Neoperlamindoroensissp. nov. (Perlidae), from Mindoro island is described. The new species is assigned to the N.recta species complex of the N.montivaga group on account of its obvious T7 and T8 with pointed processes and the presence of basolateral lobes in the everted aedeagal sac. The male adult is distinguishable by its aedeagus with a slightly raised mediodorsal lobe, fully covered with fine spinules, while the female adult has comparably small eggs (240 × 220 μm) with a punctate, chorionic surface with punctae arranged in polygonal FCIs. The life stages and sexes were assigned using COI mtDNA barcodes (2.2% maximum intraspecific genetic distance), which were compared with available barcodes of congeners, which had interspecific genetic distances varying by at least 23.5%. Biogeographic aspects, ecological habitat requirements, and suitability as potential bioindicator of the species are also briefly discussed.

Biodiversity between sand grains: Meiofauna composition across southern and western Sweden assessed by metabarcoding

The meiofauna is an important part of the marine ecosystem, but its composition and distribution patterns are relatively unexplored. Here we assessed the biodiversity and community structure of meiofauna from five locations on the Swedish western and southern coasts using a high-throughput DNA sequencing (metabarcoding) approach. The mitochondrial cytochrome oxidase 1 (COI) mini-barcode and nuclear 18S small ribosomal subunit (18S) V1-V2 region were amplified and sequenced using Illumina MiSeq technology. Our analyses revealed a higher number of species than previously found in other areas: thirteen samples comprising 6.5 dm3 sediment revealed 708 COI and 1,639 18S metazoan OTUs. Across all sites, the majority of the metazoan biodiversity was assigned to Arthropoda, Nematoda and Platyhelminthes. Alpha and beta diversity measurements showed that community composition differed significantly amongst sites. OTUs initially assigned to Acoela, Gastrotricha and the two Platyhelminthes sub-groups Macrostomorpha and Rhabdocoela were further investigated and assigned to species using a phylogeny-based taxonomy approach. Our results demonstrate that there is great potential for discovery of new meiofauna species even in some of the most extensively studied locations.

The effects of spatial and temporal replicate sampling on eDNA metabarcoding

Background

The heterogeneous nature of environmental DNA (eDNA) and its effects on species detection and community composition estimates has been highlighted in several studies in the past decades. Mostly in the context of spatial distribution over large areas, in fewer occasions looking at spatial distribution within a single body of water. Temporal variation of eDNA, similarly, has mostly been studied as seasonality, observing changes over large periods of time, and often only for small groups of organisms such as fish and amphibians.

Methods

We analyzed and compared small-scale spatial and temporal variation by sampling eDNA from two small, isolated dune lakes for 20 consecutive weeks. Metabarcoding was performed on the samples using generic COI primers. Molecular operational taxonomic unit (MOTUs) were used to assess dissimilarities between spatial and temporal replicates.

Results

Our results show large differences between samples taken within one lake at one point in time, but also expose the large differences between temporal replicates, even those taken only 1 week apart. Furthermore, between-site dissimilarities showed a linear correlation with time frame, indicating that between-site differences will be inflated when samples are taken over a period of time. We also assessed the effects of PCR replicates and processing strategies on general patterns of dissimilarity between samples. While more inclusive PCR replicate strategies lead to higher richness estimations, dissimilarity patterns between samples did not significantly change.

Conclusions

We conclude that the dissimilarity of temporal replicates at a 1 week interval is comparable to that of spatial replicate samples. It increases, however, for larger time intervals, which suggests that population turnover effects can be stronger than community heterogeneity. Spatial replicates alone may not be enough for optimal recovery of taxonomic diversity, and cross-comparisons of different locations are susceptible to inflated dissimilarities when performed over larger time intervals. Many of the observed MOTUs could be classified as either phyto- or zooplankton, two groups that have gained traction in recent years as potential novel bio-indicator species. Our results, however, indicate that these groups might be susceptible to large community shifts in relatively short periods of time, highlighting the need to take temporal variations into consideration when assessing their usability as water quality indicators.

Measuring global fish species richness with eDNA metabarcoding

Despite mounting threats to global freshwater and marine biodiversity, including climate change, habitat alteration, overharvesting and pollution, we struggle to know which species are present below the water's surface that are suffering from these stressors. However, the idea that a water sample containing environmental DNA (eDNA) can be screened using high-throughput sequencing and bioinformatics to reveal the identity of aquatic species is a revolutionary advance for studying the patterns of species extirpation, invasive species establishment and the dynamics of species richness. To date, many of the critical tests of fisheries diversity using this metabarcoding approach have been conducted in lower diversity systems (<40 fish species), but in this issue of Molecular Ecology Resources, Cilleros et al. (2018) described their eDNA application in the species-rich French Guiana fishery (>200 fish species) and showed the greater potential and some limitations of using eDNA in species-rich environments.

eDNA as a tool for identifying freshwater species in sustainable forestry: A critical review and potential future applications

Environmental DNA (eDNA) is an emerging biological monitoring tool that can aid in assessing the effects of forestry and forest manufacturing activities on biota. Monitoring taxa across broad spatial and temporal scales is necessary to ensure forest management and forest manufacturing activities meet their environmental goals of maintaining biodiversity. Our objectives are to describe potential applications of eDNA across the wood products supply chain extending from regenerating forests, harvesting, and wood transport, to manufacturing facilities, and to review the current state of the science in this context. To meet our second objective, we summarize the taxa examined with targeted (PCR, qPCR or ddPCR) or metagenomic eDNA methods (eDNA metabarcoding), evaluate how estimated species richness compares between traditional field sampling and eDNA metabarcoding approaches, and compare the geographical representation of prior eDNA studies in freshwater ecosystems to global wood baskets. Potential applications of eDNA include evaluating the effects of forestry and forest manufacturing activities on aquatic biota, delineating fish-bearing versus non fish-bearing reaches, evaluating effectiveness of constructed road crossings for freshwater organism passage, and determining the presence of at-risk species. Studies using targeted eDNA approaches focused on fish, amphibians, and invertebrates, while metagenomic studies focused on fish, invertebrates, and microorganisms. Rare, threatened, or endangered species received the least attention in targeted eDNA research, but are arguably of greatest interest to sustainable forestry and forest manufacturing that seek to preserve freshwater biodiversity. Ultimately, using eDNA methods will enable forestry and forest manufacturing managers to have data-driven prioritization for conservation actions for all freshwater species.

NGS barcoding reveals high resistance of a hyperdiverse chironomid (Diptera) swamp fauna against invasion from adjacent freshwater reservoirs

BACKGROUND

Macroinvertebrates such as non-biting midges (Chironomidae: Diptera) are important components of freshwater ecosystems. However, they are often neglected in biodiversity and conservation research because invertebrate species richness is difficult and expensive to quantify with traditional methods. We here demonstrate that Next Generation Sequencing barcodes ("NGS barcodes") can provide relief because they allow for fast and large-scale species-level sorting of large samples at low cost.

RESULTS

We used NGS barcoding to investigate the midge fauna of Singapore's swamp forest remnant (Nee Soon Swamp Forest). Based on > 14.000 barcoded specimens, we find that the swamp forest maintains an exceptionally rich fauna composed of an observed number of 289 species (estimated 336 species) in a very small area (90 ha). We furthermore barcoded the chironomids from three surrounding reservoirs that are located in close proximity. Although the swamp forest remnant is much smaller than the combined size of the freshwater reservoirs in the study (90 ha vs. > 450 ha), the latter only contains 33 (estimated 61) species. We show that the resistance of the swamp forest species assemblage is high because only 8 of the 314 species are shared despite the close proximity. Moreover, shared species are not very abundant (3% of all specimens). A redundancy analysis revealed that ~ 21% of the compositional variance of midge communities within the swamp forest was explained by a range of variables with conductivity, stream order, stream width, temperature, latitude (flow direction), and year being significant factors influencing community structure. An LME analysis demonstrates that the total species richness decreased with increasing conductivity.

CONCLUSION

Our study demonstrates that midge diversity of a swamp forest can be so high that it questions global species diversity estimates for Chironomidae, which are an important component of many freshwater ecosystems. We furthermore demonstrate that small and natural habitat remnants can have high species turnover and can be very resistant to the invasion of species from neighboring reservoirs. Lastly, the study shows how NGS barcodes can be used to integrate specimen- and species-rich invertebrate taxa in biodiversity and conservation research.

Evaluating freshwater macroinvertebrates from eDNA metabarcoding: A river Nalón case study

Rivers are a vital resource for human wellbeing. To reduce human impact on water bodies, the European Union has established an essential regulatory framework for protection and sustainable management (WFD; 2000/60/EC). In this strategy, reliable and economic bioindicators are a fundamental component. Benthic macroinvertebrates are the group most commonly used as bioindicators through all European countries. However, their conventional assessment currently entails serious cost-efficiency limitations. In this study, we have tested the reliability of metabarcoding as a tool to record river macroinvertebrates using samples from a mock community (in vitro validation) and eDNA extracted for field validation from water from six sites within a north Iberian river (River Nalón, Asturias, Spain). Two markers (V4 region within the nuclear 18S rDNA and a fragment of the mitochondrial COI gene) were amplified and sequenced using an Illumina platform. The molecular technique has proven to be more sensitive than the visual one. A cost-benefit analysis shows that the metabarcoding approach is more expensive than conventional techniques for determining macroinvertebrate communities but requires fewer sampling and identification efforts. Our results suggest metabarcoding is a useful tool for alternative assessment of freshwater quality.

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 reveals quantitative patterns of fish biodiversity in large rivers despite its downstream transportation

Despite the ecological and societal importance of large rivers, fish sampling remains costly and limited to specific habitats (e.g., river banks). Using an eDNA metabarcoding approach, we regularly sampled 500 km of a large river (Rhône River). Comparisons with long-term electrofishing surveys demonstrated the ability of eDNA metabarcoding to qualitatively and quantitatively reveal fish assemblage structures (relative species abundance) but eDNA integrated a larger space than the classical sampling location. Combination of a literature review and field data showed that eDNA behaves in the water column like fine particulate organic matter. Its detection distance varied from a few km in a small stream to more than 100 km in a large river. To our knowledge, our results are the first demonstration of the capacity of eDNA metabarcoding to describe longitudinal fish assemblage patterns in a large river, and metabarcoding appears to be a reliable, cost-effective method for future monitoring.

Sorting specimen-rich invertebrate samples with cost-effective NGS barcodes: Validating a reverse workflow for specimen processing

Biologists frequently sort specimen-rich samples to species. This process is daunting when based on morphology, and disadvantageous if performed using molecular methods that destroy vouchers (e.g., metabarcoding). An alternative is barcoding every specimen in a bulk sample and then presorting the specimens using DNA barcodes, thus mitigating downstream morphological work on presorted units. Such a "reverse workflow" is too expensive using Sanger sequencing, but we here demonstrate that is feasible with an next-generation sequencing (NGS) barcoding pipeline that allows for cost-effective high-throughput generation of short specimen-specific barcodes (313 bp of COI; laboratory cost <$0.50 per specimen) through next-generation sequencing of tagged amplicons. We applied our approach to a large sample of tropical ants, obtaining barcodes for 3,290 of 4,032 specimens (82%). NGS barcodes and their corresponding specimens were then sorted into molecular operational taxonomic units (mOTUs) based on objective clustering and Automated Barcode Gap Discovery (ABGD). High diversity of 88-90 mOTUs (4% clustering) was found and morphologically validated based on preserved vouchers. The mOTUs were overwhelmingly in agreement with morphospecies (match ratio 0.95 at 4% clustering). Because of lack of coverage in existing barcode databases, only 18 could be accurately identified to named species, but our study yielded new barcodes for 48 species, including 28 that are potentially new to science. With its low cost and technical simplicity, the NGS barcoding pipeline can be implemented by a large range of laboratories. It accelerates invertebrate species discovery, facilitates downstream taxonomic work, helps with building comprehensive barcode databases and yields precise abundance information.

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.