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.

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.

Cyanobacterial biodiversity of semiarid public drinking water supply reservoirs assessed via next-generation DNA sequencing technology.

Next-generation DNA sequencing technology was applied to generate molecular data from semiarid reservoirs during well-defined seasons. Target sequences of 16S-23S rRNA ITS and cpcBA-IGS were used to reveal the taxonomic groups of cyanobacteria present in the samples, and genes coding for cyanotoxins such as microcystins (mcyE), saxitoxins (sxtA), and cylindrospermopsins (cyrJ) were investigated. The presence of saxitoxins in the environmental samples was evaluated using ELISA kit. Taxonomic analyses of high-throughput DNA sequencing data showed the dominance of the genus Microcystis in Mundaú reservoir. Furthermore, it was the most abundant genus in the dry season in Ingazeira reservoir. In the rainy season, 16S-23S rRNA ITS analysis revealed that Cylindrospermopsis raciborskii comprised 46.8% of the cyanobacterial community in Ingazeira reservoir, while the cpcBAIGS region revealed that C. raciborskii (31.8%) was the most abundant taxon followed by Sphaerospermopsis aphanizomenoides (17.3%) and Planktothrix zahidii (16.6%). Despite the presence of other potential toxin-producing genera, the detected sxtA gene belonged to C. raciborskii, while the mcyE gene belonged to Microcystis in both reservoirs. The detected mcyE gene had good correlation with MC content, while the amplification of the sxtA gene was related to the presence of STX. The cyrJ gene was not detected in these samples. Using DNA analyses, our results showed that the cyanobacterial composition of Mundaú reservoir was similar in successive dry seasons, and it varied between seasons in Ingazeira reservoir. In addition, our data suggest that some biases of analysis influenced the cyanobacterial communities seen in the NGS output of Ingazeira reservoir.

Is qPCR a Reliable Indicator of Cyanotoxin Risk in Freshwater?

The wide distribution of cyanobacteria in aquatic environments leads to the risk of water contamination by cyanotoxins, which generate environmental and public health issues. Measurements of cell densities or pigment contents allow both the early detection of cellular growth and bloom monitoring, but these methods are not sufficiently accurate to predict actual cyanobacterial risk. To quantify cyanotoxins, analytical methods are considered the gold standards, but they are laborious, expensive, time-consuming and available in a limited number of laboratories. In cyanobacterial species with toxic potential, cyanotoxin production is restricted to some strains, and blooms can contain varying proportions of both toxic and non-toxic cells, which are morphologically indistinguishable. The sequencing of cyanobacterial genomes led to the description of gene clusters responsible for cyanotoxin production, which paved the way for the use of these genes as targets for PCR and then quantitative PCR (qPCR). Thus, the quantification of cyanotoxin genes appeared as a new method for estimating the potential toxicity of blooms. This raises a question concerning whether qPCR-based methods would be a reliable indicator of toxin concentration in the environment. Here, we review studies that report the parallel detection of microcystin genes and microcystin concentrations in natural populations and also a smaller number of studies dedicated to cylindrospermopsin and saxitoxin. We discuss the possible issues associated with the contradictory findings reported to date, present methodological limitations and consider the use of qPCR as an indicator of cyanotoxin risk.

Influence of nitrogen availability on the expression of genes involved in the biosynthesis of saxitoxin and analogs in Cylindrospermopsis raciborskii.

