Next Article in Journal
Occurrence of Zearalenone and Its Metabolites in the Blood of High-Yielding Dairy Cows at Selected Collection Sites in Various Disease States
Next Article in Special Issue
Anabaenopeptins: What We Know So Far
Previous Article in Journal
The Effects of Mixed Fusarium Mycotoxins at EU-Permitted Feed Levels on Weaned Piglets’ Tissue Lipids
Previous Article in Special Issue
Cyanobacterial Toxins and Peptides in Lake Vegoritis, Greece
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Toxins and Other Bioactive Metabolites in Deep Chlorophyll Layers Containing the Cyanobacteria Planktothrix cf. isothrix in Two Georgian Bay Embayments, Lake Huron

1
Environment and Climate Change Canada, Canada Centre for Inland Waters, 867 Lakeshore Road, Burlington, ON L7S 1A1, Canada
2
Joseph J. Zilber School of Public Health, University of Wisconsin-Milwaukee, Milwaukee, WI 53211, USA
3
Algal Taxonomy and Ecology Inc., P.O. Box 761, Stony Mountain, MB R0C 3A0, Canada
4
School of Graduate Studies, Environmental and Life Sciences, Trent University, Peterborough, ON K9L 0G2, Canada
*
Author to whom correspondence should be addressed.
Submission received: 26 March 2021 / Revised: 21 June 2021 / Accepted: 23 June 2021 / Published: 27 June 2021

Abstract

:
The understanding of deep chlorophyll layers (DCLs) in the Great Lakes—largely reported as a mix of picoplankton and mixotrophic nanoflagellates—is predominantly based on studies of deep (>30 m), offshore locations. Here, we document and characterize nearshore DCLs from two meso-oligotrophic embayments, Twelve Mile Bay (TMB) and South Bay (SB), along eastern Georgian Bay, Lake Huron (Ontario, Canada) in 2014, 2015, and 2018. Both embayments showed the annual formation of DCLs, present as dense, thin, metalimnetic plates dominated by the large, potentially toxic, and bloom-forming cyanobacteria Planktothrix cf. isothrix. The contribution of P. cf. isothrix to the deep-living total biomass (TB) increased as thermal stratification progressed over the ice-free season, reaching 40% in TMB (0.6 mg/L at 9.5 m) and 65% in South Bay (3.5 mg/L at 7.5 m) in 2015. The euphotic zone in each embayment extended down past the mixed layer, into the nutrient-enriched hypoxic hypolimnia, consistent with other studies of similar systems with DCLs. The co-occurrence of the metal-oxidizing bacteria Leptothrix spp. and bactivorous flagellates within the metalimnetic DCLs suggests that the microbial loop plays an important role in recycling nutrients within these layers, particularly phosphate (PO4) and iron (Fe). Samples taken through the water column in both embayments showed measurable concentrations of the cyanobacterial toxins microcystins (max. 0.4 µg/L) and the other bioactive metabolites anabaenopeptins (max. ~7 µg/L) and cyanopeptolins (max. 1 ng/L), along with the corresponding genes (max. in 2018). These oligopeptides are known to act as metabolic inhibitors (e.g., in chemical defence against grazers, parasites) and allow a competitive advantage. In TMB, the 2018 peaks in these oligopeptides and genes coincided with the P. cf. isothrix DCLs, suggesting this species as the main source. Our data indicate that intersecting physicochemical gradients of light and nutrient-enriched hypoxic hypolimnia are key factors in supporting DCLs in TMB and SB. Microbial activity and allelopathy may also influence DCL community structure and function, and require further investigation, particularly related to the dominance of potentially toxigenic species such as P. cf. isothrix.
Key Contribution: Deep chlorophyll layers of Planktothrix cf. isothrix in two meso-oligotrophic embayments in Georgian Bay showed measurable concentrations of the cyanobacterial toxins microcystins (max. 0.4 µg/L) and the other bioactive metabolites anabaenopeptins (max. ~7 µg/L) and cyanopeptolins (max. 1 ng/L)—along with the corresponding genes—that can act as metabolic inhibitors (e.g., in chemical defence against grazers, parasites) and allow a competitive advantage. Supporting physicochemical data provide strongs evidence that opposing and intersecting physicochemical gradients of light and nutrients from a hypoxic hypolimnion are critical in supporting the DCLs in these embayments. Microbial processes and allelopathy may also influence the DCL community structure and function and the dominance of large, toxigenic cyanobacteria such as P. cf. isothrix.

1. Introduction

Deep chlorophyll layers (DCLs) are ecologically important but often overlooked phenomena in many deep lentic ecosystems. DCLs can account for a significant fraction of pelagic primary productivity, play a key role in biogeochemical cycling of nutrients, and influence the vertical movement of micrograzers [1,2,3,4,5,6,7]. DCLs are likely far more prevalent than reported, as most sampling efforts concentrate on surface water layers. Investigations of large and deep Canadian waters such as the Great Lakes, where offshore DCLs have been observed as deep as ca. 50 m, have documented assemblages dominated by diatoms, picoplankton, and cryptophytes, along with the presence of other flagellated autotrophs/mixotrophs (e.g., chrysophytes, dinoflagellates) [8,9,10,11]. Previous studies have identified the importance of light (euphotic depth, Zeu) and thermal stratification in predicting the depth and thickness of deep-living phytoplankton communities, but less is known about the influence of chemical (e.g., nutrient) gradients and even less about the chemical ecology of this deep-living biota [4,7,12,13,14]. Several mechanisms can produce chemical gradients conducive to DCL formation and progression, including upwelling, inter-flow, groundwater, and diffusion from surficial sediments, i.e., internal loading (e.g., [15]). Internal loading often occurs in productive, stratified systems within a hypoxic hypolimnion and can provide a supply of bioavailable dissolved nutrients that are directly accessible to deep-living phytoplankton particularly in shallow waterbodies where Zmix < Zeu [16]. Studies of DCLs in smaller and shallower (<30 m) lakes in Canada report flagellated autotrophs/mixotrophs, e.g., chrysophytes as being the most prominent [1,2,6,17]. However, a more recent study documents DCLs dominated by large cyanobacteria in two nearshore embayments of eastern Georgian Bay (Lake Huron, Ontario, Canada) (e.g., [16]).
Earlier work along the eastern coast of Georgian Bay reported DCLs in a number of sheltered embayments of intermediate size and depth (<30 m), many of which regularly undergo thermal stratification and establish pronounced chemical gradients with depth, [16]. Building on this research, we carried out two nearshore surveys (2014 and 2015) of 15 embayments along this coastline to further investigate the presence of DCLs, and selected two of these waterbodies for more focused work. In 2015, detailed measurements were done in these two embayments to characterize the seasonality and vertical structure of each documented DCL and its phytoplankton composition as well as the coincident physicochemical and nutrient gradients relevant to phytoplankton physiology. We also measured a suite of bioactive metabolites that can play a role in cyanobacterial dominance (the cyanobacterial toxins microcystins, nodularin, anatoxins, saxitoxins, and cylindrospermopsins, and several other metabolic inhibitors known to act at the food web level, anabaenopeptins, cyanopeptolins, and microginins) [13]. Samples from the DCL sites were also obtained in 2018 for additional investigations of these bioactive metabolites, including genetic screening. While the ‘toxins’ have been widely measured, studies of other metabolic inhibitors have been largely limited to Europe with a few recent reports from North America [13,18,19]. To our knowledge, however, none of these other bioactive metabolites have been characterized in DCLs.

2. Results and Discussion

DCLs were documented in 2 of 15 embayments, Twelve Mile Bay (TMB) and South Bay (SB), during a spatial survey conducted in 2014 (Figure 1, Table 1). In 2015, the full survey was repeated in addition to a more detailed seasonal sampling of TMB and SB (Figure 2) where profiles and discrete depth samples were collected to characterize the vertical structure of phytoplankton composition (Table S1 in Supplementary Materials), the coincident physicochemical and nutrient gradients (Figure 3, Table 2), and a suite of bioactive metabolites (Table 3). In 2018, additional samples were collected from the two DCL sites to confirm previous results on bioactive metabolites (Table 3) and construct vertical water column profiles of corresponding genes (Figure 4).

2.1. Spatiotemporal Changes in Phytoplankton Composition in Near-Surface and Deep Chlorophyll Layers in Twelve Mile Bay

