Monitoring the ecological status of rivers with diatom eDNA metabarcoding: A comparison of taxonomic markers and analytical approaches for the inference of a molecular diatom index.

Recently, several studies demonstrated the usefulness of diatom eDNA metabarcoding as an alternative to assess the ecological quality of rivers and streams. However, the choice of the taxonomic marker as well as the methodology for data analysis differ between these studies, hampering the comparison of their results and effectiveness. The aim of this study was to compare two taxonomic markers commonly used in diatom metabarcoding and three distinct analytical approaches to infer a molecular diatom index. We used the values of classical morphological diatom index as a benchmark for this comparison. We amplified and sequenced both a fragment of the rbcL gene and the V4 region of the 18S rRNA gene for 112 epilithic samples from Swiss and French rivers. We inferred index values using three analytical approaches: by computing it directly from taxonomically assigned sequences, by calibrating de novo the ecovalues of all metabarcodes, and by using a supervised machine learning algorithm to train predictive models. In general, the values of index obtained using the two "taxonomy-free" approaches, encompassing molecular assignment and machine learning, were closer correlated to the values of the morphological index than the values based on taxonomically assigned sequences. The correlations of the three analytical approaches were higher in the case of rbcL compared to the 18S marker, highlighting the importance of the reference database which is more complete for the rbcL marker. Our study confirms the effectiveness of diatom metabarcoding as an operational tool for rivers ecological quality assessment and shows that the analytical approaches by-passing the taxonomic assignments are particularly efficient when reference databases are incomplete.

Benthic monitoring of oil and gas offshore platforms in the North Sea using environmental DNA metabarcoding.

Since 2010, considerable efforts have been undertaken to monitor the environmental status of European marine waters and ensuring the development of methodological standards for the evaluation of this status. However, the current routine biomonitoring implicates time-consuming and costly manual sorting and morphological identification of benthic macrofauna. Environmental DNA (eDNA) metabarcoding represents an alternative to the traditional monitoring method with very promising results. Here, we tested it further by performing eDNA metabarcoding of benthic eukaryotic communities in the vicinity of two offshore oil and gas platforms in the North Sea. Three different genetic markers (18S V1V2, 18S V9 and COI) were used to assess the environmental pressures induced by the platforms. All markers showed patterns of alpha and beta diversity consistent with morphology-based macrofauna analyses. In particular, the communities' structure inferred from metabarcoding and morphological data significantly changed along distance gradients from the platforms. The impact of the operational discharges was also detected by the variation of biotic index values, AMBI index showing the best correlation between morphological and eDNA data sets. Finally, the sediment physicochemical parameters were used to build a local de novo pressure index that served as benchmark to test the potential of a taxonomy-free approach. Our study demonstrates that metabarcoding approach outperforms morphology-based approach and can be used as a cost and time-saving alternative solution to the traditional morphology-based monitoring in order to monitor more efficiently the impact of industrial activities on marine biodiversity.

Ecosystems monitoring powered by environmental genomics: A review of current strategies with an implementation roadmap.

A decade after environmental scientists integrated high-throughput sequencing technologies in their toolbox, the genomics-based monitoring of anthropogenic impacts on the biodiversity and functioning of ecosystems is yet to be implemented by regulatory frameworks. Despite the broadly acknowledged potential of environmental genomics to this end, technical limitations and conceptual issues still stand in the way of its broad application by end-users. In addition, the multiplicity of potential implementation strategies may contribute to a perception that the routine application of this methodology is premature or "in development", hence restraining regulators from binding these tools into legal frameworks. Here, we review recent implementations of environmental genomics-based methods, applied to the biomonitoring of ecosystems. By taking a general overview, without narrowing our perspective to particular habitats or groups of organisms, this paper aims to compare, review and discuss the strengths and limitations of four general implementation strategies of environmental genomics for monitoring: (a) Taxonomy-based analyses focused on identification of known bioindicators or described taxa; (b) De novo bioindicator analyses; (c) Structural community metrics including inferred ecological networks; and (d) Functional community metrics (metagenomics or metatranscriptomics). We emphasise the utility of the three latter strategies to integrate meiofauna and microorganisms that are not traditionally utilised in biomonitoring because of difficult taxonomic identification. Finally, we propose a roadmap for the implementation of environmental genomics into routine monitoring programmes that leverage recent analytical advancements, while pointing out current limitations and future research needs.

Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring.

