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.

Physico-chemical and high frequency monitoring dataset from mesocosm experiments simulating extreme climate events in lakes

We present two datasets composed of high frequency sensors data, vertical in situ profiles and laboratory chemical analysis data, acquired during two different aquatic mesocosm experiments performed at the OLA ("Long-term observation and experimentation for lake ecosystems") facility at the UMR CARRTEL in Thonon les Bains, on the French shore of Lake Geneva. The DOMLAC experiment lasted 3 weeks (4-21 October 2021) and aimed to simulate predicted climate scenarios (i.e. extreme events such as storms and floods) by reproducing changes in quality and composition of lake subsidies and runoff by increased inputs of terrestrial organic matter. The PARLAC experiment lasted 3 weeks (5-23 September 2022) and aimed to simulate turbid storms by light reduction. The experimental setup consisted of nine inland polyester laminated tanks (2.1 m length, 2.1 m width and 1.1 m depth) with a total volume of approximately 4000 L and filled with water directly supplied from the lake at 4m depth. Both experimental design included three treatments each replicated three times. The DOMLAC experiment involved a control treatment (no treatment applied) and two treatments simulating allochthonous inputs from two different dissolved organic matter (DOM) extract from peat moss Sphagnum sp. (Peat-Moss treatment) and Phragmites australis (Phragmite treatment). The PARLAC experiment involved a control treatment (no treatment applied) and two treatments simulating two different intensity of light reduction. In the Medium treatment transmitted light was reduced to 70% and in the High treatment transmitted light was reduced to 15%. The datasets are composed of: 1. In situ measures from automated data loggers of temperature, conductivity, dissolved oxygen and CO₂ acquired every 5 minutes at 0.1, 0.5 and 1 m depth (DOMLAC) and 0.5m (PARLAC) for the entire period of the experiment. 2. In situ profiles (0-1 m) of temperature, conductivity, pH, dissolved oxygen (concentration and saturation) acquired twice a week during the experiment. 3. In situ measures of light spectral UV/VIS/IR irradiance (300-950 nm wavelength range) taken in the air and at 0, 0.5 and 1 m twice a week on the same day of the profiles at point 2. 4. Laboratory chemical analysis of integrated samples taken twice a week on the same day of the in situ profiles at point 2 and 3 of conductivity, pH, total alkalinity, NO₃, total and particulate nitrogen (Ntot, Npart), PO₄, total and particulate phosphorus (Ptot, Ppart), total and particulate organic carbon (TOC, POC), Ca, K, Mg, Na, Cl, SO₄ and SiO₂. Only for DOMLAC also analyses of NH₄, NO₂ and dissolved organic carbon (DOC). 5. Laboratory analysis of pigments (Chla, Chlc, carotenoids, phaeopigments) extracted from samples collected at point 4. 6. Only for DOMLAC, specific absorbance on the range 600-200nm of DOM (i.e. <0.7 µm) measured on samples collected at point 4. This dataset aims to contribute our understanding of how extreme climate events can alter lake subsidies and affect the regulation of ecosystem processes such as production, respiration, nutrient uptake and pigment composition. The data can be used for a wide range of applications as being included in meta-analysis aiming at generalising the effect of climate change on large lakes including simulating future scenarios in a broad range of geographical areas as we used different inputs of DOM leached from litters reproducing catchments characteristics typical of different latitudes, such as mostly dominated by large leaf forests and phragmites at middle latitude, and coniferous forests rich of peat mosses that spread along the water surface typical of Northern regions.

Short-Term Dynamics of Bdellovibrio and Like Organisms in Lake Geneva in Response to a Simulated Climatic Extreme Event

The short time-scale dynamics of three families of Bdellovibrio and like organisms (i.e. Bdellovibrionaceae, Peredibacteraceae, and Bacteriovoracaceae) were studied on the surface waters of Lake Geneva in summer. Using mesocosms deployed nearshore in July 2019, we simulated an extreme climatic event (an input of carbon from the watershed in response to runoff from the catchment, light reduction, and mixing in response to stormy conditions) and aimed to study the impact of both abiotic and biotic factors on their dynamics. The three families of Bdellovibrio and like organisms (BALOs) showed different dynamics during the experiment. Peredibacteraceae was the most abundant group, whereas Bacteriovoracaceae was the least abundant. Compared with the other two families, the abundance of Bdellovibrionaceae did not fluctuate, remaining relatively stable over time. Environmental variables only partially explained the dynamics of these families; in particular, temperature, pH, and chloride concentrations were positively correlated with Bacteriovoracaceae, Bdellovibrionaceae, and Peredibacteraceae abundance, respectively. Prokaryote-like particles (PLPs), such as those with high DNA content (HDNA), were strongly and positively correlated with Peredibacteraceae and Bacteriovoracaceae. In contrast, no relationships were found between Bdellovibrionaceae and PLP abundance, nor between the virus-like particles (VLPs) and the different BALOs. Overall, the experiment revealed that predation was stable in the face of the simulated climatic events. In addition, we observed that Peredibacteraceae and Bacteriovoracaceae share common traits, while Bdellovibrionaceae seems to constitute a distinct category.

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.

Spatio-temporal variability of eDNA signal and its implication for fish monitoring in lakes