Overall, TMB phytoplankton showed seasonality and a near-surface biomass, which was composed of a diverse community of flagellates, diatoms and picoplankton, typical of many meso-oligotrophic zones of the Great Lakes (e.g., [20]). While some of these near-surface taxa were also present in the deeper strata, DCLs were often composed of a distinct community, and dominated by trichomes of the cyanobacteria P. cf. isothrix, which was never observed in near-surface samples. In the summer of 2014, TMB was stratified with a distinct DCL at 10 m depth at the metalimnion, exceeding a total biomass (TB) of 1600 µg/L, almost double that in the epilimnion (Table S1). The DCL was dominated by the buoyancy-regulating cyanobacteria Planktothrix cf. isothrix (>1300 µg/L, 80% TB), with pennate diatoms (7% TB) and cryptophytes (4% TB) also abundant (Table S1). This deep-living community differed markedly from the epilimnetic assemblage, which was dominated by large flagellates—dinoflagellates (13% TB), cryptophytes (24% TB) and colonial chrysophytes (10% TB).
By autumn 2014, lake turnover brought on isothermal conditions and a more uniform vertical distribution of TB (~1000 µg/L) and composition. Both near-surface (1 m) and deep (10 m) samples had a diverse community dominated by large pennate diatoms (Fragilaria crotonensis, up to ~40% TB) and chrysophytes (up to ~40% TB; Dinobryon sertularia, Synura sp., ochromonads). The deep-living population of P. cf. isothrix seen during the summer was reduced significantly (<50 µg/L, 4% TB) and, again, this species was not observed near-surface. We also noted a small population of the nuisance flagellate Gonyostomum semens (<40 µg/L, 4% TB) near-surface. While commonplace in low numbers in softwater lakes, this raphidophyte can produce periodic blooms in boreal lakes (e.g., [21]) and outbreaks have been linked to high Fe availability [22].
During the winter ice cover of 2015, low near-surface biomass (<200 µg/L) was dominated by chrysoflagellates (55% TB), cryptophytes (24% TB) and dinoflagellates (3% TB), with low numbers of cyanobacteria also present (7% TB; Microcystis sp., Planktothrix sp.) (Figure 2, Table S1). Near-bottom (10 m) samples showed a sparse biomass (<40 µg/L), composed of mixotrophic chrysophytes (16% TB) and the large heterotrophic dinoflagellate Gymnodinium helveticum (14% TB), along with filaments of Planktothrix cf. isothrix (10% TB) and the diatom Aulacoseira subarctica (7% TB).
In spring 2015, there was a significant increase in biomass, particularly near-surface (>3300 µg/L TB) (Figure 2, Table S1). Near-surface (1 m) and deep (10 and 11.5 m) samples were all dominated by the large diatom Fragilaria crotonensis (up to 80% TB); large chrysoflagellates (up to 10% TB) were also important while cryptophytes (10% TB) were only present closest to the bottom (11.5 m). Low numbers of the colonial cyanobacteria Aphanothece minutissima and Microcystis sp. were present in the two deeper samples, but Planktothrix was not detected.
Thermal profiles in summer of 2015 revealed a stratified water column, with a mixed/epilimnetic layer (down to ~7 m) overlaying a steep thermocline to the bottom and no discernible hypolimnion (Figure 3). Epilimnetic biomass (>1700 µg/L) was again dominated by diatoms (65% TB; notably F. crotonensis), with large thecate dinoflagellates (13% TB) also being important (Table S1). Cyanobacteria constituted only a small fraction (14%) of this biomass, largely as colonial picocyanobacteria (Aphanocapsa sp., Aphanothece minutissima) and a few trichomes of the diazotrophic cyanobacteria Dolichospermum cf. fuscum. At depth, the phytoplankton community segregated into two distinct DCLs (~1 m apart), each showing strong peaks (>20 µg/L) in chlorophyll-a fluorescence (Figure 3). The more prominent DCL was concentrated within a ~1 m thick layer at the meta-hypolimnetic boundary (~9 m deep). FP mapping estimated that this DCL extended horizontally over an elliptical area of almost 0.2 km2 until the end of the first (~9 m) contour (Figure 3, Figure A1). This biomass peak (>1400 µg/L) was dominated by cryptophytes (50%TB; C. reflexa), diatoms (13% TB; possibly settling from the epilimnion) and P. cf. isothrix (8% TB) (Figure 2, Table S1). While water samples from the smaller DCL at 10 m could not be obtained, FP fluorescence indicated an assemblage dominated by “brown/golden” phytoplankton (e.g., dinoflagellates, diatoms, chrysophytes) (Figure 3).
In the autumn of 2015, the water column showed a deepened mixed/epilimnetic layer and a single metalimnetic DCL at ~9 m. Near-surface biomass had decreased to <700 µg/L and was dominated by large mixotrophic and non-mixotrophic chrysophytes (38% TB) and thecate dinoflagellates (14% TB) (Figure 2, Table S1). Small populations of colonial picocyanobacteria were also present but Planktothrix was not observed. The DCL contained double the near-surface biomass (>1400 µg/L) and was dominated by Dinobryon sp. (70% TB); Peridinium sp. (9% TB) and P. cf. isothrix (8% TB) were also important. A sample taken just above the sediments at 10 m showed a significantly lower total biomass (<700 µg/L) but an increased proportion of P. cf. isothrix (40% TB) (Figure 2, Table S1).
Although P. cf. isothrix was absent near-surface in TMB throughout all sampling, it was present in samples from both the sediment-water interface and surficial sediments throughout the seasons. Sediments—particularly those collected in winter—yielded (on qualitative inspection) viable trichomes of P. cf. isothrix under laboratory conditions simulating the summer metalimnion (Z8 media, <15 °C and <50 µmol m−2 s−1). This suggests overwintering filaments could act as a reservoir of viable cells for the metalimnetic populations of this species. Recruitment of viable, vegetative cyanobacteria cells overwintering in sediments has been observed elsewhere with bloom-forming cyanobacteria, notably Oscillatoria (syn. Planktothrix), Dolichospermum (syn. Anabaena) and Microcystis [23,24,25,26]; however, in these cases, the benthic cyanobacteria were implicated in supporting surface blooms rather than DCLs. In TMB, our observations suggest that the sediments serve as a reservoir of viable P. cf. isothrix cells which seed the DCLs, and to our knowledge remain at depth and do not appear in any detectable numbers at or near the surface. While our evidence is supportive, a more direct measurement of this recruitment process is needed to confirm our hypothesis that sediment seeding may play an important role in the origins, seasonal dynamics, and fate of cyanobacteria-dominated DCLs in TMB.

2.2. Spatiotemporal Changes in Phytoplankton Composition Near-Surface and within Deep-Chlorophyll Layers in South Bay

In summer of 2014, SB was stratified with a distinct metalimnetic DCL at 6.5 m, reaching >60 µg/L chlorophyll-a—six times that measured near-surface (1 m)—of which ~87% was derived from cyanobacteria (based on FP fluorescence). Although a detailed taxonomic identification and enumeration of the summer 2014 DCL was not possible, results from subsequent years (reported below) suggest P. cf. isothrix as a prominent and recurring member of these deep phytoplankton communities in SB. Near-surface (1 m) summer biomass in 2014 was ~1000 µg/L, and dominated by chrysophytes (35% TB), dinoflagellates (10% TB), and picocyanobacteria (7% TB) (Table S1). Large cyanobacteria were present but in very low biomass (~3% TB; notably Dolichospermum).
By autumn 2014, lake turnover brought on isothermal conditions and biomass had concentrated near-surface, increasing over three times to ~3000 µg/L; dominated by large diazotrophic cyanobacteria (Aphanizomenon flos aquae complex, 45% TB; Dolichospermum planctonicum, 4% TB), with Euglena sp. (9% TB), Aulacoseira ambigua (9% TB), and Plagioselmis nanoplanktica (8% TB) also being abundant (Table S1).
During the winter ice cover of 2015, biomass near-surface (>400 µg/L at 1 m) was dominated by large colonial, scaled chrysoflagellates (Synura sp., 77% TB), with small populations (<5% each) of dinoflagellates and picocyanobacteria (Table S1). In contrast, a much lower biomass at ~8 m depth (<200 µg/L) was composed of picocyanobacteria (Synechococcus spp., 45% TB), chrysophytes (37% TB), and dinoflagellates (9% TB).
By spring of 2015, under isothermal conditions, near-surface (1 m) biomass exceeded 1000 µg/L; dominated by a variety of chrysoflagellates (37% TB), cryptophytes (17% TB; mostly C. reflexa) and the large, horned dinoflagellate Ceratium furcoides (15% TB) (Figure 2, Table S1). At 5.5 m depth, a comparable biomass (>1000 µg/L) had an assemblage of scaled chrysophytes (36% TB; e.g., Chrysosphaerella longispina) and large dinoflagellates (Ceratium sp., 21% TB) with a few diazotrophic cyanobacteria (Dolichospermum lemmermanni, Aphanizomenon flos-aquae) also present. Near-bottom (10.5 m) samples had a less biomass (<400 µg/L), composed mostly of Cryptomonas spp. (32% TB), Chlamydomonas sp. (14% TB), and P. cf. isothrix (9% TB).
Profiles in summer of 2015 revealed a stratified water column with near-surface (1 m) biomass (~800 µg/L) dominated by filamentous diatoms (34% TB; Aulacoseira), thecate dinoflagellates (15% TB), and the diazotroph D. planktonicum (12% TB) (Figure 2, Table S1). A distinct and sizable DCL (>3500 µg/L) was concentrated within a ~1 m thick layer at the meta-hypolimnetic boundary (~7 m). FP mapping estimated that the DCL extended horizontally over an elliptical area of almost 0.03 km2, extending over an area bounded by the ~10 m depth contour (Figure 3, Figure A1). This deep biomass peak was dominated by P. cf. isothrix (65%); other cyanobacteria were also present (Aph. flos-aquae complex, 6% TB; Microcystis novacekii, 3% TB) along with dinoflagellates (7% TB) and diatoms (4% TB).
By autumn 2015, lake turnover brought on isothermal conditions and a more uniform vertical distribution of phytoplankton, i.e., no DCL was observed. A similar sized surface biomass (~600 µg/L) showed a shift in community composition from the summer assemblage towards an increased proportion of diazotrophic cyanobacteria (D. planktonicum, 32% TB; Aph. flos-aquae, 6% TB) and a decrease in the filamentous diatom Aulacoseira (16% TB) (Figure 2, Table S1). In the near-bottom sample (~9 m), a comparable biomass (~500 µg/L) showed similar composition (Aulacoseira, 17% TB; D. planktonicum, 15% TB). As seen in summer, P. cf. isothrix was again present at depth but at comparably lower proportion (9%), but other cyanobacteria (Pseudanabaena sp. 8% TB; Aph. flos-aquae complex, 7% TB) were also observed.
We note that while numerous other studies have reported DCLs that have been mixed into the epilimnion by physical processes such as turnover, storms, or artificial aeration/mixing and subsequently manifested as ephemeral surface blooms [4,18,27,28,29,30,31,32], there have been no reports of surface blooms in TMB nor SB near the locations of the DCLs, despite the high density of cottagers in the area. This suggests that the thick metalimnetic plates of Planktothrix remain essentially segregated from the surface community and are seeded and sustained by the physicochemical conditions in the bottom layers and sediments.