The last decade brought a spectacular development of so-called environmental (e)DNA studies. In general, "environmental DNA" is defined as DNA isolated from environmental samples, in contrast to genomic DNA that is extracted directly from specimens. However, the variety of different sources of eDNA and the range of taxonomic groups that are targeted by eDNA studies is large, which has led to some discussion about the breadth of the eDNA concept. In particular, there is a recent trend to restrict the use of the term "eDNA" to the DNA of macro-organisms, which are not physically present in environmental samples. In this paper, we argue that such a distinction may not be ideal, because the eDNA signal can come from organisms across the whole tree of life. Consequently, we advocate that the term "eDNA" should be used in its generic sense, as originally defined, encompassing the DNA of all organisms present in environmental samples, including microbial, meiofaunal and macrobial taxa. We first suggest specifying the environmental origin of the DNA sample, such as water eDNA, sediment eDNA or soil eDNA. A second specification would then define the taxonomic group targeted through polymerase chain reaction amplification, such as fish eDNA, invertebrate eDNA and bacterial eDNA. This terminology does also not require assumptions about the specific state of the DNA sampled (intracellular or extracellular). We hope that such terminology will help better define the scope of eDNA studies, especially for environmental managers, who use them as reference in routine biomonitoring and bioassessment.

Molecular Phylogeny and Morphology of Leannia veloxifera n. gen. et sp. Unveils a New Lineage of Monothalamous Foraminifera.

Monothalamous (single-chambered) foraminifera have long been considered as the "poor cousins" of multichambered species, which calcareous and agglutinated tests dominate in the fossil record. This view is currently changing with environmental DNA surveys showing that the monothalamids may be as diverse as hard-shelled foraminifera. Yet, the majority of numerous molecular lineages revealed by eDNA studies remain anonymous. Here, we describe a new monothalamous species and genus isolated from the sample of sea grass collected in Gulf of Eilat (Red Sea). This new species, named Leannia veloxifera, is characterized by a tiny ovoid theca (about 50-100 µm) composed of thin organic wall, with two opposite apertures. The examined individuals are multinucleated and show very active reticulopodial movement. Phylogenetic analyses of SSU rDNA, actin, and beta-tubulin (ß-tubulin) show that the species represents a novel lineage branching separately from other monothalamous foraminifera. Interestingly, the SSU rDNA sequence of the new species is very similar to an environmental foraminiferal sequence from Bahamas, suggesting that the novel lineage may represent a group of shallow-water tropical allogromiids, poorly studied until now.

Environmental Monitoring: Inferring the Diatom Index from Next-Generation Sequencing Data.

Diatoms are widely used as bioindicators for the assessment of water quality in rivers and streams. Classically, the diatom biotic indices are based on the relative abundance of morphologically identified species weighted by their autoecological value. Obtaining such indices is time-consuming, costly, and requires excellent taxonomic expertise, which is not always available. Here we tested the possibility to overcome these limitations using a next-generation sequencing (NGS) approach to identify and quantify diatoms found in environmental DNA and RNA samples. We analyzed 27 river sites in the Geneva area (Switzerland), in order to compare the values of the Swiss Diatom Index (DI-CH) computed either by microscopic quantification of diatom species or directly from NGS data. Despite gaps in the reference database and variations in relative abundance of analyzed species, the diatom index shows a significant correlation between morphological and molecular data indicating similar biological quality status for the majority of sites. This proof-of-concept study demonstrates the potential of the NGS approach for identification and quantification of diatoms in environmental samples, opening new avenues toward the routine application of genetic tools for bioassessment and biomonitoring of aquatic ecosystems.

Monitoring newt communities in urban area using eDNA metabarcoding.

