Distribution of toxigenic cyanobacteria in Alpine lakes and rivers as revealed by molecular screening.

The increasing frequency of cyanobacteria blooms in waterbodies caused by ecosystem eutrophication could endanger human health. This risk can be mitigated by effective monitoring incorporating molecular methods. To date, most molecular studies on toxigenic cyanobacteria have been limited to microcystins (MCs), disregarding other cyanotoxins, to freshwater planktic habitats while ignoring benthic habitats, and to limited geographic areas (usually one or a few specific waterbodies). In this study, we used PCR-based methods including PCR product sequencing and chemical-analytical methods (LC-MS/MS) to screen many plankton (n = 123) and biofilm samples (n = 113) originating from 29 Alpine lakes and 18 rivers for their cyanotoxin production potential. Both mcyE (indicating MC synthesis) and anaC (indicating anatoxin (ATX) synthesis) gene fragments were able to qualitatively predict MC or ATX occurrence. The abundance of mcyE gene fragments was significantly related to MC concentrations in plankton samples (R2 = 0.61). mcyE gene fragments indicative of MC synthesis were most abundant in planktic samples (65 %) and were assigned to the genera Planktothrix and Microcystis. However, mcyE rarely occurred in biofilms of lakes and rivers, i.e., 4 % and 5 %, respectively, and were assigned to Microcystis, Planktothrix, and Nostoc. In contrast, anaC gene fragments occurred frequently in planktic samples (14 % assigned to Tychonema, Phormidium (Microcoleus), and Oscillatoria), but also in biofilms of lakes (49 %) and rivers (18 %) and were assigned to the genera Phormidium, Oscillatoria, and Nostocales. The cyrJ gene fragment indicating cylindrospermopsin synthesis occurred only once in plankton (assigned to Dolichospermum), while saxitoxin synthesis potential was not detected. For plankton samples, monomictic and less eutrophic conditions were positively related to mcyE/MC occurrence frequency, while oligomictic conditions were related to anaC/ATX frequency. The anaC/ATX frequency in biofilm was related to the lake habitats generally showing higher biodiversity as revealed from metabarcoding in a parallel study.

Draft genome sequence of the anatoxin-a producing cyanobacterium Tychonema bourrellyi B0820 isolated from the epilimnion of the deep Alpine Lake Garda.

We report the draft genome sequence of strain B0820 of the cyanobacterium Tychonema bourrellyi isolated from the epilimnion of Lake Garda and assembled from a metagenome of a non-axenic culture. The strain analyzed was shown to produce anatoxin-a, a potent neurotoxin that can cause fatal intoxication in exposed organisms.

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.

Unraveling the Diversity of Eukaryotic Microplankton in a Large and Deep Perialpine Lake Using a High Throughput Sequencing Approach.

The structure of microbial communities, microalgae, heterotrophic protozoa and fungi contributes to characterize food webs and productivity and, from an anthropogenic point of view, the qualitative characteristics of water bodies. Traditionally, in freshwater environments many investigations have been directed to the study of pelagic microalgae ("phytoplankton") and periphyton (i.e., photosynthetic and mixotrophic protists) through the use of light microscopy (LM). While the number of studies on bacterioplankton communities have shown a substantial increase after the advent of high-throughput sequencing (HTS) approaches, the study of the composition, structure, and spatio-temporal patterns of microbial eukaryotes in freshwater environments was much less widespread. Moreover, the understanding of the correspondence between the relative phytoplankton abundances estimated by HTS and LM is still incomplete. Taking into account these limitations, this study examined the biodiversity and seasonality of the community of eukaryotic microplankton in the epilimnetic layer of a large and deep perialpine lake (Lake Garda) using HTS. The analyses were carried out at monthly frequency during 2014 and 2015. The results highlighted the existence of a rich and well diversified community and the presence of numerous phytoplankton taxa that were never identified by LM in previous investigations. Furthermore, the relative abundances of phytoplankton estimated by HTS and LM showed a significant relationship at different taxonomic ranks. In the 2 years of investigation, the temporal development of the whole micro-eukaryotic community showed a clear non-random and comparable distribution pattern, with the main taxonomic groups coherently distributed in the individual seasons. In perspective, the results obtained in this study highlight the importance of HTS approaches in assessing biodiversity and the relative importance of the main protist groups along environmental gradients, including those caused by anthropogenic impacts (e.g., eutrophication and climate change).

A European Multi Lake Survey dataset of environmental variables, phytoplankton pigments and cyanotoxins.

Under ongoing climate change and increasing anthropogenic activity, which continuously challenge ecosystem resilience, an in-depth understanding of ecological processes is urgently needed. Lakes, as providers of numerous ecosystem services, face multiple stressors that threaten their functioning. Harmful cyanobacterial blooms are a persistent problem resulting from nutrient pollution and climate-change induced stressors, like poor transparency, increased water temperature and enhanced stratification. Consistency in data collection and analysis methods is necessary to achieve fully comparable datasets and for statistical validity, avoiding issues linked to disparate data sources. The European Multi Lake Survey (EMLS) in summer 2015 was an initiative among scientists from 27 countries to collect and analyse lake physical, chemical and biological variables in a fully standardized manner. This database includes in-situ lake variables along with nutrient, pigment and cyanotoxin data of 369 lakes in Europe, which were centrally analysed in dedicated laboratories. Publishing the EMLS methods and dataset might inspire similar initiatives to study across large geographic areas that will contribute to better understanding lake responses in a changing environment.

