The Trait Repertoire Enabling Cyanobacteria to Bloom Assessed through Comparative Genomic Complexity and Metatranscriptomics

Insights into cyanobacterial blooms through the lens of omics

Cyanobacteria are oxygen-producing photosynthetic bacteria that convert carbon dioxide into biomass upon exposure to sunlight. However, favorable conditions cause harmful cyanobacterial blooms (HCBs), which are the dense accumulation of biomass at the water surface or subsurface, posing threats to freshwater ecosystems and human health. Understanding the mechanisms underlying cyanobacterial bloom formation is crucial for effective management. In this regard, recent advancements in omics technologies have provided valuable insights into HCBs, which have raised expectations to develop more effective control methods in the near future. This literature review aims to present the genomic architecture, adaptive mechanisms, microbial interactions, and ecological impacts of HCBs through the lens of omics. Genomic analysis indicates that the genome plasticity of cyanobacteria has enabled their resilience and effective adaptation to environmental changes. Transcriptomic investigations have revealed that cyanobacteria use various strategies for adapting to environmental stress. Additionally, metagenomic and metatranscriptomic analyses have emphasized the significant role of the microbial community in regulating HCBs. Finally, we offer perspectives on potential opportunities for further research in this field.

An evolutionary optimum amid moderate heritability in prokaryotic cell size

We investigate the distribution and evolution of prokaryotic cell size based on a compilation of 5,380 species. Size spans four orders of magnitude, from 100 nm (Mycoplasma) to more than 1 cm (Thiomargarita); however, most species congregate heavily around the mean. The distribution approximates but is distinct from log normality. Comparative phylogenetics suggests that size is heritable, yet the phylogenetic signal is moderate, and the degree of heritability is independent of taxonomic scale (i.e., fractal). Evolutionary modeling indicates the presence of an optimal cell size to which most species gravitate. The size is equivalent to a coccus of 0.70 μm in diameter. Analyses of 1,361 species with sequenced genomes show that genomic traits contribute to size evolution moderately and synergistically. Given our results, scaling theory, and empirical evidence, we discuss potential drivers that may expand or shrink cells around the optimum and propose a stability landscape model for prokaryotic cell size.

Woronichinia naegeliana: A common nontoxigenic component of temperate freshwater cyanobacterial blooms with 30% of its genome in transposons

Monitoring in the U.S. state of Washington across the period 2007-2019 showed that Woronichinia has been present in many lakes state-wide. This cyanobacterium was commonly dominant or sub-dominant in cyanobacterial blooms in the wet temperate region west of the Cascade Mountains. In these lakes, Woronichinia often co-existed with Microcystis, Dolichospermum and Aphanizomenon flos-aquae and the cyanotoxin microcystin has often been present in those blooms, although it has not been known whether Woronichinia is a toxin producer. We report the first complete genome of Woronichinia naegeliana WA131, assembled from the metagenome of a sample collected from Wiser Lake, Washington, in 2018. The genome contains no genes for cyanotoxin biosynthesis or taste-and-odor compounds, but there are biosynthetic gene clusters for other bioactive peptides, including anabaenopeptins, cyanopeptolins, microginins and ribosomally produced, post-translationally modified peptides. Genes for photosynthesis, nutrient acquisition, vitamin synthesis and buoyancy that are typical of bloom-forming cyanobacteria are present, although nitrate and nitrite reductase genes are conspicuously absent. However, the 7.9 Mbp genome is 3-4 Mbp larger than those of the above-mentioned frequently co-existing cyanobacteria. The increased genome size is largely due to an extraordinary number of insertion sequence elements (transposons), which account for 30.3% of the genome and many of which are present in multiple copies. The genome contains a relatively large number of pseudogenes, 97% of which are transposase genes. W. naegeliana WA131 thus seems to be able to limit the potentially deleterious effects of high rates of recombination and transposition to the mobilome fraction of its genome.

The Special and General Mechanism of Cyanobacterial Harmful Algal Blooms