Newts are amphibians commonly present in small ponds or garden pools in urban areas. They are protected in many countries and their presence is monitored through visual observation and/or trapping. However, newts are not easy to spot as they are small, elusive and often hidden at the bottom of water bodies. In recent years, environmental DNA (eDNA) has become a popular tool for detecting newts, with a focus on individual species using qPCR assays. Here, we assess the effectiveness of eDNA metabarcoding compared to conventional visual surveys of newt diversity in 45 ponds within urban areas of Geneva canton, Switzerland. We designed newt-specific mitochondrial 16S rRNA primers, which assign the majority of amplicons to newts, and were able to detect four species known to be present in the region, including the invasive subspecies Lissotriton vulgaris meridionalis, native to the Italian peninsula, that has been introduced in the Geneva area recently. The obtained eDNA results were congruent overall with conventional surveys, confirming the morphological observations in the majority of cases (67%). In 25% of cases, a species was only detected genetically, while in 8% of cases, the observations were not supported by eDNA metabarcoding. Our study confirms the usefulness of eDNA metabarcoding as a tool for the effective and non-invasive monitoring of newt community and suggests its broader use for the survey of newt diversity in urban area at larger scales.

Broad sampling of monothalamids (Rhizaria, Foraminifera) gives further insight into diversity of non-marine Foraminifera.

Non-marine foraminifera are among the least known groups of protists and only a handful of species have been described since the 19th century. We collected one naked and five morphologically almost identical organic-walled monothalamid species from freshwater and terrestrial environments from Germany and Austria. One of the species was identified as Lieberkuehnia wageneriClaparède and Lachmann, 1859. As its original description is ambiguous and its type specimen has been lost, a neotype is proposed. We describe four new organic-walled monothalamous foraminifera and a novel Reticulomyxa species both morphologically and genetically. Analyses of molecular data of the different isolates revealed that they are distributed across six different clades. Two new genera, Claparedellus gen. nov. and Velamentofex gen. nov., and five new monothalamous families, Lacogromiidae fam. nov., Limnogromiidae fam. nov., Lieberkuehniidae fam. nov., Edaphoallogromiidae fam. nov. and Velamentofexidae fam. nov., are established.

The Foreign Oligochaete Species Quistadrilus multisetosus (Smith, 1900) in Lake Geneva: Morphological and Molecular Characterization and Environmental Influences on Its Distribution.

The presence of the oligochaete species Quistadrilus multisetosus (Smith, 1900) originating from North America has been mentioned for several decades in Europe, the Middle East and Russia. Its distribution and abundance in Europe is still unknown but it can be considered as potentially invasive. This species was recently discovered in Lake Geneva (Switzerland/France) and three other Swiss lakes. The aims of the present work are to report its repartition and abundance in Lake Geneva, to study its ecology and to determine its invasive potential in this lake. We also provide an identification key for correctly differentiating Q. multisetosus from the closely related species Spirosperma ferox Eisen, 1879 and Embolocephalus velutinus (Grube, 1879), and study the phylogenetic position of Q. multisetosus within several Tubificinae lineages based on the cytochrome c oxidase (COI) marker. Twenty-eight sites have been monitored since 2009 in Lake Geneva. In several sites, the COI sequence corresponding to this species was also searched for in sediment samples using high-throughput sequencing. In addition, we examined specimens collected in this lake before 2009 likely to belong to Q. multisetosus and to have been misidentified. We found that Q. multisetosus was only present in the lake downstream of a wastewater treatment plant and a combined sewer overflow in the Vidy Bay (near Lausanne) and at a site located nearby. These results confirmed the high tolerance of this species to organic matter pollution. Q. multisetosus was already present in this location in 1974 (misidentified as Spirosperma ferox), which suggests that Q. multisetosus has a limited capacity to disseminate in this lake. However, we recommend continuing monitoring its presence in Lake Geneva in the future, especially in the context of warming of waters that could contribute to the expansion of this species.

High-throughput DNA barcoding of oligochaetes for abundance-based indices to assess the biological quality of sediments in streams and lakes.