2.3. Light and Pigmentation

Light is a major factor which often controls phytoplankton vertical structure and the formation of deep-living communities (e.g., [7]), and many taxa modify or supplement their effective light regimes via vertical migration, mixotrophy, and photoadaptation. The abundance of large flagellated mixotrophs in the deep layers in TMB and SB (Table S1) suggests that motility and bactivory may provide access to alternative supplies of energy and nutrients (e.g., [17]); however, these photophagotrophs rarely predominated the DCLs in these embayments. PAR measures showed that the DCLs were positioned within the euphotic zone (Table 1), and we therefore conclude that in situ light levels were sufficient to support the development of P. cf. isothrix, i.e., light availability was not a primary factor limiting these DCL communities (e.g., [10,33]). Although current literature suggests phycoerythrin-rich P. rubescens is more often the predominant Planktothrix species in deep-living maxima, the predominant species P. cf. isothrix remained green-pigmented, suggesting minimal cellular phycoerythrin (Figure 2, Figure A2). We are aware of only one other published report of a DCL dominated by green-pigmented Planktothrix [3] although broader geographical surveys are needed [9,14,34,35,36,37,38].

2.4. Nutrients and Physicochemical Conditions

In SB, DCLs dominated by P. cf. isothrix were located at the metalimnion where sufficient PAR intersected with elevated levels of major nutrients and trace metals, particularly PO4, dissolved inorganic nitrogen (DIN) and dissolved Fe (Table 1, Table 2, Figure 3). Vertical gradients of these chemical constituents increased with depth through the hypolimnion, indicative of internal loading from surficial sediments [16,39,40,41]. PO4, Fe, NH4, Mn and Co concentrations were up to three orders of magnitude higher near the sediments as conditions became hypoxic (dissolved oxygen <2 mg/L). Compared to TMB, bottom concentrations were generally higher and gradients were more pronounced in SB—which also exhibited a lower redox potential in the hypolimnion, particularly near sediments (−120 mV vs. −65 mV in TMB). Although similar gradients with depth through the hypolimnion were also evident in TMB for most measures taken, they were not as pronounced as in SB and this, along with markedly lower PO4, is consistent with the significant difference in DCL biomass between the two embayments (Figure 2 and Figure 3). Interestingly, the sheathed bacteria Leptothrix were observed in the DCL of TMB but not SB (Figure A2). These chemoorganotrophs exploit DOC-rich environments in the presence of Fe2+ and Mn2+, which they precipitate as ferrous hydroxide and manganese oxide at the interface between oxic/anoxic zones, often co-precipitating PO4 and making it inaccessible to phytoplankton (e.g., [42,43]). Their PO4-precipitating activity in TMB may, to some extent, account for the minimal change in PO4 with depth despite conditions conducive to internal loading but further study to quantify the significance of this mechanism is needed. Differences in sediment chemistry (e.g., sulphate/sulphide or aluminum) could also account for the contrasting PO4 profiles but this has not yet been characterized nor quantified in the two embayments [39,40,41].
There is a considerable body of evidence to indicate that the development of opposing and intersecting vertical gradients of light and nutrients is strongly correlated with DCL formation in stratified, meso-oligotrophic systems [4,5,7,15,34,44]. Our survey showed that stratified Georgian Bay embayments lacking these opposing and intersecting vertical gradients also lacked discernible DCLs. Some, such as Sturgeon Bay, had a hypoxic hypolimnion and associated gradient in nutrients that reached the metalimnion but did not have sufficient light penetration to support a DCL, i.e., Zeu << Zmix (Table 1). Others, such as Cognashene Lake (actually an embayment), had sufficient light at the metalimnion (Table 1) and elevated nutrients in the hypoxic hypolimnion (e.g., >1400 µg/L Fe, >600 µg/L Mn, >200 µg/L NH4) but it was unclear if these opposing gradients intersected as we lacked sufficient resolution with depth. We speculate that the lake’s greater depth and thicker hypolimnion rendered the hypolimnetic nutrients inaccessible to phytoplankton at sufficient PAR (Table 1); however, with more stable and prolonged periods of stratification, and associated hypoxia, as expected with climate change, the nutrient gradient would presumably strengthen.
Consistent with our observations, others have observed metalimnetic DCLs most commonly occuring at depths of less than 15 m, particularly those dominated by bloom-forming cyanobacteria [4,5,7,31,34,38,44,45,46]. Although significant DCLs dominated by potentially toxic cyanobacteria such as Planktothrix have been reported at depths below 15 m (e.g., [47]), this appears relatively rare [7,46]. DCLs existing deeper than 15 m (close to 30 m depth) are generally dominated by picoplankton, flagellated autotrophs/mixotrophs (e.g., large colonial chrysophytes, dinoflagellates), and diatoms [9,10,11,15]. The occurrence of these deeper DCLs has been documented in much larger (>19,000 km2) and deeper systems (>60 m) including offshore Lake Ontario [10,11], Lake Michigan, Lake Superior, and Lake Huron [8,9] and in marine systems (e.g., [15]). The significant depth of these systems and corresponding thickness of the hypolimnion suggests upwelling of bottom waters [15] and recycling from biomass [1,6,48] as more likely sources of nutrients contributing to the DCLs—this, in contrast to the direct access metalimnetic DCLs in the relatively shallower TMB and SB have to the nutrient gradient generated from internal loading at the sediment-water interface. Although opposing and intersecting physicochemical gradients in nutrients and light are considered key drivers of DCL biomass formation in general; their effects on phytoplankton community composition and activity need further investigation.

2.5. Toxins Produced by Cyanobacteria

Over the course of sampling (2014, 2015, and 2018), near-surface summer samples from 5 of 15 Georgian Bay embayments showed both detectable levels of microcystins and the genetic potential to produce these hepatotoxins (positive for mcy gene). The highest concentration observed in any sample (0.4 µg L−1 in 2018) was well below the Canadian guidelines for safe drinking and recreational water (<1.5 µg/L; <20 µg/L respectively; Health Canada 2017) (Table 3, Figure 4). MC-LA was the most frequently observed microcystin variant and also present at higher concentrations than MC-LR and MC-RR. This is consistent with other reports suggesting MC-LA is more prevalent in systems with lower trophic status (e.g., [13,49]). In both TMB and SB, both the microcystin gene (mcy) and the microcystin toxins were present at the metalimnetic DCLs and in the hypolimnion. Based on our 2018 measurements, the Planktothrix-dominated DCL in TMB coincided with peaks in microcystin concentration and corresponding gene copies (mcy) (Figure 4, Table S1) suggesting these cyanobacteria as the primary source. Other species of Planktothrix are known for the production of a wide range of bioactive metabolites including microcystins [18]; however, production by P. cf. isothrix has yet to be directly confirmed.
The presence of the mcy gene and trace amounts of microcystin in the hypolimnion of Tadenac Bay is also noteworthy (Table 3) as this embayment is relatively unimpacted (TP < 10 µg/L, TDP < 5 µg/L) with minimal shoreline development. The presence and low-level expression of the mcy gene in Tadenac Bay suggests an innate potential for toxin production despite the oligotrophic nature of this system (Table 3). None of the other toxins for which the samples were screened were detected, i.e., anatoxins, saxitoxins, cylindrospermopsins, nodularins (Table 3). In North Bay, a single sample showed the presence of low copy numbers of the saxitoxin gene (sxtA), but not the toxin.

2.6. Other Bioactive Metabolites Produced by Cyanobacteria

When associated with cyanobacteria, the term ‘toxins’ is generally applied to the small fraction of bioactive metabolites that affect humans and other large vertebrates; however, many bioactive metabolites produced by cyanobacteria have no known effect on human, pet, or livestock health, but are active towards more proximal elements of the food web including competing phytoplankton, bacteria, pathogens (e.g., viruses, chytrids), and grazers (e.g., [13]). These bioactive metabolites, which include a variety of peptides, have been reported previously from some European and US lakes and are known to inhibit metabolic processes (e.g., protease inhibitors) and function in chemical defence [13,19]. Our study represents a first-time analysis of samples from Canadian waters for both cyanobacterial toxins and these other bioactive metabolites. As seen with microcystins, several anabaenopeptins and one cyanopeptolin (1007) as well as their corresponding genes (apn, oci) were detected in both TMB and SB (Table 3); however, clear peaks coincident with the Planktothrix-dominated DCL were only evident in TMB (Figure 4). Anabaenopeptins (A, B, and F combined) peaked in 2018 at 6.6 µg L−1 in the DCL of the meso-oligotrophic TMB, which was comparable to concentrations found in eutrophic surface waters in the United States and Europe (e.g., [51,52,53]). Anabaenopeptins are thought to play a role in chemical defence via their ability to inhibit proteases used by parasites of cyanobacteria, such as chytrids, to digest their hosts [54]. This anti-parasitic activity could afford cyanobacteria in TMB an additional competitive advantage to establish dominance within the DCL.
Our results indicate that these (and potentially other) bioactive metabolites—typically reported from eutrophic systems (e.g., [51])—are also produced in meso-oligotrophic waters. Furthermore, while most studies have concentrated efforts on surface waters where visible blooms of cyanobacteria are most evident, our data demonstrate the potential importance of these semiochemicals in the often-overlooked communities in DCLs.

