Comparative genomics reveals diversified CRISPR-Cas systems of globally distributed Microcystis aeruginosa, a freshwater bloom-forming cyanobacterium

CyanoStrainChip: A Novel DNA Microarray Tool for High-Throughput Detection of Environmental Cyanobacteria at the Strain Level

Detecting cyanobacteria in environments is an important concern due to their crucial roles in ecosystems, and they can form blooms with the potential to harm humans and nonhuman entities. However, the most widely used methods for high-throughput detection of environmental cyanobacteria, such as 16S rRNA sequencing, typically provide above-species-level resolution, thereby disregarding intraspecific variation. To address this, we developed a novel DNA microarray tool, termed the CyanoStrainChip, that enables strain-level comprehensive profiling of environmental cyanobacteria. The CyanoStrainChip was designed to target 1277 strains; nearly all major groups of cyanobacteria are included by implementing 43,666 genome-wide, strain-specific probes. It demonstrated strong specificity by in vitro mock community experiments. The high correlation (Pearson's R > 0.97) between probe fluorescence intensities and the corresponding DNA amounts (ranging from 1-100 ng) indicated excellent quantitative capability. Consistent cyanobacterial profiles of field samples were observed by both the CyanoStrainChip and next-generation sequencing methods. Furthermore, CyanoStrainChip analysis of surface water samples in Lake Chaohu uncovered a high intraspecific variation of abundance change within the genus Microcystis between different severity levels of cyanobacterial blooms, highlighting two toxic Microcystis strains that are of critical concern for Lake Chaohu harmful blooms suppression. Overall, these results suggest a potential for CyanoStrainChip as a valuable tool for cyanobacterial ecological research and harmful bloom monitoring to supplement existing techniques.

Characterization of Limnospira platensis PCC 9108 R-M and CRISPR-Cas systems

The filamentous cyanobacterium Limnospira platensis, formerly known as Arthrospira platensis or spirulina, is one of the most commercially important species of microalgae. Due to its high nutritional value, pharmacological and industrial applications it is extensively cultivated on a large commercial scale. Despite its widespread use, its precise manipulation is still under development due to the lack of effective genetic protocols. Genetic transformation of Limnospira has been attempted but the methods reported have not been generally reproducible in other laboratories. Knowledge of the transformation defense mechanisms is essential for understanding its physiology and for broadening their applications. With the aim to understand more about the genetic defenses of L. platensis, in this work we have identified the restriction-modification and CRISPR-Cas systems and we have cloned and characterized thirteen methylases. In parallel, we have also characterized the methylome and orphan methyltransferases using genome-wide analysis of DNA methylation patterns and RNA-seq. The identification and characterization of these enzymes will be a valuable resource to know how this strain avoids being genetically manipulated and for further genomics studies.

Ecological Dynamics of Broad- and Narrow-Host-Range Viruses Infecting the Bloom-Forming Toxic Cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa is predicted to interact and coexist with diverse broad- and narrow-host-range viruses within a bloom; however, little is known about their affects on Microcystis population dynamics. Here, we developed a real-time PCR assay for the quantification of these viruses that have different host ranges. During the sampling period, total Microcystis abundance showed two peaks in May and August with a temporary decrease in June. The Microcystis population is largely divided into three phylotypes based on internal transcribed sequences (ITS; ITS types I to III). ITS I was the dominant phylotype (66% to 88%) except in June. Although the ITS II and III phylotypes were mostly less abundant, these phylotypes temporarily increased to approximately equivalent abundances of the ITS I population in June. During the same sampling period, the abundances of the broad-host-range virus MVGF_NODE331 increased from April to May and from July to October with a temporary decrease in June, in which its dynamics were in proportion to the increase of total Microcystis abundances regardless of changes in host ITS population composition. In contrast, the narrow-host-range viruses MVG_NODE620 and Ma-LMM01 were considerably less abundant than the broad-host-range virus and generally did not fluctuate in the environment. Considering that M. aeruginosa could increase the abundance and sustain the bloom under the prevalence of the broad-host-range virus, host abundant and diverse antiviral mechanisms might contribute to coexistence with its viruses. IMPORTANCE The bloom-forming toxic cyanobacterium Microcystis aeruginosa interacts with diverse broad- and narrow-host-range viruses. However, the dynamics of the Microcystis population (at the intraspecies level) and viruses with different host ranges remain unknown. Our real-time PCR assays unveiled that the broad-host-range virus gradually increased in abundance over the sampling period, in proportion to the increase in total Microcystis abundance regardless of changes in genotypic composition. The narrow-host-range viruses were considerably less abundant than the broad-host-range virus and did not generally fluctuate in the environment. The expansion and maintenance of the Microcystis bloom even under the increased infection by the broad-host-range virus suggested that highly abundant and diverse antiviral mechanisms allowed them to coexist with viruses under selective pressure. This paper expands our knowledge about the ecological dynamics of Microcystis viruses and provides potential insights into their coexistence with their host.

Genomic comparison of Planktothrix agardhii isolates from a Lake Erie embayment

Planktothrix agardhii is a filamentous cyanobacterial species that dominates harmful algal blooms in Sandusky Bay, Lake Erie and other freshwater basins across the world. P. agardhii isolates were obtained from early (June) blooms via single filament isolation; eight have been characterized from 2016, and 12 additional isolates have been characterized from 2018 for a total of 20 new cultures. These novel isolates were processed for genomic sequencing, where reads were used to generate scaffolds and contigs which were annotated with DIAMOND BLAST hit, Pfam, and GO. Analyses include whole genome alignment to generate phylogenetic trees and comparison of genetic rearrangements between isolates. Nitrogen acquisition and metabolism was compared across isolates. Secondary metabolite production was genetically explored including microcystins, two types of aeruginosin clusters, anabaenopeptins, cyanopeptolins, microviridins, and prenylagaramides. Two common and 4 unique CRISPR-cas islands were analyzed for similar sequences across all isolates and against the known Planktothrix-specific cyanophage, PaV-LD. Overall, the uniqueness of each genome from Planktothrix blooms sampled from the same site and at similar times belies the unexplored diversity of this genus.

Analysis of Probiotic Bacteria Genomes: Comparison of CRISPR/Cas Systems and Spacer Acquisition Diversity

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) genes constitute an adaptive (acquired) defense system of bacteria and archaea. Here 72 probiotic bacteria genomes were investigated in terms of the presence of CRISPR/Cas systems and phage/plasmid invaders through spacer analysis. 49 CRISPR/Cas systems were detected within probiotic strains, namely,17 type II-A, 10 type I-C, 8 type I-E, 5 Type I-U (I-G), 4 type III-A, 2 type I-B, 1 type I-A, 1 type IV-B, and 1 type II-C. The predicted target of spacers was determined in 25 strains and consequently, three different spacer and target patterns were revealed. The diversity of CRISPR spacers provides insight and understanding to determine strain-specific invaders of probiotic bacteria as well as their relationships between strains. CRISPR systems were clarified in many studies for genomic characterization. However, recently, endogenous genome editing with CRISPR has provided an approach for various genome editing projects. Thus, in the future, producing strain-specific phage-resistant starter cultures or probiotics by endogenous genome editing methods according to phage/plasmid survey can be utilized for industrial and pharmaceutical applications. Therefore, this study intended a comprehensive investigation of CRISPR systems of probiotic bacteria.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12088-021-00971-1.

