Assessing pollution of aquatic environments with diatoms’ DNA metabarcoding: experience and developments from France water framework directive networks

Fit for purpose? Evaluating benthic invertebrate DNA metabarcoding for ecological status class assessment in streams under the Water Framework Directive

The ecological state of aquatic ecosystems is systematically monitored using various bioindicators in many countries worldwide. In the European Union, freshwater biomonitoring is the central component of the EU Water Framework Directive (WFD, 2000/60/EC) and currently based on morpho-taxonomic methods. DNA metabarcoding is a novel approach to assess the ecological state fast and efficiently based on organismal DNA signatures and thereby support and upscale biomonitoring. However, compliance of metabarcoding with existing morpho-taxonomic methods must be ensured prior to official implementation. Thus, this study, co-designed by research institutions and environmental agencies, explored necessary key parameters and performed method intercalibration for the implementation of metabarcoding into WFD assessments of running waters. We focussed on benthic invertebrates as the most commonly used bioindicators. We analysed 170 invertebrate samples collected as part of the German federal state WFD routine stream biomonitoring, first via microscopic determination and then using metabarcoding. Our goals were to quantify overlap in i) taxonomic composition and ii) ecological status derived with both methods. For this purpose, we established data harmonisation measures to integrate invertebrate metabarcoding data into the official national WFD classification modules considering abundance and presence/absence data. Our results revealed a high (ca. 70 %) overlap of bioindicator taxa found with both methods. Metabarcoding identified significantly more small invertebrate taxa and detected similar proportions of the important bioindicator 'EPT' taxa (mayflies, stoneflies, caddisflies). Despite deviations in some detected bioindicator taxa, the derived ecological status classes were highly correlated between methods, particularly after intercalibration (R2 = 0.74, Spearman rho = 0.86). Regardless of whether we used abundance or presence/absence data, the resulting stream type classifications showed strong agreement. Thus, our study not only demonstrates the consistency of the methods for the stream types analysed but is also the first to operationalise a path to integration of metabarcoding data into the WFD assessment modules based on formal intercalibration guidelines.

Proficiency testing and cross-laboratory method comparison to support standardisation of diatom DNA metabarcoding for freshwater biomonitoring

DNA metabarcoding of benthic diatoms has been successfully applied for biomonitoring at the national scale and can now be considered technically ready for routine application. However, protocols and methods still vary between and within countries, limiting their transferability and the comparability of results. In order to overcome this, routine use of DNA metabarcoding for diatom biomonitoring requires knowledge of the sources of variability introduced by the different steps of the procedure. Here, we examine how elements of routine procedures contribute to variability between European laboratories. A set of four experiments were performed focusing on DNA extraction and PCR amplification steps to evaluate their reproducibility between different laboratories and the variability introduced by different protocols currently applied by the scientific community. Under the guidance of a reference laboratory, 17 participants from 14 countries performed DNA extraction and PCR amplification in parallel, using the same fixed protocol and their own choice of protocol. Experiments were performed by each participant on a set of standardised DNA and biofilm samples (river, lake and mock community) to investigate potential systematic and random errors. Our results revealed the successful transferability of a protocol amongst labs and a highly similar and consistent ecological assessment outcome obtained regardless of the protocols used by each participant. We propose an "all for one but prove them all" strategy, suggesting that distinct protocols can be used within the scientific community, as long as their consistency is be proven by following minimum standard requirements.

The Diatom Diversity and Ecological Status of a Tufa-Depositing River through eDNA Metabarcoding vs. a Morphological Approach-A Case Study of the Una River (Bosnia and Herzegovina)

