Presence of potential toxin-producing cyanobacteria in an oligo-mesotrophic lake in Baltic Lake District, Germany: an ecological, genetic and toxicological survey

Lake browning counteracts cyanobacteria responses to nutrients: Evidence from phytoplankton dynamics in large enclosure experiments and comprehensive observational data

Lakes worldwide are affected by multiple stressors, including climate change. This includes massive loading of both nutrients and humic substances to lakes during extreme weather events, which also may disrupt thermal stratification. Since multi-stressor effects vary widely in space and time, their combined ecological impacts remain difficult to predict. Therefore, we combined two consecutive large enclosure experiments with a comprehensive time-series and a broad-scale field survey to unravel the combined effects of storm-induced lake browning, nutrient enrichment and deep mixing on phytoplankton communities, focusing particularly on potentially toxic cyanobacterial blooms. The experimental results revealed that browning counteracted the stimulating effect of nutrients on phytoplankton and caused a shift from phototrophic cyanobacteria and chlorophytes to mixotrophic cryptophytes. Light limitation by browning was identified as the likely mechanism underlying this response. Deep-mixing increased microcystin concentrations in clear nutrient-enriched enclosures, caused by upwelling of a metalimnetic Planktothrix rubescens population. Monitoring data from a 25-year time-series of a eutrophic lake and from 588 northern European lakes corroborate the experimental results: Browning suppresses cyanobacteria in terms of both biovolume and proportion of the total phytoplankton biovolume. Both the experimental and observational results indicated a lower total phosphorus threshold for cyanobacterial bloom development in clearwater lakes (10-20 μg P L-1 ) than in humic lakes (20-30 μg P L-1 ). This finding provides management guidance for lakes receiving more nutrients and humic substances due to more frequent extreme weather events.

Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies

BACKGROUND

The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82].

RESULTS

At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%).

CONCLUSIONS

Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.

Species-Level Spatio-Temporal Dynamics of Cyanobacteria in a Hard-Water Temperate Lake in the Southern Baltics

Cyanobacteria are important primary producers in temperate freshwater ecosystems. However, studies on the seasonal and spatial distribution of cyanobacteria in deep lakes based on high-throughput DNA sequencing are still rare. In this study, we combined monthly water sampling and monitoring in 2019, amplicon sequence variants analysis (ASVs; a proxy for different species) and quantitative PCR targeting overall cyanobacteria abundance to describe the seasonal and spatial dynamics of cyanobacteria in the deep hard-water oligo-mesotrophic Lake Tiefer See, NE Germany. We observed significant seasonal variation in the cyanobacterial community composition (p < 0.05) in the epi- and metalimnion layers, but not in the hypolimnion. In winter-when the water column is mixed-picocyanobacteria (Synechococcus and Cyanobium) were dominant. With the onset of stratification in late spring, we observed potential niche specialization and coexistence among the cyanobacteria taxa driven mainly by light and nutrient dynamics. Specifically, ASVs assigned to picocyanobacteria and the genus Planktothrix were the main contributors to the formation of deep chlorophyll maxima along a light gradient. While Synechococcus and different Cyanobium ASVs were abundant in the epilimnion up to the base of the euphotic zone from spring to fall, Planktothrix mainly occurred in the metalimnetic layer below the euphotic zone where also overall cyanobacteria abundance was highest in summer. Our data revealed two potentially psychrotolerant (cold-adapted) Cyanobium species that appear to cope well under conditions of lower hypolimnetic water temperature and light as well as increasing sediment-released phosphate in the deeper waters in summer. The potential cold-adapted Cyanobium species were also dominant throughout the water column in fall and winter. Furthermore, Snowella and Microcystis-related ASVs were abundant in the water column during the onset of fall turnover. Altogether, these findings suggest previously unascertained and considerable spatiotemporal changes in the community of cyanobacteria on the species level especially within the genus Cyanobium in deep hard-water temperate lakes.

Unmasking the identity of toxigenic cyanobacteria driving a multi-toxin bloom by high-throughput sequencing of cyanotoxins genes and 16S rRNA metabarcoding