Cyanobacterial blooms in China: diversity, distribution, and cyanotoxins

Cyanobacterial blooms, which refer to the massive growth of harmful cyanobacteria, have altered the global freshwater ecosystems during the past decades. China has the largest population in the world, and it is suffering from the harmful effect of water eutrophication and cyanobacterial blooms along with rapid development of the economy and society. Research on cyanobacterial blooms and cyanotoxins in China have been overwhelmingly enhanced and emphasized during the past decades. In the present review, the research on cyanobacterial blooms in China is generally introduced, including the history of cyanobacterial bloom studies, the diversity of the bloom-forming cyanobacteria species (BFCS), and cyanotoxin studies in China. Most studies have focused on Microcystis, its blooms, and microcystins. Newly emerging blooms with the dominance of non-Microcystis BFCS have been gradually expanding to wide regions in China. Understanding the basic features of these non-Microcystis BFCS and their blooms, including their diversity, occurrence, physio-ecology, and harmful metabolites, will provide direction on future studies of cyanobacterial blooms in China.

Identification and Characterization of the CRISPR/Cas System in Staphylococcus aureus Strains From Diverse Sources

The CRISPR-Cas [clustered regularly interspaced short palindromic repeats and the CRISPR-associated genes (Cas)] system provides defense mechanisms in bacteria and archaea vs. mobile genetic elements (MGEs), such as plasmids and bacteriophages, which can either be harmful or add sequences that can provide virulence or antibiotic resistance. Staphylococcus aureus is a Gram-positive bacterium that could be the etiological agent of important soft tissue infections that can lead to bacteremia and sepsis. The role of the CRISPR-Cas system in S. aureus is not completely understood since there is a lack of knowledge about it. We analyzed 716 genomes and 1 genomic island from GENOMES-NCBI and ENA-EMBL searching for the CRISPR-Cas systems and their spacer sequences (SSs). Our bioinformatic analysis shows that only 0.83% (6/716) of the analyzed genomes harbored the CRISPR-Cas system, all of them were subtype III-A, which is characterized by the presence of the cas10/csm1 gene. Analysis of SSs showed that 91% (40/44) had no match to annotated MGEs and 9% of SSs corresponded to plasmids and bacteriophages, indicating that those phages had infected those S. aureus strains. Some of those phages have been proposed as an alternative therapy in biofilm-forming or infection with S. aureus strains, but these findings indicate that such antibiotic phage strategy would be ineffective. More research about the CRISPR/Cas system is necessary for a bigger number of S. aureus strains from different sources, so additional features can be studied.