Tufa deposits in karst rivers are unique habitats created by mutual interactions between specific environmental and biotope features and inhabited by diatoms as a highly abundant and diverse algal group. This pilot study aimed to investigate the diversity of diatom communities on tufa depositing habitats and assess the Una River's ecological status using a comparative molecular and morphological approach for diatom identification. The 312 base pairs of the rbcL gene were barcoded and analyzed using MiSeq reads and amplicon sequence variants (ASVs) obtained by the DADA2 pipeline. The reference database Diat.barcode v7 was used for taxonomic assignment. The morphological identification of the diatoms was carried out in parallel. In total, the combined dataset revealed 46 taxa identified at genus rank, 125 on the subgenus, and 145 on combined taxonomy rank. The metabarcoding approach mostly leads to a lower number of identified taxa at species rank (58 in molecular vs. 119 in optical inventory), resulting in higher values of beta diversity and heterogeneity in diatom assemblages in samples obtained by morphological approach. Despite the high percentage of taxonomically not assigned diatom ASVs to the species rank, high Shannon diversity index values and a similar number of taxa per locations compared to the morphological approach were obtained. Taxa Achnanthidium minutissimum (Kützing) Czarnecki, Achnanthidium pyrenaicum (Hustedt) H.Kobayasi, Amphora pediculus (Kützing) Grunow, Diatoma vulgaris Bory, Navicula cryptotenella Lange-Bertalot, and Navicula tripunctata (O.F.Müller) Bory were identified at all locations in both inventories. Although limited consistency in the diatom abundances between the two inventory datasets was found, a similar grouping of samples was observed connected to the river's longitudinal gradient. The data obtained using molecular approach in most sites indicated a mostly lower ecological status (good or moderate) compared to the data obtained from the morphological approach (high, good, and moderate). The potential of environmental DNA (eDNA) diatom metabarcoding for water monitoring and diversity studies is undeniable, but to fully realize the benefits of these methods in the future, it is essential to standardize protocols and expand the reference database for species found in specific habitats, such as tufa deposits.

Critical considerations for communicating environmental DNA science

The economic and methodological efficiencies of environmental DNA (eDNA) based survey approaches provide an unprecedented opportunity to assess and monitor aquatic environments. However, instances of inadequate communication from the scientific community about confidence levels, knowledge gaps, reliability, and appropriate parameters of eDNA-based methods have hindered their uptake in environmental monitoring programs and, in some cases, has created misperceptions or doubts in the management community. To help remedy this situation, scientists convened a session at the Second National Marine eDNA Workshop to discuss strategies for improving communications with managers. These include articulating the readiness of different eDNA applications, highlighting the strengths and limitations of eDNA tools for various applications or use cases, communicating uncertainties associated with specified uses transparently, and avoiding the exaggeration of exploratory and preliminary findings. Several key messages regarding implementation, limitations, and relationship to existing methods were prioritized. To be inclusive of the diverse managers, practitioners, and researchers, we and the other workshop participants propose the development of communication workflow plans, using RACI (Responsible, Accountable, Consulted, Informed) charts to clarify the roles of all pertinent individuals and parties and to minimize the chance for miscommunications. We also propose developing decision support tools such as Structured Decision-Making (SDM) to help balance the benefits of eDNA sampling with the inherent uncertainty, and developing an eDNA readiness scale to articulate the technological readiness of eDNA approaches for specific applications. These strategies will increase clarity and consistency regarding our understanding of the utility of eDNA-based methods, improve transparency, foster a common vision for confidently applying eDNA approaches, and enhance their benefit to the monitoring and assessment community.

Genetic Markers for Metabarcoding of Freshwater Microalgae: Review

The metabarcoding methods for studying the diversity of freshwater microalgae and routine biomonitoring are actively used in modern research. A lot of experience has been accumulated already, and many methodological questions have been solved (such as the influence of the methods and time of sample conservation, DNA extraction and bioinformatical processing). The reproducibility of the method has been tested and confirmed. However, one of the main problems-choosing a genetic marker for the study-still lacks a clear answer. We analyzed 70 publications and found out that studies on eukaryotic freshwater microalgae use 12 markers (different nuclear regions 18S and ITS and plastids rbcL, 23S and 16S). Each marker has its peculiarities; they amplify differently and have various levels of efficiency (variability) in different groups of algae. The V4 and V9 18S and rbcL regions are used most often. We concentrated especially on the studies that compare the results of using different markers and microscopy. We summarize the data on the primers for each region and on how the choice of a marker affects the taxonomic composition of a community.

Current Status of Omics in Biological Quality Elements for Freshwater Biomonitoring

