Nodularia spumigena is a nitrogen-fixing cyanobacterium that forms toxic blooms in the Baltic Sea each summer and the availability of phosphorous is an important factor limiting the formation of these blooms. Bioinformatic analysis identified a phosphonate degrading (phn) gene cluster in the genome of N. spumigena suggesting that this bacterium may use phosphonates as a phosphorus source. Our results show that strains of N. spumigena could grow in medium containing methylphosphonic acid (MPn) as the sole source of phosphorous and released methane when growing in medium containing MPn. We analyzed the total transcriptomes of N. spumigena UHCC 0039 grown using MPn and compared them with cultures growing in Pi-replete medium. The phnJ, phosphonate lyase gene, was upregulated when MPn was the sole source of phosphorus, suggesting that the expression of this gene could be used to indicate the presence of bioavailable phosphonates. Otherwise, growth on MPn resulted in only a minor reconstruction of the transcriptome and enabled good growth. However, N. spumigena strains were not able to utilize any of the anthropogenic phosphonates tested. The phosphonate utilizing pathway may offer N. spumigena a competitive advantage in the Pi-limited cyanobacterial blooms of the Baltic Sea.

cInorganic phosphorus (P i ) is one of the main growth-limiting factors of diazotrophic cyanobacteria. Due to human activity, the availability of P i has increased in water bodies, resulting in eutrophication and the formation of massive cyanobacterial blooms. In this study, we examined the molecular responses of the cyanobacterium Anabaena sp. strain 90 to phosphorus deprivation, aiming at the identification of candidate genes to monitor the P i status in cyanobacteria. Furthermore, this study increased the basic understanding of how phosphorus affects diazotrophic and bloom-forming cyanobacteria as a major growth-limiting factor. Based on RNA sequencing data, we identified 246 differentially expressed genes after phosphorus starvation and 823 differentially expressed genes after prolonged P i limitation, most of them related to central metabolism and cellular growth. The transcripts of the genes related to phosphorus transport and assimilation (pho regulon) were most upregulated during phosphorus depletion. One of the most increased transcripts encodes a giant protein of 1,869 amino acid residues, which contains, among others, a phytase-like domain. Our findings predict its crucial role in phosphorus starvation, but future studies are still needed. Using twodimensional difference in gel electrophoresis (2D-DIGE) and liquid chromatography-tandem mass spectrometry (LC-MS/MS), we found 43 proteins that were differentially expressed after prolonged phosphorus stress. However, correlation analysis unraveled an association only to some extent between the transcriptomic and proteomic abundances. Based on the present results, we suggest that the method used for monitoring the P i status in cyanobacterial bloom should contain wider combinations of pho regulon genes (e.g., PstABCS transport systems) in addition to the commonly used alkaline phosphatase gene alone. N itrogen and phosphorus loading into aquatic ecosystems is a global problem causing enhanced primary production of photoautotrophic organisms (1-3). Among them, cyanobacteria are major beneficiaries of turbid and eutrophic water. They employ numerous ecophysiological strategies such as the capability for using various nutrient sources, gas vacuoles, photoprotective pigments, and UV-absorbing compounds, which enable them to survive and grow under eutrophic conditions and may lead to the formation of harmful and toxic blooms (1, 4). The theory that nutrient loading is a crucial factor causing harmful cyanobacterial blooms is supported by numerous field and laboratory studies showing that toxin production and growth rate both increase under nutrient-replete conditions (2, 5-8).The diazotrophic Anabaena is one of the most common bloomforming cyanobacterial genera in freshwater and brackish-water ecosystems, alongside Microcystis and Nodularia, respectively (4, 9), and some strains are capable of producing the following toxins: microcystins, anatoxin-a, anatoxin-a(S), and saxitoxins (4, 10). The ecology of these cyanobacteria thus deserves more attention. The abil...

Salinity is an important abiotic factor controlling the distribution and abundance of Nodularia spumigena, the dominating diazotrophic and toxic phototroph, in the brackish water cyanobacterial blooms of the Baltic Sea. To expand the available genomic information for brackish water cyanobacteria, we sequenced the isolate Nodularia spumigena UHCC 0039 using an Illumina-SMRT hybrid sequencing approach, revealing a chromosome of 5,294,286 base pairs (bp) and a single plasmid of 92,326 bp. Comparative genomics in Nostocales showed pronounced genetic similarity among Nodularia spumigena strains evidencing their short evolutionary history. The studied Baltic Sea strains share similar sets of CRISPR-Cas cassettes and a higher number of insertion sequence (IS) elements compared to Nodularia spumigena CENA596 isolated from a shrimp production pond in Brazil. Nodularia spumigena UHCC 0039 proliferated similarly at three tested salinities, whereas the lack of salt inhibited its growth and triggered transcriptome remodeling, including the up-regulation of five sigma factors and the down-regulation of two other sigma factors, one of which is specific for strain UHCC 0039. Down-regulated genes additionally included a large genetic region for the synthesis of two yet unidentified natural products. Our results indicate a remarkable plasticity of the Nodularia salinity acclimation, and thus salinity strongly impacts the intensity and distribution of cyanobacterial blooms in the Baltic Sea.

The Baltic Sea is a shallow basin of brackish water in which the spatial salinity gradient is one of the most important factors contributing to species distribution. The Baltic Sea is infamous for its annual cyanobacterial blooms comprised of Nodularia spumigena, Aphanizomenon spp., and Dolichospermum spp. that cause harm, especially for recreational users. To broaden our knowledge of the cyanobacterial adaptation strategies for brackish water environments, we sequenced the entire genome of Dolichospermum sp. UHCC 0315, a species occurring not only in freshwater environments but also in brackish water. Comparative genomics analyses revealed a close association with Dolichospermum sp. UHCC 0090 isolated from a lake in Finland. The genome closure of Dolichospermum sp. UHCC 0315 unraveled a mixture of two subtypes in the original culture, and subtypes exhibited distinct buoyancy phenotypes. Salinity less than 3 g L−1 NaCl enabled proper growth of Dolichospermum sp. UHCC 0315, whereas growth was arrested at moderate salinity (6 g L−1 NaCl). The concentrations of toxins, microcystins, increased at moderate salinity, whereas RNA sequencing data implied that Dolichospermum remodeled its primary metabolism in unfavorable high salinity. Based on our results, the predicted salinity decrease in the Baltic Sea may favor toxic blooms of Dolichospermum spp.

Bacteria are major drivers of biogeochemical nutrient cycles and energy fluxes in marine environments, yet how bacterial communities respond to environmental change is not well known. Metagenomes allow examination of genetic responses of the entire microbial community to environmental change. However, it is challenging to link metagenomes directly to biogeochemical process rates. Here, we investigate metagenomic responses in natural bacterioplankton communities to simulated environmental stressors in the Baltic Sea, including increased river water input, increased nutrient concentration, and reduced oxygen level. This allowed us to identify informative prokaryotic gene markers, responding to environmental perturbation. Our results demonstrate that metagenomic and metabolic changes in bacterial communities in response to environmental stressors are influenced both by the initial community composition and by the biogeochemical factors shaping the functional response. Furthermore, the different sources of dissolved organic matter (DOM) had the largest impact on metagenomic blueprint. Most prominently, changes in DOM loads influenced specific transporter types reflecting the substrate availability and DOC assimilation and consumption pathways. The results provide new knowledge for developing models of ecosystem structure and biogeochemical cycling in future climate change scenarios and advance our exploration of the potential use of marine microorganisms as markers for environmental conditions.