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

Potential effects of the discharge of wastewater treatment plant (WWTP) effluents in benthic communities: evidence from three distinct WWTP systems

Wastewater treatment plant (WWTP) effluents can be sources of environmental contamination. In this study, we aimed to understand whether effluents of three different WWTPs may have ecological effects in riverine recipient ecosystems. To achieve this, we assessed benthic phytobenthos and macroinvertebrate communities at three different locations relative to the effluent discharge: immediately upstream, immediately downstream and 500-m downstream the effluent discharge. Two approaches were employed: the ecological status classification as defined in the Water Framework Directive (WFD) based on biological indicators; constrained multivariate analysis to disentangle the environmental drivers (physicochemical variables and contaminants, namely metals, polycyclic aromatic hydrocarbons, pharmaceuticals, and personal care products) of ecological changes across the study sites. The results showed inconsistencies between the WFD approach and the multivariate approach, as well as between the responses of macroinvertebrates and diatoms. The WWTP effluents impacted benthic communities in a single case: macroinvertebrates were negatively affected by one of the WWTP effluents, likely by the transported pharmaceuticals (other stressors are essentially homogeneous among sites). Given the findings and the scarcity of consistent evidence on ecological impacts that WWTP effluents may have in recipient ecosystems, further research is needed towards more sustainable regulation and linked environmental protection measures.

Successive accumulation of biotic assemblages at a fine spatial scale along glacier-fed waters

Glacier-fed waters create strong environmental filtering for biota, whereby different organisms may assume distinct distribution patterns. By using environmental DNA-based metabarcoding, we investigated the multi-group biodiversity distribution patterns of the Parlung No. 4 Glacier, on the Tibetan Plateau. Altogether, 642 taxa were identified from the meltwater stream and the downstream Ranwu Lake, including 125 cyanobacteria, 316 diatom, 183 invertebrate, and 18 vertebrate taxa. As the distance increased from the glacier terminus, community complexity increased via sequential occurrences of cyanobacteria, diatoms, invertebrates, and vertebrates, as well as increasing taxa numbers. The stream and lake showed different community compositions and distinct taxa. Furthermore, the correlations with environmental factors and community assembly mechanisms showed group- and habitat-specific patterns. Our results reveal the rapid spatial succession and increasing community complexity along glacial flowpaths and highlight the varying adaptivity of different organisms, while also providing insight into the ecosystem responses to global change.

Towards the fate of antibiotics and the development of related resistance genes in stream biofilms

Antibiotics are ubiquitously found in natural surface waters and cause great harm to aquatic organisms. Stream biofilm is a complex and active community composed of algae, bacteria, fungi and other microorganisms, which mainly adheres to solid substances such as rocks and sediments. The durability and diverse structural and metabolic characteristics of biofilms make them a representative of microbial life in aquatic micrecosystems and can reflect major ecosystem processes. Microorganisms and extracellular polymeric substances in biofilms can adsorb and actively accumulate antibiotics. Therefore, biofilms are excellent biological indicators for detecting antibiotic in polluted aquatic environments, but the biotransformation potential of stream biofilms for antibiotics has not been fully explored in the aquatic environment. The characteristics of stream biofilm, such as high abundance and activity of bacterial community, wide contact area with pollutants, etc., which increases the opportunity of biotransformation of antibiotics in biofilm and contribute to bioremediation to improve ecosystem health. Recent studies have demonstrated that both exposure to high and sub-minimum inhibitory concentrations of antibiotics may drive the development of antibiotic resistance genes (ARGs) in natural stream biofilms, which are susceptible to the effects of antibiotic residues, microbial communities and mobile genetic elements, etc. On the basis of peer-reviewed papers, this review explores the distribution behavior of antibiotics in stream biofilms and the contribution of biofilms to the acquisition and spread of antibiotic resistance. Considering that antibiotics and ARGs alter the structure and ecological functions of natural microbial communities and pose a threat to river organisms and human health, our research findings provide comprehensive insights into the migration, transformation, and bioavailability of antibiotics in biofilms.