Environmental DNA (eDNA) metabarcoding is revolutionizing the monitoring of aquatic biodiversity. The use of eDNA has the potential to enable non-invasive, cost-effective, time-efficient and high-sensitivity monitoring of fish assemblages. Although the capacity of eDNA metabarcoding to describe fish assemblages is recognised, research efforts are still needed to better assess the spatial and temporal variability of the eDNA signal and to ultimately design an optimal sampling strategy for eDNA monitoring. In this context, we sampled three different lakes (a dam reservoir, a shallow eutrophic lake and a deep oligotrophic lake) every 6 weeks for 1 year. We performed four types of sampling for each lake (integrative sampling of sub-surface water along transects on the left shore, the right shore and above the deepest zone, and point sampling in deeper layers near the lake bottom) to explore the spatial variability of the eDNA signal at the lake scale over a period of 1 year. A metabarcoding approach was applied to analyse the 92 eDNA samples in order to obtain fish species inventories which were compared with traditional fish monitoring methods (standardized gillnet samplings). Several species known to be present in these lakes were only detected by eDNA, confirming the higher sensitivity of this technique in comparison with gillnetting. The eDNA signal varied spatially, with shoreline samples being richer in species than the other samples. Furthermore, deep-water samplings appeared to be non-relevant for regularly mixed lakes, where the eDNA signal was homogeneously distributed. These results also demonstrate a clear temporal variability of the eDNA signal that seems to be related to species phenology, with most of the species detected in spring during the spawning period on shores, but also a peak of detection in winter for salmonid and coregonid species during their reproduction period. These results contribute to our understanding of the spatio-temporal distribution of eDNA in lakes and allow us to provide methodological recommendations regarding where and when to sample eDNA for fish monitoring in lakes.

Paleoreconstructions of ciliate communities reveal long-term ecological changes in temperate lakes

Ciliates are unicellular heterotrophic organisms that play a key role in aquatic planktonic and benthic food webs. Advances in sedimentary DNA (sed-DNA) analysis offer the possibility to integrate these bioindicators in paleoenvironmental reconstructions. In this study, we used the top–bottom paleolimnological approach and metabarcoding techniques applied to sed-DNA to compare the recent and past (i.e. prior to major anthropogenic impacts) ciliate communities of 48 lakes located along an elevation gradient. Our results show an overall decline in the β-diversity in recent time, especially in lowland lakes, which are more strongly exposed to local human pressures. Analyses of the functional groups indicate important restructuration of the food web, including the recent increase in mixotrophs. Moreover, changes in the benthic ciliates were consistent with the widespread increase in deep water anoxia. Our results provided evidence that sed-DNA can uncover information about past ciliate communities on a wide variety of lakes. Overall, our study demonstrates the potential of using ciliates as new paleoindicators, integrating information from the pelagic to the benthic zones, and providing valuable insights into ecosystem functioning through a trait-based functional community approach. As paleoindicator, they thus offer a more holistic view on the long-term changes of aquatic ecosystems.

Environmental paleomicrobiology: using DNA preserved in aquatic sediments to its full potential

In‐depth knowledge about spatial and temporal variation in microbial diversity and function is needed for a better understanding of ecological and evolutionary responses to global change. In particular, the study of microbial ancient DNA preserved in sediment archives from lakes and oceans can help us to evaluate the responses of aquatic microbes in the past and make predictions about future biodiversity change in those ecosystems. Recent advances in molecular genetic methods applied to the analysis of historically deposited DNA in sediments have not only allowed the taxonomic identification of past aquatic microbial communities but also enabled tracing their evolution and adaptation to episodic disturbances and gradual environmental change. Nevertheless, some challenges remain for scientists to take full advantage of the rapidly developing field of paleo‐genetics, including the limited ability to detect rare taxa and reconstruct complete genomes for evolutionary studies. Here, we provide a brief review of some of the recent advances in the field of environmental paleomicrobiology and discuss remaining challenges related to the application of molecular genetic methods to study microbial diversity, ecology, and evolution in sediment archives. We anticipate that, in the near future, environmental paleomicrobiology will shed new light on the processes of microbial genome evolution and microbial ecosystem responses to quaternary environmental changes at an unprecedented level of detail. This information can, for example, aid geological reconstructions of biogeochemical cycles and predict ecosystem responses to environmental perturbations, including in the context of human‐induced global changes.

Phytool, a ShinyApp to homogenise taxonomy of freshwater microalgae from DNA barcodes and microscopic observations

Methods for biomonitoring of freshwater phytoplankton are evolving rapidly with eDNA-based methods, offering great complementarity with microscopy. Metabarcoding approaches have been more commonly used over the last years, with a continuous increase in the amount of data generated. Depending on the researchers and the way they assigned barcodes to species (bioinformatic pipelines and molecular reference databases), the taxonomic assignment obtained for HTS DNA reads might vary. This is also true for traditional taxonomic studies by microscopy with regular adjustments of the classification and taxonomy.

For those reasons (leading to non-homogeneous taxonomies), gap-analyses and comparisons between studies become even more challenging and the curation processes to find potential consensus names are time-consuming. Here, we present a web-based application (Phytool), developed with ShinyApp (Rstudio), that aims to make the harmonisation of taxonomy easier and in a more efficient way, using a complete and up-to-date taxonomy reference database for freshwater microalgae. Phytool allows users to homogenise and update freshwater phytoplankton taxonomical names from sequence files and data tables directly uploaded in the application. It also gathers barcodes from curated references in a user-friendly way in which it is possible to search for specific organisms. All the data provided are downloadable with the possibility to apply filters in order to select only the required taxa and fields (e.g. specific taxonomic ranks). The main goal is to make accessible to a broad range of users the connection between microscopy and molecular biology and taxonomy through different ready-to-use functions. This study estimates that only 25% of species of freshwater phytoplankton in Phytobs are associated with a barcode. We plead for an increased effort to enrich reference databases by coupling taxonomy and molecular methods. Phytool should make this crucial work more efficient.

