Next Article in Journal
Comparative Genomics and Identification of an Enterotoxin-Bearing Pathogenicity Island, SEPI-1/SECI-1, in Staphylococcus epidermidis Pathogenic Strains
Next Article in Special Issue
Determination of the Role of Microcystis aeruginosa in Toxin Generation Based on Phosphoproteomic Profiles
Previous Article in Journal
Tolerance and Excretion of the Mycotoxins Aflatoxin B1, Zearalenone, Deoxynivalenol, and Ochratoxin A by Alphitobius diaperinus and Hermetia illucens from Contaminated Substrates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Algal Blooms and Cyanotoxins in Jordan Lake, North Carolina

1
Department of Marine, Earth and Atmospheric Sciences, North Carolina State University, Raleigh, NC 27695, USA
2
North Carolina Department of Environmental Quality, Division of Water Resources, Raleigh, NC 27699, USA
*
Author to whom correspondence should be addressed.
Toxins 2018, 10(2), 92; https://doi.org/10.3390/toxins10020092
Submission received: 31 December 2017 / Revised: 17 February 2018 / Accepted: 19 February 2018 / Published: 24 February 2018
(This article belongs to the Special Issue Harmful Algal Bloom Dynamics)

Abstract

:
The eutrophication of waterways has led to a rise in cyanobacterial, harmful algal blooms (CyanoHABs) worldwide. The deterioration of water quality due to excess algal biomass in lakes has been well documented (e.g., water clarity, hypoxic conditions), but health risks associated with cyanotoxins remain largely unexplored in the absence of toxin information. This study is the first to document the presence of dissolved microcystin, anatoxin-a, cylindrospermopsin, and β-N-methylamino-l-alanine in Jordan Lake, a major drinking water reservoir in North Carolina. Saxitoxin presence was not confirmed. Multiple toxins were detected at 86% of the tested sites and during 44% of the sampling events between 2014 and 2016. Although concentrations were low, continued exposure of organisms to multiple toxins raises some concerns. A combination of discrete sampling and in-situ tracking (Solid Phase Adsorption Toxin Tracking [SPATT]) revealed that microcystin and anatoxin were the most pervasive year-round. Between 2011 and 2016, summer and fall blooms were dominated by the same cyanobacterial genera, all of which are suggested producers of single or multiple cyanotoxins. The study’s findings provide further evidence of the ubiquitous nature of cyanotoxins, and the challenges involved in linking CyanoHAB dynamics to specific environmental forcing factors are discussed.
Key Contribution: Four cyanotoxins could be detected simultaneously in Jordan Lake. Concentrations of dissolved toxins were generally low (well below EPA thresholds) but found year-round. Passive in-situ tracking approaches proved powerful in characterizing toxin dynamics.

1. Introduction

The eutrophication of waterways causes water quality issues worldwide and these may intensify with climate change [1,2,3,4]. One issue, linked to excess nutrient input from agricultural land and urbanized areas, is Harmful Algal Blooms (HABs) [5,6,7]. In freshwater systems and estuaries, HABs are typically dominated by cyanobacteria (or blue-green algae, CyanoHABs) [6,8,9] that can have multiple adverse effects on aquatic ecosystems, from the blocking of sunlight to benthic vegetation, to oxygen depletion that may kill fish [1,10,11,12]. Global annual estimates of the socioeconomic costs of CyanoHABs are significant and range from millions to billions of dollars (e.g., water monitoring and testing, drinking water treatment, adverse impacts on recreational use and fisheries) [13,14,15].
Various environmental factors impact the initiation, peak and demise of a CyanoHAB. Increased nutrients, mainly nitrogen (N) and phosphorus (P), have long been associated with bloom development [16,17], while other triggers, such as shifts in nutrient ratios throughout a bloom’s lifecycle, may also play a role in cyanobacterial composition shifts. Decreasing N:P ratios can promote bloom-forming cyanobacterial genera, capable of N-fixation [17,18,19,20,21], and low-flow conditions within lakes or estuaries reportedly favor the growth of cyanobacteria over other algal taxa [22,23,24]. Rising water temperatures have been linked to increased bloom activity [25,26,27] and to potential shifts from non-toxic to toxic strains [28]. At least a dozen cyanobacterial genera have been implicated with toxin production, and at least eight toxin groups have been characterized, of which microcystin (MCY) has been studied most extensively [13,29]. However, not all species within a genus can produce toxins and those that can, do not do so continuously. Field and laboratory studies show that MCY concentrations tend to be positively correlated with dissolved inorganic nutrients (mainly N and P), temperature and light levels [30,31,32]. For instance, temperature optima for MCY production by Microcystis and Dolichospermum strains were reported between 18 to 25 °C, and for Dolichospermum spp., temperature seems to influence which MCY congener is produced [32,33]. While many studies report that absolute nutrient concentrations are linked to cyanotoxin presence, shortages of certain nutrients that lead to shifts in nutrient ratios may be a factor. N-limitation was linked to increased MCY and anatoxin-a (ANA) concentrations by N-fixing members of the genera Aphanizomenon, Cylindrospermopsis and Dolichospermum [30,32], while P-limitation was associated with low MCY production by Dolichospermum spp. and Microcystis spp., and with low ANA levels in Aphanizomenon spp. [30,32]. In laboratory experiments, light levels of <20 µmol photons m−2 s−1 seem to be conducive to MCY production [30,32,34]. Overall, a better understanding of the complex interplay between environmental factors, cyanobacterial growth and/or the onset of toxin production is needed in order to mitigate, and ultimately, prevent, CyanoHAB-related issues within given environments.
Cyanotoxin consumption can harm fish, livestock, pets and humans in varying ways [35,36]. Exposure to MCY and cylindrospermopsin (CYN) can impair liver function and at high doses can be lethal [35,37,38,39]. ANA and saxitoxin (STX) are both neurotoxins [29,39]. ANA causes an overstimulation in neuromuscular junctions, leading to respiratory failure [39]. STX is responsible for paralytic shellfish poisoning (PSP), a condition that can cause paralysis and death in humans [40,41,42,43]. More recently, β-N-methylamino-l-alanine (BMAA) has been investigated for its connection to neurological diseases, including amyotrophic lateral sclerosis (ALS), Alzheimer’s disease and Parkinson’s disease [44,45,46]. While an increasing number of studies are addressing potential health risks due to these substances, major knowledge gaps remain in regard to exposure pathways, concentrations in the field and environmental triggers for toxin production.
A recent US-wide survey of over 1100 lakes showed that MCY was present in 32% of the tested lakes (range = below detection (BD) to 230 µg L−1; average = 3.0 µg L−1), and at least one of the cyanotoxins (MCY, CYN, STX or ANA) could be detected in 92% of the States [47]. CYN was reported in 4.0% (range = BD to 4.4 µg L−1; average = 0.6 µg L−1), STX in 7.7% (range = BD to 0.38 µg L−1; average = 0.06 µg L−1) and ANA in approximately 0.3% of samples (range and average not given) [47]. Comprehensive toxin surveys, especially for multiple toxins, are still rare, and risk assessment by the World Health Organization (WHO) is mainly based on chlorophyll a (Chl a) levels and cyanobacterial abundance ranges. As Loftin et al. [47] demonstrate, these metrics, in contrast to actual toxin information, can lead to an overestimation of MCY risk. This overestimation is partly due to the fact that, as aforementioned, not all cyanobacteria are toxin producers, and known toxic species do not produce toxins continuously. Under the Drinking Water Protection Act, the US Environmental Protection Agency (EPA) released national 10-day health advisories for MCY of 0.3 µg L−1 for infants and 1.6 µg L−1 for adults, and for CYN, 0.7 µg L−1 for infants and 3 µg L−1 for adults, based on body weight and water intake [48,49]. In December 2016, the EPA also suggested recreational guidelines of 4 µg L−1 for MCY and 8 µg L−1 for CYN [50]. State-specific advisories do exist for ANA and STX in several states, but little information is available on BMAA, and therefore, there are no guidelines at this time [51]. In order to protect human and ecosystem health, more monitoring is essential for evaluating exposure risks, especially for these emerging substances.
For North Carolina (NC), only limited information is available on the presence of cyanotoxins. For instance, MCY was detected at low levels in 11 reservoirs across the Piedmont during summer 2002 [52], and in four lakes during 2011 and 2012 [53] (<0.8 µg L−1 MCY for all three studies). The EPA National Lakes Assessment 2007 also reported MCY and STX in NC waters, but CYN or ANA presence could not be confirmed [47]. While most previous records indicate a low MCY exposure risk based on WHO guidelines, concentrations of over 800 µg L−1, measured in Waterville Reservoir in October of 2007, are a reminder of how little is known about natural toxin ranges and their spatiotemporal dynamics in NC [54]. Most water treatment plants have procedures in place to eliminate toxins from drinking water [55,56], but the pervasiveness of cyanotoxins raises questions on chronic recreational exposure (e.g., swimming, boating, wading) [57,58] or the potential for food web poisoning via fish or shellfish consumption [10,57,59,60].
For this study, CyanoHAB dynamics were characterized in NC’s B. Everett Jordan Reservoir (henceforth, “Jordan Lake”), based on a six-year data set compiled through the North Carolina Division of Water Resources (NCDWR) Ambient Lakes Monitoring program. In addition to the continued collection of community structure data, the specific goal of this study was to test for the year-round presence of multiple cyanotoxins in an artificial reservoir that provides drinking water for nearly 300,000 people in locations such as Morrisville, Cary, and Apex. As the number of lakes and reservoirs that experience severe CyanoHAB blooms increases in the US and worldwide, newly developed approaches for measuring varying toxins are slow to be implemented in routine surveys. Although the presence of MCY was confirmed at multiple sites during summer 2002 (average 0.2 µg L−1) [47] within Jordan Lake, this is, to the best of our knowledge, the first comprehensive study to investigate five common cyanotoxins (2014–2016) using a combination of traditional and recently-developed tracking approaches. Cyanotoxin data were interpreted in relation to phytoplankton dynamics over multiple years (2011–2016) and in relation to pertinent environmental parameters throughout Jordan Lake.

2. Results

2.1. Phytoplankton Dynamics

