Coherence of Microcystis species revealed through population genomics

Metabologenomics reveals strain-level genetic and chemical diversity of Microcystis secondary metabolism

Microcystis spp. are renowned for producing the hepatotoxin microcystin in freshwater cyanobacterial harmful algal blooms around the world, threatening drinking water supplies and public and environmental health. However, Microcystis genomes also harbor numerous biosynthetic gene clusters (BGCs) encoding the biosynthesis of other secondary metabolites, including many with toxic properties. Most of these BGCs are uncharacterized and currently lack links to biosynthesis products. However, recent field studies show that many of these BGCs are abundant and transcriptionally active in natural communities, suggesting potentially important yet unknown roles in bloom ecology and water quality. Here, we analyzed 21 xenic Microcystis cultures isolated from western Lake Erie to investigate the diversity of the biosynthetic potential of this genus. Through metabologenomic and in silico approaches, we show that these Microcystis strains contain variable BGCs, previously observed in natural populations, and encode distinct metabolomes across cultures. Additionally, we find that the majority of metabolites and gene clusters are uncharacterized, highlighting our limited understanding of the chemical repertoire of Microcystis spp. Due to the complex metabolomes observed in culture, which contain a wealth of diverse congeners as well as unknown metabolites, these results underscore the need to deeply explore and identify secondary metabolites produced by Microcystis beyond microcystins to assess their impacts on human and environmental health.IMPORTANCEThe genus Microcystis forms dense cyanobacterial harmful algal blooms (cyanoHABs) and can produce the toxin microcystin, which has been responsible for drinking water crises around the world. While microcystins are of great concern, Microcystis also produces an abundance of other secondary metabolites that may be of interest due to their potential for toxicity, ecological importance, or pharmaceutical applications. In this study, we combine genomic and metabolomic approaches to study the genes responsible for the biosynthesis of secondary metabolites as well as the chemical diversity of produced metabolites in Microcystis strains from the Western Lake Erie Culture Collection. This unique collection comprises Microcystis strains that were directly isolated from western Lake Erie, which experiences substantial cyanoHAB events annually and has had negative impacts on drinking water, tourism, and industry.

Dissolved organic matter, calcium ion and extracellular polymeric substances on living associated bacteria of Microcystis colony are crucial for unicellular Microcystis to efficiently form colonies

Microcystis typically forms colonies under natural conditions, which contributes to occurrence and prevalence of algal blooms. The colonies consist of Microcystis and associated bacteria (AB), embedded in extracellular polymeric substances (EPS). Previous studies indicate that AB can induce Microcystis to form colonies, however the efficiency is generally low and results in a uniform morphotype. In this study, by using filtrated natural water, several AB strains induced unicellular M. aeruginosa to form colonies resembling several Microcystis morphotypes. The mechanisms were investigated with Methylobacterium sp. Z5. Ca2+ was necessary for Z5 to induce Microcystis to form colonies, while dissolved organic matters (DOM) facilitated AB to agglomerate Microcystis to form large colonies. EPS of living Z5, mainly the aromatic protein components, played a key role in colony induction. Z5 initially aggregated Microcystis via the bridging effects of Ca2+ and DOM, followed by the induction of EPS synthesis and secretion in Microcystis. In this process, the colony forming mode shifted from cell adhesion to a combination of cell adhesion and cell division. Intriguingly, Z5 drove the genomic rearrangement of Microcystis by upregulating some transposase genes. This study unveiled a novel mechanism about Microcystis colony formation and identified a new driver of Microcystis genomic evolution.

Synechococcus dominance induced after hydrogen peroxide treatment of Microcystis bloom in the Caloosahatchee River, Florida

Few field trials examining hydrogen peroxide as a cyanobacterial harmful algal bloom (cHAB) treatment have been conducted in subtropical and tropical regions. None have been tested in Florida, home to Lake Okeechobee and downstream waterways which periodically experience Microcystis bloom events. To investigate treatment effects in Florida, we applied a 490 μM (16.7 mg/L; 0.0015%) hydrogen peroxide spray to a minor bloom of Microcystis aeruginosa on the downstream side of Franklin Lock and Dam in the Caloosahatchee River. Although hydrogen peroxide decreased to background level one day post-treatment, succession was observed in phytoplankton community amplicon sequencing. The relative abundance of Microcystis decreased on day 3 by 86%, whereas the picocyanobacteria Synechococcus became dominant, increasing by 77% on day 3 and by 173% on day 14 to 57% of the phytoplankton community. Metatranscriptomics revealed Synechococcus likely benefitted from the antioxidant defense of upregulated peroxiredoxin, peroxidase/catalase, and rubrerythrin expressions immediately after treatment, and upregulated nitrate transport and urease to take advantage of available nitrogen. Our results indicated hydrogen peroxide induces succession of the phytoplankton community from Microcystis to non-toxic picocyanobacteria and could be used for selective suppression of harmful cyanobacteria.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

Comparative genomic analysis of Microcystis strain diversity using conserved marker genes

