Stability of toxin gene proportion in red-pigmented populations of the cyanobacterium Planktothrix during 29 years of re-oligotrophication of Lake Zürich

Deep Cyanobacteria Layers: An Overlooked Aspect of Managing Risks of Cyanobacteria

The risk of human exposure to cyanotoxins is partially influenced by the location of toxin-producing cyanobacteria in waterbodies. Cyanotoxin production can occur throughout the water column, with deep water production representing a potential public health concern, specifically for drinking water supplies. Deep cyanobacteria layers are often unreported, and it remains to be seen if lower incident rates reflect an uncommon phenomenon or a monitoring bias. Here, we examine Sunfish Lake, Ontario, Canada as a case study lake with a known deep cyanobacteria layer. Cyanotoxin and other bioactive metabolite screening revealed that the deep cyanobacteria layer was toxigenic [0.03 μg L^-1 microcystins (max) and 2.5 μg L^-1 anabaenopeptins (max)]. The deep layer was predominantly composed of Planktothrix isothrix (exhibiting a lower cyanotoxin cell quota), with Planktothrix rubescens (exhibiting a higher cyanotoxin cell quota) found at background levels. The co-occurrence of multiple toxigenic Planktothrix species underscores the importance of routine surveillance for prompt identification leading to early intervention. For instance, microcystin concentrations in Sunfish Lake are currently below national drinking water thresholds, but shifting environmental conditions (e.g., in response to climate change or nutrient modification) could fashion an environment favoring P. rubescens, creating a scenario of greater cyanotoxin production. Future work should monitor the entire water column to help build predictive capacities for identifying waterbodies at elevated risk of developing deep cyanobacteria layers to safeguard drinking water supplies.

The Red Harmful Plague in Times of Climate Change: Blooms of the Cyanobacterium Planktothrix rubescens Triggered by Stratification Dynamics and Irradiance

Planktothrix rubescens is a harmful planktonic cyanobacterium, forming concentrated metalimnetic populations in deep oligo- and mesotrophic lakes, even after successful restoration. In Lake Zurich (Switzerland), P. rubescens emerged as a keystone species with annual mass developments since the 1970s. Its success was partly attributed to effects of lake warming, such as changes in thermal stratification and seasonal deep mixing. However, recent observations based on a biweekly monitoring campaign (2009-2020) revealed two massive breakdowns and striking seasonal oscillations of the population. Here, we disentangle positive from negative consequences of secular lake warming and annual variations in weather conditions on P. rubescens dynamics: (i) despite the high survival rates of overwintering populations (up to 25%) during three consecutive winters (2014-2016) of incomplete deep convective mixing, cyanobacterial regrowth during the following stratified season was moderate and not overshooting a distinct standing stock threshold. Moreover, we recorded a negative trend for annual population maxima and total population size, pointing to a potential nutrient limitation after a series of incomplete winter mixing. Thus, the predication of steadily increasing blooms of P. rubescens could not be confirmed for the last decade. (ii) The seasonal reestablishment of P. rubescens was strongly coupled with a timely formation of a stable metalimnion structure, where the first positive net growth in the following productive summer season was observed. The trigger for the vertical positioning of filaments within the metalimnion was irradiance and not maximal water column stability. Repetitive disruptions of the vernal metalimnion owing to unstable weather conditions, as in spring 2019, went in parallel with a massive breakdown of the standing stock and marginal regrowth during thermal stratification. (iii) Driven by light intensity, P. rubescens was entrained into the turbulent epilimnion in autumn, followed by a second peak in population growth. Thus, the typical bimodal growth pattern was still intact during the last decade. Our long-term study highlights the finely tuned interplay between climate-induced changes and variability of thermal stratification dynamics and physiological traits of P. rubescens, determining its survival in a mesotrophic temperate lake.

Restoration, conservation and phytoplankton hysteresis

