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

Insights into cyanobacterial blooms through the lens of omics

Cyanobacteria are oxygen-producing photosynthetic bacteria that convert carbon dioxide into biomass upon exposure to sunlight. However, favorable conditions cause harmful cyanobacterial blooms (HCBs), which are the dense accumulation of biomass at the water surface or subsurface, posing threats to freshwater ecosystems and human health. Understanding the mechanisms underlying cyanobacterial bloom formation is crucial for effective management. In this regard, recent advancements in omics technologies have provided valuable insights into HCBs, which have raised expectations to develop more effective control methods in the near future. This literature review aims to present the genomic architecture, adaptive mechanisms, microbial interactions, and ecological impacts of HCBs through the lens of omics. Genomic analysis indicates that the genome plasticity of cyanobacteria has enabled their resilience and effective adaptation to environmental changes. Transcriptomic investigations have revealed that cyanobacteria use various strategies for adapting to environmental stress. Additionally, metagenomic and metatranscriptomic analyses have emphasized the significant role of the microbial community in regulating HCBs. Finally, we offer perspectives on potential opportunities for further research in this field.

Profiling the interplay and coevolution of Microcystis aeruginosa and cyanosiphophage Mic1

The cyanosiphophage Mic1 specifically infects the bloom-forming Microcystis aeruginosa FACHB 1339 from Lake Chaohu, China. Previous genomic analysis showed that its 92,627 bp double-stranded DNA genome consists of 98 putative open reading frames, 63% of which are of unknown function. Here, we investigated the transcriptome dynamics of Mic1 and its host using RNA sequencing. In the early, middle, and late phases of the 10 h lytic cycle, the Mic1 genes are sequentially expressed and could be further temporally grouped into two distinct clusters in each phase. Notably, six early genes, including gp49 that encodes a TnpB-like transposase, immediately reach the highest transcriptional level in half an hour, representing a pioneer cluster that rapidly regulates and redirects host metabolism toward the phage. An in-depth analysis of the host transcriptomic profile in response to Mic1 infection revealed significant upregulation of a polyketide synthase pathway and a type III-B CRISPR system, accompanied by moderate downregulation of the photosynthesis and key metabolism pathways. The constant increase of phage transcripts and relatively low replacement rate over the host transcripts indicated that Mic1 utilizes a unique strategy to gradually take over a small portion of host metabolism pathways after infection. In addition, genomic analysis of a less-infective Mic1 and a Mic1-resistant host strain further confirmed their dynamic interplay and coevolution via the frequent horizontal gene transfer. These findings provide insights into the mutual benefit and symbiosis of the highly polymorphic cyanobacteria M. aeruginosa and cyanophages.

IMPORTANCE

The highly polymorphic Microcystis aeruginosa is one of the predominant bloom-forming cyanobacteria in eutrophic freshwater bodies and is infected by diverse and abundant cyanophages. The presence of a large number of defense systems in M. aeruginosa genome suggests a dynamic interplay and coevolution with the cyanophages. In this study, we investigated the temporal gene expression pattern of Mic1 after infection and the corresponding transcriptional responses of its host. Moreover, the identification of a less-infective Mic1 and a Mic1-resistant host strain provided the evolved genes in the phage-host coevolution during the multiple-generation cultivation in the laboratory. Our findings enrich the knowledge on the interplay and coevolution of M. aeruginosa and its cyanophages and lay the foundation for the future application of cyanophage as a potential eco-friendly and bio-safe agent in controlling the succession of harmful cyanobacterial blooms.

Multi-year molecular quantification and 'omics analysis of Planktothrix-specific cyanophage sequences from Sandusky Bay, Lake Erie

Introduction

Planktothrix agardhii is a microcystin-producing cyanobacterium found in Sandusky Bay, a shallow and turbid embayment of Lake Erie. Previous work in other systems has indicated that cyanophages are an important natural control factor of harmful algal blooms. Currently, there are few cyanophages that are known to infect P. agardhii, with the best-known being PaV-LD, a tail-less cyanophage isolated from Lake Donghu, China. Presented here is a molecular characterization of Planktothrix specific cyanophages in Sandusky Bay.

Methods and Results

Putative Planktothrix-specific viral sequences from metagenomic data from the bay in 2013, 2018, and 2019 were identified by two approaches: homology to known phage PaV-LD, or through matching CRISPR spacer sequences with Planktothrix host genomes. Several contigs were identified as having viral signatures, either related to PaV-LD or potentially novel sequences. Transcriptomic data from 2015, 2018, and 2019 were also employed for the further identification of cyanophages, as well as gene expression of select viral sequences. Finally, viral quantification was tested using qPCR in 2015-2019 for PaV-LD like cyanophages to identify the relationship between presence and gene expression of these cyanophages. Notably, while PaV-LD like cyanophages were in high abundance over the course of multiple years (qPCR), transcriptomic analysis revealed only low levels of viral gene expression.

