Next Article in Journal
Overview of Pharmacokinetics and Liver Toxicities of Radix Polygoni Multiflori
Next Article in Special Issue
Kinetics of Microcystin-LR Removal in a Real Lake Water by UV/H2O2 Treatment and Analysis of Specific Energy Consumption
Previous Article in Journal
Multiplex Immunoassay Techniques for On-Site Detection of Security Sensitive Toxins
Previous Article in Special Issue
Chitosan as a Coagulant to Remove Cyanobacteria Can Cause Microcystin Release
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity Assessment of Toxic Cyanobacterial Blooms during Oxidation

1
Department of Civil, Geological, and Mining Engineering, Polytechnique Montréal, Montréal, QC H3T 1J4, Canada
2
Department of Biological Science, Université de Montréal, Montréal, QC H2V 0B3, Canada
3
Water Research Australia (WaterRA), Adelaide, SA 5001, Australia
4
BGA Innovation Hub and Water Research Centre, School of Civil and Environmental Engineering, University of New South Wales (UNSW), Sydney, NSW 2052, Australia
*
Author to whom correspondence should be addressed.
Submission received: 18 September 2020 / Revised: 28 October 2020 / Accepted: 18 November 2020 / Published: 20 November 2020
(This article belongs to the Special Issue Removal of Cyanobacteria and Cyanotoxins in Waters)

Abstract

:
Fresh-water sources of drinking water are experiencing toxic cyanobacterial blooms more frequently. Chemical oxidation is a common approach to treat cyanobacteria and their toxins. This study systematically investigates the bacterial/cyanobacterial community following chemical oxidation (Cl2, KMnO4, O3, H2O2) using high throughput sequencing. Raw water results from high throughput sequencing show that Proteobacteria, Actinobacteria, Cyanobacteria and Bacteroidetes were the most abundant phyla. Dolichospermum, Synechococcus, Microcystis and Nostoc were the most dominant genera. In terms of species, Dolichospermum sp.90 and Microcystis aeruginosa were the most abundant species at the beginning and end of the sampling, respectively. A comparison between the results of high throughput sequencing and taxonomic cell counts highlighted the robustness of high throughput sequencing to thoroughly reveal a wide diversity of bacterial and cyanobacterial communities. Principal component analysis of the oxidation samples results showed a progressive shift in the composition of bacterial/cyanobacterial communities following soft-chlorination with increasing common exposure units (CTs) (0–3.8 mg·min/L). Close cyanobacterial community composition (Dolichospermum dominant genus) was observed following low chlorine and mid-KMnO4 (287.7 mg·min/L) exposure. Our results showed that some toxin producing species may persist after oxidation whether they were dominant species or not. Relative persistence of Dolichospermum sp.90 was observed following soft-chlorination (0.2–0.6 mg/L) and permanganate (5 mg/L) oxidation with increasing oxidant exposure. Pre-oxidation using H2O2 (10 mg/L and one day contact time) caused a clear decrease in the relative abundance of all the taxa and some species including the toxin producing taxa. These observations suggest selectivity of H2O2 to provide an efficient barrier against toxin producing cyanobacteria entering a water treatment plant.
Key Contribution: This study provides a diversity assessment of cyanobacterial bloom following oxidation. Changes, shifts and oxidant persistence potential of cyanobacterial community during oxidation were quantified and reported.

1. Introduction

The occurrence of cyanobacterial blooms in fresh-water bodies has been enhanced due to eutrophication and temperature increases [1,2]. Cyanobacterial blooms may produce and release taste and odour compounds as well as cyanotoxins into water bodies. More than 40 species of cyanobacteria are known as potentially toxic species [3,4]. Microcystin (MC), anatoxin (ATX-a), saxitoxin (STX), cylindrospermopsin (CYN) and β-Methylamino-L-alanine (BMAA) are the five major groups of cyanotoxins due to their toxicity and frequent occurrence around the world [5].
Conventional drinking water treatment plants are often challenged by the removal of cyanobacteria and cyanotoxins [6,7,8,9] including both low risk (low cell numbers in the source) and high-risk water treatment plants (high cell numbers in the source). Several studies have evaluated the removal of the cyanobacteria and their harmful metabolites using different oxidants such as chlorine, ozone and, potassium permanganate [10,11,12,13,14,15]. Pre-oxidation dampens the cyanobacterial shock before entering the drinking water treatment plant and can limit the accumulation of cyanobacteria within the plant. Chlorine and ozone are also used as primary disinfectants providing an additional oxidation barrier after the filtration to remove cyanobacterial harmful metabolites (cyanotoxins) [16,17].
Dynamic and complex behaviour of cyanobacterial blooms, including cyanotoxin production and release under various environmental conditions represents a treatment challenge. The oxidation efficiency of cyanobacteria blooms may vary according to water quality parameters (such as pH, Dissolved Organic Carbon (DOC) and the presence of other bacterial communities), cyanobacterial community shape, potential agglomeration, and the growth phase [14,18,19,20]. Cyanobacteria treatment efficiency improvement requires an understanding of the cyanobacterial composition structure in response to treatment processes [21]. Molecular methods have been used to study the fate of the microbial community, including cyanobacteria, within different conditions [5,22]. Molecular methods can overcome the challenges of microscopic cell counts such as time and qualified person requirements, as well as changes in biovolumes during analyses [23,24,25].
Molecular methods such as high throughput sequencing have been deployed to study cyanobacterial communities and identify cyanotoxin biosynthesis genes [5,26,27,28,29,30,31,32,33,34]. Diversity of the cyanobacterial community and toxigenic cyanobacteria was assessed based on the Operational Taxonomy Units (OTUs) derived from 16S rRNA gene amplification (metabarcoding) [5]. Limited studies have applied high throughput sequencing to monitor the fate of cyanobacteria during the treatment processes. Xu, Pei [35] studied the microbial community of the sludge in six different drinking water treatment plants using high throughput sequencing. Results showed that cyanobacteria were the most dominant phylum in two treatment plants with a higher level of nutrients in raw water. Planktothrix, Microcystis and Cyindrospermopsis were the most abundant genera and were positively correlated with the nutrient levels in raw water. Pei, Xu [36] used 16S rRNA sequencing to study the shifts in the microbial community in clarifier sludge following coagulation by FeCl3, AlCl3 and PAFC (Polyaluminium Ferric Chloride). Results revealed selective removal of the different bacterial species, as the relative abundance of the Microcystis, Rhodobacter, Phenylobacterium and Hydrogenophaga decreased in AlCl3 sludge compare to the FeCl3 and PAFC. Lower Microcystis abundance could be related to high Al toxicity or large and high-density floc in FeCl3 and PAFC, which plays a protective role for microorganisms [36]. Lusty and Gobler (2020) used 16S rRNA to evaluate the mitigation of cyanobacterial blooms using H2O2. Results showed relative persistence of Cyanobium and Cylindrospermopsis to a moderate H2O2 dose (4 mg/L); Plankthotrix and Microcystis (abundant genus) were the most sensitive genera, respectively [37].
High throughput sequencing has been widely applied to study bacterial/cyanobacterial communities. Fewer studies focused on the diversity of bacterial/cyanobacterial communities during water treatment processes. However, no study has focused on the cyanobacterial community following chemical oxidation using high throughput sequencing. Understanding the shifts and the potential selective persistence in cyanobacterial communities following oxidation is important for choosing an efficient oxidant. Thus, the objective of this study was to assess the structural composition of the cyanobacteria community following oxidation (with Cl2, O3, KMnO4, H2O2) using high throughput metagenomic shotgun sequencing over the seasonal bloom period.

2. Results and Discussion

2.1. Cyanobacterial Bloom Characteristics Throughout Sampling

The cyanobacteria bloom samples were collected on 5 days of the bloom period (1, 13, 15, 21 and 29 August 2018) from Missisquoi Bay (Lake Champlain) close to the water intake of the drinking water treatment plant. Cyanobacterial bloom characteristics are presented in the Table 1. DOC (Dissolved Organic Carbon) and pH did not demonstrate considerable variation throughout the sampling period. Total cell counts and biovolumes follow the same decreasing and increasing trend during the sampling period. For both parameters, the highest values were found at the beginning of the sampling period (1 August), followed by a significant drop on 13 August, where cyanobacterial cell counts decreased from 3.3 × 105 to 7.8 × 104 cells/mL and the biovolumes from 30.6 mm3/L to 4.6 mm3/L. The second drop was observed between 15 August and 21 August, from 1.4 × 105 to 6.8 × 104 cells/mL for cell count and from 9.4 to 0.3 mm3/L for the biovolume. According to algal cell abundance descriptors (biovolume and cell counts) during the bloom period, the main peaks of cyanobacterial bloom occurred on 1 August and to a lower extent on 15 August. The observed cell count exceeds the alert level of the 6.5 × 104 cells/mL for drinking water treatment plants [38], except for the 29 August with 5.4 × 104 cells/mL.

2.2. Variation of the Cyanobacterial Bloom Composition

