The future of fish-based ecological assessment of European rivers: from traditional EU Water Framework Directive compliant methods to eDNA metabarcoding-based approaches

European river typologies fail to capture diatom, fish, and macrophyte community composition

Typology systems are frequently used in applied and fundamental ecology and are relevant for environmental monitoring and conservation. They aggregate ecosystems into discrete types based on biotic and abiotic variables, assuming that ecosystems of the same type are more alike than ecosystems of different types with regard to a specific property of interest. We evaluated whether this assumption is met by the Broad River Types (BRT), a recently proposed European river typology system, that classifies river segments based on abiotic variables, when it is used to group biological communities. We compiled data on the community composition of diatoms, fishes, and aquatic macrophytes throughout Europe and evaluated whether the composition is more similar in site groups with the same river type than in site groups of different river types using analysis of similarities, classification strength, typical species analysis, and the area under zeta diversity decline curves. We compared the performance of the BRT with those of four region-based typology systems, namely, Illies Freshwater Ecoregions, the Biogeographic Regions, the Freshwater Ecoregions of the World, and the Environmental Zones, as well as spatial autocorrelation (SA) classifications. All typology systems received low scores from most evaluation methods, relative to predefined thresholds and the SA classifications. The BRT often scored lowest of all typology systems. Within each typology system, community composition overlapped considerably between site groups defined by the types of the systems. The overlap tended to be the lowest for fishes and between Illies Freshwater Ecoregions. In conclusion, we found that existing broad-scale river typology systems fail to delineate site groups with distinct and compositionally homogeneous communities of diatoms, fishes, and macrophytes. A way to improve the fit between typology systems and biological communities might be to combine segment-based and region-based typology systems to simultaneously account for local environmental variation and historical distribution patterns, thus potentially improving the utility of broad-scale typology systems for freshwater biota.

eDNA testing reveals surprising findings on fish population dynamics in Thailand

COVID-19, a global health concern, has an effect on all aspects of the economy. The aquaculture and fishing industries were severely harmed as a result of the closures in multiple nations. Regular systems for inventory monitoring, production, and supply were disrupted. Cancellation of programmes for research, fieldwork, sampling, and tagging influences management-required data. For effective species management, fish dispersion assessments are indispensable. However, due to the difficulty of accessing sampling sites and the associated costs, there is frequently a lack of comprehensive information regarding the distribution and abundance of organisms. The COVID-19 prohibition made fish monitoring more problematic. Due to constant pressure, populations of the stone lapping minnow (Garra cambodgiensis), one of Thailand's overfished fish, are rapidly declining. Therefore, eDNA-based monitoring was devised and implemented to reveal the likely dispersal of the species in Thailand prior to and following the lockdown. At 28 locations within the Chao Phraya River Basin, water samples were collected. qPCR was used to determine the presence or absence of G. cambodgiensis in water samples. In 78 of 252 water samples, a wide range of computed copy numbers for G. cambodgiensis eDNA was observed. It was discovered that samples collected in 2021 (after the lockdown) contain a higher concentration of G. cambodgiensis eDNA than samples collected in 2018 or 2019 (prior to the lockdown). The closure appears to be a boon and may result in a substantial restocking of the fish we have studied. Overall, eDNA-based analysis is an extremely promising new survey instrument.

Disentangling the effects of different human disturbances on multifaceted biodiversity indices in freshwater fish

Evaluating the effects of anthropogenic pressures on several biodiversity metrics can inform the management and monitoring of biodiversity loss. However, the type of disturbances can lead to different responses in different metrics. In this study, we aimed at disentangling the effects of different types of anthropogenic disturbances on freshwater fish communities. We calculated diversity indices for 1109 stream fish communities across France by computing richness and evenness components for ecological, morphological, and phylogenetic diversity, and used null models to estimate standardized effect sizes. We used generalized linear mixed models to assess the relative effects of environmental and anthropogenic drivers in driving those diversity indices. Our results demonstrated that all diversity indices exhibited significant responses to both climatic conditions and anthropogenic disturbances. While we observed a decrease of ecological and phylogenetic richness with the intensity of disturbance, a weak increase in morphological richness and evenness was apparent. Overall, our results demonstrated the importance of disentangling various types of disturbances when assessing human-induced ecological impacts and highlighted that different facets of diversity are not impacted identically by anthropogenic disturbances in stream fish communities. This calls for further work seeking to integrate biodiversity responses to human disturbances into a multifaceted framework, and could have beneficial implications when planning conservation action in freshwater ecosystems.