Overall, twelve cyanobacterial genera were identified with Anabaenopsis, Aphanizomenon, Aphanocapsa, Aphanothece, Chroococcus, Cylindrospermopsis, Dolichospermum, Microcystis, Merismopedia, Planktolyngbya, Pseudanabaena, and Raphidiopsis across the nine sampling sites (Figure 1, Table 1 and Table S1). In addition, 48 microeukaryote phytoplankton genera/species could be distinguished belonging to the diatoms, chlorophytes, chrysophytes, cryptophytes, euglenophytes, prymnesiophytes, and dinoflagellates (Table S1). Cyanobacteria dominated the phytoplankton assemblages based on cell counts (94% of total phytoplankton; range = 5% to 100%), while they constituted, on average, 39% (range ≤ 1% to 98%) to total phytoplankton biovolume. Cyanobacterial and microeukaryote phytoplankton abundance varied slightly with season, year and across the lake (r = 0.077 to 0.360 at p < 0.0003; three-way ANOSIM; Table S2). Peak densities were reached during summer (range = 4.3 × 103 to 5.0 × 106 and 180 to 3.2 × 105, respectively) and fall months (250 to 4.3 × 106 cells mL−1 and 180 to 1.1 × 105, respectively; Figure 2 and Figure 3) and, in agreement with cell abundance, Chl a values reached their maximum during late summer/early fall, after a first initial increase typically during spring (overall range = 1 to 128 µg L−1; mean = 41 µg L−1; Figure 2).
Community structure based on Bray–Curtis similarities for both cyanobacteria and microeukaryote phytoplankton varied with season, site and year (r = 0.070 to 0.434; at p < 0.003; three-way ANOSIM; Table S3). For cyanobacteria, these similarities were highest during summer and fall (50% and 45% of the community were shared, respectively) compared to winter and spring (28% and 32%, respectively; Figure 4A). The less abundant microeukaryote assemblages consistently shared between 40% and 46% of their makeup within each of the seasons and throughout the year. Analyses of intra-annual community structure changes further revealed a recurrent pattern where the species composition followed a cyclic year-round pattern resulting in a “reset” of the assemblage by the onset of the following year (RELATE test; ρ = 0.022 to 0.639, p < 0.05). In total, 82% of the yearly datasets tested positive for cyclicity (shown for site A in 2015 in Figure 4B).

2.2. Physicochemical, Meteorological and Hydrological Parameters

Temperature, NOx (nitrate plus nitrite) and dissolved oxygen (DO) concentrations in surface waters (surface to twice Secchi depth) showed the strongest seasonal changes (r = 0.146 to 0.614, p = 0.0001; one-way ANOSIM; Table S4), with temperature maxima in summer and NOx, DO and ammonia (NH3) levels higher during colder months (Figure 5). Total Kjeldahl nitrogen (TKN = particulate and dissolved organic N plus ammonia) and TKN:TP ratios slightly increased during summer months, while no consistent seasonal trends were observed over time for total phosphorous (TP = particulate and dissolved fractions; Table S4). NH3, DO and turbidity also varied somewhat between years (Table S4). All tested environmental parameters showed some spatial variability, except for temperature, NH3 and DO (r = 0.07 to 0.385, p = 0.0001; one-way ANOSIM; Table S4). Throughout the sampling period, surface waters at the sample location with the deepest water column depth of 12.2 m (site B) had some of the lowest NOx and TP concentrations and the highest TKN:TP ratios (Table S5). NOx, TP and turbidity tended to be higher at some of the shallower sites (e.g., C, D, F and H in Figure 1, Table S5). Average TKN and TP concentrations ranged from 0.78 to 1.16 mg L−1 and 0.04 to 0.13 mg L−1, respectively (Table S5). Overall, mean surface DO levels (~0 to 3 m) ranged from 8.0 to 9.7 mg L−1 across the stations with individual measurements rarely falling beyond 4 mg L−1 (in <1% of measurements; Table S5).
Combining all of the physicochemical data for each sampling event resulted in “environmental fingerprints” which varied with month, season, and site (average r = 0.16 to 0.348 at p = 0.0001; one-way ANOSIM; Figure 6). The most notable difference was seen when comparing physicochemical settings among seasons—conditions varied little throughout the summer (49% similar) compared to fall, winter and spring (<5% similar; Figure 6). Also available were single meteorological and hydrological values to characterize conditions across the entire lake during each sampling date (Figure 7, Table S6). Weekly precipitation averaged 0.02 cm h−1 with yearly maxima occurring throughout late spring to early fall (Figure 7A, Table S6). Wind speeds showed their maxima during the winter, while Photosynthetically Active Radiation (PAR) peaked in the summer (Figure 7B,C). Finally, overall river flow (Haw River, Morgan Creek and New Hope Creek combined) varied considerably from year to year with maxima typically occurring during spring or fall (Figure 7D). Haw River flow (mean = 27.79 m3 s−1) exceeded flows for Morgan Creek and New Hope Creek (0.12 and 0.14 m3 s−1, respectively).

2.3. Cyanotoxins

Four out of the five tested toxins were detected in Jordan Lake based on discrete (grab) samples collected between August 2015 and December 2016 at stations A through G (Figure 1). Dissolved MCY was confirmed in 10 out of 65 samples, ANA in 39 out of 69, CYN in six out of 63 and BMAA in nine out of 64 samples (Figure 8 and Figure 9, Table 2). STX presence could not be confirmed (n = 40; LDL (low detection limit) = 0.015).
In addition to the grab samples, in-situ toxin tracking was employed (Solid Phase Adsorption Toxin Tracking or SPATT) to confirm the presence of dissolved MCY, ANA and CYN at stations E and G (Table 2, Figure 8A). The combined sampling approaches revealed the occurrence of MCY at multiple sites (A, B and D through G) and throughout all seasons (Figure 9). SPATTs allowed for the confirmation of MCY in 92% of samples, but the toxin was only detected in 15% of discrete samples (Table 2). Dissolved CYN was measured at sites E and G during spring, summer and fall but not during winter (Figure 8 and Figure 9). CYN presence was indicated using both SPATT (13%) and grab (10%) sampling (Figure 8, Table 2). Similar to MCY, dissolved ANA was also found at multiple sites (A through G) and during all seasons (Figure 9). SPATTs confirmed ANA in 100% and grab sampling in 57% of tests across all sites (Table 2). Finally, dissolved BMAA was found at four sites (D through G) during fall, winter and spring (14% of samples; Figure 9, Table 3). All in all, multiple toxins were present at six out of seven sites and during 30 out of 69 sampling events. As stated earlier, to minimize loss of ANA and STX in lake water with a pH outside the range of 5 to 7, the addition of a diluent is recommended (Abraxis manual). Since samples prior to October 2016 were not treated with diluent and 96% of the lake water samples during that survey period measured above a pH of 7 (range = 5.3 to 9.5; mean = 7.8), both ANA and STX presence may have been underestimated in this study.
Due to a limited number of positives for CYN, STX and BMAA, statistical analyses to examine spatiotemporal trends were limited to dissolved MCY and ANA. MCY showed higher concentrations during summer and fall at sites E and G based on in-situ tracking (r = 0.194, p = 0.024; n = 24; one-way ANOSIM Table S7), but no trend was indicated based on grab sampling (n = 65). ANA concentrations did not vary significantly over time based on SPATTS data (n = 23), and while concentrations based on discrete sampling indicated some spatial variability (r = 0.094, p = 0.018; n = 69; one-way ANOSIM; Table S7), no consistent trend was apparent across the lake. For the most commonly detected toxins, MCY and ANA, concentrations were also examined for possible linkages to cyanobacterial composition shifts. For MCY, changes in concentration could be linked to shifts in genera—Pseudanabaena, Merismopedia and Aphanothece—while changes in Raphidiopsis spp. abundance linked to variance in ANA concentrations (Bio-Env [BEST] routine; ρ = 0.283 and 0.183, respectively, at p = 0.0001).

2.4. Linkages between Environmental Factors and Phytoplankton Dynamics

Correlation analyses indicated positive relationships for Chl a, cyanobacterial and microeukaryote density and biovolume (r = 0.314 to 0.851, p < 0.05; Table 3). Increases in dissolved MCY were correlated with increases in cyanobacterial density and biovolume as well as microeukaryote phytoplankton density and dissolved ANA (r = 0.270 to 0.408, p < 0.05). ANA showed a positive relationship with dissolved MCY, microphytoplankton biovolume and TKN:TP ratios (Table 3). Chl a as well as cyanobacterial and microeukaryote abundances correlated negatively with NH3 and NOx but increased with TKN, temperature and pH (Table 3). Only the Chl a concentration was associated with increases in TP and turbidity (Table 3). There was no statistical significance when these analyses were conducted using average values for Chl a, cell densities and biovolumes across the lake (average for all stations) in regard to their relationships with meteorological and hydrological parameters (PAR, river flow, wind speed or precipitation).
Multiple regression analyses were performed and indicated that NH3, NOx, TKN and DO explained up to 68% of variance in Chl a (Table 4). NOx and TKN, combined with pH and DO, were linked to 52% of variance in cyanobacterial densities and, without DO, explained up to 41% of variance in cyanobacterial biovolumes (Table 4). Only 13% and 19% of variance in microeukaryote phytoplankton densities and biovolume could be linked to a similar suite of physicochemical parameters (Table 4). Finally, community structure patterns (Bray–Curtis similarity matrices) were matched to varying combinations of the physicochemical, meteorological and hydrological variables (environmental fingerprints) using a BEST routine [61], and these trend analyses indicated that NOx and temperature correlated most strongly with community structure patterns for cyanobacteria (ρ = 0.4 at p = 0.0001), and NOx, TKN and Morgan Creek flow data correlated with changed phytoplankton community structure for the microeukaryotes (ρ = 0.22 at p = 0.0001). No significant correlations were found when BEST routines were performed to match lake-wide community structure patterns to meteorological and hydrological variables.

3. Discussion

3.1. Cyanotoxins and Phytoplankton Dynamics in Jordan Lake