The diversity and community variation of bacterial and cyanobacterial communities during the bloom sampling period were studied using comparative metagenomics reads levels of phylum, order and genus. The number of reads for taxonomic data was normalized by relative abundance (Figure 1).
Analysis of relative abundance of the bacterial community at the phylum level, based on high throughput sequencing data, showed that Proteobacteria, Actinobacteria, Cyanobacteria, Bacteroidetes, Firmicutes and Verrucomicrobia were the six most abundant phyla throughout the sampling period (Figure 1a). As expected, Proteobacteria was by far the most abundant phylum at the beginning of the sampling period (1 August) as has been observed by Pei et al. (2017) [36]. The high relative abundance of Proteobacteria, especially in the beginning of sampling and Bacteroidetes at the end of the sampling period, may indicate contamination of the sampling point with human/animal-associated fecal markers [39,40,41]. For the rest of the sampling dates (13, 15 and 21 August), Proteobacteria remained the predominant phylum, but at a lower extent than what was observed in the first and last days of sampling. The cyanobacteria phylum accounts for 5 to 10% of total relative abundance assigned to the phylum level in all samples. The cyanobacteria relative abundance started at 5% of the phyla, followed by an increase in the middle of the sampling (13 August and 15 August) to 10%. By the end of August, the cyanobacterial contribution in the whole bacterial community decreased to 5%.
At the order level, the cyanobacterial community was dominated by members of the Chroococcales, Nostocales and Oscillatoriales during the cyanobacterial bloom period (from 1 August and 29 August) (Figure 1b). The relative abundance of the Nostocales and Chroococcales varied between the different dates and even within the same day (15 August). On the other hand, Oscillatoriales relative abundance remained steady (Figure 1b).
Analysis at the genus level showed that within the Chroococcales order, the predominant genera were Microcystis and Synechococcus. Synechococcus was the dominant genus within Chroococcales until 21 August of the sampling period, followed by Microcystis, which became the dominant genus on 29 August (Figure 1c). The predominant genera in the Oscillatoriales and Nostocales orders were Oscillatoria and Dolichospermum (formerly known as Anabaena), respectively (Figure 1c). Depending on the sampling date, Synechococcus, Microcystis and Dolichospermum were the predominant genera during the bloom period. Our observations are consistent with short/long term investigations of cyanobacterial bloom dynamics using taxonomic cell count and metagenomics [14,21,42,43,44]. Microcystis and Dolichospermum share a spatio-temporal niche during blooms and often are dominant within the cyanobacterial community, but have distinct environmental preferences; for example, they have different responses to nutrients [45,46,47]. During this study, the abundant genus shifted from Nostocales members, dominated by Dolichospermum for the first three weeks, to the Microcystis genus by the end of the sampling period (29 August). The shift to Microcystis dominance can probably be attributed to species-specific associations between Microcystis spp. and associated bacteria referred to as the epibiont phenomenon. Some species within Proteobacteria, Bacteroidetes and other phyla can shelter in the Microcystis mucilage to avoid being grazed [48]. These species play essential roles in enhancing the environmental adaptation of Microcystis within the cyanobacterial bloom, like maintaining redox balance and coping with oxidative stress [48]. Our results are in accordance with the previously reported results that Proteobacteria, Bacteriodetes tend to dominate in the Microcystis mucilage (blooms and culture-dependent studies) [49,50]. Moreover, allelopathy may influence the successional dominance of Microcystis and Dolichospermum within aquatic systems, whereby organisms produce bioactive compounds (allelochemicals) in the environment to positively or negatively influence the growth of neighboring species [51]. The effects of these allelochemicals on cyanobacterial community distribution within the environment are connected with nutrient availability and environmental conditions [52,53,54]. Competition experiments between toxic Microcystis and Dolichospermum strains, based on the lab coculture-dependent method, showed that Microcystis significantly inhibited the growth of Dolichospermum, whereas the effects of Dolichospermum on Microcystis were minimal [53,55]. Dolichospermum biovolume and biomass were sharply reduced after exposure to Microcystis strains [56]. Further investigation is required to explain the interaction and succession of Microcystis and Dolichospermum in cyanobacterial blooms.
An interesting insight was gained considering within-day changes in bacterial community composition (morning—AM and afternoon—PM) for 15 August (Figure 1). Based on the relative abundance analysis, a shift in bacterial composition at the phylum level was observed between the two samples on 15 August. In the afternoon, Proteobacteria and Actinobacteria increased, while Bacteroidetes and Cyanobacteria decreased as compared to their initial composition at the start of the sampling in the morning. At the order level, the most abundant order shifted from Nostocales in the morning to Chroococcolas in the afternoon. At the genus level, the changes were more evident for Synechococcus (increase) and Dolichospermum (decrease). These within-day changes in bacterial composition reflect the variation in stratification and mixing patterns of the cyanobacterial community in Missisquoi bay, as shown by Ndong et al. (2014) Ndong, Bird [57]. These findings highlight the importance of timing the sample collection by considering the stratification variation in the diel cycle (morning, noon, afternoon).

2.3. Impact of Oxidation on Cyanobacterial Diversity

The impact of oxidation on cyanobacterial diversity was assessed by using the impact of increased CT (the terminology used in water treatment is CT, which represents the product of the oxidant concentration (C) and contact time (T) to inactivate microorganisms) on two different samples, first on 1 August, dominated by Dolichospermum and the second on 29 August dominated by the Microcystis genus. The oxidants considered Cl2, KMnO4, O3, H2O2 differ widely in their mode of action and their persistence [58]. As a result, CT values may vary over three orders of magnitude from 0.1 mg·min/L Cl2 and a maximum of 7035 mg·min/L for H2O2.
Using non-normalized data is a source of bias in statistical analyses of the effects of oxidants on the bacterial/cyanobacterial communities. Thus, the data were normalized prior to analysis. Using the common exposure unit (CT) is not representative to compare the effects of different oxidants simultaneously. A normalized oxidant exposure (relative CT) is used to compare different oxidants among each other (Relative CT for each oxidant; max CT = 1 and min CT = 0, the exposure points in between were calculated accordingly); for example, max CT obtained for H2O2, on 1 August, is 7035 mg·min/L and is considered as relative CT = 1, the control condition is considered as relative CT = 0.
The doses (concentration) and contact time (T) were selected using the product CT, which is the foundation of disinfection and oxidation in drinking water. The choice of CT also considered prior evidence showing differences in species’ resistance to the oxidant exposure (CT) [6,10,12,14,15,16,18,42,59,60]. However, all previously reported observations were based on taxonomic cell counts and not on high throughput sequencing. The applied dosages reflect drinking water industry practices using regulated treatment processes.

2.3.1. Cyanobacterial Composition

To assess the cyanobacterial bloom composition following oxidation using KMnO4, Cl2 and H2O2, samples were taken on 1 August (Figure 2) and 29 August (Figure S1). Dolichospermum and Microcystis were the dominant genera on 1 August and 29 August, respectively.
Figure 2 shows the dissimilarity among groups of bacterial communities following oxidation at different exposures (relative CT) using Principal Component Analysis (PCA) on 1 August. Samples that appear more closely together within a PCA are assumed to be more similar in bacterial and cyanobacterial composition (Figure 2). For the bacterial community, principal axis one and principal axis two for PCA represent 51.5 and 31.4% of the variation among the samples, respectively. High-relative CT exerted clustering of two groups, the first one includes samples after Cl2 and KMnO4 oxidation and the second one encompasses KMnO4 and H2O2 oxidation. Chlorinated samples showed a clear progressive shift (as the relative CT increased) in the bacterial community.
On the other hand, the KMnO4 and H2O2 induce large shifts as relative CT increases. Bacterial and cyanobacterial composition similarity following the different oxidation on 29 August (Microcystis) is presented in Figure S1. Like the 1 August result, a large shift in the bacterial community composition following H2O2 oxidation was observed for 29 August. On the other hand, KMnO4 results showed similar bacterial composition on 29 August.
For the cyanobacterial community (1 August), principal axis one, and principal axis two for PCA represent 74.4 and 19.8% of the variation among the samples, respectively (Figure 2b). A clear trend in the cyanobacterial composition variation is observed following chlorination. Cyanobacterial composition after exposure to Cl2 (low relative CT < 0.05), KMnO4 (low and high relative CT), H2O2 (low relative CT = 0.08) oxidation clustered in the same group with the control condition. High-relative (CT = 1) using H2O2 and (mid-relative = 0.53) KMnO4 are grouped in a distinct second cluster, revealing similar cyanobacterial assemblages. This is in contrast with cyanobacterial/bacterial community trends following chlorination that display a progressive shift. Differences between the observed trends for the three oxidants were expected because of different mechanisms of actions and kinetics, persistence and selectivity. Cell count based studies show progressive shifts following oxidation [12,14]. Moreover, selective cyanobacteria oxidation has been demonstrated in the lab [11,60], in the field [37,61], and in drinking water treatment plants [8,62].
Our results enlighten the different oxidation impacts on the bacterial community between the Cl2 and H2O2 oxidation. For 29 August, which was dominated by the Microcystis genus, KMnO4 results revealed similar cyanobacterial composition. However, a large variation in cyanobacterial composition following H2O2 was observed (e.g., the samples from the 1 August, Figure 3).
Redundancy analysis (RDA) was used to explore the correlation between the relative abundance of nine dominant cyanobacterial genera (60% of the cyanobacterial community) observed in 1 August or 29 August samples and the different oxidant exposures (Figure 3). The relationship between the CT and diversity indexes (Shannon and Chao1) varies between the two dates. On 1 August, both diversity indices increase with CT, while on 29 August only the Richness index (Chao1) increased with CT. In general, the diversity of the cyanobacterial community increased due to oxidation of the dominated genus, as discussed further in detail for each oxidant.
For 1 August (Dolichospermum dominant), RDA analysis establishes an inverse correlation between the Dolichospermum genus and a wide range of chlorine exposures with relative CT > 0.04. Chlorination appears to have a lesser impact on the Dolichospermum genus than the Microcystis, possibly because of its abundance. On the other hand, no such relationship is observed between Dolichospermum and KMnO4 exposure (relative CT = 1) and high H2O2 exposure (relative CT = 1), suggesting that these oxidants at these exposure levels had a negative impact on the Dolichospermum genus persistence within the community. Furthermore, for KMnO4 at relative CT > 0.53, Microcystis, Synechococcus, and Leptolyngbya are less impacted by oxidation. Figure 3b shows the effect of the different oxidants on the cyanobacterial community of 29 August, when Microcystis genus is dominant. An immediate shift from the control is observed for all oxidants and for all relative CT. Unlike 1 August, the Shannon diversity is no longer correlated with Dolichospermum. These differences could reflect the morphological differences between the dominant genera (Microcystis vs. Dolichospermum) as a unicellular aggregate that are less resistant to oxidation as compared to filaments [63]. Chlorine and KMnO4 exposures had a low impact on Microcystis and Dolichospermum. On the other hand, any H2O2 exposure causes a reduction in Dolichospermum and Microcystis, more so on 29 August. Removal of Microcystis using H2O2 was shown by Lusty and Gobler. (2020) and is in accordance with our observations of diverging correlation between Chao1 and Microcystis [37].
Figure S2 shows the cyanobacterial composition following oxidation at genus level. Results show no significant variation in the relative abundance of Dolichospermum and Microcystis following Cl2, KMnO4 oxidation. Thus, no relative persistence to oxidation was observed at the genus level. On the other hand, high H2O2 exposure caused a decline in both Dolichospermum and Microcystis genus, which demonstrates the effect of high H2O2 exposure on cyanobacterial removal.