Quantitative monitoring of diverse fish communities on a large scale combining eDNA metabarcoding and qPCR

Environmental DNA (eDNA) metabarcoding is an effective method for studying fish communities but allows only an estimation of relative species abundance (density/biomass). Here, we combine metabarcoding with an estimation of the total abundance of eDNA amplified by our universal marker (teleo) using a quantitative (q)PCR approach to infer the absolute abundance of fish species. We carried out a 2850-km eDNA survey within the Danube catchment using a spatial integrative sampling protocol coupled with traditional electrofishing for fish biomass and density estimation. Total fish eDNA concentrations and total fish abundance were highly correlated. The correlation between eDNA concentrations per taxon and absolute specific abundance was of comparable strength when all sites were pooled and remained significant when the sites were considered separately. Furthermore, a nonlinear mixed model showed that species richness was underestimated when the amount of teleo-DNA extracted from a sample was below a threshold of 0.65 × 106 copies of eDNA. This result, combined with the decrease in teleo-DNA concentration by several orders of magnitude with river size, highlights the need to increase sampling effort in large rivers. Our results provide a comprehensive description of longitudinal changes in fish communities and underline our combined metabarcoding/qPCR approach for biomonitoring and bioassessment surveys when a rough estimate of absolute species abundance is sufficient.

Fishing for fish environmental DNA: Ecological applications, methodological considerations, surveying designs, and ways forward

Vast global declines of freshwater and marine fish diversity and population abundance pose serious threats to both ecosystem sustainability and human livelihoods. Environmental DNA (eDNA)-based biomonitoring provides robust, efficient, and cost-effective assessment of species occurrences and population trends in diverse aquatic environments. Thus, it holds great potential for improving conventional surveillance frameworks to facilitate fish conservation and fisheries management. However, the many technical considerations and rapid developments underway in the eDNA arena can overwhelm researchers and practitioners new to the field. Here, we systematically analysed 416 fish eDNA studies to summarize research trends in terms of investigated targets, research aims, and study systems, and reviewed the applications, rationales, methodological considerations, and limitations of eDNA methods with an emphasis on fish and fisheries research. We highlighted how eDNA technology may advance our knowledge of fish behaviour, species distributions, population genetics, community structures, and ecological interactions. We also synthesized the current knowledge of several important methodological concerns, including the qualitative and quantitative power eDNA has to recover fish biodiversity and abundance, and the spatial and temporal representations of eDNA with respect to its sources. To facilitate ecological applications implementing fish eDNA techniques, recent literature was summarized to generate guidelines for effective sampling in lentic, lotic, and marine habitats. Finally, we identified current gaps and limitations, and pointed out newly emerging research avenues for fish eDNA. As methodological optimization and standardization improve, eDNA technology should revolutionize fish monitoring and promote biodiversity conservation and fisheries management that transcends geographic and temporal boundaries.

A strategy for successful integration of DNA-based methods in aquatic monitoring

Recent advances in molecular biomonitoring open new horizons for aquatic ecosystem assessment. Rapid and cost-effective methods based on organismal DNA or environmental DNA (eDNA) now offer the opportunity to produce inventories of indicator taxa that can subsequently be used to assess biodiversity and ecological quality. However, the integration of these new DNA-based methods into current monitoring practices is not straightforward, and will require coordinated actions in the coming years at national and international levels.