Aquatic oligochaete communities are valuable indicators of the biological quality of sediments in streams and lakes, but identification of specimens to the species level based on morphological features requires solid expertise in taxonomy and is possible only for a fraction of specimens present in a sample. The identification of aquatic oligochaetes using DNA barcodes would facilitate their use in biomonitoring and allow a wider use of this taxonomic group for ecological diagnoses. Previous approaches based on DNA metabarcoding of samples composed of total sediments or pools of specimens have been proposed for assessing the biological quality of ecosystems, but such methods do not provide precise information on species abundance, which limits the value of resulting ecological diagnoses. Here, we tested how a DNA barcoding approach based on high-throughput sequencing of sorted and genetically tagged specimens performed to assess oligochaete species diversity and abundance and the biological quality of sediments in streams and lakes. We applied both molecular and morphological approaches at 13 sites in Swiss streams and at 7 sites in Lake Geneva. We genetically identified 33 or 66 specimens per site. For both approaches, we used the same index calculations. We found that the ecological diagnoses derived from the genetic approach matched well with those of the morphological approach and that the genetic identification of only 33 specimens per site provided enough ecological information for correctly estimating the biological quality of sediments in streams and lakes.

Diatom DNA metabarcoding for ecological assessment: Comparison among bioinformatics pipelines used in six European countries reveals the need for standardization.

Ecological assessment of lakes and rivers using benthic diatom assemblages currently requires considerable taxonomic expertise to identify species using light microscopy. This traditional approach is also time-consuming. Diatom metabarcoding is a promising alternative and there is increasing interest in using this approach for routine assessment. However, until now, analysis protocols for diatom metabarcoding have been developed and optimised by research groups working in isolation. The diversity of existing bioinformatics methods highlights the need for an assessment of the performance and comparability of results of different methods. The aim of this study was to test the correspondence of outputs from six bioinformatics pipelines currently in use for diatom metabarcoding in different European countries. Raw sequence data from 29 biofilm samples were treated by each of the bioinformatics pipelines, five of them using the same curated reference database. The outputs of the pipelines were compared in terms of sequence unit assemblages, taxonomic assignment, biotic index score and ecological assessment outcomes. The three last components were also compared to outputs from traditional light microscopy, which is currently accepted for ecological assessment of phytobenthos, as required by the Water Framework Directive. We also tested the performance of the pipelines on the two DNA markers (rbcL and 18S-V4) that are currently used by the working groups participating in this study. The sequence unit assemblages produced by different pipelines showed significant differences in terms of assigned and unassigned read numbers and sequence unit numbers. When comparing the taxonomic assignments at genus and species level, correspondence of the taxonomic assemblages between pipelines was weak. Most discrepancies were linked to differential detection or quantification of taxa, despite the use of the same reference database. Subsequent calculation of biotic index scores also showed significant differences between approaches, which were reflected in the final ecological assessment. Use of the rbcL marker always resulted in better correlation among molecular datasets and also in results closer to these generated using traditional microscopy. This study shows that decisions made in pipeline design have implications for the dataset's structure and the taxonomic assemblage, which in turn may affect biotic index calculation and ecological assessment. There is a need to define best-practice bioinformatics parameters in order to ensure the best representation of diatom assemblages. Only the use of similar parameters will ensure the compatibility of data from different working groups. The future of diatom metabarcoding for ecological assessment may also lie in the development of new metrics using, for example, presence/absence instead of relative abundance data.

Benthic foraminiferal metabarcoding and morphology-based assessment around three offshore gas platforms: Congruence and complementarity.

