Genomic and Genotypic Characterization of Cylindrospermopsis raciborskii: Toward an Intraspecific Phylogenetic Evaluation by Comparative Genomics

Coupling of cylindrospermopsin and pho-harboring Verrucomicrobia supports the formation of Raphidiopsis blooms in low-phosphorus waters

Cylindrospermopsin (CYN) can induce phytoplankton community to secrete alkaline phosphatase (ALP), which is one of the important strategies for the bloom-forming cyanobacterium Raphidiopsis to thrive in extremely low-phosphorus (P) waters. However, how bacterioplankton community, another major contributor to ALPs in waters, couples to Raphidiopsis through CYN, and the role of this coupling in supporting the dominance of Raphidiopsis in nature remain largely unknown. Here, we conducted microcosm experiments to address this knowledge gap, using a combination of differential filtration-based and metagenomics-based methods to identify the sources of ALPs. We found that, compared with algal-derived ALPs, bacteria-derived ALPs exhibited a more pronounced and sensitive response to CYN. This response to CYN was enhanced under low-P conditions. Interestingly, we found that Verrucomicrobia made the largest contribution to the total abundance of pho genes, which encode ALPs. Having high gene abundance of the CYN-sensing PI3K-AKT signaling pathway, Verrucomicrobia's proportion increased with higher concentrations of CYN under low-P conditions, thereby explaining the observed increase in pho gene abundance. Compared with other cyanobacterial genera, Raphidiopsis had a higher abundance of the pst gene. This suggests that Raphidiopsis exhibited a greater capacity to uptake the inorganic P generated by ALPs secreted by other organisms. Overall, our results reveal the mechanism of CYN-induced ALP secretion and its impact on planktonic P-cycling, and provide valuable insights into the role of CYN in supporting the formation of Raphidiopsis blooms.

Phosphonate consumers potentially contributing to methane production in Brazilian soda lakes

Oxic methane production (OMP) has been reported to significantly contribute to methane emissions from oxic surface waters. Demethylation of organic compounds, photosynthesis-associated methane production, and (bacterio)chlorophyll reduction activity are some of the investigated mechanisms as potential OMP sources related to photosynthetic organisms. Recently, cyanobacteria have often been correlated with methane accumulation and emission in freshwater, marine, and saline systems. The Brazilian Pantanal is the world's largest wetland system, with approximately 10,000 shallow lakes, most of which are highly alkaline and saline extreme environments. We initiated this study with an overall investigation using genetic markers, from which we explored metagenomic and limnological data from the Pantanal soda for five potential OMP pathways. Our results showed a strong positive correlation between dissolved methane concentrations and bloom events. Metagenomic data and nutrients, mainly orthophosphate, nitrogen, iron, and methane concentrations, suggest that the organic phosphorous demethylation pathway has the most potential to drive OMP in lakes with blooms. A specialized bacterial community was identified, including the Cyanobacteria Raphidiopsis, although the bloom does not contain the genes to carry out this process. These data showed enough evidence to infer the occurrence of an OMP pathway at Pantanal soda lakes, including the microbial sources and their relation to the 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.

Cyanobacteria: A Promising Source of Antifungal Metabolites

Cyanobacteria are a rich source of secondary metabolites, and they have received a great deal of attention due to their applicability in different industrial sectors. Some of these substances are known for their notorious ability to inhibit fungal growth. Such metabolites are very chemically and biologically diverse. They can belong to different chemical classes, including peptides, fatty acids, alkaloids, polyketides, and macrolides. Moreover, they can also target different cell components. Filamentous cyanobacteria have been the main source of these compounds. This review aims to identify the key features of these antifungal agents, as well as the sources from which they are obtained, their major targets, and the environmental factors involved when they are being produced. For the preparation of this work, a total of 642 documents dating from 1980 to 2022 were consulted, including patents, original research, review articles, and theses.

An updated classification of cyanobacterial orders and families based on phylogenomic and polyphasic analysis