Phytoplankton growth depends not only upon external factors that are not strongly altered by the presence of phytoplankton, such as temperature, but also upon factors that are strongly influenced by activity of phytoplankton, including photosynthetically active radiation, and the availability of the macronutrients carbon, nitrogen, phosphorus and, for some, silicate. Since phytoplankton therefore modify, and to an extent create, their own habitats, established phytoplankton communities can show resistance and resilience to change, including managed changes in nutrient regimes. Phytoplankton blooms and community structures can be predicted from the overall biogeochemical setting and inputs, but restorations may be influenced by the physiological responses of established phytoplankton taxa to nutrient inputs, temperature, second-order changes in illumination and nutrient recycling. In this review we discuss the contributions of phytoplankton ecophysiology to biogeochemical hysteresis and possible effects on community composition in the face of management, conservation or remediation plans.

Unusual cohabitation and competition between Planktothrix rubescens and Microcystis sp. (cyanobacteria) in a subtropical reservoir (Hammam Debagh) located in Algeria

Succession in bloom-forming cyanobacteria belonging to distant functional groups in freshwater ecosystems is currently an undescribed phenomenon. However in the Hammam Debagh reservoir (Algeria), P. rubescens and Microcystis sp. co-occur and sometimes proliferate. With the aim of identifying the main factors and processes involved in this unusual cohabitation, water samples were collected monthly from February 2013 to June 2015 at the subsurface at four sampling stations and along the entire water column at one sampling station. In addition, the composition of the cyanobacterial communities was estimated by Illumina sequencing of a 16S rRNA gene fragment from samples collected over one year (October 2013-November 2014). This molecular approach showed that the Hammam Debagh reservoir displays high species richness (89 species) but very low diversity due to the high dominance of Microcystis in this community. Furthermore, it appears that Planktothrix rubescens and Microcystis sp. coexisted (from September to January) but proliferated alternately (Spring 2015 for P. rubescens and Spring 2014 and Autumn 2014/2015 for Microcystis). The main factors and processes explaining these changes in bloom-forming species seem to be related to the variation in the depth of the lake during the mixing period and to the water temperatures during the winter prior to the bloom season in spring.

Predicting blooms of toxic cyanobacteria in eutrophic lakes with diverse cyanobacterial communities

We investigated possibility of predicting whether blooms, if they occur, would be formed of microcystin-producing cyanobacteria. DGGE analysis of 16S-ITS and mcyA genes revealed that only Planktothrix and Microcystis possessed mcy-genes and Planktothrix was the main microcystin producer. qPCR analysis revealed that the proportion of cells with mcy-genes in Planktothrix populations was almost 100%. Microcystin concentration correlated with the number of potentially toxic and total Planktothrix cells and the proportion of Planktothrix within all cyanobacteria, but not with the proportion of cells with mcy-genes in total Planktothrix. The share of Microcystis cells with mcy-genes was low and variable in time. Neither the number of mcy-possessing cells, nor the proportion of these cells in total Microcystis, correlated with the concentration of microcystins. This suggests that it is possible to predict whether the bloom in the Masurian Lakes will be toxic based on Planktothrix occurrence. Two species of toxin producing Planktothrix, P. agardhii and P. rubescens, were identified by phylogenetic analysis of 16S-ITS. Based on morphological and ecological features, the toxic Planktothrix was identified as P. agardhii. However, the very high proportion of cells with mcy-genes suggests P. rubescens. Our study reveals the need of universal primers for mcyA genes from environment.

Sedimentary DNA Reveals Cyanobacterial Community Diversity over 200 Years in Two Perialpine Lakes

