Estimating fish alpha- and beta-diversity along a small stream with environmental DNA metabarcoding

An integrated spatio-temporal view of riverine biodiversity using environmental DNA metabarcoding

Anthropogenically forced changes in global freshwater biodiversity demand more efficient monitoring approaches. Consequently, environmental DNA (eDNA) analysis is enabling ecosystem-scale biodiversity assessment, yet the appropriate spatio-temporal resolution of robust biodiversity assessment remains ambiguous. Here, using intensive, spatio-temporal eDNA sampling across space (five rivers in Europe and North America, with an upper range of 20-35 km between samples), time (19 timepoints between 2017 and 2018) and environmental conditions (river flow, pH, conductivity, temperature and rainfall), we characterise the resolution at which information on diversity across the animal kingdom can be gathered from rivers using eDNA. In space, beta diversity was mainly dictated by turnover, on a scale of tens of kilometres, highlighting that diversity measures are not confounded by eDNA from upstream. Fish communities showed nested assemblages along some rivers, coinciding with habitat use. Across time, seasonal life history events, including salmon and eel migration, were detected. Finally, effects of environmental conditions were taxon-specific, reflecting habitat filtering of communities rather than effects on DNA molecules. We conclude that riverine eDNA metabarcoding can measure biodiversity at spatio-temporal scales relevant to species and community ecology, demonstrating its utility in delivering insights into river community ecology during a time of environmental change.

Use of passive sampling in environmental DNA metabarcoding technology: Monitoring of fish diversity in the Jiangmen coastal waters

Environmental DNA (eDNA) metabarcoding technology is proving to be the most up-to-date and promising method for monitoring marine fish diversity. Fish eDNA is usually collected on a filter membrane after the filtration of water. Not only does this require the use of specialized equipment, but the amount of filtered water needed is also difficult to meet. The recently proposed passive eDNA collection method can expand the sampling scale, providing new perspectives for monitoring marine biodiversity. The role of collection methods in eDNA surveys, however, remains unclear. In this study, a low-cost custom framework with two types of filter membrane materials was used to conduct passive submersion samplings at the north and south ends of Shangchuan Island, Jiangmen, China. After defined periods of submersion, the filter membranes were recovered and eDNA extracted. Metabarcoding techniques were applied to detect fish species information in the eDNA samples. A total of 106 marine fish species from 27 orders, 53 families, and 92 genera, including one cartilaginous fish, were identified in the samples. The majority of fish detected by active filtration were also found in the passively collected samples, within the same location. Both sampling methods, therefore, showed similar species richness. Passive sampling was effective in identifying fish species diversity and provided a higher spatial resolution owing to the sample replicates. Passive sampling was also more sensitive in detecting species that differ significantly in abundance (biomarkers) between different sampling depths. When active filtration is not possible, or when large-scale sampling is the purpose of the study, passive sampling methods certainly provide a promising alternative. The findings of our study provide guidance for fish surveys and continuous bio-stereoscopic monitoring in coastal waters.

First pelagic fish biodiversity assessment of Cosmonaut Sea based on environmental DNA

The Cosmonaut Sea is a typical marginal sea in East Antarctica that has not yet been greatly impacted by climate change. As one of the least explored areas in the Southern Ocean, our knowledge regarding its fish taxonomy and diversity has been sparse. eDNA metabarcoding, as an emerging and promising tool for marine biodiversity research and monitoring, has been widely used across taxa and habitats. During the 38th Chinese Antarctic Research Expedition (CHINARE-38), we collected seawater and surface sediment samples from 38 stations in the Cosmonaut Sea and performed the first, to our knowledge, eDNA analysis of fish biodiversity in the Southern Ocean based on the molecular markers of 12S rRNA and 16S rRNA. There were 48 fish species detected by the two markers in total, with 30 and 34 species detected by the 12S rRNA and 16S rRNA marker, respectively. This was more than the trawling results (19 species) and historical survey records (16 species, "BROKE-West" cruise). With some nonsignificant differences between the Gunnerus Ridge and the Oceanic Area of Enderby Land, the Cosmonaut Sea had a richer fish biodiversity in this research compared with previous studies, and its overall composition and distribution patterns were consistent with what we know in East Antarctica. We also found that the eDNA composition of fish in the Cosmonaut Sea might be related to some environmental factors. Our study demonstrated that the use of the eDNA technique for Antarctic fish biodiversity research is likely to yield more information with less sampling effort than traditional methods. In the context of climate change, the eDNA approach will provide a novel and powerful tool that is complementary to traditional methods for polar ecology research.

