Combined Earth observations reveal the sequence of conditions leading to a large algal bloom in Lake Geneva.

Freshwater algae exhibit complex dynamics, particularly in meso-oligotrophic lakes with sudden and dramatic increases in algal biomass following long periods of low background concentration. While the fundamental prerequisites for algal blooms, namely light and nutrient availability, are well-known, their specific causation involves an intricate chain of conditions. Here we examine a recent massive Uroglena bloom in Lake Geneva (Switzerland/France). We show that a certain sequence of meteorological conditions triggered this specific algal bloom event: heavy rainfall promoting excessive organic matter and nutrients loading, followed by wind-induced coastal upwelling, and a prolonged period of warm, calm weather. The combination of satellite remote sensing, in-situ measurements, ad-hoc biogeochemical analyses, and three-dimensional modeling proved invaluable in unraveling the complex dynamics of algal blooms highlighting the substantial role of littoral-pelagic connectivities in large low-nutrient lakes. These findings underscore the advantages of state-of-the-art multidisciplinary approaches for an improved understanding of dynamic systems as a whole.

Morphology and phylogeny of Nitzschianandorii sp. nov. (Bacillariophyceae), a new small-celled lanceolate species from a post-mining reservoir.

Post-mining reservoirs are distinguished by characteristic environmental conditions where specific diatom communities can be observed. Reservoirs created as a part of the reclamation plan after human mining activities are marked by unique chemical and physical water parameters. In the course of research on the diatoms from Bogdalów reservoir, we examined the taxonomic and morphological diversity of Nitzschia taxa from the section Lanceolatae occurring in a post-mining lignite reservoir. Our study describes a new species of Nitzschia from a post-mining reservoir, Nitzschianandorii Olszynski, Zakrzewski & Zelazna-Wieczorek, sp. nov. Morphometry and morphology analyses of new species were performed with light and scanning electron microscopy. Chloroplast morphology analysis was conducted with differential interference contrast microscopy and confocal laser scanning microscopy. Molecular data from SSU 18S, rbcL and psbC sequences were obtained from cultures of this taxon. Differential diagnosis of Nitzschianandorii Olszynski, Zakrzewski & Zelazna-Wieczorek, sp. nov. with co-occurring taxa: N.lacuum and N.alpinobacillum was performed using morphological traits and nMDS analysis of the valves' morphometry.

Assessing the relevance of DNA metabarcoding compared to morphological identification for lake phytoplankton monitoring.

Phytoplankton is a key biological group used to assess the ecological status of lakes. The classical monitoring approach relies on microscopic identification and counting of phytoplankton species, which is time-consuming and requires high taxonomic expertise. High-throughput sequencing, combined with metabarcoding, has recently demonstrated its potential as an alternative approach for plankton surveys. Several studies have confirmed the relevance of the diatom metabarcoding approach to calculate biotic indices based on species ecology. However, phytoplankton communities have not yet benefited from such validation. Here, by comparing the results obtained with the two methods (molecular and microscopic counting), we evaluated the relevance of metabarcoding approach for phytoplankton monitoring by considering different metrics: alpha diversity, taxonomic composition, community structure and a phytoplankton biotic index used to assess the trophic level of lakes. For this purpose, 55 samples were collected in four large alpine lakes (Aiguebelette, Annecy, Bourget, Geneva) during the year 2021. For each sample, a metabarcoding analysis based on two genetic markers (16S and 23S rRNA) was performed, in addition to the microscopic count. Regarding the trophic level of lakes, significant differences were found between index values obtained with the two approaches. The main hypothesis to explain these differences comes from the incompleteness, particularly at the species level, of the barcode reference library for the two genetic markers. It is therefore necessary to complete reference libraries for using such species-based biotic indices with metabarcoding data. Besides this, species richness and diversity were higher in the molecular inventories than in the microscopic ones. Moreover, despite differences in taxonomic composition of the floristic lists obtained by the two approaches, their community structures were similar. These results support the possibility of using metabarcoding for phytoplankton monitoring but in a different way. We suggest exploring alternative approaches to index development, such as a taxonomy-free approach.

