Cylindrospermopsis raciborskii Virus and host: genomic characterization and ecological relevance

Tradeoffs between phage resistance and nitrogen fixation drive the evolution of genes essential for cyanobacterial heterocyst functionality

Harmful blooms caused by diazotrophic (nitrogen-fixing) Cyanobacteria are becoming increasingly frequent and negatively impact aquatic environments worldwide. Cyanophages (viruses infecting Cyanobacteria) can potentially regulate cyanobacterial blooms, yet Cyanobacteria can rapidly acquire mutations that provide protection against phage infection. Here, we provide novel insights into cyanophage:Cyanobacteria interactions by characterizing the resistance to phages in two species of diazotrophic Cyanobacteria: Nostoc sp. and Cylindrospermopsis raciborskii. Our results demonstrate that phage resistance is associated with a fitness tradeoff by which resistant Cyanobacteria have reduced ability to fix nitrogen and/or to survive nitrogen starvation. Furthermore, we use whole-genome sequence analysis of 58 Nostoc-resistant strains to identify several mutations associated with phage resistance, including in cell surface-related genes and regulatory genes involved in the development and function of heterocysts (cells specialized in nitrogen fixation). Finally, we employ phylogenetic analyses to show that most of these resistance genes are accessory genes whose evolution is impacted by lateral gene transfer events. Together, these results further our understanding of the interplay between diazotrophic Cyanobacteria and their phages and suggest that a tradeoff between phage resistance and nitrogen fixation affects the evolution of cell surface-related genes and of genes involved in heterocyst differentiation and nitrogen fixation.

Sequencing the genomes of LPP-1, the first isolated cyanophage, and its relative LPP-2 reveal different integration mechanisms in closely related phages

In the early 1960s, the first cyanophage was isolated. The description of this phage, named LPP-1, led to the extensive investigation of various cyanophages and to the study of their interactions with their cyanobacterial hosts towards controlling blooms. Here, the genomes of LPP-1 and its putative relative, LPP-2 were sequenced. Sequencing these genomes revealed that LPP-1 and LPP-2 are members of a group of short-tailed cyanophages, which are distantly related to the T7-like cyanophages. Most of the phages in this group have the ability to lysogenize their hosts. Their ability to switch between lytic and lysogenic infection may explain the formation of cyanobacterial blooms despite the persistence of their phages. This lysogenic capacity of the LPP-1-like phages occurs despite the lack of an obvious integrase gene within their genomes. Interestingly, we show that LPP-2 integrates into the host genome through an integration site in high proximity to a recombination endonuclease that may have integrase activity. Further understanding of cyanobacterial-phage relationships may provide insight into their population dynamics and suggest novel approaches for control of destructive cyanobacterial blooms.

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.

Cyanophage infections reduce photosynthetic activity and expression of CO2 fixation genes in the freshwater bloom-forming cyanobacterium Aphanizomenon flos-aquae

Cyanobacteria play a significant role in ecosystem functioning as photosynthetic and CO₂ fixing microorganisms. Whether and to what extent cyanophages alter these carbon and energy cycles in their cyanobacterial hosts is still poorly understood. In this study, we investigated changes in photosynthetic activity (PSII), expression of genes associated with the light phase of photosynthesis (psbA, petA, ndhK) and carbon metabolism (rbcL, zwf) as well as intracellular ATP and NADHP concentrations in freshwater bloom-forming filamentous cyanobacterium Aphanizomenon flos-aquae infected by cyanophage vB_AphaS-CL131. We found that PSII activity and expression level of rbcL genes, indicating potential for CO₂ fixation, had decreased in response to cyanophage adsorption and DNA injection. During the period of viral DNA replication and assembly, PSII performance and gene expression remained at this decreased level and did not change significantly, indicating lack of transcriptional shutdown by the cyanophage. Combined, these observations suggest that although there is little to no interference between cyanophage DNA replication, host transcription and cellular metabolism, A. flos-aquae underwent a physiological state-shift toward lower efficiency of carbon and energy cycling. This further suggest potential cascading effect for co-occurring non-infected members of the microbial community.

Isolation and characterization of a novel Lambda-like phage infecting the bloom-forming cyanobacteria Cylindrospermopsis raciborskii