2.3.2. Effect of Oxidation on Cyanobacterial Community Richness and Diversity

The diversity of the cyanobacterial composition following oxidation (1 August) at the species level is presented in Figure 4 (top 25 most abundant species). Dolichospermum sp.90 was the dominant species in control conditions, and its relative abundance increases following chlorination and decreases after KMnO4 and H2O2. Although no trends in relative abundance were seen at the genus level (Figure S2), a similar trend can be seen for the three species present Dolichospermum sp.90, Dolichospermum cylindrica and, Dolichospermum sp. PCC7108. Regardless of the oxidant Dolichospermum sp.90 remains the dominant species. Following the KMnO4 and H2O2 oxidation, Dolichospermum sp.90 was still the dominant species; its relative abundance declined compared to the control. The relative abundance of Microcystis aeruginosa is not impacted by chlorine but increases with KMnO4 and H2O2 exposure, confirming selective removal of Microcystis showed by Lusty and Gobler. (2020) [37].
In the 29 August samples, Microcystis aeruginosa was the most abundant species (Figure S3); its relative abundance increased after chlorination and KMnO4 oxidation (as the chlorine exposure increased) and decreased after H2O2 oxidation as compared to the control condition. Despite the very low relative abundance of Dolichospermum sp.90, Dolichospermum cylindrica and Dolichospermum sp. PCC7108, similar trends were observed on 1 August (when Dolichospermum was abundant).
The community richness and diversity indices for each treatment for 1 August and 29 August samples are illustrated in Figure 5 and Figure S4. Shannon and Chao1 show a small decline following chlorination in comparison with the control, while they increase slightly following KMnO4 exposure. The total cell numbers following KMnO4 decreased by up to 63% for high KMnO4 exposure (Figure S7). A remarkable decrease in richness is observed at high relative CT of H2O2, while the diversity (Shannon index) increases. The decline in the richness index could be the result of some less abundant species no longer identified. Indeed, total cell counts following H2O2 relative CT = 1, decreased by more than 50% (Figure S6). The same trends in the alpha diversity measurements are observed in the last week of the sampling (29 August), where the Microcysits were the most abundant genus (Figure S4).

2.4. Cyanobacterial Community Assessment Following Oxidation; Longitudinal Study

The induced changes of cyanobacterial composition structure following the oxidation (Cl2, KMnO4, O3 and, H2O2) are assessed separately. The analysis was performed at the genus and species level.

2.4.1. Chlorination (Cl2)

The chlorination experiments were conducted on 1 August (Figure 6a), and 29 August (Figure 6b). In the first chlorination trial (1 August), the abundant genera were Dolichospermum and Nostoc, representing approximatively 20% and 10% of the cyanobacterial community, respectively (Figure 6a). Figure 6 shows a limited effect of chlorination on the relative abundance of all cyanobacteria genera except Synechococcus. In addition, taxonomic cell counts show decrease (up to 30%) in total cyanobacteria cell counts for the trial on 1 August and limited variation for the 29 August trial (15% variation) (Figure S5). In terms of the species, Dolichospermum sp.90 and Dolichospermum cylindrica were dominant for the 1 August trial. For the 29 August trial, Microcystis aeruginosa was the dominant species, followed by Dolichospermum sp.90 (Figure S6). In all trials, as chlorination exposure increased, the relative abundance of the abundant species, either Dolichospermum sp.90 or Microcystis aeruginosa, increased slightly. Moreover, the relative abundance of Dolichospermum sp.90 increased as Microcystis aeruginosa did on 29 August. Chlorination results show that Dolichospermum species and Microcystis aeruginosa are relatively more persistent than the other species.

2.4.2. Potassium Permanganate (KMnO4)

The relative abundance of the different genera following oxidation using permanganate shows limited variation for both KMnO4 tests (Figure 7). Total cell counts decreased (up to 57%) in the first trial (1 August) and remained stable in the second KMnO4 trial (less than 1% variation) (29 August). Total cell counts following the first KMnO4 trial showed a decrease at 278 mg·min/L (Figure S7). The relative abundance of the Dolichospermum sp.90 increased slightly, whether it was the abundant species or not, suggesting the relative persistence of the Dolichospermum sp.90 during KMnO4 oxidation (Figure S8).

2.4.3. Ozonation (O3)

The first and second ozonation trials were performed on 15 August (Dolichospermum most abundant genus) and 21 August (Synechococcus most abundant genus). Cyanobacterial community results following ozonation in the second trial (at genus level—Figure 8) showed a decline of the relative abundance of Synechococcus, followed by an increase for Microcystis, as compared to the control condition. However, no significant variation was observed for the relative abundance of the different genera in the first ozonation trials. Furthermore, total cyanobacteria cell counts revealed no significant change for both ozonation tests (up to 15%) (Figure S9). In the control condition of the 15 August ozonation trial, Microcystis aeruginosa was the dominant species, followed by Dolichospermum sp.90 (Figure S10). Dolichospermum sp.90 was not the dominant species, but it remained intact following ozonation. Although Cyanobioum gracil and Synechococcus sp. were the dominant species in the 21 August trial control, Microcystis aeruginosa became the dominant species following ozonation. At the low dosage applied, only secondary oxidation radical by-products are likely to react with cyanobacterial cells. Under these soft ozonation conditions, unicellular Cyanobium gracil and Synechococcus sp. cells were more susceptible than Microcystis aeruginosa.

2.4.4. Hydrogen Peroxide (H2O2)

The Dolichospermum/Dolichospermum sp.90 and Microcystis/Microcystis aeruginosa were the most abundant genus/species on 1 August and 29 August H2O2 oxidation, respectively. The relative abundance of the Dolichospermum declined following H2O2 exposure of CT = 7035 mg·min/L. Microcystis relative abundance decreases by more than 10% at 8442 mg·min/L exposure of H2O2 on 29 August experiment (Figure 9). At the same H2O2 exposure, Dolichospermum decreased by 5%. Total cyanobacteria cell counts declined following H2O2 oxidation for both dates: 52% decrease on 1 August and 49% decrease on 29 August (Figure S11). The relative abundance of the dominant species (Dolichospermum sp.90 1 August and Microcystis aeruginosa 29 August) decreases after H2O2 oxidation. The relative abundance of Microcystis aeruginosa in the H2O2 trial on 1 August increased after high H2O2 exposure. Microcystis aeruginosa was susceptible to H2O2 oxidation when abundant (as the relative abundance decreased), but it persists as an abundant species. Our results are in accordance with Lusty and Gobler. (2020) [37]. Figure S12 unveils higher susceptibility of Dolichospermum species to high H2O2 exposure as compared to Microcystis aeruginosa. Our results show the effect of high H2O2 exposure on the cyanobacteria species, which is in accordance with previous studies [61,64].
Cyanobacterial composition analysis at species level exhibits relative persistence of different Dolichospermum species and Microcystis aeruginosa when it is abundant following soft-chlorination, soft ozonation and permanganate oxidation. This result is in accordance with the previous studies, which were mainly based on the cell count and the lab-cultured species [10,12,14,58,59]. Our results highlight the ability of hydrogen peroxide to decrease the relative abundance of different taxa, including the toxin-producing taxa of interest. H2O2 selectivity provides an efficient barrier against toxic cyanobacteria entering a drinking water treatment plant.