Cyanobacterial taxonomy is facing a period of rapid changes thanks to the ease of 16S rRNA gene sequencing and established workflows for description of new taxa. Since the last comprehensive review of the cyanobacterial system in 2014 until 2021, at least 273 species in 140 genera were newly described. These taxa were mainly placed into previously defined orders and families although several new families were proposed. However, the classification of most taxa still relied on hierarchical relationships inherited from the classical morphological taxonomy. Similarly, the obviously polyphyletic orders such as Synechococcales and Oscillatoriales were left unchanged. In this study, the rising number of genomic sequences of cyanobacteria and well-described reference strains allowed us to reconstruct a robust phylogenomic tree for taxonomic purposes. A less robust but better sampled 16S rRNA gene phylogeny was mapped to the phylogenomic backbone. Based on both these phylogenies, a polyphasic classification throughout the whole phylum of Cyanobacteria was created, with ten new orders and fifteen new families. The proposed system of cyanobacterial orders and families relied on a phylogenomic tree but still employed phenotypic apomorphies where possible to make it useful for professionals in the field. It was, however, confirmed that morphological convergence of phylogenetically distant taxa was a frequent phenomenon in cyanobacteria. Moreover, the limited phylogenetic informativeness of the 16S rRNA gene, resulting in ambiguous phylogenies above the genus level, emphasized the integration of genomic data as a prerequisite for the conclusive taxonomic placement of a vast number of cyanobacterial genera in the future.

Recent Advances in the Ecology of Bloom-Forming Raphidiopsis (Cylindrospermopsis) raciborskii: Expansion in China, Intraspecific Heterogeneity and Critical Factors for Invasion

Water blooms caused by the invasive cyanobacterium Raphidiopsis raciborskii occur in many reservoirs in the tropical and subtropical regions of China. In recent decades, this species has spread rapidly to temperate regions. Phenotypic plasticity and climate warming are thought to promote the worldwide dispersion of R. raciborskii. However, investigations into the genetic and phenotypic diversities of this species have revealed significant intraspecific heterogeneity. In particular, competition between R. raciborskii and Microcystis aeruginosa was highly strain dependent. Although the concept of an ecotype was proposed to explain the heterogeneity of R. raciborskii strains with different geographic origins, microevolution is more reasonable for understanding the coexistence of different phenotypes and genotypes in the same environment. It has been suggested that intraspecific heterogeneity derived from microevolution is a strong driving force for the expansion of R. raciborskii. Additionally, temperature, nutrient fluctuations, and grazer disturbance are critical environmental factors that affect the population establishment of R. raciborskii in new environments. The present review provides new insights into the ecological mechanisms underlying the invasion of R. raciborskii in Chinese freshwater ecosystems.

Microbial one‑carbon and nitrogen metabolisms are beneficial to the reservoir recovery after cyanobacterial bloom

Cyanobacterial blooms have profound effects on the structure and function of plankton communities in inland waters, but few studies have focused on the effects of microbial-based processes in one‑carbon and nitrogen cycling on water quality improvement following the bloom. Here, we compared the structure and function of the bacterial community, focusing on microbial one‑carbon and nitrogen metabolisms during and after a cyanobacterial Microcystis bloom in a deep subtropical reservoir. Our data showed that microbial one‑carbon and nitrogen cycles were closely related to different periods of the bloom, and the changes of functional genes in microbial carbon and nitrogen cycling showed the same consistent trend as that of Methylomonas sp. With the receding of the bloom, the abundance of Methylomonas as well as the functional genes of microbial one‑carbon and nitrogen cycling reached the peak and then recovered. Our results indicate that microbial one‑carbon and nitrogen metabolisms were beneficial to the recovery of water quality from the cyanobacterial bloom. This study lays a foundation for a deep understanding of the cyanobacterial decomposition mediated by microbes in one‑carbon and nitrogen cycles in inland freshwaters.

Heterocyte production, gene expression and phylogeography in Raphidiopsis ( = Cylindrospermopsis) Raciborskii

