Structural Diversity of Bacterial Communities Associated with Bloom-Forming Freshwater Cyanobacteria Differs According to the Cyanobacterial Genus

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.

Spatio-temporal connectivity of a toxic cyanobacterial community and its associated microbiome along a freshwater-marine continuum

Due to climate changes and eutrophication, blooms of predominantly toxic freshwater cyanobacteria are intensifying and are likely to colonize estuaries, thus impacting benthic organisms and shellfish farming representing a major ecological, health and economic risk. In the natural environment, Microcystis form large mucilaginous colonies that influence the development of both cyanobacterial and embedded bacterial communities. However, little is known about the fate of natural colonies of Microcystis by salinity increase. In this study, we monitored the fate of a Microcystis dominated bloom and its microbiome along a French freshwater-marine gradient at different phases of a bloom. We demonstrated changes in the cyanobacterial genotypic composition, in the production of specific metabolites (toxins and compatible solutes) and in the heterotrophic bacteria structure in response to the salinity increase. In particular M. aeruginosa and M. wesenbergii survived salinities up to 20. Based on microcystin gene abundance, the cyanobacteria became more toxic during their estuarine transfer but with no selection of specific microcystin variants. An increase in compatible solutes occurred along the continuum with extensive trehalose and betaine accumulations. Salinity structured most the heterotrophic bacteria community, with an increased in the richness and diversity along the continuum. A core microbiome in the mucilage-associated attached fraction was highly abundant suggesting a strong interaction between Microcystis and its microbiome and a likely protecting role of the mucilage against an osmotic shock. These results underline the need to better determine the interactions between the Microcystis colonies and their microbiome as a likely key to their widespread success and adaptation to various environmental conditions.

Cyanosphere Dynamic During Dolichospermum Bloom: Potential Roles in Cyanobacterial Proliferation

Under the effect of global change, management of cyanobacterial proliferation becomes increasingly pressing. Given the importance of interactions within microbial communities in aquatic ecosystems, a handful of studies explored the potential relations between cyanobacteria and their associated bacterial community (i.e., cyanosphere). Yet, most of them specifically focused on the ubiquitous cyanobacteria Microcystis, overlooking other genera. Here, based on 16s rDNA metabarcoding analysis, we confirmed the presence of cyanosphere representing up to 30% of the total bacterial community diversity, during bloom episode of another preponderant cyanobacterial genus, Dolichospermum. Moreover, we highlighted a temporal dynamic of this cyanosphere. A sPLS-DA model permits to discriminate three important dates and 220 OTUs. With their affiliations, we were able to show how these variations potentially imply a turnover in ecological functions depending on bloom phases. Although more studies are necessary to quantify the impacts of these variations, we argue that cyanosphere can have an important, yet underestimated, role in the modulation of cyanobacterial blooms.

The Raphidiopsis (= Cylindrospermopsis) raciborskii pangenome updated: Two new metagenome-assembled genomes from the South American clade

Two Raphidiopsis (=Cylindrospermopsis) raciborskii metagenome-assembled genomes (MAGs) were recovered from two freshwater metagenomic datasets sampled in 2011 and 2012 in Pampulha Lake, a hypereutrophic, artificial, shallow reservoir, located in the city of Belo Horizonte (MG), Brazil. Since the late 1970s, the lake has undergone increasing eutrophication pressure, due to wastewater input, leading to the occurrence of frequent cyanobacterial blooms. The major difference observed between PAMP2011 and PAMP2012 MAGs was the lack of the saxitoxin gene cluster in PAMP2012, which also presented a smaller genome, while PAMP2011 presented the complete sxt cluster and all essential proteins and clusters. The pangenome analysis was performed with all Raphidiopsis/Cylindrospermopsis genomes available at NCBI to date, with the addition of PAMP2011 and PAMP2012 MAGs (All33 subset), but also without the South American strains (noSA subset), and only among the South American strains (SA10 and SA8 subsets). We observed a substantial increase in the core genome size for the 'noSA' subset, in comparison to 'All33' subset, and since the core genome reflects the closeness among the pangenome members, the results strongly suggest that the conservation level of the essential gene repertoire seems to be affected by the geographic origin of the strains being analyzed, supporting the existence of a distinct SA clade. The Raphidiopsis pangenome comprised a total of 7943 orthologous protein clusters, and the two new MAGs increased the pangenome size by 11%. The pangenome based phylogenetic relationships among the 33 analyzed genomes showed that the SA genomes clustered together with 99% bootstrap support, reinforcing the metabolic particularity of the Raphidiopsis South American clade, related to its saxitoxin producing unique ability, while also indicating a different evolutionary history due to its geographic isolation.

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.

Algae drive convergent bacterial community assembly at low dilution frequency

Microbial community assembly is a complex dynamical process that determines community structure and function. The interdependence of inter-species interactions and nutrient availability presents a challenge for understanding community assembly. We sought to understand how external nutrient supply rate modulated interactions to affect the assembly process. A statistical decomposition of taxonomic structures of bacterial communities assembled with and without algae and at varying dilution frequencies allowed the separation of the effects of biotic (presence of algae) and abiotic (dilution frequency) factors on community assembly. For infrequent dilutions, the algae strongly impact community assembly, driving initially diverse bacterial consortia to converge to a common structure. Analyzing sequencing data revealed that this convergence is largely mediated by a decline in the relative abundance of specific taxa in the presence of algae. This study shows that complex phototroph-heterotroph communities can be powerful model systems for understanding assembly processes relevant to the global ecosystem functioning.

Combined impacts of algae-induced variations in water soluble organic matter and heavy metals on bacterial community structure in sediment from Chaohu Lake, a eutrophic shallow lake

Many lakes are suffering from eutrophication and heavy metals-contamination. However, the combined impacts of algae bloom and its induced variations in heavy metals on microbial community in sediment from eutrophic lakes remain unclear. In this study, we performed field experiments to investigate how algae bloom impacted water soluble organic matter (WSOM) and heavy metals in sediment from Chaohu Lake, a eutrophic shallow lake, and probed their combined impacts on sediment bacterial community structure. The results showed that algae bloom increased WSOM quantity, in particular, the soluble microbial by-product-like (SMP) and fulvic acid-like (Fa-L) components markedly enhanced by 203.70 % and 70.17 %, respectively. We also found that algae bloom redistributed the spatial patterns of heavy metals and altered their chemical species in sediment, then promoted contamination degree and potential ecological risk of heavy metals in sediment. Moreover, sediment bacterial community richness and diversity obviously decreased after algae bloom, and the variance partitioning analysis (VPA) results showed that combined impacts of algae-induced changes in WSOM and heavy metals explained 66.56 % of the variations in bacterial community structure. These findings depicted how algae bloom influence sediment WSOM and heavy metals, and revealed the combined impacts of algae-induced variations on microbial community structure in shallow eutrophic lake.

Widespread formation of intracellular calcium carbonates by the bloom-forming cyanobacterium Microcystis