2.5. Comparison of the Microscopic Cell Count vs. High Throughput Sequencing Results

The taxonomic cell count results (genus) from the chlorination on 1 August and ozonation on 15 August are presented in Figure S13. The observed genera from microscopic cell counts do not completely match with high throughput sequencing results as Aphanocapsa and Aphanothece were only reported by taxonomic cell count. Misclassification at the genus level is less common than for species level; microscopic cell counts showed that Dolichospermum spiroids, Aphanocapsa delicatissma, Aphanotheche clathrate brevis and Aphanocapsa holistca were the most abundant species following chlorination and ozonation. However, high throughput sequencing did not identify them as an abundant species. The abundant species from high throughput sequencing were Microcystis aeruginosa, Dolichospermum Sp.90 and, Dolichospermum cylindrica (Figures S6 and S10). Despite the potential of the microscopic cell count to provide absolute quantitative data such as cells/mL and biovolume, it has some drawbacks.
The differences in community composition structure retrieved from the microscopic taxonomic cell counts and high throughput sequencing could be the result of the limitations inherent to these methods. For high throughput sequencing, these limitations include incomplete DNA sequencing libraries or using different libraries to identify genus and species. In the case of microscopic taxonomic cell counts, several sources of uncertainty and error have been identified [65,66]. The morphological similarity among cyanobacteria taxa may lead to overlooking or misidentifying cyanobacteria under the microscope, especially low abundant species [5,28]. Few genomes might be available for some taxa identified by taxonomic cell count (e.g., Aphanocapsa and Aphanothece), which may result in low relative abundance of these taxa using high throughput sequencing. Additionally, oxidation may further hinder the ability to identify cells because of its impact on the cell structure. Oxidation at higher dosages can cause significant morphological deformation of the cyanobacteria species, especially for H2O2 and KMnO4 in our study [63]. Flow cytometry conducted using the method described by Moradinejad et al. (2019) showed partial membrane damage (up to 80%) under the soft-oxidation conditions performed in this study (data not shown) [63].
This study is the first to provide insights into the impact of different pre-oxidants (doses and contact time) on the diversity and cyanobacterial community composition. The experimental design focused on CTs at which the shifts occurred, in order to provide actionable results for water utilities. Additional investigation would be beneficial to extend our observations to other water bodies. In addition, using assembly/binning method is suggested for future studies to provide a more accurate view of cyanobacterial genus/species.

3. Conclusions

A comparison of the microscopic vs. high throughput sequencing results demonstrates the ability and robustness of high throughput sequencing to fully reveal a wide diversity of cyanobacterial communities in response to oxidant stress.
Results from longitudinal sampling over a bloom period of 4 weeks by high throughput sequencing highlight quick composition and abundance shifts in the cyanobacterial communities that occur within a day. High throughput sequencing revealed clearer shifts during a bloom from an initial dominance of Dolichospermum/Dolichospermum Sp.90 toward a late summer dominance by Microcystis/Microcystis aeruginosa.
Overall, pre-oxidation caused deeper changes in the diversity of whole bacterial communities, especially proteobacteria, than was observed for the cyanobacterial community. Such changes should be considered when assessing the impact of using oxidants for onsite source control.
Depending on the oxidants used, alpha diversity indexes (Shannon and Chao1) showed that oxidation resulted in different structural composition shifts within bacterial and cyanobacterial communities.
Soft-chlorination using dosages 0.2 and 0.6 mg/L that were typically used for pre-oxidation caused a progressive shift in the bacterial and cyanobacterial communities with increasing CTs (0–3.8 mg·min/L). The results using KMnO4 (5 mg/L) and H2O2 (10 mg/L) induced larger and distinct shifts in the structural composition of bacterial and cyanobacterial communities.
Regardless of the significant differences in community distribution shifts caused by different oxidants, some toxin producing species could persist after oxidation whether they were dominant species or not.
Soft-chlorination results revealed that Dolichospermum sp.90 was relatively persistent with increasing CTs whether it was the dominant species or not. Relative persistence of Dolichospermum sp.90 was also observed within KMnO4 oxidation with increasing KMnO4 exposure, regardless of the dominant species.
Soft-ozonation using dosage 0.1–0.3 mg/L results showed the relative persistence of Microcystis aeruginosa (the dominant species) with increasing ozone exposures (0–0.9 mg·min/L).
Only pre-oxidation with H2O2 (10 mg/L) caused a clear decrease in the relative abundance of all taxa and some species including the toxin-producing taxa of interest. As such, H2O2 would provide an effective first barrier against toxin producing cyanobacteria entering the drinking water treatment plant.
Selection of the most effective pre-oxidant for drinking water purposes should be made considering its impact on the cyanobacterial community structure/diversity, prevention of disinfection by-product formation and other drivers such as color and the presence of taste odor compounds.

4. Material and Methods

4.1. Sampling Site Description

The oxidation tests were conducted using natural bloom samples. Cyanobacterial bloom samples were collected and transferred to the laboratory from the Bedford water treatment plant intake (Missisquoi Bay—Lake Champlain) in southern Quebec, Canada (45°02′22.0″ N 73°04′40.5″ W). Sampling was performed during the bloom season from 1 August to 29 August 2018. Several cyanobacterial blooms have been documented for the Lake Champlain in previous years [14,21,42,43].

4.2. Chemicals and Reagents

A free chlorine stock of 2000 mg/L was freshly prepared from sodium hypochlorite (5.25%) on the day of the experiment. Free chlorine residual was measured using an N,N-diethyl-p-phenylenediamine (DPD) colorimetric method based on Standard Methods (SM) 4500-Cl G [67]. Samples were dosed at room temperature (22 °C). A stock of sodium thiosulfate 3000 mg/L was used to quench the chlorine samples at a dose of 1.1 mg/L per 1 mg/L chlorine.
A potassium permanganate (KMnO4) stock solution of 5000 mg/L was prepared by dissolving KMnO4 crystals into the ultrapure water. A DPD colorimetric method ratio of 0.891 KMnO4/Cl2 SM 4500-Cl G was used to determine the KMnO4 residual [67]. An amount of 1.2 mg/L sodium thiosulfate per 1 mg/L KMnO4 was used to quench further oxidation with KMnO4.
A bench-scale ozone generator (details in [14]) was used to prepare an ozone stock solution 50–60 mg/L. SM 4500-O3 was used to determine the stock and ozone residual concentration [14]. An ozone stock solution (50–60 mg/L) was prepared with gaseous ozone using a bench-scale ozone generator (additional details in [14]). Ozone stock concentration and residual ozone in water samples were measured using SM 4500-O3 [14]. An amount of 1.6 mg/L sodium thiosulfate per 1 mg/L ozone was used to quench the ozonation samples.
Stabilized hydrogen peroxide (30%, Sigma Aldrich, MO, USA) was used to prepare the Hydrogen peroxide (H2O2) stock (10 g/L). A colorimetric test kit (Chemetrics K-5510, Midlands, VA, USA) was performed to measure the hydrogen peroxide residual. Moreover, 1.2mg/L sodium thiosulfate per 1 mg/L H2O2 was used to quench further oxidation with H2O2.
The selected oxidant doses and contact times are presented in Table 2. The doses (concentration) and contact time have been selected based on the common pre-oxidation doses used in drinking water treatment intake.
Oxidant exposures, concentration vs. contact time, (CT) were calculated using Equation (1):
C T = 0 t [ O x i d a n t ] d t = C 0 k d e c a y ( e k d e c a y t 1 )
where kdecay (min−1) is the first-order decay rate, t (min) is the exposure time, and C0 (mg/L) is the initial concentration of oxidant at time zero. The selection of oxidant dose was based on the commonly used pre-oxidation doses in the operation of the drinking water treatment plant.

4.3. DNA Extraction, Metagenomic Preparation, Sequencing and Bioinformatic Analysis

