Functional Genomic Analyses Reveal an Open Pan-genome for the Chloroviruses and a Potential for Genetic Innovation in New Isolates.

Chloroviruses (family Phycodnaviridae) are large double-stranded DNA (dsDNA) viruses that infect unicellular green algae present in inland waters. These viruses have been isolated using three main chlorella-like green algal host cells, traditionally called NC64A, SAG, and Pbi, revealing extensive genetic diversity. In this study, we performed a functional genomic analysis on 36 chloroviruses that infected the three different hosts. Phylogenetic reconstruction based on the DNA polymerase B family gene clustered the chloroviruses into three distinct clades. The viral pan-genome consists of 1,345 clusters of orthologous groups of genes (COGs), with 126 COGs conserved in all viruses. Totals of 368, 268, and 265 COGs are found exclusively in viruses that infect NC64A, SAG, and Pbi algal hosts, respectively. Two-thirds of the COGs have no known function, constituting the "dark pan-genome" of chloroviruses, and further studies focusing on these genes may identify important novelties. The proportions of functionally characterized COGs composing the pan-genome and the core-genome are similar, but those related to transcription and RNA processing, protein metabolism, and virion morphogenesis are at least 4-fold more represented in the core genome. Bipartite network construction evidencing the COG sharing among host-specific viruses identified 270 COGs shared by at least one virus from each of the different host groups. Finally, our results reveal an open pan-genome for chloroviruses and a well-established core genome, indicating that the isolation of new chloroviruses can be a valuable source of genetic discovery. IMPORTANCE Chloroviruses are large dsDNA viruses that infect unicellular green algae distributed worldwide in freshwater environments. They comprise a genetically diverse group of viruses; however, a comprehensive investigation of the genomic evolution of these viruses is still missing. Here, we performed a functional pan-genome analysis comprising 36 chloroviruses associated with three different algal hosts in the family Chlorellaceae, referred to as zoochlorellae because of their endosymbiotic lifestyle. We identified a set of 126 highly conserved genes, most of which are related to essential functions in the viral replicative cycle. Several genes are unique to distinct isolates, resulting in an open pan-genome for chloroviruses. This profile is associated with generalist organisms, and new insights into the evolution and ecology of chloroviruses are presented. Ultimately, our results highlight the potential for genetic diversity in new isolates.

Chlorovirus PBCV-1 protein A064R has three of the transferase activities necessary to synthesize its capsid protein N-linked glycans.

Paramecium bursaria chlorella virus-1 (PBCV-1) is a large double-stranded DNA (dsDNA) virus that infects the unicellular green alga Chlorella variabilis NC64A. Unlike many other viruses, PBCV-1 encodes most, if not all, of the enzymes involved in the synthesis of the glycans attached to its major capsid protein. Importantly, these glycans differ from those reported from the three domains of life in terms of structure and asparagine location in the sequon of the protein. Previous data collected from 20 PBCV-1 spontaneous mutants (or antigenic variants) suggested that the a064r gene encodes a glycosyltransferase (GT) with three domains, each with a different function. Here, we demonstrate that: domain 1 is a ß-l-rhamnosyltransferase; domain 2 is an α-l-rhamnosyltransferase resembling only bacterial proteins of unknown function, and domain 3 is a methyltransferase that methylates the C-2 hydroxyl group of the terminal α-l-rhamnose (Rha) unit. We also establish that methylation of the C-3 hydroxyl group of the terminal α-l-Rha is achieved by another virus-encoded protein A061L, which requires an O-2 methylated substrate. This study, thus, identifies two of the glycosyltransferase activities involved in the synthesis of the N-glycan of the viral major capsid protein in PBCV-1 and establishes that a single protein A064R possesses the three activities needed to synthetize the 2-OMe-α-l-Rha-(1→2)-ß-l-Rha fragment. Remarkably, this fragment can be attached to any xylose unit.

N-glycans from Paramecium bursaria chlorella virus MA-1D: Re-evaluation of the oligosaccharide common core structure.