The development of cyanobacterial blooms in inland aquatic ecosystems is greatly promoted by nutrient availability, especially nitrogen and phosphorous. When blooms are dominated by toxigenic species the harmful effects of nutrient loading becomes particularly relevant. Among toxic species, Cylindrospermopsis raciborskii found in South American ecosystems is characterized by the production of saxitoxin and analogs (Paralytic Shellfish Poisoning, PSP), for which the factors that trigger their production have not been elucidated. In this study, the effect of nitrate availability on the relative transcript abundance of two genes (sxtU and sxtI), both involved in different steps of PSP biosynthetic pathway, was addressed in C. raciborskii MVCC19 by qPCR. The relative transcript abundance of both genes significantly increased from the beginning to the end of growth, independent of nitrate availability in the culture medium. Differences between the genes in terms of the levels of relative expression were also found, implying that during growth in nitrate-rich or nitrate-deprived conditions C. raciborskii MVCC19 has the ability to produce different kind of PSP molecules. The presence of nifH transcripts in the nitrogen-depleted treatment confirmed that in the absence of nitrate C. raciborskii fixed atmospheric N2. Moreover, after transferring filaments to nitrate-rich conditions the synthesis of nifH mRNA continued for few hours, suggesting that cell adjustments enabling the utilization of soluble nitrogen sources are not immediate. Our results show that biosynthesis of saxitoxin and analogs in C. raciborskii is not related to nitrate availability, but rather is linked to cyanobacteria growth rate.

Application of ancient DNA to the reconstruction of past microbial assemblages and for the detection of toxic cyanobacteria in subtropical freshwater ecosystems.

Ancient DNA (aDNA) analysis of lake sediments is a promising tool for detecting shifts in past microbial assemblages in response to changing environmental conditions. We examined sediment core samples from subtropical, freshwater Laguna Blanca (Uruguay), which has been severely affected by cultural eutrophication since 1960 and where cyanobacterial blooms, particularly those of the saxitoxin-producer Cylindrospermopsis raciborskii, have been reported since the 1990s. Samples corresponding to ~1846, 1852, 2000 and 2007 AD were selected to perform denaturing gradient gel electrophoresis (DGGE) analysis of the 16S-23S rRNA intergenic transcribed spacer (ribosomal ITS) to compare their prokaryotic assemblage composition. Each stratum showed different ITS patterns, but the composition of 21st century samples was clearly different than those of mid-19th century. This compositional change was correlated with shifts in sediment organic matter and chlorophyll a content, which were significantly higher in recent samples. The presence of saxitoxin-producing cyanobacteria was addressed by quantitative real-time PCR of the sxtU gene involved in toxin biosynthesis. This gene was present only in recent samples, for which clone libraries and ITS sequencing indicated the presence of Cyanobacteria. Phylogenetic analyses identified C. raciborskii only in the 2000 sample, shortly after several years when blooms were recorded in the lake. These data suggest the utility of aDNA for the reconstruction of microbial assemblage shifts in subtropical lakes, at least on centennial scales. The application of aDNA analysis to genes involved in cyanotoxin synthesis extends the applicability of molecular techniques in palaeolimnological studies to include key microbial community characteristics of great scientific and social interest.

Cylindrospermopsin and saxitoxin synthetase genes in Cylindrospermopsis raciborskii strains from Brazilian freshwater.

The Cylindrospermopsis raciborskii population from Brazilian freshwater is known to produce saxitoxin derivatives (STX), while cylindrospermopsin (CYN), which is commonly detected in isolates from Australia and Asia continents, has thus far not been detected in South American strains. However, during the investigation for the presence of cyrA, cyrB, cyrC and cyrJ CYN synthetase genes in the genomes of four laboratory-cultured C. raciborskii Brazilian strains, the almost complete cyrA gene sequences were obtained for all strains, while cyrB and cyrC gene fragments were observed in two strains. These nucleotide sequences were translated into amino acids, and the predicted protein functions and domains confirmed their identity as CYN synthetase genes. Attempts to PCR amplify cyrJ gene fragments from the four strains were unsuccessful. Phylogenetic analysis grouped the nucleotide sequences together with their homologues found in known CYN synthetase clusters of C. raciborskii strains with high bootstrap support. In addition, fragments of sxtA, sxtB and sxtI genes involved in STX production were also obtained. Extensive LC-MS analyses were unable to detect CYN in the cultured strains, whereas the production of STX and its analogues was confirmed in CENA302, CENA305 and T3. To our knowledge, this is the first study reporting the presence of cyr genes in South American strains of C. raciborskii and the presence of sxt and cyr genes in a single C. raciborskii strain. This discovery suggests a shift in the type of cyanotoxin production over time of South American strains of C. raciborskii and contributes to the reconstruction of the evolutionary history and diversification of cyanobacterial toxins.