Diatom endemism and taxonomic turnover: Assessment in high-altitude alpine lakes covering a large geographical range.

Diatoms are widely used as ecological indicators and show various degrees of endemism. Many studies that support the idea of endemic species integrate several climate zones, a variety of ecosystem types, and often focus on a global scale. Here, we investigated whether endemism could be detected when considering a homogeneous type of ecosystem in a single climate zone. We sampled stone biofilms at 40-50 cm depth in high-altitude lakes in the Alpine climate zone. A total of 149 samples were obtained from the French and Georgian mountains, two areas separated by ∼3000 km. Using Amplicon Sequence Variants derived from DNA metabarcoding, we assessed taxonomic turnover and Zeta-diversity (a measure of endemism). We ran haplotype networks and phylogenetic tests to measure geographical signal in the phylogenies of dominant taxa. The French and Georgian communities shared 51 % of species. Species that were not shared across both regions were mostly rare, and often not characteristic of lakes but of neighboring habitats instead. In contrast, at the sub-species level, 87 % of the genotypes showed restricted distributions. Whereas endemism was the rule at sub-species level, most species were shared across both French and Georgian lakes, suggesting that geographic barriers strongly limited dispersal at the sub-species level but not species level. Dominant species hosted higher levels of sub-specific diversity than rare species. In contrast to global-scale studies, we did not find any significant geographical structuring in the phylogeny of the investigated species. This could indicate ongoing dispersal at a frequency fast enough to prevent allopatric divergence, yet slow enough to prevent sharing most haplotypes between France and Georgia. These results have implications for biomonitoring: depending on the taxonomic level chosen, robust generic tools (species level) or tools dedicated to a region able to discriminate fine pressures differences (sub-species level) may be developed.

Diatom DNA metabarcoding to assess the effect of natural radioactivity in mineral springs on ASV of benthic diatom communities.

Little is still known about the low dose effects of radiation on the microbial communities in the environment. Mineral springs are ecosystems than can be affected by natural radioactivity. These extreme environments are, therefore, observatories for studying the influence of chronic radioactivity on the natural biota. In these ecosystems we find diatoms, unicellular microalgae, playing an essential role in the food chain. The present study aimed to investigate, using DNA metabarcoding, the effect of natural radioactivity in two environmental compartments (i.e. spring sediments and water) on the genetic richness, diversity and structure of diatom communities in 16 mineral springs in the Massif Central, France. Diatom biofilms were collected during October 2019, and a 312 bp region of the chloroplast gene rbcL (coding for the Ribulose Bisphosphate Carboxylase) used as a barcode for taxonomic assignation. A total of 565 amplicon sequence variants (ASV) were found. The dominant ASV were associated with Navicula sanctamargaritae, Gedaniella sp., Planothidium frequentissimum, Navicula veneta, Diploneis vacillans, Amphora copulata, Pinnularia brebissonii, Halamphora coffeaeformis, Gomphonema saprophilum, and Nitzschia vitrea, but some of the ASVs could not be assigned at the species level. Pearson correlation failed to show a correlation between ASV' richness and radioactivity parameters. Non-parametric MANOVA analysis based on ASVs occurrence or abundances revealed that geographical location was the main factor influencing ASVs distribution. Interestingly, 238U was the second factor that explained diatom ASV structure. Among the ASVs in the mineral springs monitored, ASV associated with one of the genetic variants of Planothidium frequentissimum was well represented in the springs and with higher levels of 238U, suggesting its high tolerance to this particular radionuclide. This diatom species may therefore represent a bio-indicator of high natural levels of uranium.

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.

DNA sequence and taxonomic gap analyses to quantify the coverage of aquatic cyanobacteria and eukaryotic microalgae in reference databases: Results of a survey in the Alpine region.