Raphidiopsis ( = Cylindrospermopsis) raciborskii was described as a subtropical-tropical cyanobacterium, later reported expanding into temperate regions. Heterocyte presence used to distinguish Cylindrospermopsis from the very similar Raphidiopsis, but recently the two genera were recognized as one and unified. This study aimed to investigate how heterocyte production is related to nitrogen (N) limitation in heterocytous and non-heterocytous strains of R.raciborskii. High N-concentrations did not inhibit heterocyte development in some strains, while prolonged N-starvation periods never stimulated production in others. RT-qPCR was used to examine the genetic background, through the expression patterns of nifH, ntcA and hetR. While gene expression increased under N-restriction, N-sufficiency did not suppress nifH transcripts as previously observed in other diazotrophyc cyanobacteria, suggesting that heterocyte production in R. raciborskii is not regulated by N-availability. Heterocytous and non-heterocytous strains were genotypically characterized to assess their phylogenetic relationships,. In the phylogenetic tree, clusters were intermixed and confirmed Raphidiopsis and Cylindrospermopsis as the same genus. The tree supported previous findings of earlier splitting of American strains, while contesting the African origin hypothesis. The existence of two lines of Chinese strains, with distinct evolutionary patterns, is a significant addition that could lead to new hypotheses of the species biogeography.

Phosphorus deficiency stimulates dominance of Cylindrospermopsis through facilitating cylindrospermopsin-induced alkaline phosphatase secretion: Integrating field and laboratory-based evidences

Potentially toxic Cylindrospermopsis raciborskii blooms are of emerging concerns, as its scale is spreading from tropical regions to high latitudes, increasing the risk of aquatic biota being exposed to cylindrospermopsin (CYN). So far, CYN-producing C. raciborskii strains have only been reported in tropical waters which are commonly phosphorus (P)-deficient, where they can dominate phytoplankton communities. However, the influence of CYN on phytoplankton communities under different P status remains unclear. In this study, we first analyzed the summer observations of 120 tropical reservoirs in Guangdong Province. The proportion of potential CYN-producers was significantly higher in P-deficient and CYN-present reservoirs than that in P-sufficient or CYN-absent ones. This suggested that in P-deficient condition, the potential CYN producers might gain more advantages by the help of CYN. Then, in laboratory experiments we found that upon P deprivation, CYN did not inhibit the cell growth of other algal cells, but significantly stimulates them to secret more alkaline phosphatase (ALP) than in P-sufficient condition. Through transcriptomics, we further revealed that under such P-deficient condition, CYN remarkably induced intracellular nitrogen allocation and protein export system by activating the PIK3/Akt-cGMP/PKG signaling pathways in Scenedesmus bijugatus, thus enhancing its ALP secretion. Our study implies that CYN-induced ALP secretion is facilitated upon P deficiency, thus supporting the dominance of its producers C. raciborskii.

Genome Streamlining, Plasticity, and Metabolic Versatility Distinguish Co-occurring Toxic and Nontoxic Cyanobacterial Strains of Microcoleus

Harmful cyanobacterial bloom occurrences have increased worldwide due to climate change and eutrophication, causing nuisance and animal deaths. Species from the benthic cyanobacterial genus Microcoleus are ubiquitous and form thick mats in freshwater systems, such as rivers, that are sometimes toxic due to the production of potent neurotoxins (anatoxins). Anatoxin-producing (toxic) strains typically coexist with non-anatoxin-producing (nontoxic) strains in mats, although the reason for this is unclear. To determine the genetic mechanisms differentiating toxic and nontoxic Microcoleus, we sequenced and assembled genomes from 11 cultures and compared these to another 31 Microcoleus genomes. Average nucleotide identities (ANI) indicate that toxic and nontoxic strains are distinct species (ANI, <95%), and only 6% of genes are shared across all 42 genomes, suggesting a high level of genetic divergence among Microcoleus strains. Comparative genomics showed substantial genome streamlining in toxic strains and a potential dependency on external sources for thiamine and sucrose. Toxic and nontoxic strains are further differentiated by an additional set of putative nitrate transporter (nitrogen uptake) and cyanophycin (carbon and nitrogen storage) genes, respectively. These genes likely confer distinct competitive advantages based on nutrient availability and suggest nontoxic strains are more robust to nutrient fluctuations. Nontoxic strains also possess twice as many transposable elements, potentially facilitating greater genetic adaptation to environmental changes. Our results offer insights into the divergent evolution of Microcoleus strains and the potential for cooperative and competitive interactions that contribute to the co-occurrence of toxic and nontoxic species within mats. IMPORTANCE Microcoleus autumnalis, and closely related Microcoleus species, compose a geographically widespread group of freshwater benthic cyanobacteria. Canine deaths due to anatoxin-a poisoning, following exposure to toxic proliferations, have been reported globally. While Microcoleus proliferations are on the rise, the mechanisms underpinning competition between, or coexistence of, toxic and nontoxic strains are unknown. This study identifies substantial genetic differences between anatoxin-producing and non-anatoxin-producing strains, pointing to reduced metabolic flexibility in toxic strains, and potential dependence on cohabiting nontoxic strains. Results provide insights into the metabolic and evolutionary differences between toxic and nontoxic Microcoleus, which may assist in predicting and managing aquatic proliferations.