The formation of intracellular amorphous calcium carbonates (iACC) has been recently observed in a few cultured strains of Microcystis, a potentially toxic bloom-forming cyanobacterium found worldwide in freshwater ecosystems. If iACC-forming Microcystis are abundant within blooms, they may represent a significant amount of particulate Ca. Here, we investigate the significance of iACC biomineralization by Microcystis. First, the presence of iACC-forming Microcystis cells has been detected in several eutrophic lakes, indicating that this phenomenon occurs under environmental conditions. Second, some genotypic (presence/absence of ccyA, a marker gene of iACC biomineralization) and phenotypic (presence/absence of iACC) diversity have been detected within a collection of strains isolated from one single lake. This illustrates that this trait is frequent but also variable within Microcystis even at a single locality. Finally, one-third of publicly available genomes of Microcystis were shown to contain the ccyA gene, revealing a wide geographic and phylogenetic distribution within the genus. Overall, the present work shows that the formation of iACC by Microcystis is common under environmental conditions. While its biological function remains undetermined, this process should be further considered regarding the biology of Microcystis and implications on the Ca geochemical cycle in freshwater environments.

Coexistence of Synechococcus and Microcystis Blooms in a Tropical Urban Reservoir and Their Links with Microbiomes

Bacteria play a crucial role in driving ecological processes in aquatic ecosystems. Studies have shown that bacteria-cyanobacteria interactions contributed significantly to phytoplankton dynamics. However, information on the contribution of bacterial communities to blooms remains scarce. Here, we tracked changes in the bacterial community during the development of a cyanobacterial bloom in an equatorial estuarine reservoir. Two forms of blooms were observed simultaneously corresponding to the lotic and lentic characteristics of the sampling sites where significant spatial variabilities in physicochemical water quality, cyanobacterial biomass, secondary metabolites, and cyanobacterial/bacterial compositions were detected. Microcystis dominated the upstream sites during peak periods and were succeeded by Synechococcus when the bloom subsided. For the main body of the reservoir, a mixed bloom featuring coccoid and filamentous cyanobacteria (Microcystis, Synechococcus, Planktothricoides, Nodosilinea, Raphidiopsis, and Prochlorothrix) was observed. Concentrations of the picocyanobacteria Synechococcus remained high throughout the study, and their positive correlations with cylindrospermopsin and anatoxin-a suggested that they could produce cyanotoxins, which pose more damaging impacts than previously supposed. Succession of different cyanobacteria (Synechococcus and Microcystis) following changes in nutrient composition and ionic strength was demonstrated. The microbiomes associated with blooms were unique to the dominant cyanobacteria. Generic and specialized bloom biomarkers for the Microcystis and downstream mixed blooms were also identified. Microscillaceae, Chthoniobacteraceae, and Roseomonas were the major heterotrophic bacteria associated with Microcystis bloom, whereas Phycisphaeraceae and Methylacidiphilaceae were the most prominent groups for the Synechococcus bloom. Collectively, bacterial community can be greatly deviated by the geological condition, monsoon season, cyanobacterial density, and dominant cyanobacteria.

Multi-omics analyses from a single sample: prior metabolite extraction does not alter the 16S rRNA-based characterization of prokaryotic community in a diversity of sample types

Massive sequencing of the 16S rRNA gene has become a standard first step to describe and compare microbial communities from various samples. Parallel analysis of high numbers of samples makes it relevant to the statistical testing of the influence of natural or experimental factors and variables. However, these descriptions fail to document changes in community or ecosystem functioning. Nontargeted metabolomics are a suitable tool to bridge this gap, yet extraction protocols are different. In this study, prokaryotic community compositions are documented by 16S rRNA gene sequencing after direct DNA extraction or after metabolites extraction followed by DNA extraction. Results obtained using the V3-V4 region on nonaxenic cultures of cyanobacteria, lake water column, biofilm, and gut of wild and lab-reared fish indicate that prior extraction of metabolites does not influence the obtained image of prokaryotic communities. This validates sequential extraction of metabolites followed by DNA as a way to combine 16S rRNA sequencing with metabolome characterization from a single sample. This approach has the potential to complement community structure characterization with a proxy of their functioning, without the uncertainties associated with the use of separate samples.

Metagenome-Based Exploration of Bacterial Communities Associated with Cyanobacteria Strains Isolated from Thermal Muds

Cyanobacteria constitute a pioneer colonizer of specific environments for whom settlement in new biotopes precedes the establishment of composite microbial consortia. Some heterotrophic bacteria constitute cyanobacterial partners that are considered as their cyanosphere, being potentially involved in mutualistic relationships through the exchange and recycling of key nutrients and the sharing of common goods. Several non-axenic cyanobacterial strains have been recently isolated, along with their associated cyanospheres, from the thermal mud of Balaruc-les-Bains (France) and the biofilms of the retention basin where they develop. The community structure and relationships among the members of the isolated cyanobacterial strains were characterized using a metagenomic approach combined with taxonomic and microscopic descriptions of the microbial consortia. The results provided insights into the potential role and metabolic capabilities of the microorganisms of thermal mud-associated cyanobacterial biofilms. Thus, the physical proximity, host-specificity, and genetic potential functions advocate for their complementarity between cyanobacteria and their associated microbiota. Besides these findings, our results also highlighted the great influence of the reference protein database chosen for performing functional annotation of the metagenomes from organisms of the cyanosphere and the difficulty of selecting one unique database that appropriately covers both autotroph and heterotroph metabolic specificities.

Temporal and functional interrelationships between bacterioplankton communities and the development of a toxigenic Microcystis bloom in a lowland European reservoir

The cyanobacteria-associated microbiome is constantly reshaped by bloom development. However, the synergistic-antagonistic nature of the relationships between Microcystis and its microbiome still remains unclear. Therefore, temporal changes of bacterioplankton communities and their functional potential through different developing stages of a Microcystis toxigenic bloom were investigated, considering bacterioplankton assemblages as particle-attached (PAB) and free-living (FLB) bacteria. 16S rRNA sequencing revealed that PAB were represented by Proteobacteria and Cyanobacteria, while FLB by Proteobacteria and Actinobacteria. Network and ordination analyses indicated that PAB inter-relationships were more complex-numerous connections between taxa with stronger correlations, than FLB-rather influenced by physico-chemical parameters. PAB in pre-summer was diverse with Proteobacteria containing potential taxa involved in nitrogen-transforming processes. In mid-summer, PAB presented a mix-bloom dominated by Snowella, Aphanizomenon, and Microcystis, which were succeeded by toxigenic Microcystis in post-summer. Both periods were associated to potential taxa with parasitic/predatory lifestyles against cyanobacteria. In post-summer, Sutterellaceae were recognized as poor water quality indicators, and their strong association with Microcystis could have represented an increased threat for that period. Microcystis was a major factor significantly reducing PAB diversity and evenness, suggesting that it negatively influenced bacterioplankton assemblages, probably also altering the overall community functional potential.

The role of organic nutrients in structuring freshwater phytoplankton communities in a rapidly changing world