Total nucleic acid was extracted from the frozen filters (after filtration, samples were quickly transferred to −80 °C before DNA extraction) using an RNeasy PowerWater Kit (Qiagen Group, Germantwon, MD, USA) with modification. Before the extraction, 200 µL of nuclease-free water and 5 µL of TATAA Universal DNA spike II (TATAA Biocenter AB) were added to the filters to evaluate DNA extraction yields using RT-qPCR. RNeasy PowerWater Isolation kit solution PM1 was used to lyse the cells along with Dithiothreitol (DTT), which prevents disulfide bonds forming residues of proteins. A total volume of 60 µL nuclease-free water provided with the kit was used to elute the total nucleic acid, of which 30 µL of DNA (with minimum 1 ng of DNA) extracts were stored at −20 °C. DNA was subsequently purified with the Zymo Kit (Zymo Research, Irvine, CA, USA) according to the manufacturer’s instructions. Each DNA sample was resuspended in 60 µL of nuclease-free water and quantified with a Qubit v.2.0 fluorometer (Life Technologies, Burlington, ON, Canada). A volume of 30 µL DNA was sent for pyrosequencing (Roche 454 FLX instrumentation with Titanium chemistry) to the Genome Quebec.
An Illumina NovaSeq 6000 platform using S4 flow cells was applied to sequence DNA libraries. A home-made bioinformatic pipeline was used for further analysis of Paired-end raw reads of 150 base pairs (bp) as follows. First, raw reads trimming quality was performed using the SolexaQA v3.1.7.1, default parameters [68]. Further analyses were carried out on the trimmed reads shorter than 75 nt. An in-house script, based on the screening of identical leading 20 bp, was used to remove artificial duplicates. Gene fragments were predicted using FragGeneScan-Plus v3.0 based on the trimmed high-quality reads [69]. Then, predicted fragments of protein were clustered at 90% similarity level using cd-hit v4.8.1 [70]. A diamond engine was used for similarity search on the M5nr database based on a representative of each cluster (https://github.com/MG-RAST/myM5NR). Best hits (minimal e-value of 1 × 10−5) combined with the last common ancestor approach were used to assess the taxonomic affiliation of protein fragments. It should be mentioned that the annotation process uses a read-mapping process of small gene fragments of encoding proteins on a large database of proteins. One gene fragments encoding protein in the database could match with multiple species or strains.

4.4. DOC and Cyanobacteria Cell Count

Pre-rinsed 0.45 µm membrane filters (Supor 45 µm, 47 m, PES PALL, Port Washington, NY, USA) and carbon-free glass vials were used for Dissolved Organic Carbon (DOC) samples. DOC measurements were performed via a 5310 total organic carbon analyzer (Sievers Analytical Instruments, Boulder, CO, USA). An inverted microscope with 20× magnification was used to perform cell count samples preserved with Lugol’s Iodine [71,72].

4.5. Statistical Analysis

All analyses were performed using a custom bioinformatics pipeline implemented in R (v.3.6,2, RStudio, Inc., Boston, MA, USA), phyloseq (V.1.28.0) to visualize the community composition at phylum (all bacteria reads), order, and genus (cyanobacteria reads) [73]. The twenty-five most abundant cyanobacteria species were visualized using pheatmap (v.1.0.12) [74]. Then, the alpha diversity metrics were estimated using phyloseq’s estimate richness function (Shannon and Chao1). Taxonomic data were normalized by the centred log-ratio transformation using easy CODA (v.0.31.1) [75]. The beta-diversity was analyzed using the vegan package (v.2.5-6), where the similarity matrices were calculated based on the Euclidean distance [76]. The homogeneity of variances of normalized data related to each oxidant was analyzed before building the model. A Redundancy Analysis (RDA) constrained ordination to each oxidant applied to the cyanobacteria genus and tested by the permutation test (>95% significance).

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/12/11/728/s1. Figure S1: Principal components analysis (PCA) of the normalized relative abundance of comparative metagenomics reads in 29 August 2018 sample. Data are plotted following the genus-level classification (a) PCA analysis of bacterial community following oxidation using different CT (b) PCA of the cyanobacterial community following oxidation using different CT. Figure S2: Relative abundance of the most abundant genus following the oxidation using Cl2, KMnO4, H2O2 (1 August 2018 abundant: Dolichospermum). Figure S3: Cyanobacterial Species heat map following the oxidation using Cl2, KMnO4, H2O2 (29 August 2018 abundant: Microcystis). Figure S4: Alpha diversity measures of cyanobacterial community following oxidation Cl2, KMnO4, H2O2 (29 August 2018, abundant genus: Microcystis). Figure S5: Total cyanobacteria cell counts following chlorination for 1 August 2018 trial and 29 August 2018 trial. Figure S6: Cyanobacterial species heat map following the chlorination (a) 1 August 2018 trial, (b) 29 August 2018 trial. Figure S7: Total cyanobacteria cell counts following the permanganate oxidation for 1 August 2018 trial and 29 August 2018 trial. Figure S8: Cyanobacterial Species heat map following the oxidation using KMnO4 (a) 1 August 2018 trial, (b) 29 August 2018 trial. Figure S9: Total cyanobacteria cell counts following the O3 oxidation for 15 August 2018 trial and 21 August 2018 trial. Figure S10: Cyanobacterial Species heat map following the oxidation using O3 (a) 15 August 2018 trial, (b) 21 August 2018 trial. Figure S11: Total cyanobacteria cell counts following the H2O2 oxidation for 1 August 2018 trial and 29 August 2018 trial. Figure S12: Cyanobacterial Species heat map following the oxidation using H2O2 (a) 1 August 2018 trial (b) 29 August 2018 trial. Figure S13: Relative abundance of cyanobacteria species (via light Microscopy) following oxidation (a) O3 second trial (15 August 2018) (b) Cl2 first trial (1 August 2018).

Author Contributions

Conceptualization, S.M., H.T., A.Z., S.D. and M.P.; Formal analysis, S.M., H.T., J.F.G.M. and Y.T.; Funding acquisition, J.S., S.D. and M.P.; Methodology, S.M., H.T., J.F.G.M., J.S., Y.T., A.Z., S.D. and M.P.; Software, J.F.G.M.; Supervision, A.Z., S.D. and M.P.; Visualization, S.M., H.T., J.F.G.M., S.D. and M.P.; Writing—original draft, S.M., H.T. and J.F.G.M.; Writing—review and editing, S.M., H.T., J.F.G.M., J.S., Y.T., A.Z., S.D. and M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Genome Canada and Genome Quebec: Algal Blooms, Treatment, Risk Assessment, Prediction and Prevention through Genomics (ATRAPP) Project, Grant number Genome Canada/UM RQ000607 and the APC was funded by Genome Canada and Genome Quebec (ATRAPP project).

Acknowledgments

The authors acknowledge support from Algal Blooms, Treatment, Risk Assessment, Prediction and Prevention through Genomics (ATRAPP), the authors thank the staff at NSERC Industrial Chair on Drinking Water at Polytechnique Montreal. The authors thank the staff at Microbial Evolutionary Genomics (Shapiro lab) at Universite de Montreal.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wells, M.L.; Trainer, V.L.; Smayda, T.J.; Karlson, B.S.; Trick, C.G.; Kudela, R.M.; Ishikawa, A.; Bernard, S.; Wulff, A.; Anderson, D.M.; et al. Harmful algal blooms and climate change: Learning from the past and present to forecast the future. Harmful Algae 2015, 49, 68–93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Paerl, H.W.; Paul, V.J. Climate change: Links to global expansion of harmful cyanobacteria. Water Res. 2012, 46, 1349–1363. [Google Scholar] [CrossRef]
  3. Al-Sammak, M.A.; Hoagland, K.D.; Cassada, D.; Snow, D.D. Co-occurrence of the cyanotoxins BMAA, DABA and anatoxin-a in Nebraska reservoirs, fish, and aquatic plants. Toxins 2014, 6, 488–508. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Westrick, J.A.; Szlag, D.C.; Southwell, B.J.; Sinclair, J. A review of cyanobacteria and cyanotoxins removal/inactivation in drinking water treatment. Anal. Bioanal. Chem. 2010, 397, 1705–1714. [Google Scholar] [CrossRef]
  5. Casero, M.C.; Velázquez, D.; Medina-Cobo, M.; Quesada, A.; Cirés, S. Unmasking the identity of toxigenic cyanobacteria driving a multi-toxin bloom by high-throughput sequencing of cyanotoxins genes and 16S rRNA metabarcoding. Sci. Total Environ. 2019, 665, 367–378. [Google Scholar] [CrossRef] [PubMed]
  6. Zamyadi, A.; Ho, L.; Newcombe, G.; Bustamante, H.; Prévost, M. Fate of toxic cyanobacterial cells and disinfection by-products formation after chlorination. Water Res. 2012, 46, 1524–1535. [Google Scholar] [CrossRef] [PubMed]
  7. Pazouki, P.; Prévost, M.; McQuaid, N.; Barbeau, B.; De Boutray, M.-L.; Zamyadi, A.; Dorner, S. Breakthrough of cyanobacteria in bank filtration. Water Res. 2016, 102, 170–179. [Google Scholar] [CrossRef] [PubMed]
  8. Zamyadi, A.; Dorner, S.; Ndong, M.; Ellis, D.; Bolduc, A.; Bastien, C.; Prévost, M. Low-risk cyanobacterial bloom sources: Cell accumulation within full-scale treatment plants. J. Am. Water Work. Assoc. 2013, 105, E651–E663. [Google Scholar] [CrossRef] [Green Version]
  9. Almuhtaram, H.; Cui, Y.; Zamyadi, A.; Hofmann, R. Cyanotoxins and Cyanobacteria Cell Accumulations in Drinking Water Treatment Plants with a Low Risk of Bloom Formation at the Source. Toxins 2018, 10, 430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Fan, J.; Daly, R.; Hobson, P.; Ho, L.; Brookes, J. Impact of potassium permanganate on cyanobacterial cell integrity and toxin release and degradation. Chemosphere 2013, 92, 529–534. [Google Scholar] [CrossRef] [PubMed]
  11. Fan, J.; Ho, L.; Hobson, P.; Daly, R.; Brookes, J.D. Application of Various Oxidants for Cyanobacteria Control and Cyanotoxin Removal in Wastewater Treatment. J. Environ. Eng. 2014, 140, 04014022. [Google Scholar] [CrossRef]
  12. Coral, L.A.; Zamyadi, A.; Barbeau, B.; Bassetti, F.J.; Lapolli, F.R.; Prevost, M. Oxidation of M. aeruginosa and A. flos-aquae by ozone: Impacts on cell integrity and chlorination by-product formation. Water Res. 2013, 47, 2983–2994. [Google Scholar] [CrossRef] [PubMed]
  13. Ding, J.; Shi, H.; Timmons, T.; Adams, C. Release and Removal of Microcystins from Microcystis during Oxidative-, Physical-, and UV-Based Disinfection. J. Environ. Eng. 2010, 136, 2–11. [Google Scholar] [CrossRef]
  14. Zamyadi, A.; Coral, L.A.; Barbeau, B.; Dorner, S.; Lapolli, F.R.; Prévost, M. Fate of toxic cyanobacterial genera from natural bloom events during ozonation. Water Res. 2015, 73, 204–215. [Google Scholar] [CrossRef]
  15. Zamyadi, A.; Greenstein, K.E.; Glover, C.M.; Adams, C.; Rosenfeldt, E.; Wert, E.C. Impact of Hydrogen Peroxide and Copper Sulfate on the Delayed Release of Microcystin. Water 2020, 12, 1105. [Google Scholar] [CrossRef]
  16. Zamyadi, A.; Ho, L.; Newcombe, G.; Daly, R.I.; Burch, M.; Baker, P.; Prévost, M. Release and Oxidation of Cell-Bound Saxitoxins during Chlorination of Anabaena circinalis Cells. Environ. Sci. Technol. 2010, 44, 9055–9061.e9. [Google Scholar] [CrossRef]
  17. Vlad, S.; Anderson, W.B.; Peldszus, S.; Huck, P.M. Removal of the cyanotoxin anatoxin-a by drinking water treatment processes: A review. J. Water Health 2014, 12, 601–617. [Google Scholar] [CrossRef]
  18. He, X.; Wert, E.C. Colonial cell disaggregation and intracellular microcystin release following chlorination of naturally occurring Microcystis. Water Res. 2016, 101, 10–16. [Google Scholar] [CrossRef]
  19. Wert, E.C.; Dong, M.M.; Rosario-Ortiz, F.L. Using digital flow cytometry to assess the degradation of three cyanobacteria species after oxidation processes. Water Res. 2013, 47, 3752–3761. [Google Scholar] [CrossRef]
  20. Merel, S.; Walker, D.; Chicana, R.; Snyder, S.A.; Baurès, E.; Thomas, O. State of knowledge and concerns on cyanobacterial blooms and cyanotoxins. Environ. Int. 2013, 59, 303–327. [Google Scholar] [CrossRef]
  21. Zamyadi, A.; Dorner, S.; Sauvé, S.; Ellis, D.; Bolduc, A.; Bastien, C.; Prévost, M. Species-dependence of cyanobacteria removal efficiency by different drinking water treatment processes. Water Res. 2013, 47, 2689–2700. [Google Scholar] [CrossRef]
  22. Zhu, L.; Zuo, J.; Song, L.; Gan, N. Microcystin-degrading bacteria affect mcyD expression and microcystin synthesis in Microcystis spp. J. Environ. Sci. 2016, 41, 195–201. [Google Scholar] [CrossRef] [PubMed]
  23. Zamyadi, A.; Romanis, C.; Mills, T.; Neilan, B.; Choo, F.; Coral, L.A.; Gale, D.; Newcombe, G.; Crosbie, N.D.; Stuetz, R.M.; et al. Diagnosing water treatment critical control points for cyanobacterial removal: Exploring benefits of combined microscopy, next-generation sequencing, and cell integrity methods. Water Res. 2019, 152, 96–105. [Google Scholar] [CrossRef] [PubMed]
  24. Hawkins, P.R.; Holliday, J.; Kathuria, A.; Bowling, L. Change in cyanobacterial biovolume due to preservation by Lugol’s Iodine. Harmful Algae 2005, 4, 1033–1043. [Google Scholar] [CrossRef]
  25. American Water Works Association (AWWA). M57-Algae Source to Treatment; American Water Works Association: Denver, CO, USA, 2010. [Google Scholar]
  26. Berry, M.A.; Davis, T.W.; Cory, R.M.; Duhaime, M.B.; Johengen, T.H.; Kling, G.W.; Marino, J.A.; Uyl, P.A.D.; Gossiaux, D.; Dick, G.J.; et al. Cyanobacterial harmful algal blooms are a biological disturbance to Western Lake Erie bacterial communities. Environ. Microbiol. 2017, 19, 1149–1162. [Google Scholar] [CrossRef] [PubMed]
  27. Lezcano, M.Á.; Velázquez, D.; Quesada, A.; El-Shehawy, R. Diversity and temporal shifts of the bacterial community associated with a toxic cyanobacterial bloom: An interplay between microcystin producers and degraders. Water Res. 2017, 125, 52–61. [Google Scholar] [CrossRef] [PubMed]
  28. Kim, K.H.; Yoon, Y.; Hong, W.-Y.; Kim, J.; Cho, Y.-C.; Hwang, S.-J. Application of metagenome analysis to characterize the molecular diversity and saxitoxin-producing potentials of a cyanobacterial community: A case study in the North Han River, Korea. Appl. Biol. Chem. 2018, 61, 153–161. [Google Scholar] [CrossRef] [Green Version]
  29. Eldridge, S.L.C.; Wood, T.M. Annual variations in microcystin occurrence in Upper Klamath Lake, Oregon, based on high-throughput DNA sequencing, qPCR, and environmental parameters. Lake Reserv. Manag. 2019, 36, 31–44. [Google Scholar] [CrossRef]
  30. Casero, M.C.; Ballot, A.; Agha, R.; Quesada, A.; Cirés, S. Characterization of saxitoxin production and release and phylogeny of sxt genes in paralytic shellfish poisoning toxin-producing Aphanizomenon gracile. Harmful Algae 2014, 37, 28–37. [Google Scholar] [CrossRef]
  31. Woodhouse, J.N.; Kinsela, A.S.; Collins, R.N.; Bowling, L.C.; Honeyman, G.L.; Holliday, J.K.; Neilan, B.A. Microbial communities reflect temporal changes in cyanobacterial composition in a shallow ephemeral freshwater lake. ISME J. 2016, 10, 1337–1351. [Google Scholar] [CrossRef]
  32. Pessi, I.S.; Maalouf, P.D.C.; Laughinghouse, H.D.; Baurain, D.; Wilmotte, A. On the use of high-throughput sequencing for the study of cyanobacterial diversity in Antarctic aquatic mats. J. Phycol. 2016, 52, 356–368. [Google Scholar] [CrossRef] [PubMed]
  33. Scherer, P.I.; Millard, A.D.; Miller, A.; Schoen, R.; Raeder, U.; Geist, J.; Zwirglmaier, K. Temporal Dynamics of the Microbial Community Composition with a Focus on Toxic Cyanobacteria and Toxin Presence during Harmful Algal Blooms in Two South German Lakes. Front. Microbiol. 2017, 8, 2387. [Google Scholar] [CrossRef] [PubMed]
  34. Willis, A.; Woodhouse, J.N. Defining Cyanobacterial Species: Diversity and Description Through Genomics. Crit. Rev. Plant Sci. 2020, 39, 101–124. [Google Scholar] [CrossRef]
  35. Xu, H.; Pei, H.; Jin, Y.; Ma, C.; Wang, Y.; Sun, J.; Li, H. High-throughput sequencing reveals microbial communities in drinking water treatment sludge from six geographically distributed plants, including potentially toxic cyanobacteria and pathogens. Sci. Total Environ. 2018, 634, 769–779. [Google Scholar] [CrossRef] [PubMed]
  36. Pei, H.; Xu, H.; Wang, J.; Jin, Y.; Xiao, H.; Ma, C.; Sun, J.; Li, H. 16S rRNA Gene Amplicon Sequencing Reveals Significant Changes in Microbial Compositions during Cyanobacteria-Laden Drinking Water Sludge Storage. Environ. Sci. Technol. 2017, 51, 12774–12783. [Google Scholar] [CrossRef]
  37. Lusty, M.W.; Gobler, C.J. The Efficacy of Hydrogen Peroxide in Mitigating Cyanobacterial Blooms and Altering Microbial Communities across Four Lakes in NY, USA. Toxins 2020, 12, 428. [Google Scholar] [CrossRef]
  38. Newcombe, G.; House, J.; Ho, L.; Baker, P.; Burch, M. Management Strategies for Cyanobacteria (Blue-Green Algae): A Guide for Water Utilities; The Cooperative Research Centre for Water Quality and Treatment: Adelaïde, Australia, 2010; p. 112. [Google Scholar]
  39. Lee, C.; Marion, J.; Cheung, M.; Lee, C.S.; Lee, J. Associations among Human-Associated Fecal Contamination, Microcystis aeruginosa, and Microcystin at Lake Erie Beaches. Int. J. Environ. Res. Public Health 2015, 12, 11466–11485. [Google Scholar] [CrossRef] [Green Version]
  40. Vadde, K.K.; Feng, Q.; Wang, J.; McCarthy, A.J.; Sekar, R. Next-generation sequencing reveals fecal contamination and potentially pathogenic bacteria in a major inflow river of Taihu Lake. Environ. Pollut. 2019, 254, 113108. [Google Scholar] [CrossRef]
  41. Ballesté, E.; Blanch, A.R. Persistence of Bacteroides Species Populations in a River as Measured by Molecular and Culture Techniques. Appl. Environ. Microbiol. 2010, 76, 7608–7616. [Google Scholar] [CrossRef] [Green Version]
  42. Zamyadi, A.; MacLeod, S.L.; Fan, Y.; McQuaid, N.; Dorner, S.; Sauvé, S.; Prévost, M. Toxic cyanobacterial breakthrough and accumulation in a drinking water plant: A monitoring and treatment challenge. Water Res. 2012, 46, 1511–1523. [Google Scholar] [CrossRef]
  43. McQuaid, N.; Zamyadi, A.; Prévost, M.; Bird, D.F.; Dorner, S. Use of in vivophycocyanin fluorescence to monitor potential microcystin-producing cyanobacterial biovolume in a drinkingwater source. J. Environ. Monit. 2011, 13, 455–463. [Google Scholar] [CrossRef] [PubMed]
  44. Tromas, N.; Fortin, N.; Bedrani, L.; Terrat, Y.; Cardoso, P.; Bird, D.; Greer, C.W.; Shapiro, B.J. Characterising and predicting cyanobacterial blooms in an 8-year amplicon sequencing time course. ISME J. 2017, 11, 1746–1763. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Andersson, A.; Höglander, H.; Karlsson, C.; Huseby, S. Key role of phosphorus and nitrogen in regulating cyanobacterial community composition in the northern Baltic Sea. Estuar. Coast. Shelf Sci. 2015, 164, 161–171. [Google Scholar] [CrossRef] [Green Version]
  46. Fortin, N.; Munoz-Ramos, V.; Bird, D.; Lévesque, B.; Whyte, L.G.; Greer, C.W. Toxic Cyanobacterial Bloom Triggers in Missisquoi Bay, Lake Champlain, as Determined by Next-Generation Sequencing and Quantitative PCR. Life 2015, 5, 1346–1380. [Google Scholar] [CrossRef] [Green Version]
  47. Harke, M.J.; Davis, T.W.; Watson, S.B.; Gobler, C.J. Nutrient-Controlled Niche Differentiation of Western Lake Erie Cyanobacterial Populations Revealed via Metatranscriptomic Surveys. Environ. Sci. Technol. 2016, 50, 604–615. [Google Scholar] [CrossRef]
  48. Li, Q.; Lin, F.; Yang, C.; Wang, J.; Lin, Y.; Shen, M.; Park, M.S.; Li, T.; Zhao, J. A Large-Scale Comparative Metagenomic Study Reveals the Functional Interactions in Six Bloom-Forming Microcystis-Epibiont Communities. Front. Microbiol. 2018, 9, 746. [Google Scholar] [CrossRef]
  49. Shi, L.; Cai, Y.; Li, P.; Yang, H.; Liu, Z.; Kong, L.; Yü, Y.; Kong, F. Molecular Identification of the Colony-Associated Cultivable Bacteria of the CyanobacteriumMicrocystis aeruginosaand Their Effects on Algal Growth. J. Freshw. Ecol. 2009, 24, 211–218. [Google Scholar] [CrossRef] [Green Version]
  50. Maruyama, T.; Kato, K.; Yokoyama, A.; Tanaka, T.; Hiraishi, A.; Park, H.-D. Dynamics of microcystin-degrading bacteria in mucilage of Microcystis. Microb. Ecol. 2003, 46, 279–288. [Google Scholar] [CrossRef] [Green Version]
  51. Wacklin, P.; Hoffmann, L.; Komarek, J. Nomenclatural validation of the genetically revised cyanobacterial genus Dolichospermum (RALFS ex BORNET et FLAHAULT) comb. nova. Fottea 2009, 9, 59–64. [Google Scholar] [CrossRef] [Green Version]
  52. Van Wichelen, J.; Vanormelingen, P.; Codd, G.A.; Vyverman, W. The common bloom-forming cyanobacterium Microcystis is prone to a wide array of microbial antagonists. Harmful Algae 2016, 55, 97–111. [Google Scholar] [CrossRef]
  53. Li, Y.; Li, D. Competition between toxic Microcystis aeruginosa and nontoxic Microcystis wesenbergii with Anabaena PCC. Environ. Boil. Fishes 2012, 24, 69–78. [Google Scholar]
  54. Kardinaal, W.; Janse, I.; Agterveld, M.K.-V.; Meima, M.; Snoek, J.; Mur, L.; Huisman, J.; Zwart, G.; Visser, P. Microcystis genotype succession in relation to microcystin concentrations in freshwater lakes. Aquat. Microb. Ecol. 2007, 48, 1–12. [Google Scholar] [CrossRef] [Green Version]
  55. Zhang, X.-W.; Fu, J.; Song, S.; Zhang, P.; Yang, X.-H.; Zhang, L.-R.; Luo, Y.; Liu, C.-H.; Zhu, H.-L. Interspecific competition between Microcystis aeruginosa and Anabaena flos-aquae from Taihu Lake, China. Z. für Nat. C 2014, 69, 53–60. [Google Scholar] [CrossRef]
  56. Chia, M.A.; Jankowiak, J.G.; Kramer, B.J.; Goleski, J.A.; Huang, I.-S.; Zimba, P.V.; Bittencourt-Oliveira, M.D.C.; Gobler, C.J. Succession and toxicity of Microcystis and Anabaena (Dolichospermum) blooms are controlled by nutrient-dependent allelopathic interactions. Harmful Algae 2018, 74, 67–77. [Google Scholar] [CrossRef] [PubMed]
  57. Ndong, M.; Bird, D.; Nguyen-Quang, T.; De Boutray, M.-L.; Zamyadi, A.; Vinçon-Leite, B.; Lemaire, B.J.; Prévost, M.; Dorner, S. Estimating the risk of cyanobacterial occurrence using an index integrating meteorological factors: Application to drinking water production. Water Res. 2014, 56, 98–108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Fan, J.; Hobson, P.; Ho, L.; Daly, R.; Brookes, J.D. The effects of various control and water treatment processes on the membrane integrity and toxin fate of cyanobacteria. J. Hazard. Mater. 2014, 264, 313–322. [Google Scholar] [CrossRef]
  59. Zamyadi, A.; Fan, Y.; Daly, R.I.; Prévost, M. Chlorination of Microcystis aeruginosa: Toxin release and oxidation, cellular chlorine demand and disinfection by-products formation. Water Res. 2013, 47, 1080–1090. [Google Scholar] [CrossRef]
  60. Fan, J.; Ho, L.; Hobson, P.; Brookes, J.D. Evaluating the effectiveness of copper sulphate, chlorine, potassium permanganate, hydrogen peroxide and ozone on cyanobacterial cell integrity. Water Res. 2013, 47, 5153–5164. [Google Scholar] [CrossRef]
  61. Matthijs, H.C.; Visser, P.M.; Reeze, B.; Meeuse, J.; Slot, P.C.; Wijn, G.; Talens, R.; Huisman, J. Selective suppression of harmful cyanobacteria in an entire lake with hydrogen peroxide. Water Res. 2012, 46, 1460–1472. [Google Scholar] [CrossRef]
  62. Zamyadi, A.; Henderson, R.K.; Stuetz, R.; Newcombe, G.; Newtown, K.; Gladman, B. Cyanobacterial management in full-scale water treatment and recycling processes: Reactive dosing following intensive monitoring. Environ. Sci. Water Res. Technol. 2016, 2, 362–375. [Google Scholar] [CrossRef]
  63. Moradinejad, S.; Glover, C.M.; Mailly, J.; Seighalani, T.Z.; Peldszus, S.; Barbeau, B.; Dorner, S.; Prévost, M.; Zamyadi, A. Using Advanced Spectroscopy and Organic Matter Characterization to Evaluate the Impact of Oxidation on Cyanobacteria. Toxins 2019, 11, 278. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Zhou, S.; Shao, Y.; Gao, N.; Deng, Y.; Qiao, J.; Ou, H.; Deng, J. Effects of different algaecides on the photosynthetic capacity, cell integrity and microcystin-LR release of Microcystis aeruginosa. Sci. Total Environ. 2013, 463, 111–119. [Google Scholar] [CrossRef] [PubMed]
  65. Park, J.; Kim, Y.; Kim, M.; Lee, W.H. A novel method for cell counting of Microcystis colonies in water resources using a digital imaging flow cytometer and microscope. Environ. Eng. Res. 2018, 24, 397–403. [Google Scholar] [CrossRef]
  66. Xiao, X.; Sogge, H.; Lagesen, K.; Tooming-Klunderud, A.; Jakobsen, K.S.; Rohrlack, T. Use of High Throughput Sequencing and Light Microscopy Show Contrasting Results in a Study of Phytoplankton Occurrence in a Freshwater Environment. PLoS ONE 2014, 9, e106510. [Google Scholar]
  67. American Public Health Association (APHA); American Water Works Association (AWWA); Water Environment Federation (WEF). Standard Methods for the Examination of Water and Wastewater; American Public Health Association: Washington, DC, USA, 2012; Volume 5, p. 1360. [Google Scholar]
  68. Cox, M.P.; Peterson, D.A.; Biggs, P.J. SolexaQA: At-a-glance quality assessment of Illumina second-generation sequencing data. BMC Bioinform. 2010, 11, 1–6. [Google Scholar] [CrossRef] [Green Version]
  69. Kim, D.; Hahn, A.S.; Wu, S.-J.; Hanson, N.W.; Konwar, K.M.; Hallam, S.J. FragGeneScan-Plus for Scalable High-Throughput Short-Read Open Reading Frame Prediction; IEEE CIBCB Conference: Niagara Falls, ON, Canada, 2015; pp. 1–8. [Google Scholar]
  70. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 2012, 28, 3150–3152. [Google Scholar] [CrossRef]
  71. Lund, J.W.G. A Simple Counting Chamber for Nannoplankton. Limnol. Oceanogr. 1959, 4, 57–65. [Google Scholar] [CrossRef]
  72. Planas, D.; Desrosiers, M.; Groulx, S.; Paquet, S.; Carignan, R. Pelagic and benthic algal responses in eastern Canadian Boreal Shield lakes following harvesting and wildfires. Can. J. Fish. Aquat. Sci. 2000, 57, 136–145. [Google Scholar] [CrossRef]
  73. McMurdie, P.J.; Holmes, S.P. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar]
  74. Raivo, K. Pheatmap: Pretty Heatmaps, R package version 1.0.12; GitHub, Inc.: San Francisco, CA, USA, 2019. [Google Scholar]
  75. Greenacre, M. Compositional Data Analysis in Practice, 1st ed.; CRC Press: Cleveland, OH, USA, 2018. [Google Scholar]
  76. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.R.; O’hara, R.B.; Simpson, G.L.; Solymos, P.; Stevens, M.H.H.; Wagner, H. Community Ecology Package, R Package ‘Vegan’, version 2.5-0; GitHub, Inc.: San Francisco, CA, USA, 2016; p. 285. [Google Scholar]