Paramecium bursaria chlorella virus MA-1D is a chlorovirus that infects Chlorella variabilis strain NC64A, a symbiont of the protozoan Paramecium bursaria. MA-1D has a 339-kb genome encoding ca. 366 proteins and 11 tRNAs. Like other chloroviruses, its major capsid protein (MCP) is decorated with N-glycans, whose structures have been solved in this work by using nuclear magnetic spectroscopy and matrix-assisted laser desorption ionization-time of flight mass spectrometry along with MS/MS experiments. This analysis identified three N-linked oligosaccharides that differ in the nonstoichiometric presence of three monosaccharides, with the largest oligosaccharide composed of eight residues organized in a highly branched fashion. The N-glycans described here share several features with those of the other chloroviruses except that they lack a distal xylose unit that was believed to be part of a conserved core region for all the chloroviruses. Examination of the MA-1D genome detected a gene with strong homology to the putative xylosyltransferase in the reference chlorovirus PBCV-1 and in virus NY-2A, albeit mutated with a premature stop codon. This discovery means that we need to reconsider the essential features of the common core glycan region in the chloroviruses.

Structure of the chlorovirus PBCV-1 major capsid glycoprotein determined by combining crystallographic and carbohydrate molecular modeling approaches.

The glycans of the major capsid protein (Vp54) of Paramecium bursaria chlorella virus (PBCV-1) were recently described and found to be unusual. This prompted a reexamination of the previously reported Vp54 X-ray structure. A detailed description of the complete glycoprotein was achieved by combining crystallographic data with molecular modeling. The crystallographic data identified most of the monosaccharides located close to the protein backbone, but failed to detect those further from the glycosylation sites. Molecular modeling complemented this model by adding the missing monosaccharides and examined the conformational preference of the whole molecule, alone or within the crystallographic environment. Thus, combining X-ray crystallography with carbohydrate molecular modeling resulted in determining the complete glycosylated structure of a glycoprotein. In this case, it is the chlorovirus PBCV-1 major capsid protein.

Predators catalyze an increase in chloroviruses by foraging on the symbiotic hosts of zoochlorellae.

Virus population growth depends on contacts between viruses and their hosts. It is often unclear how sufficient contacts are made between viruses and their specific hosts to generate spikes in viral abundance. Here, we show that copepods, acting as predators, can bring aquatic viruses and their algal hosts into contact. Specifically, predation of the protist Paramecium bursaria by copepods resulted in a >100-fold increase in the number of chloroviruses in 1 d. Copepod predation can be seen as an ecological "catalyst" by increasing contacts between chloroviruses and their hosts, zoochlorellae (endosymbiotic algae that live within paramecia), thereby facilitating viral population growth. When feeding, copepods passed P. bursaria through their digestive tract only partially digested, releasing endosymbiotic algae that still supported viral reproduction and resulting in a virus population spike. A simple predator-prey model parameterized for copepods consuming protists generates cycle periods for viruses consistent with those observed in natural ponds. Food webs are replete with similar symbiotic organisms, and we suspect the predator catalyst mechanism is capable of generating blooms for other endosymbiont-targeting viruses.

Size-dependent Catalysis of Chlorovirus Population Growth by A Messy Feeding Predator.

Many chloroviruses replicate in endosymbiotic zoochlorellae that are protected from infection by their symbiotic host. To reach the high virus concentrations that often occur in natural systems, a mechanism is needed to release zoochlorellae from their hosts. We demonstrate that the ciliate predator Didinium nasutum foraging on zoochlorellae-bearing Paramecium bursaria can release live zoochlorellae from the ruptured prey cell that can then be infected by chloroviruses. The catalysis process is very effective, yielding roughly 95% of the theoretical infectious virus yield as determined by sonication of P. bursaria. Chlorovirus activation is more effective with smaller Didinia, as larger Didinia typically consume entire P. bursaria cells without rupturing them, precluding the release of zoochlorellae. We also show that the timing of Chlorovirus growth is tightly linked to the predator-prey cycle between Didinium and Paramecium, with the most rapid increase in chloroviruses temporally linked to the peak foraging rate of Didinium, supporting the idea that predator-prey cycles can drive cycles of Chlorovirus abundance.

The Autonomous Glycosylation of Large DNA Viruses.

