Environmental Studies of Cyanobacterial Harmful Algal Blooms Should Include Interactions with the Dynamic Microbiome

Microbial diversity, genomics, and phage-host interactions of cyanobacterial harmful algal blooms

The occurrence of cyanobacterial harmful algal blooms (cyanoHABs) is related to their physical and chemical environment. However, less is known about their associated microbial interactions and processes. In this study, cyanoHABs were analyzed as a microbial ecosystem, using 1 year of 16S rRNA sequencing and 70 metagenomes collected during the bloom season from Lake Okeechobee (Florida, USA). Biogeographical patterns observed in microbial community composition and function reflected ecological zones distinct in their physical and chemical parameters that resulted in bloom "hotspots" near major lake inflows. Changes in relative abundances of taxa within multiple phyla followed increasing bloom severity. Functional pathways that correlated with increasing bloom severity encoded organic nitrogen and phosphorus utilization, storage of nutrients, exchange of genetic material, phage defense, and protection against oxidative stress, suggesting that microbial interactions may promote cyanoHAB resilience. Cyanobacterial communities were highly diverse, with picocyanobacteria ubiquitous and oftentimes most abundant, especially in the absence of blooms. The identification of novel bloom-forming cyanobacteria and genomic comparisons indicated a functionally diverse cyanobacterial community with differences in its capability to store nitrogen using cyanophycin and to defend against phage using CRISPR and restriction-modification systems. Considering blooms in the context of a microbial ecosystem and their interactions in nature, physiologies and interactions supporting the proliferation and stability of cyanoHABs are proposed, including a role for phage infection of picocyanobacteria. This study displayed the power of "-omics" to reveal important biological processes that could support the effective management and prediction of cyanoHABs.

IMPORTANCE

Cyanobacterial harmful algal blooms pose a significant threat to aquatic ecosystems and human health. Although physical and chemical conditions in aquatic systems that facilitate bloom development are well studied, there are fundamental gaps in the biological understanding of the microbial ecosystem that makes a cyanobacterial bloom. High-throughput sequencing was used to determine the drivers of cyanobacteria blooms in nature. Multiple functions and interactions important to consider in cyanobacterial bloom ecology were identified. The microbial biodiversity of blooms revealed microbial functions, genomic characteristics, and interactions between cyanobacterial populations that could be involved in bloom stability and more coherently define cyanobacteria blooms. Our results highlight the importance of considering cyanobacterial blooms as a microbial ecosystem to predict, prevent, and mitigate them.

Characterizing the Impact of Cyanobacterial Blooms on the Photoreactivity of Surface Waters from New York Lakes: A Combined Statewide Survey and Laboratory Investigation

Cyanobacterial blooms introduce autochthonous dissolved organic matter (DOM) into aquatic environments, but their impact on surface water photoreactivity has not been investigated through collaborative field sampling with comparative laboratory assessments. In this work, we quantified the apparent quantum yields (Φapp,RI) of reactive intermediates (RIs), including excited triplet states of dissolved organic matter (3DOM*), singlet oxygen (1O2), and hydroxyl radicals (•OH), for whole water samples collected by citizen volunteers from more than 100 New York lakes. Multiple comparisons tests and orthogonal partial least-squares analysis identified the level of cyanobacterial chlorophyll a as a key factor in explaining the enhanced photoreactivity of whole water samples sourced from bloom-impacted lakes. Laboratory recultivation of bloom samples in bloom-free lake water demonstrated that apparent increases in Φapp,RI during cyanobacterial growth were likely driven by the production of photoreactive moieties through the heterotrophic transformation of freshly produced labile bloom exudates. Cyanobacterial proliferation also altered the energy distribution of 3DOM* and contributed to the accelerated transformation of protriptyline, a model organic micropollutant susceptible to photosensitized reactions, under simulated sunlight conditions. Overall, our study provides insights into the relationship between the photoreactivity of surface waters and the limnological characteristics and trophic state of lakes and highlights the relevance of cyanobacterial abundance in predicting the photoreactivity of bloom-impacted surface waters.

Indole-3-acetic acid promotes growth in bloom-forming Microcystis via an antioxidant response