Insights in Pharmaceutical Pollution: The Prospective Role of eDNA Metabarcoding

Environmental pollution is a growing threat to natural ecosystems and one of the world's most pressing concerns. The increasing worldwide use of pharmaceuticals has elevated their status as significant emerging contaminants. Pharmaceuticals enter aquatic environments through multiple pathways related to anthropogenic activity. Their high consumption, insufficient waste treatment, and the incapacity of organisms to completely metabolize them contribute to their accumulation in aquatic environments, posing a threat to all life forms. Various analytical methods have been used to quantify pharmaceuticals. Biotechnology advancements based on next-generation sequencing (NGS) techniques, like eDNA metabarcoding, have enabled the development of new methods for assessing and monitoring the ecotoxicological effects of pharmaceuticals. eDNA metabarcoding is a valuable biomonitoring tool for pharmaceutical pollution because it (a) provides an efficient method to assess and predict pollution status, (b) identifies pollution sources, (c) tracks changes in pharmaceutical pollution levels over time, (d) assesses the ecological impact of pharmaceutical pollution, (e) helps prioritize cleanup and mitigation efforts, and (f) offers insights into the diversity and composition of microbial and other bioindicator communities. This review highlights the issue of aquatic pharmaceutical pollution while emphasizing the importance of using modern NGS-based biomonitoring actions to assess its environmental effects more consistently and effectively.

Sub-MIC antibiotics influence the microbiome, resistome and structure of riverine biofilm communities

The effects of sub-minimum inhibitory concentrations (sub-MICs) of antibiotics on aquatic environments is not yet fully understood. Here, we explore these effects by employing a replicated microcosm system fed with river water where biofilm communities were continuously exposed over an eight-week period to sub-MIC exposure (1/10, 1/50, and 1/100 MIC) to a mix of common antibiotics (ciprofloxacin, streptomycin, and oxytetracycline). Biofilms were examined using a structure-function approach entailing microscopy and metagenomic techniques, revealing details on the microbiome, resistome, virulome, and functional prediction. A comparison of three commonly used microbiome and resistome databases was also performed. Differences in biofilm architecture were observed between sub-MIC antibiotic treatments, with an overall reduction of extracellular polymeric substances and autotroph (algal and cyanobacteria) and protozoan biomass, particularly at the 1/10 sub-MIC condition. While metagenomic analyses demonstrated that microbial diversity was lowest at the sub-MIC 1/10 antibiotic treatment, resistome diversity was highest at sub-MIC 1/50. This study also notes the importance of benchmarking analysis tools and careful selection of reference databases, given the disparity in detected antimicrobial resistance genes (ARGs) identity and abundance across methods. Ultimately, the most detected ARGs in sub-MICs exposed biofilms were those that conferred resistance to aminoglycosides, tetracyclines, β-lactams, sulfonamides, and trimethoprim. Co-occurrence of microbiome and resistome features consistently showed a relationship between Proteobacteria genera and aminoglycoside ARGs. Our results support the hypothesis that constant exposure to sub-MICs antibiotics facilitate the transmission and promote prevalence of antibiotic resistance in riverine biofilms communities, and additionally shift overall microbial community metabolic function.

Emerging challenges of the impacts of pharmaceuticals on aquatic ecosystems: A diatom perspective