Since the 1960 s, there has been a rapid expansion of drilling activities in the central and northern Adriatic Sea to meet the increasing global energy demand. The discharges of organic and inorganic pollutants, as well as the alteration of the sediment substrate, are among the main impacts associated with these activities. In the present study, we evaluate the response of benthic foraminifera to the activities of three gas platforms in the northwestern Adriatic Sea, with a special focus on the Armida A platform for which extensive geochemical data (organic matter, trace elements, polycyclic aromatic hydrocarbons, other hydrocarbons, and volatile organic compounds) are available. The response to disturbance is assessed by analyzing the foraminiferal diversity using the traditional morphology-based approach and by 18S rDNA-based metabarcoding. The two methods give congruent results, showing relatively lower foraminiferal diversity and higher dominance values at stations closer to the platforms (<50 m). The taxonomic compositions of the morphological and metabarcoding datasets are very different, the latter being dominated by monothalamous, mainly soft-walled species. However, compositional changes consistently occur at 50 m from the platform and can be related to variations in sediment grain-size variation and higher concentrations of Ni, Zn, Ba, hydrocarbons and total organic carbon. Additionally, several morphospecies and Molecular Operational Taxonomic Units (MOTUs) show strong correlations with distance from the platform and with environmental parameters extracted from BIOENV analysis. Some of these MOTUs have the potential to become new bioindicators, complementing the assemblage of hard-shelled foraminiferal species detected through microscopic analyses. The congruence and complementarity between metabarcoding and morphological approaches support the application of foraminiferal metabarcoding in routine biomonitoring surveys as a reliable, time- and cost-effective methodology to assess the environmental impacts of marine industries.

Embracing Environmental Genomics and Machine Learning for Routine Biomonitoring.

Genomics is fast becoming a routine tool in medical diagnostics and cutting-edge biotechnologies. Yet, its use for environmental biomonitoring is still considered a futuristic ideal. Until now, environmental genomics was mainly used as a replacement of the burdensome morphological identification, to screen known morphologically distinguishable bioindicator taxa. While prokaryotic and eukaryotic microbial diversity is of key importance in ecosystem functioning, its implementation in biomonitoring programs is still largely unappreciated, mainly because of difficulties in identifying microbes and limited knowledge of their ecological functions. Here, we argue that the combination of massive environmental genomics microbial data with machine learning algorithms can be extremely powerful for biomonitoring programs and pave the way to fill important gaps in our understanding of microbial ecology.

Multi-marker eDNA metabarcoding survey to assess the environmental impact of three offshore gas platforms in the North Adriatic Sea (Italy).

The environmental DNA (eDNA) metabarcoding represents a new promising tool for biomonitoring and environmental impact assessment. One of the main advantages of eDNA metabarcoding, compared to the traditional morphotaxonomy-based methods, is to provide a more holistic biodiversity information that includes inconspicuous morphologically non-identifiable taxa. Here, we use eDNA metabarcoding to survey marine biodiversity in the vicinity of the three offshore gas platforms in North Adriatic Sea (Italy). We isolated eDNA from 576 water and sediment samples collected at 32 sampling sites situated along four axes at increasing distances from the gas platforms. We obtained about 46 million eDNA sequences for 5 markers from nuclear 18S V1V2, 18S V4, 18S 37F and mitochondrial 16S and COI genes that cover a wide diversity of benthic and planktonic eukaryotes. Our results showed some impact of platform activities on benthic and pelagic communities at very close distance (<50 m), while communities for intermediate (125 m, 250 m, 500 m) and reference (1000 m, 2000 m) sites did not show any particular biodiversity changes that could be related to platforms activities. The most significant community change along the distance gradient was obtained with the 18S V1V2 marker targeting benthic eukaryotes, even though other markers showed similar trends, but to a lesser extent. These results were congruent with the AMBI index inferred from the eDNA sequences assigned to benthic macrofauna. We finally explored the relation between various physicochemical parameters, including hydrocarbons, on benthic community in the case of one of the platforms. Our results showed that these communities were not significantly impacted by most of hydrocarbons, but rather by macro-elements and sediment texture.

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.

Development and implementation of eco-genomic tools for aquatic ecosystem biomonitoring: the SYNAQUA French-Swiss program.