Microcystis-dominated cyanobacterial harmful algal blooms (cyanoHABs) have a global impact on freshwater environments, affecting both wildlife and human health. Microcystis diversity and function in field samples and laboratory cultures can be determined by sequencing whole genomes of cultured isolates or natural populations, but these methods remain computationally and financially expensive. Amplicon sequencing of marker genes is a lower cost and higher throughput alternative to characterize strain composition and diversity in mixed samples. However, the selection of appropriate marker gene region(s) and primers requires prior understanding of the relationship between single gene genotype, whole genome content, and phenotype. To identify phylogenetic markers of Microcystis strain diversity, we compared phylogenetic trees built from each of 2,351 individual core genes to an established phylogeny and assessed the ability of these core genes to predict whole genome content and bioactive compound genotypes. We identified single-copy core genes better able to resolve Microcystis phylogenies than previously identified marker genes. We developed primers suitable for current Illumina-based amplicon sequencing with near-complete coverage of available Microcystis genomes and demonstrate that they outperform existing options for assessing Microcystis strain composition. Results showed that genetic markers can be used to infer Microcystis gene content and phenotypes such as potential production of bioactive compounds , although marker performance varies by bioactive compound gene and sequence similarity. Finally, we demonstrate that these markers can be used to characterize the Microcystis strain composition of laboratory or field samples like those collected for surveillance and modeling of Microcystis-dominated cyanobacterial harmful algal blooms.

Genomes of nine biofilm-forming filamentous strains of Cyanobacteria (genera Jaaginema, Scytonema, and Karukerafilum gen. nov.) isolated from mangrove habitats of Guadeloupe (Lesser Antilles)

Biofilm-forming cyanobacteria are abundant in mangrove ecosystems, colonizing various niches including sediment surface and periphyton where they can cover large areas, yet have received limited attention. Several filamentous isolates were recently isolated from Guadeloupe, illustrating the diversity and novelty present in these biofilms. In this study, nine strains belonging to three novel lineages found abundantly in Guadeloupe biofilms are characterized by genome sequencing, morphological and ultrastructural examination, metabolome fingerprinting and searched for secondary metabolites biosynthesis pathways. Assignation of two lineages to known genera is confirmed, namely Scytonema and Jaaginema. The third lineage corresponds to a new Coleofasciculales genus herein described as Karukerafilum gen. nov. The four strains belonging to this genus group into two subclades, one of which displays genes necessary for nitrogen fixation as well as the complete pathway for geosmin production. This study gives new insights into the diversity of mangrove biofilm-forming cyanobacteria, including genome-based description of a new genus and the first genome sequence available for the genus Jaaginema.

Variation in resource competition traits among Microcystis strains is affected by their microbiomes

Freshwater harmful algal blooms are often dominated by Microcystis, a phylogenetically cohesive group of cyanobacteria marked by extensive genetic and physiological diversity. We have previously shown that this genetic diversity and the presence of a microbiome of heterotrophic bacteria influences competitive interactions with eukaryotic phytoplankton. In this study, we sought to explain these observations by characterizing Monod equation parameters for resource usage (maximum growth rate μ max, half-saturation value for growth K s, and quota) as a function of N and P levels for four strains (NIES-843, PCC 9701, PCC 7806 [WT], and PCC 7806 ΔmcyB) in presence and absence of a microbiome derived from Microcystis isolated from Lake Erie. Results indicated limited differences in maximum growth rates but more pronounced differences in half-saturation values among Microcystis strains. The largest impact of the microbiome was reducing the minimal nitrogen concentration sustaining growth and reducing half saturation values, with variable results depending on the Microcystis strain. Microcystis strains also differed from each other in their N and P quotas and the extent to which microbiome presence affected them. Our data highlight the importance of the microbiome in altering Microcystis-intrinsic traits, strain competitive hierarchies, and thus bloom dynamics. As quota, μ max, and K s are commonly used in models for harmful algal blooms, our data suggest that model improvement may be possible by incorporating genotype dependencies of resource-use parameters.

Trade-off in genome turnover events leading to adaptive evolution of Microcystis aeruginosa species complex

BACKGROUND

Numerous studies in the past have expanded our understanding of the genetic differences of global distributed cyanobacteria that originated around billions of years ago, however, unraveling how gene gain and loss drive the genetic evolution of cyanobacterial species, and the trade-off of these evolutionary forces are still the central but poorly understood issues.

RESULTS

To delineate the contribution of gene flow in mediating the hereditary differentiation and shaping the microbial evolution, a global genome-wide study of bloom-forming cyanobacterium, Microcystis aeruginosa species complex, provided robust evidence for genetic diversity, reflected by enormous variation in gene repertoire among various strains. Mathematical extrapolation showed an 'open' microbial pan-genome of M. aeruginosa species, since novel genes were predicted to be introduced after new genomes were sequenced. Identification of numerous horizontal gene transfer's signatures in genome regions of interest suggested that genome expansion via transformation and phage-mediated transduction across bacterial lineage as an evolutionary route may contribute to the differentiation of Microcystis functions (e.g., carbohydrate metabolism, amino acid metabolism, and energy metabolism). Meanwhile, the selective loss of some dispensable genes at the cost of metabolic versatility is as a mean of adaptive evolution that has the potential to increase the biological fitness.

CONCLUSIONS

Now that the recruitment of novel genes was accompanied by a parallel loss of some other ones, a trade-off in gene content may drive the divergent differentiation of M. aeruginosa genomes. Our study provides a genetic framework for the evolution of M. aeruginosa species and illustrates their possible evolutionary patterns.