Freshwater ecosystems have been experiencing various forms of threats, mainly since the last century. The severity of this adverse scenario presents unprecedented challenges to human health, water supply, agriculture, forestry, ecological systems, and biodiversity, among other areas. Despite the progress made in various biomonitoring techniques tailored to specific countries and biotic communities, significant constraints exist, particularly in assessing and quantifying biodiversity and its interplay with detrimental factors. Incorporating modern techniques into biomonitoring methodologies presents a challenging topic with multiple perspectives and assertions. This review aims to present a comprehensive overview of the contemporary advancements in freshwater biomonitoring, specifically by utilizing omics methodologies such as genomics, metagenomics, transcriptomics, proteomics, metabolomics, and multi-omics. The present study aims to elucidate the rationale behind the imperative need for modernization in this field. This will be achieved by presenting case studies, examining the diverse range of organisms that have been studied, and evaluating the potential benefits and drawbacks associated with the utilization of these methodologies. The utilization of advanced high-throughput bioinformatics techniques represents a sophisticated approach that necessitates a significant departure from the conventional practices of contemporary freshwater biomonitoring. The significant contributions of omics techniques in the context of biological quality elements (BQEs) and their interpretations in ecological problems are crucial for biomonitoring programs. Such contributions are primarily attributed to the previously overlooked identification of interactions between different levels of biological organization and their responses, isolated and combined, to specific critical conditions.

Zeta diversity patterns in metabarcoded lotic algal assemblages as a tool for bioassessment

Assessments of the ecological health of algal assemblages in streams typically focus on measures of their local diversity and classify individuals by morphotaxonomy. Such assemblages are often connected through various ecological processes, such as dispersal, and may be more accurately assessed as components of regional-, rather than local-scale assemblages. With recent declines in the costs of sequencing and computation, it has also become increasingly feasible to use metabarcoding to more accurately classify algal species and perform regional-scale bioassessments. Recently, zeta diversity has been explored as a novel method of constructing regional bioassessments for groups of streams. Here, we model the use of zeta diversity to investigate whether stream health can be determined by the landscape diversity of algal assemblages. We also compare the use of DNA metabarcoding and morphotaxonomy classifications in these zeta diversity-based bioassessments of regional stream health. From 96 stream samples in California, we used various orders of zeta diversity to construct models of biotic integrity for multiple assemblages of diatoms, as well as hybrid assemblages of diatoms in combination with soft-bodied algae, using taxonomy data generated with both DNA sequencing as well as traditional morphotaxonomic approaches. We compared our ability to evaluate the ecological health of streams with the performance of multiple algal indices of biological condition. Our zeta diversity-based models of regional biotic integrity were more strongly correlated with existing indices for algal assemblages classified using metabarcoding compared to morphotaxonomy. Metabarcoding for diatoms and hybrid algal assemblages involved rbcL and 18S V9 primers, respectively. Importantly, we also found that these algal assemblages, independent of the classification method, are more likely to be assembled under a process of niche differentiation rather than stochastically. Taken together, these results suggest the potential for zeta diversity patterns of algal assemblages classified using metabarcoding to inform stream bioassessments.

New recorded diatoms in Holocene sediment cores from the Gonggeom-ji Wetland in Korea

The Gonggeom-ji reservoir is an agricultural one built for rice farming during the Proto-Three Kingdoms period and was designated as Gyeongsangbuk-do monument No. 121 because of its high historical value. The Nakdonggang National Institute of Biological Resources has been conducting paleontological and paleoenvironmental studies on major wetlands from Korea since 2016, as well as diatom, geological, and depth distribution analyses on the sedimentary soil of Gonggeom-ji. This study summarized the description and ecological characteristics of six newly recorded diatoms (Gomphonema lacusrankala, Pinnularia diandae, P. gibba var. hyaline, P. lacunarum, Sellaphora labda var. nipponica, Stauroneis angustilancea) found in samples collected through drilling in Gonggeom-ji in 2019.

Environmental filtering and mass effect are two important processes driving lake benthic diatoms: Results of a DNA metabarcoding study in a large lake

Environmental filtering is often found to dominate assembly rules in diatoms. These microalgae are diverse, especially at subspecies level, and tend to exhibit a niche phylogenetic conservatism. Therefore, other rules, such as competition or mass effects, should be detectable when environmental gradients and dispersal barriers are limited. We used metabarcoding to analyse benthic littoral diatom communities in 153 sites in a large lake (Geneva) exhibiting weak geographical barriers and weak environmental gradients outside river estuaries. We assessed assembly rules using variance partitioning, phylogenetic and source tracking analyses. No phylogenetic over-dispersion of communities, indicative of exclusive competition, was detected. Instead, we found these communities to be phylogenetically over-clustered, indicating environmental filtering, which was even stronger near river estuaries where environmental gradients are stronger. Finally, using a Bayesian method (SourceTracker), we found that rivers flowing into the lake bring communities that settle, especially in sites close to estuaries. Rivers with the highest discharges are primarily responsible for immigration, explaining 27% of lake composition. Therefore, despite favourable conditions to observe other rules, our results support that diatom communities are prominently assembled by environmental filtering and immigration processes, in particular from rivers. However, this does not exclude that other assembly rules may be at play at a finer spatial, temporal and/or phylogenetic scale.