Cell Death in Cyanobacteria: Current Understanding and Recommendations for a Consensus on Its Nomenclature

Cyanobacteria are globally widespread photosynthetic prokaryotes and are major contributors to global biogeochemical cycles. One of the most critical processes determining cyanobacterial eco-physiology is cellular death. Evidence supports the existence of controlled cellular demise in cyanobacteria, and various forms of cell death have been described as a response to biotic and abiotic stresses. However, cell death research in this phylogenetic group is a relatively young field and understanding of the underlying mechanisms and molecular machinery underpinning this fundamental process remains largely elusive. Furthermore, no systematic classification of modes of cell death has yet been established for cyanobacteria. In this work, we analyzed the state of knowledge in the field of cyanobacterial cell death. Based on that, we propose unified criterion for the definition of accidental, regulated, and programmed forms of cell death in cyanobacteria based on molecular, biochemical, and morphologic aspects following the directions of the Nomenclature Committee on Cell Death (NCCD). With this, we aim to provide a guide to standardize the nomenclature related to this topic in a precise and consistent manner, which will facilitate further ecological, evolutionary, and applied research in the field of cyanobacterial cell death.

Viruses of freshwater bloom-forming cyanobacteria: genomic features, infection strategies and coexistence with the host

Freshwater bloom-forming cyanobacteria densely grow in the aquatic environments, leading to an increase in the viral-contact rate. They possess numerous antiviral genes, as well as cell differentiation- and physiological performance-related genes, owing to genome expansion. Their genomic features and unique lifestyles suggest that they coexist with cyanoviruses in ways different from marine cyanobacteria. Furthermore, genome contents of isolated freshwater bloom-forming cyanobacterial viruses have little in common with those of marine cyanoviruses studied to date. They lack the marine cyanoviral hallmark genes that sustain photosynthetic activity and redirect host metabolism to viral reproduction; therefore, they are predicted to share metabolisms and precursor pools with host cyanobacteria to ensure efficient viral reproduction and avoid nutrient deficiencies and antiviral response. Additionally, cyanovirus-cyanobacteria coexistence strategies may change as bloom density increases. Diverse genotypic populations of cyanoviruses and hosts coexist and fluctuate under high viral-contact rate conditions, leading to their rapid coevolution through antiviral responses. The ancestral and newly evolved genotypes coexist, thereby expanding the diversity levels of host and viral populations. Bottleneck events occurring due to season-related decreases in bloom-forming species abundance provide each genotype within cyanobacterial population an equal chance to increase in prevalence during the next bloom and enhance further diversification.

Genomic Characteristics of the Toxic Bloom-Forming Cyanobacterium Microcystis aeruginosa NIES-102

Microcystis aeruginosa, a bloom-forming cyanobacterium distributed mainly in freshwater environments, can be divided into at least 12 groups (A-K and X) based on multi-locus phylogenetic analyses. In this study, we characterized the genome of microcystin-producing M. aeruginosa NIES-102, assigned to group A, isolated from Lake Kasumigaura, Japan. The complete genome sequence of M. aeruginosa NIES-102 comprised a 5.87-Mbp circular chromosome containing 5,330 coding sequences. The genome was the largest among all sequenced genomes for the species. In a comparison with the genome of M. aeruginosa NIES-843, which belongs to the same group, the microcystin biosynthetic gene cluster and CRISPR-Cas locus were highly similar. A synteny analysis revealed small-scale rearrangements between the two genomes. Genes encoding transposases were more abundant in these two genomes than in other Microcystis genomes. Our results improve our understanding of structural genomic changes and adaptation to a changing environment in the species.

Cooccurrence of Broad- and Narrow-Host-Range Viruses Infecting the Bloom-Forming Toxic Cyanobacterium Microcystis aeruginosa