Glycosylation of surface molecules is a key feature of several eukaryotic viruses, which use the host endoplasmic reticulum/Golgi apparatus to add carbohydrates to their nascent glycoproteins. In recent years, a newly discovered group of eukaryotic viruses, belonging to the Nucleo-Cytoplasmic Large DNA Virus (NCLDV) group, was shown to have several features that are typical of cellular organisms, including the presence of components of the glycosylation machinery. Starting from initial observations with the chlorovirus PBCV-1, enzymes for glycan biosynthesis have been later identified in other viruses; in particular in members of the Mimiviridae family. They include both the glycosyltransferases and other carbohydrate-modifying enzymes and the pathways for the biosynthesis of the rare monosaccharides that are found in the viral glycan structures. These findings, together with genome analysis of the newly-identified giant DNA viruses, indicate that the presence of glycogenes is widespread in several NCLDV families. The identification of autonomous viral glycosylation machinery leads to many questions about the origin of these pathways, the mechanisms of glycan production, and eventually their function in the viral replication cycle. The scope of this review is to highlight some of the recent results that have been obtained on the glycosylation systems of the large DNA viruses, with a special focus on the enzymes involved in nucleotide-sugar production.

Viral DNA Accumulation Regulates Replication Efficiency of Chlorovirus OSy-NE5 in Two Closely Related Chlorella variabilis Strains.

Many chloroviruses replicate in Chlorella variabilis algal strains that are ex-endosymbionts isolated from the protozoan Paramecium bursaria, including the NC64A and Syngen 2-3 strains. We noticed that indigenous water samples produced a higher number of plaque-forming viruses on C. variabilis Syngen 2-3 lawns than on C. variabilis NC64A lawns. These observed differences led to the discovery of viruses that replicate exclusively in Syngen 2-3 cells, named Only Syngen (OSy) viruses. Here, we demonstrate that OSy viruses initiate infection in the restricted host NC64A by synthesizing some early virus gene products and that approximately 20% of the cells produce a small number of empty virus capsids. However, the infected cells did not produce infectious viruses because the cells were unable to replicate the viral genome. This is interesting because all previous attempts to isolate host cells resistant to chlorovirus infection were due to changes in the host receptor for the virus.

Catalysis of Chlorovirus Production by the Foraging of Bursaria truncatella on Paramecia bursaria Containing Endosymbiotic Algae.

Chloroviruses are large viruses that replicate in chlorella-like green algae and normally exist as mutualistic endosymbionts (referred to as zoochlorellae) in protists such as Paramecium bursaria. Chlorovirus populations rise and fall in indigenous waters through time; however, the factors involved in these virus fluctuations are still under investigation. Chloroviruses attach to the surface of P. bursaria but cannot infect their zoochlorellae hosts because the viruses cannot reach the zoochlorellae as long as they are in the symbiotic phase. Predators of P. bursaria, such as copepods and didinia, can bring chloroviruses into contact with zoochlorellae by disrupting the paramecia, which results in an increase in virus titers in microcosm experiments. Here, we report that another predator of P. bursaria, Bursaria truncatella, can also increase chlorovirus titers. After two days of foraging on P. bursaria, B. truncatella increased infectious chlorovirus abundance about 20 times above the controls. Shorter term foraging (3 h) resulted in a small increase of chlorovirus titers over the controls and more foraging generated more chloroviruses. Considering that B. truncatella does not release viable zoochlorellae either during foraging or through fecal pellets, where zoochlorellae could be infected by chlorovirus, we suggest a third pathway of predator virus catalysis. By engulfing the entire protist and digesting it slowly, virus replication can occur within the predator and some of the virus is passed out through a waste vacuole. These results provide additional support for the hypothesis that predators of P. bursaria are important drivers of chlorovirus population sizes and dynamics.

Chlorovirus PBCV-1 Multidomain Protein A111/114R Has Three Glycosyltransferase Functions Involved in the Synthesis of Atypical N-Glycans.