The application is available at https://caninuzzo.shinyapps.io/phytool_v1/

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.

Seasonal Dynamics of Abundance, Structure, and Diversity of Methanogens and Methanotrophs in Lake Sediments

Understanding temporal and spatial microbial community abundance and diversity variations is necessary to assess the functional roles played by microbial actors in the environment. In this study, we investigated spatial variability and temporal dynamics of two functional microbial sediment communities, methanogenic Archaea and methanotrophic bacteria, in Lake Bourget, France. Microbial communities were studied from 3 sites sampled 4 times over a year, with one core sampled at each site and date, and 5 sediment layers per core were considered. Microbial abundance in the sediment were determined using flow cytometry. Methanogens and methanotrophs community structures, diversity, and abundance were assessed using T-RFLP, sequencing, and real-time PCR targeting mcrA and pmoA genes, respectively. Changes both in structure and abundance were detected mainly at the water-sediment interface in relation to the lake seasonal oxygenation dynamics. Methanogen diversity was dominated by Methanomicrobiales (mainly Methanoregula) members, followed by Methanosarcinales and Methanobacteriales. For methanotrophs, diversity was dominated by Methylobacter in the deeper area and by Methylococcus in the shallow area. Organic matter appeared to be the main environmental parameter controlling methanogens, while oxygen availability influenced both the structure and abundance of the methanotrophic community.

Physico-chemical dataset from an in situ mesocosm experiment simulating extreme climate events in Lake Geneva (MESOLAC)