Ecophysiological Aspects and sxt Genes Expression Underlying Induced Chemical Defense in STX-Producing Raphidiopsis raciborskii (Cyanobacteria) against the Zooplankter Daphnia gessneri

Cyanobacteria stand out among phytoplankton when they form massive blooms and produce toxins. Because cyanotoxin genes date to the origin of metazoans, the hypothesis that cyanotoxins function as a defense against herbivory is still debated. Although their primary cellular function might vary, these metabolites could have evolved as an anti-predator response. Here we evaluated the physiological and molecular responses of a saxitoxin-producing Raphidiopsis raciborskii to infochemicals released by the grazer Daphnia gessneri. Induced chemical defenses were evidenced in R. raciborskii as a significant increase in the transcription level of sxt genes, followed by an increase in saxitoxin content when exposed to predator cues. Moreover, cyanobacterial growth decreased, and no significant effects on photosynthesis or morphology were observed. Overall, the induced defense response was accompanied by a trade-off between toxin production and growth. These results shed light on the mechanisms underlying zooplankton-cyanobacteria interactions in aquatic food webs. The widespread occurrence of the cyanobacterium R. raciborskii in freshwater bodies has been attributed to its phenotypic plasticity. Assessing the potential of this species to thrive over interaction filters such as zooplankton grazing pressure can enhance our understanding of its adaptive success.

Genetic Identity and Herbivory Drive the Invasion of a Common Aquatic Microbial Invader

Despite the increasing number of species invasions, the factors driving invasiveness are still under debate. This is particularly the case for “invisible” invasions by aquatic microbial species. Since in many cases only a few individuals or propagules enter a new habitat, their genetic variation is low and might limit their invasion success, known as the genetic bottleneck. Thus, a key question is, how genetic identity and diversity of invading species influences their invasion success and, subsequently, affect the resident community. We conducted invader-addition experiments using genetically different strains of the globally invasive, aquatic cyanobacterium Raphidiopsis raciborskii (formerly: Cylindrospermopsis raciborskii) to determine the role of invader identity and genetic diversity (strain richness) at four levels of herbivory. We tested the invasion success of solitary single strain invasions against the invader genetic diversity, which was experimentally increased up to ten strains (multi-strain populations). By using amplicon sequencing we determined the strain-specific invasion success in the multi-strain treatments and compared those with the success of these strains in the single-strain treatments. Furthermore, we tested for the invasion success under different herbivore pressures. We showed that high grazing pressure by a generalist herbivore prevented invasion, whereas a specialist herbivore enabled coexistence of consumer and invader. We found a weak effect of diversity on invasion success only under highly competitive conditions. When invasions were successful, the magnitude of this success was strain-specific and consistent among invasions performed with single-strain or multi-strain populations. A strain-specific effect was also observed on the resident phytoplankton community composition, highlighting the strong role of invader genetic identity. Our results point to a strong effect of the genetic identity on the invasion success under low predation pressure. The genetic diversity of the invader population, however, had little effect on invasion success in our study, in contrast to most previous findings. Instead, it is the interaction between the consumer abundance and type together with the strain identity of the invader that defined invasion success. This study underlines the importance of strain choice in invasion research and in ecological studies in general.

Intra-population strain variation in phosphorus storage strategies of the freshwater cyanobacterium Raphidiopsis raciborskii