3. Conclusions

We observed the annual formation of persistent DCLs in TMB and SB during the stratification period, with a predominance of the large cyanobacteria Planktothrix cf isothrix. These deep-living layers were formed near the hypoxic hypolimnion at the intersection of opposing vertical gradients of PAR and nutrients, notably Fe, Mn and PO4. Similar DCLs were not observed in the other (13) embayments surveyed along the same coastline, where analogous gradients of PAR and water chemistry were not detected. We also observed the Fe- and PO4-precipitating bacteria Leptothrix within the DCL of TMB but not SB, suggesting that microbial processes may contribute to differences in their metalimnetic chemical gradients and DCL composition and biomass. In both embayments, the Planktothrix-dominated DCLs coincided with measurable levels of cyanobacterial toxins and other bioactive metabolites (most notably microcystins and anabaenopeptins), along with the associated genes. This is the first report of these compounds in DCLs from meso-oligotrophic Canadian lakes, and it merits more focused work to understand if and how these compounds function in the establishment and maintenance of these deep-living communities. The presence of significant numbers of live Planktothrix in surficial sediments and overlying water suggest these act as important seed populations to the DCLs, meriting further investigation. Overall, the prevalence of these significant DCLs dominated by large, potentially harmful cyanobacteria across different waterbodies is unknown, and their formation in apparently oligotrophic systems requires more research to understand the driving factors and track their change with changing environmental conditions and anthropogenic development in respective regions.

4. Materials and Methods

4.1. Study Sites

Fifteen embayments along the eastern shore of Georgian Bay, Lake Huron (Ontario, Canada) were surveyed, with a focus on sheltered embayments of intermediate size and depth (<30 m), particularly those that are reported to regularly undergo stratification and establish pronounced chemical gradients facilitated by hypolimnetic hypoxia (Figure 1) [16,55]. All 15 embayments were sampled during summer (July/August) and autumn (September/October) stratification of 2014 and 2015 (Figure 1 and Figure A1, Table A1). Twelve Mile Bay (TMB) and South Bay (SB) were selected, based on the presence of DCLs, for additional and more detailed sampling and analysis in 2015 in winter (February/March), spring (June), summer (August), and autumn (September/October). Based on the earlier results, TMB and SB were sampled again during stratification in the summer of 2018 (July/August) and analyzed for a broader suite of bioactive metabolites beyond the cyanobacterial toxins commonly measured (Table 3).

4.2. Physicochemical Profiles of Water Column

At each sampling, photosynthetically active radiation (PAR; 400–700 nm) was measured with a LI-193 spherical underwater quantum sensor (LI-COR Biosciences, NE, USA) at 0.5 m depth intervals. Measurements were corrected for variance in incident irradiance using a LI-190R quantum sensor (LI-COR Biosciences, NE, USA), and euphotic zone depth calculated (i.e., <1% sub-surface irradiance). A YSI multi-parameter sonde (Xylem Inc., New York, NY, USA) was used to obtain depth profiles of temperature, dissolved oxygen, pH, specific conductivity, redox, and turbidity at each site, while a FluoroProbe (FP) (bbe Moldaenke GmbH, Schwentinental, Germany) was used to measure depth profiles of in situ, fluorescence-based, phytoplankton pigment class-specific chlorophyll (‘green algae’, ‘cyanobacteria’, ’brown’ algae’ (diatoms, chrysophytes and dinoflagellates), and ‘cryptophytes’).

4.3. Water Sampling

At each site and date, whole water samples were collected from near-surface (1 m) and bottom (1 m from sediments) using a 10 L Niskin water sampler, subsampled into acid washed polyethylene bottles and kept in the dark at 4 °C until processed (within 24 h). During stratification, samples were also collected from the middle of the thermocline (metalimnion), the epilimnion (1 m), and hypolimnion (1 m from sediments). Samples were also taken at additional depths where DCLs were detected in the FP profiles.

4.4. Sediment Sampling

Surficial sediments were sampled using a modified gravity corer (Uwitech, Austria). Overlying water was siphoned and surficial 1.0 cm of the core was extruded on-site, placed in pre-labeled Whirl-Pak® bags, and immediately transported to the laboratory on ice in a dark cooler for microscopic imaging before being frozen in the dark at −20 °C.

4.5. Water Quality Analysis

Dissolved and particulate components were analyzed using filtrate from a cellulose acetate filter (0.45 μm pore size, 47 mm diameter) and material collected on a Whatman GF/C filter (1.2 μm nominal pore size, 47 mm diameter), respectively. These included major nutrients (dissolved inorganic and organic carbon (DIC, DOC), particulate organic carbon (POC), nitrate/nitrite (NO2/3), ammonium (NH4), total dissolved Kjeldahl nitrogen (TKN), particulate organic nitrogen (PON), phosphorus (P), dissolved phosphorus (DP), phosphate (PO4)), dissolved silica, and extracted chlorophyll-a, all of which were analyzed at the National Laboratory for Environmental Testing (NLET, Burlington, Ontario) using standard methods [56]. Dissolved organic nitrogen (DON) was calculated by subtracting NH41+ from TKN.

4.6. Extraction and Analysis of Cyanobacterial Bioactive Metabolites

Whole water samples were concentrated onto Whatman GF/C filters (1.2 µm nominal pore size, 47 mm diameter) and stored at −80 °C in the dark until extraction. Cyanobacterial toxins and other bioactive metabolites were extracted in 10 mL of analytical grade aqueous methanol (1:1 v/v) amended with analytical grade formic acid (0.1%) using probe sonication (three 30 sec pulses) (Fisher Scientific Co. Qsonica Sonicator Q500, Ontario, Canada). Samples were then centrifuged at 3000 rpm for 15 min to pellet debris. A PTFE syringe filter (1.0 µm pore size, 30 mm diameter) attached to a 10 mL gas-tight glass syringe was used to filter the resulting supernatant into a glass vial, which was then evaporated to dryness using nitrogen gas-flow and heat (30 °C). The resulting residue was reconstituted with 1 mL of analytical grade aqueous methanol (1:1 v/v) and vortexed. A final filtration was done using a PTFE syringe filter (0.45 µm pore size, 15 mm diameter) attached to a gas-tight glass syringe into a 1.5 mL HPLC amber glass vial and stored at −80 °C in the dark until analysis.
Nineteen cyanobacterial peptides and five alkaloids were analyzed by high performance liquid chromatography tandem mass spectrometry (HPLC-MS/MS), as described previously [51]. Certified reference standards of microcystin-LR, -[Dha7]LR and nodularin (NOD-R) were purchased from the National Research Council (NRC) of Canada Biotoxins program (Nova Scotia, Canada). Microcystin-RR, -LA, -LF, -YR, -WR, -LY, -LW, -HtyR, and -HilR were purchased from Enzo Life Sciences (New York, NY, USA). Anabaenopeptin-A, -B, and -F; cyanopeptolin-1007, -1021, and -1040; as well as microginin-690 were purchased from MARBIONC (NC, USA). Anatoxin-a fumarate was purchased from Tocris Bioscience (MN, USA) as a racemic mixture. Homo-anatoxin-a was purchased from Abraxis. Cylindrospermopsin was purchased from Enzo Life Sciences. Saxitoxin and neosaxtoxin were purchased from the NRC. Each analyte exceeded 95% purity as per the manufacturer’s certification.
HPLC-MS/MS analysis was performed on a Sciex 4000 QTRAP (AB Sciex, USA) tandem mass spectrometer equipped with a Shimadzu Prominence HPLC. Peptides were chromatographically separated in 20 µL injections of extracts on a Luna C8 column (Phenomenex, CA, USA) using gradient elution where the mobile phase consisted of A (0.1% formic acid and 5 mM ammonium acetate in HPLC grade water) and B (0.1% formic acid and 5 mM ammonium acetate in 95% acetonitrile). The gradient began at 30% B for 3 min, increasing over a linear gradient to 95% B at 9 min, and held at 95% B until 15 min at which point B was returned to the starting condition for 5 min. The mass spectrometer was operated in positive mode using a scheduled multiple reaction monitoring method. Alkaloids were separated by hydrophilic interaction liquid chromatography (HILIC) (SeQuant®, 5 μm, 150 × 2.1 mm I.D., EMD Millipore Corporation, MA, USA) with mobile phases of (A) HPLC water with 60 mM formic acid and (B) 100% acetonitrile. Isocratic elution (60% B) was used for the HILIC method. The concentration of target analytes was determined based on a linear regression model of peak area relative to the known concentration of an 8-point calibration curve prepared in 1:1 methanol:water. Regression coefficients of R > 0.98 were accepted at an accuracy of >90% at each calibration level.

4.7. Taxonomic Identification and Enumeration

Unfiltered sample aliquots of 100 mL were preserved in Lugol’s iodine solution (2% v/v) for later taxonomic identification, abundance, and biomass using the Utermöhl technique [57].

4.8. DNA Extraction and PCR Amplification