We reconstructed cyanobacterial community structure and phylogeny using DNA that was isolated from layers of stratified sediments spanning 200 years of lake history in the perialpine lakes Greifensee and Lake Zurich (Switzerland). Community analysis based on amplification and sequencing of a 400-nucleotide (nt)-long 16S rRNA fragment specific to Cyanobacteria revealed operational taxonomic units (OTUs) capturing the whole phylum, including representatives of a newly characterized clade termed Melainabacteria, which shares common ancestry with Cyanobacteria and has not been previously described in lakes. The reconstruction of cyanobacterial richness and phylogenetic structure was validated using a data set consisting of 40 years of pelagic microscopic counts from each lake. We identified the OTUs assigned to common taxa known to be present in Greifensee and Lake Zurich and found a strong and significant relationship (adjusted R² = 0.89; P < 0.001) between pelagic species richness in water and OTU richness in the sediments. The water-sediment richness relationship varied between cyanobacterial orders, indicating that the richness of Chroococcales and Synechococcales may be underestimated by microscopy. PCR detection of the microcystin synthetase gene mcyA confirmed the presence of potentially toxic cyanobacterial taxa over recent years in Greifensee and throughout the last century in Lake Zurich. The approach presented in this study demonstrates that it is possible to reconstruct past pelagic cyanobacterial communities in lakes where the integrity of the sedimentary archive is well preserved and to explore changes in phylogenetic and functional diversity over decade-to-century timescales.

IMPORTANCE

Cyanobacterial blooms can produce toxins that affect water quality, especially under eutrophic conditions, which are a consequence of human-induced climate warming and increased nutrient availability. Lakes worldwide have suffered from regular cyanobacterial blooms over the last century. The lack of long-term data limits our understanding of how these blooms form. We successfully reconstructed the past diversity of whole cyanobacterial communities over two hundred years by sequencing genes preserved in the sediments of two perialpine lakes in Switzerland. We identified changes in diversity over time and validated our results using existing data collected in the same two lakes over the past 40 years. This work shows the potential of our approach for addressing important ecological questions about the effects of a changing environment on lake ecology.

Role of toxic and bioactive secondary metabolites in colonization and bloom formation by filamentous cyanobacteria Planktothrix

Bloom-forming cyanobacteria Planktothrix agardhii and P. rubescens are regularly involved in the occurrence of cyanotoxin in lakes and reservoirs. Besides microcystins (MCs), which inhibit eukaryotic protein phosphatase 1 and 2A, several families of bioactive peptides are produced, thereby resulting in impressive secondary metabolite structural diversity. This review will focus on the current knowledge of the phylogeny, morphology, and ecophysiological adaptations of Planktothrix as well as the toxins and bioactive peptides produced. The relatively well studied ecophysiological adaptations (buoyancy, shade tolerance, nutrient storage capacity) can partly explain the invasiveness of this group of cyanobacteria that bloom within short periods (weeks to months). The more recent elucidation of the genetic basis of toxin and bioactive peptide synthesis paved the way for investigating its regulation both in the laboratory using cell cultures as well as under field conditions. The high frequency of several toxin and bioactive peptide synthesis genes observed within P. agardhii and P. rubescens, but not for other Planktothrix species (e.g. P. pseudagardhii), suggests a potential functional linkage between bioactive peptide production and the colonization potential and possible dominance in habitats. It is hypothesized that, through toxin and bioactive peptide production, Planktothrix act as a niche constructor at the ecosystem scale, possibly resulting in an even higher ability to monopolize resources, positive feedback loops, and resilience under stable environmental conditions. Thus, refocusing harmful algal bloom management by integrating ecological and phylogenetic factors acting on toxin and bioactive peptide synthesis gene distribution and concentrations could increase the predictability of the risks originating from Planktothrix blooms.

Cyanobacteria biennal dynamic in a volcanic mesotrophic lake in central Italy: Strategies to prevent dangerous human exposures to cyanotoxins