Recommendations for the preservation of environmental samples in diatom metabarcoding studies

Implementation of DNA metabarcoding for diatoms for environmental monitoring is now moving from a research to an operational phase, requiring rigorous guidelines and standards. In particular, the first steps of the diatom metabarcoding process, which consist of sampling and storage, have been addressed in various ways in scientific and pilot studies and now need to be rationalised. The objective of this study was to compare three currently applied preservation protocols through different storage durations (ranging from one day to one year) for phytobenthos and phytoplankton samples intended for diatom DNA metabarcoding analysis. The experimental design used samples from four freshwater and two marine sites of diverse ecological characteristics. The impact of the sample preservation and storage duration was assessed through diatom metabarcoding endpoints: DNA quality and quantity, diversity and richness, diatom assemblage composition and ecological index values (for freshwater samples). The yield and quality of extracted DNA only decreased for freshwater phytobenthos samples preserved with ethanol. Diatom diversity was not affected and their taxonomic composition predominantly reflected the site origin. Only rare taxa (< 100 reads) differed among preservation methods and storage durations. For biomonitoring purposes, freshwater ecological index values were not affected by the preservation method and storage duration tested (including ethanol preservation), all treatments returning the same ecological status for a site. This study contributes to consolidating diatom metabarcoding. Thus, accompanied by operational standards, the method will be ready to be confidently deployed and prescribed in future regulatory monitoring.

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.

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.

How can integrated morphotaxonomy- and metabarcoding-based diatom assemblage analyses best contribute to the ecological assessment of streams?

Environmental conditions, such as nutrient concentrations, salinity, elevation etc., shape diatom assemblages of periphytic biofilms. These assemblages respond rapidly to environmental changes, a fact which makes diatoms valuable bioindicators. Hence, freshwater biomonitoring programmes currently use diatom indices (e.g. EU Water Framework Directive - WFD). To date, microscopy-based assessments require high taxonomic expertise for diatom identification at the species level. High-throughput technologies now provide cost-effective identification approaches that are promising, complementary or alternative tools for bioassessment. The suitability of the metabarcoding method is evaluated for the first time in the Cyprus streams WFD monitoring network, an eastern Mediterranean country with many endemic species and results are compared to the results acquired from the morphotaxonomic analysis. Morphotaxonomic identification was conducted microscopically, using the most updated taxonomic concepts, literature and online resources. At the same time, DNA metabarcoding involved the use of the rbcL 312 bp barcode, high-throughput sequencing and bioinformatic analysis. The ecological status was calculated using the IPS Index. Results show a positive correlation between morpho-taxonomic and molecular IPS scores. Discrepancies between the two methodologies are related to the limitations of both techniques. This study confirmed that Fistulifera saprophila can have a crucial role in key differences observed, as it negatively influences IPS scores and microscopy methods frequently overlook it. Importantly, gaps in the DNA barcoding reference databases lead to a positive overestimation in IPS scores. Overall, we conclude that DNA metabarcoding offsets the morphotaxonomic methodology for the ecological quality assessment of freshwaters.

Potential for cross-contamination of diatom DNA samples when using toothbrushes

The use of toothbrushes and similar devices for sampling diatoms from hard surfaces is a well-established approach. Toothbrushes are routinely cleaned and reused when sampling for analysis by light microscopy. This paper looks at the scale of contamination encountered when this technique is used to sample diatoms for metabarcoding analyses, as well as at the scale of contamination to be expected if stream, rather than distilled water, is used to wash diatoms from stones. Although some contamination attributable to toothbrushes was detected, read numbers were low and had no effect on index calculation or ecological status estimates. However, if the primary focus of a study is to thoroughly document diversity in a sample, then even this small level of contamination may be unacceptable and more stringent measures may be required.

Epiphytic Diatom-Based Biomonitoring in Mediterranean Ponds: Traditional Microscopy versus Metabarcoding Approaches