Cylindrospermopsis raciborskii is a central bloom‐forming cyanobacteria. However, despite its ecological significance, little is known of its interactions with the phages that infect it. Currently, only a single sequenced genome of a Cylindrospermopsis‐infecting phage is publicly available. Here we describe the isolation and characterization of Cr‐LKS3, a second phage infecting Cylindrospermopsis. Cr‐LKS3 is a siphovirus with a higher genome similarity to prophages within heterotrophic bacteria genomes than to any other cyanophage/cyano‐prophage, suggesting that it represents a novel cyanophage group. The function, order and orientation of the 72 genes in the Cr‐LKS3 genome are highly similar to those of Escherichia virus Lambda (hereafter Lambda), despite the very low sequence similarity between these phages, showing high evolutionary convergence despite the substantial difference in host characteristics. Similarly to Lambda, the genome of Cr‐LKS3 contains various genes that are known to be central to lysogeny, suggesting it can enter a lysogenic cycle. Cr‐LKS3 has a unique ability to infect a host with a dramatically different GC content, without carrying any tRNA genes to compensate for this difference. This ability, together with its potential lysogenic lifestyle shed light on the complex interactions between C. raciborskii and its phages.

The combined effect of clethodim (herbicide) and nitrogen variation on allelopathic interactions between Microcystis aeruginosa and Raphidiopsis raciborskii

The large-scale use of herbicides deteriorates water quality and threatens aquatic biodiversity. Unfortunately, there are few studies on the ecological effects of herbicides on toxin-producing strains of cyanobacteria under changing nutrient conditions. The objective of the present study was to investigate the effects of the herbicide clethodim and nitrogen variation on the allelopathic interactions and toxin production of Microcystis aeruginosa BCCUSP232 and Raphidiopsis raciborskii (formerly known as Cylindrospermopsis raciborskii) ITEPA1. M. aeruginosa had increased cell density when exposed to the clethodim (H +) (23.55 mg/L), whereas the highest cell density of R. raciborskii was observed in the treatment with clethodim plus limited nitrogen. Also, the cell-free exudate of R. raciborskii significantly stimulated the growth of M. aeruginosa on day 3 of the experiment. The concentration of chlorophyll-a in M. aeruginosa cultures generally increased in all the treatments, while in R. raciborskii cultures, the opposite occurred. Total microcystins (MCs) content of M. aeruginosa in the mixed cultures was 68% higher in nitrogen-enriched conditions than the control. A similar increase in MC content occurred in M. aeruginosa unialgal culture treated with R. raciborskii exudate. Total saxitoxin concentration was 81% higher in mixed cultures of R. raciborskii simultaneously exposed to high nitrogen and clethodim. Similarly, unialgal cultures of R. raciborskii exposed to either high nitrogen or clethodim had higher saxitoxins concentrations than the control. The intracellular H₂O₂ content of M. aeruginosa cultures decreased, whereas, in R. raciborskii cultures, it increased during exposure to high nitrogen and clethodim. Only R. raciborskii had a significant variation in peroxidase activity. The activities of glutathione S-transferase of both strains were higher in the presence of clethodim. These results revealed that nitrogen enrichment and the presence of clethodim might lead to the excessive proliferation of M. aeruginosa and R. raciborskii and increased production of cyanotoxins in aquatic environments.

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.

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 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.

The Landscape of Genetic Content in the Gut and Oral Human Microbiome

Despite substantial interest in the species diversity of the human microbiome and its role in disease, the scale of its genetic diversity, which is fundamental to deciphering human-microbe interactions, has not been quantified. Here, we conducted a cross-study meta-analysis of metagenomes from two human body niches, the mouth and gut, covering 3,655 samples from 13 studies. We found staggering genetic heterogeneity in the dataset, identifying a total of 45,666,334 non-redundant genes (23,961,508 oral and 22,254,436 gut) at the 95% identity level. Fifty percent of all genes were "singletons," or unique to a single metagenomic sample. Singletons were enriched for different functions (compared with non-singletons) and arose from sub-population-specific microbial strains. Overall, these results provide potential bases for the unexplained heterogeneity observed in microbiome-derived human phenotypes. One the basis of these data, we built a resource, which can be accessed at https://microbial-genes.bio.