Figure 1. Identity of major detected bloom-associated cyanobacterial community members during the sampling period: (a) relative abundance of the different phylum, (b) the relative abundance of orders belonging to cyanobacterial phylum, (c) relative abundance of genera belonging to the Nostocales, Chrooccocales and Oscillatoriales orders.
Figure 1. Identity of major detected bloom-associated cyanobacterial community members during the sampling period: (a) relative abundance of the different phylum, (b) the relative abundance of orders belonging to cyanobacterial phylum, (c) relative abundance of genera belonging to the Nostocales, Chrooccocales and Oscillatoriales orders.
Toxins 12 00728 g001
Figure 2. Principal components analysis (PCA) of the normalized relative abundance of comparative metagenomics reads in 1 August sample. Data are plotted following the genus-level classification (a) PCA analysis of bacterial community following oxidation using different common exposure units (CT), (b) PCA of the cyanobacterial community following oxidation using different CT.
Figure 2. Principal components analysis (PCA) of the normalized relative abundance of comparative metagenomics reads in 1 August sample. Data are plotted following the genus-level classification (a) PCA analysis of bacterial community following oxidation using different common exposure units (CT), (b) PCA of the cyanobacterial community following oxidation using different CT.
Toxins 12 00728 g002
Figure 3. Redundancy analysis (RDA) of oxidant effect on cyanobacterial diversity and the cyanobacterial community at genus level Cl2 (0.6 mg/L), KMnO4 (5 mg/L), H2O2 (10 mg/L) (a) 1 August 2018 (b) 29 August 2018.
Figure 3. Redundancy analysis (RDA) of oxidant effect on cyanobacterial diversity and the cyanobacterial community at genus level Cl2 (0.6 mg/L), KMnO4 (5 mg/L), H2O2 (10 mg/L) (a) 1 August 2018 (b) 29 August 2018.
Toxins 12 00728 g003
Figure 4. Cyanobacterial species heat map following the oxidation using Cl2 (0.6 mg/L), KMnO4 (5 mg/L), H2O2 (10 mg/L) (1 August 2018).
Figure 4. Cyanobacterial species heat map following the oxidation using Cl2 (0.6 mg/L), KMnO4 (5 mg/L), H2O2 (10 mg/L) (1 August 2018).
Toxins 12 00728 g004
Figure 5. Alpha diversity measures of the cyanobacterial community following oxidation Cl2 (0.6 mg/L), KMnO4 (5 mg/L), H2O2 (10 mg/L) (1 August 2018).
Figure 5. Alpha diversity measures of the cyanobacterial community following oxidation Cl2 (0.6 mg/L), KMnO4 (5 mg/L), H2O2 (10 mg/L) (1 August 2018).
Toxins 12 00728 g005
Figure 6. The relative abundance of the most abundant genus following chlorination (0.6 mg/L and 0.2 mg/L) (a) 1 August 2018 (Dolichospermum genus abundant) (b) 29 August 2018 (Microcystis genus abundant).
Figure 6. The relative abundance of the most abundant genus following chlorination (0.6 mg/L and 0.2 mg/L) (a) 1 August 2018 (Dolichospermum genus abundant) (b) 29 August 2018 (Microcystis genus abundant).
Toxins 12 00728 g006
Figure 7. Relative abundance of the most abundant genus following KMnO4 (5 mg/L) oxidation (a) 1 August 2018 (Dolichospermum genus abundant) (b) 29 August 2018 (Microcystis genus abundant).
Figure 7. Relative abundance of the most abundant genus following KMnO4 (5 mg/L) oxidation (a) 1 August 2018 (Dolichospermum genus abundant) (b) 29 August 2018 (Microcystis genus abundant).
Toxins 12 00728 g007
Figure 8. Relative abundance of the most abundant genus following O3 (0.3 mg/L and 0.1 mg/L) oxidation (a) 15 August 2018 (Dolichospermum genus Abundant), (b) 21 August 2018 (Synechococcus genus abundant).
Figure 8. Relative abundance of the most abundant genus following O3 (0.3 mg/L and 0.1 mg/L) oxidation (a) 15 August 2018 (Dolichospermum genus Abundant), (b) 21 August 2018 (Synechococcus genus abundant).
Toxins 12 00728 g008
Figure 9. Relative abundance of the most abundant genus following H2O2 (10 mg/L) oxidation on (a) 1 August 2018 (Dolichospermum genus abundant) (b) 29 August 2018 (Microcystis genus abundant).
Figure 9. Relative abundance of the most abundant genus following H2O2 (10 mg/L) oxidation on (a) 1 August 2018 (Dolichospermum genus abundant) (b) 29 August 2018 (Microcystis genus abundant).
Toxins 12 00728 g009
Table 1. Cyanobacterial bloom characteristics.
Table 1. Cyanobacterial bloom characteristics.
Sampling DateDOC (mg/L)pHCell Count (cells/mL)Biovolume (mm3/L)
1 August 20185.97.63.3 × 10530.6
13 August 20185.87.37.8 × 1044.6
15 August 20185.57.41.4 × 1059.4
21 August 20184.97.46.8 × 1040.3
29 August 20185.67.55.4 × 1042.1
Table 2. Experimental plan.
Table 2. Experimental plan.
OxidantWater TypeOxidant Dose (mg/L)Contact Time
Cl2Real Bloom, Missisquoi Bay, Quebec, Canada0.2

0.6
1 min
2 min
5 min
10 min
60 min
120 min
O30.1

0.3
1 min
2 min
5 min
10 min
KMnO4530 min
60 min
120 min
H2O2106 h
24 h
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Moradinejad, S.; Trigui, H.; Guerra Maldonado, J.F.; Shapiro, J.; Terrat, Y.; Zamyadi, A.; Dorner, S.; Prévost, M. Diversity Assessment of Toxic Cyanobacterial Blooms during Oxidation. Toxins 2020, 12, 728. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12110728

AMA Style

Moradinejad S, Trigui H, Guerra Maldonado JF, Shapiro J, Terrat Y, Zamyadi A, Dorner S, Prévost M. Diversity Assessment of Toxic Cyanobacterial Blooms during Oxidation. Toxins. 2020; 12(11):728. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12110728

Chicago/Turabian Style

Moradinejad, Saber, Hana Trigui, Juan Francisco Guerra Maldonado, Jesse Shapiro, Yves Terrat, Arash Zamyadi, Sarah Dorner, and Michèle Prévost. 2020. "Diversity Assessment of Toxic Cyanobacterial Blooms during Oxidation" Toxins 12, no. 11: 728. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12110728

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