Discussion

This work aims to provide a broader understanding of Planktothrix cyanophage diversity with the goals of teasing apart the role of cyanophages in the control and regulation of harmful algal blooms and designing monitoring methodology for potential toxin-releasing lysis events.

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.

Phylogenomics of five Pseudanabaena cyanophages and evolutionary traces of horizontal gene transfer

BACKGROUND

Along with the fast development and urbanization in developing countries, the waterbodies aside the growing cities become heavily polluted and highly eutrophic, thus leading to the seasonal outbreak of cyanobacterial bloom. Systematic isolation and characterization of freshwater cyanophages might provide a biological solution to control the awful blooms. However, genomic sequences and related investigations on the freshwater cyanophages remain very limited to date.

RESULTS

Following our recently reported five cyanophages Pam1~Pam5 from Lake Chaohu in China, here we isolated another five cyanophages, termed Pan1~Pan5, which infect the cyanobacterium Pseudanabaena sp. Chao 1811. Whole-genome sequencing showed that they all contain a double-stranded DNA genome of 37.2 to 72.0 kb in length, with less than half of the putative open reading frames annotated with known functions. Remarkably, the siphophage Pan1 encodes an auxiliary metabolic gene phoH and constitutes, together with the host, a complete queuosine modification pathway. Proteomic analyses revealed that although Pan1~Pan5 are distinct from each other in evolution, Pan1 and Pan3 are somewhat similar to our previously identified cyanophages Pam3 and Pam1 at the genomic level, respectively. Moreover, phylogenetic analyses suggested that Pan1 resembles the α-proteobacterial phage vB_DshS-R5C, revealing direct evidence for phage-mediated horizontal gene transfer between cyanobacteria and α-proteobacteria.

CONCLUSION

In addition to the previous reports of Pam1~Pam5, the present findings on Pan1~Pan5 largely enrich the library of reference freshwater cyanophages. The abundant genomic information provides a pool to identify novel genes and proteins of unknown function. Moreover, we found for the first time the evolutionary traces in the cyanophage that horizontal gene transfer might occur at the level of not only inter-species, but even inter-phylum. It indicates that the bacteriophage or cyanophage could be developed as a powerful tool for gene manipulation among various species or phyla.

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.

Comparative genomic analysis of five freshwater cyanophages and reference-guided metagenomic data mining

BACKGROUND

As important producers using photosynthesis on Earth, cyanobacteria contribute to the oxygenation of atmosphere and the primary production of biosphere. However, due to the eutrophication of urban waterbodies and global warming, uncontrollable growth of cyanobacteria usually leads to the seasonal outbreak of cyanobacterial blooms. Cyanophages, a group of viruses that specifically infect and lyse cyanobacteria, are considered as potential environment-friendly agents to control the harmful blooms. Compared to the marine counterparts, only a few freshwater cyanophages have been isolated and genome sequenced to date, largely limiting their characterizations and applications.

RESULTS

Here, we isolated five freshwater cyanophages varying in tail morphology, termed Pam1~Pam5, all of which infect the cyanobacterium Pseudanabaena mucicola Chao 1806 that was isolated from the bloom-suffering Lake Chaohu in Anhui, China. The whole-genome sequencing showed that cyanophages Pam1~Pam5 all contain a dsDNA genome, varying in size from 36 to 142 Kb. Phylogenetic analyses suggested that Pam1~Pam5 possess different DNA packaging mechanisms and are evolutionarily distinct from each other. Notably, Pam1 and Pam5 have lysogeny-associated gene clusters, whereas Pam2 possesses 9 punctuated DNA segments identical to the CRISPR spacers in the host genome. Metagenomic data-based calculation of the relative abundance of Pam1~Pam5 at the Nanfei estuary towards the Lake Chaohu revealed that the short-tailed Pam1 and Pam5 account for the majority of the five cyanophages. Moreover, comparative analyses of the reference genomes of Pam1~Pam5 and previously reported cyanophages enabled us to identify three circular and seven linear contigs of virtual freshwater cyanophages from the metagenomic data of the Lake Chaohu.

CONCLUSIONS

We propose a high-throughput strategy to systematically identify cyanophages based on the currently available metagenomic data and the very limited reference genomes of experimentally isolated cyanophages. This strategy could be applied to mine the complete or partial genomes of unculturable bacteriophages and viruses. Transformation of the synthesized whole genomes of these virtual phages/viruses to proper hosts will enable the rescue of bona fide viral particles and eventually enrich the library of microorganisms that exist on Earth. Video abstract.