To plan and stimulate such an integration, the European network DNAqua-Net (COST Action CA15219) brought together international experts from academia, as well as key environmental biomonitoring stakeholders from different European countries. Together, this transdisciplinary consortium developed a roadmap for implementing DNA-based methods with a focus on inland waters assessed by the EU Water Framework Directive (2000/60/EC). This was done through a series of online workshops held in April 2020, which included fifty participants, followed by extensive synthesis work.

The roadmap is organised around six objectives: 1) to highlight the effectiveness and benefits of DNA-based methods, 2) develop an adaptive approach for the implementation of new methods, 3) provide guidelines and standards for best practice, 4) engage stakeholders and ensure effective knowledge transfer, 5) support the environmental biomonitoring sector to achieve the required changes, 6) steer the process and harmonise efforts at the European level.

This paper provides an overview of the forum discussions and the common European views that have emerged from them, while reflecting the diversity of situations in different countries. It highlights important actions required for a successful implementation of DNA-based biomonitoring of aquatic ecosystems by 2030.

Strategies for sample labelling and library preparation in DNA metabarcoding studies

Metabarcoding of DNA extracted from environmental or bulk specimen samples is increasingly used to profile biota in basic and applied biodiversity research because of its targeted nature that allows sequencing of genetic markers from many samples in parallel. To achieve this, PCR amplification is carried out with primers designed to target a taxonomically informative marker within a taxonomic group, and sample-specific nucleotide identifiers are added to the amplicons prior to sequencing. The latter enables assignment of the sequences back to the samples they originated from. Nucleotide identifiers can be added during the metabarcoding PCR and during "library preparation", that is, when amplicons are prepared for sequencing. Different strategies to achieve this labelling exist. All have advantages, challenges and limitations, some of which can lead to misleading results, and in the worst case compromise the fidelity of the metabarcoding data. Given the range of questions addressed using metabarcoding, ensuring that data generation is robust and fit for the chosen purpose is critically important for practitioners seeking to employ metabarcoding for biodiversity assessments. Here, we present an overview of the three main workflows for sample-specific labelling and library preparation in metabarcoding studies on Illumina sequencing platforms; one-step PCR, two-step PCR, and tagged PCR. Further, we distill the key considerations for researchers seeking to select an appropriate metabarcoding strategy for their specific study. Ultimately, by gaining insights into the consequences of different metabarcoding workflows, we hope to further consolidate the power of metabarcoding as a tool to assess biodiversity across a range of applications.

A Comparison of Fish Diversity in Rocky Reef Habitats by Multi-Mesh Gillnets and Environmental DNA Metabarcoding

This study was undertaken in order to explore the practical effectiveness of the environmental DNA (eDNA) metabarcoding approach in evaluating fish composition and diversity in a high heterogeneous rocky reef habitat. We assessed the fish composition and diversity characteristics of the rocky reef habitat at Dachen Islands, Taizhou and the Zhejiang Province in China in November 2020 by comparing two methods: multi-mesh gillnets and eDNA. A comparative analysis was carried out on the fish composition and diversity characteristics gained under the two methods by using taxonomy, ecotypes and diversity indices. The results showed that there were 28 species of fish collected through gillnets, distributed under 24 genera, 19 families, 6 orders and one class. Among them, 4, 18, and 6 species of near-surface, near groundfish and groundfish were found, respectively, with Thryssa mystax, Johnius belangerii, and Sebastiscus marmoratus being the dominant species in each water layer. A total of 81 species of fish detected by eDNA metabarcoding belonging to 67 genera, 46 families, 15 orders and 2 classes. The near-surface, near groundfish and groundfish species were 17, 42, and 22, with Thryssa vitrirostris, Benthosema pterotum, Harpadon nehereus, and Dasyatis akajei being the dominant species in each water layer. Twenty species (71.4%) and 41 species (50.6%) of reef fish were counted by gillnets and eDNA, respectively. The results showed that multi-mesh gillnets can accurately obtain information on fish composition in rocky reef habitats, but with some selectivity. The eDNA technology can detect species not collected by gillnets, but the number of species detected in areas with fast water velocity is significantly less than other eDNA stations where the water velocity is slow. In summary, the combination of traditional nets and eDNA will provide more information on taxonomic diversity and population biomass, transforming natural resource management and ecological studies of fish communities on a larger spatial and temporal scale.