Viruses play important roles in regulating the abundance and composition of bacterial populations in aquatic ecosystems. The bloom-forming toxic cyanobacterium Microcystis aeruginosa is predicted to interact with diverse cyanoviruses, resulting in Microcystis population diversification. However, current knowledge of the genomes from these viruses and their infection programs is limited to those of Microcystis virus Ma-LMM01. Here, we performed a time series sampling at a small pond in Japan during a Microcystis bloom and then investigated the genomic information and transcriptional dynamics of Microcystis-interacting viruses using metagenomic and metatranscriptomic approaches. We identified 15 viral genomic fragments classified into three groups, groups I (including Ma-LMM01), II (high abundance and transcriptional activity), and III (new lineages). According to the phylogenetic distribution of Microcystis strains possessing spacers against each viral group, the group II-original viruses interacted with all three phylogenetically distinct Microcystis population types (phylotypes), whereas the groups I and III-original viruses interacted with only one or two phylotypes, indicating the cooccurrence of broad- (group II) and narrow (groups I and III)-host-range viruses in the bloom. These viral fragments showed the highest transcriptional levels during daytime regardless of their genomic differences. Interestingly, M. aeruginosa expressed antiviral defense genes in the environment, unlike what was seen with an Ma-LMM01 infection in a previous culture experiment. Given that broad-host-range viruses often induce antiviral responses within alternative hosts, our findings suggest that such antiviral responses might inhibit viral multiplication, mainly that of broad-host-range viruses like those in group II.IMPORTANCE The bloom-forming toxic cyanobacterium Microcystis aeruginosa is thought to have diversified its population through the interactions between host and viruses in antiviral defense systems. However, current knowledge of viral genomes and infection programs is limited to those of Microcystis virus Ma-LMM01, which was a narrow host range in which it can escape from the highly abundant host defense systems. Our metagenomic approaches unveiled the cooccurrence of narrow- and broad-host-range Microcystis viruses, which included fifteen viral genomic fragments from Microcystis blooms that were classified into three groups. Interestingly, Microcystis antiviral defense genes were expressed against viral infection in the environment, unlike what was seen in a culture experiment with Ma-LMM01. Given that viruses with a broad host range often induce antiviral responses within alternative hosts, our findings suggest that antiviral responses inhibit viral reproduction, especially that of broad-range viruses like those in group II. This paper augments our understanding of the interactions between M. aeruginosa and its viruses and fills an important knowledge gap.

Interaction between cyanophage MaMV-DC and eight Microcystis strains, revealed by genetic defense systems

Cyanophage MaMV-DC is a member of Myoviridae that was reported to specifically infect and lyse Microcystis aeruginosa FACHB-524 among 21 selected cyanobacterial strains. We reidentified the infection specificity of MaMV-DC among seven other Microcystis strains of different species. In our experiments, MaMV-DC infected three Microcystis strains but did not form plaque in Microcystis lawns. This indicated that MaMV-DC is at least a genus- rather than strain-specific virus. Cyanophage MaMV-DC genes were transcribed in M. aeruginosa FACHB-524, M. flos-aquae TF09, M. aeruginosa TA09 and M. wesenbergii DW09, and the growth of these Microcystis strains was inhibited by the addition of MaMV-DC. The predicted defense of eight Microcystis strains by CRISPR-Cas systems has shown mixed consistency with the infection experiment results, suggesting other defense or anti-defense systems play roles during infection process. Restriction-modification (RM) system analysis revealed an abundance of four types of RM proteins that may play roles in defense against cyanophages.

Insertions within the Saxitoxin Biosynthetic Gene Cluster Result in Differential Toxin Profiles

The neurotoxin saxitoxin and related paralytic shellfish toxins are produced by multiple species of cyanobacteria and dinoflagellates. This study investigates the two saxitoxin-producing strains of Scytonema crispum, CAWBG524 and CAWBG72, isolated in New Zealand. Each strain was previously reported to have a distinct paralytic shellfish toxin profile, a rare observation between strains within the same species. Sequencing of the saxitoxin biosynthetic clusters ( sxt) from S. crispum CAWBG524 and S. crispum CAWBG72 revealed the largest sxt gene clusters described to date. The distinct toxin profiles of each strain were correlated to genetic differences in sxt tailoring enzymes, specifically the open-reading frame disruption of the N-21 sulfotransferase sxtN, adenylylsulfate kinase sxtO, and the C-11 dioxygenase sxtDIOX within S. crispum CAWBG524 via genetic insertions. Heterologous overexpression of SxtN allowed for the proposal of saxitoxin and 3'-phosphoadenosine 5'-phosphosulfate as substrate and cofactor, respectively, using florescence binding assays. Further, catalytic activity of SxtN was confirmed by the in vitro conversion of saxitoxin to the N-21 sulfonated analog gonyautoxin 5, making this the first known report to biochemically confirm the function of a sxt tailoring enzyme. Further, SxtN could not convert neosaxitoxin to its N-21 sulfonated analog gonyautoxin 6, indicating paralytic shellfish toxin biosynthesis most likely occurs along a predefined route. In this study, we identified key steps toward the biosynthetic conversation of saxitoxin to other paralytic shellfish toxins.