Several cyanobacteria, including diazotrophic Raphidiopsis raciborskii, can form harmful blooms when dissolved inorganic phosphorus concentrations are very low. We hypothesized that R. raciborskii strains would vary in phosphorus (P) allocations to cell growth and storage, providing resilience of populations to continuously low or variable P supplies. We tested this hypothesis using six toxic strains (producing cylindrospermopsins) isolated from a field population using batch monocultures with and without P and dissolved inorganic nitrogen (DIN). Treatments replete with DIN, irrespective of P addition, had similar exponential growth rates for individual strains. P storage capacity varied 4-fold among strains and was significantly higher in DIN-free treatments than in replete treatments. P was stored by all R. raciborskii strains, in preference to allocation to increase growth rates. P stores decreased with increased growth rate across strains, but weeere not related to the time to P starvation in P-free treatments. The storage capacity of R. raciborskii, combined with strategies to efficiently uptake P, means that P controls may not control R. raciborskii populations in the short term. Intra-population strain variation in P storage capacity will need to be reflected in process-based models to predict blooms of R. raciborskii and other cyanobacteria adapted to low-P conditions.

Biogeography of the cyanobacterium Raphidiopsis (Cylindrospermopsis) raciborskii: Integrating genomics, phylogenetic and toxicity data

Raphidiopsis (Cylindrospermopsis) raciborskii, a globally distributed bloom-forming cyanobacterium, produces either the cytotoxin cylindrospermopsin (CYL) in Oceania, Asia and Europe or the neurotoxin saxitoxin (STX) and analogues (paralytic shellfish poison, PSP) in South America (encoded by sxt genetic cluster) and none of them in Africa. Nevertheless, this particular geographic pattern is usually overlooked in current hypotheses about the species dispersal routes. Here, we combined genomics, phylogenetic analyses, toxicity data and a literature survey to unveil the evolutionary history and spread of the species. Phylogenies based on 354 orthologous genes from all the available genomes and ribosomal ITS sequences of the taxon showed two well-defined clades: the American, having the PSP producers; and the Oceania/Europe/Asia, including the CYL producers. We propose central Africa as the original dispersion center (non-toxic populations), reaching North Africa and North America (in former Laurasia continent). The ability to produce CYL probably took place in populations that advanced to sub-Saharan Africa and then to Oceania and South America. According to the genomic context of the sxt cluster found in PSP-producer strains, this trait was acquired once by horizontal transfer in South America, where the ability to produce CYL was lost.

Phylogenomic Analysis of Secondary Metabolism in the Toxic Cyanobacterial Genera Anabaena, Dolichospermum and Aphanizomenon

Cyanobacteria produce an array of toxins that pose serious health risks to humans and animals. The closely related diazotrophic genera, Anabaena, Dolichospermum, and Aphanizomenon, frequently form poisonous blooms in lakes and brackish waters around the world. These genera form a complex now termed the Anabaena, Dolichospermum, and Aphanizomenon (ADA) clade and produce a greater array of toxins than any other cyanobacteria group. However, taxonomic confusion masks the distribution of toxin biosynthetic pathways in cyanobacteria. Here we obtained 11 new draft genomes to improve the understanding of toxin production in these genera. Comparison of secondary metabolite pathways in all available 31 genomes for these three genera suggests that the ability to produce microcystin, anatoxin-a, and saxitoxin is associated with specific subgroups. Each toxin gene cluster was concentrated or even limited to a certain subgroup within the ADA clade. Our results indicate that members of the ADA clade encode a variety of secondary metabolites following the phylogenetic clustering of constituent species. The newly sequenced members of the ADA clade show that phylogenetic separation of planktonic Dolichospermum and benthic Anabaena is not complete. This underscores the importance of taxonomic revision of Anabaena, Dolichospermum, and Aphanizomenon genera to reflect current phylogenomic understanding.

Phylogenetic relationships and genetic divergence of paralytic shellfish toxin- and cylindrospermopsin- producing Cylindrospermopsis and raphidiopsis