The Western Lake Erie culture collection: A promising resource for evaluating the physiological and genetic diversity of Microcystis and its associated microbiome

Cyanobacterial harmful algal blooms (cyanoHABs) dominated by Microcystis spp. have significant public health and economic implications in freshwater bodies around the world. These blooms are capable of producing a variety of cyanotoxins, including microcystins, that affect fishing and tourism industries, human and environmental health, and access to drinking water. In this study, we isolated and sequenced the genomes of 21 primarily unialgal Microcystis cultures collected from western Lake Erie between 2017 and 2019. While some cultures isolated in different years have a high degree of genetic similarity (genomic Average Nucleotide Identity >99%), genomic data show that these cultures also represent much of the breadth of known Microcystis diversity in natural populations. Only five isolates contained all the genes required for microcystin biosynthesis while two isolates contained a previously described partial mcy operon. Microcystin production within cultures was also assessed using Enzyme-Linked Immunosorbent Assay (ELISA) and supported genomic results with high concentrations (up to 900 μg L⁻¹) in cultures with complete mcy operons and no or low toxin detected otherwise. These xenic cultures also contained a substantial diversity of bacteria associated with Microcystis, which has become increasingly recognized as an essential component of cyanoHAB community dynamics. These results highlight the genomic diversity among Microcystis strains and associated bacteria in Lake Erie, and their potential impacts on bloom development, toxin production, and toxin degradation. This culture collection significantly increases the availability of environmentally relevant Microcystis strains from temperate North America.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Intra-population genomic diversity of the bloom-forming cyanobacterium, Aphanizomenon gracile, at low spatial scale

Cyanobacteria are oxygenic photosynthetic bacteria that perform a substantial part of the global primary production. Some species are responsible for catastrophic environmental events, called blooms, which have become increasingly common in lakes and freshwater bodies as a consequence of global changes. Genotypic diversity is considered essential for marine cyanobacterial population, allowing it to cope with spatio-temporal environmental variations and to adapt to specific micro-niches in the ecosystem. This aspect is underestimated in the study of bloom development, however, and given little notice in studies of the ecology of harmful cyanobacteria. Here we compared the genomes of four strains of Aphanizomenon gracile, a species of filamentous toxinogenic cyanobacteria (Nostocales) found worldwide in fresh and brackish water. Millimeter-sized fascicles were isolated from a single water sample and have been maintained in culture since 2010. A comparative study revealed extensive heterogeneity in gene contents, despite similar genome size and high similarity indices. These variations were mainly associated with mobile genetic elements and biosynthetic gene clusters. For some of the latter, metabolomic analysis confirmed the production of related secondary metabolites, such as cyanotoxins and carotenoids, which are thought to play a fundamental role in the cyanobacterial fitness. Altogether, these results demonstrated that an A. gracile bloom could be a highly diverse population at low spatial scale and raised questions about potential exchanges of essential metabolites between individuals.

Expression of Microcystis Biosynthetic Gene Clusters in Natural Populations Suggests Temporally Dynamic Synthesis of Novel and Known Secondary Metabolites in Western Lake Erie

Microcystis spp. produce diverse secondary metabolites within freshwater cyanobacterial harmful algal blooms (cyanoHABs) around the world. In addition to the biosynthetic gene clusters (BGCs) encoding known compounds, Microcystis genomes harbor numerous BGCs of unknown function, indicating a poorly understood chemical repertoire. While recent studies show that Microcystis produces several metabolites in the lab and field, little work has focused on analyzing the abundance and expression of its broader suite of BGCs during cyanoHAB events. Here, we use metagenomic and metatranscriptomic approaches to track the relative abundance of Microcystis BGCs and their transcripts throughout the 2014 western Lake Erie cyanoHAB. The results indicate the presence of several transcriptionally active BGCs that are predicted to synthesize both known and novel secondary metabolites. The abundance and expression of these BGCs shifted throughout the bloom, with transcript abundance levels correlating with temperature, nitrate, and phosphorus concentrations and the abundance of co-occurring predatory and competitive eukaryotic microorganisms, suggesting the importance of both abiotic and biotic controls in regulating expression. This work highlights the need for understanding the chemical ecology and potential risks to human and environmental health posed by secondary metabolites that are produced but often unmonitored. It also indicates the prospects for identifying pharmaceutical-like molecules from cyanoHAB-derived BGCs. IMPORTANCE Microcystis spp. dominate cyanobacterial harmful algal blooms (cyanoHABs) worldwide and pose significant threats to water quality through the production of secondary metabolites, many of which are toxic. While the toxicity and biochemistry of microcystins and several other compounds have been studied, the broader suite of secondary metabolites produced by Microcystis remains poorly understood, leaving gaps in our understanding of their impacts on human and ecosystem health. We used community DNA and RNA sequences to track the diversity of genes encoding synthesis of secondary metabolites in natural Microcystis populations and assess patterns of transcription in western Lake Erie cyanoHABs. Our results reveal the presence of both known gene clusters that encode toxic secondary metabolites as well as novel ones that may encode cryptic compounds. This research highlights the need for targeted studies of the secondary metabolite diversity in western Lake Erie, a vital freshwater source to the United States and Canada.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km² region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

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.

Effect of hydrological modification on the potential toxicity of Microcystis aeruginosa complex in Salto Grande reservoir, Uruguay