Microbial Diversity and Toxin Risk in Tropical Freshwater Reservoirs of Cape Verde

The Cape Verde islands are part of the African Sahelian arid belt that possesses an erratic rain pattern prompting the need for water reservoirs, which are now critical for the country’s sustainability. Worldwide, freshwater cyanobacterial blooms are increasing in frequency due to global climate change and the eutrophication of water bodies, particularly in reservoirs. To date, there have been no risk assessments of cyanobacterial toxin production in these man-made structures. We evaluated this potential risk using 16S rRNA gene amplicon sequencing and full metagenome sequencing in freshwater reservoirs of Cape Verde. Our analysis revealed the presence of several potentially toxic cyanobacterial genera in all sampled reservoirs. Faveta potentially toxic and bloom-forming Microcystis sp., dominated our samples, while a Cryptomonas green algae and Gammaproteobacteria dominated Saquinho and Poilão reservoirs. We reconstructed and assembled the Microcystis genome, extracted from the metagenome of bulk DNA from Faveta water. Phylogenetic analysis of Microcystis cf. aeruginosa CV01’s genome revealed its close relationship with other Microcystis genomes, as well as clustering with other continental African strains, suggesting geographical coherency. In addition, it revealed several clusters of known toxin-producing genes. This survey reinforces the need to better understand the country’s microbial ecology as a whole of water reservoirs on the rise.

The highly heterogeneous methylated genomes and diverse restriction-modification systems of bloom-forming Microcystis

The occurrence of harmful Microcystis blooms is increasing in frequency in a myriad of freshwater ecosystems. Despite considerable research pertaining to the cause and nature of these blooms, the molecular mechanisms behind the cosmopolitan distribution and phenotypic diversity in Microcystis are still unclear. We compared the patterns and extent of DNA methylation in three strains of Microcystis, PCC 7806SL, NIES-2549 and FACHB-1757, using Single Molecule Real-Time (SMRT) sequencing technology. Intact restriction-modification (R-M) systems were identified from the genomes of these strains, and from two previously sequenced strains of Microcystis, NIES-843 and TAIHU98. A large number of methylation motifs and R-M genes were identified in these strains, which differ substantially among different strains. Of the 35 motifs identified, eighteen had not previously been reported. Strain NIES-843 contains a larger number of total putative methyltransferase genes than have been reported previously from any bacterial genome. Genomic comparisons reveal that methyltransferases (some partial) may have been acquired from the environment through horizontal gene transfer.

Comparative Genomics of the Baltic Sea Toxic Cyanobacteria Nodularia spumigena UHCC 0039 and Its Response to Varying Salinity

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.

Incomplete Selective Sweeps of Microcystis Population Detected by the Leader-End CRISPR Fragment Analysis in a Natural Pond

The freshwater cyanobacterium Microcystis aeruginosa frequently forms toxic massive blooms and exists in an arms race with its infectious phages in aquatic natural environments, and as a result, has evolved extremely diverse and elaborate antiviral defense systems, including the clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated genes (Cas) system. Here, to assess Microcystis population dynamics associated with exogenous mobile genetic elements such as phages and plasmids, we examined the temporal variation in CRISPR genotypes (CTs) by analyzing spacer sequences detected in a natural pond between June and October 2013 when a cyanobacterial bloom occurred. A total of 463,954 high-quality leader-end CRISPR sequences were obtained and the sequences containing spacers were classified into 31 previously reported CTs and 68 new CTs based on the shared order of the leader-end spacers. CT19 was the most dominant genotype (32%) among the 16 most common CTs, followed by CT52 (14%) and CT58 (9%). Spacer repertoires of CT19 showed mainly two different types; CT19_origin, which was identical to the CT19 spacer repertoire of previously isolated strains, and CT19_new+, which contained a new spacer at the leader-end of the CRISPR region of CT19_origin, which were present in almost equal abundance, accounting for up to 99.94% of CT19 sequences. Surprisingly, we observed the spacer repertoires of the second to tenth spacers of CT19_origin at the most leader-end of proto-genotype sequences of CT19_origin. These were observed during the sampling in this study and our previous study at the same ecosystem in 2010 and 2011, suggesting these CTs persisted from 2011 to 2013 in spite of phage pressure. The leader-end variants were observed in other CT genotypes. These findings indicated an incomplete selective sweep of Microcystis populations. We explained the phenomenon as follow; the abundance of Microcystis varied seasonally and drastically, resulting that Microcystis populations experience a bottleneck once a year, and thereby founder effects following a bottleneck mean that older CTs have an equal chance of increasing in prevalence as the CTs generated following acquisition of newer spacers.