Pharmaceuticals are a ubiquitous group of emerging pollutants of considerable importance due to their biological potency and potential to elicit effects in wildlife and humans. Pharmaceuticals have been quantified in terrestrial, marine, fresh, and transitional waters, as well as the fauna and macro-flora that inhabit them. Pharmaceuticals can enter water ways through different human and veterinary pathways with traditional wastewater treatment, unable to completely remove pharmaceuticals, discharging often unknown quantities to aquatic ecosystems. However, there is a paucity of available information regarding the effects of pharmaceuticals on species at the base of aquatic food webs, especially on phytoplankton, with research typically focussing on fish and aquatic invertebrates. Diatoms are one of the main classes of phytoplankton and are some of the most abundant and important organisms in aquatic systems. As primary producers, diatoms generate ∼40 % of the world's oxygen and are a vital food source for primary consumers. Diatoms can also be used for bioremediation of polluted water bodies but perhaps are best known as bio-indicators for water quality studies. However, this keystone, non-target group is often ignored during ecotoxicological studies to assess the effects of pollutants of concern. Observed effects of pharmaceuticals on diatoms have the potential to be used as an indicator of pharmaceutical-induced impacts on higher trophic level organisms and wider ecosystem effects. The aim of this review is to present a synthesis of research on pharmaceutical exposure to diatoms, considering ecotoxicity, bioremediation and the role of diatoms as bio-indicators. We highlight significant omissions and knowledge gaps which need addressing to realise the potential role of diatoms in future risk assessment approaches and help evaluate the impacts of pharmaceuticals in the aquatic environment at local and global scales.

Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities

Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae, and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable, and highly reproducible in terms of microbial composition, function, and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation.

Wastewater microorganisms impact microbial diversity and important ecological functions of stream periphyton

Effluents of wastewater treatment plants can impact microbial communities in the receiving streams. However, little is known about the role of microorganisms in wastewater as opposed to other wastewater constituents, such as nutrients and micropollutants. We aimed therefore at determining the impact of wastewater microorganisms on the microbial diversity and function of periphyton, key microbial communities in streams. We used a flow-through channel system to grow periphyton upon exposure to a mixture of stream water and unfiltered or ultra-filtered wastewater. Impacts were assessed on periphyton biomass, activities and tolerance to micropollutants, as well as on microbial diversity. Our results showed that wastewater microorganisms colonized periphyton and modified its community composition, resulting for instance in an increased abundance of Chloroflexi and a decreased abundance of diatoms and green algae. This led to shifts towards heterotrophy, as suggested by the changes in nutrient stoichiometry and the increased mineralization potential of carbon substrates. An increased tolerance towards micropollutants was only found for periphyton exposed to unfiltered wastewater but not to ultra-filtered wastewater, suggesting that wastewater microorganisms were responsible for this increased tolerance. Overall, our results highlight the need to consider the role of wastewater microorganisms when studying potential impacts of wastewater on the receiving water body.

Community reassemblies of eukaryotes, prokaryotes, and viruses in the hexabromocyclododecanes-contaminated microcosms

The microbial community in seriously contaminated environment were not well known. This research investigated the community reassemblies in microcosms made of two distinct mangrove sediments amended with high levels of hexabromocyclododecanes (HBCDs). After eight months of contamination, the transformation of HBCDs yielded various lower brominated products and resulted in acidification (pH ~2). Therefore, the degraders and dehalogenase homologous genes involved in transformation of HBCDs only presented in low abundance to avoid further deterioration of the habitats. Moreover, in these deteriorated habitats, 1344 bacterial, 969 archaeal, 599 eukaryotic (excluded fungi), 187 fungal OTUs, and 10 viral genera, were reduced compared with controls. Specifically, in two groups of microcosms, Zetaproteobacteria, Deinococcus-Thermus, Spirochaetes, Bacteroidetes, Euryarchaeota, and Ascomycota, were positively responding taxa to HBCDs. Caloneis (Bacillariophyta) and Ascomycota turned to the dominant eukaryotic and fungal taxa. Most of predominant taxa were related to the contamination of brominated flame retardants (BFRs). Microbial communities were reassembled in divergent and sediment-dependent manner. The long-term contamination of HBCDs leaded to the change of relations between many taxa, included some of the environmental viruses and their known hosts. This research highlight the importance of monitoring the ecological effects around plants producing or processing halogenated compounds.