This dataset complement a previously published dataset [1] and corresponds to the physico-chemical parameters data series produced during the MESOLAC experimental project [2]. The presented dataset is composed of: 1. In situ profiles (0-3m) of temperature, conductivity, pH, dissolved oxygen (concentration and saturation). 2. In situ measurements of light spectral UV/VIS/IR irradiance (300-950 nm wavelength range) taken at 0, 0.25, 0.5, 1, 1.5, 2 and 2.5m. 3. Laboratory chemical analysis of samples collected at 0 and 2 m (conductivity, pH, total alkalinity, NH4, NO2, NO3, total and particulate nitrogen (Ntot, Npart), PO4, total and particulate phosphorus (Ptot, Ppart), total, organic particulate and total particulate carbon (Ctot, Cpart-org, Cpart-tot), Cl, SO4, SiO2. 4. Laboratory analysis of pigments extracted from samples collected at 0 and 2 m (Chla, Chlc, carotenoids, phaeopigments). The experimental design is the same as in Tran-Khac et al [1]. Briefly, it consisted of nine pelagic mesocosms (about 3000 L, 3m depth) deployed in July 2019 in Lake Geneva near the shore of Thonon les Bains (France) aiming to simulate predicted climate scenarios (i.e. extreme events) and assess the response of planktonic communities, ecosystem functioning and resilience. During the experiment, physical parameters were measured twice a week. At the same time, samples were collected at 0 and 2m of depth for subsequent chemical laboratory analyses. These data are presented in the dataset file, ordered by sampling event (numbered from S1 to S8), treatment (Control-C, High-H and Medium-M) and replicates (1 to 3). For each sampling point the measured parameters are listed in columns, missing data and values below the detection limit are marked as NA (not available). This data set aims to contribute to the understanding of the effect of environmental forcing on lake physico-chemical characteristics (such as temperature, oxygen and nutrient concentration) under simulated intense weather events. To a broader extent, the presented data can be used for a wide variety of applications, including monitoring of a large peri-alpine lake functioning under environmental stress and being included in further meta-analysis to generalise the effect of climate change on large lakes. The two complementary dataset differ in the acquired data and methods, temporal and spatial resolution. They complete each other in terms of physico-chemical characterization of the experimental treatments and together can allow comparison of the two different monitoring strategies (continuous vs punctual) during in situ experimental manipulations.

Landscape Setting Drives the Microbial Eukaryotic Community Structure in Four Swedish Mountain Lakes over the Holocene

On the annual and interannual scales, lake microbial communities are known to be heavily influenced by environmental conditions both in the lake and in its terrestrial surroundings. However, the influence of landscape setting and environmental change on shaping these communities over a longer (millennial) timescale is rarely studied. Here, we applied an 18S metabarcoding approach to DNA preserved in Holocene sediment records from two pairs of co-located Swedish mountain lakes. Our data revealed that the microbial eukaryotic communities were strongly influenced by catchment characteristics rather than location. More precisely, the microbial communities from the two bedrock lakes were largely dominated by unclassified Alveolata, while the peatland lakes showed a more diverse microbial community, with Ciliophora, Chlorophyta and Chytrids among the more predominant groups. Furthermore, for the two bedrock-dominated lakes-where the oldest DNA samples are dated to only a few hundred years after the lake formation-certain Alveolata, Chlorophytes, Stramenopiles and Rhizaria taxa were found prevalent throughout all the sediment profiles. Our work highlights the importance of species sorting due to landscape setting and the persistence of microbial eukaryotic diversity over millennial timescales in shaping modern lake microbial communities.

Protist taxonomic and functional diversity in soil, freshwater and marine ecosystems

Protists dominate eukaryotic diversity and play key functional roles in all ecosystems, particularly by catalyzing carbon and nutrient cycling. To date, however, a comparative analysis of their taxonomic and functional diversity that compares the major ecosystems on Earth (soil, freshwater and marine systems) is missing. Here, we present a comparison of protist diversity based on standardized high throughput 18S rRNA gene sequencing of soil, freshwater and marine environmental DNA. Soil and freshwater protist communities were more similar to each other than to marine protist communities, with virtually no overlap of Operational Taxonomic Units (OTUs) between terrestrial and marine habitats. Soil protists showed higher γ diversity than aquatic samples. Differences in taxonomic composition of the communities led to changes in a functional diversity among ecosystems, as expressed in relative abundance of consumers, phototrophs and parasites. Phototrophs (eukaryotic algae) dominated freshwater systems (49% of the sequences) and consumers soil and marine ecosystems (59% and 48%, respectively). The individual functional groups were composed of ecosystem- specific taxonomic groups. Parasites were equally common in all ecosystems, yet, terrestrial systems hosted more OTUs assigned to parasites of macro-organisms while aquatic systems contained mostly microbial parasitoids. Together, we show biogeographic patterns of protist diversity across major ecosystems on Earth, preparing the way for more focused studies that will help understanding the multiple roles of protists in the biosphere.

In situ pelagic dataset from continuous monitoring: A mesocosm experiment in Lake Geneva (MESOLAC)

This dataset corresponds to a data series produced from automated data loggers during the MESOLAC experimental project. Nine pelagic mesocosms (about 3000 L, 3 m depth) were deployed in July 2019 in Lake Geneva near the shore of Thonon les Bains (France), simulating predicted climate scenarios (i.e. intense weather events) by applying a combination of forcing. The design consisted of three treatments each replicated three times: a control treatment (named C - no treatment applied) and two different treatments simulating different intensities of weather events. The high intensity treatment (named H) aimed to reproduce short and intense weather events such as violent storms. It consisted of a short-term stress applied during the first week, with high pulse of dissolved organic carbon (5x increased concentration, i.e. total DOC ∼ 6 mg L-1), transmitted light reduced to 15% and water column manual mixing. The medium intensity treatment (named M) simulated less intense and more prolonged exposures such as during flood events. It was maintained during the 4 weeks of the experiment and consisted of 1.5x increased concentration of dissolved organic carbon (i.e. total DOC ∼ 2 mg L-1), 70% transmitted light and water column manual mixing. Automated data loggers were placed for the entire period of the experiment in the mesocosms and in the lake for comparison with natural conditions. Temperature, conductivity, dissolved oxygen and CO2 were monitored every 15 min at different depths (0.15, 0.25, 1 and 2 m). This data set aims to contribute our understanding of the effect of environmental forcing on lake ecosystem processes (such as production, respiration and CO2 exchange) under simulated intense weather events and the ability of the planktonic community to recover after perturbation. To a broader extent, the presented data can be used for a wide variety of applications, including monitoring of lake community functioning during a period of high productivity on a large peri-alpine lake and being included in further meta-analysis aiming at generalising the effect of climate change on large lakes.

Assessing the response of micro-eukaryotic diversity to the Great Acceleration using lake sedimentary DNA

Long-term time series have provided evidence that anthropogenic pressures can threaten lakes. Yet it remains unclear how and the extent to which lake biodiversity has changed during the Anthropocene, in particular for microbes. Here, we used DNA preserved in sediments to compare modern micro-eukaryotic communities with those from the end of the 19th century, i.e., before acceleration of the human imprint on ecosystems. Our results obtained for 48 lakes indicate drastic changes in the composition of microbial communities, coupled with a homogenization of their diversity between lakes. Remote high elevation lakes were globally less impacted than lowland lakes affected by local human activity. All functional groups (micro-algae, parasites, saprotrophs and consumers) underwent significant changes in diversity. However, we show that the effects of anthropogenic changes have benefited in particular phototrophic and mixotrophic species, which is consistent with the hypothesis of a global increase of primary productivity in lakes.

Global distribution of Trebouxiophyceae diversity explored by high-throughput sequencing and phylogenetic approaches

Trebouxiophyceae are a ubiquitous class of Chlorophyta encountered in aquatic and terrestrial environments. Most taxa are photosynthetic, and many acts as photobionts in symbiotic relationships, while others are free-living. Trebouxiophyceae have also been widely investigated for their use for biotechnological applications. In this work, we aimed at obtaining a comprehensive image of their diversity by compiling the information of 435 freshwater, soil and marine environmental DNA samples surveyed with Illumina sequencing technology in order to search for the most relevant environments for bioprospecting. Freshwater and soil were most diverse and shared more than half of all operational taxonomic units (OTUs), however, their communities were significantly distinct. Oceans hosted the highest genetic novelty, and did not share any OTUs with the other environments; also, marine samples host more diversity in warm waters. Symbiotic genera usually found in lichens such as Trebouxia, Myrmecia and Symbiochloris were also abundantly detected in the ocean, suggesting either free-living lifestyles or unknown symbiotic relationships with marine planktonic organisms. Altogether, our study opens the way to new prospection for trebouxiophycean strains, especially in understudied environments like the ocean.

Diversity, spatial distribution and activity of fungi in freshwater ecosystems

High-throughput sequencing has given new insights into aquatic fungal community ecology over the last 10 years. Based on 18S ribosomal RNA gene sequences publicly available, we investigated fungal richness and taxonomic composition among 25 lakes and four rivers. We used a single pipeline to process the reads from raw data to the taxonomic affiliation. In addition, we studied, for a subset of lakes, the active fraction of fungi through the 18S rRNA transcripts level. These results revealed a high diversity of fungi that can be captured by 18S rRNA primers. The most OTU-rich groups were Dikarya (47%), represented by putative filamentous fungi more diverse and abundant in freshwater habitats than previous studies have suggested, followed by Cryptomycota (17.6%) and Chytridiomycota (15.4%). The active fraction of the community showed the same dominant groups as those observed at the 18S rRNA genes level. On average 13.25% of the fungal OTUs were active. The small number of OTUs shared among aquatic ecosystems may result from the low abundances of those microorganisms and/or they constitute allochthonous fungi coming from other habitats (e.g., sediment or catchment areas). The richness estimates suggest that fungi have been overlooked and undersampled in freshwater ecosystems, especially rivers, though they play key roles in ecosystem functioning as saprophytes and parasites.

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

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

The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems

The bioassessment of aquatic ecosystems is currently based on various biotic indices that use the occurrence and/or abundance of selected taxonomic groups to define ecological status. These conventional indices have some limitations, often related to difficulties in morphological identification of bioindicator taxa. Recent development of DNA barcoding and metabarcoding could potentially alleviate some of these limitations, by using DNA sequences instead of morphology to identify organisms and to characterize a given ecosystem. In this paper, we review the structure of conventional biotic indices, and we present the results of pilot metabarcoding studies using environmental DNA to infer biotic indices. We discuss the main advantages and pitfalls of metabarcoding approaches to assess parameters such as richness, abundance, taxonomic composition and species ecological values, to be used for calculation of biotic indices. We present some future developments to fully exploit the potential of metabarcoding data and improve the accuracy and precision of their analysis. We also propose some recommendations for the future integration of DNA metabarcoding to routine biomonitoring programs.

Trophic relationships between metazooplankton communities and their plankton food sources in the Iles Eparses (Western Indian Ocean)

Coral reef and atoll lagoons are among the most diversified marine ecosystems but also the most affected by the combined effects of climate change and human activities. The Iles Eparses (Scattered Islands) in the Western Indian Ocean have been little affected by human pressure and can be considered to be "pristine" ecosystems. Metazooplankton plays a major role in the functioning and productivity of aquatic ecosystems, and this study was undertaken: (i) to determine the spatial abundance, distribution and species composition of metazooplankton, (ii) to assess the effect of metazooplankton grazing on pico- and nanophytoplankton and (iii) to analyze the trophic positions of metazooplankton by using the stable isotope signatures of a wide variety of taxa and particulate organic matter from the Iles Eparses and Mayotte. Tromelin Island (which is not located in the Mozambique Channel) had the lowest metazooplankton abundance with no cyanobacteria Trichodesmium spp. or mollusks (pteropods) presence, and with δ(15)N signatures of organisms that were higher than for the islands in the Mozambique Channel. Trichodesmium spp. was found in the Mozambique Channel and the plankton food web was probably based preferentially on these cyanobacteria with lower δ(15)N signatures indicating direct or indirect trophic transfer of diazotrophic nitrogen to metazooplankton. Three of the islands were distinct: Europa had the highest proportion of copepods, with oithonids being dominant, which is typical of rich mangrove systems, while Juan de Nova and Mayotte seemed to be the sites most affected by human activity with a high abundance of appendicularians and distinct particulate organic matter ∂(13)C signatures. Grazing experiments showed that food could be a limiting factor for metazooplankton in the Iles Eparses. However, the effect of metazooplankton grazing on phytoplankton appeared to be very low (0.01-2.32% of the total phytoplankton per day).

Diversity and Dynamics of Active Small Microbial Eukaryotes in the Anoxic Zone of a Freshwater Meromictic Lake (Pavin, France)

Microbial eukaryotes play a crucial role in ecosystem functioning and oxygen is considered to be one of the strongest barriers against their local dispersal. However, diversity of microbial eukaryotes in freshwater habitats with oxygen gradients has previously received very little attention. We applied high-throughput sequencing (V4 region of the 18S rRNA gene) in conjunction with quantitative PCR (DNA and RNA) and fluorescent in situ hybridization (FISH) analyses, to provide an unique spatio-temporal analysis of microbial eukaryotes diversity and potential activity in a meromictic freshwater lake (lake Pavin). This study revealed a high genetic diversity of unicellular eukaryotes in the permanent anoxic zone of lake Pavin and allowed the discrimination of active vs. inactive components. Forty-two percent of the OTUs (Operational Taxonomic Units) are exclusively present in the monimolimnion, where Alveolata (Ciliophora and Dinophyceae) and Fungi (Dikarya and Chytrids) are the most active phyla and are probably represented by species capable of anaerobic metabolism. Pigmented eukaryotes (Haptophyceae and Chlorophyceae) are also present and active in this zone, which opens up questions regarding their metabolism.

Is Planktonic Diversity Well Recorded in Sedimentary DNA? Toward the Reconstruction of Past Protistan Diversity

Studies based on the coupling of a paleolimnological approach and molecular tools (e.g., sequencing of sedimentary DNA) present a promising opportunity to obtain long-term data on past lacustrine biodiversity. However, certain validations are still required, such as the evaluation of DNA preservation in sediments for various planktonic taxa that do not leave any morphological diagnostic features. In this study, we focused on the diversity of planktonic unicellular eukaryotes and verified the presence of their DNA in sediment archives. We compared the molecular inventories (high-throughput sequencing of 18S ribosomal DNA) obtained from monitoring the water column with those obtained for DNA archived in the first 30 cm of sediment. Seventy-one percent of taxonomic units found in the water samples were detected in sediment samples, including pigmented taxa, such as Chlorophyta, Dinophyceae, and Chrysophyceae, phagotrophic taxa, such as Ciliophora, parasitic taxa, such as Apicomplexa and Chytridiomycota, and saprotrophs, such as Cryptomycota. Parallel analysis of 18S ribosomal RNA (rRNA) transcripts revealed the presence of living eukaryotic taxa only in the top 2 cm of sediment; although some limits exist in using RNA/DNA ratio as indicator of microbial activity, these results suggested that the sedimentary DNA mostly represented DNA from past and inactive communities. Only the diversity of a few groups, such as Cryptophyta and Haptophyta, seemed to be poorly preserved in sediments. Our overall results showed that the application of sequencing techniques to sedimentary DNA could be used to reconstruct past diversity for numerous planktonic eukaryotic groups.

Temporal Dynamics of Active Prokaryotic Nitrifiers and Archaeal Communities from River to Sea

To test if different niches for potential nitrifiers exist in estuarine systems, we assessed by pyrosequencing the diversity of archaeal gene transcript markers for taxonomy (16S ribosomal RNA (rRNA)) during an entire year along a salinity gradient in surface waters of the Charente estuary (Atlantic coast, France). We further investigated the potential for estuarine prokaryotes to oxidize ammonia and hydrolyze urea by quantifying thaumarchaeal amoA and ureC and bacterial amoA transcripts. Our results showed a succession of different nitrifiers from river to sea with bacterial amoA transcripts dominating in the freshwater station while archaeal transcripts were predominant in the marine station. The 16S rRNA sequence analysis revealed that Thaumarchaeota marine group I (MGI) were the most abundant overall but other archaeal groups like Methanosaeta were also potentially active in winter (December-March) and Euryarchaeota marine group II (MGII) were dominant in seawater in summer (April-August). Each station also contained different Thaumarchaeota MGI phylogenetic clusters, and the clusters' microdiversity was associated to specific environmental conditions suggesting the presence of ecotypes adapted to distinct ecological niches. The amoA and ureC transcript dynamics further indicated that some of the Thaumarchaeota MGI subclusters were involved in ammonia oxidation through the hydrolysis of urea. Our findings show that ammonia-oxidizing Archaea and Bacteria were adapted to contrasted conditions and that the Thaumarchaeota MGI diversity probably corresponds to distinct metabolisms or life strategies.

Temporal dynamics of active Archaea in oxygen-depleted zones of two deep lakes

Deep lakes are of specific interest in the study of archaeal assemblages as chemical stratification in the water column allows niche differentiation and distinct community structure. Active archaeal community and potential nitrifiers were investigated monthly over 1 year by pyrosequencing 16S rRNA transcripts and genes, and by quantification of archaeal amoA genes in two deep lakes. Our results showed that the active archaeal community patterns of spatial and temporal distribution were different between these lakes. The meromictic lake characterized by a stable redox gradient but variability in nutrient concentrations exhibited large temporal rearrangements of the dominant euryarchaeal phylotypes, suggesting a variety of ecological niches and dynamic archaeal communities in the hypolimnion of this lake. Conversely, Thaumarchaeota Marine Group I (MGI) largely dominated in the second lake where deeper water layers exhibited only short periods of complete anoxia and constant low ammonia concentrations. Investigations conducted on archaeal amoA transcripts abundance suggested that not all lacustrine Thaumarchaeota conduct the process of nitrification. A high number of 16S rRNA transcripts associated to crenarchaeal group C3 or the Miscellaneous Euryarchaeotic Group indicates the potential for these uncharacterized groups to contribute to nutrient cycling in lakes.

Acceleration of cyanobacterial dominance in north temperate-subarctic lakes during the Anthropocene

Increases in atmospheric temperature and nutrients from land are thought to be promoting the expansion of harmful cyanobacteria in lakes worldwide, yet to date there has been no quantitative synthesis of long-term trends. To test whether cyanobacteria have increased in abundance over the past ~ 200 years and evaluate the relative influence of potential causal mechanisms, we synthesised 108 highly resolved sedimentary time series and 18 decadal-scale monitoring records from north temperate-subarctic lakes. We demonstrate that: (1) cyanobacteria have increased significantly since c. 1800 ce, (2) they have increased disproportionately relative to other phytoplankton, and (3) cyanobacteria increased more rapidly post c. 1945 ce. Variation among lakes in the rates of increase was explained best by nutrient concentration (phosphorus and nitrogen), and temperature was of secondary importance. Although cyanobacterial biomass has declined in some managed lakes with reduced nutrient influx, the larger spatio-temporal scale of sedimentary records show continued increases in cyanobacteria throughout the north temperate-subarctic regions.

The need for ecological monitoring of freshwaters in a changing world: a case study of Lakes Annecy, Bourget, and Geneva

Lakes Annecy, Bourget, and Geneva are large, deep carbonated peri-alpine lakes in eastern France. They are located in the same ecoregion but have been subject to differing degrees of anthropogenic pressure over the past decades. A comparative analysis of these ecosystems can therefore provide valuable information on how the lakes have responded to changes in phosphorus runoff, fish management practices, and global warming. Each of these lakes has undergone a restoration process, and changes in water quality and trophic state, as measured using parameters like transparency, chlorophyll a, nutrient concentrations, and phytoplankton biomass and structure, can be used to evaluate efforts made to preserve these ecosystems. Our results reveal that (1) peri-alpine lakes are exemplary cases of restoration in the world where freshwater eutrophication is on the increase, and (2) efforts must be maintained because of the new context of climate change, the effects of which on the quality and the ecological functioning of lakes are still poorly understood.

When is sampling complete? The effects of geographical range and marker choice on perceived diversity in Nitzschia palea (Bacillariophyta)

DNA barcoding, being developed for biomonitoring, requires a database of reference sequences and knowledge of how much sequences can deviate before they are assigned to separate species. The molecular hunt for hidden species also raises the question of species definitions. We examined whether there are objective criteria for sequence-based species delimitation in diatoms, using Nitzschia palea, an important monophyletic indicator species already known to contain cryptic diversity. Strains from a wide geographical range were sequenced for 28S rRNA, COI and rbcL. Homogeneity indices and the Chao index failed to objectively select a precise number of species existing in N. palea as well as an evolutionary method based on coalescence theory. COI always gave higher diversity estimations than 28S rRNA or rbcL. Mating data did not provide a precise calibration of molecular species thresholds. Rarefaction curves indicated that further MOTUs would be detected with more isolates than we sampled (81 clones, 42 localities). Although some genotypes had intercontinental distributions, there was a positive relationship between genetic and geographical distance, suggesting even higher richness than we assessed, given that many regions were not sampled. Overall, no objective criteria were found for species separation; instead barcoding will need a consensual approach to molecular species limits.

Trophic transfer of microcystins through the lake pelagic food web: evidence for the role of zooplankton as a vector in fish contamination

An in situ study was performed to investigate the role of zooplankton as a vector of microcystins (MCs) from Planktothrix rubescens filaments to fish during a metalimnic bloom of P. rubescens in Lake Hallwil (Switzerland). The concentrations of MCs in P. rubescens and various zooplanktonic taxa (filter-feeders and predators) were assessed in different water strata (epi-, meta- and hypolimnion) using replicated sampling over a 24-hour survey. The presence of P. rubescens in the gut content of various zooplanktonic taxa (Daphnia, Bosmina and Chaoborus) was verified by targeting the cyanobacterial nucleic acids (DNA). These results highlighted that cyanobacterial cells constitute a part of food resource for herbivorous zooplanktonic taxa during metalimnic bloom periods. Furthermore, presence of MCs in Chaoborus larvae highlighted the trophic transfer of MCs between herbivorous zooplankton and their invertebrate predators. Our results suggest that zooplanktonic herbivores by diel vertical migration (DVM) act as vectors of MCs by encapsulating grazed cyanobacteria. As a consequence, they largely contribute to the contamination of zooplanktonic predators, and in fine of zooplanktivorous whitefish. Indeed, we estimated the relative contribution of three preys of the whitefish (i.e. Daphnia, Bosmina and Chaoborus) to diet contamination. We showed that Chaoborus and Daphnia were the highest contributor as MC vectors in the whitefish diet (74.6 and 20.5% of MC-LR equivalent concentrations, respectively). The transfer of MCs across the different trophic compartments follows complex trophic pathways involving various trophic levels whose relative importance in fish contamination might vary at daily and seasonal scale.

Phylogenetic affiliation of SSU rRNA genes generated by massively parallel sequencing: new insights into the freshwater protist diversity

Recent advances in next-generation sequencing (NGS) technologies spur progress in determining the microbial diversity in various ecosystems by highlighting, for example, the rare biosphere. Currently, high-throughput pyrotag sequencing of PCR-amplified SSU rRNA gene regions is mainly used to characterize bacterial and archaeal communities, and rarely to characterize protist communities. In addition, although taxonomic assessment through phylogeny is considered as the most robust approach, similarity and probabilistic approaches remain the most commonly used for taxonomic affiliation. In a first part of this work, a tree-based method was compared with different approaches of taxonomic affiliation (BLAST and RDP) of 18S rRNA gene sequences and was shown to be the most accurate for near full-length sequences and for 400 bp amplicons, with the exception of amplicons covering the V5-V6 region. Secondly, the applicability of this method was tested by running a full scale test using an original pyrosequencing dataset of 18S rRNA genes of small lacustrine protists (0.2-5 µm) from eight freshwater ecosystems. Our results revealed that i) fewer than 5% of the operational taxonomic units (OTUs) identified through clustering and phylogenetic affiliation had been previously detected in lakes, based on comparison to sequence in public databases; ii) the sequencing depth provided by the NGS coupled with a phylogenetic approach allowed to shed light on clades of freshwater protists rarely or never detected with classical molecular ecology approaches; and iii) phylogenetic methods are more robust in describing the structuring of under-studied or highly divergent populations. More precisely, new putative clades belonging to Mamiellophyceae, Foraminifera, Dictyochophyceae and Euglenida were detected. Beyond the study of protists, these results illustrate that the tree-based approach for NGS based diversity characterization allows an in-depth description of microbial communities including taxonomic profiling, community structuring and the description of clades of any microorganisms (protists, Bacteria and Archaea).