Carbon, nitrogen, and phosphorus are critical macroelements in freshwater systems. Historically, researchers and managers have focused on inorganic forms, based on the premise that the organic pool was not available for direct uptake by phytoplankton. We now know that phytoplankton can tap the organic nutrient pool through a number of mechanisms including direct uptake, enzymatic hydrolysis, mixotrophy, and through symbiotic relationships with microbial communities. In this review, we explore these mechanisms considering current and projected future anthropogenically-driven changes to freshwater systems. In particular, we focus on how naturally- and anthropogenically- derived organic nutrients can influence phytoplankton community structure. We also synthesize knowledge gaps regarding phytoplankton physiology and the potential challenges of nutrient management in an organically dynamic and anthropogenically modified world. Our review provides a basis for exploring these topics and suggests several avenues for future work on the relation between organic nutrients and eutrophication and their ecological implications in freshwater systems.

Temporal heterogeneity of bacterial communities and their responses to Raphidiopsis raciborskii blooms

To elucidate the interspecies connectivity between cyanobacteria and other bacteria (noncyanobacteria), microbial diversity and composition were investigated through high-throughput sequencing (HTS) in a drinking water reservoir in Chongqing city, Southwest China, during Raphidiopsis raciborskii blooms. Significant temporal changes were observed in microbial community composition during the sampling period, primarily reflected by variations in relative bacterial abundance. The modularity analysis of the network demonstrated that the bacterial community forms co-occurrence/exclusion patterns in response to variations in environmental factors. Moreover, five modules involved in the dynamic phases of the R. raciborskii bloom were categorized into the Pre-Bloom, Bloom, Post-Bloom, and Non-Bloom Groups. The reservoir was eutrophic (i.e., the average concentrations of total nitrogen (TN) and total phosphorus (TP) were 2.32 and 0.07 mg L^-1, respectively) during the investigation; however, Pearson's correlation coefficient showed that R. raciborskii was not significantly correlated with nitrogen and phosphorus. However, other environmental factors, such as water temperature, pH, and the permanganate index, were positively correlated with R. raciborskii. Importantly, Proteobacteria (α-, γ-Proteobacteria), Acidobacteria, Chloroflexi, and Firmicutes were preferentially associated with increased R. raciborskii blooms. These results suggested that the transition of R. raciborskii bloom-related microbial modules and their keystone species could be crucial in the development and collapse of R. raciborskii blooms and could provide a fundamental basis for understanding the linkage between the structure and function of the microbial community during bloom dynamics.

Metagenomics Analysis to Investigate the Microbial Communities and Their Functional Profile During Cyanobacterial Blooms in Lake Varese

Toxic cyanobacterial blooms represent a natural phenomenon caused by a mass proliferation of photosynthetic prokaryotic microorganisms in water environments. Bloom events have been increasingly reported worldwide and their occurrence can pose serious threats to aquatic organisms and human health. In this study, we assessed the microbial composition, with a focus on Cyanobacteria, in Lake Varese, a eutrophic lake located in northern Italy. Water samples were collected and used for obtaining a 16S-based taxonomic profile and performing a shotgun sequencing analysis. The phyla found to exhibit the greatest relative abundance in the lake included Proteobacteria, Cyanobacteria, Actinobacteriota and Bacteroidota. In the epilimnion and at 2.5 × Secchi depth, Cyanobacteria were found to be more abundant compared to the low levels detected at greater depths. The blooms appear to be dominated mainly by the species Lyngbya robusta, and a specific functional profile was identified, suggesting that distinct metabolic processes characterized the bacterial population along the water column. Finally, analysis of the shotgun data also indicated the presence of a large and diverse phage population.

An ensemble approach to the structure-function problem in microbial communities

The metabolic activity of microbial communities plays a primary role in the flow of essential nutrients throughout the biosphere. Molecular genetics has revealed the metabolic pathways that model organisms utilize to generate energy and biomass, but we understand little about how the metabolism of diverse, natural communities emerges from the collective action of its constituents. We propose that quantifying and mapping metabolic fluxes to sequencing measurements of genomic, taxonomic, or transcriptional variation across an ensemble of diverse communities, either in the laboratory or in the wild, can reveal low-dimensional descriptions of community structure that can explain or predict their emergent metabolic activity. We survey the types of communities for which this approach might be best suited, review the analytical techniques available for quantifying metabolite fluxes in communities, and discuss what types of data analysis approaches might be lucrative for learning the structure-function mapping in communities from these data.

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.

Metagenomic analysis of Raphidiopsisraciborskii microbiome: beyond the individual

Raphidiopsisraciborskii is a toxic, invasive bacteria with a defined biogeographic pattern attributed to the generation of ecotypes subjected to local environmental filters and to phenotypic plasticity. The interactions taking place between the cyanobacterium and the other bacteria inhabiting the external polysaccharide-rich matrix surrounding the cells, or phycosphere, may be ecotype-specific and would have different influence on the carbon and nutrient cycling in the ecosystem. Here, we describe the bacterial community or microbiome (assessed by 16S rRNA metagenomics) associated to two R.raciborskii strains that have been described as different ecotypes: the saxitoxin-producer MVCC19 and the non-toxic LB2897. Our results showed that both ecotypes share 50% of their microbiomes and differ in their dominant taxa. The taxon having the highest abundance in the microbiome of MVCC19 was Neorhizobium (22.5% relative abundance), while the dominant taxon in LB2897 was the PlanctomycetesSM1A02 (26.2% relative abundance). These groups exhibit different metabolic capabilities regarding nitrogen acquisition (symbiotic nitrogen-fixing in Neorhizobium vs. anammox in SM1A02), suggesting the existence of ecotype-specific microbiomes that play a relevant role in cyanobacterial niche-adaptation. In addition, as saxitoxin and analogues are nitrogen-rich (7 atoms per molecule), we hypothesise that saxitoxin-producing R.raciborskii benefits from external sources of nitrogen provided by the microbiome bacteria. Based on these findings, we propose that the mechanisms involved in the assembly of the cyanobacterial microbiome community are ecotype-dependent.

The Future Is Big-and Small: Remote Sensing Enables Cross-Scale Comparisons of Microbiome Dynamics and Ecological Consequences

Coupling remote sensing with microbial omics-based approaches provides a promising new frontier for scientists to scale microbial interactions across space and time. These data-rich, interdisciplinary methods allow us to better understand interactions between microbial communities and their environments and, in turn, their impact on ecosystem structure and function. Here, we highlight current and novel examples of applying remote sensing, machine learning, spatial statistics, and omics data approaches to marine, aquatic, and terrestrial systems. We emphasize the importance of integrating biochemical and spatiotemporal environmental data to move toward a predictive framework of microbiome interactions and their ecosystem-level effects. Finally, we emphasize lessons learned from our collaborative research with recommendations to foster productive and interdisciplinary teamwork.

Spatial distribution of sediment archaeal and bacterial communities relates to the source of organic matter and hypoxia - a biogeographical study on Lake Remoray (France)

Bottom waters hypoxia spreads in many lakes worldwide causing severe consequences on whole lakes trophic network. Here, we aimed at understanding the origin of organic matter stored in the sediment compartment and the related diversity of sediment microbial communities in a lake with deoxygenated deep water layers. We used a geostatistical approach to map and compare both the variation of organic matter and microbial communities in sediment. Spatialisation of C/N ratio and δ13C signature of sediment organic matter suggested that Lake Remoray was characterized by an algal overproduction which could be related to an excess of nutrient due to the close lake-watershed connectivity. Three spatial patterns were observed for sediment microbial communities after the hypoxic event, each characterized by specific genetic structure, microbial diversity and composition. The relative abundance variation of dominant microbial groups across Lake Remoray such as Cyanobacteria, Gammaproteobacteria, Deltaproteobacteria and Chloroflexi provided us important information on the lake areas where hypoxia occurs. The presence of methanogenic species in the deeper part of the lake suggests important methane production during hypoxia period. Taken together, our results provide an extensive picture of microbial communities' distribution related to quantity and quality of organic matter in a seasonally hypoxic lake.