The structures of the four N-linked glycans from the prototype chlorovirus PBCV-1 major capsid protein do not resemble any other glycans in the three domains of life. All known chloroviruses and antigenic variants (or mutants) share a unique conserved central glycan core consisting of five sugars, except for antigenic mutant virus P1L6, which has four of the five sugars. A combination of genetic and structural analyses indicates that the protein coded by PBCV-1 gene a111/114r, conserved in all chloroviruses, is a glycosyltransferase with three putative domains of approximately 300 amino acids each. Here, in addition to in silico sequence analysis and protein modeling, we measured the hydrolytic activity of protein A111/114R. The results suggest that domain 1 is a galactosyltransferase, domain 2 is a xylosyltransferase and domain 3 is a fucosyltransferase. Thus, A111/114R is the protein likely responsible for the attachment of three of the five conserved residues of the core region of this complex glycan, and, if biochemically corroborated, it would be the second three-domain protein coded by PBCV-1 that is involved in glycan synthesis. Importantly, these findings provide additional support that the chloroviruses do not use the canonical host endoplasmic reticulum-Golgi glycosylation pathway to glycosylate their glycoproteins; instead, they perform glycosylation independent of cellular organelles using virus-encoded enzymes.

Genetic Diversity of Potassium Ion Channel Proteins Encoded by Chloroviruses That Infect Chlorella heliozoae.

Chloroviruses are large, plaque-forming, dsDNA viruses that infect chlorella-like green algae that live in a symbiotic relationship with protists. Chloroviruses have genomes from 290 to 370 kb, and they encode as many as 400 proteins. One interesting feature of chloroviruses is that they encode a potassium ion (K+) channel protein named Kcv. The Kcv protein encoded by SAG chlorovirus ATCV-1 is one of the smallest known functional K+ channel proteins consisting of 82 amino acids. The KcvATCV-1 protein has similarities to the family of two transmembrane domain K+ channel proteins; it consists of two transmembrane α-helixes with a pore region in the middle, making it an ideal model for studying K+ channels. To assess their genetic diversity, kcv genes were sequenced from 103 geographically distinct SAG chlorovirus isolates. Of the 103 kcv genes, there were 42 unique DNA sequences that translated into 26 new Kcv channels. The new predicted Kcv proteins differed from KcvATCV-1 by 1 to 55 amino acids. The most conserved region of the Kcv protein was the filter, the turret and the pore helix were fairly well conserved, and the outer and the inner transmembrane domains of the protein were the most variable. Two of the new predicted channels were shown to be functional K+ channels.

Diversity of tRNA Clusters in the Chloroviruses.

Viruses rely on their host's translation machinery for the synthesis of their own proteins. Problems belie viral translation when the host has a codon usage bias (CUB) that is different from an infecting virus due to differences in the GC content between the host and virus genomes. Here, we examine the hypothesis that chloroviruses adapted to host CUB by acquisition and selection of tRNAs that at least partially favor their own CUB. The genomes of 41 chloroviruses comprising three clades, each infecting a different algal host, have been sequenced, assembled and annotated. All 41 viruses not only encode tRNAs, but their tRNA genes are located in clusters. While differences were observed between clades and even within clades, seven tRNA genes were common to all three clades of chloroviruses, including the tRNAArg gene, which was found in all 41 chloroviruses. By comparing the codon usage of one chlorovirus algal host, in which the genome has been sequenced and annotated (67% GC content), to that of two of its viruses (40% GC content), we found that the viruses were able to at least partially overcome the host's CUB by encoding tRNAs that recognize AU-rich codons. Evidence presented herein supports the hypothesis that a chlorovirus tRNA cluster was present in the most recent common ancestor (MRCA) prior to divergence into three clades. In addition, the MRCA encoded a putative isoleucine lysidine synthase (TilS) that remains in 39/41 chloroviruses examined herein, suggesting a strong evolutionary pressure to retain the gene. TilS alters the anticodon of tRNAMet that normally recognizes AUG to then recognize AUA, a codon for isoleucine. This is advantageous to the chloroviruses because the AUA codon is 12-13 times more common in the chloroviruses than their host, further helping the chloroviruses to overcome CUB. Among large DNA viruses infecting eukaryotes, the presence of tRNA genes and tRNA clusters appear to be most common in the Phycodnaviridae and, to a lesser extent, in the Mimiviridae.

SMRT Sequencing of Paramecium Bursaria Chlorella Virus-1 Reveals Diverse Methylation Stability in Adenines Targeted by Restriction Modification Systems.