Dynamics of ammonia-oxidizing Archaea and Bacteria in contrasted freshwater ecosystems

Thaumarchaeota have been recognized as the main drivers of aerobic ammonia oxidation in many ecosystems. However, little is known about the role of ammonia-oxidizing Archaea (AOA) and Bacteria (AOB) in lacustrine ecosystems. In this study, the photic zone of three contrasted freshwater ecosystems located in France was sampled during two periods: winter homothermy (H) and summer thermal stratification (TS), to investigate the distribution of planktonic AOA and AOB. We showed that AOB were predominant in nutrient-rich ecosystems, whereas AOA dominated when ammonia concentrations were the lowest and during winter, which could provide a favorable environment for their growth. Moreover, analyses of archaeal libraries revealed the ubiquity of the thaumarchaeal I.1a clade associated with higher diversity of AOA in the most nutrient-poor lake. More generally, this work assesses the presence of AOA in lakes, but also highlights the existence of clades typically associated with lacustrine and hot spring ecosystems and specific ecological niches occupied by these microorganisms.

Impact of toxic cyanobacterial blooms on Eurasian perch (Perca fluviatilis): experimental study and in situ observations in a peri-alpine lake

Due to the importance of young-of-the-year (YOY) perch in the peri-alpine regions where they are consumed, the microcystin (MC) contamination of YOY perch was analysed both in field (Lake Bourget, France) and experimentally using force-feeding protocols with pure MCs. In-situ, schools of YOY perch present in the epilimnion of the lake were never found in direct contact with the P. rubescens blooms that were present in the metalimnion. However, MCs were detected in the muscles and liver of the fish and were thus assumed to reach YOY perch through dietary routes, particularly via the consumption of MC-containing Daphnia. Force-feeding experiment demonstrates the existence of MC detoxification/excretion processes and suggests that in situ, YOY perch could partly detoxify and excrete ingested MCs, thereby limiting the potential negative effects on perch populations under bloom conditions. However, because of chronic exposure these processes could not allow for the complete elimination of MCs. In both experimental and in situ studies, no histological change was observed in YOY perch, indicating that MC concentrations that occurred in Lake Bourget in 2009 were too low to cause histological damage prone to induce mortality. However, Deoxyribonucleic acid (DNA) damages were observed for both the high and low experimental MC doses, suggesting that similar effects could occur in situ and potentially result in perch population disturbance during cyanobacterial blooms. Our results indicate the presence of MCs in wild perch, the consumption of this species coming from Lake Bourget is not contested but more analyses are needed to quantify the risk.