Fish eDNA metabarcoding from aquatic biofilm samples: Methodological aspects

Fish eDNA metabarcoding is usually performed from filtered water samples. The volume of filtered water depends on the study scope and can rapidly become time consuming according to the number of samples that have to be processed. To avoid time allocated to filtration, passive DNA samplers have been used to recover fish eDNA from marine environments faster. In freshwater ecosystems, aquatic biofilms were used to catch eDNA from macroinvertebrates. Here, we test the capacity of aquatic biofilms to entrap fish eDNA in a large lake and, therefore, the possibility to perform fish eDNA metabarcoding from this matrix compared to the traditional fish eDNA approach from filtered water samples. Methodological aspects of the use of aquatic biofilms for fish eDNA metabarcoding (e.g. PCR replicates, biological replicates, bioinformatics pipeline, reference database and taxonomic assignment) were validated against a mock community. When using biofilms from habitats sheltered from wind and waves, biofilm and water approach provided similar inventories. Richness and diversity were comparable between both approaches. Approaches differed only for rare taxa. Our results illustrate the capacity of aquatic biofilms to act as passive eDNA samplers of fish eDNA and, therefore, the possibility to use biofilms to monitor fish communities efficiently from biofilms. Furthermore, our results open up avenues of research to study a diversity of biological groups (among which bioindicators as diatoms, macroinvertebrates and fish) from eDNA isolated from a single environmental matrix reducing sampling efforts, analysis time and costs.

The particle size distribution of environmental DNA varies with species and degradation

Environmental DNA (eDNA) analysis is frequently used as a non-invasive method to investigate species and biodiversity in ecosystems. However, such eDNA may represent both organisms currently present as well as species that released their DNA some point in the past, thereby representing a mix of current and historic biodiversity. This may lead to a false-positive detection of organisms' presence. As the eDNA particle size distribution (PSD) changes along with the decay process, it may facilitate solving the above problem. Here, we set up tank experiments with snails, zebrafish and daphnids, respectively, to monitor the change in eDNA PSD and eDNA degradation through time after removing organisms. We found that zebrafish eDNA decays more slowly for larger particle sizes. Across all species tested, the percentage of large size ranges tended to increase over time while the smaller sizes showed relatively fast decay rates. As a result, PSD changed consistently with eDNA decay, although initial PSD varied between species. In combination, we propose that eDNA PSD can be used to assess the current prevalence of organisms at an eDNA sampling location while avoiding false-positives on the presence of species. Our findings expand the applicability of eDNA for monitoring target species in freshwater ecosystems.

Highly comparable metabarcoding results from MGI-Tech and Illumina sequencing platforms

With the developments in DNA nanoball sequencing technologies and the emergence of new platforms, there is an increasing interest in their performance in comparison with the widely used sequencing-by-synthesis methods. Here, we test the consistency of metabarcoding results from DNBSEQ-G400RS (DNA nanoball sequencing platform by MGI-Tech) and NovaSeq 6000 (sequencing-by-synthesis platform by Illumina) platforms using technical replicates of DNA libraries that consist of COI gene amplicons from 120 soil DNA samples. By subjecting raw sequencing data from both platforms to a uniform bioinformatics processing, we found that the proportion of high-quality reads passing through the filtering steps was similar in both datasets. Per-sample operational taxonomic unit (OTU) and amplicon sequence variant (ASV) richness patterns were highly correlated, but sequencing data from DNBSEQ-G400RS harbored a higher number of OTUs. This may be related to the lower dominance of most common OTUs in DNBSEQ data set (thus revealing higher richness by detecting rare taxa) and/or to a lower effective read quality leading to generation of spurious OTUs. However, there was no statistical difference in the ASV and post-clustered ASV richness between platforms, suggesting that additional denoising step in the ASV workflow had effectively removed the 'noisy' reads. Both OTU-based and ASV-based composition were strongly correlated between the sequencing platforms, with essentially interchangeable results. Therefore, we conclude that DNBSEQ-G400RS and NovaSeq 6000 are both equally efficient high-throughput sequencing platforms to be utilized in studies aiming to apply the metabarcoding approach, but the main benefit of the former is related to lower sequencing cost.