Draft Genome Sequences of Novel Pseudomonas, Flavobacterium, and Sediminibacterium [corrected] Strains from a Freshwater Ecosystem.

Freshwater ecosystems represent 0.01% of the water on Earth, but they support 6% of global biodiversity that is still mostly uncharacterized. Here, we describe the genome sequences of three strains belonging to novel species in the Pseudomonas, Flavobacterium, and Sediminibacterium genera recovered from a water sample of Lake Garda, Italy.

Diversity and Cyclical Seasonal Transitions in the Bacterial Community in a Large and Deep Perialpine Lake.

High-throughput sequencing (HTS) was used to analyze the seasonal variations in the bacterioplankton community composition (BCC) in the euphotic layer of a large and deep lake south of the Alps (Lake Garda). The BCC was analyzed throughout two annual cycles by monthly samplings using the amplification and sequencing of the V3-V4 hypervariable region of the 16S rRNA gene by the MiSeq Illumina platform. The dominant and most diverse bacterioplankton phyla were among the more frequently reported in freshwater ecosystems, including the Proteobacteria, Cyanobacteria, Bacteroidetes, Verrucomicrobia, Actinobacteria, and Planctomycetes. As a distinctive feature, the development of the BCC showed a cyclical temporal pattern in the two analyzed years and throughout the euphotic layer. The recurring temporal development was controlled by the strong seasonality in water temperature and thermal stratification, and by cyclical temporal changes in nutrients and, possibly, by the remarkable annual cyclical development of cyanobacteria and eukaryotic phytoplankton hosting bacterioplankton that characterizes Lake Garda. Further downstream analyses of operational taxonomic units associated to cyanobacteria allowed confirming the presence of the most abundant taxa previously identified by microscopy and/or phylogenetic analyses, as well as the presence of other small Synechococcales/Chroococcales and rare Nostocales never identified so far in the deep lakes south of the Alps. The implications of the high diversity and strong seasonality are relevant, opening perspectives for the definition of common and discriminating patterns characterizing the temporal and spatial distribution in the BCC, and for the application of the new sequencing technologies in the monitoring of water quality in large and deep lakes.

Draft Genome Sequence of the Planktic Cyanobacterium Tychonema bourrellyi, Isolated from Alpine Lentic Freshwater.

We describe here the draft genome sequence of the cyanobacterium Tychonema bourrellyi, assembled from a metagenome of a nonaxenic culture. The strain (FEM_GT703) was isolated from a freshwater sample taken from Lake Garda, Italy. The draft genome sequence represents the first assembled T. bourrellyi strain.

Planktic Tychonema (Cyanobacteria) in the large lakes south of the Alps: phylogenetic assessment and toxigenic potential.

This work allowed assessing a widespread occurrence of Tychonema bourrellyi in the largest lakes south of the Alps (Garda, Iseo, Como and Maggiore). The taxonomy of the species was confirmed adopting a polyphasic approach, which included microscopic examinations, molecular (16S rRNA and rbcLX sequences) and (Lake Garda) ecological characterisations. Over 70% of the 36 isolates of Tychonema sampled from the four lakes tested positive for the presence of genes implicated in the biosynthesis of anatoxins (anaF and/or anaC) and for the production of anatoxin-a (ATX) and homoanatoxin-a (HTX). A detailed analysis carried out in Lake Garda showed strong ongoing changes in the cyanobacterial community, with populations of Tychonema developing with higher biovolumes compared to the microcystins (MCs) producer Planktothrix rubescens Moreover, the time × depth distribution of Tychonema was paralleled by a comparable distribution of ATX and HTX. The increasing importance of Tychonema in Lake Garda was also suggested by the opposite trends of ATX and MCs observed since 2009. These results suggest that radical changes are occurring in the largest lakes south of the Alps. Their verification and implications will require to be assessed by extending a complete experimental work to the other large perialpine lakes.

Variability of microcystin cell quota in metapopulations of Planktothrix rubescens: causes and implications for water management.

In this study, we investigated the relationships between microcystin (MCs) concentrations and the biovolumes of Planktothrix rubescens (BPr) in 2 natural lakes (Pusiano and Garda) and 2 artificially dammed reservoirs (Occhito and Ledro) in Italy. In all the considered water bodies, P. rubescens was the dominant cyanobacterium. All the lakes were characterized by significant relationships between MCs and BPr, with limited variability in the MC quota (the content of MCs per unit of biovolume) within each water body compared with the variability between sites. The results were consistent with the development of specific MC-genotypes, with moderate seasonal and spatial changes in the proportion between toxic and non-toxic strains. The MC cell quota obtained in our work (ECQ, Environmental Cell Quota) were in the same range of values computed on the basis of analyses made on environmental samples dominated by P. rubescens or Planktothrix agardhii, and on isolates of the same two species (<1 to over 10 µg mm(-3)). Besides the usual ordinary least square regressions, models have been evaluated by using quantile regression, a method that allows estimating the conditional median or other quantiles of the response variable. We showed that the use of quantile regressions has different advantages, which included the computation of MC quota based on the whole range of available data, the robustness against outliers, and the ability to estimate models also in cases where there is no or only weak relationships. The highest ECQ values estimated from 95% quantile regressions in specific water bodies might be used to estimate the worst-case MC concentrations from algal abundances. Nevertheless, it was stressed that a realistic assessment of toxicity and potential adverse health effects necessarily should take into account the toxicity potential of the more abundant MC-congeners produced by specific cyanobacteria populations.