Water samples were filtered onto 0.22 µm pore size, 47 mm polycarbonate filters and stored at −80 °C for DNA extraction. Samples were extracted using DNeasy PowerWater kit (Qiagen, Canada) following the manufacturer’s protocol.
The availability of a commercial multiplex, quantitative PCR (qPCR) kit for genes of cylindrospermopsin (cyrA), microcystin/nodularin (mcyE/ndA), and saxitoxin (sxtA) facilitated the simultaneous measurement of respective gene copy numbers. The qPCR assays were performed using the PhytoxigeneTM kit (Diagnostic Technology, Australia) according to the manufacturer’s protocol. Briefly, the lyophilized master mix was first spun down and then reconstituted in PCR-grade water. Each qPCR reaction consisted of 20 µL of the reconstituted master mix and 5 µL of template DNA (10–50 ng) and was carried out in a Bio-Rad CFX96 cycler (Bio-Rad, USA). The cycling conditions consisted of an initial denaturation at 95 °C for 2 min, followed by 40 cycles of denaturation at 95 °C for 15 sec and annealing-extension at 60 °C for 30 s. Gene copies per sample were calculated using a standard curve (target gene copy number vs. Ct) determined for each target gene. Standard curves for all target genes were constructed using standards purchased from the same manufacturer (correlation coefficient and efficiency: mcyE/ndA: R2 = 0.999; E = 100.7%; cyrA: R2 = 0.999; E = 100.7%; sxtA: R2 = 1.000; E = 102.1%).
All samples were screened by conventional endpoint PCR for the presence of anatoxin-a synthesis gene and several synthesis genes of oligopeptides know to be bioactive metabolites: aeruginoside, anabaenopeptin, cyanopeptolin, microcystin, microginin, and prenylagaramide. PCR was performed in a total volume of 25 µL composed of 5 µL of GoTaq 10X buffer (Promega), 0.5 µL of dNTPs (10 mM each), 0.5 µL each forward and reverse primers (10 µM), 0.125 µL GoTaq Polymerase (5 units µL−1), 2 µL of template DNA (10–50 ng), and 16 µL of nuclease-free water (see Table A2 for primer sequences). The thermal cycling condition consisted of an initial denaturation step at 95° C for 3 min followed by 35 cycles of denaturation at 95° C for 1 min, annealing for 30 sec (variable annealing temperature; see Table A2), and extension at 72° for 30 s. PCR products were evaluated on a 2% agarose gel.
Anabaenopeptin gene copies were quantified using digital PCR (dPCR). Each dPCR mixture consisted of 7.5 µL of QuantStudio 3D digital PCR master mix v.2 (Thermo Fisher Scientific, Canada), 1.4 µL of each 900 nM forward and reverse primer, 0.75 µL of 250 nM probe, 2 µL of DNA extract, and 195 µL nuclease-free water. The reaction mixtures were loaded into a QuantStudio 3D digital PCR 20K chip v2 (Thermo Fisher Scientific, Canada) using a QuantStudio 3D chip loader (Thermo Fisher Scientific, Canada). The reactions were carried out in a ProFlex thermocycler (Thermo Fisher Scientific, Canada) using the following cycling conditions: 96 °C for 10 min, 40 cycles of 60 °C for 2 min and 98 °C for 30 s, and a final step of 60 °C for 2 min. The chips were read on a QuantStudio 3D digital PCR instrument and the results analyzed using the QuantStudio AnalysisSuite software (V 3.2, 2019) (Thermo Fisher Scientific, Canada).

Supplementary Materials

The following are available at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/toxins13070445/s1, Table S1: Taxonomic identification and enumeration of phytoplankton in Twelve Mile Bay and South Bay embayments of Georgian Bay, Lake Huron (Ontario, Canada).

Author Contributions

Conceptualization, A.Z. and S.B.W.; methodology, T.R.M., L.C.W. and H.K.; validation, A.Z.; formal analysis, A.Z. and S.B.W.; investigation, A.Z. and S.B.W.; resources, A.Z. and S.B.W.; data curation, A.Z., S.B.W., T.R.M., L.C.W. and H.K.; writing—original draft preparation, A.Z.; writing—review and editing, A.Z., S.B.W., T.R.M., L.C.W. and H.K; visualization, A.Z.; supervision, A.Z. and S.B.W.; project administration, A.Z. and S.B.W.; funding acquisition, A.Z. and S.B.W. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by Environment and Climate Change Canada’s Lake Simcoe/South-eastern Georgian Bay Clean-Up Fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or supplementary material.

Acknowledgments

We would like to extend our sincerest gratitude to the hardworking students and staff, Kaitlin Burek, Shelby Grassick, Flora (Zi Wan) Dong, Jay Guo, Joseph Gzik, Anqi Liang, and the skilled technical team of the Research Support Section, Adam Morden, Daniel Abbey, and Corey Treen, within the Watershed Hydrology and Ecology Research Division at Environment and Climate Change Canada for their assistance in the field and laboratory. Thank you also to Patricia Chow-Fraser and Stuart Campbell at McMaster University as well as Aisha Chiandet at the Severn Sound Environmental Association for their support in the field as well as for their knowledge, guidance, and navigation in planning and execution of our science in this complex watershed. A thank you also to Bob Rowsell, David Depew and Jacqui Milne for their support, encouragement, and engagement in the early stages of the work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Location within (A) Twelve Mile Bay (TMB), scale 52.16 m and (B) South Bay (SB), scale 0.105 km, in Georgian Bay, Lake Huron, Ontario, Canada where each of the deep-chlorophyll layers were sampled, denoted by a filled, red circle. Three-dimensional mapping by Fluoroprobe© revealed that each DCL extended horizontally to the end of the first contour demarcated on the map (~9 m in TMB and ~10 m in SB, approximate extent indicated by red lines) and therefore each of these contours were used to estimate the elliptical area over which each of the DCLs existed. Maps reproduced from Navionics Web API v2, www.navionics.com accessed on January 2021.
Figure A1. Location within (A) Twelve Mile Bay (TMB), scale 52.16 m and (B) South Bay (SB), scale 0.105 km, in Georgian Bay, Lake Huron, Ontario, Canada where each of the deep-chlorophyll layers were sampled, denoted by a filled, red circle. Three-dimensional mapping by Fluoroprobe© revealed that each DCL extended horizontally to the end of the first contour demarcated on the map (~9 m in TMB and ~10 m in SB, approximate extent indicated by red lines) and therefore each of these contours were used to estimate the elliptical area over which each of the DCLs existed. Maps reproduced from Navionics Web API v2, www.navionics.com accessed on January 2021.
Toxins 13 00445 g0a1
Table A1. Coordinates and depths of sites under study in Georgian Bay, Lake Huron (Ontario, Canada). All 15 embayments were sampled during summer (July/August) and autumn (September/October) stratification of 2014 and 2015. Twelve Mile Bay and South Bay were selected (based on presence of DCLs) for additional and more detailed sampling and analysis in 2015 in winter (February/March), spring (June), summer (August), and autumn (September/October). TMB and SB were again sampled during stratification in the summer of 2018 (July/August) and analyzed for a broader suite of bioactive metabolites beyond the cyanobacterial toxins commonly measured (Table 3).
Table A1. Coordinates and depths of sites under study in Georgian Bay, Lake Huron (Ontario, Canada). All 15 embayments were sampled during summer (July/August) and autumn (September/October) stratification of 2014 and 2015. Twelve Mile Bay and South Bay were selected (based on presence of DCLs) for additional and more detailed sampling and analysis in 2015 in winter (February/March), spring (June), summer (August), and autumn (September/October). TMB and SB were again sampled during stratification in the summer of 2018 (July/August) and analyzed for a broader suite of bioactive metabolites beyond the cyanobacterial toxins commonly measured (Table 3).
SiteLatitudeLongitudeDepth (m)
Byng Inlet45.77056−80.568363
Sturgeon Bay North45.6133−80.432514
Sturgeon Bay South45.6055−80.40966
Deep Bay45.3953−80.223615
Woods Bay45.13785−79.98948
Twelve Mile Bay45.0838−79.94612
Tadenac Bay45.0588−79.977524
Go Home Bay44.9913−79.937610
Longuissa Bay44.9613−79.88848
Cognashene Lake44.9510−79.918616
Musqwash Bay West44.9520−79.878730
Musqwash Bay East44.9483−79.8500620
North Bay44.8919−79.792518
South Bay44.8762−79.785710
Honey Harbour44.8770−79.82637
Figure A2. Microscopic imaging of water sample at the DCL (~9 m) in Twelve Mile Bay showing Planktothrix cf. isothrix (black arrow) and Leptothrix with Fe oxide precipitates (five red arrows). Images are not to scale and are only for representative purpose.
Figure A2. Microscopic imaging of water sample at the DCL (~9 m) in Twelve Mile Bay showing Planktothrix cf. isothrix (black arrow) and Leptothrix with Fe oxide precipitates (five red arrows). Images are not to scale and are only for representative purpose.
Toxins 13 00445 g0a2
Table A2. Primers used to amplify regions of genes coding for cyanobacterial toxins and other bioactive metabolites (oligopeptides).
Table A2. Primers used to amplify regions of genes coding for cyanobacterial toxins and other bioactive metabolites (oligopeptides).
Gene LociPrimer 1Sequence (5′-3′) 2Amplicon Length (bp)Annealing
Temperature (°C)
Aeruginoside
aerDaerD-FGAAACCAGTAGTGAACAGACCTTAAATTATC49360
aerD-RGACCAACTCATCTCAATTTTCCC
Anabaenopeptin
apnDapnDTe-FGACACGCCTTCTATTTTCGAGA58160
apnDTe-RCGCGAAATAAAACAATAGGGG
apnCapnNMT-FCGTGCAGATGATGACCTATCCA47055
apnNMT-RAAGGTTCGCAATACTTCAGGGTT
Anatoxin 3
anaCanaC-gen-F2TCTGGTATTCAGTMCCCTCYAT36658
anaC-gen-R2CCCAATARCCTGTCATCAA
Cyanopeptolin
ociCcptDTe-FGATCTCTATCAACAGTTTGGAGCAA69061
cptDTe-RACTGTTCGGCTAACACTTGAACAT
ociBociB-FTGGTTTTTAGATCAATTTGAGTCCG51266
ociB-RCCACTGTTTTTGCCAAAGAGTG
Microcystin
mcyCmcyCTe-FTTACAAGCGATGAATCTCATGG50360
mcyCTe-RGGGATTTAATAAGAAACCATCAACC
mcyEHep-FTTTGGGGTTAACTTTTTTGGGCATAGTC46660
Hep-RAATTCTTGAGGCTGTAAATCGGGTTT
Microginin
micDmgnTe-FTGGTCAATGGGAGGAGTGATAG51460
mgnTe-RCTGTAGTGATCTCCACTAATCCATTG
micAmcnA-FAAACCCTTTGATTTGAGCCAG42860
mcnA-RCAGCAAGGTACAGCCCTGTT
Prenylagaramide
pagAacyA-FCCCTGGAAAAGATATTTTAGGAGC48960
acyA-RATTTGAAGTACGGTTTGACATGG
1 Primers for cylindrospermopsin, nodularin, and saxitoxin are not listed here since these toxins were only measured using the PhytoxigeneTM kit (Diagnostic Technology, Australia) and the primers are included in the kit’s master mix. 2 [18]. 3 [50].