Interactions between bacteria and phytoplankton in the phycosphere facilitate and constrain biogeochemical cycling in aquatic ecosystems. Indole-3-acetic acid (IAA) is a bacterially produced chemical signal that promotes growth of phytoplankton and plants. Here, we explored the impact of IAA on bloom-forming cyanobacteria and their associated bacteria. Exposure to IAA and its precursor, tryptophan, resulted in a strong growth response in a bloom of the freshwater cyanobacterium, Microcystis. Metatranscriptome analysis revealed the induction of an antioxidant response in Microcystis upon exposure to IAA, potentially allowing populations to increase photosynthetic rate and overcome internally generated reactive oxygen. Our data reveal that co-occurring bacteria within the phycosphere microbiome exhibit a division of labor for supportive functions, such as nutrient mineralization and transport, vitamin synthesis, and reactive oxygen neutralization. These complex dynamics within the Microcystis phycosphere microbiome are an example of interactions within a microenvironment that can have ecosystem-scale consequences.

Stringent Response of Cyanobacteria and Other Bacterioplankton during Different Stages of a Harmful Cyanobacterial Bloom

We conducted a field study to investigate the role of stringent response in cyanobacteria and coexisting bacterioplankton during nutrient-deprived periods at various stages of bloom in a freshwater lake (Utah Lake) for the first time. Using metagenomics and metatranscriptomics analyses, we examined the cyanobacterial ecology and expression of important functional genes related to stringent response, N and P metabolism, and regulation. Our findings mark a significant advancement in understanding the mechanisms by which toxic cyanobacteria survive and proliferate during nitrogen (N) and phosphorus (P) limitations. We successfully identified and analyzed the metagenome-assembled genomes (MAGs) of the dominant bloom-forming cyanobacteria, namely, Dolichospermum circinale, Aphanizomenon flos-aquae UKL13-PB, Planktothrix agardhii, and Microcystis aeruginosa. By mapping RNA-seq data to the coding sequences of the MAGs, we observed that these four prevalent cyanobacteria species activated multiple functions to adapt to the depletion of inorganic nutrients. During and after the blooms, the four dominant cyanobacteria species expressed high levels of transcripts related to toxin production, such as microcystins (mcy), anatoxins (ana), and cylindrospermopsins (cyr). Additionally, genes associated with polyphosphate (poly-P) storage and the stringent response alarmone (p)ppGpp synthesis/hydrolysis, including ppk, relA, and spoT, were highly activated in both cyanobacteria and bacterioplankton. Under N deficiency, the main N pathways shifted from denitrification and dissimilatory nitrate reduction in bacterioplankton toward N2-fixing and assimilatory nitrate reduction in certain cyanobacteria with a corresponding shift in the community composition. P deprivation triggered a stringent response mediated by spoT-dependent (p)ppGpp accumulation and activation of the Pho regulon in both cyanobacteria and bacterioplankton, facilitating inorganic and organic P uptake. The dominant cyanobacterial MAGs exhibited the presence of multiple alkaline phosphatase (APase) transcripts (e.g., phoA in Dolichospermum, phoX in Planktothrix, and Microcystis), suggesting their ability to synthesize and release APase enzymes to convert ambient organic P into bioavailable forms. Conversely, transcripts associated with bacterioplankton-dominated pathways like denitrification were low and did not align with the occurrence of intense cyanoHABs. The strong correlations observed among N, P, stringent response metabolisms and the succession of blooms caused by dominant cyanobacterial species provide evidence that the stringent response, induced by nutrient limitation, may activate unique N and P functions in toxin-producing cyanobacteria, thereby sustaining cyanoHABs.

Spatiotemporal diversity and community structure of cyanobacteria and associated bacteria in the large shallow subtropical Lake Okeechobee (Florida, United States)

