The influence of macroinvertebrate abundance on the assessment of freshwater quality in The Netherlands

From DNA to diagnostics: A case study using macroinvertebrate metabarcoding to assess the effectiveness of restoration measures in a Dutch stream

Stream ecosystems are under pressure due to multiple stressors. Restoration measures can halt further degradation and improve their ecological status. However, assessment of the effectiveness of the implemented measures is often insufficient because of logistic and financial constraints. DNA-metabarcoding has been proposed to scale up sample processing, although its application as a diagnostic tool has received less attention. The aim of our study was to evaluate if DNA-metabarcoding of stream macroinvertebrates can be used to compute a stressor-specific index to assess the effectiveness of a stream restoration project. For this purpose, we sampled the upstream, restored, and downstream section of a recently restored lowland stream in the Netherlands. At each site, we applied three different methods of macroinvertebrate identification: morphological identification of bulk samples (morphology), DNA-metabarcoding of the same bulk samples (DNA) and metabarcoding of eDNA extracted from the water (eDNA). First, we compared the community composition identified by each method. The communities identified by morphology and DNA were highly similar, whereas the communities generated by the eDNA differed. Second, we analysed whether the identification methods could be used to assess the effectiveness of the restoration project, focussing on a stressor-specific index for flow as the restoration measures aimed at improving flow conditions. Both the morphology and bulk DNA samples indicated improved flow conditions in the restored section of the stream (i.e., less stress from the reduction or absence of flow than in the unrestored sections). Contrary, the eDNA-water samples did not differentiate the amount of stress throughout the catchment, although applying recent developments in eDNA sampling could lead to more robust results. In conclusion, this study forms proof of concept that DNA from bulk samples can be utilized to assess the effectiveness of restoration measures, showing the added value of this approach for water managers.

Testing the Influence of Incomplete DNA Barcode Libraries on Ecological Status Assessment of Mediterranean Transitional Waters

Simple Summary

The biodiversity and ecological status assessment of transitional water ecosystems by benthic macroinvertebrates investigation could use DNA barcode tools for more rapid and efficient outputs. The principal limits of this application are the incompleteness of DNA barcode databases, the identification of optimal primers set, and the gap in the species sequences. The influence of the incompleteness of DNA barcode libraries on species diversity indices, ecological indicators, and ecological status assessment in transitional waters of the southeast Mediterranean were analysed, underlying the importance to implement DNA barcode libraries and to put an effort toward specific species at a local level.

Abstract

The ecological assessment of European aquatic ecosystems is regulated under the framework directives on strategy for water and marine environments. Benthic macroinvertebrates are the most used biological quality element for ecological assessment of rivers, coastal-marines, and transitional waters. The morphological identification of benthic macroinvertebrates is the current tool for their assessment. Recently, DNA-based tools have been proposed as effective alternatives. The main current limits of DNA-based applications include the incompleteness of species recorded in the DNA barcode reference libraries and the primers bias. Here, we analysed the influence of the incompleteness of DNA barcode databases on species diversity indices, ecological indicators, and ecological assessment in transitional waters of the southeast Mediterranean, taking into account the availability of commonly sequenced and deposited genomic regions for listed species. The ecological quality status assigned through the potential application of both approaches to the analysed transitional water ecosystems was different in 27% of sites. We also analysed the inter-specific genetic distances to evaluate the potential application of the DNA metabarcoding method. Overall, this work highlights the importance to expand the barcode databases and to analyse, at the regional level, the gaps in the DNA barcodes.

All together now: Limitations and recommendations for the simultaneous analysis of all eukaryotic soil sequences

The soil environment contains a large, but historically underexplored, reservoir of biodiversity. Sequencing prokaryotic marker genes has become commonplace for the discovery and characterization of soil bacteria and archaea. Increasingly, this approach is also applied to eukaryotic marker genes to characterize the diversity and distribution of soil eukaryotes. However, understanding the properties and limitations of eukaryotic marker sequences is essential for correctly analysing, interpreting, and synthesizing the resulting data. Here, we illustrate several biases from sequencing data that affect measurements of biodiversity that arise from variation in morphology, taxonomy and phylogeny between organisms, as well as from sampling designs. We recommend analytical approaches to overcome these limitations, and outline how the benchmarking and standardization of sequencing protocols may improve the comparability of the data.

Exploring the capacity of aquatic biofilms to act as environmental DNA samplers: Test on macroinvertebrate communities in rivers

Aquatic biofilms are heterogeneous assemblages of microorganisms surrounded by a matrix of extracellular polymeric substances (EPS). Recent studies suggest that aquatic biofilms can physically act as sorptive sponges of DNA. We took the opportunity from already available samples of stone biofilms and macroinvertebrates specimens collected in parallel at the same sites to test the capacity of biofilms to act as DNA samplers of macroinvertebrate communities in streams. Macroinvertebrate communities are usually studied with metabarcoding using the DNA extracted from their bodies bulk samples, which remains a time-consuming approach and involves the destruction of all individual specimens from the samples. The ability of biofilms to capture DNA was explored on 19 rivers sites of a tropical island (Mayotte Island, France). First, macroinvertebrate specimens were identified based on their morphological characteristics. Second, DNA was extracted from biofilms, and macroinvertebrate communities were targeted using a standard COI barcode. The resulting morphological and molecular inventories were compared. They provided comparable structures and diversities for macroinvertebrate communities when one is working with the unassigned OTU data. After taxonomic assignment of the OTU data, diversity and richness were no longer correlated. The ecological assessment derived from morphological bulk samples was conserved by the biofilms samples. We also showed that the biofilm method allows to detect a higher diversity for some organisms (Cnidaria), that is hardly accessible with the morphological method. The results of this study exploring the DNA signal captured by natural biofilms are encouraging. However, a more detailed study integrating more replicates and comparing the biodiversity signal based on both morphological and molecular bulk macroinvertebrate samples to the one captured by biofilms will be necessary. Better understanding how the DNA signal captured by natural biofilms represents the biodiversity of a given sampling site is necessary before considering its use for bioassessment applications.

TaxonTableTools: A comprehensive, platform-independent graphical user interface software to explore and visualise DNA metabarcoding data