References

  1. Pick, F.R.; Nalewajko, C.; Lean, D.R.S. The origin of a metalimnetic chrysophyte peak. Limnol. Oceanogr. 1984, 29, 125–134. [Google Scholar] [CrossRef]
  2. Fee, E.J. The vertical and seasonal distribution of chlorophyll in lakes of the Experimental Lakes Area, northwestern Ontario: Implications for primary production estimates. Limnol. Oceanogr. 1976, 21, 767–783. [Google Scholar] [CrossRef]
  3. Pannard, A.; Planas, D.; Noac’h, P.; Bormans, M.; Jourdain, M.; Beisner, B. Contribution of the deep chlorophyll maximum to primary production, phytoplankton assemblages and diversity in a small stratified lake. J. Plankton Res. 2020, 42, 630–649. [Google Scholar] [CrossRef]
  4. Pannard, A.; Planas, D.; Beisner, B. Macrozooplankton and the persistence of the deep chlorophyll maximum in a stratified lake. Freshw. Biol. 2015, 60, 1717–1733. [Google Scholar] [CrossRef]
  5. Lofton, M.; Leach, T.; Beisner, B.; Carey, C. Relative importance of top-down vs. bottom-up control of lake phytoplankton vertical distributions varies among fluorescence-based spectral groups. Limnol. Oceanogr. 2020, 65, 2485–2501. [Google Scholar] [CrossRef]
  6. Pick, F.R.; Lean, D.R.S.; Nalewajko, C. Nutrient status of metalimnetic phytoplankton peaks. Limnol. Oceanogr. 1984, 29, 960–971. [Google Scholar] [CrossRef]
  7. Leach, T.H.; Beisner, B.E.; Carey, C.C.; Pernica, P.; Rose, K.C.; Huot, Y.; Brentrup, J.A.; Domaizon, I.; Grossart, H.-P.; Ibelings, B.W.; et al. Patterns and drivers of deep chlorophyll maxima structure in 100 lakes: The relative importance of light and thermal stratification. Limnol. Oceanogr. 2018, 63, 628–646. [Google Scholar] [CrossRef] [Green Version]
  8. Barbiero, R.P.; Tuchman, M.L. The Deep chlorophyll maximum in Lake Superior. J. Great Lakes Res. 2004, 30, 256–268. [Google Scholar] [CrossRef]
  9. Bramburger, A.J.; Reavie, E.D. A comparison of phytoplankton communities of the deep chlorophyll layers and epilimnia of the Laurentian Great Lakes. J. Great Lakes Res. 2016, 42, 1016–1025. [Google Scholar] [CrossRef]
  10. Twiss, M.; Ulrich, C.; Zastepa, A.; Pick, F. On phytoplankton growth and loss rates to microzooplankton in the epilimnion and metalimnion of Lake Ontario in mid-summer. J. Great Lakes Res. 2012, 38, 146–153. [Google Scholar] [CrossRef]
  11. Scofield, A.E.; Watkins, J.M.; Weidel, B.C.; Luckey, F.J.; Rudstam, L.G. The deep chlorophyll layer in Lake Ontario: Extent, mechanisms of formation, and abiotic predictors. J. Great Lakes Res. 2017, 43, 782–794. [Google Scholar] [CrossRef]
  12. Legrand, C.; Rengefors, K.; Fistarol, G.; Granéli, E. Allelopathy in phytoplankton—Biochemical, ecological and evolutionary aspects. Phycologia 2003, 42, 406–419. [Google Scholar] [CrossRef] [Green Version]
  13. Janssen, E.M.L. Cyanobacterial peptides beyond microcystins—A review on co-occurrence, toxicity, and challenges for risk assessment. Water Res. 2019, 151, 488–499. [Google Scholar] [CrossRef]
  14. Camacho, A. On the occurrence and ecological features of deep chlorophyll maxima (DCM) in Spanish stratified lakes. Limnetica 2006, 25, 453–478. [Google Scholar] [CrossRef]
  15. Silsbe, G.; Malkin, S. Where light and nutrients collide: The global distribution and activity of subsurface chlorophyll maximum layers. In Aquatic Microbial Ecology and Biogeochemistry: A Dual Perspective; Glibert, P.M., Kana, T.M., Eds.; Springer: Cham, Switzerland, 2016; pp. 141–152. [Google Scholar]
  16. Verschoor, M.; Powe, C.; McQuay, E.; Schiff, S.; Venkiteswaran, J.; Li, J.; Molot, L. Internal iron loading and warm temperatures are pre-conditions for cyanobacterial dominance in embayments along Georgian Bay, Great Lakes. Can. J. Fish. Aquat. Sci. 2017, 74, 1439–1453. [Google Scholar] [CrossRef]
  17. Bird, D.F.; Kalff, J. Phagotrophic sustenance of a metalimnetic phytoplankton peak. Limnol. Oceanogr. 1989, 34, 155–162. [Google Scholar] [CrossRef]
  18. Kurmayer, R.; Deng, L.; Entfellner, E. Role of toxic and bioactive secondary metabolites in colonization and bloom formation by filamentous cyanobacteria Planktothrix. Harmful Algae 2016, 54, 69–86. [Google Scholar] [CrossRef] [Green Version]
  19. Śliwińska-Wilczewska, S.; Wisniewska, K.; Konarzewska, Z.; Cieszyńska, A.; Felpeto, A.B.; Lewandowska, A.; Latała, A. The current state of knowledge on taxonomy, modulating factors, ecological roles, and mode of action of phytoplankton allelochemicals. Sci. Total. Environ. 2021, 773, 145681. [Google Scholar] [CrossRef]
  20. Reavie, E.D.; Barbiero, R.P.; Allinger, L.E.; Warren, G.J. Phytoplankton trends in the Great Lakes, 2001–2011. J. Great Lakes Res. 2014, 40, 618–639. [Google Scholar] [CrossRef]
  21. Findlay, D.L.; Paterson, J.J.; Hendzel, L.L.; Kling, H.J. Factors influencing Gonyostomum semen blooms in a small boreal reservoir lake. Hydrobiologia 2005, 533, 243–252. [Google Scholar] [CrossRef]
  22. Lebret, K.; Östman, Ö.; Langenheder, S.; Drakare, S.; Guillemette, F.; Lindström, E.S. High abundances of the nuisance raphidophyte Gonyostomum semen in brown water lakes are associated with high concentrations of iron. Sci. Rep. 2018, 8, 13463. [Google Scholar] [CrossRef]
  23. Head, R.M.; Jones, R.I.; Bailey-Watts, A.E. An assessment of the influence of recruitment from the sediment on the development of planktonic populations of cyanobacteria in a temperate mesotrophic lake. Freshw. Biol. 1999, 41, 759–769. [Google Scholar] [CrossRef]
  24. Brunberg, A.-K.; Blomqvist, P. Recruitment of Microcystis (Cyanophyceae) from lake sediments: The importance of littoral inocula. J. Phycol. 2003, 39, 58–63. [Google Scholar] [CrossRef]
  25. Holland, D.P.; Walsby, A.E. Viability of the cyanobacterium Planktothrix rubescens in the cold and dark, related to over-winter survival and summer recruitment in Lake Zürich. Eur. J. Phycol. 2008, 43, 179–184. [Google Scholar] [CrossRef]
  26. Kravchuk, E.; Ivanova, E.; Gladyshev, M. Spatial distribution of resting stages (akinetes) of the cyanobacteria Anabaena flos-aquae in sediments and its influence on pelagic populations. Mar. Freshw. Res. 2011, 62, 450–461. [Google Scholar] [CrossRef]
  27. Nürnberg, G.K.; LaZerte, B.D.; Olding, D.D. An Artificially induced Planktothrix rubescens surface bloom in a small kettle lake in Southern Ontario compared to blooms world-wide. Lake Reserv. Manag. 2003, 19, 307–322. [Google Scholar] [CrossRef] [Green Version]
  28. Grach-Pogrebinsky, O.; Sedmak, B.; Carmeli, S. Seco[D-Asp3]microcystin-RR and [D-Asp3,D-Glu(OMe)6]microcystin-RR, two new microcystins from a toxic water bloom of the cyanobacterium Planktothrix rubescens. J. Nat. Prod. 2004, 67, 337–342. [Google Scholar] [CrossRef]
  29. Cantin, A.; Beisner, B.E.; Gunn, J.M.; Prairie, Y.T.; Winter, J.G. Effects of thermocline deepening on lake plankton communities. Can. J. Fish. Aquat. Sci. 2011, 68, 260–276. [Google Scholar] [CrossRef]
  30. Salmaso, N.; Boscaini, A.; Shams, S.; Cerasino, L. Strict coupling between the development of Planktothrix rubescens and microcystin content in two nearby lakes south of the Alps (lakes Garda and Ledro). Ann. Limnol. Int. J. Limnol. 2013, 49, 309–318. [Google Scholar] [CrossRef] [Green Version]
  31. Kasprzak, P.; Shatwell, T.; Gessner, M.; Gonsiorczyk, T.; Kirillin, G.; Selmeczy, G.; Padisak, J.; Engelhardt, C. Extreme weatherevent triggers cascade towards extreme turbidity in a clear-water Lake. Ecosystems 2017, 20, 1407–1420. [Google Scholar] [CrossRef]
  32. Giling, D.P.; Nejstgaard, J.C.; Berger, S.A.; Grossart, H.-P.; Kirillin, G.; Penske, A.; Lentz, M.; Casper, P.; Sareyka, J.; Gessner, M.O. Thermocline deepening boosts ecosystem metabolism: Evidence from a large-scale lake enclosure experiment simulating a summer storm. Glob. Chang. Biol. 2017, 23, 1448–1462. [Google Scholar] [CrossRef]
  33. de Marsac, N.T. Occurrence and nature of chromatic adaptation in cyanobacteria. J. Bacteriol. 1977, 130, 82. [Google Scholar] [CrossRef] [Green Version]
  34. Micheletti, S.; Schanz, F.; Walsby, A.E. The daily integral of photosynthesis by Planktothrix rubescens during summer stratification and autumnal mixing in Lake Zürich. New Phytol. 1998, 139, 233–246. [Google Scholar] [CrossRef]
  35. Davis, P.A.; Beard, S.J.; Walsby, A.E. Variation in filament width and gas vesicles of red and green isolates of Planktothrix spp. Arch. Hydrobiol. Suppl. 2003, 146, 15–29. [Google Scholar]
  36. Legnani, E.; Copetti, D.; Oggioni, A.; Tartari, G.; Maria, T.; Palumbo, G.; Morabito, G. Planktothrix rubescens’ seasonal dynamics and vertical distribution in Lake Pusiano (North Italy). J. Limnol. 2005, 64, 61–73. [Google Scholar] [CrossRef]
  37. Copetti, D.; Tartari, G.; Morabito, G.; Oggioni, A.; Legnani, E.; Imberger, J. A biogeochemical model of Lake Pusiano (North Italy) and its use in the predictability of phytoplankton blooms: First preliminary results. J. Limnol. 2006, 65, 59–64. [Google Scholar] [CrossRef] [Green Version]
  38. Halstvedt, C.B.; Rohrlack, T.; Andersen, T.; Skulberg, O.; Edvardsen, B. Seasonal dynamics and depth distribution of Planktothrix spp. in Lake Steinsfjorden (Norway) related to environmental factors. J. Plankton Res. 2007, 29, 471–482. [Google Scholar] [CrossRef] [Green Version]
  39. Boström, B.; Andersen, J.M.; Fleischer, S.; Jansson, M. Exchange of phosphorus across the sediment-water interface. Hydrobiologia 1988, 170, 229–244. [Google Scholar] [CrossRef]
  40. Katsev, S.; Dittrich, M. Modeling of decadal scale phosphorus retention in lake sediment under varying redox conditions. Ecol. Model. 2013, 251, 246–259. [Google Scholar] [CrossRef]
  41. Loh, P.S.; Molot, L.A.; Nowak, E.; Nürnberg, G.K.; Watson, S.B.; Ginn, B. Evaluating relationships between sediment chemistry and anoxic phosphorus and iron release across three different water bodies. Inland Waters 2013, 3, 105–118. [Google Scholar] [CrossRef] [Green Version]
  42. Rentz, J.; Turner, I.; Ullman, J. Removal of phosphorus from solution using biogenic iron oxides. Water Res. 2009, 43, 2029–2035. [Google Scholar] [CrossRef] [PubMed]
  43. Buliauskaitė, R.; Wilfert, P.; Kumar, P.S.; de Vet, W.W.J.M.; Witkamp, G.-J.; Korving, L.; van Loosdrecht, M.C.M. Biogenic iron oxides for phosphate removal. Environ. Technol. 2020, 41, 260–266. [Google Scholar] [CrossRef]
  44. Hamre, K.D.; Lofton, M.E.; McClure, R.P.; Munger, Z.W.; Doubek, J.P.; Gerling, A.B.; Schreiber, M.E.; Carey, C.C. In situ fluorometry reveals a persistent, perennial hypolimnetic cyanobacterial bloom in a seasonally anoxic reservoir. Freshw. Sci. 2018, 37, 483–495. [Google Scholar] [CrossRef]
  45. Padisak, J.; Barbosa, F.; Koschel, R.; Krienitz, L. Deep layer cyanoprokaryota maxima in temperate and tropical lakes. Arch. Hydrobiol. 2003, 58, 175–199. [Google Scholar]
  46. Molot, L.A.; Watson, S.B.; Creed, I.F.; Trick, C.G.; McCabe, S.K.; Verschoor, M.J.; Sorichetti, R.J.; Powe, C.; Venkiteswaran, J.J.; Schiff, S.L. A novel model for cyanobacteria bloom formation: The critical role of anoxia and ferrous iron. Freshw. Biol. 2014, 59, 1323–1340. [Google Scholar] [CrossRef]
  47. Selmeczy, G.; Tapolczai, K.; Casper, P.; Krienitz, L.; Padisak, J. Spatial- and niche segregation of DCM-forming cyanobacteria in Lake Stechlin (Germany). Hydrobiologia 2015, 764, 229–240. [Google Scholar] [CrossRef]
  48. Padisak, J.; Krienitz, L.; Koschel, R.; Nedoma, J. Deep-layer autotrophic picoplankton maximum in the oligotrophic Lake Stechlin, Germany: Origin, activity, development and erosion. Eur. J. Phycol. 1997, 32, 403–416. [Google Scholar] [CrossRef]
  49. Taranu, Z.E.; Zastepa, A.; Creed, I.; Pick, F.R. Beyond total microcystins, the importance of examining different microcystin variants. LakeLine 2017, 37, 14–18. [Google Scholar]
  50. Rantala-Ylinen, A.; Känä, S.; Wang, H.; Rouhiainen, L.; Wahlsten, M.; Rizzi, E.; Berg, K.; Gugger, M.; Sivonen, K. Anatoxin-a synthetase gene cluster of the cyanobacterium Anabaena sp. Strain 37 and molecular methods to detect potential producers. Appl. Environ. Microbiol. 2011, 77, 7271–7278. [Google Scholar] [CrossRef] [Green Version]
  51. Beversdorf, L.J.; Weirich, C.A.; Bartlett, S.L.; Miller, T.R. Variable cyanobacterial toxin and metabolite profiles across six eutrophic lakes of differing physiochemical characteristics. Toxins 2017, 9, 62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Bartlett, S.; Brunner, S.; Klump, J.; Houghton, E.M.; Miller, T.R. Spatial analysis of toxic or otherwise bioactive cyanobacterial peptides in Green Bay, Lake Michigan. J. Great Lakes Res. 2018, 44, 924–933. [Google Scholar] [CrossRef] [PubMed]
  53. Chorus, I.; Sivonen, K.; Codd, G.A.; Börner, T.; Von Doehren, H.; Welker, M.; Dittmann, E.; Claussner, Y.; Christopffersen, K.; Schober, E.; et al. Toxic and Bioactive Peptides in Cyanobacteria. 2006. Available online: https://ci.nii.ac.jp/naid/10030356501/ (accessed on 26 June 2021).
  54. Rohrlack, T.; Christiansen, G.; Kurmayer, R. Putative antiparasite defensive system involving ribosomal and nonribosomal oligopeptides in cyanobacteria of the genus Planktothrix. Appl. Environ. Microbiol. 2013, 79, 2642–2647. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Charlton, M.; Mayne, G. Science and Monitoring Synthesis for South-Eastern Georgian Bay; Report prepared for Environment Canada. 2013. Available online: https://www.stateofthebay.ca/wp-content/uploads/2019/03/Science-and-Monitoring-Synthesis-for-South-Eastern-Georgian-Bay_Charlton-and-Mayne_2013.pdf (accessed on 26 June 2021).
  56. NLET (National Laboratory for Environmental Testing). Schedule of Services. Environment and Climate Change Canada, C.C.f.I.W.; NLET (National Laboratory for Environmental Testing): Burlington, ON, Canada, 1997. [Google Scholar]
  57. Findlay, D.; Kling, H. Protocols for Measuring Biodiversity: Phytoplankton in Freshwater. 2003. Available online: https://www.researchgate.net/publication/264881321_Protocols_for_measuring_biodiversity_Phytoplankton_in_freshwater (accessed on 26 June 2021).