Cyanobacterial harmful algal blooms (CyanoHABs) are longstanding aquatic hazards worldwide, of which the mechanism is not yet fully understood, i.e., the process in which cyanobacteria establish dominance over coexisting algae in the same eutrophic waters. The dominance of CyanoHABs represents a deviation from their low abundance under conventional evolution in the oligotrophic state, which has been the case since the origin of cyanobacteria on early Earth. To piece together a comprehensive mechanism of CyanoHABs, we revisit the origin and adaptive radiation of cyanobacteria in oligotrophic Earth, demonstrating ubiquitous adaptive radiation enabled by corresponding biological functions under various oligotrophic conditions. Next, we summarize the biological functions (ecophysiology) which drive CyanoHABs and ecological evidence to synthesize a working mechanism at the population level (the special mechanism) for CyanoHABs: CyanoHABs are the consequence of the synergistic interaction between superior cyanobacterial ecophysiology and elevated nutrients. Interestingly, these biological functions are not a result of positive selection by water eutrophication, but an adaptation to a longstanding oligotrophic state as all the genes in cyanobacteria are under strong negative selection. Last, to address the relative dominance of cyanobacteria over coexisting algae, we postulate a "general" mechanism of CyanoHABs at the community level from an energy and matter perspective: cyanobacteria are simpler life forms and thus have lower per capita nutrient demand for growth than coexisting eukaryotic algae. We prove this by comparing cyanobacteria and eukaryotic algae in cell size and structure, genome size, size of genome-scale metabolic networks, cell content, and finally the golden standard-field studies with nutrient supplementation in the same waters. To sum up, the comprehensive mechanism of CyanoHABs comprises a necessary condition, which is the general mechanism, and a sufficient condition, which is the special mechanism. One prominent prediction based on this tentative comprehensive mechanism is that eukaryotic algal blooms will coexist with or replace CyanoHABs if eutrophication continues and goes over the threshold nutrient levels for eukaryotic algae. This two-fold comprehensive mechanism awaits further theoretic and experimental testing and provides an important guide to control blooms of all algal species.

Cyanobacterial Blooms Are Not a Result of Positive Selection by Freshwater Eutrophication

Long-standing cyanobacterial harmful algal blooms (CyanoHABs) are known to result from synergistic interaction between elevated nutrients and superior ecophysiology of cyanobacteria. However, it remains to be determined whether CyanoHABs are a result of positive selection by eutrophic waters. To address this, we conducted molecular evolutionary analyses on the genomes of 9 bloom-forming cyanobacteria, combined with pangenomics and metatranscriptomics. The results showed no positive selection by water eutrophication. Instead, all homologous genes in the species are under strong purifying selection based on the ratio of divergence at nonsynonymous and synonymous sites (dN/dS) and phylogeny. The dN/dS < 0.85 (median = 0.3) for all homologous genes are similar between the genes in the pathways driving CyanoHABs and housekeeping functions. Phylogenetic support for non-positive selection comes from the mixed clustering of strains: strains of the same species from diverse geographic origins form the same clusters, while strains from the same origins form different clusters. Further support lies in the codon adaptation index (CAI) and single nucleotide polymorphism (SNP). The CAI ranged from 0.42 to 0.9 (mean = 0.75), which indicates high-level codon usage bias; the pathways for CyanoHABs and housekeeping functions showed a similar CAI. Interestingly, CAI was negatively correlated with gene expression in 3 metatranscriptomes. The numbers of SNPs were concentrated around 5 to 50. As the SNP number increases, the gene expression level decreases. These negative correlations agree with the population-level dN/dS and phylogeny in supporting purifying selection in bloom-forming cyanobacteria. In summary, superior ecophysiology appears to be acquired prior to water eutrophication. IMPORTANCE CyanoHABs are global environmental hazards, and their mechanisms of action are being intensively investigated. On an ecological scale, CyanoHABs are consequences of synergistic interactions between biological functions and elevated nutrients in eutrophic waters. On an evolutionary scale, one important question is how bloom-forming cyanobacteria acquire these superior biological functions. There are several possibilities, including adaptive evolution and horizontal gene transfer. Here, we explored the possibility of positive selection. We reasoned that there are two possible periods for cyanobacteria to acquire these functions: before the onset of water eutrophication or during water eutrophication. Either way, there should be molecular signatures in protein sequences for positive selection. Interestingly, we found no positive selection by water eutrophication, but strong purifying selection instead on nearly all the genes, suggesting these superior functions aiding CyanoHABs are acquired prior to water eutrophication.

Transcriptomic responses to phosphorus in an invasive cyanobacterium, Raphidiopsis raciborskii: Implications for nutrient management