Cyanobacterial harmful algal blooms (CyanoHABs) are complex communities that include coexisting toxic and non-toxic strains only distinguishable by genetic methods. This study shows a water-management oriented use of next generation sequencing (NGS) to specifically pinpoint toxigenic cyanobacteria within a bloom simultaneously containing three of the most widespread cyanotoxins (the hepatotoxins microcystins, MCs; and the neurotoxins anatoxin-a, ATX, and saxitoxins, STXs). The 2013 summer bloom in Rosarito reservoir (Spain) comprised 33 cyanobacterial OTUs based on 16S rRNA metabarcoding, 7 of which accounted for as much as 96.6% of the community. Cyanotoxins and their respective biosynthesis genes were concurrently present throughout the entire bloom event including: MCs and mcyE gene; ATX and anaF gene; and STXs and sxtI gene. NGS applied to amplicons of cyanotoxin-biosynthesis genes unveiled 6 toxigenic OTUs, comprising 3 involved in MCs production (Planktothrix agardhii and 2 Microcystis spp.), 2 in ATX production (Cuspidothrix issatschenkoi and Phormidium/Tychonema spp.) and 1 in STXs production (Aphanizomenon gracile). These toxigenic taxa were also present in 16S rRNA OTUs list and their relative abundance was positively correlated with the respective toxin concentrations. Our results point at MC-producing P. agardhii and ATX-producing C. issatschenkoi as the main contributors to the moderate toxin concentrations observed, and suggest that their distribution in Southern Europe is broader than previously thought. Our findings also stress the need for monitoring low-abundance cyanobacteria (<1% relative abundance) in cyanotoxicity studies, and provide novel data on the presence of picocyanobacteria and potentially ATX-producing benthic taxa (e.g., Phormidium) in deep thermally-stratified water bodies. This study showcases a straightforward use of amplicon metagenomics of cyanotoxin biosynthesis genes in a multi-toxin bloom thus illustrating the broad applicability of NGS for water management in risk-oriented monitoring of CyanoHABs.

High-throughput DNA sequencing reveals the dominance of pico- and other filamentous cyanobacteria in an urban freshwater Lake

The current study presents findings related to algal blooms in a fresh water lake, which has been experiencing severe cyanobacterial blooms (CyanoHABs). Primarily, picocyanobacteria belonging to the genus Synechococcus and filamentous cyanobacterial group belonging to Aphanizomenon and Dolichospermum dominated top water column during non-bloom and bloom periods respectively. The dominance of Synechococcus in early summer informs that blooming in Utah Lake starts in early summer and then later is taken over by other bloom-forming cyanobacteria, such as species belonging to the genus Aphanizomenon. A strong negative correlation (r = -0.9, p < 0.001) was found between the occurrence of Aphanizomenon and Synechococcus which correlates very well with the fact that the blooms of these two different cyanobacteria never coexisted. The predominance of cyanobacteria in 2017 was attributed more to temperature (r = 0.18, p < 0.001). The Actinobacteria was negatively correlated with primary production and high chlorophyll a concentration. Flavobacterium and Limnohabitans were the main phytoplankton colonizers and predators detected that could secrete extracellular enzymes to degrade algal exudates (such as proteins and polysaccharides). Additionally, cyanotoxins producers Microcystis aeruginosa and Planktothrix accounted for up to 12.43% and 7.04% of total cyanobacteria abundance during blooms. The relative abundance of chloroplast reads was overall lower than the cyanobacteria reads, except for the May 5th sampling in 2017. There was inter-annual variability in the bloom-associated heterotrophic bacterial populations, but these populations were consistent with bloom-associated bacterial populations found in other lakes. Community diversity analysis for both Shannon and Simpson indices indicated lower community diversity during the bloom period. The beta diversity conducted by PCoA and UPGMA trees suggested the significant temporal rather than spatial impacts on shaping the phytoplankton community structures during the summer season.