Figure 1. Locations of study sites in Georgian Bay, Lake Huron (Ontario, Canada) denoted by red circles.
Figure 1. Locations of study sites in Georgian Bay, Lake Huron (Ontario, Canada) denoted by red circles.
Toxins 13 00445 g001
Figure 2. Seasonal (W, winter; Sp, spring; Su, summer; Au, autumn) changes in total phytoplankton biomass (line series) and extracted chlorophyll-a concentration (bars) in 2015. Near-surface (1 m) water column measurements are in panel (A) and at the depth of the recurring metalimnetic peaks, i.e., deep chlorophyll layers (DCLs) (~9 m in Twelve Mile Bay (TMB) and ~7 m in South Bay (SB), see Table 1 and Figure 3) are in panel (B). In the absence of DCLs during winter (W) and spring (Sp), i.e., during isothermal conditions, the sample was taken at the depth corresponding to summer DCL peak for comparison, i.e., ~9 m in TMB and ~7 m in SB. See Materials and Methods for sampling dates.
Figure 2. Seasonal (W, winter; Sp, spring; Su, summer; Au, autumn) changes in total phytoplankton biomass (line series) and extracted chlorophyll-a concentration (bars) in 2015. Near-surface (1 m) water column measurements are in panel (A) and at the depth of the recurring metalimnetic peaks, i.e., deep chlorophyll layers (DCLs) (~9 m in Twelve Mile Bay (TMB) and ~7 m in South Bay (SB), see Table 1 and Figure 3) are in panel (B). In the absence of DCLs during winter (W) and spring (Sp), i.e., during isothermal conditions, the sample was taken at the depth corresponding to summer DCL peak for comparison, i.e., ~9 m in TMB and ~7 m in SB. See Materials and Methods for sampling dates.
Toxins 13 00445 g002
Figure 3. Water column profiles in Twelve Mile Bay (panel A) and South Bay (panel B) for: phosphate (PO4), iron (Fe), manganese, ammonium (NH4) (●), nitrite/nitrate (NO2/3) (▲), temperature, chlorophyll-a, dissolved oxygen, pH, and redox in summer (August) of 2015. See Materials and Methods section for measurement and analytical protocols used.
Figure 3. Water column profiles in Twelve Mile Bay (panel A) and South Bay (panel B) for: phosphate (PO4), iron (Fe), manganese, ammonium (NH4) (●), nitrite/nitrate (NO2/3) (▲), temperature, chlorophyll-a, dissolved oxygen, pH, and redox in summer (August) of 2015. See Materials and Methods section for measurement and analytical protocols used.
Toxins 13 00445 g003
Figure 4. Vertical water column gradients of bioactive metabolites (microcystins and anabaenopeptins) produced by cyanobacteria and corresponding genes (mcy and apnDTe, respectively) in Twelve Mile Bay in summer (August) of 2018. Note that higher concentrations of these bioactive metabolites (and their genes) co-occur with biomass peaks, notably of Planktothrix isothrix, at depth (Table S1, 2018 data).
Figure 4. Vertical water column gradients of bioactive metabolites (microcystins and anabaenopeptins) produced by cyanobacteria and corresponding genes (mcy and apnDTe, respectively) in Twelve Mile Bay in summer (August) of 2018. Note that higher concentrations of these bioactive metabolites (and their genes) co-occur with biomass peaks, notably of Planktothrix isothrix, at depth (Table S1, 2018 data).
Toxins 13 00445 g004
Table 1. Depth of mixed (Zmix) and euphotic zones (Zeu) during mixed and stratified conditions and presence (Y: Yes/N: No) of deep-chlorophyll layer (DCL) in different embayments of Georgian Bay in the spring (June) and summer (August) of 2015. During mixed conditions in spring (June), Zmix was equivalent to the site depth (Zsite).
Table 1. Depth of mixed (Zmix) and euphotic zones (Zeu) during mixed and stratified conditions and presence (Y: Yes/N: No) of deep-chlorophyll layer (DCL) in different embayments of Georgian Bay in the spring (June) and summer (August) of 2015. During mixed conditions in spring (June), Zmix was equivalent to the site depth (Zsite).
MixedStratified
SiteZsite (m)Zeu (m)Zeu (m)Zmix (m)Zeu:ZmixDCL (m)
Byng Inlet53.0-5-N
Sturgeon Bay North142.44.580.6N
Sturgeon Bay South6----N
Deep Bay152.73.780.5N a
Woods Bay8----N
Twelve Mile Bay123.49.171.3Y (~9)
Tadenac Bay246.55.3110.5N
Go Home Bay103.6---N
Longuissa Bay93.84.780.6N
Cognashene Lake163.99.571.4N
Musqwash Bay West30----N
Musqwash Bay East20----N
North Bay183.26.480.8N
South Bay103.95.761.0Y (~7)
Honey Harbour74.63.670.5N
a Not detected during this study but reported in 2012 and 2014 by [16].
Table 2. Nutrient concentrations (µg/L) and depth ratios in the hypolimnion (Hypo) and epilimnion (Epi) in each embayment (DON = dissolved organic nitrogen). Data are based on the summer 2015 sampling expedition during stratification.
Table 2. Nutrient concentrations (µg/L) and depth ratios in the hypolimnion (Hypo) and epilimnion (Epi) in each embayment (DON = dissolved organic nitrogen). Data are based on the summer 2015 sampling expedition during stratification.
NutrientTwelve Mile BaySouth Bay
HypoEpiHypo: EpiHypoEpiHypo:Epi
P151412321318
Fe127816416010408
NO2/31191680.8
Si2.51.32842
Mn2100.366010200.52040
Co0.20.01200.30.0126
NH420545895.5107
SO4101112.37.70.3
DON21219613672991
PO40.40.80.51720.7246
Table 3. Major bioactive metabolites (right column) 1 and corresponding genes (left column) 2 (detect “+”; non-detect “-”; not measured “NM”; or maximum concentration (µg/L), any depth) in water samples from embayments along eastern Georgian Bay during 2014, 2015, and 2018. MC: Microcystins, NOD: Nodularin, CYN: Cylindrospermopsins, STX: Saxitoxins, ATX: Anatoxins, APT: Anabaenopeptins, CPT: Cyanopeptolins, MG: Microginins. NOD, CYN, and ATX were not detected in any of the samples.
Table 3. Major bioactive metabolites (right column) 1 and corresponding genes (left column) 2 (detect “+”; non-detect “-”; not measured “NM”; or maximum concentration (µg/L), any depth) in water samples from embayments along eastern Georgian Bay during 2014, 2015, and 2018. MC: Microcystins, NOD: Nodularin, CYN: Cylindrospermopsins, STX: Saxitoxins, ATX: Anatoxins, APT: Anabaenopeptins, CPT: Cyanopeptolins, MG: Microginins. NOD, CYN, and ATX were not detected in any of the samples.
SitesMCsSTXsAPTsCPTsMGs
Byng Inlet----------
Sturgeon Bay North+0.01--------
Sturgeon Bay SouthNMNMNMNMNMNMNMNMNMNM
Deep Bay----------
Woods BayNMNMNMNMNMNMNMNMNMNM
Twelve Mile Bay 3+0.4--+6.6+0.001--
Tadenac Bay+0.01--------
Go Home Bay----------
Longuissa Bay----------
Cognashene Lake----------
Musqwash Bay WestNMNMNMNMNMNMNMNMNMNM
Musqwash Bay EastNMNMNMNMNMNMNMNMNMNM
North Bay--+-------
South Bay+0.03--+0.01+0.001+-
Honey Harbour+0.04--------
1 ELISA and PPIA results verified by LC-MS/MS. 2 PCR primers listed in Table A2 (adapted from [18,50]. 3 Peaks in Twelve Mile Bay were measured in 2018 and the highest concentrations of these bioactive metabolites (and their genes) co-occurred with increasing biomass, notably of Planktothrix cf. isothrix, at the depth of the deep chlorophyll layer (Figure 4, Table S1, 2018 data).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zastepa, A.; Miller, T.R.; Watson, L.C.; Kling, H.; Watson, S.B. Toxins and Other Bioactive Metabolites in Deep Chlorophyll Layers Containing the Cyanobacteria Planktothrix cf. isothrix in Two Georgian Bay Embayments, Lake Huron. Toxins 2021, 13, 445. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13070445

AMA Style

Zastepa A, Miller TR, Watson LC, Kling H, Watson SB. Toxins and Other Bioactive Metabolites in Deep Chlorophyll Layers Containing the Cyanobacteria Planktothrix cf. isothrix in Two Georgian Bay Embayments, Lake Huron. Toxins. 2021; 13(7):445. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13070445

Chicago/Turabian Style

Zastepa, Arthur, Todd R. Miller, L. Cynthia Watson, Hedy Kling, and Susan B. Watson. 2021. "Toxins and Other Bioactive Metabolites in Deep Chlorophyll Layers Containing the Cyanobacteria Planktothrix cf. isothrix in Two Georgian Bay Embayments, Lake Huron" Toxins 13, no. 7: 445. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13070445

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop