Environmental DNA metabarcoding reveals comparable responses to agricultural stressors on different trophic levels of a freshwater community.

Freshwater habitats are under stress from agricultural land use, most notably the influx of neonicotinoid pesticides and increased nutrient pressure from fertilizer. Traditional studies investigating the effects of stressors on freshwater systems are often limited to a narrow range of taxa, depending heavily on morphological expertise. Additionally, disentanglement of multiple simultaneous stressors can be difficult in field studies, whereas controlled laboratory conditions do not accurately reflect natural conditions and food webs. To overcome these drawbacks, we investigated the impacts of two agricultural stressors (the neonicotinoid insecticide thiacloprid and fertilizer) in full-factorial design in a semi-natural research site, using environmental DNA sampling to study three different taxonomic groups representing three trophic levels: bacteria (decomposers), phytoplankton (primary producers), and chironomids (consumers). The results show considerable impact of both stressors across trophic levels, with an additive effect of fertilizer and thiacloprid on community composition at all levels. These findings suggest that agricultural stressors affect the entire food web, either directly or through cascade reactions. They are also consistent with morphological assessments that were performed in the same study site, even at a lower number of replicates. The study presented shows that the use of multimarker environmental DNA provides a more comprehensive assessment of stressor impacts across multiple trophic levels, at a higher taxonomic resolution than traditional surveys. Additionally, many putative novel bioindicators for both agricultural stressors were discovered. We encourage further investigations into stressors impacts at different trophic levels, which will lead to more effective monitoring and management of freshwater systems.

The mitochondrial genome of Nemalecium lighti (Hydrozoa, Leptothecata).

The hydrozoan species Nemalecium lighti (Hargitt, 1924) is widely distributed in tropical marine waters around the world. Here we report the complete linear mitochondrial genome of N. lighti from Sint Eustatius (Lesser Antilles). The mitochondrial genome with a length of 14,320 bp encodes for 13 protein-coding genes, two tRNA genes, and two rRNA genes. Gene arrangement differs from that found in other species of the same taxonomic order and a phylogenetic analysis shows that based on mitochondrial genes, N. lighti clusters outside of the Leptothecata, rendering the order paraphyletic.

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.

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.

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.

Uncovering bacterial and functional diversity in macroinvertebrate mitochondrial-metagenomic datasets by differential centrifugation.

PCR-free techniques such as meta-mitogenomics (MMG) can recover taxonomic composition of macroinvertebrate communities, but suffer from low efficiency, as >90% of sequencing data is mostly uninformative due to the great abundance of nuclear DNA that cannot be identified with current reference databases. Current MMG studies do not routinely check data for information on macroinvertebrate-associated bacteria and gene functions. However, this could greatly increase the efficiency of MMG studies by revealing yet overlooked diversity within ecosystems and making currently unused data available for ecological studies. By analysing six 'mock' communities, each containing three macroinvertebrate taxa, we tested whether this additional data on bacterial taxa and functional potential of communities can be extracted from MMG datasets. Further, we tested whether differential centrifugation, which is known to greatly increase efficiency of macroinvertebrate MMG studies by enriching for mitochondria, impacts on the inferred bacterial community composition. Our results show that macroinvertebrate MMG datasets contain a high number of mostly endosymbiont bacterial taxa and associated gene functions. Centrifugation reduced both the absolute and relative abundance of highly abundant Gammaproteobacteria, thereby facilitating detection of rare taxa and functions. When analysing both taxa and gene functions, the number of features obtained from the MMG dataset increased 31-fold ('enriched') respectively 234-fold ('not enriched'). We conclude that analysing MMG datasets for bacteria and gene functions greatly increases the amount of information available and facilitates the use of shotgun metagenomic techniques for future studies on biodiversity.

DNA barcode reference libraries for the monitoring of aquatic biota in Europe: Gap-analysis and recommendations for future work.

Effective identification of species using short DNA fragments (DNA barcoding and DNA metabarcoding) requires reliable sequence reference libraries of known taxa. Both taxonomically comprehensive coverage and content quality are important for sufficient accuracy. For aquatic ecosystems in Europe, reliable barcode reference libraries are particularly important if molecular identification tools are to be implemented in biomonitoring and reports in the context of the EU Water Framework Directive (WFD) and the Marine Strategy Framework Directive (MSFD). We analysed gaps in the two most important reference databases, Barcode of Life Data Systems (BOLD) and NCBI GenBank, with a focus on the taxa most frequently used in WFD and MSFD. Our analyses show that coverage varies strongly among taxonomic groups, and among geographic regions. In general, groups that were actively targeted in barcode projects (e.g. fish, true bugs, caddisflies and vascular plants) are well represented in the barcode libraries, while others have fewer records (e.g. marine molluscs, ascidians, and freshwater diatoms). We also found that species monitored in several countries often are represented by barcodes in reference libraries, while species monitored in a single country frequently lack sequence records. A large proportion of species (up to 50%) in several taxonomic groups are only represented by private data in BOLD. Our results have implications for the future strategy to fill existing gaps in barcode libraries, especially if DNA metabarcoding is to be used in the monitoring of European aquatic biota under the WFD and MSFD. For example, missing species relevant to monitoring in multiple countries should be prioritized for future collaborative programs. We also discuss why a strategy for quality control and quality assurance of barcode reference libraries is needed and recommend future steps to ensure full utilisation of metabarcoding in aquatic biomonitoring.

The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems.

The bioassessment of aquatic ecosystems is currently based on various biotic indices that use the occurrence and/or abundance of selected taxonomic groups to define ecological status. These conventional indices have some limitations, often related to difficulties in morphological identification of bioindicator taxa. Recent development of DNA barcoding and metabarcoding could potentially alleviate some of these limitations, by using DNA sequences instead of morphology to identify organisms and to characterize a given ecosystem. In this paper, we review the structure of conventional biotic indices, and we present the results of pilot metabarcoding studies using environmental DNA to infer biotic indices. We discuss the main advantages and pitfalls of metabarcoding approaches to assess parameters such as richness, abundance, taxonomic composition and species ecological values, to be used for calculation of biotic indices. We present some future developments to fully exploit the potential of metabarcoding data and improve the accuracy and precision of their analysis. We also propose some recommendations for the future integration of DNA metabarcoding to routine biomonitoring programs.