Global co-occurrence of methanogenic archaea and methanotrophic bacteria in Microcystis aggregates

Global warming and eutrophication contribute to the worldwide increase in cyanobacterial blooms, and the level of cyanobacterial biomass is strongly associated with rises in methane emissions from surface lake waters. Hence, methane-metabolizing microorganisms may be important for modulating carbon flow in cyanobacterial blooms. Here, we surveyed methanogenic and methanotrophic communities associated with floating Microcystis aggregates in 10 lakes spanning four continents, through sequencing of 16S rRNA and functional marker genes. Methanogenic archaea (mainly Methanoregula and Methanosaeta) were detectable in 5 of the 10 lakes and constituted the majority (~50%-90%) of the archaeal community in these lakes. Three of the 10 lakes contained relatively more abundant methanotrophs than the other seven lakes, with the methanotrophic genera Methyloparacoccus, Crenothrix, and an uncultured species related to Methylobacter dominating and nearly exclusively found in each of those three lakes. These three are among the five lakes in which methanogens were observed. Operational taxonomic unit (OTU) richness and abundance of methanotrophs were strongly positively correlated with those of methanogens, suggesting that their activities may be coupled. These Microcystis-aggregate-associated methanotrophs may be responsible for a hitherto overlooked sink for methane in surface freshwaters, and their co-occurrence with methanogens sheds light on the methane cycle in cyanobacterial aggregates.

Reviewing Interspecies Interactions as a Driving Force Affecting the Community Structure in Lakes via Cyanotoxins

The escalating occurrence of toxic cyanobacterial blooms worldwide is a matter of concern. Global warming and eutrophication play a major role in the regularity of cyanobacterial blooms, which has noticeably shifted towards the predomination of toxic populations. Therefore, understanding the effects of cyanobacterial toxins in aquatic ecosystems and their advantages to the producers are of growing interest. In this paper, the current literature is critically reviewed to provide further insights into the ecological contribution of cyanotoxins in the variation of the lake community diversity and structure through interspecies interplay. The most commonly detected and studied cyanobacterial toxins, namely the microcystins, anatoxins, saxitoxins, cylindrospermopsins and β-N-methylamino-L-alanine, and their ecotoxicity on various trophic levels are discussed. This work addresses the environmental characterization of pure toxins, toxin-containing crude extracts and filtrates of single and mixed cultures in interspecies interactions by inducing different physiological and metabolic responses. More data on these interactions under natural conditions and laboratory-based studies using direct co-cultivation approaches will provide more substantial information on the consequences of cyanotoxins in the natural ecosystem. This review is beneficial for understanding cyanotoxin-mediated interspecies interactions, developing bloom mitigation technologies and robustly assessing the hazards posed by toxin-producing cyanobacteria to humans and other organisms.

Insights into the cyanosphere: capturing the respective metabolisms of cyanobacteria and chemotrophic bacteria in natural conditions?

Specific interactions have been highlighted between cyanobacteria and chemotrophic bacteria within the cyanosphere, suggesting that nutrients recycling could be optimized by cyanobacteria/bacteria exchanges. In order to determine the respective metabolic roles of the cyanobacterial and bacterial consortia (microbiome), a day-night metatranscriptomic analysis was performed on Dolichospermum sp. (N₂ -fixer) and Microcystis sp. (non N₂ -fixer) natural blooms occurring successively within a French peri-urban lake. The taxonomical and functional analysis of the metatranscriptoms have highlighted specific association of bacteria within the cyanosphere, driven by the cyanobacteria identity, without strongly modifying the functional composition of the microbiomes, suggesting functional redundancy within the cyanosphere. Moreover, the functional composition of these active communities was driven by the living mode. During the two successive bloom events, it appeared that NH₄ ⁺ (newly fixed and/or allochthonous) was preferentially transformed into amino acids for the both the microbiome and the cyanobacteria, while phosphate metabolism was enhanced, suggesting that due to a high cellular growth, P limitation might take place within the cyanosphere consortium.

Individual Microcystis colonies harbour distinct bacterial communities that differ by Microcystis oligotype and with time

Interactions between bacteria and phytoplankton in the phycosphere have impacts at the scale of whole ecosystems, including the development of harmful algal blooms. The cyanobacterium Microcystis causes toxic blooms that threaten freshwater ecosystems and human health globally. Microcystis grows in colonies that harbour dense assemblages of other bacteria, yet the taxonomic composition of these phycosphere communities and the nature of their interactions with Microcystis are not well characterized. To identify the taxa and compositional variance within Microcystis phycosphere communities, we performed 16S rRNA V4 region amplicon sequencing on individual Microcystis colonies collected biweekly via high-throughput droplet encapsulation during a western Lake Erie cyanobacterial bloom. The Microcystis phycosphere communities were distinct from microbial communities in whole water and bulk phytoplankton seston in western Lake Erie but lacked 'core' taxa found across all colonies. However, dissimilarity in phycosphere community composition correlated with sampling date and the Microcystis 16S rRNA oligotype. Several taxa in the phycosphere were specific to and conserved with Microcystis of a single oligotype or sampling date. Together, this suggests that physiological differences between Microcystis strains, temporal changes in strain phenotypes, and the composition of seeding communities may impact community composition of the Microcystis phycosphere.

Shotgun metagenomic sequencing reveals the full taxonomic, trophic, and functional diversity of a coral reef benthic cyanobacterial mat from Bonaire, Caribbean Netherlands

Anthropogenic forcing is spurring cyanobacterial proliferation in aquatic ecosystems worldwide. While planktonic cyanobacterial blooms have received substantial research attention, benthic blooms of mat-forming cyanobacteria have received considerably less attention, especially benthic mat blooms on coral reefs. Resultingly, numerous aspects of coral reef benthic cyanobacterial bloom ecology remain unknown, including underlying biodiversity in the mat communities. Most previous characterizations of coral reef cyanobacterial mat composition have only considered the cyanobacterial component. Without an unbiased characterization of full community diversity, we cannot predict whole-community response to anthropogenic inputs or effectively determine appropriate mitigation strategies. Here, we advocate for the implementation of shotgun sequencing techniques to study coral reef cyanobacterial mats worldwide, utilizing a case study of a coral reef benthic cyanobacterial mat sampled from the island of Bonaire, Caribbean Netherlands. Read-based taxonomic profiling revealed that Cyanobacteria was present at only 47.57% relative abundance in a coral reef cyanobacterial mat, with non-cyanobacterial members of the sampled mat community, including diatoms (0.78%), fungi (0.25%), Archaea (0.34%), viruses (0.08%), and other bacteria (45.78%), co-dominating the community. We found numerous gene families for regulatory systems and for functional pathways (both aerobic and anaerobic). These gene families were involved in community coordination; photosynthesis; nutrient scavenging; and the cycling of sulfur, nitrogen, phosphorous, and iron. We also report bacteriophage (including prophage) sequences associated with this subtidal coral reef cyanobacterial mat, which could contribute to intra-mat nutrient cycling and bloom dynamics. Overall, our results suggest that Cyanobacteria-focused analysis of coral reef cyanobacterial mats underestimates mat diversity and fails to capture community members possessing broad metabolic potential for intra-mat nutrient scavenging, recycling, and retention that likely contribute to the contemporary success of cyanobacterial mats on reefs. We advocate for increased collaboration between microbiologists and coral reef ecologists to unite insights from each discipline and improve efforts to understand mat ecology.