A Review of Cyanophage–Host Relationships: Highlighting Cyanophages as a Potential Cyanobacteria Control Strategy

Harmful algal blooms (HABs) are naturally occurring phenomena, and cyanobacteria are the most commonly occurring HABs in freshwater systems. Cyanobacteria HABs (cyanoHABs) negatively affect ecosystems and drinking water resources through the production of potent toxins. Furthermore, the frequency, duration, and distribution of cyanoHABs are increasing, and conditions that favor cyanobacteria growth are predicted to increase in the coming years. Current methods for mitigating cyanoHABs are generally short-lived and resource-intensive, and have negative impacts on non-target species. Cyanophages (viruses that specifically target cyanobacteria) have the potential to provide a highly specific control strategy with minimal impacts on non-target species and propagation in the environment. A detailed review (primarily up to 2020) of cyanophage lifecycle, diversity, and factors influencing infectivity is provided in this paper, along with a discussion of cyanophage and host cyanobacteria relationships for seven prominent cyanoHAB-forming genera in North America, including: Synechococcus, Microcystis, Dolichospermum, Aphanizomenon, Cylindrospermopsis, Planktothrix, and Lyngbya. Lastly, factors affecting the potential application of cyanophages as a cyanoHAB control strategy are discussed, including efficacy considerations, optimization, and scalability for large-scale applications.

A comparative study of metatranscriptomic assessment methods to characterize Microcystis blooms

Harmful algal blooms are increasing in duration and severity globally, resulting in increased research interest. The use of genetic sequencing technologies has provided a wealth of opportunity to advance knowledge, but also poses a risk to that knowledge if handled incorrectly. The vast numbers of sequence processing tools and protocols provide a method to test nearly every hypothesis, but each method has inherent strengths and weaknesses. Here, we tested six methods to classify and quantify metatranscriptomic activity from a harmful algal bloom dominated by Microcystis spp. Three online tools were evaluated (Kaiju, MG-RAST, and GhostKOALA) in addition to three local tools that included a command line BLASTx approach, recruitment of reads to individual Microcystis genomes, and recruitment to a combined Microcystis composite genome generated from sequenced isolates with complete, closed genomes. Based on the analysis of each tool presented in this study, two recommendations are made that are dependent on the hypothesis to be tested. For researchers only interested in the function and physiology of Microcystis spp., read recruitments to the composite genome, referred to as "Frankenstein's Microcystis", provided the highest total estimates of transcript expression. However, for researchers interested in the entire bloom microbiome, the online GhostKOALA annotation tool, followed by subsequent read recruitments, provided functional and taxonomic characterization, in addition to transcript expression estimates. This study highlights the critical need for careful evaluation of methods before data analysis.

Tracing the active genetic diversity of Microcystis and Microcystis phage through a temporal survey of Taihu

Harmful algal blooms are commonly thought to be dominated by a single genus, but they are not homogenous communities. Current approaches, both molecular and culture-based, often overlook fine-scale variations in community composition that can influence bloom dynamics. We combined homology-based searches (BLASTX) and phylogenetics to distinguish and quantify Microcystis host and phage members across a summer season during a 2014 Microcystis- dominated bloom that occurred in Lake Tai (Taihu), China. We found 47 different genotypes of the Microcystis-specific DNA-dependent RNA polymerase (rpoB), which included several morphospecies. Microcystis flos-aquae and Microcystis wesenbergii accounted for ~86% of total Microcystis transcripts, while the more commonly studied Microcystis aeruginosa only accounted for ~7%. Microcystis genotypes were classified into three temporal groups according to their expression patterns across the course of the bloom: early, constant and late. All Microcystis morphospecies were present in each group, indicating that expression patterns were likely dictated by competition driven by environmental factors, not phylogeny. We identified three primary Microcystis-infecting phages based on the viral terminase, including a novel Siphoviridae phage that may be capable of lysogeny. Within our dataset, Myoviridae phages consistent with those infecting Microcystis in a lytic manner were positively correlated to the early host genotypes, while the Siphoviridae phages were positively correlated to the late host genotypes, when the Myoviridae phages express putative genetic markers for lysogeny. The expression of genes in the microcystin-encoding mcy cassette was estimated using mcyA, which revealed 24 Microcystis-specific genotypes that were negatively correlated to the early host genotypes. Of all environmental factors measured, pH best described the temporal shift in the Microcystis community genotypic composition, promoting hypotheses regarding carbon concentration mechanisms and oxidative stress. Our work expounds on the complexity of HAB events, using a well-studied dataset to highlight the need for increased resolution of community dynamics.