DNA metabarcoding is increasingly used as a tool to assess biodiversity in research and environmental management. Powerful analysis software exists to process raw data. However, the translation of sequence read data into biological information and downstream analyses may be difficult for end users with limited expertise in bioinformatics. Thus, there is a growing need for easy-to-use, graphical user interface (GUI) software to analyse and visualise DNA metabarcoding data. Here, we present TaxonTableTools (TTT), a new platform-independent GUI that aims to fill this gap by providing simple, reproducible analysis and visualisation workflows. At its base, TTT uses a "TaXon table", which is a data format that can be generated easily within TTT from two input files: a read table and a taxonomy table obtained using various published metabarcoding pipelines. TTT analysis and visualisation modules include Venn diagrams to compare taxon overlap among replicates, samples, or analysis methods. TTT analyses and visualises basic statistics, such as read proportion per taxon, as well as more sophisticated visualisations, such as interactive Krona charts for taxonomic data exploration. Various ecological analyses can be produced directly, including alpha or beta diversity estimates, and rarefaction analysis ordination plots. Metabarcoding data can be converted into formats required for traditional, taxonomy-based analyses performed by regulatory bioassessment programs. In addition, TTT is able to produce html-based interactive graphics that can be analysed in any web browser. The software comes with a manual and tutorial, is free and publicly available through GitHub (https://github.com/TillMacher/TaxonTableTools) or the Python package index (https://pypi.org/project/taxontabletools/).

Connecting high-throughput biodiversity inventories: Opportunities for a site-based genomic framework for global integration and synthesis

High-throughput sequencing (HTS) is increasingly being used for the characterization and monitoring of biodiversity. If applied in a structured way, across broad geographical scales, it offers the potential for a much deeper understanding of global biodiversity through the integration of massive quantities of molecular inventory data generated independently at local, regional and global scales. The universality, reliability and efficiency of HTS data can potentially facilitate the seamless linking of data among species assemblages from different sites, at different hierarchical levels of diversity, for any taxonomic group and regardless of prior taxonomic knowledge. However, collective international efforts are required to optimally exploit the potential of site-based HTS data for global integration and synthesis, efforts that at present are limited to the microbial domain. To contribute to the development of an analogous strategy for the nonmicrobial terrestrial domain, an international symposium entitled "Next Generation Biodiversity Monitoring" was held in November 2019 in Nicosia (Cyprus). The symposium brought together evolutionary geneticists, ecologists and biodiversity scientists involved in diverse regional and global initiatives using HTS as a core tool for biodiversity assessment. In this review, we summarize the consensus that emerged from the 3-day symposium. We converged on the opinion that an effective terrestrial Genomic Observatories network for global biodiversity integration and synthesis should be spatially led and strategically united under the umbrella of the metabarcoding approach. Subsequently, we outline an HTS-based strategy to collectively build an integrative framework for site-based biodiversity data generation.

Functional traits explain crayfish invasive success in the Netherlands

Biological invasions by nonindigenous species can have negative effects on economies and ecosystems. To limit this impact, current research on biological invasions uses functional traits to facilitate a mechanistic understanding of theoretical and applied questions. Here we aimed to assess the role of functional traits in the progression of crayfish species through different stages of invasion and determine the traits associated with invasive success. A dataset of thirteen functional traits of 15 species currently occurring or available for sale in the Netherlands was evaluated. Six of these crayfish appeared invasive. Important traits distinguishing successful from unsuccessful invaders were a temperate climate in the native range, a medium to high egg count and producing more than one egg clutch per year. The most successful invaders had different functional trait combinations: Procambarus clarkii has a higher reproductive output, can migrate over longer distances and possesses a higher aggression level; Faxonius limosus is adapted to a colder climate, can reproduce parthenogetically and has broader environmental tolerances. Using a suit of functional traits to analyse invasive potential can help risk management and prevention. For example, based on our data Procambarus virginalis is predicted to become the next successful invasive crayfish in the Netherlands.

Decision-making and best practices for taxonomy-free environmental DNA metabarcoding in biomonitoring using Hill numbers

Environmental DNA (eDNA) metabarcoding is raising expectations for biomonitoring of organisms that have hitherto been neglected. To bypass current limitations in taxonomic assignments due to incomplete or erroneous reference databases, taxonomy-free approaches are proposed for biomonitoring at the level of operational taxonomic units (OTUs). This is challenging, because OTUs cannot be annotated and directly compared against classically derived taxonomic data. The application of good stringency treatments to infer the validity of OTUs and clear understanding of the consequences of such treatments is especially relevant for biodiversity assessments. We investigated how common practices of stringency filtering affect eDNA diversity estimates in the statistical framework of Hill numbers. We collected water eDNA samples at 61 sites across a 740-km² river catchment, reflecting a spatially realistic scenario in biomonitoring. After bioinformatic processing of the data, we studied how different stringency treatments affect conclusions with respect to biodiversity at the catchment and site levels. The applied stringency treatments were based on the consistent appearance of OTUs across filter replicates, a relative abundance cut-off and rarefaction. We detected large differences in diversity estimates when accounting for presence/absence only, such that detected diversity at the catchment scale differed by an order of magnitude between the treatments. These differences disappeared when using stringency treatments with increasing weighting of the OTU abundances. Our study demonstrated the usefulness of Hill numbers for biodiversity analyses and comparisons of eDNA data sets that strongly differ in diversity. We recommend best practice for data stringency filtering for biomonitoring using eDNA.

Can metabarcoding resolve intraspecific genetic diversity changes to environmental stressors? A test case using river macrozoobenthos

Genetic diversity is the most basal level of biodiversity and determines the evolutionary capacity of species to adapt to changing environments, yet it is typically neglected in routine biomonitoring and stressor impact assessment. For a comprehensive analysis of stressor impacts on genetic diversity, it is necessary to assess genetic variants simultaneously in many individuals and species. Such an assessment is not as straightforward and usually limited to one or few focal species. However, nowadays species diversity can be assessed by analysing thousands of individuals of a community simultaneously with DNA metabarcoding. Recent bioinformatic advances also allow for the extraction of exact sequence variants (ESVs or haplotypes) in addition to Operational Taxonomic Units (OTUs). By using this new capability, we here evaluated if the analysis of intraspecific mitochondrial diversity in addition to species diversity can provide insights into responses of stream macrozoobenthic communities to environmental stressors. For this purpose, we analysed macroinvertebrate bulk samples of three German river systems with different stressor levels using DNA metabarcoding. While OTU and haplotype number were negatively correlated with stressor impact, this association was not as clear when studying haplotype diversity across all taxa. However, stressor responses were found for sensitive EPT (Ephemeroptera, Plecoptera, Trichoptera) taxa and those exceedingly resistant to organic stress. An increase in haplotype number per OTU and haplotype diversity of sensitive taxa was observed with an increase in ecosystem quality and stability, while the opposite pattern was detected for pollution resistant taxa. However, this pattern was less prominent than expected based on the strong differences in stressor intensity between sites. To compare genetic diversity among communities in river systems, we focussed on OTUs, which were present in all systems. As OTU composition differed strongly between rivers, this led to the exclusion of a high number of OTUs, especially in diverse river systems of good quality, which potentially diminished the increase in intraspecific diversity. To better understand responses of intraspecific genetic diversity to environmental stressors, for example in river ecosystems, it would be important to increase OTU overlap between compared sites, e.g. by sampling a narrower stressor gradient, and to perform calibrated studies controlling for the number of individuals and their haplotypes. However, this pioneer study shows that the extraction of haplotypes from DNA metabarcoding datasets is a promising source of information to simultaneously assess intraspecific diversity changes in response to environmental impacts for a metacommunity.

Analysis of 13,312 benthic invertebrate samples from German streams reveals minor deviations in ecological status class between abundance and presence/absence data

Benthic invertebrates are the most commonly used organisms used to assess ecological status as required by the EU Water Framework Directive (WFD). For WFD-compliant assessments, benthic invertebrate communities are sampled, identified and counted. Taxa × abundance matrices are used to calculate indices and the resulting scores are compared to reference values to determine the ecological status class. DNA-based tools, such as DNA metabarcoding, provide a new and precise method for species identification but cannot deliver robust abundance data. To evaluate the applicability of DNA-based tools to ecological status assessment, we evaluated whether the results derived from presence/absence data are comparable to those derived from abundance data. We analysed benthic invertebrate community data obtained from 13,312 WFD assessments of German streams. Broken down to 30 official stream types, we compared assessment results based on abundance and presence/absence data for the assessment modules “organic pollution” (i.e., the saprobic index) and “general degradation” (a multimetric index) as well as their underlying metrics.

In 76.6% of cases, the ecological status class did not change after transforming abundance data to presence/absence data. In 12% of cases, the status class was reduced by one (e.g., from good to moderate), and in 11.2% of cases, the class increased by one. In only 0.2% of cases, the status shifted by two classes. Systematic stream type-specific deviations were found and differences between abundance and presence/absence data were most prominent for stream types where abundance information contributed directly to one or several metrics of the general degradation module. For a single stream type, these deviations led to a systematic shift in status from ‘good’ to ‘moderate’ (n = 201; with only n = 3 increasing). The systematic decrease in scores was observed, even when considering simulated confidence intervals for abundance data. Our analysis suggests that presence/absence data can yield similar assessment results to those for abundance-based data, despite type-specific deviations. For most metrics, it should be possible to intercalibrate the two data types without substantial efforts. Thus, benthic invertebrate taxon lists generated by standardised DNA-based methods should be further considered as a complementary approach.

Increased performance of DNA metabarcoding of macroinvertebrates by taxonomic sorting

DNA-based identification through the use of metabarcoding has been proposed as the next step in the monitoring of biological communities, such as those assessed under the Water Framework Directive (WFD). Advances have been made in the field of metabarcoding, but challenges remain when using complex samples. Uneven biomass distributions, preferential amplification and reference database deficiencies can all lead to discrepancies between morphological and DNA-based taxa lists. The effects of different taxonomic groups on these issues remain understudied. By metabarcoding WFD monitoring samples, we analyzed six different taxonomic groups of freshwater organisms, both separately and combined. Identifications based on metabarcoding data were compared directly to morphological assessments performed under the WFD. The diversity of taxa for both morphological and DNA-based assessments was similar, although large differences were observed in some samples. The overlap between the two taxon lists was 56.8% on average across all taxa, and was highest for Crustacea, Heteroptera, and Coleoptera, and lowest for Annelida and Mollusca. Taxonomic sorting in six basic groups before DNA extraction and amplification improved taxon recovery by 46.5%. The impact on ecological quality ratio (EQR) scoring was considerable when replacing morphology with DNA-based identifications, but there was a high correlation when only replacing a single taxonomic group with molecular data. Different taxonomic groups provide their own challenges and benefits. Some groups might benefit from a more consistent and robust method of identification. Others present difficulties in molecular processing, due to uneven biomass distributions, large genetic diversity or shortcomings of the reference database. Sorting samples into basic taxonomic groups that require little taxonomic knowledge greatly improves the recovery of taxa with metabarcoding. Current standards for EQR monitoring may not be easily replaced completely with molecular strategies, but the effectiveness of molecular methods opens up the way for a paradigm shift in biomonitoring.

Advancing the use of molecular methods for routine freshwater macroinvertebrate biomonitoring – the need for calibration experiments

Over the last decade, steady advancements have been made in the use of DNA-based methods for detection of species in a wide range of ecosystems. This progress has culminated in molecular monitoring methods being employed for the detection of several species for enforceable management purposes of endangered, invasive, and illegally harvested species worldwide. However, the routine application of DNA-based methods to monitor whole communities (typically a metabarcoding approach) in order to assess the status of ecosystems continues to be limited. In aquatic ecosystems, the limited use is particularly true for macroinvertebrate communities. As part of the DNAqua-Net consortium, a structured discussion was initiated with the aim to identify potential molecular methods for freshwater macroinvertebrate community assessment and identify important knowledge gaps for their routine application. We focus on three complementary DNA sources that can be metabarcoded: 1) DNA from homogenised samples (bulk DNA), 2) DNA extracted from sample preservative (fixative DNA), and 3) environmental DNA (eDNA) from water or sediment. We provide a brief overview of metabarcoding macroinvertebrate communities from each DNA source and identify challenges for their application to routine monitoring. To advance the utilisation of DNA-based monitoring for macroinvertebrates, we propose an experimental design template for a series of methodological calibration tests. The template compares sources of DNA with the goal of identifying the effects of molecular processing steps on precision and accuracy. Furthermore, the same samples will be morphologically analysed, which will enable the benchmarking of molecular to traditional processing approaches. In doing so we hope to highlight pathways for the development of DNA-based methods for the monitoring of freshwater macroinvertebrates.