Conditions that promote the formation of black bloom in aquatic microcosms and its effects on sediment bacteria related to iron and sulfur cycling

Black bloom occurs frequently in eutrophic waters. We investigated the conditions promoted the formation of black bloom via in-situ measurement in two aquatic microcosms and the effects of black bloom on the bacterial community composition. Although larger changes in dissolved oxygen (DO) were detected in the Hydrilla verticillata-dominated microcosm over the 90-day simulation, black bloom occurred more readily in the phytoplankton-dominated than macrophyte-dominated microcosm under conditions of O₂ depletion and temperature above 30 °C. The sediment bacterial community composition shifted after black bloom; the relative abundance of Thiobacillus and Sideroxydans, which oxidize iron (Fe) and sulfur (S), decreased by 47% and 48%, respectively, in the phytoplankton-dominated microcosm and by 18% and 20% in the macrophyte-dominated microcosm. By contrast, Desulfatiglans increased by 13% and 19%, respectively, after black bloom. Furthermore, inter-taxa correlations remarkably changed according to co-occurrence network analysis. Thirty-six different taxa from the phylum to the genus level were identified as biomarkers of sediments collected before and after the black bloom event. Most of these biomarkers are related to Fe/S cycling in aquatic ecosystems.

Alteration of dominant cyanobacteria in different bloom periods caused by abiotic factors and species interactions

Freshwater cyanobacterial blooms have drawn public attention because they threaten the safety of water resources and human health worldwide. Heavy cyanobacterial blooms outbreak in Lake Taihu in summer annually and vanish in other months. To find out the factors impacting the cyanobacterial blooms, the present study measured the physicochemical parameters of water and investigated the composition of microbial community using the 16S rRNA gene and internal transcribed spacer amplicon sequencing in the months with or without bloom. The most interesting finding is that two major cyanobacteria, Planktothrix and Microcystis, dramatically alternated during a cyanobacterial bloom in 2016, which is less mentioned in previous studies. When the temperature of the water began increasing in July, Planktothrix appeared first and showed as a superior competitor for M. aeruginosa in NO₃^--rich conditions. Microcystis became the dominant genus when the water temperature increased further in August. Laboratory experiments confirmed the influence of temperature and the total dissolved nitrogen (TDN) form on the growth of Planktothrix and Microcystis in a co-culture system. Besides, species interactions between cyanobacteria and non-cyanobacterial microorganisms, especially the prokaryotes, also played a key role in the alteration of Planktothrix and Microcystis. The present study exhibited the alteration of two dominant cyanobacteria in the different bloom periods caused by the temperature, TDN forms as well as the species interactions. These results helped the better understanding of cyanobacterial blooms and the factors which contribute to them.

Stability issues of microcystins, anabaenopeptins, anatoxins, and cylindrospermopsin during short-term and long-term storage of surface water and drinking water samples

Reproducible analytical procedures and rigorous quality control are imperative for an accurate monitoring of cyanobacterial toxins in environmental water samples. In this study, the short-term and long-term storage stability of diverse cyanotoxins (anatoxins, cylindrospermopsin, anabaenopeptins, and 12 microcystins) was evaluated in water samples, under different scenarios. Transport controls were performed at three monitoring sites in spiked ultrapure water and lake water to investigate short-term stability issues. Medium-term storage stability was evaluated for up to 14-28 days in ultrapure water, chlorine-treated drinking water (amended with reductant), and surface water (filtered and unfiltered) stored at different temperatures (20 °C, 4 °C, and -20 °C). Substantial decreases of cylindrospermopsin and anabaenopeptins were observed in tap water (20 °C) and unfiltered surface water (20 °C or 4 °C). Regardless of matrix type, cyanotoxin recoveries generally remained within an 80-120% range when the water samples were kept frozen. After a prolonged storage duration of 365 days at -20 °C, most cyanotoxins experienced decreases in the range of 10-20%. The notable exception was for the tryptophan-containing MC-LW and MC-WR, with more substantial variations (30% to 50% decrease) and conversion to N-formylkynurenine analogs. Reanalysis of field-collected surface waters after long-term storage at -20 °C also indicated significantly decreasing trends of cyanotoxins (between 6% and 23% decrease). In view of the above, short sample hold times should be favored as recommended in EPA methods.

Cyanobacterial blooms contribute to the diversity of antibiotic-resistance genes in aquatic ecosystems

Cyanobacterial blooms are a global ecological problem that directly threatens human health and crop safety. Cyanobacteria have toxic effects on aquatic microorganisms, which could drive the selection for resistance genes. The effect of cyanobacterial blooms on the dispersal and abundance of antibiotic-resistance genes (ARGs) of concern to human health remains poorly known. We herein investigated the effect of cyanobacterial blooms on ARG composition in Lake Taihu, China. The numbers and relative abundances of total ARGs increased obviously during a Planktothrix bloom. More pathogenic microorganisms were present during this bloom than during a Planktothrix bloom or during the non-bloom period. Microcosmic experiments using additional aquatic ecosystems (an urban river and Lake West) found that a coculture of Microcystis aeruginosa and Planktothrix agardhii increased the richness of the bacterial community, because its phycosphere provided a richer microniche for bacterial colonization and growth. Antibiotic-resistance bacteria were naturally in a rich position, successfully increasing the momentum for the emergence and spread of ARGs. These results demonstrate that cyanobacterial blooms are a crucial driver of ARG diffusion and enrichment in freshwater, thus providing a reference for the ecology and evolution of ARGs and ARBs and for better assessing and managing water quality.

Seasonal dynamics of the bacterial communities associated with cyanobacterial blooms in the Han River