Phosphorus (P) is a vital macronutrient associated with the growth and proliferation of Raphidiopsis raciborskii, an invasive and notorious bloom-forming cyanobacterium. However, the molecular mechanisms involved in P acclimation remain largely unexplored for Raphidiopsis raciborskii. Here, transcriptome sequencing of Raphidiopsis raciborskii was conducted to reveal multifaceted mechanisms involved in mimicking dipotassium phosphate (DIP), β-glycerol phosphate (Gly), 2-aminoethylphosphonic acid (AEP), and P-free conditions (NP). Chlorophyll a fluorescence parameters showed significant differences in the NP and AEP groups compared with the DIP and Gly-groups. Expression levels of genes related to phosphate transportation and uptake, organic P utilization, nitrogen metabolism, urea cycling, carbon fixation, amino acid metabolism, environmental information, the ATP-synthesis process in glycolysis, the tricarboxylic acid (TCA) cycle, and the pentose phosphate pathway were remarkably upregulated, while those related to photosynthesis, phycobiliproteins, respiration, oxidative phosphorylation, sulfur metabolism, and genetic information were markedly downregulated in the NP group relative to the DIP group. However, the expression of genes involved in organic P utilization, the urea cycle, and genetic information in the Gly-group, and carbon-phosphorus lyase, genetic information and environmental information in the AEP group were significantly increased compared to the DIP group. Together, these results indicate that Raphidiopsis raciborskii exhibits the evolution of coordination of multiple metabolic pathways and certain key genes to adapt to ambient P changes, which implies that if P is reduced to control Raphidiopsis raciborskii bloom, there is a risk that external nutrients (such as nitrogen, amino acids, and urea) will stimulate the growth or metabolism of Raphidiopsis.

Cyanobacterial dominance and succession: Factors, mechanisms, predictions, and managements

Eutrophication of natural water bodies worldwide has led to cyanobacteria becoming the dominant species in phytoplankton communities, causing serious harm environmentally and economically. Cyanobacterial succession makes effective treatment of cyanobacterial blooms a challenge. Although there are many studies about cyanobacterial dominance and succession, it is still lack of relevant review summarizing the advances on this topic. To control cyanobacterial blooms and manage water quality effectively, we conducted a critical review and drew the following conclusions: (1) cyanobacterial dominance and succession occur from spring to summer, with changes of multiple environmental factors dominated by temperature and nutrients conditions; (2) the cyanobacterial dominance and succession are inherently attributed to the distinctive traits of cyanobacteria including colony formation, gas vesicles, toxin release, and nitrogen fixation; (3) given the current meta-omics explorations on mechanisms of cyanobacterial succession, how to combine the extensive data to draw general conclusions is a challenge in the future; (4) the dominant niche of high temperature-adapted cyanobacteria genera will be further reinforced with global warming and elevated carbon dioxide in the future; (5) considering the causes and future developments of cyanobacterial blooms, the management strategies for controlling cyanobacterial blooms include reducing external nutrient input and removing internal nutrient in sediment, artificial mixing waters to decrease buoyancy of cyanobacteria, and biological control using allelopathy of aquatic plants and/or enhancing zooplankton feeding.

Nutrient Loading and Viral Memory Drive Accumulation of Restriction Modification Systems in Bloom-Forming Cyanobacteria

The mechanisms driving cyanobacterial harmful algal blooms (HABs) like those caused by Microcystis aeruginosa remain elusive, but improved defense against viral predation has been implicated for success in eutrophic environments. Our genus-level analyses of 139,023 genomes revealed that HAB-forming cyanobacteria carry vastly more restriction modification systems per genome (RMPG) than nearly all other prokaryotic genera, suggesting that viral defense is a cornerstone of their ecological success. In contrast, picocyanobacteria that numerically dominate nutrient-poor systems have the fewest RMPG within the phylum Cyanobacteria. We used classic resource competition models to explore the hypothesis that nutrient enrichments drive ecological selection for high RMPG due to increased host-phage contact rate. These classic models, agnostic to the mechanism of defense, explain how nutrient loading can select for increased RMPG but, importantly, fail to explain the extreme accumulation of these defense systems. However, extreme accumulation of RMPG can be achieved in a novel “memory” model that accounts for a unique activity of restriction modification systems: the accidental methylation of viral DNA by the methyltransferase. The methylated virus “remembers” the RM defenses of its former host and can evade these defenses if they are present in the next host. This viral memory leads to continual RM system devaluation; RMs accumulate extensively because the benefit of each addition is diminished. Our modeling leads to the hypothesis that nutrient loading and virion methylation drive the extreme accumulation of RMPG in HAB-forming cyanobacteria. Finally, our models suggest that hosts with different RMPG values can coexist when hosts have unique sets of RM systems.