Chloroviruses (family Phycodnaviridae) infect eukaryotic, freshwater, unicellular green algae. A unique feature of these viruses is an abundance of DNA methyltransferases, with isolates dedicating up to 4.5% of their protein coding potential to these genes. This diversity highlights just one of the long-standing values of the chlorovirus model system; where group-wide epigenomic characterization might begin to elucidate the function(s) of DNA methylation in large dsDNA viruses. We characterized DNA modifications in the prototype chlorovirus, PBCV-1, using single-molecule real time (SMRT) sequencing (aka PacBio). Results were compared to total available sites predicted in silico based on DNA sequence alone. SMRT-software detected N6-methyl-adenine (m6A) at GATC and CATG recognition sites, motifs previously shown to be targeted by PBCV-1 DNA methyltransferases M.CviAI and M. CviAII, respectively. At the same time, PacBio analyses indicated that 10.9% of the PBCV-1 genome had large interpulse duration ratio (ipdRatio) values, the primary metric for DNA modification identification. These events represent 20.6x more sites than can be accounted for by all available adenines in GATC and CATG motifs, suggesting base or backbone modifications other than methylation might be present. To define methylation stability, we cross-compared methylation status of each GATC and CATG sequence in three biological replicates and found ∼81% of sites were stably methylated, while ∼2% consistently lack methylation. The remaining 17% of sites were stochastically methylated. When methylation status was analyzed for both strands of each target, we show that palindromes existed in completely non-methylated states, fully-methylated states, or hemi-methylated states, though GATC sites more often lack methylation than CATG sequences. Given that both sequences are targeted by not just methyltransferases, but by restriction endonucleases that are together encoded by PBCV-1 as virus-originating restriction modification (RM) systems, there is strong selective pressure to modify all target sites. The finding that most instances of non-methylation are associated with hemi-methylation is congruent with observations that hemi-methylated palindromes are resistant to cleavage by restriction endonucleases. However, sites where hemi-methylation is conserved might represent a unique regulatory function for PBCV-1. This study serves as a baseline for future investigation into the epigenomics of chloroviruses and their giant virus relatives.

Chloroviruses.

Chloroviruses are large dsDNA, plaque-forming viruses that infect certain chlorella-like green algae; the algae are normally mutualistic endosymbionts of protists and metazoans and are often referred to as zoochlorellae. The viruses are ubiquitous in inland aqueous environments throughout the world and occasionally single types reach titers of thousands of plaque-forming units per ml of native water. The viruses are icosahedral in shape with a spike structure located at one of the vertices. They contain an internal membrane that is required for infectivity. The viral genomes are 290 to 370 kb in size, which encode up to 16 tRNAs and 330 to ~415 proteins, including many not previously seen in viruses. Examples include genes encoding DNA restriction and modification enzymes, hyaluronan and chitin biosynthetic enzymes, polyamine biosynthetic enzymes, ion channel and transport proteins, and enzymes involved in the glycan synthesis of the virus major capsid glycoproteins. The proteins encoded by many of these viruses are often the smallest or among the smallest proteins of their class. Consequently, some of the viral proteins are the subject of intensive biochemical and structural investigation.

Chitin synthesis by Chlorella cells infected by chloroviruses: Enhancement by adopting a slow-growing virus and treatment with aphidicolin.

Chlorella viruses or chloroviruses contain a gene that encodes an enzyme that catalyzes chitin synthesis. This gene is expressed early in viral infections to produce chitin on the outside of the Chlorella cell wall. Interestingly, chitin synthesis by microalgal Chlorella cells in combination with chloroviruses represents a unique eco-friendly process for converting solar energy and CO2 into useful materials. However, during the final viral infection stage, the host cells are completely lysed, so chitin should be harvested before cells lyse. To increase chitin yields, slow-growing chlorovirus isolates were adopted and the viral replication process was modified with an inhibitor of DNA synthesis. The accumulation of chitin on the surface of Chlorella cells infected with one of nine chlorovirus isolates carrying the chitin synthase gene was compared with that of CVK2 (a standard virus)-infected cells. Chlorella cells infected with CVNF-1 (a slow-growing virus) accumulated chitin over the entire cell surface within 15 min post-infection (p.i.), and chitin continued to accumulate for up to 8 h p.i. before cells lysed. This was 2-fold longer than the chitin-accumulation period for cells infected with CVK2. The addition of aphidicolin delayed the progression of the virus replication cycle and extended the chitin-accumulation period of CVNF-1-infected cells to 12 h p.i. before cells lysed. Additionally, chitin production in the aphidicolin-treated CVNF-1-infected cells was approximately 6-fold higher than in CVK2-infected cells not treated with aphidicolin. Thus, chitin synthesis in a Chlorella-virus system may be prolonged by using slow-growing viral isolates treated with aphidicolin.

