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.

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.

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.

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.