Wastewater microorganisms impact the micropollutant biotransformation potential of natural stream biofilms

Biotransformation is the most important process removing manmade chemicals from the environment, yet mechanisms governing this essential ecosystem function are underexplored. To understand these mechanisms, we conducted experiments in flow-through systems, by colonizing stream biofilms under different conditions of mixing river water with treated (and ultrafiltered) wastewater. We performed biotransformation experiments with those biofilms, using a set of 75 micropollutants, and could disentangle potential mechanisms determining the biotransformation potential of stream biofilms. We showed that the increased biotransformation potential downstream of wastewater treatment plants that we observed for specific micropollutants contained in household wastewaters (downstream effect) is caused by microorganisms released with the treated effluent, rather than by the in-stream exposure to those micropollutants. Complementary data from 16S rRNA amplicon-sequencing revealed 146 amplicon sequence variants (ASVs) that followed the observed biotransformation patterns. Our results align with findings for community tolerance, and provide clear experimental evidence that microorganisms released with treated wastewater integrate into downstream biofilms and impact crucial ecosystem functions.

Environmental diagnosis with Raman Spectroscopy applied to diatoms

Freshwater quality has been changing due to the ever greater use of water resources and the contamination load resulting from human activities. Management of these systems, thus, requires constant diagnose of water quality with fast and efficient methodologies. The conventional methods adopted are, however, time-consuming, often very expensive, and require specialised expertise. Raman spectroscopy (RS) is a simple, fast and label-free technique that can be applied to environmental diagnosis using diatoms. Here, we developed a diagnostic method based on Raman spectroscopy applied to freshwater diatoms. For this, Raman spectra were recorded from diatoms of three lakes of a natural city park. The data acquired was analysed by chemometrics methods to describe the data (Partial Least Squares Regression), infer relationships in the dataset (Cluster Analysis) and produce classification models (Artificial Neural Network). The classification models developed diagnosed the lakes with excellent accuracy (89%) without requiring taxonomic information about the diatom species recorded. This study provides a proof-of-concept for the application of diatom Raman spectroscopy to diagnosing water quality, laying an important foundation for future environmental studies aiming at assessing freshwater systems, to be replicated at larger scales and to varied geographic settings.

Wastewater constituents impact biofilm microbial community in receiving streams

Microbial life in natural biofilms is dominated by prokaryotes and microscopic eukaryotes living in dense association. In stream ecosystems, microbial biofilms influence primary production, elemental cycles, food web interactions as well as water quality. Understanding how biofilm communities respond to anthropogenic impacts, such as wastewater treatment plant (WWTP) effluent, is important given the key role of biofilms in stream ecosystem function. Here, we implemented 16S and 18S rRNA gene sequencing of stream biofilms upstream (US) and downstream (DS) of WWTP effluents in four Swiss streams to test how bacterial and eukaryotic communities respond to wastewater constituents. Stream biofilm composition was strongly affected by geographic location - particularly for bacteria. However, the abundance of certain microbial community members was related to micropollutants in the wastewater - among bacteria, micropollutant-associated members were found e.g. in Alphaproteobacteria, and among eukaryotes e.g. in Bacillariophyta (algal diatoms). This study corroborates several previously characterized responses (e.g. as seen in diatoms), but also reveals previously unknown community responses - such as seen in Alphaproteobacteria. This study advances our understanding of the ecological impact of the current wastewater treatment practices and provides information about potential new marker organisms to assess ecological change in stream biofilms.

Finding needles in a haystack-Extensive diversity in the eustigmatophyceae revealed by community metabarcode analysis targeting the rbcL gene using lineage-directed primers