Short-term uptake of microcystin-LR by Coregonus lavaretus: GST activity and genotoxicity

In the present study, juvenile whitefish weighing 2 g were exposed by force-feeding to two ecologically relevant doses (0.05 and 0.5 μg per fish) of microcystin-LR (MC-LR). Then over 96 h the MC uptake in fish liver and muscle was measured, as the activity of the detoxification enzyme glutathione S-transferase (GST) in the liver, and the genotoxicity impact on red blood cells. Results show that (1) the MC-LR equivalent concentrations increased for both doses and in both organs of whitefish with approximately threefold lower concentrations for the low dose compared to the high dose in both organs and threefold lower concentrations in the muscle compared to the liver for each dose (2) the liver GST activity increased during the first 48 h of exposure with fivefold higher GST activity for the highest dose at 48 h compared to control and (3) MC-LR leads to deoxyribonucleic acid strand breaks that were detected by the comet assay and shown to be partially repaired. This work demonstrates that European whitefish could be impacted by cyanobacteria toxins due to rapid microcystin uptake, especially in the context of chronic contamination, which can occur during long bloom episodes.

Short-term responses of unicellular planktonic eukaryotes to increases in temperature and UVB radiation

Background

Small size eukaryotes play a fundamental role in the functioning of coastal ecosystems, however, the way in which these micro-organisms respond to combined effects of water temperature, UVB radiations (UVBR) and nutrient availability is still poorly investigated.