Vico Lake, a volcanic meso-eutrophic lake in Central Italy, whose water is used for drinking and recreational activities, experienced the presence of the microcystins (MC) producing cyanobacterium Planktothrix rubescens. In order to assess the human health risks and to provide the local health authorities with a scientific basis for planning tailored monitoring activities, we studied P. rubescens ecology and toxicity for two years. P. rubescens generally dominated the phytoplankton community, alternating with Limnothrix redekei, potentially toxic. P. rubescens was distributed throughout the water column during winter; in summer it produced intense blooms where drinking water is collected (-20 m); here MC were detected all year round (0.5-5 μg/L), with implications for drinking water quality. In surface waters, MC posed no risk for recreational activities in summer, while in winter surface blooms and foams (containing up to 56 μg MC/L) can represent a risk for people and children practicing water sports and for animals consuming raw water. Total phosphorus, phosphate and inorganic nitrogen were not relevant to predict densities nor toxicity; however, a strong correlation between P. rubescens density and aminopeptidase ectoenzymatic activity, an enzyme involved in protein degradation, suggested a role of organic nitrogen for this species. The fraction of potentially toxic population, determined both as mcyB(+)/16SrDNA (10-100%) and as the MC/mcyB(+) cells (0.03-0.79 pg MC/cell), was much more variable than usually observed for P. rubescens. Differently from other Italian and European lakes, the correlation between cell density or the mcyB(+) cells and MC explained only ∼50 and 30% of MC variability, respectively: for Vico Lake, monitoring only cell or the mcyB(+) cell density is not sufficient to predict MC concentrations, and consequently to protect population health. Finally, during a winter bloom one site has been sampled weekly, showing that monthly sampling during such a phase could greatly underestimate the 'hazard'. Our results highlight the need to adopt a stepwise monitoring activity, considering the lake and the cyanobacteria specific features. This activity should be complemented with communication to the public and involvement of stakeholders.

Identification of compounds bound to suspended solids causing sub-lethal toxic effects in Daphnia magna. A field study on re-suspended particles during river floods in Ebro River

Identifying chemicals causing adverse effects in organisms present in water remains a challenge in environmental risk assessment. This study aimed to assess and identify toxic compounds bound to suspended solids re-suspended during a prolonged period of flushing flows in the lower part of Ebro River (NE, Spain). This area is contaminated with high amounts of organochlorine and mercury sediment wastes. Chemical characterization of suspended material was performed by solid phase extraction using a battery of non-polar and polar solvents and analyzed by GC-MS/MS and LC-MS/MS. Mercury content was also determined for all sites. Post-exposure feeding rates of Daphnia magna were used to assess toxic effects of whole and filtered water samples and of re-constituted laboratory water with re-suspended solid fractions. Organochlorine and mercury residues in the water samples increased from upstream to downstream locations. Conversely, toxic effects were greater at the upstream site than downstream of the superfund Flix reservoir. A further analysis of the suspended solid fraction identified a toxic component eluted within the 80:20 methanol:water fraction. Characterization of that toxic component fraction by LC-MS/MS identified the phytotoxin anatoxin-a, whose residue levels were correlated with observed feeding inhibition responses. Further feeding inhibition assays conducted in the lab using anatoxin-a produced from Planktothrix agardhii, a filamentous cyanobacteria, confirmed field results. This study provides evidence that in real field situation measured contaminant residues do not always agree with toxic effects.

Integrating phylogeny, geographic niche partitioning and secondary metabolite synthesis in bloom-forming Planktothrix