The effectiveness of environmental protection measures is based on the early identification and diagnosis of anthropogenic pressures. Similarly, restoration actions require precise monitoring of changes in the ecological quality of ecosystems, in order to highlight their effectiveness. Monitoring the ecological quality relies on bioindicators, which are organisms revealing the pressures exerted on the environment through the composition of their communities. Their implementation, based on the morphological identification of species, is expensive because it requires time and experts in taxonomy. Recent genomic tools should provide access to reliable and high-throughput environmental monitoring by directly inferring the composition of bioindicators' communities from their DNA (metabarcoding). The French-Swiss program SYNAQUA (INTERREG France-Switzerland 2017-2019) proposes to use and validate the tools of environmental genomic for biomonitoring and aims ultimately at their implementation in the regulatory bio-surveillance. SYNAQUA will test the metabarcoding approach focusing on two bioindicators, diatoms, and aquatic oligochaetes, which are used in freshwater biomonitoring in France and Switzerland. To go towards the renewal of current biomonitoring practices, SYNAQUA will (1) bring together different actors: scientists, environmental managers, consulting firms, and biotechnological companies, (2) apply this approach on a large scale to demonstrate its relevance, (3) propose robust and reliable tools, and (4) raise public awareness and train the various actors likely to use these new tools. Biomonitoring approaches based on such environmental genomic tools should address the European need for reliable, higher-throughput monitoring to improve the protection of aquatic environments under multiple pressures, guide their restoration, and follow their evolution.

Taxonomic revision of freshwater foraminifera with the description of two new agglutinated species and genera.

Most foraminifera inhabit marine habitats, but some species of monothalamids have been described from freshwater environments, mainly from Swiss water bodies over 100 years ago. Recent environmental DNA surveys revealed the presence of four major phylogenetic clades of freshwater foraminifera. However, until now only one of them (clade 2) has been associated to a morphologically described taxon-the family Reticulomyxidae. Here, we present morphological and molecular data for the genera representing the three remaining clades. We describe two new agglutinated freshwater genera from China and the Netherlands, Lacogromia and Limnogromia, which represent clades 3 and 4, respectively. We also report the first ribosomal DNA sequences of the genus Lieberkuehnia, which place this genus within clade 1. Our study provides the first morphotaxonomic documentation of molecular clades of freshwater foraminifera, showing that the environmental DNA sequences correspond to the agglutinated monothalamous species, morphologically similar to those described 100 years ago.

Taxonomy-free molecular diatom index for high-throughput eDNA biomonitoring.

Current biodiversity assessment and biomonitoring are largely based on the morphological identification of selected bioindicator taxa. Recently, several attempts have been made to use eDNA metabarcoding as an alternative tool. However, until now, most applied metabarcoding studies have been based on the taxonomic assignment of sequences that provides reference to morphospecies ecology. Usually, only a small portion of metabarcoding data can be used due to a limited reference database and a lack of phylogenetic resolution. Here, we investigate the possibility to overcome these limitations using a taxonomy-free approach that allows the computing of a molecular index directly from eDNA data without any reference to morphotaxonomy. As a case study, we use the benthic diatoms index, commonly used for monitoring the biological quality of rivers and streams. We analysed 87 epilithic samples from Swiss rivers, the ecological status of which was established based on the microscopic identification of diatom species. We compared the diatom index derived from eDNA data obtained with or without taxonomic assignment. Our taxonomy-free approach yields promising results by providing a correct assessment for 77% of examined sites. The main advantage of this method is that almost 95% of OTUs could be used for index calculation, compared to 35% in the case of the taxonomic assignment approach. Its main limitations are under-sampling and the need to calibrate the index based on the microscopic assessment of diatoms communities. However, once calibrated, the taxonomy-free molecular index can be easily standardized and applied in routine biomonitoring, as a complementary tool allowing fast and cost-effective assessment of the biological quality of watercourses.

Arnoldiellina fluorescens gen. et sp. nov.--a new green autofluorescent foraminifer from the Gulf of Eilat (Israel).

A new monothalamous (single-chambered) soft-walled foraminiferal species, Arnoldiellina fluorescens gen. et sp. nov., was isolated from samples collected in the Gulf of Eilat, Israel. The species is characterized by a small elongate organic theca with a single aperture of allogromiids. It is characterized by the emission of green autofluorescence (GAF) that has so far not been reported from foraminifera. Phylogenetic analysis of a fragment of the 18S rDNA indicates that the species is related to a group of monothalamous foraminiferans classified as clade I. Although the morphology of the new species is very different compared to the other members of this clade, a specific helix in 18S rRNA secondary structure strongly supports this position.