Lake Okeechobee is a large eutrophic, shallow, subtropical lake in south Florida, United States. Due to decades of nutrient loading and phosphorus rich sediments, the lake is eutrophic and frequently experiences cyanobacterial harmful algal blooms (cyanoHABs). In the past, surveys of the phytoplankton community structure in the lake have been conducted by morphological studies, whereas molecular based studies have been seldom employed. With increased frequency of cyanoHABs in Lake Okeechobee (e.g., 2016 and 2018 Microcystis-dominated blooms), it is imperative to determine the diversity of cyanobacterial taxa that exist within the lake and the limnological parameters that drive bloom-forming genera. A spatiotemporal study of the lake was conducted over the course of 1 year to characterize the (cyano)bacterial community structure, using 16S rRNA metabarcoding, with coincident collection of limnological parameters (e.g., nutrients, water temperature, major ions), and cyanotoxins. The objectives of this study were to elucidate spatiotemporal trends of community structure, identify drivers of community structure, and examine cyanobacteria-bacterial relationships within the lake. Results indicated that cyanobacterial communities within the lake were significantly different between the wet and dry season, but not between periods of nitrogen limitation and co-nutrient limitation. Throughout the year, the lake was primarily dominated by the picocyanobacterium Cyanobium. The bloom-forming genera Cuspidothrix, Dolichospermum, Microcystis, and Raphidiopsis were highly abundant throughout the lake and had disparate nutrient requirements and niches within the lake. Anatoxin-a, microcystins, and nodularins were detected throughout the lake across both seasons. There were no correlated (cyano)bacteria shared between the common bloom-forming cyanobacteria Dolichospermum, Microcystis, and Raphidiopsis. This study is the first of its kind to use molecular based methods to assess the cyanobacterial community structure within the lake. These data greatly improve our understanding of the cyanobacterial community structure within the lake and the physiochemical parameters which may drive the bloom-forming taxa within Lake Okeechobee.

The phyto-bacterioplankton couple in a shallow freshwater ecosystem: Who leads the dance?

Bloom-forming phytoplankton dynamics are still unpredictable, even though it is known that several abiotic factors, such as nutrient availability and temperature, are key factors for bloom development. We investigated whether biotic factors, i.e. the bacterioplankton composition (via 16SrDNA metabarcoding), were correlated with phytoplankton dynamics, through a weekly monitoring of a shallow lake known to host recurrent cyanobacterial blooms. We detected concomitant changes in both bacterial and phytoplankton community biomass and diversity. During the bloom event, a significant decrease in phytoplankton diversity, was detected, with a first co-dominance of Ceratium, Microcystis and Aphanizomenon, followed by a co-dominance of the two cyanobacterial genera. In the same time, we observed a decrease of the particle-associated (PA) bacterial richness and the emergence of a specific bacterial consortium that was potentially better adapted to the new nutritional niche. Unexpectedly, changes in PA bacterial communities occurred just before the development the emergence of the phytoplanktonic bloom and the associated modification of the phytoplanktonic community composition, suggesting that changes in environmental conditions leading to the bloom, were first sensed by the bacterial PA community. This last was quite stable throughout the bloom event, even though there were changes in the blooming species, suggesting that the association between cyanobacterial species and bacterial communities may not be as tight as previously described for monospecific blooming communities. Finally, the dynamics of the free-living (FL) bacterial communities displayed a different trajectory from those of the PA and phytoplankton communities. This FL communities can be viewed as a reservoir for bacterial recruitment for the PA fraction. Altogether, these data also highlight s that the spatial organization within these different microenvironments in the water column is a relevant factor in the structuring of these communities.

Abundance trade-offs and dominant taxa maintain the stability of the bacterioplankton community underlying Microcystis blooms

Microcystis blooms are an intractable global environmental problem that pollute water and compromise ecosystem functioning. Closed-lake management practices keep lakes free of sewage and harmful algae invasions and have succeeded in controlling local Microcystis blooms; however, there is little understanding of how the bacterioplankton communities associated with Microcystis have changed. Here, based on metagenomic sequencing, the phyla, genera, functional genes and metabolic functions of the bacterioplankton communities were compared between open lakes (underlying Microcystis blooms) and closed lakes (no Microcystis blooms). Water properties and zooplankton density were investigated and measured as factors influencing blooms. The results showed that (1) the water quality of closed lakes was improved, and the nitrogen and phosphorus concentrations were significantly reduced. (2) The stability of open vs. closed-managed lakes differed notably at the species and genus levels (p < 0.01), but no significant variations were identified at the phylum and functional genes levels (p > 0.05). (3) The relative abundance of Microcystis (Cyanobacteria) increased dramatically in the open lakes (proportions from 1.44 to 41.76%), whereas the relative abundance of several other dominant genera of Cyanobacteria experienced a trade-off and decreased with increasing Microcystis relative abundance. (4) The main functions of the bacterioplankton communities were primarily related to dominant genera of Proteobacteria and had no significant relationship with Microcystis. Overall, the closed-lake management practices significantly reduced nutrients and prevented Microcystis blooms, but the taxonomic and functional structures of bacterioplankton communities remained stable overall.