Environmental DNA gives comparable results to morphology-based indices of macroinvertebrates in a large-scale ecological assessment

Anthropogenic activities are changing the state of ecosystems worldwide, affecting community composition and often resulting in loss of biodiversity. Rivers are among the most impacted ecosystems. Recording their current state with regular biomonitoring is important to assess the future trajectory of biodiversity. Traditional monitoring methods for ecological assessments are costly and time-intensive. Here, we compared monitoring of macroinvertebrates based on environmental DNA (eDNA) sampling with monitoring based on traditional kick-net sampling to assess biodiversity patterns at 92 river sites covering all major Swiss river catchments. From the kick-net community data, a biotic index (IBCH) based on 145 indicator taxa had been established. The index was matched by the taxonomically annotated eDNA data by using a machine learning approach. Our comparison of diversity patterns only uses the zero-radius Operational Taxonomic Units assigned to the indicator taxa. Overall, we found a strong congruence between both methods for the assessment of the total indicator community composition (gamma diversity). However, when assessing biodiversity at the site level (alpha diversity), the methods were less consistent and gave complementary data on composition. Specifically, environmental DNA retrieved significantly fewer indicator taxa per site than the kick-net approach. Importantly, however, the subsequent ecological classification of rivers based on the detected indicators resulted in similar biotic index scores for the kick-net and the eDNA data that was classified using a random forest approach. The majority of the predictions (72%) from the random forest classification resulted in the same river status categories as the kick-net approach. Thus, environmental DNA validly detected indicator communities and, combined with machine learning, provided reliable classifications of the ecological state of rivers. Overall, while environmental DNA gives complementary data on the macroinvertebrate community composition compared to the kick-net approach, the subsequently calculated indices for the ecological classification of river sites are nevertheless directly comparable and consistent.

Wet grinding of invertebrate bulk samples – a scalable and cost-efficient protocol for metabarcoding and metagenomics

Most metabarcoding protocols for invertebrate bulk samples start with sample homogenisation, followed by DNA extraction, amplification of a specific marker region, and sequencing. Many of the above-mentioned laboratory steps have been verified thoroughly and best practice strategies exist, yet, no clear recommendation for the basis of almost all metabarcoding studies exists: the homogenisation of samples itself. Two different categories of devices are typically used for homogenisation: bead mills or blenders. Both have upsides and downsides. Bead mills rely on single-use plastics and therefore produce a lot of waste and are expensive. In addition to that, processing times can go up to 30 minutes making them unsuitable for large-scale studies. Blenders can handle larger sample volumes in a shorter time, and be cleaned – yet suffer from an increased risk of cross-contamination. We aimed to develop a fast, robust, cheap, and reliable sample homogenisation protocol that overcomes limitations of both approaches, i.e. does not produce difficult to discard waste and avoid single-use plastics while reducing overall costs. We tested the performance of the new protocol using six size-sorted Malaise trap samples and six unsorted stream macroinvertebrate kick-net samples. We used 14 replicates per sample and included many negative controls at different steps of the protocol to quantify the impacts of i) insufficient homogenisation and ii) cross-contamination. Our results show that 3-min homogenisation is sufficient to recover about 80% of OTUs per sample in each replicate and that a non-hazardous DIY cleaning solution provides an effective and efficient way of cleaning. The improvements of the protocol in terms of speed, ease of handling, an overall reduction of costs as well as the documented reliability and robustness make it an important candidate for sample homogenisation after sampling in particular for large-scale and regulatory metabarcoding but also metagenomics biodiversity assessments and monitoring.

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.

