Kristin J.
Painter
*a,
Jason J.
Venkiteswaran
b,
Dana F.
Simon
c,
Sung
Vo Duy
c,
Sébastien
Sauvé
c and
Helen M.
Baulch
a
aSchool of Environment and Sustainability, Global Institute for Water Security, University of Saskatchewan, Saskatoon, SK S7N 5C8, Canada. E-mail: kristin.painter@usask.ca
bDepartment of Geography and Environmental Studies, Wilfrid Laurier University, Waterloo, ON N2L 3C5, Canada
cDepartment of Chemistry, Université de Montréal, Montréal, QC H2V 0B3, Canada
First published on 12th July 2022
Cyanobacterial blooms present challenges for water treatment, especially in regions like the Canadian prairies where poor water quality intensifies water treatment issues. Buoyant cyanobacteria that resist sedimentation present a challenge as water treatment operators attempt to balance pre-treatment and toxic disinfection by-products. Here, we used microscopy to identify and describe the succession of cyanobacterial species in Buffalo Pound Lake, a key drinking water supply. We used indicator species analysis to identify temporal grouping structures throughout two sampling seasons from May to October 2018 and 2019. Our findings highlight two key cyanobacterial bloom phases – a mid-summer diazotrophic bloom of Dolichospermum spp. and an autumn Planktothrix agardhii bloom. Dolichospermum crassa and Woronichinia compacta served as indicators of the mid-summer and autumn bloom phases, respectively. Different cyanobacterial metabolites were associated with the distinct bloom phases in both years: toxic microcystins were associated with the mid-summer Dolichospermum bloom and some newly monitored cyanopeptides (anabaenopeptin A and B) with the autumn Planktothrix bloom. Despite forming a significant proportion of the autumn phytoplankton biomass (>60%), the Planktothrix bloom had previously not been detected by sensor or laboratory-derived chlorophyll-a. Our results demonstrate the power of targeted taxonomic identification of key species as a tool for managers of bloom-prone systems. Moreover, we describe an autumn Planktothrix agardhii bloom that has the potential to disrupt water treatment due to its evasion of detection. Our findings highlight the importance of identifying this autumn bloom given the expectation that warmer temperatures and a longer ice-free season will become the norm.
Environmental significanceCyanobacterial blooms threaten water supplies, especially in the vast, yet understudied prairie regions where pressures on water security mean water quality is often considered secondarily to water quantity. Understanding the structure, composition, and metabolites of cyanobacterial communities is essential for water treatment operators who must balance risks of high cyanobacterial biomass with those of toxic disinfection by-products in systems with inherently poor and declining water quality. Our work demonstrated two cyanobacterial bloom periods: a toxic mid-summer bloom of nitrogen-fixers, and a secondary bloom of the genera Planktothrix that, despite in-lake monitoring with high-frequency sensors, was previously undetected. This autumn bloom highlights a knowledge gap that must be considered by researchers and practitioners as open-water seasons lengthen with climate warming. |
The impacts of cyanobacterial species are of particular concern due to their ability to produce toxic metabolites, thereby presenting potential human health risks. Shallow prairie lakes appear to be especially susceptible to toxin-producing cyanobacterial blooms, for example, Loewen et al.14 found a negative correlation between depth and toxin production in Alberta lakes, a link they attributed to the polymictic nature of many shallow prairie lakes whereby frequent mixing provides phosphorus (P) enrichment from deeper waters to phytoplankton. Plentiful P and limited inflow of nitrogen (N) from the watershed throughout the summer months in many prairie systems favours the growth of N-fixing genera such as Dolichospermum and Aphanizomenon known to produce toxins.7,15 Shallow prairie lakes are also more susceptible to warming,16 a driver that has also been linked to blooms17 thus concerns also remain about non-diazotrophic species (e.g., Planktothrix, Microcystis) that have the potential to proliferate and contribute to toxin risk. Consequently, improved understanding of the species and associated toxins is needed to guide more in-depth analysis and management of these systems.
The physicochemical properties of many toxin-producing species can also impede their removal during drinking water treatment. For example, filamentous species (e.g., Aphanizomenon and Planktothrix spp.) can foul filtering infrastructure while others, e.g., Microcystis genera, are resistant to coagulation making their removal during sedimentation processes challenging.18–20 High phytoplankton biomass can be particularly difficult to deal with for water 2treatment operators who must manage particulate loads and balance the use of disinfection processes such as pre-chlorination with the risk of toxic disinfection by-products. When present in large quantities, phytoplankton biomass may contribute to rising floc – aggregates of solid material that cannot be effectively overcome by coagulation (and subsequent sedimentation) alone. Pre-chlorination at the outset of the treatment process may reduce problems associated with bloom impacts on processes such as rising flocs by killing buoyant algae,21 however, when chlorine reacts with naturally occurring organic matter, toxic trihalomethanes (THMs) are produced.22 Operators must meet THM limits while maintaining adequate drinking water supplies. Moreover, some toxin producing species may persist after oxidation whether they were dominant species or not,23 adding to the complexity of managing bloom-prone water supplies.
Though cyanobacterial blooms are common in shallow prairie lakes, the intra- and interannual variation inherent to the prairie region contributes to the difficulty in predicting the severity of blooms.24 Treatment operations may need to be altered quickly in response and the availability of additional tools with high ease-of-use to managers is imperative. Non-specific indicators of algal biomass, such as chlorophyll-a (chl-a), have enabled managers to establish links between in-lake algal blooms and water treatment issues;25,26 however, such measurements provide limited information about the organisms responsible for treatment disruptions. Operators of municipal water supplies in Canada also frequently employ routine in-house microscopy to estimate cyanobacterial phytoplankton cell counts as part of source water monitoring.27 Thus, additional identification of the presence of key species indicative of cyanobacterial blooms could provide managers with an early indicator of those blooms that have the potential to affect water treatment operations, something not possible with non-specific biomass indicators.
Here, we examine the structure of the taxonomically identified algal community throughout two ice-free seasons in a shallow, polymictic prairie lake that serves as a major drinking water supply. We use a combination of clustering, indicator species analysis, and ordination to describe the succession of algal community abundance over time and identify potential indicator species common across years and time periods. We also provide new information about the types and concentrations of cyanotoxins present and describe their correlation to the cyanobacterial community. Finally, we contrast the common algal biomass indicator chl-a and the findings of our taxonomic analysis. The objective of this paper is to identify the seasonal evolution of harmful cyanobacteria and their toxins in a crucial drinking water supply and demonstrate the utility of targeted taxonomic identification of key species as a tool for managers.
![]() | ||
Fig. 1 Map of Buffalo Pound Lake, Saskatchewan, Canada, and its gross (grey) and effective (yellow) drainage area, including upstream Lake Diefenbaker. The inset map (red outline) indicates the location of the Buffalo Pound Water Treatment Plant (green WTP symbol) and the sampling site (star symbol; 50°35′8.8′′ N, 105°23′0.24′′ W) near its intakes. The Buffalo Pound Dam, which controls lake water level, is shown at the southeast end of the lake (black and white dam symbol). Geospatial data used to create this map was downloaded from CanVec32 and ArcGIS Online.33 |
Buffalo Pound Lake serves as an important water supply in a relatively water scarce, semi-arid to subhumid region where water quality is affected by eutrophication and annual cyanobacterial blooms.8,34 The lake is an important drinking water source, supplying drinking water to a quarter of the province's population, approximately 260000 people, in addition to substantive water withdrawals to support industry including local mines. The annual range of known concentrations for chl-a is approx. 5.0 to 170 μg L−1 and <2.0 to 340 μg L−1 for total P.35 Additionally, dissolved organic matter (DOM) in Buffalo Pound is relatively high for a drinking water source. The poor source water quality creates major challenges for water treatment, for example in 2015 water shortages occurred resulting from early bloom conditions co-occurring with extreme stratification.36 The drinking water treatment plant uses a process which includes pre-chlorination, coagulation (with alum or polyaluminum chloride), flocculation, clarification, filtration, and then final disinfection before the water is distributed to the cities of Moose Jaw and Regina, Saskatchewan. Cascade degasification and granular activated carbon is used seasonally, near the start and end of the process train respectively while powdered activated carbon is used on an as-needed basis to due reduce excess taste and odour compounds. While pre-chlorination was briefly discontinued to manage disinfection by-products (THMs); issues with flocculation and clarification have occurred periodically associated with bloom periods, and vigilance is required to manage the concurrent issues of blooms of buoyant cyanobacteria, transient lake stratification, and high DOM.
Water samples collected for intra- and extracellular cyanobacterial metabolites were filtered (100 ml) on the day of collection through a 0.45 μm hydrophilic polypropylene (GHP) membrane filter,37 without adjusting sample pH (mean pH = 8.29 ± 0.31 in 2018 and 8.64 ± 0.37 in 2019). Filters were enclosed in a Petri dish and frozen until time of analysis for intracellular metabolites. The filtrate phase was collected in PETG amber bottles for dissolved (extracellular) toxin analysis and frozen at −20 °C until analysis.38 Samples were shipped to Université de Montréal in Montréal, Canada for analysis of total microcystins and 17 cyanobacterial metabolites (see also Texts S1–S2 and Table S1†). The analysis of filters proceeded as follows. Following addition of 5 mL of 90:
10 MeOH/water, the filters were ultrasonicated (15 min), centrifuged (10 min, 6000 rpm), and the liquid phase transferred to a fresh extraction tube. The extraction process was repeated once (two extraction cycles), and the combined solution evaporated to dryness prior to reconstitution in 5 mL of HPLC-water. Following filtration through GHP and addition of internal standards, aqueous solutions were analyzed by on-line solid-phase extraction (SPE) ultra-high-performance liquid chromatography (UHPLC) tandem mass spectrometry (MS/MS) or high-resolution mass spectrometry (HRMS), using the same instrumental methods as for dissolved samples.
For total microcystins, a 5 mL aliquot of water was filtered through a GHP; a sub-aliquot was amended with MMPB-d3 prior to permanganate oxidation (to generate the 2-methyl-3-methoxy-4-phenylbutyric acid (MMPB) moiety), and subsequently analyzed by on-line SPE – UHPLC-MS/MS (Thermo TSQ Quantiva).39,40 For individual cyanotoxins (anabaenopeptin A, anabaenopeptin B, anatoxin-a, homoanatoxin-a, cylindrospermopsin, and microcystins -dm-RR, -RR, -YR, -HtyR, -LR, -dm-LR, -HilR, -WR, -LA, -LY, -LF, and -LW), a 5 mL aliquot of water was filtered through a GHP; a sub-aliquot was amended with an internal standard mix and analyzed by on-line SPE – UHPLC-HRMS (Thermo Orbitrap Q-Exactive).37 Matrix-matched calibration (additions to blank lake water, calibration range of LOQ-1000 ng L−1) with internal standard correction was used to quantify cyanotoxins, with accuracies within 81–113%.37 All samples were analyzed in duplicate. If specific samples surpassed the high-end calibration level, a new sample preparation was performed by using a diluted aliquot of the sample for analysis; concentration in the initial sample was then derived applying the dilution factor.
Phytoplankton samples were preserved with Lugol's iodine upon collection, concentrated, and shipped to Plankton R Us in Winnipeg, Canada for enumeration and identification by D. L. Findlay. Larsen et al.41 provide a detailed description of enumeration methods. In brief, samples were enumerated and identified to the species level using an inverted microscope following the methods described in Findlay and Kling.42 In addition to viable vegetative cells (i.e., with intact chloroplast), cyanobacterial heterocytes were also enumerated if present. Cell counts were converted to wet weight biomass by approximating cell volume using measurements of individual cells and applying the geometric formula best fitted to the shape of the cell per Vollenweider.43
Chl-a samples were analyzed at the Saskatchewan Water Chemistry and Ecology laboratory at the University of Saskatchewan Global Institute for Water Security in Saskatoon, Canada. Samples were prepared by filtering 100–300 ml of sampled lake water through glass fiber filters (GF/Fs) and freezing filters until analysis. Prior to analysis, GF/Fs were digested in 95% ethanol for 24 hours44 and absorbance of the resulting solution was measured at 649, 655 and 750 nm on a Shimadzu UV-1601PC spectrophotometer.
IndVals were calculated in two ways, first considering only the associations between species and each group, and second considering possible associations between species and single groups plus group combinations. Group combinations were included because indicator species associated to two or more groups allow for interpretation of similarity between groupings because they characterize a higher-level grouping structure than would otherwise be indicated by single groups only.47
Prior to analysis of individual intracellular cyanotoxin data, values below LOD were replaced with zero. In the rare instance that values (n = 3 occurrences for individual MC congeners) were greater than LOD but less than level of quantification (LOQ), they were replaced with their respective LOQs (see Table S1†). A squareroot transformation was then applied to the dataset to reduce skewness and the envfit function in R package vegan50 was used to fit toxin vectors onto a non-metric multidimensional scaling (NMDS) ordination of the cyanobacterial phytoplankton community using permutation tests. Finally, non-parametric Spearman rank correlations (rs) were used to assess relationships between the biomass of individual cyanobacterial species and the most predominant toxin concentrations (total microcystins, anabaenopeptins A and B).
Despite differences in the overall number of species detected, partitioning placed the phytoplankton community into three similar temporally defined groups in both years: Group A consisting largely of early-season algal species; Group B contained the early to mid-summer blooming diazotrophic cyanophyte species; and Group C was dominated by the late-season non-diazotrophic cyanophyte species (Fig. 2). In 2018 (Fig. 2a), Group A contained only the first two sampling days (late May and early June) and was dominated by chrysophytes (73% of Group A biomass), especially Dinobryon sertularia (43%) and Chrysochromulina parva (11%). Cyanophytes accounted for 75% and 70% of 2018 Group B and C biomass, respectively. Group B spanned from mid-June to mid-July (DOY 170 to 193) and also included August 8th (DOY 220). Three cyanophyte species accounted for much of the biomass: Dolichospermum flos-aquae and D. crassa accounted for 44% of Group B biomass while Planktothrix agardhii contributed 29%. Group C, which included mid-July to mid-September (DOY 197 to 260 but excluding 220), was associated with a shift in the cyanophyte community from Dolichospermum spp. dominant to P. agardhii dominant. P. agardhii contributed 63% of the total biomass in Group C while the Dolichospermum species contributed less than 5% collectively.
In 2019 (Fig. 2b), Group A spanned the first four sampling days, extending later into June (DOY 170) than in 2018. Cryptophytes and chrysophytes accounted for nearly half of the Group A biomass in 2019 (31% and 16% respectively), however, cyanophytes contributed a much larger proportion (43%) than was observed for Group A in 2018 (5%), largely due to D. flos-aquae (41% of Group A biomass). Like in 2018, in 2019 cyanophytes accounted for much of the biomass for the remainder of the year: 68% of Group B biomass and 74% of Group C biomass. Group B spanned from early July to late August (DOY 186 to 242) and was again dominated by diazotrophs including D. crassa (33%), D. flos-aquae (16%) and Aphanizomenon flos-aquae (14%). Unlike 2018, P. agardhii only accounted for 2% of Group B biomass in 2019; however, it again became the wholly dominant species in Group C (DOY 264 to 276) contributing 63% to 88% of the total biomass from mid-September to October.
Individual indicator species were also similar between years. For example, the chrysophyte species Dinobryon sertularia was a significant Group A indicator in both 2018 (IndVal = 0.963, p = 0.001) and 2019 (IndVal = 1.00, p = 0.001). The diazotrophs D. crassa (IndVal = 0.931, p = 0.001) and A. flos-aquae (IndVal = 0.997, p = 0.001) were significant indicator species of Group B in 2018 and 2019, respectively. D. crassa was an indicator of the combination of Group A and B in 2019 (IndVal = 1.00, p = 0.001). The cyanophyte species Woronichinia compacta was a common Group C indicator species in both years (IndVal > 0.960, p = 0.001). Additionally, non-diazotrophic cyanophytes P. agardhii (IndVal = 0.990, p = 0.001) and Merismopedia tenuissima were Group C indicators (IndVal = 0.834, p = 0.017) in 2019 while P. agardhii was an indicator of the Group B and C combination in 2018 (IndVal = 1.00, p = 0.009). A complete list of indicator species by group and group combination and their indicator values can be found in Table S2.†
The seasonal patterns apparent in both years indicate a shift from a community dominated by diazotrophic species (i.e., Group B) to non-diazotrophic species (i.e., Group C) composed largely of Dolichospermum spp. and Planktothrix spp. (Fig. 3a and b and S1†). N-fixation during the diazotrophic period was demonstrated by the presence of heterocytes belonging to diazotrophic taxa (sensu Boyer 2021) including D. crassa, D. flos-aquae, and D. solitaria in 2018 (Fig. 3c) and D. crassa, D. flos-aquae and Aphanizomenon flos-aquae in 2019 (Fig. 3d). The majority of heterocyte abundance was contributed by D. crassa and D. flos-aquae (Fig. 3c and d). Both Dolichospermum spp. and heterocyte biomass was greater in 2019 (max = 11848 and 942 mg m−3, respectively) than 2018 (max = 4400 and 357 mg m−3, respectively), particularly that of D. crassa.
![]() | ||
Fig. 4 Concentrations (ng L−1) of screened intracellular (grey line) and extracellular (black line) cyanobacterial metabolites present in water from the 0.8 m depth of Buffalo Pound Lake in 2018 (a) and 2019 (c) including individual metabolites and total microcystins (“Total MC”; b, d). Only those metabolites that were present above the limit of detection (see Table S1†) in each year are shown. |
Of the 12 individual metabolites present in 2019 intracellular samples, nine were microcystin congeners (MC-LR, dmMC-LR, MC-HiIR, MC-LA, MC-LF, MC-LW, MC-LY, MC-RR, and MC-YR), and the remainder included anabaenopeptin A and B, and anatoxin-A (Fig. 4c). Only MC-LR and dmMC-LR were present in extracellular samples, both below 50 ng L−1. As in 2018, the most prevalent microcystin congener in intracellular samples was MC-LR (from <10 to 7664 ng L−1); MC-LR and the other microcystin congeners present peaked in mid-August (DOY 227), approximately one month later than in 2018 but again corresponding with Group B. Anatoxin-A was present in very small concentrations (below its LOQ, max = 6 ng L−1) in September. Both anabaenopeptin A and B (max = 459 and 303 ng L−1, respectively) peaked in September on DOY 260. Like in 2018, anabaenopeptin B concentrations also appeared to peak twice, the first smaller peak (202 ng L−1) co-occurring with the microcystin peak and the second (303 ng L−1) occurring on DOY 260. However, anabaenopeptin B was present at far lower concentrations in 2019 than in the previous year.
![]() | ||
Fig. 5 Non-metric multidimensional scaling plots of cyanobacterial phytoplankton community composition throughout the 2018 (stress = 0.091; a, c) and 2019 (stress = 0.072; b, d) sampling seasons. Points indicate day of year in 2018 (a, c) and 2019 (b, d). Abbreviated species name codes using the first three letters of genus and species (gen.spe; shown in green) represent cyanobacterial species. Vectors represent toxins (see Table S1† for abbreviations) for which permutation tests were significant (p > 0.05) in 2018 (c) and 2019 (d) with vector length indicating strength of correlation. |
For example, in 2018, significant vectors showed an association between total microcystins (r2 = 0.39, p = 0.005) and several microcystin congeners (dmMC-LR, r2 = 0.30, p = 0.035; and MC-YR, r2 = 0.35, p = 0.016) and the Group B (DOYs 184 to 212) cyanobacterial community dominated by diazotrophic taxa (e.g., D. crassa and D. flos-aquae). The strongest associations were those between anabaenopeptin A (r2 = 0.48, p = 0.002) and anabaenopeptin B (r2 = 0.68, p = 0.001) and the late-season Group C (e.g., DOYs 235 to 260) cyanobacterial community dominated by Planktothrix agardhii.
Similarly, in 2019, microcystins were also correlated with the mid-season (DOYs 186 to 227, or Group B) bloom comprised largely of the N-fixing taxa (Dolichospermum flos-aquae, D. crassa, A. flos-aquae) and which also included Microcystis aeruginosa. Total microcystins (r2 = 0.45, p = 0.014) and four microcystin congeners, MC-LR (r2 = 0.47, p = 0.01), dmMC-LR (r2 = 0.49, p = 0.008), MC-LY (r2 = 0.37, p = 0.039), and MC-YR (r2 = 0.39, p = 0.036) were significantly correlated to the mid-season bloom cyanobacterial community. Anabaenopeptins A (r2 = 0.81, p = 0.001) and B (r2 = 0.63, p = 0.001) were most strongly correlated with the cyanobacterial community present after DOY 246 (Group C) which was largely dominated by P. agardhii, but also included Woronichinia compacta, Rhabdogloea linearis, and Merismopedia tenuissima.
Spearman rank correlations between the biomass of individual cyanobacterial species and the most prevalent intracellular cyanotoxins (total microcystins, anabaenopeptins A and B) in 2018 and 2019 were again suggestive of links between the Group B species and microcystins, and non-diazotrophic Group C species and anabaenopeptins. In 2018, correlations were strongest between total microcystins and diazotrophic D. crassa (rs = 0.68, p < 0.001) as well as D. flos-aquae (rs = 0.57, p = 0.003) and non-diazotrophic M. aeruginosa (rs = 0.52, p = 0.009). Anabaenopeptin A was most strongly correlated with W. compacta biomass (rs = 0.78, p < 0.001) and P. agardhii (rs = 0.49, p = 0.01) but also negatively correlated with D. crassa (rs = −0.53, p = 0.007) and D. flos-aquae (rs = −0.49, p = 0.02). Anabaenopeptin B was most strongly correlated with P. agardhii (rs = 0.76, p < 0.001) and W. compacta (rs = 0.57, p = 0.003).
Of the Spearman rank correlations between the biomass of cyanobacterial species present in 2019 and intracellular toxins, the strongest was between total microcystins and A. flos-aquae (rs = 0.92, p < 0.001). Significant positive correlations also existed between total microcystins and D. crassa (rs = 0.74, p = 0.001), D. flos-aquae (rs = 0.72, p = 0.001), M. aeruginosa (rs = 0.69, p = 0.002), and R. lineare (rs = 0.56, p = 0.02). Anabaenopeptin A was most strongly correlated with P. agardhii (rs = 0.77, p < 0.001) and W. compacta biomass (rs = 0.78, p < 0.001) as well as M. tenuissima (rs = 0.67, p = 0.002). Anabaenopeptin B was correlated with the same species (P. agardhii, rs = 0.65, p = 0.004; W. compacta, rs = 0.64, p = 0.006; M. tenuissima, rs = 0.61, p = 0.008) although these correlations were slightly weaker; however, it should be stated that correlations between individual species biomass and cyanopeptides are not indicative of a causal relationship between a particular species and the production of a cyanopeptide.
When the results of the Spearman rank correlations are taken into consideration alongside the results of the indicator species analysis, some key similarities are apparent. For example, D. crassa was an indicator for Group B in 2018, and the combination Group B + C in 2019, and A. flos-aquae was an indicator of Group B in 2018. Both species were strongly correlated with microcystins suggesting that they serve as important indicators of potentially toxic mid-summer bloom formation. W. compacta was an indicator of Group C in both years and correlated with the anabaenopeptins. Its presence appears to indicate a shift from a primarily microcystin-producing N-fixing bloom to an anabaenopeptin-producing bloom comprised largely of filamentous Planktothrix.
Our results showed that phytoplankton community structure near the Buffalo Pound Water Treatment Plant's raw water intake changed seasonally but did so in a consistent pattern during the two study years. Cryptophytes and chrysophytes (i.e., Group A) were succeeded by N-fixing cyanobacteria, especially Dolichospermum spp. (i.e., Group B). By September much of the community was composed of non-diazotrophic Planktothrix agardhii (i.e., Group C). Though this general pattern was consistent between years, the timing of succession between groups differed as did their abundance. This succession from non-cyanobacterial species to a diazotroph-dominated community in the warm summer months is well established for many temperate eutrophic lakes.52,53 In the prairie region, diazotrophic blooms have received much research attention by limnologists11,54,55 given highly productive prairie systems tend towards N-limitation, a demonstrated driver of diazotrophic abundance.56 The secondary late-season Planktothrix bloom exhibited in our study has received comparatively little research in prairie systems. We noted active N-fixation verified by the presence of heterocytes during the diazotrophic bloom phase but the mechanism leading to diazotrophic collapse and replacement with Planktothrix remains unclear given low inorganic N (e.g., median September ammonium = ∼30 μg N/L and nitrate = <LOD; Fig. S3†) typically persists late into the ice-free season.
Accompanying the above pattern of species succession was a shift in the cyanometabolites identified at our study site from microcystins in mid-summer to anabaenopeptins in early autumn. Cyanometabolites are detected frequently when cyanobacterial species are present and many species, including those identified here, can produce commonly detected compounds including microcystins and anabaenopeptins – the production of which may co-occur with each other or other cyanometabolites.2,57 For example, the Microcystis, Dolichospermum, and Planktothrix species we identified are known producers of microcystins58 and Microcystis and Planktothrix are also commonly shown to produce anabaenopeptins.57,59 However, organisms that appear morphologically the same may in fact possess unique biosynthetic gene clusters differentiating their toxin-producing abilities;60,61 therefore, causal attribution of cyanotoxin production to any one species is beyond the scope of the work presented here, requiring metatranscriptomics to verify activation of genes responsible for cyanopeptide synthesis.60,62 Instead, we have identified patterns of co-occurrence with indicator species to establish the utility of such indicators for enhanced biomonitoring of bloom-prone systems.
The observed production of microcystins during the mid-season bloom was expected. Increased probabilities of microcystin concentrations exceeding public health thresholds in mid-summer have been described in Buffalo Pound and other prairie lakes previously.12 Indeed, well-known microcystin producers (D. crassa, D. flos-aquae, M. aeruginosa, and in 2018 P. agardhii) occurred in Buffalo Pound in mid-summer, one of which, Dolichospermum crassa, served well as a specific indicator of that bloom phase in both years. Notably, during the mid-summer bloom period, total microcystins in extracellular samples were often higher than the sum of the 12 common individual congeners targeted by our methods – a finding suggestive of the presence of unknown microcystins, of which there are more than 200 known microcystin congeners.63 Another possibility is that transformation products were present containing the ADDA moiety detected by our total microcystin analysis. Such transformation products may contribute to potential toxicity of drinking/recreational waters,58 thus their inclusion with the total microcystin concentrations reported here is warranted to assess risk.
The shift to an anabaenopeptin-producing community in early autumn dominated by Planktothrix agardhii was suggestive of either (1) a change in the metabolite producing species, or (2) some shift in metabolite production owing to environmental factors; however, confirmation of this finding with genomic tools is needed. Detection and subsequent study of anabaenopeptins in surface waters is relatively new;64 however, many of the species we identified in our study, including the abundant Planktothrix agardhii59 and our Group C indicator, Woronichinia compacta,65 may produce anabaenopeptins. In both years, the detected concentrations of intracellular anabaenopeptins in Buffalo Pound increased as Planktothrix biomass increased. Moreover, anti-correlation was observed between Dolichospermum species (D. crassa and D. flos-aquae) biomass and anabaenopeptin A in 2018 while no correlations were observed between Dolichospermum biomass and anabaenopeptins A or B otherwise. While the pathway leading to this observation is unresolved and definitive attribution of anabaenopeptin production to Planktothrix or any other species is beyond the scope of the work presented here, our observation of a large Planktothrix bloom, indicated by W. compacta and accompanied by the presence of anabaenopeptins late in the year is notable given this secondary bloom has been relatively understudied and potential risks associated with anabaenopeptins remain unclear.64
Environmental conditions are likely to play a role in bloom succession from diazotrophs to Planktothrix and should be another focus of future research examining the understudied autumn bloom. Planktothrix agardhii has been shown to thrive and form large blooms under very low to non-detectable inorganic N concentrations66 therefore a shift in nutrient availability is unlikely to explain why Planktothrix formed blooms under low inorganic N conditions in autumn (Fig. S3†). It is more likely the proliferation of Planktothrix is related to its tolerance of seasonal cooling and decreased solar irradiance.59 Day length measured at the nearest weather station (Environment and Climate Change Canada station at Moose Jaw, SK) decreases from 16.35 hours on July 1 (∼DOY 182) to 11.63 hours on October 1 (∼DOY 270).67 The late-season water temperatures in Buffalo Pound of just above 16 °C in early September (DOY ∼245) to 5–8 °C in early October (DOY ∼270, Fig. S4†) are well below the ideal range for Dolichospermum spp.68 but still tolerable for Planktothrix agardhii which has been shown to persist as the dominant species even under ice.69 For example, a similar study of a temperate lagoon by Pilkaitytė et al.70 has shown dominance of N-fixing species (A. flos-aquae and Dolichospermum spp.) during the hot summer months while Planktothrix agardhii and Woronichinia spp. co-occurred in the cool autumn season, often achieving their highest biomass at that time – a finding consistent with the work presented here.
In addition to tolerating low temperatures, Planktothrix spp. are capable of photoheterotrophy and can energize slow filament growth using in-lake carbon sources even at very low irradiance.71 Notably, cyanometabolites, such as microcystins and anabaenopeptins, are rich in C and N and may serve as major nutrient sources to heterotrophic bacteria that occur concurrently with phytoplankton blooms.72,73 These nutrient-rich peptides, as well as other complex molecules produced by cyanobacteria, are taken up and subsequently re-mineralized by heterotrophs – a pathway that could serve as an important labile nutrient source to Planktothrix during photoheterotrophy.72,74Planktothrix may also use these cyanopeptides as a defense mechanism in response to microbial antagonists as illustrated by Rohrlack et al.75 or to displace other bloom-forming cyanobacterial species via niche construction.59,72,76 However, these compounds were largely detected in intracellular samples (e.g., anabaenopeptins) thus their purpose, and the species responsible for their production, is unclear. It is also possible that rapid breakdown or uptake by heterotrophs in surface waters limits extracellular detection of these compounds.
We found that chl-a concentrations correlated well with Dolichospermum but did not capture the later Planktothrix bloom, suggesting reliance on photosynthesis is diminished later in the ice-free season. This finding may support the hypothesis that Planktothrix in Buffalo Pound are photoheterotrophic because less light is required to energize photoheterotrophy than photosynthesis;71,77 thus, likely resulting in reduced chl-a production late in the season when irradiance is low. However, it is also possible that commonly used methods to measure chl-a may not always detect Planktothrix spp. biomass. For example, Wentzky et al.78 reported this lack of association for P. rubescens due to the red pigment phycoerythrin which is detected at different wavelengths (e.g., 525 and 570 nm (ref. 79)) than those commonly used to detect chl-a (e.g., 659, 655, and 750 nm; see methods). Planktothrix in Buffalo Pound was identified as P. agardhii, or “green” Planktothrix, which typically lacks phycoerythrin.80 However, morphology alone may not be sufficient to identify similar species within the same genera.81 Recent work has shown P. agardhii and P. rubescens can have significant genetic overlap, including horizontal gene transfer of the phycoerythrin gene cluster between the two species thus resulting in both red and green phenotypes.82,83 It is conceivable that P. rubescens existed in Buffalo Pound or, more likely, its deeper upstream reservoir previously, therefore an extremely low light adapted “red” phenotype of P. agardhii is a possibility; however algal filaments with red coloration have not been noted in raw intake water at the water treatment plant (Blair Kardash, BPWTP, pers. comm.).
Our results indicate that, under current conditions, the Planktothrix bloom associated with anabaenopeptins poses a potential but unknown risk to human and wildlife health.84Planktothrix was not correlated with toxic microcystins like the earlier Dolichospermum bloom. However, anabaenopeptin toxicity is less studied than microcystin toxicity.72 Yet despite the correlations we have found, without deeper analysis (e.g., at the metatranscriptomic level) we cannot conclude that a particular species has a causal relationship with any individual cyanobacterial metabolite. Besides, toxicogenic strains of P. agardhii have been widely observed elsewhere in temperate Canadian systems (e.g., in Lake Erie85) and environmental drivers may in fact mediate toxin-producing gene activity;59,62 thus, it is possible that toxin expression by the late season Planktothrix bloom in Buffalo Pound may change depending on other factors not explored here (e.g., changing climatic conditions or nutrient sources). As such, continued monitoring and study of this important late season bloom is needed to establish possible toxin risk.
The largest amount of biomass during the sampling period in both years was recorded during the Planktothrix bloom. Given that Planktothrix is both filamentous and can potentially produce toxins, the ability to detect it is critical for drinking water operations. However, our results show that the most cost-effective and easy-to-use chl-a sensors and laboratory methods may not provide operators with reliable indication of increasing Planktothrix biomass as solar irradiance declines later in the year. We found the presence of Woronichinia compacta to be a good indicator of the conditions in which large Planktothrix blooms occurred, but, given the abundance of Planktothrix compared to W. compacta, drinking water operators who routinely use microscopy are likely to find Planktothrix itself easier to target and identify in their weekly samples. Such information can be used to supplement chl-a concentrations and algal cell counts to allow for better preparation for potential water treatment challenges posed by large amounts of buoyant, filamentous biomass. However, identification of indicators may require technicians to go beyond the basic algal identification outlined in standard methods. Given the scale and urgency of the problems posed by cyanobacterial blooms, updated guidance for the identification of cyanobacterial species should therefore be a priority for regulators.
Our finding of common seasonal patterns and indicators each year suggests there may be underlying stability to community structure albeit the time period presented here is limited thus additional years of study would be needed to make such conclusions. We hypothesize that environmental variables are likely driving the differences in magnitude and timing of blooms and their associated cyanotoxins observed in the two years. For example, succession between community groups occurred earlier in 2018 than in 2019 and the cyanobacterial community was, for the most part, Planktothrix-dominated after mid-July. However, in 2019, the bloom phases occurred two weeks to one month later, with the Planktothrix bloom only taking hold at the end of August. Moreover, N-fixation and diazotrophic abundance was also greater in 2019. These findings point to differences in nutrient availability and climatic conditions (e.g., temperature, precipitation, wind, and runoff) potentially acting to modulate bloom formation and collapse. Future work should explore the environmental drivers of the bloom phylogeny illustrated here to better enable the prediction of these distinct bloom phases.
The similarity in patterns observed over two years is certainly useful to understand this complex system but it is not certain that those two years are representative of current or future conditions in Buffalo Pound Lake. With mounting threats to water security in the Saskatchewan River Basin,34 and large-scale irrigation development planned86 directly upstream, understanding and protection of this key drinking water resource in the face of change is increasingly critical. Climate models for the prairie region continue to point to increased future warming87 and there is a growing body of evidence suggesting harmful blooms and their toxins will increase in shallow, polymictic prairie lakes like Buffalo Pound.12 The impacts of this warming may be an increasingly toxic mid-season N-fixing bloom which could exacerbate the challenges for toxin removal. The potential for increased landscape runoff from more severe precipitation events could bring further water treatment challenges related to disinfection by-products if delivery of DOM to the lake increases. Autumn irradiance will remain low regardless, thus with longer ice-free seasons induced by climate warming, filamentous Planktothrix is likely to proliferate in deeper areas of the lake (e.g., near BPWTP intake) later into the year—as suggested by recent reports of P. agardhii blooms (and toxins) persisting through the winter under ice.69 Improved understanding of the autumn Planktothrix bloom, its persistence, and its potential to produce toxins will be increasingly important, thus we recommend targeted study of this key bloom phase, and its potential indicators, in the future.
Footnote |
† Electronic supplementary information (ESI) available. See https://doi.org/10.1039/d2em00078d |
This journal is © The Royal Society of Chemistry 2022 |