A Systematic Literature Review for Evidence of Aphanizomenon flos-aquae Toxigenicity in Recreational Waters and Toxicity of Dietary Supplements: 2000⁻2017

Previous studies of recreational waters and blue-green algae supplements (BGAS) demonstrated co-occurrence of Aphanizomenon flos-aquae (AFA) and cyanotoxins, presenting exposure risk. The authors conducted a systematic literature review using a GRADE PRISMA-p 27-item checklist to assess the evidence for toxigenicity of AFA in both fresh waters and BGAS. Studies have shown AFA can produce significant levels of cylindrospermopsin and saxitoxin in fresh waters. Toxicity studies evaluating AFA-based BGAS found some products carried the mcyE gene and tested positive for microcystins at levels ≤ 1 μg microcystin (MC)-LR equivalents/g dry weight. Further analysis discovered BGAS samples had cyanotoxins levels exceeding tolerable daily intake values. There is evidence that Aphanizomenon spp. are toxin producers and AFA has toxigenic genes such as mcyE that could lead to the production of MC under the right environmental conditions. Regardless of this ability, AFA commonly co-occur with known MC producers, which may contaminate BGAS. Toxin production by cyanobacteria is a health concern for both recreational water users and BGAS consumers. Recommendations include: limit harvesting of AFA to months when toxicity is lowest, include AFA in cell counts during visible blooms, and properly identify cyanobacteria species using 16S rRNA methods when toxicity levels are higher than advisory levels.

Recovering Genomics Clusters of Secondary Metabolites from Lakes Using Genome-Resolved Metagenomics

Metagenomic approaches became increasingly popular in the past decades due to decreasing costs of DNA sequencing and bioinformatics development. So far, however, the recovery of long genes coding for secondary metabolites still represents a big challenge. Often, the quality of metagenome assemblies is poor, especially in environments with a high microbial diversity where sequence coverage is low and complexity of natural communities high. Recently, new and improved algorithms for binning environmental reads and contigs have been developed to overcome such limitations. Some of these algorithms use a similarity detection approach to classify the obtained reads into taxonomical units and to assemble draft genomes. This approach, however, is quite limited since it can classify exclusively sequences similar to those available (and well classified) in the databases. In this work, we used draft genomes from Lake Stechlin, north-eastern Germany, recovered by MetaBat, an efficient binning tool that integrates empirical probabilistic distances of genome abundance, and tetranucleotide frequency for accurate metagenome binning. These genomes were screened for secondary metabolism genes, such as polyketide synthases (PKS) and non-ribosomal peptide synthases (NRPS), using the Anti-SMASH and NAPDOS workflows. With this approach we were able to identify 243 secondary metabolite clusters from 121 genomes recovered from our lake samples. A total of 18 NRPS, 19 PKS, and 3 hybrid PKS/NRPS clusters were found. In addition, it was possible to predict the partial structure of several secondary metabolite clusters allowing for taxonomical classifications and phylogenetic inferences. Our approach revealed a high potential to recover and study secondary metabolites genes from any aquatic ecosystem.

Temporal Dynamics of the Microbial Community Composition with a Focus on Toxic Cyanobacteria and Toxin Presence during Harmful Algal Blooms in Two South German Lakes