Genome variation in nine co-occurring toxic Cylindrospermopsis raciborskii strains

Cyanobacteria form harmful algal blooms and are highly adapted to a range of habitats, in part due to their phenotype plasticity. This plasticity is partially the result of co-existence of multiple strains within a single population. The toxic cyanobacterium Cylindrospermopsis raciborskii has remarkable phenotypic plasticity, strain variation and environmental adaptation resulting in an expansion of its global range. To understand the genetic basis of the high level of plasticity within a C. raciborskii population, the genomes of nine co-occurring strains were compared. The strains differed in morphology, toxin cell quotas and physiology, despite being obtained from a single water sample. Comparative genomics showed that three coiled strains were 3.9 Mbp in size, with 3544 ± 11 genes, while straight strains were 3.8 Mbp in size, with 3485 ± 20 genes. The core proteome comprised 86% of the genome and consisted of 2891 orthologous groups (OGs), whereas the variable genome comprised ∼14% (847 OGs), and the strain specific genome only ∼1% (433 OGs).There was a high proportion of variable strain-specific genes for the very closely related strains, which may underpin strain differentiation. The variable genes were associated with environmental responses and adaptation, particularly phage defence, DNA repair, membrane transport, and stress, illustrative of the adaptability of the strains in response to environmental and biological stressors. This study shows that high genomic variability exists between co-occurring strains and may be the basis of strain phenotypic differences and plasticity of populations. Therefore management and prediction of blooms of this harmful species requires different approaches to capture this strain variability.

Metatranscriptomics analysis of cyanobacterial aggregates during cyanobacterial bloom period in Lake Taihu, China

Molecular mechanism of interaction between the bloom-forming cyanobacterial species and attached microbios within cyanobacterial aggregates has not been elucidated yet and understanding of which would help to unravel the cyanobacteria bloom-forming mechanism. In this study, we profiled the metabolically active community by high-throughput metatranscriptome sequencing from cyanobacterial aggregates during cyanobacterial bloom period in Lake Taihu, China. A total of 308 million sequences were obtained using the HiSeq 2500 sequencing platform, which provided a great sequence coverage to carry out the in-depth taxonomic classification, functional classification, and metabolic pathway analysis of the cyanobacterial aggregates. The results show that bacteria dominated in cyanobacterial aggregates, accounting for more than 96.66% of total sequences. Microcystis was the most abundant genus, accounted for 26.80% of total assigned sequences at the genus level in cyanobacterial aggregates community; however, Proteobacteria (46.20%) was found to be as the most abundant active bacterial populations at the phylum level. More importantly, nitrogen, phosphonate, and phosphinate metabolism which associated with eutrophication were found in this study. Especially, the enzymes and organisms relating to denitrification and anammox of nitrogen metabolism, which reduced nitrogen concentration by reducing nitrate to nitrogen to inhibit the eutrophication, were first discovered in Lake Taihu during cyanobacterial bloom period. The present study provides a snapshot of metatranscriptome for cyanobacterial aggregates in Lake Taihu and demonstrates that cyanobacterial aggregates could play a key role in the nitrogen cycle in eutrophic water.

Metagenomic analysis reveals potential interactions in an artificial coculture

Disentangling the interactions between cyanobacteria and associated bacterial community is important for understanding the mechanisms that mediate the formation of cyanobacterial blooms in freshwater ecosystems. Despite the fact that a metagenomic approach enables researchers to profile the structure of microbial communities associated with cyanobacteria, reconstructing genome sequences for all members remains inefficient, due to the inherent enormous microbial diversity. Here, we have established a stable coculture system under high salinity, originally from a mixture of an axenic cyanobacterium Synechococcus sp. PCC 7002 and a non-axenic bloom-forming cyanobacterium Microcystis colony. Metagenomic analysis showed that the coculture consists of S. sp. PCC 7002 and two heterotrophic bacteria, designated as Pseudomonas stutzeri TAIHU and Mesorhizobium sp. TAIHU, respectively. And near-complete genome sequences of both bacteria were reconstructed from the metagenomic dataset with an average completeness of 99.8%. Genome-wide pathway analysis revealed that M. sp. TAIHU carried all the genes involved in the de novo biosynthesis of cobalamin, which is required by S. sp. PCC 7002 for growth. To cope with the high salinity in the coculture, experimental evidence demonstrated that S. sp. PCC 7002 would synthesize the compatible solutes including sucrose and glucosylglycerol, which are supposed to be exploited by both heterotrophic bacteria as potential carbon and/or nitrogen sources. Furthermore, the genes encoding for the biosynthesis of the ectoine, another common osmolyte are found exclusively in P. stutzeri TAIHU, while the genes responsible for the catabolism of ectoine and its derives are present only in M. sp. TAIHU. These genomic evidence indicates beneficial interaction between three members in the coculture. Establishment of the coculture system with relative simplicity provides a useful model system for investigating the interspecies interactions, and genome sequences of both bacteria associated with Microcystis bloom described here will facilitate the researcher to elucidate the role of these heterotrophic bacteria in the formation and maintenance of cyanobacterial bloom in freshwater ecosystem.