Characterization of a new chlorovirus type with permissive and non-permissive features on phylogenetically related algal strains.

A previous report indicated that prototype chlorovirus PBCV-1 replicated in two Chlorella variabilis algal strains, NC64A and Syngen 2-3, that are ex-endosymbionts isolated from the protozoan Paramecium bursaria. Surprisingly, plaque-forming viruses on Syngen 2-3 lawns were often higher than on NC64A lawns from indigenous water samples. These differences led to the discovery of viruses that exclusively replicate in Syngen 2-3 cells, named Only Syngen (OSy) viruses. OSy-NE5, the prototype virus for the proposed new species, had a linear dsDNA genome of 327kb with 44-nucleotide-long, incompletely base-paired, covalently closed hairpin ends. Each hairpin structure was followed by an identical 2612 base-paired inverted sequence after which the DNA sequence diverged. OSy-NE5 encoded 357 predicted CDSs and 13 tRNAs. Interestingly, OSy-NE5 attached to and initiated infection in NC64A cells but infectious progeny viruses were not produced; thus OSy-NE5 replication in NC64A is blocked at some later stage of replication.

Chloroviruses Have a Sweet Tooth.

Chloroviruses are large double-stranded DNA (dsDNA) viruses that infect certain isolates of chlorella-like green algae. They contain up to approximately 400 protein-encoding genes and 16 transfer RNA (tRNA) genes. This review summarizes the unexpected finding that many of the chlorovirus genes encode proteins involved in manipulating carbohydrates. These include enzymes involved in making extracellular polysaccharides, such as hyaluronan and chitin, enzymes that make nucleotide sugars, such as GDP-L-fucose and GDP-D-rhamnose and enzymes involved in the synthesis of glycans attached to the virus major capsid proteins. This latter process differs from that of all other glycoprotein containing viruses that traditionally use the host endoplasmic reticulum and Golgi machinery to synthesize and transfer the glycans.

Large dsDNA chloroviruses encode diverse membrane transport proteins.

Many large DNA viruses that infect certain isolates of chlorella-like green algae (chloroviruses) are unusual because they often encode a diverse set of membrane transport proteins, including functional K(+) channels and aquaglyceroporins as well as K(+) transporters and calcium transporting ATPases. Some chloroviruses also encode putative ligand-gated-like channel proteins. No one protein is present in all of the chloroviruses that have been sequenced, but the K(+) channel is the most common as only two chloroviruses have been isolated that lack this complete protein. This review describes the properties of these membrane-transporting proteins and suggests possible physiological functions and evolutionary histories for some of them.

Dynamic attachment of Chlorovirus PBCV-1 to Chlorella variabilis.

Chloroviruses infect their hosts by specifically binding to and degrading the cell wall of their algal hosts at the site of attachment, using an intrinsic digesting enzyme(s). Chlorovirus PBCV-1 stored as a lysate survived longer than virus alone, suggesting virus attachment to cellular debris may be reversible. Ghost cells (algal cells extracted with methanol) were used as a model to study reversibility of PBCV-1 attachment because ghost cells are as susceptible to attachment and wall digestion as are live cells. Reversibility of attachment to ghost cells was examined by releasing attached virions with a cell wall degrading enzyme extract. The majority of the released virions retained infectivity even after re-incubating the released virions with ghost cells two times. Thus the chloroviruses appear to have a dynamic attachment strategy that may be beneficial in indigenous environments where cell wall debris can act as a refuge until appropriate host cells are available.