Benchmarking bioinformatic tools for fast and accurate eDNA metabarcoding species identification

Bioinformatic analysis of eDNA metabarcoding data is a crucial step toward rigorously assessing biodiversity. Many programs are now available for each step of the required analyses, but their relative abilities at providing fast and accurate species lists have seldom been evaluated. We used simulated mock communities and real fish eDNA metabarcoding data to evaluate the performance of 13 bioinformatic programs and pipelines to retrieve fish occurrence and read abundance using the 12S mt rRNA gene marker. We used four indices to compare the outputs of each program with the simulated samples: sensitivity, F-measure, root-mean-square error (RMSE) on read relative abundances, and execution time. We found marked differences among programs only for the taxonomic assignment step, both in terms of sensitivity, F-measure and RMSE. Running time was highly different between programs for each step. The fastest programs with best indices for each step were assembled into a pipeline. We compared this pipeline to pipelines constructed from existing toolboxes (OBITools, Barque, and QIIME 2). Our pipeline and Barque obtained the best performance for all indices and appear to be better alternatives to highly used pipelines for analysing fish eDNA metabarcoding data when a complete reference database is available. Analysis on real eDNA metabarcoding data also indicated differences for taxonomic assignment and execution time only. This study reveals major differences between programs during the taxonomic assignment step. The choice of algorithm for the taxonomic assignment can have a significant impact on diversity estimates and should be made according to the objectives of the study.

Can we manage fisheries with the inherent uncertainty from eDNA?

Environmental (e)DNA, as a general approach in aquatic systems, seeks to connect the presence of species' genetic material in the water and hence to infer the species' physical presence. However, fisheries managers face making decisions with risk and uncertainty when eDNA indicates a fish is present but traditional methods fail to capture the fish. In comparison with traditional methods such as nets, electrofishing and piscicides, eDNA approaches have more sources of underlying error that could give rise to false positives. This has resulted in some managers to question whether eDNA can be used to make management decisions because there is no fish in hand. As a relatively new approach, the methods and techniques have quickly evolved to improve confidence in eDNA. By evaluating an eDNA based research programmes through the pattern of the eDNA signal, assay design, experimental design, quality assurance and quality control checks, data analyses and concurrent search for fish using traditional gears, the evidence for fish presence can be evaluated to build confidence in the eDNA approach. The benefits for fisheries management from adopting an eDNA approach are numerous but include cost effectiveness, broader geographic coverage of habitat occupancy, early detection of invasive species, non-lethal stock assessments, exploration of previously inaccessible aquatic environments and discovery of new species hidden beneath the water's surface. At a time when global freshwater and marine fisheries are facing growing threats from over-harvest, pollution and climate change, we anticipate that growing confidence in eDNA will overcome the inherent uncertainty of not having a fish in hand and will empower the informed management actions necessary to protect and restore our fisheries.

Translational Molecular Ecology in practice: Linking DNA-based methods to actionable marine environmental management

Molecular-based approaches can provide timely biodiversity assessments, showing an immense potential to facilitate decision-making in marine environmental management. However, the uptake of molecular data into environmental policy remains minimal. Here, we showcase a selection of local to global scale studies applying molecular-based methodologies for environmental management at various stages of implementation. Drawing upon lessons learned from these case-studies, we provide a roadmap to facilitate applications of DNA-based methods to marine policies and to overcome the existing challenges. The main impediment identified is the need for standardized protocols to guarantee data comparison across spatial and temporal scales. Adoption of Translational Molecular Ecology - the sustained collaboration between molecular ecologists and stakeholders, will enhance consensus with regards to the objectives, methods, and outcomes of environmental management projects. Establishing a sustained dialogue among stakeholders is key to accelerating the adoption of molecular-based approaches for marine monitoring and assessment.