Genome sequences of lower Great Lakes Microcystis sp. reveal strain-specific genes that are present and expressed in western Lake Erie blooms

Blooms of the potentially toxic cyanobacterium Microcystis are increasing worldwide. In the Laurentian Great Lakes they pose major socioeconomic, ecological, and human health threats, particularly in western Lake Erie. However, the interpretation of "omics" data is constrained by the highly variable genome of Microcystis and the small number of reference genome sequences from strains isolated from the Great Lakes. To address this, we sequenced two Microcystis isolates from Lake Erie (Microcystis aeruginosa LE3 and M. wesenbergii LE013-01) and one from upstream Lake St. Clair (M. cf aeruginosa LSC13-02), and compared these data to the genomes of seventeen Microcystis spp. from across the globe as well as one metagenome and seven metatranscriptomes from a 2014 Lake Erie Microcystis bloom. For the publically available strains analyzed, the core genome is ~1900 genes, representing ~11% of total genes in the pan-genome and ~45% of each strain's genome. The flexible genome content was related to Microcystis subclades defined by phylogenetic analysis of both housekeeping genes and total core genes. To our knowledge this is the first evidence that the flexible genome is linked to the core genome of the Microcystis species complex. The majority of strain-specific genes were present and expressed in bloom communities in Lake Erie. Roughly 8% of these genes from the lower Great Lakes are involved in genome plasticity (rapid gain, loss, or rearrangement of genes) and resistance to foreign genetic elements (such as CRISPR-Cas systems). Intriguingly, strain-specific genes from Microcystis cultured from around the world were also present and expressed in the Lake Erie blooms, suggesting that the Microcystis pangenome is truly global. The presence and expression of flexible genes, including strain-specific genes, suggests that strain-level genomic diversity may be important in maintaining Microcystis abundance during bloom events.

Genome analysis and gene nblA identification of Microcystis aeruginosa myovirus (MaMV-DC) reveal the evidence for horizontal gene transfer events between cyanomyovirus and host

The genome sequence, genetic characterization and nblA gene function of Microcystis aeruginosa myovirus isolated from Lake Dianchi in China (MaMV-DC) have been analysed. The genome DNA is 169 223 bp long, with 170 predicted protein-coding genes (001L–170L) and a tRNA gene. About one-sixth of these genes have homologues in the host cyanobacteria M. aeruginosa. The genome carries a gene homologous to host nblA, which encodes a protein involved in the degradation of cyanobacterial phycobilisome. Its expression during MaMV-DC infection was confirmed by reverse transcriptase PCR and Western blot detection and abundant expression was companied by the significant decline of phycocyanin content and massive release of progeny MaMV-DC. In addition, expressing MaMV-DC nblA reduced the phycocyanin peak and the phycocyanin to chlorophyll ratio in model cyanobacteria. These results confirm that horizontal gene transfer events have occurred between cyanobacterial host and cyanomyovirus and suggest that MaMV-DC carrying host-derived genes (such as 005L, that codes for NblA) is responsible for more efficient expression of cyanophage genes and release of progeny cyanophage. This study provides novel insight into the horizontal gene transfer in cyanophage and the interactions between cyanophage and their host.

A review of the global ecology, genomics, and biogeography of the toxic cyanobacterium, Microcystis spp