Toxic freshwater cyanobacteria form harmful algal blooms that can cause acute toxicity to humans and livestock. Globally distributed, bloom-forming cyanobacteria Planktothrix either retain or lose the mcy gene cluster (encoding the synthesis of the secondary metabolite hepatotoxin microcystin or MC), resulting in a variable spatial/temporal distribution of (non)toxic genotypes. Despite their importance to human well-being, such genotype diversity is not being mapped at scales relevant to nature. We aimed to reveal the factors influencing the dispersal of those genotypes by analyzing 138 strains (from Europe, Russia, North America and East Africa) for their (i) mcy gene cluster composition, (ii) phylogeny and adaptation to their habitat and (iii) ribosomally and nonribosomally synthesized oligopeptide products. Although all the strains from different species contained at least remnants of the mcy gene cluster, various phylogenetic lineages evolved and adapted to rather specific ecological niches (for example, through pigmentation and gas vesicle protein size). No evidence for an increased abundance of specific peptides in the absence of MC was found. MC and peptide distribution rather depended on phylogeny, ecophysiological adaptation and geographic distance. Together, these findings provide evidence that MC and peptide production are primarily related to speciation processes, while within a phylogenetic lineage the probability that strains differ in peptide composition increases with geographic distance.

CYANOCHIP: an antibody microarray for high-taxonomical-resolution cyanobacterial monitoring

Cyanobacteria are Gram-negative photosynthetic prokaryotes that are widespread on Earth. Eutrophication and global warming make some aquatic ecosystems behave as bioreactors that trigger rapid and massive cyanobacterial growth with remarkable economic and health consequences. Rapid and efficient early warning systems are required to support decisions by water body authorities. We have produced 17 specific antibodies to the most frequent cyanobacterial strains blooming in freshwater ecosystems, some of which are toxin producers. A sandwich-type antibody microarray immunoassay (CYANOCHIP) was developed for the simultaneous testing of any of the 17 strains, or other closely related strains, in field samples from different habitats (water, rocks, and sediments). We titrated and tested all of the antibodies in succession using a fluorescent sandwich microarray immunoassay. Although most showed high specificity, we applied a deconvolution method based on graph theory to disentangle the few existing cross-reactions. The CYANOCHIP sensitivity ranged from 10(2) to 10(4) cells mL(-1), with most antibodies detecting approximately 10(2) cells mL(-1). We validated the system by testing multiple isolates and crude natural samples from freshwater reservoirs and rocks, both in the laboratory and by in situ testing in the field. The results demonstrated that CYANOCHIP is a valuable tool for the sensitive and reliable detection of cyanobacteria for early warning and research purposes.

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.

Genetic variability of microcystin biosynthesis genes in Planktothrix as elucidated from samples preserved by heat desiccation during three decades

Historic samples of phytoplankton can provide information on the abundance of the toxigenic genotypes of cyanobacteria in dependence on increased or decreased eutrophication. The analysis of a time-series from preserved phytoplankton samples by quantitative PCR (qPCR) extends observation periods considerably. The analysis of DNA from heat-desiccated samples by qPCR can be aggravated by point substitutions or the fragmentation of DNA introduced by the high temperature. In this study, we analyzed whether the heat desiccation of the cellular material of the cyanobacterium Planktothrix sp. introduced potential errors to the template DNA that is used for qPCR within (i) 16S rDNA and phycocyanin genes and (ii) the mcyA gene indicative of the incorporation of either dehydrobutyrine (Dhb) or N-methyl-dehydroalanine (Mdha) in position 7, and (ii) the mcyB gene, which is indicative of homotyrosine (Hty) in position 2 of the microcystin (MC) molecule. Due to high temperature desiccation, the deterioration of the DNA template quality was rather due to fragmentation than due to nucleotide substitutions. By using the heat-desiccated samples of Lake Zürich, Switzerland the abundance of the Dhb, Mdha and Hty genotypes was determined during three decades (1977-2008). Despite major changes in the trophic state of the lake resulting in a major increase of the total Planktothrix population density, the proportion of these genotypes encoding the synthesis of different MC congeners showed high stability. Nevertheless, a decline of the most abundant mcyA genotype indicative of the synthesis of Dhb in position 7 of the MC molecule was observed. This decline could be related to the gradual incline in the proportion of a mutant genotype carrying a 1.8kbp deletion of this gene region. The increase of this mcyA (Dhb) gene deletion mutant has been minor so far, however, and likely did not affect the overall toxicity of the population.