CyanoHABs are a worldwide problem that has resulted in the development of WHO guidelines to assess risks (low, moderate and high) from MCY exposure based on toxin concentration, Chl a and cyanobacterial density [1,2,3,12]. However, applying these three metrics, a water body can be at risk based on one, but not all, of these criteria. For instance, for over 1100 lakes in the US, agreement for risk assessment based on all three parameters was only observed for 27% of the systems [47]. Given this discrepancy and the fact that most monitoring programs routinely measure Chl a and cyanobacterial density, but rarely employ approaches to measure toxins, complicates the tasks of water resource managers to protect designated lake uses and human health. Jordan Lake has been known for water quality issues due to eutrophication and recurrent CyanoHABs since its impoundment in the early 1980s. Despite its importance as a drinking water source for nearly 300,000 people and its recreational use by over a million visitors annually (NC Department of Natural and Cultural Resources [62], potential health risks from cyanotoxin presence had remained largely unexplored. Over a 2-year study period, from 2014 to 2016, a total of 36% of the examined samples from Jordan Lake tested positive for MCY but only one discrete sample (1.98 µg MCY L−1, site G on 21 June 2016) exceeded WHO guidelines for drinking water, with 1 µg MCY L−1, while values never reached those for EPA recreational guidelines of 4 µg L−1 [50,63]. The average concentration of dissolved MCY (0.06 µg L−1) was within the range of values observed in other NC freshwater systems (0.05 to 0.54 µg L−1) and across the US (BDL to 230 µg L−1) [47]. Notably, levels remained well below concentrations in CyanoHAB-prone systems, such as Lake Erie, where typical annual maxima peak at ~200 µg MCY L−1 and, in one instance, exceeded 1200 µg L−1 [64]. Following the aforementioned official WHO guidelines [47], Jordan Lake would be generally categorized as high risk based on its Chl a and cyanobacterial density; however, based on this study, only a low risk was observed for both MCY and CYN from 2014 to 2016. These findings further corroborate how critical toxin information is for the refinement of health risk metrics that directly inform lake-specific management decisions but also help shape national and international guidelines.
CyanoHABs may consist of multiple forms of toxins, but limited data is currently available on where and when toxins co-occur and under what environmental conditions [47,65,66,67]. To our knowledge, this study is unique in providing a year-round and multi-year record that allows the confirmation of the presence and co-occurrence of four cyanotoxins (dissolved MCY, CYN, ANA and BMAA) in a US freshwater body (Table 2 and Table 5). In contrast to this study, previous state-wide surveys tested a small number of samples (~seven or less) and these were typically collected during one season. Such limited temporal coverage is common and increases the probability of missing toxic events in any water body, due to the ephemeral nature of CyanoHABs. In NC, for instance, the detection of an unprecedented high MCY level of over 800 µg L−1 in Waterville Reservoir in 2007 raises questions on how well natural toxin ranges and spatiotemporal dynamics have been assessed based on traditional grab sampling and existing monitoring frequencies [54]. CYN had been tested for, but was not detected in, any major NC freshwater body [47,53]. ANA genes were found in six lakes (City Lake, Oak Hollow Lake, Randleman Reservoir, Lake Brandt, Lake Mackintosh, and Belews Lake) but the toxin itself was not confirmed [47,53]. Little to no information is currently available on BMAA presence or STX in NC lakes and reservoirs. US-wide, STX was confirmed in 7.7% of lakes during the National Lakes Assessment 2007, and this survey included one NC lake: Lake Rhodhiss (Table 5) [47]. Universally, more comprehensive datasets are needed to aid the development of risk thresholds for newly emerging cyanotoxins (e.g., no national guidelines currently exist for ANA or BMAA [51,68], to allow comparisons across freshwater systems and regions and to begin to inform epidemiological studies on the possible synergistic effects of multiple toxins.
A lack of information on co-occurring toxins typically goes hand-in-hand with limited data on year-round toxin dynamics. In this study, discrete toxin sampling that provided momentary snapshots of conditions was complemented by year-round in-situ tracking (SPATTs approach) at two sites in Jordan Lake, from 2014 to 2016. The advantages of employing SPATTs come from their higher sensitivity in detecting low toxin levels via a time-integrative signal, their use in freshwater to marine environments, the facilitation for testing multiple toxins, and their easy deployment and recovery [69]. A major limitation of using passive samplers, especially as a sole toxin detection approach, comes from the semi-quantitative nature of the data that can currently not be linked to regulatory limits and, hence, makes risk determination in systems difficult [69]. SPATT-based average MCY concentrations in Jordan Lake (36.2 ng (g resin)−1 d−1) fell within the ranges reported in several California studies (19.6–137.7 ng (g resin)−1 d−1) [70,71], and, similar to those reports, the SPATT method proved more sensitive for MCY detection compared to grab sampling. In Jordan Lake, MCY was present year-round, with 92% of the tested samples based on in-situ tracking compared to only 15% based on grab sampling (Table 2). This supports the effectiveness of in-situ tracking approaches in addressing emerging concerns in regard to the potential impacts of chronic or subacute exposure for wildlife and humans [72,73]. Using SPATTs for the detection of cyanotoxins other than MCY requires careful consideration of resin type [74]. As such, the hydrophobic HP-20 resin used for this study was thoroughly tested for its efficiency in detecting MCY but has not been fully evaluated for its efficiency in adsorbing other toxin types (e.g., CYN or ANA) [74]. For instance, only a total of three samples tested positive for CYN in this study (Figure 9B, Table 2). While both SPATTs and grab samples allowed for consistent detection of ANA, relative concentrations based on SPATTs were relatively low compared to MCY levels based on in-situ accumulation. This difference could have been a direct consequence of ANA being less prevalent throughout deployment periods, which would lower accumulation potential, but was likely also an artifact of the toxin not being efficiently adsorbed and/or retained during prolonged deployment [74,75]. An increasing number of studies have been conducted to test resins for the detection of algal toxins, to better evaluate the potential of in-situ passive samplers, to inform future health risk assessments and management decisions (review in [74]).
Phytoplankton assemblages in Jordan Lake were dominated by cyanobacteria (~94% based on cell density) with microeukaryote phytoplankton only rarely outnumbering the prokaryotes during non-bloom months. Changes in overall community structure followed consistent intra-annual patterns for both cyanobacteria and less abundant microphytoplankton. All of the six most abundant cyanobacterial taxa, identified via microscopy, were potential producers of single or multiple toxins, which included Pseudanabaena spp. (MCY), Cylindrospermopsis raciborskii (CYN, ANA, STX, BMAA), Aphanocapsa delicatissima and A. pulchra (MCY), Chroococcus spp. (MCY), several species of Dolichospermum (MCY, CYN, ANA, STX, BMAA) and Microcystis aeruginosa and M. firma (MCY, BMAA) [76,77]. Of these main genera, Aphanocapsa, Cylindrospermopsis and Pseudanabaena occurred in 2-year dominance shifts, a pattern that could not be linked to any of the physicochemical or hydrological factors tested in this study. Exploring the relationships between toxin presence and cyanobacterial community data indicated that relative changes in the abundance of Pseudanabaena spp., Merismopedia punctata and Aphanothece saxicola were linked to shifts in MCY, and Raphidiopsis spp. abundances were associated with changes in dissolved ANA [78,79,80,81]. Whether these taxa were truly responsible for toxin production remains unconfirmed and would have required further taxonomic resolution on the species and strain levels, since toxicity is not a genus-specific trait, nor is toxin production continuous. Combining field studies, such as this, with culture-based trials using isolates will allow us to verify taxonomic affiliations based on genomics, explore gene expression and tie findings to meta-omics profiles for natural cyanobacterial communities [28,53].

3.2. Environmental Factors in Relation to Phytoplankton and Toxin Dynamics

Jordan Lake has been consistently rated as eutrophic or hyper-eutrophic, and nutrient input from urban (26%) and agricultural (16%) land uses upstream (the remaining 58% are forested) serves as important stimulant for phytoplankton growth [82]. NOx and NH3 concentrations were positively correlated with overall river flow in this study, and the availability of both nitrogen sources subsequently declined with increasing algal biomass, cell densities and total TKN. This overall shift in TKN was most likely attributed to the incorporation of N into algal biomass. A suite of environmental factors, including NOx, TKN, NH3, DO and pH, was linked to 68% of the variance in Chl a and 52% in cyanobacterial density (Table 4). In contrast, TKN and NOx were associated with only 19% of the variance observed for microeukaryote phytoplankton, indicating that these main algal groups flourish under different environmental conditions. In agreement with studies elsewhere, changes in temperature, together with nutrient availability (i.e., NOx), were linked to shifts in overall cyanobacterial community structure [9,25,28,83]. Only a weak correlative relationship was observed between microphytoplankton composition and a combination of NOx, TKN, and river flow (Morgan Creek). Information on additional key environmental factors is needed to further characterize the significance of these potential forcing factors for microphytoplankton but also for cyanobacterial and toxin dynamics in Jordan Lake. For instance, urea has been reported to specifically stimulate cyanobacteria [16,84], and changes in the availability of both urea and inorganic P have been linked to increased abundances of toxic species [85,86]. Additionally, shifts in N:P ratios have been suggested to promote N-fixing cyanobacteria [17,87], a group also represented in Jordan Lake (i.e., genera Cylindrospermopsis, Dolichospermum and Pseudanabaena) [88,89,90]. Examining potential linkages for the two most commonly detected toxins, ANA and MCY, only revealed a positive correlation between TKN:TP ratios and dissolved ANA, giving some indication that P might have been less readily available relative to N. However, this remains speculative since TKN and TP estimates included varying dissolved and particulate fractions, and no separate information was obtained on the availability of dissolved P to further explore relationships between toxin and dissolved versus cell-bound nutrients. As cyanotoxin production may be tied to a complex array of environmental conditions, the collection of larger toxin datasets seems imperative for deciphering these linkages.
Finally, the role that river flow and water retention time play for CyanoHAB dynamics in Jordan Lake warrants closer examination in future studies. The lower arm of the lake, where river input is higher due to the Haw River (90% of water inflow into the lake), has estimated retention times of about 5 days, while they reportedly exceed 400 days in the upper arms (New Hope and Morgan Creeks) of Jordan Lake (NCDWR). Previous studies have shown that low flow conditions tend to result in larger blooms and shifts in species composition [22,24,91,92]. However, examining correlations between algal abundances and/or community structure among individual stations or regions (grouped stations adjacent to river outflows) did not yield significant differences. Moreover, and in contrast to previous studies that reported drought conditions favorable to bloom activity [23,52,93], no correlative relationships were seen between precipitation, algal abundances and community structure over the investigated study period.

3.3. Conclusions and Recommendations

This study is the first to show that cyanobacterial communities in Jordan Lake are linked with the recurrence of multiple cyanotoxins throughout the year. These findings fall in line with an increasing number of studies that have confirmed the ubiquitous nature of cyanotoxins, their simultaneous presence in varying environments and the need for further research to characterize the conditions that favor toxin production. The continued development and employment of highly sensitive toxin-tracking approaches (e.g., SPATTs), together with an expanding tool-kit for genomic and transcriptomic testing, will be essential for examining cause–effect relationships and providing the knowledge needed to predict the likelihood for toxin exposure via varying exposure pathways, be it to single or multiple agents. The presented study approach can inform efforts in similar water bodies where continued issues with eutrophication due to increasing population growth threaten water quality. The study’s findings serve as a baseline to better characterize CyanoHAB events in Jordan Lake and guide continued testing for selected toxins as part of routine water quality monitoring to protect the lake’s dedicated uses (i.e., drinking water and recreation).

4. Materials and Methods

4.1. Study Area and Data Collection

Jordan Lake is a 56 km2 artificial reservoir in central NC in Chatham County, constructed between 1967 and 1983 by the U.S. Army Corps of Engineers. The lake is filled by the Haw River in the south and Morgan Creek and New Hope Creek in the north (Figure 1). The Haw River is the largest of the three inflows and accounts for 70–90% of the total water volume entering the lake [94]. The reservoir has an average depth of 4.9 m and volume of 265 × 106 m3 [52,95]. Sampling by NCDWR was conducted via small boats, as outlined in the DWR manual for standard operating procedures [96,97]. Briefly, Secchi depth was measured, and physical data (temperature, DO, and pH) were collected from the surface to depth at approximately 1 m intervals throughout the water column, using either Hydrolab (Hach Environmental, Loveland, CO, USA) or YSI (Yellow Springs Instrument Co., Yellow Springs, OH, USA) sondes. Grab samples for community structure and chemical analyses (NH3, NOx, TKN, TP and turbidity) were collected via a depth-integrated sampler from the surface to twice Secchi depth. Transport and processing for phytoplankton, nutrients, Chl a and turbidity followed standard protocols detailed in DWR’s standard operating manual [97]. Drought index measurements were obtained from DWR’s online Drought Monitor History database based on weekly drought averages for drought conditions by percent area for Chatham County (representative of the upper Cape Fear River watershed). Hourly meteorological parameters (wind speed, precipitation, and PAR) were obtained from the Reedy Creek Field Laboratory (State Climate Office of North Carolina), located approximately 27 km from Jordan Lake (35.807° N, 78.744° W). Daily hydrological data, obtained from the US Geological Survey (USGS), included stream discharge (flow) for Morgan Creek (near Chapel Hill, NC at 35.89333° N, 79.01972° W, site ID 02097517), New Hope Creek (near Blands, NC at 35.885° N, 78.96528° W, site ID 02097314) and Haw River (near Bynum, NC at 35.76528° N, 79.13583° W, site ID 02096960).

4.2. Phytoplankton Data

Phytoplankton community analyses were conducted microscopically by NCDWR using Leitz inverted microscopes and Utermöhl counting chambers [98,99]. Briefly, samples were preserved with Lugol’s solution (0.4% final concentration) upon collection, and a 5 mL subsample was settled for 24 h. Samples were analyzed within 14 days of collection. Samples were counted until 100 units (single cells, colonies or filaments, depending on the specific taxon) of the most dominant taxa were recorded. Taxonomic identification [100] was established to at least genus level. Biovolumes were calculated using cell densities (cells mL−1) multiplied by reference values [99].

4.3. Toxin Analyses

Discrete (grab) samples were collected at approximately 0.5 m depth using pre-cleaned (acid washed followed by three Milli-Q [MQ] water rinses) polyethylene terephthalate glycol (PETG) bottles and after pre-rinsing the bottles with lake water. The bottles were chilled on ice in a cooler for transport to the lab, where 50 mL aliquots were filtered through GF/F filters with a nominal 0.7 µm pore size (Whatman grade, GE Healthcare Life Sciences, Chicago, IL, USA). The filtrate was collected in glass scintillation vials for analysis of dissolved (extracellular) toxins and stored frozen at −20 °C until analysis using commercially available ELISAs (Abraxis Inc., Warminster, PA, USA; see details below). Samples were analyzed using a BioTek ELx800 Absorbance Microplate Reader (BioTek, Winooski, VT, USA). Dissolved samples for ANA and STX were pretreated with a diluent to prevent toxin loss, following the manufacturer’s guidelines (Abraxis), except for samples collected prior to October 2016. Toxin analyses were conducted for sites A through G from August 2015 through December 2016.
In addition to grab samples, Solid Phase Adsorption Toxin Tracking (SPATT) [71,101] units were used to determine in-situ toxin accumulation over approximately monthly intervals (average deployment time was 28 days). SPATTs were deployed at 2 sites at 0.5 m depth (Figure 1, sites E and G) from August 2015 to December 2016. Construction, deployment and extraction procedures for SPATTs followed previously published guidelines [101]. Briefly, 3 g of HP-20 resin (Sigma–Aldrich, St. Louis, MO, USA) was activated in 100% methanol for 30 min, then rinsed with three equivalent volumes of MQ-water and sonicated for 45 s at 50% amplitude with a sonic dismembrator (Fisher Scientific, Hampton, NH, USA, Model FB120). After sonication, activated bags were stored in chilled MQ water in the refrigerator until deployment [101]. Buoys with a weighted rope were deployed with a mesh bag containing two SPATTs attached at sites E and G. Retrieved units were kept out of direct sunlight, put on ice for transport to the lab and transferred into a −80 °C freezer within ~2 h of collection. The resin from the SPATT bags was extracted according to previously published protocols [70] with the following modifications: samples were vortexed before each of the three extractions and extracts two and three were combined for analysis, while extract one was run separately. All cyanotoxin analyses for SPATT extracts and dissolved samples were conducted using ELISAs (Abraxis Inc., Warminster, PA, USA): MCY-ADDA (Product #520011; sensitive to MCY-LR, -YR, -LF, -RR, LW, and nodularin; LDL = 0.10 µg L−1), CYN (Product #522011; sensitive to CYN and deoxy-CYN; LDL = 0.04 µg L−1), ANA (Product #520060; sensitive to anatoxin-a and homoanatoxin-a; LDL = 0.1 µg L−1), STX (Product #52255B; sensitive to STX and other paralytic shellfish poison [PSP] toxins; LDL = 0.015 µg L−1), and BMAA (Product #520040; sensitive to BMAA and other amino acids; limit of quantitation = 4 ng mL−1). Immediately prior to analysis, SPATT samples were diluted with the sample diluent provided with each ELISA kit to avoid methanol interference during assays. Final methanol concentrations were <5% for MCY, <20% for CYN and <2.5% for ANA (ELISA manuals and Abraxis recommendations). In addition, diluted SPATT extracts were centrifuged for 2 min at 13,000 rpm at room temperature (Eppendorf 5424 R Microcentrifuge) to remove any particulate matter. SPATT results were normalized as nanograms toxin per gram resin per day (ng toxin (g resin)−1 d−1). STX and BMAA were not analyzed using the SPATT approach due to adsorption bias when using HP-20 resin (see further details in the discussion) [74].

4.4. Chlorophyll Analyses

Chl a concentration (µg L−1) was determined using the EPA method 445.0 via fluorescence [102]. Briefly, 50–100 mL aliquots of lake water were concentrated onto 0.7 µm GF/F filters, extracted using acetone and measured fluorometrically (Turner Designs Model 10 Series fluorometer).

4.5. Statistical Analyses

Statistical analyses were performed using the PRIMER v7 [61] and STATISTICA 13 (TIBCO Software) statistical software packages. Community data (cell densities and biovolumes) were square-root transformed and compared based on Bray–Curtis similarity values, while physical environmental parameters (averaged over the upper water column from the surface to twice Secchi depth) were log-transformed, normalized (mean subtracted from each value and divided by the standard deviation) and compared after the computation of Euclidean distance resemblance matrices [61]. Three-way ANOSIM (analysis of similarity) tests were computed to examine temporal (seasonal [spring: March–May; summer: June–August; fall: September–November; winter: December–February] and yearly) as well as spatial trends (crossed design; 9999 permutations). This resulted in r values which represent a measure of distinction between groups. For instance, r values of 0 indicated that groups were similar, while an r value of 1, or close to 1, implied that groups were dissimilar. Similarity patterns over temporal or spatial scales were further illustrated using non-metric multidimensional scaling (nMDS) plots where more closely clustered data points represented higher similarity. Stress values were calculated for MDS plots to reflect the level of distortion that results from representing similarity rankings between multiple samples in a two-dimensional space. Generally, a stress value of <0.2 indicates an accurate representation of similarity rankings [61]. Additionally, one-way ANOSIM tests were conducted to examine whether toxin concentrations varied over time and location. Cyanobacteria and phototrophic microeukaryotes were analyzed separately and combined (total phytoplankton), and temporal or spatial differences were determined based on abundance data (cells mL−1) as well as biovolume (mm3 m−3). In contrast to the environmental data, phytoplankton data was square-root transformed and compared using Bray–Curtis similarity indices [61].
Chl a, cell densities, biovolumes and toxin concentrations were examined for their relationships with physical (temperature, DO and pH) and chemical (NH3, NOx, TKN, TP and turbidity) parameters using correlation and regression analyses (Pearson’s product-moment correlations, r; adjusted coefficient of determination, R2; STATISTICA 13, TIBCO Software). The same analyses were conducted using average values for Chl a, cell densities and biovolumes across the lake (all stations combined per sampling date) and related to meteorological and hydrological parameters (PAR, river flow, wind speed or precipitation). A BEST routine (PRIMER v7) was used to establish matches between similarities in site-specific physicochemical data (Euclidean distance-based matrix) and community structure information (Bray–Curtis similarity matrix) using Spearman’s rank correlations (rho, ρ) [61]. This routine was also repeated to link community data across the lake (all stations combined per sampling date) to meteorological and hydrological parameters. The BEST routine, unlike multiple regression analyses, cannot differentiate among positive or negative relationships, but identifies similarities between the two matrices. The RELATE test was used to test a cyclical model for an annual resetting of algal and microeukaryote assemblages at differing stations and in each of the 6 years [61].

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6651/10/2/92/s1, Table S1: Most abundant cyanobacterial and microeukaryote phytoplankton groups, identified to genus or, where possible, species level using microscopy, Table S2: Results of three-way ANOSIM tests comparing cyanobacteria, microeukaryote phytoplankton and total phytoplankton abundances across season, year and site (crossed design), Table S3: Results of three-way ANOSIM tests comparing cyanobacteria, microeukaryote phytoplankton and total phytoplankton community structure across season, year and site (crossed design), Table S4: Results of one-way ANOSIM comparing environmental, meteorological and hydrological parameters across month, season, year and site, Table S5: Minimum (Min), maximum (Max) and average (Ave) values for parameters at each site, Table S6: Minimum (Min), maximum (Max) and average (Ave) values for meteorological and hydrological parameters, Table S7: Results of one-way ANOSIM comparing MCY and ANA concentrations based on SPATTs data and ANA levels based on grab sampling across month, season, year and site.

Acknowledgments

Thank you to the DWR Ambient Monitoring Team for support during field surveys. Support for in-situ tracking deployments was provided by Tarek Aziz, Dan Obenour and Jeremy Smithheart (NCSU). Assistance with toxin analyses was provided by Dave Deardorff and Justin Carpino (Abraxis) and for PRIMER statistical analyses by Bob Clarke. This study was supported through funds from the Urban Water Consortium (UWC)/Water Resources Research Institute (WRRI) (#16-11-U) and North Carolina Sea Grant (#5104348).

Author Contributions

D.W. and A.S. conceived the study and wrote the paper; J.G. oversaw field sampling, D.W., J.G., M.V.B. and E.F. participated in sample collection and analyses.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Selman, M.; Greenhalgh, S.; Diaz, R.; Sugg, Z. Eutrophication and hypoxia in coastal areas: A global assessment of the state of knowledge. World Resour. Inst. 2008, 284, 1–6. [Google Scholar]
  2. Kundzewicz, Z.W.; Mata, L.J.; Arnell, N.W.; Döll, P.; Jimenez, B.; Miller, K.; Oki, T.; Sen, Z.; Shiklomanov, I. The implications of projected climate change for freshwater resources and their management. Hydrol. Sci. J. 2008, 53, 3–10. [Google Scholar] [CrossRef]
  3. Rabalais, N.N.; Turner, R.E.; Díaz, R.J.; Justic, D. Global change and eutrophication of coastal waters. ICES J. Mar. Sci. 2009, 66, 1528–1537. [Google Scholar] [CrossRef]
  4. Woodward, G.; Perkins, D.M.; Brown, L.E. Climate change and freshwater ecosystems: Impacts across multiple levels of organization. Philos. Trans. R. Soc. B 2010, 365, 2093–2106. [Google Scholar] [CrossRef] [PubMed]
  5. Paerl, H.W.; Huisman, J. Climate change: A catalyst for global expansion of harmful cyanobacterial blooms. Environ. Microbiol. Rep. 2009, 1, 27–37. [Google Scholar] [CrossRef] [PubMed]
  6. O’Neil, J.M.; Davis, T.W.; Burford, M.A.; Gobler, C.J. The rise of harmful cyanobacteria blooms: The potential roles of eutrophication and climate change. Harmful Algae 2012, 14, 313–334. [Google Scholar] [CrossRef]
  7. Paerl, H.W.; Gardner, W.S.; Havens, K.E.; Joyner, A.R.; McCarthy, M.J.; Newell, S.E.; Qin, B.; Scott, J.T. Mitigating cyanobacterial harmful algal blooms in aquatic ecosystems impacted by climate change and anthropogenic nutrients. Harmful Algae 2016, 54, 213–222. [Google Scholar] [CrossRef] [PubMed]
  8. Dokulil, M.T.; Teubner, K. Cyanobacterial dominance in lakes. Hydrobiologia 2000, 438, 1–12. [Google Scholar] [CrossRef]
  9. Paerl, H.W.; Fulton, R.S., 3rd; Moisander, P.H.; Dyble, J. Harmful freshwater algal blooms, with an emphasis on cyanobacteria. Sci. World J. 2001, 1, 76–113. [Google Scholar] [CrossRef] [PubMed]
  10. Ibelings, B.W.; Chorus, I. Accumulation of cyanobacterial toxins in freshwater “seafood” and its consequences for public health: A review. Environ. Pollut. 2007, 150, 177–192. [Google Scholar] [CrossRef] [PubMed]
  11. Paerl, H.W.; Pinckney, J.L.; Fear, J.M.; Peierls, B.L. Ecosystem responses to internal and watershed organic matter loading: Consequences for hypoxia in the eutrophying Neuse River Estuary, North Carolina, USA. Mar. Ecol. Prog. Ser. 1998, 166, 17–25. [Google Scholar] [CrossRef]
  12. Watson, S.B.; Miller, C.; Arhonditsis, G.; Boyer, G.L.; Carmichael, W.; Charlton, M.N.; Confesor, R.; Depew, D.C.; Höök, T.O.; Ludsin, S.A. The re-eutrophication of Lake Erie: Harmful algal blooms and hypoxia. Harmful Algae 2016, 56, 44–66. [Google Scholar] [CrossRef] [PubMed]
  13. Steffensen, D.A. Economic cost of cyanobacterial blooms. In Cyanobacterial Harmful Algal Blooms: State of the Science and Research Needs; Hudnell, H.K., Ed.; Springer: New York, NY, USA, 2008; pp. 855–865. [Google Scholar]
  14. Le, C.; Zha, Y.; Li, Y.; Sun, D.; Lu, H.; Yin, B. Eutrophication of lake waters in China: Cost, causes, and control. Environ. Manag. 2010, 45, 662–668. [Google Scholar] [CrossRef] [PubMed]
  15. Hamilton, D.P.; Wood, S.A.; Dietrich, D.R.; Puddick, J. Costs of harmful blooms of freshwater cyanobacteria. In Cyanobacteria: An Economic Perspective; Sharma, N.K., Rai, A.K., Stal, L.J., Eds.; Wiley: Chichester, UK, 2014; pp. 245–256. [Google Scholar]
  16. Kratz, W.A.; Myers, J. Nutrition and growth of several blue-green algae. Am. J. Bot. 1955, 42, 282–287. [Google Scholar] [CrossRef]
  17. Schindler, D.W. Evolution of phosphorus limitation in lakes. Science 1977, 195, 260–262. [Google Scholar] [CrossRef] [PubMed]
  18. Havens, K.E.; James, R.T.; East, T.L.; Smith, V.H. N:P ratios, light limitation, and cyanobacterial dominance in a subtropical lake impacted by non-point source nutrient pollution. Environ. Pollut. 2003, 122, 379–390. [Google Scholar] [CrossRef]
  19. Paerl, H. Nutrient and other environmental controls of harmful cyanobacterial blooms along the freshwater–marine continuum. In Cyanobacterial Harmful Algal Blooms: State of the Science and Research Needs; Hudnell, H.K., Ed.; Springer: New York, NY, USA, 2008; pp. 217–237. [Google Scholar]
  20. Liu, X.; Lu, X.; Chen, Y. The effects of temperature and nutrient ratios on Microcystis blooms in Lake Taihu, China: An 11-year investigation. Harmful Algae 2011, 10, 337–343. [Google Scholar] [CrossRef]
  21. Lee, T.A.; Rollwagen-Bollens, G.; Bollens, S.M. The influence of water quality variables on cyanobacterial blooms and phytoplankton community composition in a shallow temperate lake. Environ. Monit. Assess. 2015, 187, 315. [Google Scholar] [CrossRef] [PubMed]
  22. Cloern, J.E. Annual variations in river flow and primary production in the South San Francisco Bay Estuary (USA). Estuar. Coast. 1991, 19, 91–96. [Google Scholar]
  23. Van Vliet, M.T.H.; Zwolsman, J.J.G. Impact of summer droughts on the water quality of the Meuse river. J. Hydrol. 2008, 353, 1–17. [Google Scholar] [CrossRef]
  24. Michalak, A.M.; Anderson, E.J.; Beletsky, D.; Boland, S.; Bosch, N.S.; Bridgeman, T.B.; Chaffin, J.D.; Cho, K.; Confesor, R.; Daloglu, I.; et al. Record-setting algal bloom in Lake Erie caused by agricultural and meteorological trends consistent with expected future conditions. Proc. Natl. Acad. Sci. USA 2013, 110, 6448–6452. [Google Scholar] [CrossRef] [PubMed]
  25. Konopka, A.; Brock, T.D. Effect of temperature on blue-green algae (cyanobacteria) in Lake Mendota. Appl. Environ. Microbiol. 1978, 36, 572–576. [Google Scholar] [PubMed]
  26. Paerl, H.W.; Huisman, J. Blooms like it hot. Science 2008, 320, 57–58. [Google Scholar] [CrossRef] [PubMed]
  27. Wu, Y.; Li, L.; Zheng, L.; Dai, G.; Ma, H.; Shan, K.; Wu, H.; Zhou, Q.; Song, L. Patterns of succession between bloom-forming cyanobacteria Aphanizomenon flos-aquae and Microcystis and related environmental factors in large, shallow Dianchi Lake, China. Hydrobiologia 2016, 765, 1–13. [Google Scholar] [CrossRef]
  28. Davis, T.W.; Berry, D.L.; Boyer, G.L.; Gobler, C.J. The effects of temperature and nutrients on the growth and dynamics of toxic and non-toxic strains of Microcystis during cyanobacteria blooms. Harmful Algae 2009, 8, 715–725. [Google Scholar] [CrossRef]
  29. Cheung, M.Y.; Liang, S.; Lee, J. Toxin-producing cyanobacteria in freshwater: A review of the problems, impact on drinking water safety, and efforts for protecting public health. J. Microbiol. 2013, 51, 1–10. [Google Scholar] [CrossRef] [PubMed]
  30. Sivonen, K. Cyanobacterial Toxins. In Encyclopedia of Microbiology, 3rd ed.; Schaechter, M., Ed.; Elsevier: Oxford, UK, 2009. [Google Scholar]
  31. Orr, P.T.; Jones, G.J. Relationship between microcystin production and cell division rates in nitrogen-limited Microcystis aeruginosa cultures. Limnol. Oceanogr. 1998, 43, 1604–1614. [Google Scholar] [CrossRef]
  32. Kaebernick, M.; Neilan, B.A. Ecological and molecular investigations of cyanotoxin production. FEMS Microbiol. Ecol. 2001, 35, 1–9. [Google Scholar] [CrossRef] [PubMed]
  33. Rapala, J.; Sivonen, K.; Lyra, C.; Niemela, S.I. Variation of microcystins, cyanobacterial hepatotoxins, in Anabaena spp. as a function of growth stimuli. Appl. Environ. Microbiol. 1997, 63, 2206–2212. [Google Scholar] [PubMed]
  34. Gao, K.; Yu, H.; Brown, M.T. Solar PAR and UV radiation affects the physiology and morphology of the cyanobacterium Anabaena sp. PCC 7120. J. Photochem. Photobiol. B 2007, 89, 117–124. [Google Scholar] [CrossRef] [PubMed]
  35. Råbergh, C.M.I.; Bylund, G.; Eriksson, J.E. Histopathological effects of microcystin-LR, a cyclic peptide toxin from the cyanobacterium (blue-green alga) Microcystis aeruginosa on common carp (Cyprinus carpio L.). Aquat. Toxicol. 1991, 20, 131–145. [Google Scholar] [CrossRef]
  36. Carmichael, W.W.; Boyer, G.L. Health impacts from cyanobacteria harmful algae blooms: Implications for the North American Great Lakes. Harmful Algae 2016, 54, 194–212. [Google Scholar] [CrossRef] [PubMed]
  37. Hawkins, P.R.; Runnegar, M.T.; Jackson, A.R.; Falconer, I.R. Severe hepatotoxicity caused by the tropical cyanobacterium (blue-green alga) Cylindrospermopsis raciborskii (Woloszynska) Seenaya and Subba Raju isolated from a domestic water supply reservoir. Appl. Environ. Microbiol. 1985, 50, 1292–1295. [Google Scholar] [PubMed]
  38. Falconer, I.R.; Humpage, A.R. Cyanobacterial (blue-green algal) toxins in water supplies: Cylindrospermopsins. Environ. Toxicol. 2006, 21, 299–304. [Google Scholar] [CrossRef] [PubMed]
  39. Falconer, I.R. Health effects associated with controlled exposures to cyanobacterial toxins. In Cyanobacterial Harmful Algal Blooms: State of the Science and Research Needs; Hudnell, H.K., Ed.; Springer: New York, NY, USA, 2008; pp. 607–612. [Google Scholar]
  40. Acres, J.; Gray, J. Paralytic shellfish poisoning. Can. Med. Assoc. J. 1978, 119, 1195–1197. [Google Scholar] [PubMed]
  41. Halstead, B.W.; Schantz, E.J. Paralytic Shellfish Poisoning; Report prepared for World Health Organization; World Health Organization: Geneva, Switzerland, 1984. [Google Scholar]
  42. Kaas, H.; Henriksen, P. Saxitoxins (PSP toxins) in Danish lakes. Water Res. 2000, 34, 2089–2097. [Google Scholar] [CrossRef]
  43. Wiese, M.; D’Agostino, P.M.; Mihali, T.K.; Moffitt, M.C.; Neilan, B.A. Neurotoxic alkaloids: Saxitoxin and its analogs. Mar. Drugs 2010, 8, 2185–2211. [Google Scholar] [CrossRef] [PubMed]
  44. Hirano, A.; Kurland, L.T.; Krooth, R.S.; Lessell, S. Parkinsonism-dementia complex, an endemic disease on the island of Guam: I. Clinical features. Brain 1961, 84, 642–661. [Google Scholar] [CrossRef] [PubMed]
  45. Murch, S.J.; Cox, P.A.; Banack, S.A.; Steele, J.C.; Sacks, O.W. Occurrence of β-methylamino-l-alanine (BMAA) in ALS/PDC patients from Guam. Acta Neurol. Scand. 2004, 110, 267–269. [Google Scholar] [CrossRef] [PubMed]
  46. Banack, S.A.; Caller, T.A.; Stommel, E.W. The cyanobacteria derived toxin beta-N-methylamino-l-alanine and amyotrophic lateral sclerosis. Toxins 2010, 2, 2837–2850. [Google Scholar] [CrossRef] [PubMed]
  47. Loftin, K.A.; Graham, J.L.; Hilborn, E.D.; Lehmann, S.C.; Meyer, M.T.; Dietze, J.E.; Griffith, C.B. Cyanotoxins in inland lakes of the United States: Occurrence and potential recreational health risks in the EPA National Lakes Assessment 2007. Harmful Algae 2016, 56, 77–90. [Google Scholar] [CrossRef] [PubMed]
  48. Drinking Water Health Advisory for the Cyanobacterial Microcystin Toxins; Report Prepared for Environmental Protection Agency; Environmental Protection Agency: Washington, DC, USA, 2015.
  49. Drinking Water Health Advisory for the Cyanobacterial Toxin Cylindrospermopsin; Report Prepared for Environmental Protection Agency; Environmental Protection Agency: Washington, DC, USA, 2015.
  50. Human Health Recreational Ambient Water Quality Criteria or Swimming Advisories for Microcystins and Cylindrospermopsins—Draft; Report Prepared for Environmental Protection Agency; Environmental Protection Agency: Washington, DC, USA, 2016.
  51. Farrer, D.; Counter, M.; Hillwig, R.; Cude, C. Health-based cyanotoxin guideline values allow for cyanotoxin-based monitoring and efficient public health response to cyanobacterial blooms. Toxins 2015, 7, 457–477. [Google Scholar] [CrossRef] [PubMed]
  52. Touchette, B.W.; Burkholder, J.M.; Allen, E.H.; Alexander, J.L.; Kinder, C.A.; Brownie, C.; James, J.; Britton, C.H. Eutrophication and cyanobacteria blooms in run-of-river impoundments in North Carolina, USA. Lake Reserv. Manag. 2007, 23, 179–192. [Google Scholar] [CrossRef]
  53. Grubbs, L.F. Quantification of Select Cyanobacteria and Cyanotoxins in Piedmont North Carolina Lakes Using Real-Time PCR. Master’s Thesis, The University of North Carolina at Greensboro, Greensboro, NC, USA, 2014. [Google Scholar]
  54. Allender, C.J.; LeCleir, G.R.; Rinta-Kanto, J.M.; Small, R.L.; Satchwell, M.F.; Boyer, G.L.; Wilhelm, S.W. Identifying the source of unknown microcystin genes and predicting microcystin variants by comparing genes within uncultured cyanobacterial cells. Appl. Environ. Microbiol. 2009, 75, 3598–3604. [Google Scholar] [CrossRef] [PubMed]
  55. Himberg, K.; Keijola, A.M.; Hiisvirta, L.; Pyysalo, H.; Sivonen, K. The effect of water treatment processes on the removal of hepatotoxins from Microcystis and Oscillatoria cyanobacteria: A laboratory study. Water Res. 1989, 23, 979–984. [Google Scholar] [CrossRef]
  56. Hitzfeld, B.C.; Hoger, S.J.; Dietrich, D.R. Cyanobacterial toxins: Removal during drinking water treatment, and human risk assessment. Environ. Health Perspect. 2000, 108, 113–122. [Google Scholar] [CrossRef] [PubMed]
  57. Stewart, I.; Webb, P.M.; Schluter, P.J.; Fleming, L.E.; Burns, J.W.; Gantar, M.; Backer, L.C.; Shaw, G.R. Epidemiology of recreational exposure to freshwater cyanobacteria—An international prospective cohort study. BMC Public Health 2006, 6, 93. [Google Scholar] [CrossRef] [PubMed]
  58. Backer, L.C.; McNeel, S.V.; Barber, T.; Kirkpatrick, B.; Williams, C.; Irvin, M.; Zhou, Y.; Johnson, T.B.; Nierenberg, K.; Aubel, M. Recreational exposure to microcystins during algal blooms in two California lakes. Toxicon 2010, 55, 909–921. [Google Scholar] [CrossRef] [PubMed]
  59. Garcia, A.C.; Bargu, S.; Dash, P.; Rabalais, N.N.; Sutor, M.; Morrison, W.; Walker, N.D. Evaluating the potential risk of microcystins to blue crab (Callinectes sapidus) fisheries and human health in a eutrophic estuary. Harmful Algae 2010, 9, 134–143. [Google Scholar] [CrossRef]
  60. Lehman, P.W.; Teh, S.J.; Boyer, G.L.; Nobriga, M.L.; Bass, E.; Hogle, C. Initial impacts of Microcystis aeruginosa blooms on the aquatic food web in the San Francisco Estuary. Hydrobiologia 2010, 637, 229–248. [Google Scholar] [CrossRef]
  61. Clarke, K.R.; Warwick, R.N. Change in Marine Communities: An Approach to Statistical Analysis and Interpretation, 3rd ed.; PRIMER-E: Plymouth, UK, 2001. [Google Scholar]
  62. NC Department of Natural and Cultural Resources. Available online: https://www.ncdcr.gov/ (accessed on 22 February 2018).
  63. Cyanobacterial Toxins: Microcystin-LR in Drinking Water; Background Document for Preparation of WHO Guidelines for Drinking-Water Quality; Report Prepared for World Health Organizaition; World Health Organization: Geneva, Switzerland, 2003.
  64. Steffen, M.M.; Belisle, B.S.; Watson, S.B.; Boyer, G.L.; Wilhelm, S.W. Status, causes and controls of cyanobacterial blooms in Lake Erie. J. Gt. Lakes Res. 2014, 40, 215–225. [Google Scholar] [CrossRef]
  65. Howard, M.D.A.; Nagoda, C.; Kudela, R.M.; Hayashi, K.; Tatters, A.; Caron, D.A.; Busse, L.; Brown, J.; Sutula, M.; Stein, E.D. Microcystin prevalence throughout lentic waterbodies in coastal Southern California. Toxins 2017, 9, 231. [Google Scholar] [CrossRef] [PubMed]
  66. Graham, J.L.; Loftin, K.A.; Meyer, M.T.; Ziegler, A.C. Cyanotoxin mixtures and taste-and-odor compounds in cyanobacterial blooms from the Midwestern United States. Environ. Sci. Technol. 2010, 44, 7361–7368. [Google Scholar] [CrossRef] [PubMed]
  67. Al-Sammak, M.A.; Hoagland, K.D.; Snow, D.D.; Cassada, D. Methods for simultaneous detection of the cyanotoxins BMAA, DABA, and anatoxin-a in environmental samples. Toxicon 2013, 76, 316–325. [Google Scholar] [CrossRef] [PubMed]
  68. Health Effects Support Document for the Cyanobacterial Toxin Anatoxin-A; Report Prepared for Environmental Protection Agency; Environmental Protection Agency: Washington, DC, USA, 2015.
  69. MacKenzie, L.; Beuzenberg, V.; Holland, P.; McNabb, P.; Selwood, A. Solid Phase Adsorption Toxin Tracking (SPATT): A new monitoring tool that simulates the biotoxin contamination of filter feeding bivalves. Toxicon 2004, 44, 901–918. [Google Scholar] [CrossRef] [PubMed]
  70. Kudela, R.M. Characterization and deployment of Solid Phase Adsorption Toxin Tracking (SPATT) resin for monitoring of microcystins in fresh and saltwater. Harmful Algae 2011, 11, 117–125. [Google Scholar] [CrossRef]
  71. Gibble, C.M.; Kudela, R.M. Detection of persistent microcystin toxins at the land–sea interface in Monterey Bay, California. Harmful Algae 2014, 39, 146–153. [Google Scholar] [CrossRef]
  72. White, S.H.; Duivenvoorden, L.; Fabbro, L.D. A decision-making framework for ecological impacts associated with the accumulation of cyanotoxins (cylindrospermopsin and microcystin). Lakes Reserv. Res. Manag. 2005, 10, 25–37. [Google Scholar] [CrossRef]
  73. Zanchett, G.; Oliveira-Filho, E.C. Cyanobacteria and cyanotoxins: From impacts on aquatic ecosystems and human health to anticarcinogenic effects. Toxins 2013, 5, 1896–1917. [Google Scholar] [CrossRef] [PubMed]
  74. Kudela, R.M. Passive sampling for freshwater and marine algal toxins. Compr. Anal. Chem. 2017, 78, 379–409. [Google Scholar]
  75. Kudela, R.M.; University of California Santa Cruz, Santa Cruz, CA, USA. Personal communication, 2017.
  76. Pearson, L.; Mihali, T.; Moffitt, M.; Kellmann, R.; Neilan, B. On the chemistry, toxicology and genetics of the cyanobacterial toxins, microcystin, nodularin, saxitoxin and cylindrospermopsin. Mar. Drugs 2010, 8, 1650–1680. [Google Scholar] [CrossRef] [PubMed]
  77. Otten, T.G.; Paerl, H.W. Health effects of toxic cyanobacteria in US drinking and recreational waters: Our current understanding and proposed direction. Curr. Environ. Health Rep. 2015, 2, 75–84. [Google Scholar] [CrossRef] [PubMed]
  78. Namikoshi, M.; Murakami, T.; Watanabe, M.F.; Oda, T.; Yamada, J.; Tsujimura, S.; Nagai, H.; Oishi, S. Simultaneous production of homoanatoxin-a, anatoxin-a, and a new non-toxic 4-hydroxyhomoanatoxin-a by the cyanobacterium Raphidiopsis mediterranea Skuja. Toxicon 2003, 42, 533–538. [Google Scholar] [CrossRef]
  79. Carmichael, W.W.; Li, R. Cyanobacteria toxins in the Salton Sea. Saline Syst. 2006, 2, 5. [Google Scholar] [CrossRef] [PubMed]
  80. Furtado, A.L.F.F.; do Carmo Calijuri, M.; Lorenzi, A.S.; Honda, R.Y.; Genuário, D.B.; Fiore, M.F. Morphological and molecular characterization of cyanobacteria from a Brazilian facultative wastewater stabilization pond and evaluation of microcystin production. Hydrobiologia 2009, 627, 195–209. [Google Scholar] [CrossRef]
  81. Paerl, H.W.; Otten, T.G. Harmful cyanobacterial blooms: Causes, consequences, and controls. Microb. Ecol. 2013, 65, 995–1010. [Google Scholar] [CrossRef] [PubMed]
  82. Lake, B. Everett Jordan Watershed Model Report. Available online: ftp://ftp.tjcog.org/pub/planning/water/JordanAllocationModel/TTDataFiles/Documents/Model_Development_Reports/Jordan_Watershed_Model_Report_July2014final.pdf (accessed on 10 December 2017).
  83. Bullerjahn, G.S.; McKay, R.M.; Davis, T.W.; Baker, D.B.; Boyer, G.L.; D’Anglada, L.V.; Doucette, G.J.; Ho, J.C.; Irwin, E.G.; Kling, C.L.; et al. Global solutions to regional problems: Collecting global expertise to address the problem of harmful cyanobacterial blooms. A Lake Erie case study. Harmful Algae 2016, 54, 223–238. [Google Scholar] [CrossRef] [PubMed]
  84. Valladares, A.; Montesinos, M.L.; Herrero, A.; Flores, E. An ABC-type, high-affinity urea permease identified in cyanobacteria. Mol. Microbiol. 2002, 43, 703–715. [Google Scholar] [CrossRef] [PubMed]
  85. Davis, T.W.; Harke, M.J.; Marcoval, M.A.; Goleski, J.; Orano-Dawson, C.; Berry, D.L.; Gobler, C.J. Effects of nitrogenous compounds and phosphorus on the growth of toxic and non-toxic strains of Microcystis during cyanobacterial blooms. Aquat. Microb. Ecol. 2010, 61, 149–162. [Google Scholar] [CrossRef]
  86. Huang, W.; Bi, Y.; Hu, Z. Effects of fertilizer-urea on growth, photosynthetic activity and microcystins production of Microcystis aeruginosa isolated from Dianchi lake. Bull. Environ. Contam. Toxicol. 2014, 92, 514–519. [Google Scholar] [CrossRef] [PubMed]
  87. Paerl, H.W.; Xu, H.; Hall, N.S.; Zhu, G.; Qin, B.; Wu, Y.; Rossignol, K.L.; Dong, L.; McCarthy, M.J.; Joyner, A.R. Controlling cyanobacterial blooms in hypertrophic Lake Taihu, China: Will nitrogen reductions cause replacement of non-N2 fixing by N2 fixing taxa? PLoS ONE 2014, 9, e113123. [Google Scholar] [CrossRef] [PubMed]
  88. Rippka, R.; Waterbury, J.B. The synthesis of nitrogenase by non-heterocystous cyanobacteria. FEMS Microbiol. Lett. 1977, 2, 83–86. [Google Scholar] [CrossRef]
  89. Bergman, B.; Gallon, J.R.; Rai, A.N.; Stal, L.J. N2 fixation by non-heterocystous cyanobacteria. FEMS Microbiol. Rev. 1997, 19, 139–185. [Google Scholar] [CrossRef]
  90. Farnelid, H.; Bentzon-Tilia, M.; Andersson, A.F.; Bertilsson, S.; Jost, G.; Labrenz, M.; Jürgens, K.; Riemann, L. Active nitrogen-fixing heterotrophic bacteria at and below the chemocline of the central Baltic Sea. ISME J. 2013, 7, 1413–1423. [Google Scholar] [CrossRef] [PubMed]
  91. Robson, B.J.; Hamilton, D.P. Summer flow event induces a cyanobacterial bloom in a seasonal Western Australian estuary. Mar. Freshw. Res. 2003, 54, 139–151. [Google Scholar] [CrossRef]
  92. Mitrovic, S.M.; Hardwick, L.; Dorani, F. Use of flow management to mitigate cyanobacterial blooms in the Lower Darling River, Australia. J. Plankton Res. 2010, 33, 229–241. [Google Scholar] [CrossRef]
  93. Reichwaldt, E.S.; Ghadouani, A. Effects of rainfall patterns on toxic cyanobacterial blooms in a changing climate: Between simplistic scenarios and complex dynamics. Water Res. 2012, 46, 1372–1393. [Google Scholar] [CrossRef] [PubMed]
  94. Chattopadhyay, S.; Jha, M.K. Hydrological response due to projected climate variability in Haw River watershed, North Carolina, USA. Hydrol. Sci. J. 2016, 61, 495–506. [Google Scholar] [CrossRef]
  95. Cape Fear River Basinwide Water Quality Plan; Report Prepared for North Carolina Division of Water Resources; North Carolina Division of Water Resources: Raleigh, NC, USA, 2000.
  96. DWR Manual for Standard Operating Procedures. Available online: https://deq.nc.gov/about/divisions/water-resources/water-resources-data/water-sciences-home-page/intensive-survey-branch/ambient-lakes-monitoring (accessed on 22 February 2018).
  97. Intensive Survey Branch Standard Operating Procedures Manual: Physical and Chemical Monitoring; Report Prepared for North Carolina Division of Water Resources; North Carolina Division of Water Resources: Raleigh, NC, USA, 2013.
  98. Utermöhl, H. Zur Vervollkommnung der quantitativen Phytoplankton-Methodik. Mitt. Int. Ver. Theor. Angew. Limnol. 1958, 9, 1–38. [Google Scholar] [CrossRef]
  99. Standard Operating Procedures for the Collection and Analysis of Algae; Report Prepared for North Carolina Division of Water Resources; North Carolina Division of Water Resources: Raleigh, NC, USA, 2016.
  100. Wehr, J.D.; Sheath, R.G.; Kociolek, J.P. Freshwater Ecology of North America: Ecology and Classification, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2015. [Google Scholar]
  101. Lane, J.Q.; Roddam, C.; Langlois, G.W.; Kudela, R.M. Application of Solid Phase Adsorption Toxin Tracking (SPATT) for field detection of the hydrophilic phycotoxins domoic acid and saxitoxin in coastal California. Limnol. Oceanogr. Methods 2010, 8, 645–660. [Google Scholar] [CrossRef]
  102. Method 445.0 In Vitro Determination of Chlorophyll a and Pheophytin a in Marine and Freshwater Algae by Fluorescence; Report Prepared for Environmental Protection Agency; Environmental Protection Agency: Washington, DC, USA, 1997.
Figure 1. Map of Jordan Lake sampling sites. Biological, chemical and physical data were analyzed over a 6-year period (2011 to 2016) for sites A, B and C (red circles). Information over approximately 2 years (2014 to 2016) was available for an additional six sites (sites D through I, blue circles). Arrows indicate the three main rivers flowing into the lake. Map from snazzymaps.com.
Figure 1. Map of Jordan Lake sampling sites. Biological, chemical and physical data were analyzed over a 6-year period (2011 to 2016) for sites A, B and C (red circles). Information over approximately 2 years (2014 to 2016) was available for an additional six sites (sites D through I, blue circles). Arrows indicate the three main rivers flowing into the lake. Map from snazzymaps.com.
Toxins 10 00092 g001
Figure 2. Average cell densities for cyanobacteria (black circles) and microeukaryote phytoplankton (white triangles) and for Chl a concentration (white squares). Before October 2014, averages were calculated for sites A through C. After October 2014, averages were calculated for all nine sites. Note: all axes are log-transformed.
Figure 2. Average cell densities for cyanobacteria (black circles) and microeukaryote phytoplankton (white triangles) and for Chl a concentration (white squares). Before October 2014, averages were calculated for sites A through C. After October 2014, averages were calculated for all nine sites. Note: all axes are log-transformed.
Toxins 10 00092 g002
Figure 3. Cyanobacterial cell densities at each of the sampling locations (sites A through I are shown as panels A through I). Colors depict the six most abundant genera, and less abundant taxa are grouped as “Other”. Long-term sites A through C were sampled from January 2011 to December 2016. D, F, H and I were monitored from October 2014 to June 2016, while monitoring at E and G continued through December 2016. Vertical dashed lines separate years. Note, there are differing scales on the y-axes for A, C and E.
Figure 3. Cyanobacterial cell densities at each of the sampling locations (sites A through I are shown as panels A through I). Colors depict the six most abundant genera, and less abundant taxa are grouped as “Other”. Long-term sites A through C were sampled from January 2011 to December 2016. D, F, H and I were monitored from October 2014 to June 2016, while monitoring at E and G continued through December 2016. Vertical dashed lines separate years. Note, there are differing scales on the y-axes for A, C and E.
Toxins 10 00092 g003
Figure 4. (A) MDS plot based on Bray–Curtis similarities for cyanobacterial communities as a non-metric multi-dimensional scaling (nMDS) plot by season (data from all years and stations combined). (B) Relative changes in cyanobacterial community composition shown along a month-to-month trajectory (site A in 2015). Stress values are reported in the top right corner of each plot.
Figure 4. (A) MDS plot based on Bray–Curtis similarities for cyanobacterial communities as a non-metric multi-dimensional scaling (nMDS) plot by season (data from all years and stations combined). (B) Relative changes in cyanobacterial community composition shown along a month-to-month trajectory (site A in 2015). Stress values are reported in the top right corner of each plot.
Toxins 10 00092 g004
Figure 5. Changes in (A) temperature, (B) NOx, (C) NH3 concentration and (D) Total Kjeldahl nitrogen (TKN):TP ratio averaged for each sampling event. Standard error bars are included. Vertical dashed lines separate years.
Figure 5. Changes in (A) temperature, (B) NOx, (C) NH3 concentration and (D) Total Kjeldahl nitrogen (TKN):TP ratio averaged for each sampling event. Standard error bars are included. Vertical dashed lines separate years.
Toxins 10 00092 g005
Figure 6. MDS plot showing Euclidian distances for environmental fingerprints by season. Parameters included in these analyses were temperature, NOx, NH3, TKN, TP and DO concentrations, TKN:TP ratios, pH levels and turbidity. The stress value is reported in the top right corner.
Figure 6. MDS plot showing Euclidian distances for environmental fingerprints by season. Parameters included in these analyses were temperature, NOx, NH3, TKN, TP and DO concentrations, TKN:TP ratios, pH levels and turbidity. The stress value is reported in the top right corner.
Toxins 10 00092 g006
Figure 7. Weekly averages of meteorological and hydrological parameters: (A) precipitation; (B) wind speed; (C) PAR; (D) overall river flow. Vertical dashed lines separate years.
Figure 7. Weekly averages of meteorological and hydrological parameters: (A) precipitation; (B) wind speed; (C) PAR; (D) overall river flow. Vertical dashed lines separate years.
Toxins 10 00092 g007
Figure 8. SPATT toxin values (columns) and toxin concentrations based on grab sampling (symbols) for (A) MCY, (B) CYN and (C) ANA for site E. SPATT toxin concentrations in ng toxin (g resin−1) d−1 are shown at the half-point of each deployment period. Grab samples are represented as µg toxin L−1 (filled symbols). Empty symbols along the x-axis indicate when toxin values fell below the LDL for each ELISA kit (LDLs shown as horizontal dashed lines originating from the secondary y-axes).
Figure 8. SPATT toxin values (columns) and toxin concentrations based on grab sampling (symbols) for (A) MCY, (B) CYN and (C) ANA for site E. SPATT toxin concentrations in ng toxin (g resin−1) d−1 are shown at the half-point of each deployment period. Grab samples are represented as µg toxin L−1 (filled symbols). Empty symbols along the x-axis indicate when toxin values fell below the LDL for each ELISA kit (LDLs shown as horizontal dashed lines originating from the secondary y-axes).
Toxins 10 00092 g008
Figure 9. Seasonal averages for dissolved toxin concentrations based on (A) discrete sample analyses for sites A through G and (B) in-situ tracking (SPATTs) of MCY, ANA and CYN for sites E and G. Standard error bars are included whenever multiple samples tested positive. Note there are differing log-scales on the y-axes.
Figure 9. Seasonal averages for dissolved toxin concentrations based on (A) discrete sample analyses for sites A through G and (B) in-situ tracking (SPATTs) of MCY, ANA and CYN for sites E and G. Standard error bars are included whenever multiple samples tested positive. Note there are differing log-scales on the y-axes.
Toxins 10 00092 g009
Table 1. Latitude (Lat) and longitude (Long) for the nine sampling sites across Jordan Lake. Sampling was conducted on a monthly basis with more frequent biweekly monitoring during months with higher bloom activity (May through September). Included are the Site ID and Division of Water Resources (DWR) site names. n = number of sampling time points per site. Depth describes the average water column depth.
Table 1. Latitude (Lat) and longitude (Long) for the nine sampling sites across Jordan Lake. Sampling was conducted on a monthly basis with more frequent biweekly monitoring during months with higher bloom activity (May through September). Included are the Site ID and Division of Water Resources (DWR) site names. n = number of sampling time points per site. Depth describes the average water column depth.
Site IDDWR IDLat (° N)Long (° W)FromTonDepth (m)
ACPF086C35.79479.004January 2011December 2016985.06
BCPF087D35.74279.021January 2011December 2016967.67
CCPF055C35.68779.083January 2011December 2016955.87
DCPF086CUPS35.83779.001October 2014June 2016261.47
ECPF086C35.82578.998October 2014December 2016342.92
FCPF081A1B35.83678.976October 2014June 2016251.82
GCPF081A1C35.81578.983October 2014December 2016323.16
HCPF055C135.69979.082October 2014June 2016272.25
ICPF055C635.68279.078October 2014June 2016278.50
Table 2. Percentage (%) of samples that tested positive for varying toxins using discrete sampling and the Solid Phase Adsorption Toxin Tracking (SPATTs) approach. Average concentrations (Ave) for dissolved (Diss) toxins are shown as µg L−1 (values below LDL were not included when calculating the average for each toxin), and for SPATT as ng toxin (g resin)1 d−1. n = number of samples tested. LDL = low detection limit of Enzyme-Linked Immunosorbent Assay (ELISA) detection method; BDL = below detection limit of ELISA test; MCY = microcystin; ANA = anatoxin-a; CYN = cylindrospermopsin; BMAA = β-N-methylamino-l-alanine; STX = saxitoxin.
Table 2. Percentage (%) of samples that tested positive for varying toxins using discrete sampling and the Solid Phase Adsorption Toxin Tracking (SPATTs) approach. Average concentrations (Ave) for dissolved (Diss) toxins are shown as µg L−1 (values below LDL were not included when calculating the average for each toxin), and for SPATT as ng toxin (g resin)1 d−1. n = number of samples tested. LDL = low detection limit of Enzyme-Linked Immunosorbent Assay (ELISA) detection method; BDL = below detection limit of ELISA test; MCY = microcystin; ANA = anatoxin-a; CYN = cylindrospermopsin; BMAA = β-N-methylamino-l-alanine; STX = saxitoxin.
ToxinSample TypeAveRangePositive (%)nLDL
MCYDiss0.37BDL—1.9815650.10
SPATT39.49BDL—347.459224
ANADiss0.2BDL—0.6857690.10
SPATT3.970.31—13.2810023
CYNDiss0.27BDL—0.8310630.04
SPATT0.05BDL—0.051324
BMAADiss10.75BDL—23.4514644.00
STXDissBDLBDL0400.015
Table 3. Results from correlation analyses. Pearson’s correlation coefficients (r) shown in bold are significant at p < 0.05. Cyano = cyanobacteria; Microphyto = microeukaryote phytoplankton; Diss = dissolved; DO = dissolved oxygen; Chl a = chlorophyll a.
Table 3. Results from correlation analyses. Pearson’s correlation coefficients (r) shown in bold are significant at p < 0.05. Cyano = cyanobacteria; Microphyto = microeukaryote phytoplankton; Diss = dissolved; DO = dissolved oxygen; Chl a = chlorophyll a.
Chl a (µg L−1)Cyano (Cells mL−1)Cyano (mm3 m−3)Microphyto(Cells mL−1)Microphyto(mm3 m−3)MCY (ng mL−1)ANA (ng mL−1)
Cyano (cells mL−1)0.555
Cyano (mm3 m−3)0.5150.851
Microphyto (cells mL−1)0.3580.3670.314
Microphyto (mm3 m−3)0.4470.4690.3640.329
Diss MCY (ng mL−1)0.1010.2720.3130.4080.055
Diss ANA (ng mL−1)−0.1660.1250.132−0.0080.2770.270
NH3 (mg L−1)0.3520.2570.2290.1540.210−0.044−0.135
NOx (mg L−1)0.5400.5860.4990.2950.349−0.102−0.106
TKN (mg L−1)0.7500.4980.4920.3140.2820.147−0.016
TP (mg L−1)0.217−0.0220.0610.049−0.031−0.023−0.183
TKN:TP0.1100.3060.2290.0890.1580.2000.326
Turbidity (NTU)0.1640.1120.0020.0680.105−0.075−0.230
Temp (°C)0.3220.5240.4700.1750.2380.1500.150
DO (mg L−1)−0.0180.3090.229−0.0810.101−0.017−0.125
pH0.4410.4580.4370.2000.2310.2550.043
Table 4. Results from multiple regression analyses. Adj. R2 = adjusted coefficient of determination at p < 0.05. Individual t-statistics are listed in parenthesis. Turb = turbidity. Degrees of freedom for each regression = 441.
Table 4. Results from multiple regression analyses. Adj. R2 = adjusted coefficient of determination at p < 0.05. Individual t-statistics are listed in parenthesis. Turb = turbidity. Degrees of freedom for each regression = 441.
Adj. R2Parameter 1Parameter 2Parameter 3Parameter 4n
Chl a0.680NH3 (−8.61)NOx (−10.50)TKN (19.99)DO (5.66)427
Cyano (cells mL−1)0.521NOx (−9.10)TKN (5.43)DO (−5.56)pH (8.62)446
Cyano (mm3 m−3)0.406NOx (−8.79)TKN (6.31)pH (6.19) 452
Microphyto (cells mL−1)0.132NOx (−4.42)TKN (5.05) 454
Microphyto (mm3 m−3)0.189NH3 (−3.12)NOx (−3.99)TKN (5.34)Turb (−3.71)454
Table 5. Summary of reported cyanotoxins in North Carolina (NC) water bodies (based on discrete sampling). Average (Ave) concentrations denoted with an asterisk are approximated from published figures. BDL = below detection limit. ND = no data provided. n = replicates. Res = Reservoir. * = average value approximated from graph.
Table 5. Summary of reported cyanotoxins in North Carolina (NC) water bodies (based on discrete sampling). Average (Ave) concentrations denoted with an asterisk are approximated from published figures. BDL = below detection limit. ND = no data provided. n = replicates. Res = Reservoir. * = average value approximated from graph.
LocationMonth/Year Water BodyToxinAve (µg L−1)Range (µg L−1)nMethodReference
ApexAugust 2015–December 2016Jordan LakeMCY0.06BDL—1.9865ELISAThis Study
CYN0.02BDL—0.8363
ANA0.11BDL—0.6869
BMAA9.32BDL—23.4564
PiedmontJune 2002–August 2002Jordan LakeMCY0.20 *ND6ELISA[52]
Kerr Scott ResMCY0.30 *ND6
Tuckertown ResMCY0.12 *ND6
Oak Hollow LakeMCY0.10 *ND6
Falls LakeMCY0.22 *ND6
Narrows ResMCY0.15 *ND6
Lake RhodhissMCY0.25 *ND6
Lake MichieMCY0.15 *ND6
High Rock LakeMCY0.05 *ND6
Lake TilleryMCY0.35 *ND6
High Point LakeMCY0.12 *ND6
PiedmontJune 2011–September 2012City LakeMCY0.22BDL—0.316ELISA[53]
Oak Hollow LakeMCY0.21BDL—0.264
Randleman ResMCY0.17BDL—0.187
Lake MackintoshMCY0.17BDL—0.175
WatervilleOctober 2007Waterville ResMCY824.3ND3LC-MS[54]
StatewideJune 2007–July 2007Lake LeeMCY0.210.17—0.242ELISA[47]
Lake RhodhissMCY0.14BDL—0.142
STX0.03BDL—0.032
Lake OrangeMCY0.540.541
Lake FisherMCY0.170.1471
High Rock LakeMCY0.520.521
Lake TownsendMCY0.160.161
Falls LakeMCY0.280.281
Lake HickoryMCY0.160.161
Beaverdam LakeMCY0.230.231
Graham-Mebane LakeMCY0.110.111

Share and Cite

MDPI and ACS Style

Wiltsie, D.; Schnetzer, A.; Green, J.; Vander Borgh, M.; Fensin, E. Algal Blooms and Cyanotoxins in Jordan Lake, North Carolina. Toxins 2018, 10, 92. https://doi.org/10.3390/toxins10020092

AMA Style

Wiltsie D, Schnetzer A, Green J, Vander Borgh M, Fensin E. Algal Blooms and Cyanotoxins in Jordan Lake, North Carolina. Toxins. 2018; 10(2):92. https://doi.org/10.3390/toxins10020092

Chicago/Turabian Style

Wiltsie, Daniel, Astrid Schnetzer, Jason Green, Mark Vander Borgh, and Elizabeth Fensin. 2018. "Algal Blooms and Cyanotoxins in Jordan Lake, North Carolina" Toxins 10, no. 2: 92. https://doi.org/10.3390/toxins10020092

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