DNA-based analyses of bacterial communities were performed to identify the bacteria co-occurring with cyanobacterial blooms in samples collected at a single site over 2 years. Microcystis aeruginosa was the most predominant species (81% in 2018, and 94% in 2019) within the phylum Cyanobacteria, and microcystins were detected during all cyanobacterial blooms. The stereo microscope and scanning electron microscope observations showed bacterial associations on and around the aggregated M. aeruginosa cells. Culture-independent analyses of filtered bacterial communities showed that the Flavobacterium species in phylum Bacteroidetes (19%) was dominant in the cyanobacterial phycosphere, followed by the Limnohabitans species in Betaproteobacteria (11%). Using principal component analysis, major bacterial genus, including Microcystis and Flavobacterium species, were clustered during cyanobacterial blooms in both years. To identify key bacterial species that develop long-term symbiosis with M. aeruginosa, another culture-independent analysis was performed after the environmental sample had been serially subcultured for 1 year. Interestingly, Brevundimonas (14%) was the most dominant species, followed by Porphyrobacter (7%) and Rhodobacter (3.5%) within the Alphaproteobacteria. Screening of 100 colonies from cyanobacterial bloom samples revealed that the majority of culturable bacteria belonged to Gammaproteobacteria (28%) and Betaproteobacteria (57%), including Pseudomonas, Curvibacter, and Paucibacter species. Several isolates of Brevundimonas, Curvibacter, and Pseudomonas species could promote the growth of axenic M. aeruginosa PCC7806. The sensitivity of M. aeruginosa PCC7806 cells to different environmental conditions was monitored in bacteria-free pristine freshwater, indicating that nitrogen addition promotes the growth of M. aeruginosa.

Covariation patterns of phytoplankton and bacterioplankton in hypertrophic shallow lakes

The aim of this work was to assess the temporal patterns in the community composition of phytoplankton (PCC) and bacterioplankton (BCC) in two interconnected and hypertrophic Pampean shallow lakes in Argentina. Factors shaping their community dynamics and community temporal covariations were also analysed. We performed 4 years of seasonal samplings (2012-2016) and communities were studied by the Utermöhl approach (PCC) and Illumina MiSeq sequencing (BCC). We found marked seasonal variations in both communities and inter-annual variations with decreasing microbial community similarities during the study. We also observed covariation in community-level dynamics among PCC and BCC within and between shallow lakes. The within-lake covariations remained positive and significant, while controlling for the effects of intrinsic (environmental) and extrinsic (temporal and meteorological) factors, suggesting a community coupling mediated by intrinsic biotic interactions. Algal-bacterial associations between different taxa of phytoplankton and bacterioplankton within each lake were also found. PCC was mainly explained by pure regional extrinsic (17-21%) and intrinsic environmental (8-9%) factors, while BCC was explained by environmental (8-10%) and biotic interactions with phytoplankton (7-8%). Our results reveal that the influence of extrinsic regional factors can be channeled to bacterioplankton through both environmental (i.e. water temperature) and phytoplankton effects.

Impact of algal organic matter on the performance, cyanotoxin removal, and biofilms of biologically-active filtration systems

The occurrence of harmful algal blooms dominated by toxic cyanobacteria has induced continuous loadings of algal organic matter (AOM) and toxins in drinking water treatment plants. However, the impact of AOM on the active biofilms and microbial community structures of biologically-active filtration (BAF), which directly affects the contaminant removal, is not well understood. In this study, we systematically examined the effects of AOM on BAF performance and bacterial biofilm formation over 240 days, tracing the removal of specific AOM components, a cyanotoxin [microcystin-LR (MC-LR)], and microbial community responses. The component analysis (excitation and emission matrix analysis) results for AOM revealed that terrestrial humic-like substances showed the highest removal among all the identified components and were strongly correlated to MC-LR removal. In addition, reduced empty bed contact time and deactivation of biofilms significantly decreased BAF performances for both AOM and MC-LR. The active biofilm, bacterial community structure, and mlrA gene (involved in microcystin degradation) abundance demonstrated that bacterial biofilm composition responded to AOM and MC-LR, in which Rhodocyclaceae, Saprospiraceae, and Comamonadaceae were dominant. In addition, MC-LR biodegradation appeared to be more active at the top than at the bottom layer in BAF. Overall, this study provides deeper insights into the role of biofilms and filter operation on the fate of AOM and MC-LR in BAF.

The Composition and Function of Microbiomes Within Microcystis Colonies Are Significantly Different Than Native Bacterial Assemblages in Two North American Lakes

The toxic cyanobacterium Microcystis is one of the most pervasive harmful algal bloom (HAB) genera and naturally occurs in large colonies known to harbor diverse heterotrophic bacterial assemblages. While colony-associated microbiomes may influence Microcystis blooms, there remains a limited understanding of the structure and functional potential of these communities and how they may be shaped by changing environmental conditions. To address this gap, we compared the dynamics of Microcystis-attached (MCA), free-living (FL), and whole water (W) microbiomes during Microcystis blooms using next-generation amplicon sequencing (16S rRNA), a predictive metagenome software, and other bioinformatic approaches. Microbiomes were monitored through high resolution spatial-temporal surveys across two North American lakes, Lake Erie (LE) and Lake Agawam (LA; Long Island, NY, United States) in 2017, providing the largest dataset of these fractions to date. Sequencing of 126 samples generated 7,922,628 sequences that clustered into 7,447 amplicon sequence variants (ASVs) with 100% sequence identity. Across lakes, the MCA microbiomes were significantly different than the FL and W fractions being significantly enriched in Gemmatimonadetes, Burkholderiaceae, Rhizobiales, and Cytophagales and depleted of Actinobacteria. Further, although MCA communities harbored > 900 unique ASVs, they were significantly less diverse than the other fractions with diversity inversely related to bloom intensity, suggesting increased selection pressure on microbial communities as blooms intensified. Despite taxonomic differences between lakes, predicted metagenomes revealed conserved functional potential among MCA microbiomes. MCA communities were significantly enriched in pathways involved in N and P cycling and microcystin-degradation. Taxa potentially capable of N2-fixation were significantly enriched (p < 0.05) and up to four-fold more abundant within the MCA faction relative to other fractions, potentially aiding in the proliferation of Microcystis blooms during low N conditions. The MCA predicted metagenomes were conserved over 8 months of seasonal changes in temperature and N availability despite strong temporal succession in microbiome composition. Collectively, these findings indicate that Microcystis colonies harbor a statistically distinct microbiome with a conserved functional potential that may help facilitate bloom persistence under environmentally unfavorable conditions.

Network analysis reveals succession of Microcystis genotypes accompanying distinctive microbial modules with recurrent patterns

Every member of the ecological community is connected via a network of vital and complex relationships, called the web of life. To elucidate the ecological network and interactions among producers, consumers, and decomposers in the Daechung Reservoir, Korea, during cyanobacterial harmful algal blooms (cyanoHAB), especially those involving Microcystis, we investigated the diversity and compositions of the cyanobacterial (16S rRNA gene), including the genotypes of Microcystis (cpcBA-IGS gene), non-cyanobacterial (16S), and eukaryotic (18S) communities through high-throughput sequencing. Microcystis blooms were divided into the Summer Major Bloom and Autumn Minor Bloom with different dominant genotypes of Microcystis. Network analysis demonstrated that the modules involved in the different phases of the Microcystis blooms were categorized into the Pre-Bloom, Bloom, Post-Bloom, and Non-Bloom Groups at all sampling stations. In addition, the non-cyanobacterial components of each Group were classified, while the same Group showed similarity across all stations, suggesting that Microcystis and other microbes were highly interdependent and organized into cyanoHAB-related module units. Importantly, the Microcystis genotype-based sub-network uncovered that Pirellula, Pseudanabaena, and Vampirovibrionales preferred to interact with specific Microcystis genotypes in the Summer Major Bloom than with other genotypes in the Autumn Minor Bloom, while the copepod Skistodiaptomus exhibited the opposite pattern. In conclusion, the transition patterns of cyanoHAB-related modules and their key components could be crucial in the succession of Microcystis genotypes and to enhance the understanding of microbial ecology in an aquatic environment.