It is widely known that the environmental conditions caused by the construction of reservoirs favor the proliferation of toxic cyanobacteria and the formation of blooms due to the high residence time of the water, low turbidity, temperature regimes, among others. Microcystin-producing cyanobacteria such as those from the Microcystis aeruginosa complex (MAC) are the most frequently found organisms in reservoirs worldwide, being the role of the environment on microcystin production poorly understood. Here, we addressed the community dynamics and potential toxicity of MAC cyanobacteria in a subtropical reservoir (Salto Grande) located in the low Uruguay river. Samples were taken from five different sites (upstream, inside the reservoir and downstream) during contrasting seasons (summer and winter) to analyze: (i) the MAC community structure by amplicon sequencing of the phycocyanin gene spacer, (ii) the genotype diversity of microcystin-producing MAC by high resolution melting analysis of the mcyJ gene, and (iii) the abundance and mcy transcription activity of the microcystin-producing (toxic) fraction. We found that MAC diversity decreased from summer to winter but, despite the observed changes in MAC community structure, the abundance of toxic organisms and the transcription of mcy genes were always higher inside the reservoir, regardless of the season. Two different genotypes of toxic MAC were detected inside the reservoir, one associated with low water temperature (15 °C) and one thriving at high water temperature (31 °C). These findings indicate that the environmental conditions inside the reservoir reduce community diversity while promoting the proliferation of toxic genotypes that actively transcribe mcy genes, whose relative abundance will depend on the water temperature.

Decoding geobiological evolution from microbiomes

Genomic records of genetic recombination and mutation rates indicate that freshwater ammonia-oxidizing archaea have evolved through paleoclimate and geohydrological history.

Microcystis pangenome reveals cryptic diversity within and across morphospecies

Microcystis, a common harmful algal bloom (HAB) taxon, threatens water supplies and human health, yet species delimitation is contentious in this taxon, leading to challenges in research and management of this threat. Historical and common morphology-based classifications recognize multiple morphospecies, most with variable and diverse ecologies, while DNA sequence-based classifications indicate a single species with multiple ecotypes. To better delimit Microcystis species, we conducted a pangenome analysis of 122 genomes. Core- and non-core gene phylogenetic analyses placed 113 genomes into 23 monophyletic clusters containing at least two genomes. Overall, genome-related indices revealed that Microcystis contains at least 16 putative genospecies. Fifteen genospecies included at least one Microcystis aeruginosa morphospecies, and 10 genospecies included two or more morphospecies. This classification system will enable consistent taxonomic identification of Microcystis and thereby aid in resolving some of the complexities and controversies that have long characterized eco-evolutionary research and management of this important HAB taxon.

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.

High genomic differentiation and limited gene flow indicate recent cryptic speciation within the genus Laspinema (cyanobacteria)

The sympatric occurrence of closely related lineages displaying conserved morphological and ecological traits is often characteristic of free-living microbes. Gene flow, recombination, selection, and mutations govern the genetic variability between these cryptic lineages and drive their differentiation. However, sequencing conservative molecular markers (e.g., 16S rRNA) coupled with insufficient population-level sampling hindered the study of intra-species genetic diversity and speciation in cyanobacteria. We used phylogenomics and a population genomic approach to investigate the extent of local genomic diversity and the mechanisms underlying sympatric speciation of Laspinema thermale. We found two cryptic lineages of Laspinema. The lineages were highly genetically diverse, with recombination occurring more frequently within than between them. That suggests the existence of a barrier to gene flow, which further maintains divergence. Genomic regions of high population differentiation harbored genes associated with possible adaptations to high/low light conditions and stress stimuli, although with a weak diversifying selection. Overall, the diversification of Laspinema species might have been affected by both genomic and ecological processes.

Early and late cyanobacterial bloomers in a shallow, eutrophic lake

Cyanobacterial blooms present challenges for water treatment, especially in regions like the Canadian prairies where poor water quality intensifies water treatment issues. Buoyant cyanobacteria that resist sedimentation present a challenge as water treatment operators attempt to balance pre-treatment and toxic disinfection by-products. Here, we used microscopy to identify and describe the succession of cyanobacterial species in Buffalo Pound Lake, a key drinking water supply. We used indicator species analysis to identify temporal grouping structures throughout two sampling seasons from May to October 2018 and 2019. Our findings highlight two key cyanobacterial bloom phases - a mid-summer diazotrophic bloom of Dolichospermum spp. and an autumn Planktothrix agardhii bloom. Dolichospermum crassa and Woronichinia compacta served as indicators of the mid-summer and autumn bloom phases, respectively. Different cyanobacterial metabolites were associated with the distinct bloom phases in both years: toxic microcystins were associated with the mid-summer Dolichospermum bloom and some newly monitored cyanopeptides (anabaenopeptin A and B) with the autumn Planktothrix bloom. Despite forming a significant proportion of the autumn phytoplankton biomass (>60%), the Planktothrix bloom had previously not been detected by sensor or laboratory-derived chlorophyll-a. Our results demonstrate the power of targeted taxonomic identification of key species as a tool for managers of bloom-prone systems. Moreover, we describe an autumn Planktothrix agardhii bloom that has the potential to disrupt water treatment due to its evasion of detection. Our findings highlight the importance of identifying this autumn bloom given the expectation that warmer temperatures and a longer ice-free season will become the norm.