Sequences from the Stramenopile class Eustigmatophyceae are rarely reported in metabarcoding studies, and when they have been reported, there are very few haplotypes. We hypothesized that the paucity of eustigmatophyte species detected in these studies may be a result of the metabarcoding techniques used, which have primarily employed universal ribosomal RNA gene regions. In this study, we examined environmental DNA samples from 22 sites in southwestern Virginia, some of which had previously been studied using ribosomal RNA analysis. We used metabarcoding techniques targeting the plastid rbcL gene with new primers designed to produce a 370 bp amplicon from all lineages of the Eustigmatophyceae in a reference collection. The amplicons were then analyzed with DADA2 to produce amplicon sequence variants (ASVs). Our results revealed 184 rbcL haplotypes that can be tentatively assigned to the Eustigmatophyceae from these sites, representing much higher diversity than has been detected by ribosomal DNA-based studies. The techniques employed can be used for future studies of population structure, ecology, distribution, and diversity of this class. With these techniques, it should be possible to make realistic estimates of the species-level diversity of the Eustigmatophyceae on local, regional, and perhaps global scales.

Impact of wastewater on the microbial diversity of periphyton and its tolerance to micropollutants in an engineered flow-through channel system

Wastewater treatment plants (WWTPs) play an important role in retaining organic matter and nutrients but to a lesser extent micropollutants. Therefore, treated wastewater is recognized as a major source of multiple stressors, including complex mixtures of micropollutants. These can potentially affect microbial communities in the receiving water bodies and the ecological functions they provide. In this study, we evaluated in flow-through channels the consequences of an exposure to a mixture of stream water and different percentages of urban WWTP effluent, ranging from 0% to 80%, on the microbial diversity and function of periphyton communities. Assuming that micropollutants exert a selective pressure for tolerant microorganisms within communities, we further examined the periphyton sensitivity to a micropollutant mixture extracted from passive samplers that were immersed in the wastewater effluent. As well, micropollutants in water and in periphyton were comprehensively quantified. Our results show that micropollutants detected in periphyton differed from those found in water, both in term of concentration and composition. Especially photosystem II inhibitors accumulated in periphyton more than other pesticides. Although effects of other substances cannot be excluded, this accumulation may have contributed to the observed higher tolerance of phototrophic communities to micropollutants upon exposure to 30% and 80% of wastewater. On the contrary, no difference in tolerance was observed for heterotrophic communities. Exposure to the gradient of wastewater led to structural differences in both prokaryotic and eukaryotic communities. For instance, the relative abundance of cyanobacteria was higher with increasing percentage of wastewater effluent, whereas the opposite was observed for diatoms. Such results could indicate that differences in community structure do not necessarily lead to higher tolerance. This highlights the need to consider other wastewater constituents such as nutrients and wastewater-derived microorganisms that can modulate community structure and tolerance. By using engineered flow-through channels that mimic to some extent the required field conditions for the development of tolerance in periphyton, our study constitutes a base to investigate the mechanisms underlying the increased tolerance, such as the potential role of microorganisms originating from wastewater effluents, and different treatment options to reduce the micropollutant load in effluents.

Dumpster diving for diatom plastid 16S rRNA genes

High throughput sequencing is improving the efficiency of monitoring diatoms, which inhabit and support aquatic ecosystems across the globe. In this study, we explored the potential of a standard V4 515F-806RB primer pair in recovering diatom plastid 16S rRNA sequences. We used PhytoREF to classify the 16S reads from our freshwater biofilm field sampling from three stream segments across two streams in south-eastern Australia and retrieved diatom community data from other, publicly deposited, Australian 16S amplicon datasets. When these diatom operational taxonomic units (OTUs) were traced using the default RDPII and NCBI databases, 68% were characterized as uncultured cyanobacteria. We analysed the 16S rRNA sequences from 72 stream biofilm samples, separated the chloroplast OTUs, and classified them using the PhytoREF database. After filtering the reads attributed to Bacillariophyta (relative abundance >1%), 71 diatom OTUs comprising more than 90% of the diatom reads in each stream biofilm sample were identified. Beta-diversity analyses demonstrated significantly different diatom assemblages and discrimination among river segments. To further test the approach, the diatom OTUs from our biofilm sampling were used as reference sequences to identify diatom reads from other Australian 16S rRNA datasets in the NCBI-SRA database. Across the three selected public datasets, 67 of our 71 diatom OTUs were detected in other Australian ecosystems. Our results show that diatom plastid 16S rRNA genes are readily amplified with existing 515F-806RB primer sets. Therefore, the volume of existing 16S rRNA amplicon datasets initially generated for microbial community profiling can also be used to detect, characterize, and map diatom distribution to inform phylogeny and ecological health assessments, and can be extended into a range of ecological and industrial applications. To our knowledge, this study represents the first attempt to classify freshwater samples using this approach and the first application of PhytoREF in Australia.