Evaluation of environmental factors and microbial community structure in an important drinking-water reservoir across seasons

The composition of microbial communities varies in water and sediments, and changes in environmental factors have major effects on microbiomes. Here, we characterized variations in microbial communities and physicochemical factors at two sites in a large subtropical drinking water reservoir in southern China. The microbiomes of all sites, including the diversity and abundance of microbial species, were determined via metagenomics, and the relationships between microbiomes and physicochemical factors were determined via redundancy analysis. The dominant species in sediment and water samples differed; Dinobryon sp. LO226KS and Dinobryon divergens were dominant in sediment samples, whereas Candidatus Fonsibacter ubiquis and Microcystis elabens were dominant in water. The diversity was also significantly different in microbial alpha diversity between water and sediment habitats (p < 0.01). The trophic level index (TLI) was the major factor affecting the microbial community in water samples; Mycolicibacterium litorale and Mycolicibacterium phlei were significantly positively related to TLI. Furthermore, we also studied the distribution of algal toxin-encoding genes and antibiotic-resistant genes (ARGs) in the reservoir. It found that water samples contained more phycotoxin genes, with the cylindrospermopsin gene cluster most abundant. We found three genera highly related to cylindrospermopsin and explored a new cyanobacteria Aphanocapsa montana that may produce cylindrospermopsin based on the correlation through network analysis. The multidrug resistance gene was the most abundant ARG, while the relationship between ARGs and bacteria in sediment samples was more complicated than in water. The results of this study enhance our understanding of the effects of environmental factors on microbiomes. In conclusion, research on the properties, including profiles of algal toxin-encoding genes and ARGs, and microbial communities can aid water quality monitoring and conservation.

Spatio-temporal connectivity of the aquatic microbiome associated with cyanobacterial blooms along a Great Lake riverine-lacustrine continuum

Lake Erie is subject to recurring events of cyanobacterial harmful algal blooms (cHABs), but measures of nutrients and total phytoplankton biomass seem to be poor predictors of cHABs when taken individually. A more integrated approach at the watershed scale may improve our understanding of the conditions that lead to bloom formation, such as assessing the physico-chemical and biological factors that influence the lake microbial community, as well as identifying the linkages between Lake Erie and the surrounding watershed. Within the scope of the Government of Canada's Genomics Research and Development Initiative (GRDI) Ecobiomics project, we used high-throughput sequencing of the 16S rRNA gene to characterize the spatio-temporal variability of the aquatic microbiome in the Thames River-Lake St. Clair-Detroit River-Lake Erie aquatic corridor. We found that the aquatic microbiome was structured along the flow path and influenced mainly by higher nutrient concentrations in the Thames River, and higher temperature and pH downstream in Lake St. Clair and Lake Erie. The same dominant bacterial phyla were detected along the water continuum, changing only in relative abundance. At finer taxonomical level, however, there was a clear shift in the cyanobacterial community, with Planktothrix dominating in the Thames River and Microcystis and Synechococcus in Lake St. Clair and Lake Erie. Mantel correlations highlighted the importance of geographic distance in shaping the microbial community structure. The fact that a high proportion of microbial sequences found in the Western Basin of Lake Erie were also identified in the Thames River, indicated a high degree of connectivity and dispersal within the system, where mass effect induced by passive transport play an important role in microbial community assembly. Nevertheless, some cyanobacterial amplicon sequence variants (ASVs) related to Microcystis, representing less than 0.1% of relative abundance in the upstream Thames River, became dominant in Lake St. Clair and Erie, suggesting selection of those ASVs based on the lake conditions. Their extremely low relative abundances in the Thames suggest additional sources are likely to contribute to the rapid development of summer and fall blooms in the Western Basin of Lake Erie. Collectively, these results, which can be applied to other watersheds, improve our understanding of the factors influencing aquatic microbial community assembly and provide new perspectives on how to better understand the occurrence of cHABs in Lake Erie and elsewhere.