The Self-Bleaching Process of Microcystis aeruginosa is Delayed by a Symbiotic Bacterium Pseudomonas sp. MAE1-K and Promoted by Methionine Deficiency

Various interactions between marine cyanobacteria and heterotrophic bacteria have been known, but the symbiotic relationships between Microcystis and heterotrophic bacteria remain unclear. An axenic M. aeruginosa culture (NIES-298) was quickly bleached after exponential growth, whereas a xenic M. aeruginosa culture (KW) showed a normal growth curve, suggesting that some symbiotic bacteria may delay this bleaching. The bleaching process of M. aeruginosa was distinguished from the phenomena of previously proposed chlorosis and programmed cell death in various characteristics. Bleached cultures of NIES-298 quickly bleached actively growing M. aeruginosa cultures, suggesting that M. aeruginosa itself produces bleach-causing compounds. Pseudomonas sp. MAE1-K delaying the bleaching of NIES-298 cultures was isolated from the KW culture. Bleached cultures of NIES-298 treated with strain MAE1-K lost their bleaching ability, suggesting that strain MAE1-K rescues M. aeruginosa from bleaching via inactivation of bleaching compounds. From Tn5 transposon mutant screening, a metZ mutant of strain MAE1-K (F-D3) unable to synthesize methionine, promoting the bleaching of NIES-298 cultures but capable of inactivating bleaching compounds, was obtained. The bleaching process of NIES-298 cultures was promoted with the coculture of mutant F-D3 and delayed by methionine supplementation, suggesting that the bleaching process of M. aeruginosa is promoted by methionine deficiency.

IMPORTANCE Cyanobacterial blooms in freshwaters represent serious global concerns for the ecosystem and human health. In this study, we found that one of the major species in cyanobacterial blooms, Microcystis aeruginosa, was quickly collapsed after exponential growth by producing self-bleaching compounds and that a symbiotic bacterium, Pseudomonas sp. MAE1-K delayed the bleaching process via the inactivation of bleaching compounds. In addition, we found that a metZ mutant of strain MAE1-K (F-D3) causing methionine deficiency promoted the bleaching process of M. aeruginosa, suggesting that methionine deficiency may induce the production of bleaching compounds. These results will provide insights into the symbiotic relationships between M. aeruginosa and heterotrophic bacteria that will contribute to developing novel strategies to control cyanobacterial blooms.

Metagenomic and Metatranscriptomic Insights into Population Diversity of Microcystis Blooms: Spatial and Temporal Dynamics of mcy Genotypes, Including a Partial Operon That Can Be Abundant and Expressed

Cyanobacterial harmful algal blooms (cyanoHABs) degrade freshwater ecosystems globally. Microcystis aeruginosa often dominates cyanoHABs and produces microcystin (MC), a class of hepatotoxins that poses threats to human and animal health. Microcystin toxicity is influenced by distinct structural elements across a diversity of related molecules encoded by variant mcy operons. However, the composition and distribution of mcy operon variants in natural blooms remain poorly understood. Here, we characterized the variant composition of mcy genes in western Lake Erie Microcystis blooms from 2014 and 2018. Sampling was conducted across several spatial and temporal scales, including different bloom phases within 2014, extensive spatial coverage on the same day (2018), and frequent, autonomous sampling over a 2-week period (2018). Mapping of metagenomic and metatranscriptomic sequences to reference sequences revealed three Microcystis mcy genotypes: complete (all genes present [mcyA-J]), partial (truncated mcyA, complete mcyBC, and missing mcyD-J), and absent (no mcy genes). We also detected two different variants of mcyB that may influence the production of microcystin congeners. The relative abundance of these genotypes was correlated with pH and nitrate concentrations. Metatranscriptomic analysis revealed that partial operons were, at times, the most abundant genotype and expressed in situ, suggesting the potential biosynthesis of truncated products. Quantification of genetic divergence between genotypes suggests that the observed strains are the result of preexisting heterogeneity rather than de novo mutation during the sampling period. Overall, our results show that natural Microcystis populations contain several cooccurring mcy genotypes that dynamically shift in abundance spatiotemporally via strain succession and likely influence the observed diversity of the produced congeners. IMPORTANCE Cyanobacteria are responsible for producing microcystins (MCs), a class of potent and structurally diverse toxins, in freshwater systems around the world. While microcystins have been studied for over 50 years, the diversity of their chemical forms and how this variation is encoded at the genetic level remain poorly understood, especially within natural populations of cyanobacterial harmful algal blooms (cyanoHABs). Here, we leverage community DNA and RNA sequences to track shifts in mcy genes responsible for producing microcystin, uncovering the relative abundance, expression, and variation of these genes. We studied this phenomenon in western Lake Erie, which suffers annually from cyanoHAB events, with impacts on drinking water, recreation, tourism, and commercial fishing.

Comprehensive and functional analyses reveal the genomic diversity and potential toxicity of Microcystis