Results

We coupled molecular tools (18S rRNA gene sequencing and fingerprinting) with microscope-based identification and counting to experimentally investigate the short-term responses of small eukaryotes (<6 μm; from a coastal Mediterranean lagoon) to a warming treatment (+3°C) and UVB radiation increases (+20%) at two different nutrient levels. Interestingly, the increase in temperature resulted in higher pigmented eukaryotes abundances and in community structure changes clearly illustrated by molecular analyses. For most of the phylogenetic groups, some rearrangements occurred at the OTUs level even when their relative proportion (microscope counting) did not change significantly. Temperature explained almost 20% of the total variance of the small eukaryote community structure (while UVB explained only 8.4%). However, complex cumulative effects were detected. Some antagonistic or non additive effects were detected between temperature and nutrients, especially for Dinophyceae and Cryptophyceae.

Conclusions

This multifactorial experiment highlights the potential impacts, over short time scales, of changing environmental factors on the structure of various functional groups like small primary producers, parasites and saprotrophs which, in response, can modify energy flow in the planktonic food webs.

Quantitative PCR enumeration of total/toxic Planktothrix rubescens and total cyanobacteria in preserved DNA isolated from lake sediments

The variability of spatial distribution and the determinism of cyanobacterial blooms, as well as their impact at the lake scale, are still not understood, partly due to the lack of long-term climatic and environmental monitoring data. The paucity of these data can be alleviated by the use of proxy data from high-resolution sampling of sediments. Coupling paleolimnological and molecular tools and using biomarkers such as preserved DNA are promising approaches, although they have not been performed often enough so far. In our study, a quantitative PCR (qPCR) technique was applied to enumerate total cyanobacterial and total and toxic Planktothrix communities in preserved DNA derived from sediments of three lakes located in the French Alps (Lake Geneva, Lake Bourget, and Lake Annecy), containing a wide range of cyanobacterial species. Preserved DNA from lake sediments was analyzed to assess its quality, quantity, and integrity, with further application for qPCR. We applied the qPCR assay to enumerate the total cyanobacterial community, and multiplex qPCR assays were applied to quantify total and microcystin-producing Planktothrix populations in a single reaction tube. These methods were optimized, calibrated, and applied to sediment samples, and the specificity and reproducibility of qPCR enumeration were tested. Accurate estimation of potential inhibition within sediment samples was performed to assess the sensitivity of such enumeration by qPCR. Some precautions needed for interpreting qPCR results in the context of paleolimnological approaches are discussed. We concluded that the qPCR assay can be used successfully for the analysis of lake sediments when DNA is well preserved in order to assess the presence and dominance of cyanobacterial and Planktothrix communities.