Benthic diatoms have traditionally been used as bioindicators of aquatic ecosystems. Because diatom-based monitoring of water quality is required by European legislation, molecular-based methods had emerged as useful alternatives to classical methods based on morphological identification using light microscopy. The aim of this study was to test the reliability of DNA metabarcoding combined with High-Throughput Sequencing (HTS) techniques in the bioassessment of the trophic status of 22 Mediterranean shallow ponds in NW Spain. For each pond, the Trophic Diatom Index (TDI) was calculated from inventories obtained by identification using light microscopy (LM) followed by high-throughput sequencing (HTS) at the molecular level. Ponds were subsequently classified into five water quality classes. The results showed a good correspondence between both methods, especially after applying a correction factor that depended on the biovolume of the cells. This correspondence led to the assignment to the same quality class in 59% of the ponds. The determination and quantification of valves or DNA sequences was one of the main pitfalls, which mainly included those related to the variability in the relative abundances of some species. Accordingly, ponds with similar relative abundances for the dominant species were assigned to the same quality class. Moreover, other difficulties leading the discrepancies were the misidentification of some species due to the presence of semi-cryptic taxa, the incompleteness of the reference database and the bioinformatic protocol. Thus, the validation of DNA-based methods for the identification of freshwater diatoms represents an important goal, as an alternative to using traditional methods in Mediterranean shallow ponds.

DNA metabarcoding effectively quantifies diatom responses to nutrients in streams

Nutrient pollution from human activities remains a common problem facing stream ecosystems. Identifying ecological responses to phosphorus and nitrogen can inform decisions affecting the protection and management of streams and their watersheds. Diatoms are particularly useful because they are a highly diverse group of unicellular algae found in nearly all aquatic environments and are sensitive responders to increased nutrient concentrations. Here, we used DNA metabarcoding of stream diatoms as an approach to quantifying effects of total phosphorus (TP) and total nitrogen (TN). Threshold indicator taxa analysis (TITAN) identified operational taxonomic units (OTUs) that increased or decreased along TP and TN gradients along with nutrient concentrations at which assemblages had substantial changes in the occurrences and relative abundances of OTUs. Boosted regression trees showed that relative abundances of gene sequence reads for OTUs identified by TITAN as low P, high P, low N, or high N diatoms had strong relationships with nutrient concentrations, which provided support for potentially using these groups of diatoms as metrics in monitoring programs. Gradient forest analysis provided complementary information by characterizing multi-taxa assemblage change using multiple predictors and results from random forest models for each OTU. Collectively, these analyses showed that notable changes in diatom assemblage structure and OTUs began around 20 µg TP/L, low P diatoms decreased substantially and community change points occurred from 75 to 150 µg/L, and high P diatoms became increasingly dominant from 150 to 300 µg/L. Diatoms also responded to TN with large decreases in low N diatoms occurring from 280 to 525 µg TN/L and a transition to dominance by high N diatoms from 525-850 µg/L. These diatom responses to TP and TN could be used to inform protection efforts (i.e., anti-degradation) and management goals (i.e., nutrient reduction) in streams and watersheds. Our results add to the growing support for using diatom metabarcoding in monitoring programs.

Substratum-associated microbiota

Highlights of new, interesting, and emerging research findings on substratum-associated microbiota covered from a survey of 2019 literature from primarily freshwaters provide insight into research trends of interest to the Water Environment Federation and others interested in benthic, aquatic environments. Coverage of topics on bottom-associated or attached algae and cyanobacteria, though not comprehensive, includes new methods, taxa new-to-science, nutrient dynamics, auto- and heterotrophic interactions, grazers, bioassessment, herbicides and other pollutants, metal contaminants, and nuisance, and bloom-forming and harmful algae. Coverage of bacteria, also not comprehensive, focuses on the ecology of benthic biofilms and microbial communities, along with the ecology of microbes like Caulobacter crescentus, Rhodobacter, and other freshwater microbial species. Bacterial topics covered also include metagenomics and metatranscriptomics, toxins and pollutants, bacterial pathogens and bacteriophages, and bacterial physiology. Readers may use this literature review to learn about or renew their interest in the recent advances and discoveries regarding substratum-associated microbiota. PRACTITIONER POINTS: This review of literature from 2019 on substratum-associated microbiota presents highlights of findings on algae, cyanobacteria, and bacteria from primarily freshwaters. Coverage of algae and cyanobacteria includes findings on new methods, taxa new to science, nutrient dynamics, auto- and heterotrophic interactions, grazers, bioassessment, herbicides and other pollutants, metal contaminants, and nuisance, bloom-forming and harmful algae. Coverage of bacteria includes findings on ecology of benthic biofilms and microbial communities, the ecology of microbes, metagenomics and metatranscriptomics, toxins and pollutants, bacterial pathogens and bacteriophages, and bacterial physiology. Highlights of new, noteworthy and emerging topics build on those from 2018 and will be of relevance to the Water Environment Federation and others interested in benthic, aquatic environments.