Cylindrospermopsis and Raphidiopsis (C/R group) are closely related species responsible for cyanobacterial blooms worldwide. Paralytic shellfish toxins (PSTs) and cylindrospermopsins (CYNs) have been identified in different C/R group strains. However, the evolutionary relationship between PST- and CYN-producing strains has not been systematically evaluated. In this study, C/R group strains and their toxin biosynthesis genes were evaluated by phylogenetic analysis and sequence comparison. None of the tested strains are able to produce PSTs and CYNs simultaneously. The C/R group strains were clustered into five clades, including two non-toxic, two CYN-producing and one PST-producing clades. A high degree of similarity was observed for rpoC1 (> 96%) and ITS-L (> 97%) sequences within each clade with the exception of the ITS-L (87% to 100%) region in CYN-producing R. curvata, which has been shown to contain variable sequence insertions. Genomic analysis revealed that sxtY and sxtZ could be found in both toxic and non-toxic strains. The transposase gene IS4 was only observed in strains from the PST-producing clade. The sxt and cyr gene clusters share five gene families with similar functions. The amino acid sequences of the adenylyl-sulfate kinase genes, sxtO and cyrN, are more similar (45% to 81%) than other pairs of genes (8.0% to 40%). SxtO and CyrN proteins from C/R group strains forms an independent clade on the phylogenetic tree with a high degree of sequence similarity (78% to 100%). In conclusion, PST- and CYN- producing C/R group species can be classified into different clades based on their phylogenetic profile. The sxtO and cyrN genes have probably diverged from a single ancestral adenylyl-sulfate kinase gene, and may be specifically used for toxin biosynthesis in C/R group species.

Intraspecific variability in response to phosphorus depleted conditions in the cyanobacteria Microcystis aeruginosa and Raphidiopsis raciborskii

Phosphorus loading plays an important role in the occurrence of cyanobacterial blooms and understanding how this nutrient affects the physiology of cyanobacteria is imperative to manage these phenomena. Microcystis aeruginosa and Raphidiopsis raciborskii are cyanobacterial species that form potentially toxic blooms in freshwater ecosystems worldwide. Blooms comprise numerous strains with high trait variability, which can contribute to the widespread distribution of these species. Here, we explored the intraspecific variability in response to phosphorus depleted conditions (P-) testing five strains of each species. Strains could be differentiated by cell volume or genetic profiles except for those of the same species, sampling location and date, though these presented differences in their response to (P-). Although differently affected by (P-) over 10 days, all strains were able to grow and maintain photosynthetic activity. For most M. aeruginosa and R. raciborskii strains growth rates were not significantly different comparing (P+) and (P-) conditions. After ten days in (P-), only one M. aeruginosa strain and two R. raciborskii strains showed reduction in biovolume yield as compared to (P+) but in most strains chlorophyll-a concentrations were lower in (P-) than in (P+). Reduced photosystem II efficiency was found for only one R. raciborskii strain while all M. aeruginosa strains were affected. Only two M. aeruginosa and one R. raciborskii strain increased alkaline phosphatase activity under (P-) as compared to (P+). Variation in P-uptake was also observed but comparison among strains yielded homogeneous groups comprised of representatives of both species. Comparing the response of each species as a whole, the (P-) condition affected growth rate, biovolume yield and chlorophyll yield. However, these parameters revealed variation among strains of the same species to the extent that differences between M. aeruginosa and R. raciborskii were not significant. Taken together, these results do not support the idea that R. raciborskii, as a species, can withstand phosphorus limitation better than M. aeruginosa and also point that the level of intraspecific variation may preclude generalizations based on studies that use only one or few strains.

A closely-related clade of globally distributed bloom-forming cyanobacteria within the Nostocales

In order to better understand the relationships among current Nostocales cyanobacterial blooms, eight genomes were sequenced from cultured isolates or from environmental metagenomes of recent planktonic Nostocales blooms. Phylogenomic analysis of publicly available sequences placed the new genomes among a group of 15 genomes from four continents in a distinct ADA clade (Anabaena/Dolichospermum/Aphanizomenon) within the Nostocales. This clade contains four species-level groups, two of which include members with both Anabaena-like and Aphanizomenon flos-aquae-like morphology. The genomes contain many repetitive genetic elements and a sizable pangenome, in which ABC-type transporters are highly represented. Alongside common core genes for photosynthesis, the differentiation of N₂-fixing heterocysts, and the uptake and incorporation of the major nutrients P, N and S, we identified several gene pathways in the pangenome that may contribute to niche partitioning. Genes for problematic secondary metabolites-cyanotoxins and taste-and-odor compounds-were sporadically present, as were other polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) gene clusters. By contrast, genes predicted to encode the ribosomally generated bacteriocin peptides were found in all genomes.