Broad and Fine Scale Variability in Bacterial Diversity and Cyanotoxin Quotas in Benthic Cyanobacterial Mats

Benthic proliferations of Microcoleus autumnalis (basionym Phormidium autumnale) and closely related taxa are being reported with increasing frequency in streams and rivers worldwide. This species commonly produces the potent neurotoxin anatoxin, and exposure to this has resulted in animal fatalities and human health concerns. Bacterial communities within cyanobacterial assemblages can facilitate processes such as nutrient cycling and are posited to influence cyanobacterial growth and function. However, there is limited knowledge on spatial variability of bacterial communities associated with benthic cyanobacteria and anatoxin content and quotas. In this study, M. autumnalis-dominated mat samples were collected from six sites in two New Zealand streams. Associated bacterial communities were characterized using 16S rRNA metabarcoding, anatoxin content by liquid chromatography-mass spectrometry and anaC copies using droplet digital PCR. Bacterial assemblages differed significantly when amplicon sequence variants were compared between streams and most sites within streams. These differences were associated with conductivity, DRP, DIN, temperature, anatoxin concentration, and quota. Despite the differences in bacterial community composition; at phyla, class and order levels there was high similarity across spatial scales, with Bacteroidetes (ca. 67%) and Proteobacteria (ca. 25%) dominant. There was significant variability in total anatoxin concentrations between sites in both streams (p < 0.001). When the data were converted to anatoxin quotas variability was reduced, suggesting that the relative abundance of toxic genotypes is a key driver of total anatoxin concentrations in mats. This study demonstrates the complexity of microbial communities within M. autumnalis-dominated mats and highlights their likely important role in within-mat nutrient cycling processes.

Drivers and ecological consequences of dominance in periurban phytoplankton communities using networks approaches

Evaluating the causes and consequences of dominance by a limited number of taxa in phytoplankton communities is of huge importance in the current context of increasing anthropogenic pressures on natural ecosystems. This is of particular concern in densely populated urban areas where usages and impacts of human populations on water ecosystems are strongly interconnected. Microbial biodiversity is commonly used as a bioindicator of environmental quality and ecosystem functioning, but there are few studies at the regional scale that integrate the drivers of dominance in phytoplankton communities and their consequences on the structure and functioning of these communities. Here, we studied the causes and consequences of phytoplankton dominance in 50 environmentally contrasted waterbodies, sampled over four summer campaigns in the highly-populated Île-de-France region (IDF). Phytoplankton dominance was observed in 32-52% of the communities and most cases were attributed to Chlorophyta (35.5-40.6% of cases) and Cyanobacteria (30.3-36.5%). The best predictors of dominance were identified using multinomial logistic regression and included waterbody features (surface, depth and connection to the hydrological network) and water column characteristics (total N, TN:TP ratio, water temperature and stratification). The consequences of dominance were dependent on the identity of the dominant organisms and included modifications of biological attributes (richness, cohesion) and functioning (biomass, RUE) of phytoplankton communities. We constructed co-occurrence networks using high resolution phytoplankton biomass and demonstrated that networks under dominance by Chlorophyta and Cyanobacteria exhibited significantly different structure compared with networks without dominance. Furthermore, dominance by Cyanobacteria was associated with more profound network modifications (e.g. cohesion, size, density, efficiency and proportion of negative links), suggesting a stronger disruption of the structure and functioning of phytoplankton communities in the conditions in which this group dominates. Finally, we provide a synthesis on the relationships between environmental drivers, dominance status, community attributes and network structure.

Characterization of Distinct CyanoHABs-Related Modules in Microbial Recurrent Association Network

To elucidate the interspecies connectivity between cyanobacteria and other bacteria (non-cyanobacteria) during cyanobacterial harmful algal blooms (cyanoHABs), samples were collected from the Nakdong River, Korea, from June 2016 to August 2017, and microbial recurrent association network (MRAN) analysis was performed to overcome the limitations of conventional network analysis. Microcystis blooms were tightly linked with Pseudanabaena in summer and were accompanied by significant changes in the non-cyanobacterial community composition (nCCC) compared to non-bloom period. Riverine bacterial communities could be clearly separated into modules that were involved in the formation, maintenance, and decomposition of cyanoHABs. Roseomonas and Herbaspirillum were directly linked with major cyanobacteria and assigned to connector and module hub in cyanoHABs-related modules, respectively. The functional profiles of the cyanoHABs-related modules suggested that nitrate reduction, aerobic ammonia oxidation, fermentation, and hydrocarbon degradation could be increased during the Microcystis bloom periods. In conclusion, MRAN analysis revealed that specific bacteria belonging to cyanoHABs-related module, including connectors and module hubs, appeared to contribute to the development and collapse of cyanoHABs. Therefore, to understand cyanoHABs, a modular microbial perspective may be more helpful than a single bacterial species perspective.

Phosphate availability affects fixed nitrogen transfer from diazotrophs to their epibionts

Dinitrogen (N₂) fixation is a major source of external nitrogen (N) to aquatic ecosystems and therefore exerts control over productivity. Studies have shown that N₂ -fixers release freshly fixed N into the environment, but the causes for this N release are largely unclear. Here, we show that the availability of phosphate can directly affect the transfer of freshly fixed N to epibionts in filamentous, diazotrophic cyanobacteria. Stable-isotope incubations coupled to single-cell analyses showed that <1% and ~15% of freshly fixed N was transferred to epibionts of Aphanizomenon and Nodularia, respectively, at phosphate scarcity during a summer bloom in the Baltic Sea. When phosphate was added, the transfer of freshly fixed N to epibionts dropped to about half for Nodularia, whereas the release from Aphanizomenon increased slightly. At the same time, the growth rate of Nodularia roughly doubled, indicating that less freshly fixed N was released and was used for biomass production instead. Phosphate scarcity and the resulting release of freshly fixed N could explain the heavy colonization of Nodularia filaments by microorganisms during summer blooms. As such, the availability of phosphate may directly affect the partitioning of fixed N₂ in colonies of diazotrophic cyanobacteria and may impact the interactions with their microbiome.

Impacts of microbial assemblage and environmental conditions on the distribution of anatoxin-a producing cyanobacteria within a river network