A pilot experiment to assess the efficiency of pharmaceutical plant wastewater treatment and the decreasing effluent toxicity to periphytic biofilms

Pharmaceutical industry effluents are complex and highly variable in time. Assessing the efficiency of a pharmaceutical industry wastewater treatment plant (WWTP) and the resulting decrease in effluent toxicity and ecological risk is thus not straightforward. We set up an original in situ pilot directly connected to a pharmaceutical WWTP to monitor the chronic toxicity of successive effluents using natural periphytic biofilms. Their structural and functional responses to effluent exposure were assessed by combining (i) a molecular approach to characterize the bacterial and diatom diversity and (ii) functional measurements of photosynthetic and enzyme activities. Effluent contamination by pharmaceuticals strongly decreased after the quaternary treatment (activated carbon). Most of the structural biological characteristics improved with cumulative WWTP treatment (bacterial diversity, microbial genetic structure, and biological diatom index), showing community recovery along the treatment process. However, functional parameters did not show clear links with treatment steps, suggesting that microbial activities were not solely driven by pharmaceuticals produced during the experimental period. Operationally, this type of pilot system offers a useful tool for biomonitoring approaches and offers new approaches for industrial managers to assess the ecological risk of production effluents in receiving water.

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.

Survivability of Alcaligenaceae and Chromatiaceae as palm oil mill effluent pollution bioindicators under fluctuations of temperature, pH and total suspended solid

Alcaligenaceae and Chromatiaceae were previously reported as the specific pollution bioindicators in the receiving river water contaminated by palm oil mill effluent (POME) final discharge. Considering the inevitable sensitivity of bacteria under environmental stresses, it is crucial to assess the survivability of both bacteria in the fluctuated environmental factors, proving their credibility as POME pollution bioindicators in the environment. In this study, the survivability of Alcaligenaceae and Chromatiaceae from facultative pond, algae (aerobic) pond and final discharge were evaluated under varying sets of temperature (25-40°C), pH (pH 7-9) and low/high total suspended solid (TSS) contents of POME collected during low/high crop seasons of oil palm, respectively. Following treatment, the viability status and compositions of the bacterial community were assessed using flow cytometry-based assay and high-throughput Illumina MiSeq, respectively, in correlation with the changes of physicochemical properties. The changes in temperature, pH and TSS indeed changed the physicochemical properties of POME. The functionality of bacterial cells was also shifted where the viable cells and high nucleic acid contents reduced at elevated levels of temperature and pH but increased at high TSS content. Interestingly, the Alcaligenaceae and Chromatiaceae continuously detected in the samples which accounted for more than 0.5% of relative abundance, with a positive correlation with biological oxygen demand (BOD₅) concentration. Therefore, either Alcaligenaceae or Chromatiaceae or both could be regarded as the reliable and specific bacterial indicators to indicate the pollution in river water due to POME final discharge despite the fluctuations in temperature, pH and TSS.

Lake Ecosystem Robustness and Resilience Inferred from a Climate-Stressed Protistan Plankton Network