This review summarizes the present state of knowledge regarding the toxic, bloom-forming cyanobacterium, Microcystis, with a specific focus on its geographic distribution, toxins, genomics, phylogeny, and ecology. A global analysis found documentation suggesting geographic expansion of Microcystis, with recorded blooms in at least 108 countries, 79 of which have also reported the hepatatoxin microcystin. The production of microcystins (originally "Fast-Death Factor") by Microcystis and factors that control synthesis of this toxin are reviewed, as well as the putative ecophysiological roles of this metabolite. Molecular biological analyses have provided significant insight into the ecology and physiology of Microcystis, as well as revealed the highly dynamic, and potentially unstable, nature of its genome. A genetic sequence analysis of 27 Microcystis species, including 15 complete/draft genomes are presented. Using the strictest biological definition of what constitutes a bacterial species, these analyses indicate that all Microcystis species warrant placement into the same species complex since the average nucleotide identity values were above 95%, 16S rRNA nucleotide identity scores exceeded 99%, and DNA-DNA hybridization was consistently greater than 70%. The review further provides evidence from around the globe for the key role that both nitrogen and phosphorus play in controlling Microcystis bloom dynamics, and the effect of elevated temperature on bloom intensification. Finally, highlighted is the ability of Microcystis assemblages to minimize their mortality losses by resisting grazing by zooplankton and bivalves, as well as viral lysis, and discuss factors facilitating assemblage resilience.

Orthocaspase and toxin-antitoxin loci rubbing shoulders in the genome of Microcystis aeruginosa PCC 7806

Programmed cell death in multicellular organisms is a coordinated and precisely regulated process. On the other hand, in bacteria we have little clue about the network of interacting molecules that result in the death of a single cell within a population or the death of almost complete population, such as often observed in cyanobacterial blooms. With the recent discovery that orthocaspase MaOC1 of the cyanobacterium Microcystis aeruginosa is an active proteolytic enzyme, we have gained a possible hint about at least one step in the process, but the picture is far from complete. Interestingly, the genomic context of MaOC1 revealed the presence of multiple copies of genes that belong to toxin-antitoxin modules. It has been speculated that these also play a role in bacterial programmed cell death. The discovery of two components linked to cell death within the same genomic region could open new ways to deciphering the underlying mechanisms of cyanobacterial cell death.

Metagenomic Analysis Reveals Symbiotic Relationship among Bacteria in Microcystis-Dominated Community

Microcystis bloom, a cyanobacterial mass occurrence often found in eutrophicated water bodies, is one of the most serious threats to freshwater ecosystems worldwide. In nature, Microcystis forms aggregates or colonies that contain heterotrophic bacteria. The Microcystis-bacteria colonies were persistent even when they were maintained in lab culture for a long period. The relationship between Microcystis and the associated bacteria was investigated by a metagenomic approach in this study. We developed a visualization-guided method of binning for genome assembly after total colony DNA sequencing. We found that the method was effective in grouping sequences and it did not require reference genome sequence. Individual genomes of the colony bacteria were obtained and they provided valuable insights into microbial community structures. Analysis of metabolic pathways based on these genomes revealed that while all heterotrophic bacteria were dependent upon Microcystis for carbon and energy, Vitamin B12 biosynthesis, which is required for growth by Microcystis, was accomplished in a cooperative fashion among the bacteria. Our analysis also suggests that individual bacteria in the colony community contributed a complete pathway for degradation of benzoate, which is inhibitory to the cyanobacterial growth, and its ecological implication for Microcystis bloom is discussed.

Next-generation sequencing (NGS) for assessment of microbial water quality: current progress, challenges, and future opportunities

Water quality is an emergent property of a complex system comprised of interacting microbial populations and introduced microbial and chemical contaminants. Studies leveraging next-generation sequencing (NGS) technologies are providing new insights into the ecology of microbially mediated processes that influence fresh water quality such as algal blooms, contaminant biodegradation, and pathogen dissemination. In addition, sequencing methods targeting small subunit (SSU) rRNA hypervariable regions have allowed identification of signature microbial species that serve as bioindicators for sewage contamination in these environments. Beyond amplicon sequencing, metagenomic and metatranscriptomic analyses of microbial communities in fresh water environments reveal the genetic capabilities and interplay of waterborne microorganisms, shedding light on the mechanisms for production and biodegradation of toxins and other contaminants. This review discusses the challenges and benefits of applying NGS-based methods to water quality research and assessment. We will consider the suitability and biases inherent in the application of NGS as a screening tool for assessment of biological risks and discuss the potential and limitations for direct quantitative interpretation of NGS data. Secondly, we will examine case studies from recent literature where NGS based methods have been applied to topics in water quality assessment, including development of bioindicators for sewage pollution and microbial source tracking, characterizing the distribution of toxin and antibiotic resistance genes in water samples, and investigating mechanisms of biodegradation of harmful pollutants that threaten water quality. Finally, we provide a short review of emerging NGS platforms and their potential applications to the next generation of water quality assessment tools.