Evaluation and sensitivity analysis of diatom DNA metabarcoding for WFD bioassessment of Mediterranean rivers

Our study of 164 diatom samples from Catalonia (NE Spain) is the first to evaluate the applicability of DNA metabarcoding, based on high throughput sequencing (HTS) using a 312-bp rbcL marker, for biomonitoring Mediterranean rivers. For this, we compared the values of a biotic index (IPS) and the ecological status classes derived from them, between light microscope-based (LM) and HTS methods. Very good correspondence between methods gives encouraging results concerning the applicability of DNA metabarcoding for Catalan rivers for the EU Water Framework Directive (WFD). However, in 10 sites, the ecological status class was downgraded from "Good"/"High" obtained by LM to "Moderate"/"Poor"/"Bad" by HTS; these "critical" sites are especially important, because the WFD requires remedial action by water managers for any river with Moderate or lower status. We investigated the contribution of each species to the IPS using a "leave-one-out" sensitivity analysis, paying special attention to critical sites. Discrepancies in IPS between LM and HTS were mainly due to the misidentification and overlooking in LM of a few species, which were better recovered by HTS. This bias was particularly important in the case of Fistulifera saprophila, whose clear underrepresentation in LM was important for explaining 8 out of the 10 critical sites and probably reflected destruction of weakly-silicified frustules during sample preparation. Differences between species in the rbcL copy number per cell affected the relative abundance obtained by HTS for Achnanthidium minutissimum, Nitzschia inconspicua and Ulnaria ulna, which were also identified by the sensitivity analysis as important for the WFD. Only minor IPS discrepancies were attributed to the incompleteness of the reference library, as most of the abundant and influential species (to the IPS) were well represented there. Finally, we propose that leave-one-out analysis is a good method for identifying priority species for isolation and barcoding.

Quantifying stream periphyton assemblage responses to nutrient amendments with a molecular approach

Nutrient (nitrogen [N] and phosphorus [P]) pollution is a pervasive water quality issue in the USA for small streams and rivers. The effect of nutrients on the biotic condition of streams is often evaluated with biological indicators such as macroinvertebrate assemblages or periphyton assemblages, particularly diatoms. Molecular approaches facilitate the use of periphyton assemblages as bioindicators because periphyton is diverse and its composition as a whole, rather than just diatoms, soft-bodied algae, or any single group, may convey additional information about responses to nutrients. To further develop the concept that a taxonomically-broad evaluation of periphyton assemblages could be useful for developing stream bioindicators, we examined microbial assemblage composition with both 16S and 18S rRNA genes, enabling us to evaluate composition in 3 domains. We measured otherwise unknown nutrient responses of different periphyton groups in situ with experiments that used glass fiber filters to allow diffusion of amended nutrients into a stream. We deployed these experimental setups in 2 streams that differ in the extent of agricultural land-use in their catchments in the southeastern USA. Experiments consisted of controls, N amendments, P amendments, and both N and P amendments. Periphyton assemblages that grew on the filters differed significantly by stream, date or season, and nutrient treatment. Assemblage differences across treatments were more consistent among Bacteria and Archaea than among eukaryotes. Effects of nutrient amendments were more pronounced in the stream with less agricultural land use and, therefore, lower nutrient loading than in the stream with more agricultural land use and higher nutrient loading. Combined N and P amendments decreased species richness and evenness for Bacteria and Archaea by ∼36 and ∼9%, respectively, compared with controls. Indicator species analysis revealed that specific clades varied in their response to treatments. Indicators based on the responses of these indicator clades were related to nutrient treatments across sites and seasons. Analyses that included all the taxa in a domain did not resolve differences in responses to N vs P. Instead, better resolution was achieved with an analysis focused on diatoms, which responded more strongly to P than N. Overall, our results showed that in situ nutrient-diffusing substrate experiments are a useful approach for describing assemblage responses to nutrients in streams. This type of molecular approach may be useful to environmental agencies and stakeholders responsible for assessing and managing stream water quality and biotic condition.