Integrating invasive species risk assessment into environmental DNA metabarcoding reference libraries

Environmental DNA (eDNA) metabarcoding has shown promise as a tool for estimating biodiversity and early detection of invasive species. In aquatic systems, advantages of this method include the ability to concurrently monitor biodiversity and detect incipient invasions simply through the collection and analysis of water samples. However, depending on the molecular markers chosen for a given study, reference libraries containing target sequences from present species may limit the usefulness of eDNA metabarcoding. To explore the extent of this issue and how it may be resolved to aid biodiversity and invasive species early detection goals, we focus on fishes in the well-studied Laurentian Great Lakes region. First, we provide a synthesis of species currently known from the region and of non-indigenous species identified as threats by international, national, regional, and introduction pathway-specific fish risk assessments. With these species lists, we then evaluate 23 primer pairs commonly used in fish eDNA metabarcoding with available databases of sequence coverage and species specificity. Finally, we identify established and potentially invasive non-indigenous fish that should be prioritized for genetic sequencing to ensure robust eDNA metabarcoding for the region. Our results should increase confidence in using eDNA metabarcoding for fisheries conservation and management in the Great Lakes region and help prioritize reference sequencing efforts. The ultimate utility of eDNA metabarcoding approaches will come when conservation management of existing fish communities is integrated with early detection efforts for invasive species surveillance to assess total fish biodiversity.

Environmental DNA (eDNA) Metabarcoding in the Fish Market and Nearby Seafood Restaurants in Taiwan Reveals the Underestimation of Fish Species Diversity in Seafood

Seafood, especially the traditional one in Taiwan, is rarely sourced from a fixed species and routinely from similar species depending on their availability. Hence, the species composition of seafood can be complicated. While a DNA-based approach has been routinely utilized for species identification, a large scale of seafood identification in fish markets and restaurants could be challenging (e.g., elevated cost and time-consuming only for a limited number of species identification). In the present study, we aimed to identify the majority of fish species potentially consumed in fish markets and nearby seafood restaurants using environmental DNA (eDNA) metabarcoding. Four eDNA samplings from a local fish market and nearby seafood restaurants were conducted using Sterivex cartridges. Nineteen universal primers previously validated for fish species identification were utilized to amplify the fragments of mitochondrial DNA (12S, COI, ND5) of species in eDNA samples and sequenced with NovaSeq 6000 sequencing. A total of 153 fish species have been identified based on 417 fish related operational taxonomic units (OTUs) generated from 50,534,995 reads. Principal Coordinate Analysis (PCoA) further showed the differences in fish species between the sampling times and sampling sites. Of these fish species, 22 chondrichthyan fish, 14 Anguilliformes species, and 15 Serranidae species were respectively associated with smoked sharks, braised moray eels, and grouper fish soups. To our best knowledge, this work represents the first study to demonstrate the feasibility of a large scale of seafood identification using eDNA metabarcoding approach. Our findings also imply the species diversity in traditional seafood might be seriously underestimated and crucial for the conservation and management of marine resources.

Plant biodiversity assessment through pollen DNA metabarcoding in Natura 2000 habitats (Italian Alps)

Monitoring biodiversity is of increasing importance in natural ecosystems. Metabarcoding can be used as a powerful molecular tool to complement traditional biodiversity monitoring, as total environmental DNA can be analyzed from complex samples containing DNA of different origin. The aim of this research was to demonstrate the potential of pollen DNA metabarcoding using the chloroplast trnL partial gene sequencing to characterize plant biodiversity. Collecting airborne biological particles with gravimetric Tauber traps in four Natura 2000 habitats within the Natural Park of Paneveggio Pale di San Martino (Italian Alps), at three-time intervals in 1 year, metabarcoding identified 68 taxa belonging to 32 local plant families. Metabarcoding could identify with finer taxonomic resolution almost all non-rare families found by conventional light microscopy concurrently applied. However, compared to microscopy quantitative results, Poaceae, Betulaceae, and Oleaceae were found to contribute to a lesser extent to the plant biodiversity and Pinaceae were more represented. Temporal changes detected by metabarcoding matched the features of each pollen season, as defined by aerobiological studies running in parallel, and spatial heterogeneity was revealed between sites. Our results showcase that pollen metabarcoding is a promising approach in detecting plant species composition which could provide support to continuous monitoring required in Natura 2000 habitats for biodiversity conservation.