Microcystis is a cyanobacteria that is widely distributed across the world. It has attracted great attention because it produces the hepatotoxin microcystin (MC) that can inhibit eukaryotic protein phosphatases and pose a great risk to animal and human health. Due to the high diversity of morphospecies and genomes, it is still difficult to classify Microcystis species. In this study, we investigated the pangenome of 23 Microcystis strains to detect the genetic diversity and evolutionary dynamics. Microcystis revealed an open pangenome containing 22,009 gene families and exhibited different functional constraints. The core-genome phylogenetic analysis accurately differentiated the toxic and nontoxic strains and could be used as a taxonomic standard at the genetic level. We also investigated the functions of HGT events, of which were mostly conferred from cyanobacteria and closely related species. In order to detect the potential toxicity of Microcystis, we searched and characterized MC biosynthetic gene clusters and other secondary metabolite gene clusters. Our work provides insights into the genetic diversity, evolutionary dynamics, and potential toxicity of Microcystis, which could benefit the species classification and development of new methods for drinking water quality control and management of bloom formation in the future.

Amplicon-based and metagenomic approaches provide insights into toxigenic potential in understudied Atlantic Canadian lakes

Cyanobacterial blooms and their toxigenic potential threaten freshwater resources worldwide. In Atlantic Canada, despite an increase of cyanobacterial blooms in the last decade, little is known about the toxigenic potential and the taxonomic affiliation of bloom-forming cyanobacteria. In this study, we employed polymerase chain reaction (PCR) and metagenomic approaches to assess the potential for cyanotoxin and other bioactive metabolite production in Harvey Lake (oligotrophic) and Washademoak Lake (mesotrophic) in New Brunswick, Canada, during summer and early fall months. The PCR survey detected the potential for microcystin (hepatotoxin) and anatoxin-a (neurotoxin) production in both lakes, despite a cyanobacterial bloom only being visible in Washademoak. Genus-specific PCR associated microcystin production potential with the presence of Microcystis in both lakes. The metagenomic strategy provided insight into temporal variations in the microbial communities of both lakes. It also permitted the recovery of a near-complete Microcystis aeruginosa genome with the genetic complement to produce microcystin and other bioactive metabolites such as piricyclamide, micropeptin/cyanopeptolin, and aeruginosin. Our approaches demonstrate the potential for production of a diverse complement of bioactive compounds and establish important baseline data for future studies of understudied lakes, which are frequently affected by cyanobacterial blooms.

Genotyping and multivariate regression trees reveal ecological diversification within the Microcystis aeruginosa complex along a wide environmental gradient

Addressing the ecological and evolutionary processes underlying biodiversity patterns is essential to identify the mechanisms shaping community structure and function. In bacteria, the formation of new ecologically distinct populations (ecotypes) is proposed as one of the main drivers of diversification. New ecotypes arise when mutations in key functional genes or acquisition of new metabolic pathways by horizontal gene transfer allow the population to exploit new resources, permitting their coexistence with the parental population. We previously reported the presence of microcystin-producing organisms of the Microcystis aeruginosa complex (toxic MAC) through an 800 km environmental gradient ranging from freshwater to estuarine-marine waters in South America. We hypothesize that the success of toxic MAC in such a gradient is due to the existence of very closely related populations that are ecologically distinct (ecotypes), each specialized to a specific arrangement of environmental variables. Here, we analyzed toxic MAC genetic diversity through qPCR and high-resolution melting analysis (HRMA) of a functional gene (mcyJ, microcystin synthetase cluster). We explored the variability of the mcyJ gene along the environmental gradient by multivariate classification and regression trees (mCART). Six groups of mcyJ genotypes were distinguished and associated with different combinations of water temperature, conductivity and turbidity. We propose that each mcyJ variant associated to a defined environmental condition is an ecotype (or species) whose relative abundances vary according to their fitness in the local environment. This mechanism would explain the success of toxic MAC in such a wide array of environmental conditions. Importance Organisms of the Microcystis aeruginosa Complex form harmful algal blooms (HABs) in nutrient-rich water bodies worldwide. MAC HABs are difficult to manage owing to the production of potent toxins (microcystins) that resist water treatment. Besides, the role of microcystins in the ecology of MAC organisms is still elusive, meaning that the environmental conditions driving the toxicity of the bloom are not clear. Furthermore, the lack of coherence between morphology-based and genomic-based species classification makes it difficult to draw sound conclusions about when and where each member species of the MAC will dominate the bloom. Here, we propose that the diversification process and success of toxic MAC in a wide range of waterbodies involves the generation of ecotypes, each specialized in a particular niche, whose relative abundance varies according to its fitness in the local environment. This knowledge can improve the generation of accurate prediction models of MAC growth and toxicity, helping to prevent human and animal intoxication.

Rapid Cyanobacteria Species Identification with High Sensitivity Using Native Mass Spectrometry

Cyanobacteria have evolved over billions of years to adapt and survive in diverse climates. Environmentally, this presents a huge challenge because cyanobacteria can now rapidly form algae blooms that are detrimental to aquatic life. In addition, many cyanobacteria produce toxins, making them hazardous to animals and humans that they encounter. Rapid identification of cyanobacteria is essential to monitor and prevent toxic algae blooms. Here, we show for the first time how native mass spectrometry can quickly and precisely identify cyanobacteria from diverse aquatic environments. By monitoring phycobiliproteins, abundant protein complexes within cyanobacteria, simple, easy-to-understand mass spectral "fingerprints" were created that were unique to each species. Moreover, our method is 10-fold more sensitive than the current MALDI-TOF mass spectrometric methods, meaning that cyanobacteria can be monitored using this technology prior to bloom formation. Together, the data show great promise for the simultaneous detection and identification of co-existing cyanobacteria in situ.