Blooms of planktonic cyanobacteria have long been of concern in lakes, but more recently, harmful impacts of riverine benthic cyanobacterial mats been recognized. As yet, we know little about how various benthic cyanobacteria are distributed in river networks, or how environmental conditions or other associated microbes in their consortia affect their biosynthetic capacities. We performed metagenomic sequencing for 22 Oscillatoriales-dominated (Cyanobacteria) microbial mats collected across the Eel River network in Northern California and investigated factors associated with anatoxin-a producing cyanobacteria. All microbial communities were dominated by one or two cyanobacterial species, so the key mat metabolisms involve oxygenic photosynthesis and carbon oxidation. Only a few metabolisms fueled the growth of the mat communities, with little evidence for anaerobic metabolic pathways. We genomically defined four cyanobacterial species, all which shared <96% average nucleotide identity with reference Oscillatoriales genomes and are potentially novel species in the genus Microcoleus. One of the Microcoleus species contained the anatoxin-a biosynthesis genes, and we describe the first anatoxin-a gene cluster from the Microcoleus clade within Oscillatoriales. Occurrence of these four Microcoleus species in the watershed was correlated with total dissolved nitrogen and phosphorus concentrations, and the species that contains the anatoxin-a gene cluster was found in sites with higher nitrogen concentrations. Microbial assemblages in mat samples with the anatoxin-a gene cluster consistently had a lower abundance of Burkholderiales (Betaproteobacteria) species than did mats without the anatoxin-producing genes. The associations of water nutrient concentrations and certain co-occurring microbes with anatoxin-a producing Microcoleus motivate further exploration for their roles as potential controls on the distributions of toxigenic benthic cyanobacteria in river networks.

Light Modulates the Physiology of Nonphototrophic Actinobacteria

Light is a source of energy and an environmental cue that is available in excess in most surface environments. In prokaryotic systems, conversion of light to energy by photoautotrophs and photoheterotrophs is well understood, but the conversion of light to information and the cellular response to that information have been characterized in only a few species. Our goal was to explore the response of freshwater Actinobacteria, which are ubiquitous in illuminated aquatic environments, to light. We found that Actinobacteria without functional photosystems grow faster in the light, likely because sugar transport and metabolism are upregulated in the light. Based on the action spectrum of the growth effect and comparisons of the genomes of three Actinobacteria with this growth rate phenotype, we propose that the photosensor in these strains is a putative CryB-type cryptochrome. The ability to sense light and upregulate carbohydrate transport during the day could allow these cells to coordinate their time of maximum organic carbon uptake with the time of maximum organic carbon release by primary producers.IMPORTANCE Sunlight provides information about both place and time. In sunlit aquatic environments, primary producers release organic carbon and nitrogen along with other growth factors during the day. The ability of Actinobacteria to coordinate organic carbon uptake and utilization with production of photosynthate enables them to grow more efficiently in the daytime, and it potentially gives them a competitive advantage over heterotrophs that constitutively produce carbohydrate transporters, which is energetically costly, or produce transporters only after detection of the substrate(s), which delays their response. Understanding how light cues the transport of organic carbon and its conversion to biomass is key to understanding biochemical mechanisms within the carbon cycle, the fluxes through it, and the variety of mechanisms by which light enhances growth.

High-throughput DNA sequencing reveals the dominance of pico- and other filamentous cyanobacteria in an urban freshwater Lake

The current study presents findings related to algal blooms in a fresh water lake, which has been experiencing severe cyanobacterial blooms (CyanoHABs). Primarily, picocyanobacteria belonging to the genus Synechococcus and filamentous cyanobacterial group belonging to Aphanizomenon and Dolichospermum dominated top water column during non-bloom and bloom periods respectively. The dominance of Synechococcus in early summer informs that blooming in Utah Lake starts in early summer and then later is taken over by other bloom-forming cyanobacteria, such as species belonging to the genus Aphanizomenon. A strong negative correlation (r = -0.9, p < 0.001) was found between the occurrence of Aphanizomenon and Synechococcus which correlates very well with the fact that the blooms of these two different cyanobacteria never coexisted. The predominance of cyanobacteria in 2017 was attributed more to temperature (r = 0.18, p < 0.001). The Actinobacteria was negatively correlated with primary production and high chlorophyll a concentration. Flavobacterium and Limnohabitans were the main phytoplankton colonizers and predators detected that could secrete extracellular enzymes to degrade algal exudates (such as proteins and polysaccharides). Additionally, cyanotoxins producers Microcystis aeruginosa and Planktothrix accounted for up to 12.43% and 7.04% of total cyanobacteria abundance during blooms. The relative abundance of chloroplast reads was overall lower than the cyanobacteria reads, except for the May 5th sampling in 2017. There was inter-annual variability in the bloom-associated heterotrophic bacterial populations, but these populations were consistent with bloom-associated bacterial populations found in other lakes. Community diversity analysis for both Shannon and Simpson indices indicated lower community diversity during the bloom period. The beta diversity conducted by PCoA and UPGMA trees suggested the significant temporal rather than spatial impacts on shaping the phytoplankton community structures during the summer season.

Diversity of Bacteria in Lakes with Different Chlorophyll Content and Investigation of Their Respiratory Activity through a Long-Term Microcosm Experiment

Bacterial community structure and metabolism are critical factors for ecosystem functioning since they affect remineralization of nutrients and carbon flow. We used Illumina sequencing of 16SrRNA V3-V4 regions to investigate whether bacterial assemblage composition differs between four samples from two lakes in the geographic region of Epirus (Greece) characterized by distinct oligotrophic to eutrophic/hypereutrophic conditions as revealed by chlorophyll-a values. We found high similarity (>60%) for bacterial assemblages recovered from the two lakes when eutrophic/hypereutrophic conditions prevailed. Distinct bacterial communities appeared in oligotrophic and mesotrophic waters. Low temperature was occasionally an important factor in shaping the bacterial community. In parallel, microcosm experiments were performed to estimate respiration rates of bacterioplankton at in situ temperature and under a 2 °C temperature increase scenario. Differently assembled communities were found to display similar rates except under hypereutrophic conditions when respiration increased significantly, leading to hypoxic conditions. Temperature increase did not affect respiration rates. Overall this study indicated a clear differentiation of bacterial communities between sites of different trophic state. However, different communities responded similarly under a specific range of chlorophyll-a values and resisted small scale temperature perturbations. Different results were found for hypereutrophic conditions and this has implications for ecosystems functioning, given the increasing occurrence of eutrophication events.

The Characteristics and Dynamics of Cyanobacteria-Heterotrophic Bacteria Between Two Estuarine Reservoirs - Tropical Versus Sub-Tropical Regions

In this study, Illumina MiSeq sequencing technique was employed to explore the characteristics and dynamics of cyanobacteria–heterotrophic bacteria between two estuarine reservoirs in sub-tropical (reservoir A in Shanghai) and tropical (reservoir B in Singapore) regions. The results indicated that significant differences in bacterial community composition were found between two estuarine reservoirs, which influenced by varied environmental variables. The environmental heterogeneity in reservoir A was much higher, which indicated that the composition of bacterial community in reservoir A was more complex. In contrast, reservoir B provided a suitable and temperate water environment conditions for bacterial growth, which resulted in higher community diversity and less co-exclusion correlations. The molecular ecological network indicated that the presence of dominant bacterial community in each of the reservoir were significant different. These differences mainly reflected the responses of bacterial community to the variations of environmental variables. Although Synechococcus was the dominant cyanobacterial species in both reservoirs, it exhibited co-occurrence patterns with different heterotrophic bacteria between reservoirs. In addition, the cyanobacteria–heterotrophic bacteria interaction exhibited highly dynamic variations, which was affected by nutrition and survive space. Also, the co-occurrence of Microcystis and Pseudanabaena found in reservoir B implied that the non-N-fixing Microcystis accompanied with N-fixing Pseudanabeana occurrence in freshwater lakes, so as to better meet the demand for nitrogen source.