The taxonomic identification of organisms based on the amplification of specific genetic markers (metabarcoding) implicitly requires adequate discriminatory information and taxonomic coverage of environmental DNA sequences in taxonomic databases. These requirements were quantitatively examined by comparing the determination of cyanobacteria and microalgae obtained by metabarcoding and light microscopy. We used planktic and biofilm samples collected in 37 lakes and 22 rivers across the Alpine region. We focused on two of the most used and best represented genetic markers in the reference databases, namely the 16S rRNA and 18S rRNA genes. A sequence gap analysis using blastn showed that, in the identity range of 99-100%, approximately 30% (plankton) and 60% (biofilm) of the sequences did not find any close counterpart in the reference databases (NCBI GenBank). Similarly, a taxonomic gap analysis showed that approximately 50% of the cyanobacterial and eukaryotic microalgal species identified by light microscopy were not represented in the reference databases. In both cases, the magnitude of the gaps differed between the major taxonomic groups. Even considering the species determined under the microscope and represented in the reference databases, 22% and 26% were still not included in the results obtained by the blastn at percentage levels of identity ≥95% and ≥97%, respectively. The main causes were the absence of matching sequences due to amplification and/or sequencing failure and potential misidentification in the microscopy step. Our results quantitatively demonstrated that in metabarcoding the main obstacles in the classification of 16S rRNA and 18S rRNA sequences and interpretation of high-throughput sequencing biomonitoring data were due to the existence of important gaps in the taxonomic completeness of the reference databases and the short length of reads. The study focused on the Alpine region, but the extent of the gaps could be much greater in other less investigated geographic areas.

Fish eDNA metabarcoding from aquatic biofilm samples: Methodological aspects.

Fish eDNA metabarcoding is usually performed from filtered water samples. The volume of filtered water depends on the study scope and can rapidly become time consuming according to the number of samples that have to be processed. To avoid time allocated to filtration, passive DNA samplers have been used to recover fish eDNA from marine environments faster. In freshwater ecosystems, aquatic biofilms were used to catch eDNA from macroinvertebrates. Here, we test the capacity of aquatic biofilms to entrap fish eDNA in a large lake and, therefore, the possibility to perform fish eDNA metabarcoding from this matrix compared to the traditional fish eDNA approach from filtered water samples. Methodological aspects of the use of aquatic biofilms for fish eDNA metabarcoding (e.g. PCR replicates, biological replicates, bioinformatics pipeline, reference database and taxonomic assignment) were validated against a mock community. When using biofilms from habitats sheltered from wind and waves, biofilm and water approach provided similar inventories. Richness and diversity were comparable between both approaches. Approaches differed only for rare taxa. Our results illustrate the capacity of aquatic biofilms to act as passive eDNA samplers of fish eDNA and, therefore, the possibility to use biofilms to monitor fish communities efficiently from biofilms. Furthermore, our results open up avenues of research to study a diversity of biological groups (among which bioindicators as diatoms, macroinvertebrates and fish) from eDNA isolated from a single environmental matrix reducing sampling efforts, analysis time and costs.

DNA metabarcoding reveals differences in distribution patterns and ecological preferences among genetic variants within some key freshwater diatom species.

Our study evaluates differences in the distribution and ecology of genetic variants within several ecologically important diatom species that are also key for Water Framework Directive monitoring of European rivers: Fistulifera saprophila (FSAP), Achnanthidium minutissimum (ADMI), Nitzschia inconspicua (NINC) and Nitzschia soratensis (NSTS). We used DADA2 to infer amplicon sequence variants (ASVs) of a short rbcL barcode in 531 environmental samples from biomonitoring campaigns in Catalonia and France. ASVs within each species showed different distribution patterns. Threshold Indicator Taxa ANalysis revealed three ecological groupings of ASVs in both ADMI and FSAP. Two of these in each species were separated by opposite responses to calcium and conductivity. Boosted regression trees additionally showed that both variables greatly influenced the occurrence of these groupings. A third grouping in FSAP was characterized by a negative response to total organic carbon and hence was better represented in waters with higher ecological status than the other FSAP ASVs, contrasting with what is generally assumed for the species. In the two Nitzschia species, our analyses confirmed earlier studies: NINC preferred higher levels of calcium and conductivity. Our findings suggest that the broad ecological tolerance of some diatom species results from overlapping preferences among genetic variants, which individually show much more restricted preferences and distributions. This work shows the importance of studying the ecological preferences of genetic variants within species complexes, now possible with DNA metabarcoding. The results will help reveal and understand biogeographical distributions and facilitate the development of more accurate biological indexes for biomonitoring programmes.