Single-colony sequencing reveals microbe-by-microbiome phylosymbiosis between the cyanobacterium Microcystis and its associated bacteria

BACKGROUND

Cyanobacteria from the genus Microcystis can form large mucilaginous colonies with attached heterotrophic bacteria-their microbiome. However, the nature of the relationship between Microcystis and its microbiome remains unclear. Is it a long-term, evolutionarily stable association? Which partners benefit? Here we report the genomic diversity of 109 individual Microcystis colonies-including cyanobacteria and associated bacterial genomes-isolated in situ and without culture from Lake Champlain, Canada and Pampulha Reservoir, Brazil.

RESULTS

We identified 14 distinct Microcystis genotypes from Canada, of which only two have been previously reported, and four genotypes specific to Brazil. Microcystis genetic diversity was much greater between than within colonies, consistent with colony growth by clonal expansion rather than aggregation of Microcystis cells. We also identified 72 bacterial species in the microbiome. Each Microcystis genotype had a distinct microbiome composition, and more closely related genotypes had more similar microbiomes. This pattern of phylosymbiosis could be explained by co-phylogeny in only two out of the nine most prevalent associated bacterial genera, Roseomonas and Rhodobacter. These phylogenetically associated genera could enrich the metabolic repertoire of Microcystis, for example by encoding the biosynthesis of complementary carotenoid molecules. In contrast, other colony-associated bacteria showed weaker signals of co-phylogeny, but stronger evidence of horizontal gene transfer with Microcystis. These observations suggest that acquired genes are more likely to be retained in both partners (Microcystis and members of its microbiome) when they are loosely associated, whereas one gene copy is sufficient when the association is physically tight and evolutionarily long-lasting.

CONCLUSIONS

We have introduced a method for culture-free isolation of single colonies from nature followed by metagenomic sequencing, which could be applied to other types of microbes. Together, our results expand the known genetic diversity of both Microcystis and its microbiome in natural settings, and support their long-term, specific, and potentially beneficial associations. Video Abstract.

Genome Evolution of Filamentous Cyanobacterium Nostoc Species: From Facultative Symbiosis to Free Living

In contrast to obligate bacteria, facultative symbiotic bacteria are mainly characterized by genome enlargement. However, the underlying relationship of this feature with adaptations to various habitats remains unclear. In this study, we used the global genome data of Nostoc strains, including 10 novel genomes sequenced in this study and 26 genomes available from public databases, and analyzed their evolutionary history. The evolutionary boundary of the real clade of Nostoc species was identified and was found to be consistent with the results of polyphasic taxonomy. The initial ancestral species of Nostoc was demonstrated to be consistent with a facultative symbiotic population. Further analyses revealed that Nostoc strains tended to shift from facultative symbiosis to a free-living one, along with an increase in genome sizes during the dispersal of each exterior branch. Intracellular symbiosis was proved to be essentially related to Nostoc evolution, and the adaptation of its members to free-living environments was coupled with a large preference for gene acquisition involved in gene repair and recombination. These findings provided unique evidence of genomic mechanisms by which homologous microbes adapt to distinct life manners and external environments.

The genetic and ecophysiological diversity of Microcystis

Microcystis is a cyanobacterium that forms toxic blooms in freshwater ecosystems around the world. Biological variation among taxa within the genus is apparent through genetic and phenotypic differences between strains and via the spatial and temporal distribution of strains in the environment, and this fine-scale diversity exerts strong influence over bloom toxicity. Yet we do not know how varying traits of Microcystis strains govern their environmental distribution, the tradeoffs and links between these traits, or how they are encoded at the genomic level. Here we synthesize current knowledge on the importance of diversity within Microcystis and on the genes and traits that likely underpin ecological differentiation of taxa. We briefly review spatial and environmental patterns of Microcystis diversity in the field and genetic evidence for cohesive groups within Microcystis. We then compile data on strain-level diversity regarding growth responses to environmental conditions and explore evidence for variation of community interactions across Microcystis strains. Potential links and tradeoffs between traits are identified and discussed. The resulting picture, while incomplete, highlights key knowledge gaps that need to be filled to enable new models for predicting strain-level dynamics, which influence the development, toxicity and cosmopolitan nature of Microcystis blooms.

Molecular responses to inorganic and organic phosphorus sources in the growth and toxin formation of Microcystis aeruginosa

Toxic cyanobacteria bloom is a ubiquitous phenomenon worldwide in eutrophic lakes or reservoirs. Microcystis, is a cosmopolitan genus in cyanobacteria and exists in many different forms. Microcystis aeruginosa (M. aeruginosa) can produce microcystins (MCs) with strong liver toxicity during its growth and decomposition. Phosphorus (P) is a typical growth limiting factor of M. aeruginosa. Though different forms and concentrations of P are common in natural water, the molecular responses in the growth and MCs formation of M. aeruginosa remain unclear. In this study, laboratory experiments were conducted to determine the uptake of P, cell activity, MCs release, and related gene expression under different concentrations of dissolved inorganic phosphorus (DIP) and dissolved organic phosphorus (DOP). We found that the growth of M. aeruginosa was promoted by increasing DIP concentration but coerced under high concentration (0.6 and 1.0 mg P/L) of DOP after P starvation. The growth stress was not related to the alkaline phosphatase activity (APA). Although alkaline phosphatase (AP) could convert DOP into algae absorbable DIP, the growth status of M. aeruginosa mainly depended on the response mechanism of phosphate transporter expression to the extracellular P concentration. High-concentration DIP promoted MCs production in M. aeruginosa, while high-concentration DOP triggered the release of intracellular MCs rather than affecting MCs production. Our study revealed the molecular responses of algal growth and toxin formation under different P sources, and provided a theoretical basis and novel idea for risk management of eutrophic lakes and reservoirs.

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.