Beyond fish eDNA metabarcoding: Field replicates disproportionately improve the detection of stream associated vertebrate species

Fast, reliable, and comprehensive biodiversity monitoring data are needed for environmental decision making and management. Recent work on fish environmental DNA (eDNA) metabarcoding shows that aquatic diversity can be captured fast, reliably, and non-invasively at moderate costs. Because water in a catchment flows to the lowest point in the landscape, often a stream, it can collect traces of terrestrial species via surface or subsurface runoff along its way or when specimens come into direct contact with water (e.g., when drinking). Thus, fish eDNA metabarcoding data can provide information on fish but also on other vertebrate species that live in riparian habitats. This additional data may offer a much more comprehensive approach for assessing vertebrate diversity at no additional costs. Studies on how the sampling strategy affects species detection especially of stream-associated communities, however, are scarce. We therefore performed an analysis on the effects of biological replication on both fish as well as (semi-)terrestrial species detection. Along a 2 km stretch of the river Mulde (Germany), we collected 18 1-L water samples and analyzed the relation of detected species richness and quantity of biological replicates taken. We detected 58 vertebrate species, of which 25 were fish and lamprey, 18 mammals, and 15 birds, which account for 50%, 22.2%, and 7.4% of all native species to the German federal state of Saxony-Anhalt. However, while increasing the number of biological replicates resulted in only 24.8% more detected fish and lamprey species, mammal, and bird species richness increased disproportionately by 68.9% and 77.3%, respectively. Contrary, PCR replicates showed little stochasticity. We thus emphasize to increase the number of biological replicates when the aim is to improve general species detections. This holds especially true when the focus is on rare aquatic taxa or on (semi-)terrestrial species, the so-called ‘bycatch’. As a clear advantage, this information can be obtained without any additional sampling or laboratory effort when the sampling strategy is chosen carefully. With the increased use of eDNA metabarcoding as part of national fish bioassessment and monitoring programs, the complimentary information provided on bycatch can be used for biodiversity monitoring and conservation on a much broader scale.

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.

Composition and distribution of fish environmental DNA in an Adirondack watershed

Background

Environmental DNA (eDNA) surveys are appealing options for monitoring aquatic biodiversity. While factors affecting eDNA persistence, capture and amplification have been heavily studied, watershed-scale surveys of fish communities and our confidence in such need further exploration.

Methods

We characterized fish eDNA compositions using rapid, low-volume filtering with replicate and control samples scaled for a single Illumina MiSeq flow cell, using the mitochondrial 12S ribosomal RNA locus for taxonomic profiling. Our goals were to determine: (1) spatiotemporal variation in eDNA abundance, (2) the filtrate needed to achieve strong sequencing libraries, (3) the taxonomic resolution of 12S ribosomal sequences in the study environment, (4) the portion of the expected fish community detectable by 12S sequencing, (5) biases in species recovery, (6) correlations between eDNA compositions and catch per unit effort (CPUE) and (7) the extent that eDNA profiles reflect major watershed features. Our bioinformatic approach included (1) estimation of sequencing error from unambiguous mappings and simulation of taxonomic assignment error under various mapping criteria; (2) binning of species based on inferred assignment error rather than by taxonomic rank; and (3) visualization of mismatch distributions to facilitate discovery of distinct haplotypes attributed to the same reference. Our approach was implemented within the St. Regis River, NY, USA, which supports tribal and recreational fisheries and has been a target of restoration activities. We used a large record of St. Regis-specific observations to validate our assignments.

Results