Climate change and the aquatic continuum: A cyanobacterial comeback story

Billions of years ago, the Earth's waters were dominated by cyanobacteria. These microbes amassed to such formidable numbers, they ushered in a new era-starting with the Great Oxidation Event-fuelled by oxygenic photosynthesis. Throughout the following eon, cyanobacteria ceded portions of their global aerobic power to new photoautotrophs with the rise of eukaryotes (i.e. algae and higher plants), which co-existed with cyanobacteria in aquatic ecosystems. Yet while cyanobacteria's ecological success story is one of the most notorious within our planet's biogeochemical history, scientists to this day still seek to unlock the secrets of their triumph. Now, the Anthropocene has ushered in a new era fuelled by excessive nutrient inputs and greenhouse gas emissions, which are again reshaping the Earth's biomes. In response, we are experiencing an increase in global cyanobacterial bloom distribution, duration, and frequency, leading to unbalanced, and in many instances degraded, ecosystems. A critical component of the cyanobacterial resurgence is the freshwater-marine continuum: which serves to transport blooms, and the toxins they produce, on the premise that "water flows downhill". Here, we identify drivers contributing to the cyanobacterial comeback and discuss future implications in the context of environmental and human health along the aquatic continuum. This Minireview addresses the overlooked problem of the freshwater to marine continuum and the effects of nutrients and toxic cyanobacterial blooms moving along these waters. Marine and freshwater research have historically been conducted in isolation and independently of one another. Yet, this approach fails to account for the interchangeable transit of nutrients and biology through and between these freshwater and marine systems, a phenomenon that is becoming a major problem around the globe. This Minireview highlights what we know and the challenges that lie ahead.

Comparative metabolomic analysis of exudates of microcystin-producing and microcystin-free Microcystis aeruginosa strains

Cyanobacterial harmful algal blooms (cHABs) dominated by Microcystis aeruginosa threaten the ecological integrity and beneficial uses of lakes globally. In addition to producing hepatotoxic microcystins (MC), M. aeruginosa exudates (MaE) contain various compounds with demonstrated toxicity to aquatic biota. Previously, we found that the ecotoxicity of MaE differed between MC-producing and MC-free strains at exponential (E-phase) and stationary (S-phase) growth phases. However, the components in these exudates and their specific harmful effects were unclear. In this study, we performed untargeted metabolomics based on liquid chromatography-mass spectrometry to reveal the constituents in MaE of a MC-producing and a MC-free strain at both E-phase and S-phase. A total of 409 metabolites were identified and quantified based on their relative abundance. These compounds included lipids, organoheterocyclic compounds, organic acid, benzenoids and organic oxygen compounds. Multivariate analysis revealed that strains and growth phases significantly influenced the metabolite profile. The MC-producing strain had greater total metabolites abundance than the MC-free strain at S-phase, whereas the MC-free strain released higher concentrations of benzenoids, lipids, organic oxygen, organic nitrogen and organoheterocyclic compounds than the MC-producing strain at E-phase. Total metabolites had higher abundance in S-phase than in E- phase in both strains. Analysis of differential metabolites (DMs) and pathways suggest that lipids metabolism and biosynthesis of secondary metabolites were more tightly coupled to growth phases than to strains. Abundance of some toxic lipids and benzenoids DMs were significantly higher in the MC-free strain than the MC-producing one. This study builds on the understanding of MaE chemicals and their biotoxicity, and adds to evidence that non-MC-producing strains of cyanobacteria may also pose a threat to ecosystem health.

Changes in Microbiome Activity and Sporadic Viral Infection Help Explain Observed Variability in Microcosm Studies