In silico analysis of putative paralytic shellfish poisoning toxins export proteins in cyanobacteria.

Paralytic shellfish poisoning toxins (PSTs) are a family of more than 30 natural alkaloids synthesized by dinoflagellates and cyanobacteria whose toxicity in animals is mediated by voltage-gated Na(+) channel blocking. The export of PST analogues may be through SxtF and SxtM, two putative MATE (multidrug and toxic compound extrusion) family transporters encoded in PSTs biosynthetic gene cluster (sxt). sxtM is present in every sxt cluster analyzed; however, sxtF is only present in the Cylindrospermopsis-Raphidiopsis clade. These transporters are energetically coupled with an electrochemical gradient of proton (H(+)) or sodium (Na(+)) ions across membranes. Because the functional role of PSTs remains unknown and methods for genetic manipulation in PST-producing organisms have not yet been developed, protein structure analyses will allow us to understand their function. By analyzing the sxt cluster of eight PST-producing cyanobacteria, we found no correlation between the presence of sxtF or sxtM and a specific PSTs profile. Phylogenetic analyses of SxtF/M showed a high conservation of SxtF in the Cylindrospermopsis-Raphidiopsis clade, suggesting conserved substrate affinity. Two domains involved in Na(+) and drug recognition from NorM proteins (MATE family) of Vibrio parahaemolyticus and V. cholerae are present in SxtF/M. The Na(+) recognition domain was conserved in both SxtF/M, indicating that Na(+) can maintain the role as a cation anti-transporter. Consensus motifs for toxin binding differed between SxtF and SxtM implying differential substrate binding. Through protein modeling and docking analysis, we found that there is no marked affinity between the recognition domain and a specific PST analogue. This agrees with our previous results of PST export in R. brookii D9, where we observed that the response to Na(+) incubation was similar to different analogues. These results reassert the hypothesis regarding the involvement of Na(+) in toxin export, as well as the motifs L(398)XGLQD(403) (SxtM) and L(390)VGLRD(395) (SxtF) in toxin recognition.

Toxicity phenotype does not correlate with phylogeny of Cylindrospermopsis raciborskii strains.

Cylindrospermopsis raciborskii is a species of freshwater, bloom-forming cyanobacterium. C. raciborskii produces toxins, including cylindrospermopsin (hepatotoxin) and saxitoxin (neurotoxin), although non toxin-producing strains are also observed. In spite of differences in toxicity, C. raciborskii strains comprise a monophyletic group, based upon 16S rRNA gene sequence identities (greater than 99%). We performed phylogenetic analyses; 16S rRNA gene and 16S-23S rRNA gene internally transcribed spacer (ITS-1) sequence comparisons, and genomic DNA restriction fragment length polymorphism (RFLP), resolved by pulsed-field gel electrophoresis (PFGE), of strains of C. raciborskii, obtained mainly from the Australian phylogeographic cluster. Our results showed no correlation between toxic phenotype and phylogenetic association in the Australian strains. Analyses of the 16S rRNA gene and the respective ITS-1 sequences (long L, and short S) showed an independent evolution of each ribosomal operon. The genes putatively involved in the cylindrospermopsin biosynthetic pathway were present in one locus and only in the hepatotoxic strains, demonstrating a common genomic organization for these genes and the absence of mutated or inactivated biosynthetic genes in the non toxic strains. In summary, our results support the hypothesis that the genes involved in toxicity may have been transferred as an island by processes of gene lateral transfer, rather than convergent evolution.