Cyprus' diatom diversity and the association of environmental and anthropogenic influences for ecological assessment of rivers using DNA metabarcoding.

Human activities are the leading cause of environmental impairments. Appropriate biomonitoring of ecosystems is needed to assess these activities effectively. In freshwater ecosystems, periphytic and epilithic biofilms have diatom assemblages. These assemblages respond rapidly to environmental changes, making diatoms valuable bioindicators. For this reason, freshwater biomonitoring programs are currently using diatoms (e.g., Water Framework Directive). In the past ten years, DNA metabarcoding coupled with next-generation sequencing and bioinformatics represents a complementary approach for diatom biomonitoring. In this study, this approach is used for the first time in Cyprus by considering the association of environmental and anthropogenic pressures to diatom assemblages. Statistical analysis was then applied to identify the environmental (i.e., river types, geo-morphological) and anthropogenic (i.e., physicochemical, human land-use pressures) variables' role in the observed diatom diversity. Results indicate differences in diatom assemblages between intermittent and perennial rivers. Achnanthidium minutissimum was more abundant in intermittent rivers; whereas Amphora pediculus and Planothidium caputium in perennial ones. Additionally, we could demonstrate the correlation between nutrients (e.g., nitrogen, phosphorus), stations' local characteristics (e.g., elevation), and land use activities on the observed differences in diatom diversity. Finally, we conclude that multi-stressors and anthropogenic pressures together as multiple stressors have a significant statistical relationship to the observed diatom diversity and play a pivotal role in determining Cyprus' rivers' ecological status.

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.

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.

Ripe for reassessment: A synthesis of available molecular data for the speciose diatom family Bacillariaceae.

The Bacillariaceae is a very species-rich family of raphid diatoms and includes the large and taxonomically difficult genus Nitzschia, whose species are often small-celled and finely structured and have few discrete morphological characters visible in the light microscope. The classification of Nitzschia is still mostly based on one developed in the second half of the 19th century by Grunow, who separated the genus into a series of sections largely on cell shape and symmetry, the position of the raphe, transverse extension of the fibulae, and folding of the valve. We assembled and analysed single-gene and concatenated alignments of nSSU, nLSU, rbcL, psbC and cox1 to test Grunow's and subsequent classifications and to examine selected morphological characters for their potential to help define monophyletic groups. The maximum likelihood trees were equivocal as to monophyly of the family itself but showed good support for each of eight main clades of Bacillariaceae, three of which corresponded more or less to existing genera (Hantzschia, Cylindrotheca and Bacillaria). The other five main clades and some subclades comprised groups of Nitzschia species or assemblies of Nitzschia species with other genera (Pseudo-nitzschia, Fragilariopsis, Neodenticula, Tryblionella, Psammodictyon). Relationships between most of the eight main clades were not resolved robustly but all analyses recovered Nitzschia as non-monophyletic. The Grunowian classification of Nitzschia into sections was not supported, though in some respects (e.g. treatment of sigmoid species) it is better than subsequent reclassifications. Several of the main clades and subclades are cryptic (lacking morphological synapomorphies) and homoplasy is common in both light microscopical and ultrastructural characters (to the extent that organisms initially assigned to the same species sometimes prove to belong to a different main clade). Nevertheless, some characters, including the structure of the raphe canal and girdle, seem to be sufficiently conservative evolutionarily to give a provisional estimate of relationships if molecular data are unavailable. No new formal classifications are proposed but various options are explored and research needs identified.

Environmental filtering and limiting similarity as main forces driving diatom community structure in Mediterranean and continental temporary and perennial streams.