Genetic, Genomics, and Responses to Stresses in Cyanobacteria: Biotechnological Implications

Cyanobacteria are widely-diverse, environmentally crucial photosynthetic prokaryotes of great interests for basic and applied science. Work to date has focused mostly on the three non-nitrogen fixing unicellular species Synechocystis PCC 6803, Synechococcus PCC 7942, and Synechococcus PCC 7002, which have been selected for their genetic and physiological interests summarized in this review. Extensive "omics" data sets have been generated, and genome-scale models (GSM) have been developed for the rational engineering of these cyanobacteria for biotechnological purposes. We presently discuss what should be done to improve our understanding of the genotype-phenotype relationships of these models and generate robust and predictive models of their metabolism. Furthermore, we also emphasize that because Synechocystis PCC 6803, Synechococcus PCC 7942, and Synechococcus PCC 7002 represent only a limited part of the wide biodiversity of cyanobacteria, other species distantly related to these three models, should be studied. Finally, we highlight the need to strengthen the communication between academic researchers, who know well cyanobacteria and can engineer them for biotechnological purposes, but have a limited access to large photobioreactors, and industrial partners who attempt to use natural or engineered cyanobacteria to produce interesting chemicals at reasonable costs, but may lack knowledge on cyanobacterial physiology and metabolism.

Complete sequence and structure of the genome of the harmful algal bloom-forming cyanobacterium Planktothrix agardhii NIES-204T and detailed analysis of secondary metabolite gene clusters

Planktothrix species are distributed worldwide, and these prevalent cyanobacteria occasionally form potentially devastating toxic blooms. Given the ecological and taxonomic importance of Planktothrix agardhii as a bloom species, we set out to determine the complete genome sequence of the type strain Planktothrix agardhii NIES-204. Remarkably, we found that the 5S ribosomal RNA genes are not adjacent to the 16S and 23S ribosomal RNA genes. The genomic structure of P. agardhii NIES-204 is highly similar to that of another P. agardhii strain isolated from a geographically distant site, although they differ distinctly by a large inversion. We identified numerous gene clusters that encode the components of the metabolic pathways that generate secondary metabolites. We found that the aeruginosin biosynthetic gene cluster was more similar to that of another toxic bloom-forming cyanobacterium Microcystis aeruginosa than to that of other strains of Planktothrix, suggesting horizontal gene transfer. Prenyltransferases encoded in the prenylagaramide gene cluster of Planktothrix strains were classified into two phylogenetically distinct types, suggesting a functional difference. In addition to the secondary metabolite gene clusters, we identified genes for inorganic nitrogen and phosphate uptake components and gas vesicles. Our findings contribute to further understanding of the ecologically important genus Planktothrix.

CyanoPATH: a knowledgebase of genome-scale functional repertoire for toxic cyanobacterial blooms

CyanoPATH is a database that curates and analyzes the common genomic functional repertoire for cyanobacteria harmful algal blooms (CyanoHABs) in eutrophic waters. Based on the literature of empirical studies and genome/protein databases, it summarizes four types of information: common biological functions (pathways) driving CyanoHABs, customized pathway maps, classification of blooming type based on databases and the genomes of cyanobacteria. A total of 19 pathways are reconstructed, which are involved in the utilization of macronutrients (e.g. carbon, nitrogen, phosphorus and sulfur), micronutrients (e.g. zinc, magnesium, iron, etc.) and other resources (e.g. light and vitamins) and in stress resistance (e.g. lead and copper). These pathways, comprised of both transport and biochemical reactions, are reconstructed with proteins from NCBI and reactions from KEGG and visualized with self-created transport/reaction maps. The pathways are hierarchical and consist of subpathways, protein/enzyme complexes and constituent proteins. New cyanobacterial genomes can be annotated and visualized for these pathways and compared with existing species. This set of genomic functional repertoire is useful in analyzing aquatic metagenomes and metatranscriptomes in CyanoHAB research. Most importantly, it establishes a link between genome and ecology. All these reference proteins, pathways and maps and genomes are free to download at http://www.csbg-jlu.info/CyanoPATH.