Network analyses of biological communities allow for identifying potential consequences of climate change on the resilience of ecosystems and their robustness to resist stressors. Using DNA metabarcoding datasets from a three-year-sampling (73 samples), we constructed the protistan plankton co-occurrence network of Lake Zurich, a model lake ecosystem subjected to climate change. Despite several documentations of dramatic lake warming in Lake Zurich, our study provides an unprecedented perspective by linking changes in biotic association patterns to climate stress. Water temperature belonged to the strongest environmental parameters splitting the data into two distinct seasonal networks (October–April; May–September). The expected ecological niche of phytoplankton, weakened through nutrient depletion because of permanent thermal stratification and through parasitic fungi, was occupied by the cyanobacterium Planktothrix rubescens and mixotrophic nanoflagellates. Instead of phytoplankton, bacteria and nanoflagellates were the main prey organisms associated with key predators (ciliates), which contrasts traditional views of biological associations in lake plankton. In a species extinction scenario, the warm season network emerged as more vulnerable than the cold season network, indicating a time-lagged effect of warmer winter temperatures on the communities. We conclude that climate stressors compromise lake ecosystem robustness and resilience through species replacement, richness differences, and succession as indicated by key network properties.

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.

Freshwater diatom biomonitoring through benthic kick-net metabarcoding

Biomonitoring is an essential tool for assessing ecological conditions and informing management strategies. The application of DNA metabarcoding and high throughput sequencing has improved data quantity and resolution for biomonitoring of taxa such as macroinvertebrates, yet, there remains the need to optimise these methods for other taxonomic groups. Diatoms have a longstanding history in freshwater biomonitoring as bioindicators of water quality status. However, multi-substrate periphyton collection, a common diatom sampling practice, is time-consuming and thus costly in terms of labour. This study examined whether the benthic kick-net technique used for macroinvertebrate biomonitoring could be applied to bulk-sample diatoms for metabarcoding. To test this approach, we collected samples using both conventional multi-substrate microhabitat periphyton collections and bulk-tissue kick-net methodologies in parallel from replicated sites with different habitat status (good/fair). We found there was no significant difference in community assemblages between conventional periphyton collection and kick-net methodologies or site status, but there was significant difference between diatom communities depending on site (P = 0.042). These results show the diatom taxonomic coverage achieved through DNA metabarcoding of kick-net is suitable for ecological biomonitoring applications. The shift to a more robust sampling approach and capturing diatoms as well as macroinvertebrates in a single sampling event has the potential to significantly improve efficiency of biomonitoring programmes that currently only use the kick-net technique to sample macroinvertebrates.

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.

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: R² ≈ 0.7, interception ≈ 0.8) and sensitive to global anthropogenic disturbance (Rs² > 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.

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

Ecological status assessment of watercourses is based on the calculation of quality indices using pollution sensitivity of targeted biological groups, including diatoms. The determination and quantification of diatom species is generally based on microscopic morphological identification, which requires expertise and is time-consuming and costly. In Europe, this morphological approach is legally imposed by standards and regulatory decrees by the Water Framework Directive (WFD). Over the past decade, a DNA-based molecular biology approach has newly been developed to identify species based on genetic criteria rather than morphological ones (i.e. DNA metabarcoding). In combination with high throughput sequencing technologies, metabarcoding makes it possible both to identify all species present in an environmental sample and to process several hundred samples in parallel. This article presents the results of two recent studies carried out on the WFD networks of rivers of Mayotte (2013–2018) and metropolitan France (2016–2018). These studies aimed at testing the potential application of metabarcoding for biomonitoring in the context of the WFD. We discuss the various methodological developments and optimisations that have been made to make the taxonomic inventories of diatoms produced by metabarcoding more reliable, particularly in terms of species quantification. We present the results of the application of this DNA approach on more than 500 river sites, comparing them with those obtained using the standardised morphological method. Finally, we discuss the potential of metabarcoding for routine application, its limits of application and propose some recommendations for future implementation in WFD.