Genotype and host microbiome alter competitive interactions between Microcystis aeruginosa and Chlorella sorokiniana

Cyanobacterial harmful algal blooms (cyanoHABs) continue to increase in frequency and magnitude, threatening global freshwater ecosystems and services. In north-temperate lakes cyanobacteria appear in early summer, succeeding green algae as the dominant phytoplankton group, a pattern thought to be mediated by changes in temperature and bioavailable nutrients. To understand additional drivers of this successional pattern our study used reciprocal invasion experiments to examine the competitive interaction between Microcystis aeruginosa, a dominant contributor to cyanoHABs, and the green alga Chlorella sorokiniana. We considered two factors that may impact these interactions: (1) strain variation, with a specific emphasis on the presence or absence of the gene for the hepatotoxin microcystin, and (2) host-associated bacteria. We used toxic M. aeruginosa PCC 7806 (microcystin producing strain), a non-toxic mutant of PCC 7806, non-toxic M. aeruginosa PCC 9701 (non-microcystin producing strain), and C. sorokiniana. Each organism was available free of all bacteria (i.e., axenic) and with a re-introduced defined bacterial community to generate their xenic counterparts. Competitive interactions were assessed with reciprocal invasion experiments between paired xenic and paired axenic populations of C. sorokiniana and one of the two Microcystis strains, each assessed separately. Flow cytometry and random forest models were used to rapidly discriminate and quantify phytoplankton population densities with 99% accuracy. We found that M. aeruginosa PCC 7806, but not strain PCC 9701, could proliferate from low abundance in a steady-state population of C. sorokiniana. Further, the presence of bacteria allowed M. aeruginosa PCC 7806 to grow to a higher population density into an established C. sorokiniana population than when grown axenic. Conversely, when M. aeruginosa was dominant, C. sorokiniana was only able to proliferate from low density into the PCC 9701 strain, and only when axenic. The mutant of PCC 7806 lacking the ability to produce microcystin behaved similarly to the toxic wild-type, implying microcystin is not responsible for the difference in competitive abilities observed between the two wild-type strains. Quantification of microcystins (MCs) when PCC 7806 M. aeruginosa was introduced into the C. sorokiniana culture showed two-fold more MCs per cell when host-associated bacteria were absent compared to present in both species cultures. Our results show that the ability of M. aeruginosa to compete with C. sorokiniana is determined by genomic differences beyond genes involved in microcystin toxin generation and indicate an important role of host-associated bacteria in mediating phytoplankton interspecies interactions. These results expand our understanding of the key drivers of phytoplankton succession and the establishment and persistence of freshwater harmful cyanobacterial blooms.

Comparative genomics analysis of c-di-GMP metabolism and regulation in Microcystis aeruginosa

BACKGROUND

Cyanobacteria are of special concern because they proliferate in eutrophic water bodies worldwide and affect water quality. As an ancient photosynthetic microorganism, cyanobacteria can survive in ecologically diverse habitats because of their capacity to rapidly respond to environmental changes through a web of complex signaling networks, including using second messengers to regulate physiology or metabolism. A ubiquitous second messenger, bis-(3',5')-cyclic-dimeric-guanosine monophosphate (c-di-GMP), has been found to regulate essential behaviors in a few cyanobacteria but not Microcystis, which are the most dominant species in cyanobacterial blooms. In this study, comparative genomics analysis was performed to explore the genomic basis of c-di-GMP signaling in Microcystis aeruginosa.

RESULTS

Proteins involved in c-di-GMP metabolism and regulation, such as diguanylate cyclases, phosphodiesterases, and PilZ-containing proteins, were encoded in M. aeruginosa genomes. However, the number of identified protein domains involved in c-di-GMP signaling was not proportional to the size of M. aeruginosa genomes (4.97 Mb in average). Pan-genome analysis showed that genes involved in c-di-GMP metabolism and regulation are conservative in M. aeruginosa strains. Phylogenetic analysis showed good congruence between the two types of phylogenetic trees based on 31 highly conserved protein-coding genes and sensor domain-coding genes. Propensity for gene loss analysis revealed that most of genes involved in c-di-GMP signaling are stable in M. aeruginosa strains. Moreover, bioinformatics and structure analysis of c-di-GMP signal-related GGDEF and EAL domains revealed that they all possess essential conserved amino acid residues that bind the substrate. In addition, it was also found that all selected M. aeruginosa genomes encode PilZ domain containing proteins.

CONCLUSIONS

Comparative genomics analysis of c-di-GMP metabolism and regulation in M. aeruginosa strains helped elucidating the genetic basis of c-di-GMP signaling pathways in M. aeruginosa. Knowledge of c-di-GMP metabolism and relevant signal regulatory processes in cyanobacteria can enhance our understanding of their adaptability to various environments and bloom-forming mechanism.