Dissolved Microcystin Release Coincident with Lysis of a Bloom Dominated by Microcystis spp. in Western Lake Erie Attributed to a Novel Cyanophage

Western Lake Erie (Laurentian Great Lakes) is prone to annual cyanobacterial harmful algal blooms (cHABs) dominated by Microcystis spp. that often yield microcystin toxin concentrations exceeding the federal EPA recreational contact advisory of 8 μg liter-1 In August 2014, microcystin levels were detected in finished drinking water above the World Health Organization 1.0 μg liter-1 threshold for consumption, leading to a 2-day disruption in the supply of drinking water for >400,000 residents of Toledo, Ohio (USA). Subsequent metatranscriptomic analysis of the 2014 bloom event provided evidence that release of toxin into the water supply was likely caused by cyanophage lysis that transformed a portion of the intracellular microcystin pool into the dissolved fraction, rendering it more difficult to eliminate during treatment. In August 2019, a similar increase in dissolved microcystins at the Toledo water intake was coincident with a viral lytic event caused by a phage consortium different in composition from what was detected following the 2014 Toledo water crisis. The most abundant viral sequence in metagenomic data sets was a scaffold from a putative member of the Siphoviridae, distinct from the Ma-LMM01-like Myoviridae that are typically documented to occur in western Lake Erie. This study provides further evidence that viral activity in western Lake Erie plays a significant role in transformation of microcystins from the particulate to the dissolved fraction and therefore requires monitoring efforts from local water treatment plants. Additionally, identification of multiple lytic cyanophages will enable the development of a quantitative PCR toolbox to assess viral activity during cHABs.IMPORTANCE Viral attack on cHABs may contribute to changes in community composition during blooms, as well as bloom decline, yet loss of bloom biomass does not eliminate the threat of cHAB toxicity. Rather, it may increase risks to the public by delivering a pool of dissolved toxin directly into water treatment utilities when the dominating Microcystis spp. are capable of producing microcystins. Detecting, characterizing, and quantifying the major cyanophages involved in lytic events will assist water treatment plant operators in making rapid decisions regarding the pool of microcystins entering the plant and the corresponding best practices to neutralize the toxin.

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.

Predetermined clockwork microbial worlds: Current understanding of aquatic microbial diel response from model systems to complex environments

In the photic zone of aquatic ecosystems, microorganisms with different metabolisms and their viruses form complex interactions and food webs. Within these interactions, phototrophic microorganisms such as eukaryotic microalgae and cyanobacteria interact directly with sunlight, and thereby generate circadian rhythms. Diel cycling originally generated in microbial phototrophs is directly transmitted toward heterotrophic microorganisms utilizing the photosynthetic products as they are excreted or exuded. Such diel cycling seems to be indirectly propagated toward heterotrophs as a result of complex biotic interactions. For example, cell death of phototrophic microorganisms induced by viral lysis and protistan grazing provides additional resources of dissolved organic matter to the microbial community, and so generates diel cycling in other heterotrophs with different nutrient dependencies. Likewise, differences in the diel transmitting pathway via complex interactions among heterotrophs, and between heterotrophs and their viruses, may also generate higher variation and time lag diel rhythms in different heterotrophic taxa. Thus, sunlight and photosynthesis not only contribute energy and carbon supply, but also directly or indirectly control diel cycling of the microbial community through complex interactions in the photic zone of aquatic ecosystems.

Genomic Analysis of Mic1 Reveals a Novel Freshwater Long-Tailed Cyanophage

Lake Chaohu, one of the five largest freshwater lakes in China, has been suffering from severe cyanobacterial blooms in the summer for many years. Cyanophages, the viruses that specifically infect cyanobacteria, play a key role in modulating cyanobacterial population, and thus regulate the emergence and decline of cyanobacterial blooms. Here we report a long-tailed cyanophage isolated from Lake Chaohu, termed Mic1, which specifically infects the cyanobacterium Microcystis aeruginosa. Mic1 has an icosahedral head of 88 nm in diameter and a long flexible tail of 400 nm. It possesses a circular genome of 92,627 bp, which contains 98 putative open reading frames. Genome sequence analysis enabled us to define a novel terminase large subunit that consists of two types of intein, indicating that the genome packaging of Mic1 is under fine control via posttranslational maturation of the terminase. Moreover, phylogenetic analysis suggested Mic1 and mitochondria share a common evolutionary origin of DNA polymerase γ gene. All together, these findings provided a start-point for investigating the co-evolution of cyanophages and its cyanobacterial hosts.