We found that 300 mL drawn through 25-mm cellulose nitrate filters yielded greater than 5 ng/µL DNA at most sites in summer, which was an approximate threshold for generating strong sequencing libraries in our hands. Using inferred sequence error rates, we binned 12S references for 110 species on a state checklist into 85 single-species bins and seven multispecies bins. Of 48 bins observed by capture survey in the St. Regis, we detected eDNA consistent with 40, with an additional four detections flagged as potential contaminants. Sixteen unobserved species detected by eDNA ranged from plausible to implausible based on distributional data, whereas six observed species had no 12S reference sequence. Summed log-ratio compositions of eDNA-detected taxa correlated with log(CPUE) (Pearson’s R = 0.655, P < 0.001). Shifts in eDNA composition of several taxa and a genotypic shift in channel catfish (Ictalurus punctatus) coincided with the Hogansburg Dam, NY, USA. In summary, a simple filtering apparatus operated by field crews without prior expertise gave useful summaries of eDNA composition with minimal evidence of field contamination. 12S sequencing achieved useful taxonomic resolution despite the short marker length, and data exploration with standard bioinformatic tools clarified taxonomic uncertainty and sources of error.

Space-time dynamics in monitoring neotropical fish communities using eDNA metabarcoding

The biodiverse Neotropical ecoregion remains insufficiently assessed, poorly managed, and threatened by unregulated human activities. Novel, rapid and cost-effective DNA-based approaches are valuable to improve understanding of the biological communities and for biomonitoring in remote areas. Here, we evaluate the potential of environmental DNA (eDNA) metabarcoding for assessing the structure and distribution of fish communities by analysing water and sediment from 11 locations along the Jequitinhonha River catchment (Brazil). Each site was sampled twice, before and after a major rain event in a five-week period and fish diversity was estimated using high-throughput sequencing of 12S rRNA amplicons. In total, 252 Molecular Operational Taxonomic Units (MOTUs) and 34 fish species were recovered, including endemic, introduced, and previously unrecorded species for this basin. Spatio-temporal variation of eDNA from fish assemblages was observed and species richness was nearly twice as high before the major rain event compared to afterwards. Yet, peaks of diversity were primarily associated with only four of the locations. No correlation between β-diversity and longitudinal distance or presence of dams was detected, but low species richness observed at sites located near dams might that these anthropogenic barriers may have an impact on local fish diversity. Unexpectedly high α-diversity levels recorded at the river mouth suggest that these sections should be further evaluated as putative "eDNA reservoirs" for rapid monitoring. By uncovering spatio-temporal changes, unrecorded biodiversity components, and putative anthropogenic impacts on fish assemblages, we further strengthen the potential of eDNA metabarcoding as a biomonitoring tool, especially in regions often neglected or difficult to access.

Gradients in richness and turnover of a forest passerine's diet prior to breeding: A mixed model approach applied to faecal metabarcoding data

Rather little is known about the dietary richness and variation of generalist insectivorous species, including birds, due primarily to difficulties in prey identification. Using faecal metabarcoding, we provide the most comprehensive analysis of a passerine's diet to date, identifying the relative magnitudes of biogeographic, habitat and temporal trends in the richness and turnover in diet of Cyanistes caeruleus (blue tit) along a 39 site and 2° latitudinal transect in Scotland. Faecal samples were collected in 2014-2015 from adult birds roosting in nestboxes prior to nest building. DNA was extracted from 793 samples and we amplified COI and 16S minibarcodes. We identified 432 molecular operational taxonomic units that correspond to putative dietary items. Most dietary items were rare, with Lepidoptera being the most abundant and taxon-rich prey order. Here, we present a statistical approach for estimation of gradients and intersample variation in taxonomic richness and turnover using a generalised linear mixed model. We discuss the merits of this approach over existing tools and present methods for model-based estimation of repeatability, taxon richness and Jaccard indices. We found that dietary richness increases significantly as spring advances, but changes little with elevation, latitude or local tree composition. In comparison, dietary composition exhibits significant turnover along temporal and spatial gradients and among sites. Our study shows the promise of faecal metabarcoding for inferring the macroecology of food webs, but we also highlight the challenge posed by contamination and make recommendations of laboratory and statistical practices to minimise its impact on inference.

Standards for Methods Utilizing Environmental DNA for Detection of Fish Species