Bacterioplankton plays an essential role in aquatic ecosystems, and cyanobacteria are an influential part of the microbiome in many water bodies. In freshwaters used for recreational activities or drinking water, toxic cyanobacteria cause concerns due to the risk of intoxication with cyanotoxins, such as microcystins. In this study, we aimed to unmask relationships between toxicity, cyanobacterial community composition, and environmental factors. At the same time, we assessed the correlation of a genetic marker with microcystin concentration and aimed to identify the main microcystin producer. We used Illumina MiSeq sequencing to study the bacterioplankton in two recreational lakes in South Germany. We quantified a microcystin biosynthesis gene (mcyB) using qPCR and linked this information with microcystin concentration to assess toxicity. Microcystin biosynthesis gene (mcyE)-clone libraries were used to determine the origin of microcystin biosynthesis genes. Bloom toxicity did not alter the bacterial community composition, which was highly dynamic at the lowest taxonomic level for some phyla such as Cyanobacteria. At the OTU level, we found distinctly different degrees of temporal variation between major bacteria phyla. Cyanobacteria and Bacteroidetes showed drastic temporal changes in their community compositions, while the composition of Actinobacteria remained rather stable in both lakes. The bacterial community composition of Alpha- and Beta-proteobacteria remained stable over time in Lake Klostersee, but it showed temporal variations in Lake Bergknappweiher. The presence of potential microcystin degraders and potential algicidal bacteria amongst prevalent Bacteroidetes and Alphaproteobacteria implied a role of those co-occurring heterotrophic bacteria in cyanobacterial bloom dynamics. Comparison of both lakes studied revealed a large shared microbiome, which was shaped toward the lake specific community composition by environmental factors. Microcystin variants detected were microcystin-LR, -RR, and -YR. The maximum microcystin concentrations measured was 6.7 μg/L, a value still acceptable for recreational waters but not drinking water. Microcystin concentration correlated positively with total phosphorus and mcyB copy number. We identified low abundant Microcystis sp. as the only microcystin producer in both lakes. Therefore, risk assessment efforts need to take into account the fact that non-dominant species may cause toxicity of the blooms observed.

IMPROVING CYANOBACTERIA AND CYANOTOXIN MONITORING IN SURFACE WATERS FOR DRINKING WATER SUPPLY

Cyanobacteria in fresh water can cause serious threats to drinking water supplies. Managing cyanobacterial blooms particularly at small drinking water treatment plants is challenging. Because large amount of cyanobacteria may cause clogging in the treatment process and various cyanotoxins are hard to remove, while they may cause severe health problems. There is lack of instructions of what cyanobacteria/toxin amount should trigger what kind of actions for drink-ing water management except for Microcystins. This demands a Cyanobacteria Management Tool (CMT) to help regula-tors/operators to improve cyanobacteria/cyanotoxin monitoring in surface waters for drinking water supply. This project proposes a CMT tool, including selecting proper indicators for quick cyanobacteria monitoring and verifying quick analysis methods for cyanobacteria and cyanotoxin. This tool is suggested for raw water management regarding cyano-bacteria monitoring in lakes, especially in boreal forest climate. In addition, it applies to regions that apply international WHO standards for water management. In Swedish context, drinking water producers which use raw water from lakes that experience cyanobacterial blooms, need to create a monitoring routine for cyanobacteria/cyanotoxin and to monitor beyond such as Anatoxins, Cylindrospermopsins and Saxitoxins. Using the proposed CMT tool will increase water safety at surface water treatment plants substantially by introducing three alerting points for actions. CMT design for each local condition should integrate adaptive monitoring program.

Could the presence of larger fractions of non-cyanobacterial species be used as a predictor of microcystin production under variable nutrient regimes?

The occurrence of cyanobacteria and microcystin is highly dynamic in natural environments and poses one of the biggest challenges to water resource management. While a number of drivers are known to be responsible for the occurrence of cyanobacterial blooms, the drivers of microcystin production are not adequately known. This study aims to quantify the effects of the changes in the structures of phytoplankton and cyanobacterial communities on the dynamics of microcystin production under highly variable nutrient concentration. In our study, nutrient variability could explain 64% of the variability in microcystin production. When changes in the fractions of non-cyanobacteria versus cyanobacteria genera were additionally included, 80% of the variability in microcystin production could be explained; under high nutrient concentrations, non-cyanobacterial phytoplankton groups were dominant over cyanobacteria and cyanobacteria produced more toxins. In contrast, changes in the cyanobacterial community structures could only explain a further 4% of the dynamics of microcystin production. As such, the dominance of non-cyanobacterial groups appears to be a useful factor to explain microcystin occurrence in addition to traditionally used factors such as absolute cyanobacterial cell numbers, especially when the nutrient regime is taken into account. This information could help to further refine the risk assessment frameworks which are currently used to manage the risk posed by cyanobacterial blooms.