The environmental conditions experienced by microbial communities are rarely fully simulated in the laboratory. Researchers use experimental containers ("bottles"), where natural samples can be manipulated and evaluated. However, container-based methods are subject to "bottle effects": changes that occur when enclosing the plankton community that are often times unexplained by standard measures like pigment and nutrient concentrations. We noted variability in a short-term, nutrient amendment experiment during a 2019 Lake Erie, Microcystis spp. bloom. We observed changes in heterotrophic bacteria activity (transcription) on a time-frame consistent with a response to experimental changes in nutrient availability, demonstrating how the often overlooked microbiome of cyanobacterial blooms can be altered. Samples processed at the time of collection (T0) contained abundant transcripts from Bacteroidetes, which reduced in abundance during incubation in all bottles, including controls. Significant biological variability in the expression of Microcystis-infecting phage was observed between replicates, with phosphate-amended treatments showing a 10-fold variation. The expression patterns of Microcystis-infecting phage were significantly correlated with ∼35% of Microcystis-specific functional genes and ∼45% of the cellular-metabolites measured across the entire microbial community, suggesting phage activity not only influenced Microcystis dynamics, but the biochemistry of the microbiome. Our observations demonstrate how natural heterogeneity among replicates can be harnessed to provide further insight on virus and host ecology.

Dynamic Responses of Endosymbiotic Microbial Communities Within Microcystis Colonies in North American Lakes to Altered Nitrogen, Phosphorus, and Temperature Levels

The toxic cyanobacterium, Microcystis, is a pervasive cyanobacterial harmful algal bloom (CHAB) - forming genus that naturally occurs in colonies that harbor diverse microbiomes of heterotrophic bacteria. While the effects of nutrient loading and climatic warming on CHABs are well-known, little is known regarding how these environmental drivers alter the structural and functional potential of the microbial assemblages associated with blooms that, in turn, may impact cyanobacterial growth. Here, we used next-generation sequencing of 16S ribosomal rRNA genes to characterize the dynamics of the bacterial assemblages within Microcystis colonies in two temperate North American lakes: Lake Erie and Lake Agawam (NY, United States) and quantified their responses to experimentally increased levels of nitrogen (N), phosphorus (P) and temperature. Across experiments, Microcystis populations were consistently and significantly promoted by N and, to a lesser extent, elevated temperature (p < 0.05). In contrast, bacterial assemblages within Microcystis colonies were more resilient to environmental perturbations, with the relative abundance of 7–16% of amplicon sequence variants changing and several individual taxa displaying significant (p < 0.05) increases and decreases in relative abundance, primarily in response to elevated temperature and to a lesser extent, N. In contrast to individual taxa, community diversity was not significantly altered by individual treatments during experiments but rather was inversely correlated with the intensity of Microcystis blooms (p < 0.001). While predicted metabolic function was even less impacted by environmental drivers than microbial diversity, the predicted abundance of nitrogenase (nifH), alkaline phosphatase (phoX), and urease (ure) genes significantly increased in response to N but decreased in response to increased temperature (p < 0.05). Collectively, the resilience of microbial community structure and function within colonies suggests they may support the ability of Microcystis to persist through short-term fluctuations in environmental conditions by supplying essential nutrients.

A comparative study of metatranscriptomic assessment methods to characterize Microcystis blooms

Harmful algal blooms are increasing in duration and severity globally, resulting in increased research interest. The use of genetic sequencing technologies has provided a wealth of opportunity to advance knowledge, but also poses a risk to that knowledge if handled incorrectly. The vast numbers of sequence processing tools and protocols provide a method to test nearly every hypothesis, but each method has inherent strengths and weaknesses. Here, we tested six methods to classify and quantify metatranscriptomic activity from a harmful algal bloom dominated by Microcystis spp. Three online tools were evaluated (Kaiju, MG-RAST, and GhostKOALA) in addition to three local tools that included a command line BLASTx approach, recruitment of reads to individual Microcystis genomes, and recruitment to a combined Microcystis composite genome generated from sequenced isolates with complete, closed genomes. Based on the analysis of each tool presented in this study, two recommendations are made that are dependent on the hypothesis to be tested. For researchers only interested in the function and physiology of Microcystis spp., read recruitments to the composite genome, referred to as "Frankenstein's Microcystis", provided the highest total estimates of transcript expression. However, for researchers interested in the entire bloom microbiome, the online GhostKOALA annotation tool, followed by subsequent read recruitments, provided functional and taxonomic characterization, in addition to transcript expression estimates. This study highlights the critical need for careful evaluation of methods before data analysis.