Climatic extreme events such as droughts (unpredictable), dry periods (predictable) or even flush floods, threaten freshwater ecosystems worldwide. The filtering mechanisms of these events and their strength on communities, however, can be different among regions. While time-for-adaptation theory defines whether or not water scarcity can be considered as disturbance, the stress-dominance theory predicts an increase in importance of environmental filtering and a decrease in the role of biotic interactions in communities with increasing environmental stress. Here, we tested whether environmental filtering (leading to trait convergence) or limiting similarity (leading to trait divergence) is the main assembly rule shaping the structure and trait composition of benthic diatom assemblages in Mediterranean (Portuguese) and continental (Hungarian) temporary and perennial streams. We assumed that the trait composition of diatom assemblages in the two stream types would be less different in the Mediterranean than in the continental region (addressed to time-for-adaptation theory). We also hypothesized that trait composition would be shaped by environmental filtering in the Hungarian streams while by biotic interactions in Portuguese streams (addressed to stress-dominance theory). Our results supported our first hypothesis since traits, which associated primarily to temporary streams were found only in the continental region. Our findings, however, only partially proved the stress-dominance hypothesis. In the continental region, where drying up of streams were induced by unpredictable droughts, biotic interactions were the main assembly rules shaping community structure. In contrast, environmental filtering was nearly as important as limiting similarity in structuring trait composition in the Mediterranean region during the predictable dry phase with no superficial flow. These analyses also highlighted that drought events (both predictable and unpredictable ones) have a complex and strong influence on benthic diatom assemblages resulting even in irreversible changes in trait composition and thereby in ecosystem functioning.

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.

A taxonomy-free approach based on machine learning to assess the quality of rivers with diatoms.

Diatoms are a compulsory biological quality element in the ecological assessment of rivers according to the Water Framework Directive. The application of current official indices requires the identification of individuals to species or lower rank under a microscope based on the valve morphology. This is a highly time-consuming task, often susceptible of disagreements among analysts. In alternative, the use of DNA metabarcoding combined with High-Throughput Sequencing (HTS) has been proposed. The sequences obtained from environmental DNA are clustered into Operational Taxonomic Units (OTUs), which can be assigned to a taxon using reference databases, and from there calculate biotic indices. However, there is still a high percentage of unassigned OTUs to species due to the incompleteness of reference libraries. Alternatively, we tested a new taxonomy-free approach based on diatom community samples to assess rivers. A combination of three machine learning techniques is used to build models that predict diatom OTUs expected in test sites, under reference conditions, from environmental data. The Observed/Expected OTUs ratio indicates the deviation from reference condition and is converted into a quality class. This approach was never used with diatoms neither with OTUs data. To evaluate its efficiency, we built a model based on OTUs lists (HYDGEN) and another based on taxa lists from morphological identification (HYDMORPH), and also calculated a biotic index (IPS). The models were trained and tested with data from 81 sites (44 reference sites) from central Portugal. Both models were considered accurate (linear regression for Observed and Expected richness: R2 ≈ 0.7, interception ≈ 0.8) and sensitive to global anthropogenic disturbance (Rs2 > 0.30 p < 0.006 for global disturbance). Yet, the HYDGEN model based on molecular data was sensitive to more types of pressures (such as, changes in land use and habitat quality), which gives promising insights to its use for bioassessment of rivers.

Diat.barcode, an open-access curated barcode library for diatoms.

Diatoms (Bacillariophyta) are ubiquitous microalgae which produce a siliceous exoskeleton and which make a major contribution to the productivity of oceans and freshwaters. They display a huge diversity, which makes them excellent ecological indicators of aquatic ecosystems. Usually, diatoms are identified using characteristics of their exoskeleton morphology. DNA-barcoding is an alternative to this and the use of High-Throughput-Sequencing enables the rapid analysis of many environmental samples at a lower cost than analyses under microscope. However, to identify environmental sequences correctly, an expertly curated reference library is needed. Several curated libraries for protists exists; none, however are dedicated to diatoms. Diat.barcode is an open-access library dedicated to diatoms which has been maintained since 2012. Data come from two sources (1) the NCBI nucleotide database and (2) unpublished sequencing data of culture collections. Since 2017, several experts have collaborated to curate this library for rbcL, a chloroplast marker suitable for species-level identification of diatoms. For the latest version of the database (version 7), 605 of the 3482 taxonomical names originally assigned by the authors of the rbcL sequences were modified after curation. The database is accessible at https://www6.inra.fr/carrtel-collection_eng/Barcoding-database .

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.