Environmental DNA (eDNA) techniques are gaining attention as cost-effective, non-invasive strategies for acquiring information on fish and other aquatic organisms from water samples. Currently, eDNA approaches are used to detect specific fish species and determine fish community diversity. Various protocols used with eDNA methods for aquatic organism detection have been reported in different eDNA studies, but there are no general recommendations for fish detection. Herein, we reviewed 168 papers to supplement and highlight the key criteria for each step of eDNA technology in fish detection and provide general suggestions for eliminating detection errors. Although there is no unified recommendation for the application of diverse eDNA in detecting fish species, in most cases, 1 or 2 L surface water collection and eDNA capture on 0.7-μm glass fiber filters followed by extraction with a DNeasy Blood and Tissue Kit or PowerWater DNA Isolation Kit are useful for obtaining high-quality eDNA. Subsequently, species-specific quantitative polymerase chain reaction (qPCR) assays based on mitochondrial cytochrome b gene markers or eDNA metabarcoding based on both 12S and 16S rRNA markers via high-throughput sequencing can effectively detect target DNA or estimate species richness. Furthermore, detection errors can be minimized by mitigating contamination, negative control, PCR replication, and using multiple genetic markers. Our aim is to provide a useful strategy for fish eDNA technology that can be applied by researchers, advisors, and managers.

Particle Size Distribution of Environmental DNA from the Nuclei of Marine Fish

Environmental DNA (eDNA) analyses have enabled a more efficient surveillance of species distribution and composition than conventional methods. However, the characteristics and dynamics of eDNA (e.g., origin, state, transport, and fate) remain unknown. This is especially limited for the eDNA derived from nuclei (nu-eDNA), which has recently been used in eDNA analyses. Here, we compared the particle size distribution (PSD) of nu-eDNA from Japanese Jack Mackerel (Trachurus japonicus) with that of mt-eDNA (eDNA derived from mitochondria) reported in previous studies. We repeatedly sampled rearing water from the tanks under multiple temperatures and fish biomass levels, and quantified the copy numbers of size-fractioned nu-eDNA. We found that the concentration of nu-eDNA was higher than that of mt-eDNA at 3-10 μm size fraction. Moreover, at the 0.8-3 μm and 0.4-0.8 μm size fractions, eDNA concentrations of both types increased with higher temperature and their degradation tended to be suppressed. These results imply that the production of eDNA from large to small size fractions could buffer the degradation of small-sized eDNA, which could improve its persistence in water. Our findings will contribute to refine the difference between nu- and mt-eDNA properties, and assist eDNA analyses as an efficient tool for the conservation of aquatic species.

Understanding PCR Processes to Draw Meaningful Conclusions from Environmental DNA Studies

As environmental DNA (eDNA) studies have grown in popularity for use in ecological applications, it has become clear that their results differ in significant ways from those of traditional, non-PCR-based surveys. In general, eDNA studies that rely on amplicon sequencing may detect hundreds of species present in a sampled environment, but the resulting species composition can be idiosyncratic, reflecting species' true biomass abundances poorly or not at all. Here, we use a set of simulations to develop a mechanistic understanding of the processes leading to the kinds of results common in mixed-template PCR-based (metabarcoding) studies. In particular, we focus on the effects of PCR cycle number and primer amplification efficiency on the results of diversity metrics in sequencing studies. We then show that proportional indices of amplicon reads capture trends in taxon biomass with high accuracy, particularly where amplification efficiency is high (median correlation up to 0.97). Our results explain much of the observed behavior of PCR-based studies, and lead to recommendations for best practices in the field.

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.

Monitoring riverine fish communities through eDNA metabarcoding: determining optimal sampling strategies along an altitudinal and biodiversity gradient

Monitoring aquatic biodiversity through DNA extracted from environmental samples (eDNA) combined with high-throughput sequencing, commonly referred to as eDNA metabarcoding, is increasing in popularity within the scientific community. However, sampling strategies, laboratory protocols and analytical pipelines can influence the results of eDNA metabarcoding surveys. While the impact of laboratory protocols and analytical pipelines have been extensively studied, the importance of sampling strategies on eDNA metabarcoding surveys has not received the same attention. To avoid underestimating local biodiversity, adequate sampling strategies (i.e. sampling intensity and spatial sampling replication) need to be implemented. This study evaluated the impact of sampling strategies along an altitudinal and biodiversity gradient in the upper section of the Murrumbidgee River (Murray-Darling Basin, Australia). An eDNA metabarcoding survey was used to determine the local fish biodiversity and evaluate the influence of sampling intensity and spatial sampling replication on the biodiversity estimates. The results show that optimal eDNA sampling strategies varied between sites and indicate that river morphology, species richness and species abundance affect the optimal sampling intensity and spatial sampling replication needed to accurately assess the fish biodiversity. While the generality of the patterns will need to be confirmed through future studies, these findings provide a basis to guide future eDNA metabarcoding surveys in river systems.