Connecting the morphological and molecular species concepts to facilitate species identification within the genus Fragilaria (Bacillariophyta).

This paper explores the diversity and taxonomy of species within Fragilaria sensu stricto, an abundant and ecologically important diatom genus, taking advantage of cultured and DNA-barcoded material. The goal is to facilitate the identification of European taxa within this complex, providing a unified view on morphological and molecular diversity. There is a general agreement that the separation of species within the group of Fragilaria is difficult because morphological descriptions of species are not consistent between authorities, ongoing taxonomic revisions have resulted in species described with standards of the late 20th and 21st centuries alongside descriptions based on 19th century (light microscopical) criteria, and because not all diagnostic characters can be seen in all specimens encountered in routine analyses. Consequent confusion could blur potentially important ecological distinctions between species. Our study demonstrated that some species defined on morphological criteria could be confirmed using the rbcL chloroplast gene as a genetic marker, for example, Fragilaria gracilis, Fragilaria tenera, Fragilaria perminuta, and Fragilaria subconstricta. However, even for those species, preliminary identifications based on morphology often differed from identifications based on phylogenetic clustering combined with detailed morphological study. Clades were well-defined by rbcL, but based on morphology, the terminal taxa of these clades did not match the currently described Fragilaria species. To clarify recognition of these taxa, we describe three new species: Fragilaria agnesiae, Fragilaria heatherae, and Fragilaria joachimii.

Benthic Diatom Communities in an Alpine River Impacted by Waste Water Treatment Effluents as Revealed Using DNA Metabarcoding.

Freshwater ecosystems are continuously affected by anthropogenic pressure. One of the main sources of contamination comes from wastewater treatment plant (WWTP) effluents that contain wide range of micro- and macropollutants. Chemical composition, toxicity levels and impact of treated effluents (TEs) on the recipient aquatic ecosystems may strongly differ depending on the wastewater origin. Compared to urban TEs, hospital ones may contain more active pharmaceutical substances. Benthic diatoms are relevant ecological indicators because of their high species and ecological diversity and rapid response to human pressure. They are routinely used for water quality monitoring. However, there is a knowledge gap on diatom communities' development and behavior in treated wastewater in relation to prevailing micro- and macropollutants. In this study, we aim to (1) investigate the response of diatom communities to urban and hospital TEs, and (2) evaluate TEs effect on communities in the recipient river. Environmental biofilms were colonized in TEs and the recipient river up- and downstream from the WWTP output to study benthic diatoms using DNA metabarcoding combined with high-throughput sequencing (HTS). In parallel, concentrations of nutrients, pharmaceuticals and seasonal conditions were recorded. Diatom metabarcoding showed that benthic communities differed strongly in their diversity and structure depending on the habitat. TE sites were generally dominated by few genera with polysaprobic preferences belonging to the motile guild, while river sites favored diverse communities from oligotrophic and oligosaprobic groups. Seasonal changes were visible to lower extent. To categorize parameters important for diatom changes we performed redundancy analysis which suggested that communities within TE sites were associated to higher concentrations of beta-blockers and non-steroidal anti-inflammatory drugs in urban effluents vs. antibiotics and orthophosphate in hospital effluents. Furthermore, indicator species analysis showed that 27% of OTUs detected in river downstream communities were indicator for urban or hospital TE sites and were absent in the river upstream. Finally, biological diatom index (BDI) calculated to evaluate the ecological status of the recipient river suggested water quality decrease linked to the release of TEs. Thus, in-depth assessment of diatom community composition using DNA metabarcoding is proposed as a promising technique to highlight the disturbing effect of pollutants in Alpine rivers.