Genome of brown tide virus (AaV), the little giant of the Megaviridae, elucidates NCLDV genome expansion and host-virus coevolution

Kratosvirus quantuckense: the history and novelty of an algal bloom disrupting virus and a model for giant virus research

Since the discovery of the first "giant virus," particular attention has been paid toward isolating and culturing these large DNA viruses through Acanthamoeba spp. bait systems. While this method has allowed for the discovery of plenty novel viruses in the Nucleocytoviricota, environmental -omics-based analyses have shown that there is a wealth of diversity among this phylum, particularly in marine datasets. The prevalence of these viruses in metatranscriptomes points toward their ecological importance in nutrient turnover in our oceans and as such, in depth study into non-amoebal Nucleocytoviricota should be considered a focal point in viral ecology. In this review, we report on Kratosvirus quantuckense (née Aureococcus anophagefferens Virus), an algae-infecting virus of the Imitervirales. Current systems for study in the Nucleocytoviricota differ significantly from this virus and its relatives, and a litany of trade-offs within physiology, coding potential, and ecology compared to these other viruses reveal the importance of K. quantuckense. Herein, we review the research that has been performed on this virus as well as its potential as a model system for algal-virus interactions.

Genomic adaptation of giant viruses in polar oceans

Despite being perennially frigid, polar oceans form an ecosystem hosting high and unique biodiversity. Various organisms show different adaptive strategies in this habitat, but how viruses adapt to this environment is largely unknown. Viruses of phyla Nucleocytoviricota and Mirusviricota are groups of eukaryote-infecting large and giant DNA viruses with genomes encoding a variety of functions. Here, by leveraging the Global Ocean Eukaryotic Viral database, we investigate the biogeography and functional repertoire of these viruses at a global scale. We first confirm the existence of an ecological barrier that clearly separates polar and nonpolar viral communities, and then demonstrate that temperature drives dramatic changes in the virus-host network at the polar-nonpolar boundary. Ancestral niche reconstruction suggests that adaptation of these viruses to polar conditions has occurred repeatedly over the course of evolution, with polar-adapted viruses in the modern ocean being scattered across their phylogeny. Numerous viral genes are specifically associated with polar adaptation, although most of their homologues are not identified as polar-adaptive genes in eukaryotes. These results suggest that giant viruses adapt to cold environments by changing their functional repertoire, and this viral evolutionary strategy is distinct from the polar adaptation strategy of their hosts.

Giant Virus Infection Signatures Are Modulated by Euphotic Zone Depth Strata and Iron Regimes of the Subantarctic Southern Ocean

Viruses can alter the abundance, evolution, and metabolism of microorganisms in the ocean, playing a key role in water column biogeochemistry and global carbon cycles. Large efforts to measure the contribution of eukaryotic microorganisms (e.g., protists) to the marine food web have been made, yet the in situ activities of the ecologically relevant viruses that infect these organisms are not well characterized. Viruses within the phylum Nucleocytoviricota ("giant viruses") are known to infect a diverse range of ecologically relevant marine protists, yet how these viruses are influenced by environmental conditions remains under-characterized. By employing metatranscriptomic analyses of in situ microbial communities along a temporal and depth-resolved gradient, we describe the diversity of giant viruses at the Southern Ocean Time Series (SOTS), a site within the subpolar Southern Ocean. Using a phylogeny-guided taxonomic assessment of detected giant virus genomes and metagenome-assembled genomes, we observed depth-dependent structuring of divergent giant virus families mirroring dynamic physicochemical gradients in the stratified euphotic zone. Analyses of transcribed metabolic genes from giant viruses suggest viral metabolic reprogramming of hosts from the surface to a 200-m depth. Lastly, using on-deck incubations reflecting a gradient of iron availability, we show that modulating iron regimes influences the activity of giant viruses in the field. Specifically, we show enhanced infection signatures of giant viruses under both iron-replete and iron-limited conditions. Collectively, these results expand our understanding of how the water column's vertical biogeography and chemical surroundings affect an important group of viruses within the Southern Ocean. IMPORTANCE The biology and ecology of marine microbial eukaryotes is known to be constrained by oceanic conditions. In contrast, how viruses that infect this important group of organisms respond to environmental change is less well known, despite viruses being recognized as key microbial community members. Here, we address this gap in our understanding by characterizing the diversity and activity of "giant" viruses within an important region in the sub-Antarctic Southern Ocean. Giant viruses are double-stranded DNA (dsDNA) viruses of the phylum Nucleocytoviricota and are known to infect a wide range of eukaryotic hosts. By employing a metatranscriptomics approach using both in situ samples and microcosm manipulations, we illuminated both the vertical biogeography and how changing iron availability affects this primarily uncultivated group of protist-infecting viruses. These results serve as a foundation for our understanding of how the open ocean water column structures the viral community, which can be used to guide models of the viral impact on marine and global biogeochemical cycling.

Diversity and genomics of giant viruses in the North Pacific Subtropical Gyre

Large double-stranded DNA viruses of the phylum Nucleocytoviricota, often referred to as "giant viruses," are ubiquitous members of marine ecosystems that are important agents of mortality for eukaryotic plankton. Although giant viruses are known to be prevalent in marine systems, their activities in oligotrophic ocean waters remain unclear. Oligotrophic gyres constitute the majority of the ocean and assessing viral activities in these regions is therefore critical for understanding overall marine microbial processes. In this study, we generated 11 metagenome-assembled genomes (MAGs) of giant viruses from samples previously collected from Station ALOHA in the North Pacific Subtropical Gyre. Phylogenetic analyses revealed that they belong to the orders Imitervirales (n = 6), Algavirales (n = 4), and Pimascovirales (n = 1). Genome sizes ranged from ~119-574 kbp, and several of the genomes encoded predicted TCA cycle components, cytoskeletal proteins, collagen, rhodopsins, and proteins potentially involved in other cellular processes. Comparison with other marine metagenomes revealed that several have broad distribution across ocean basins and represent abundant viral constituents of pelagic surface waters. Our work sheds light on the diversity of giant viruses present in oligotrophic ocean waters across the globe.

Viral Complexity

Although traditionally viewed as streamlined and simple, discoveries over the last century have revealed that viruses can exhibit surprisingly complex physical structures, genomic organization, ecological interactions, and evolutionary histories. Viruses can have physical dimensions and genome lengths that exceed many cellular lineages, and their infection strategies can involve a remarkable level of physiological remodeling of their host cells. Virus–virus communication and widespread forms of hyperparasitism have been shown to be common in the virosphere, demonstrating that dynamic ecological interactions often shape their success. And the evolutionary histories of viruses are often fraught with complexities, with chimeric genomes including genes derived from numerous distinct sources or evolved de novo. Here we will discuss many aspects of this viral complexity, with particular emphasis on large DNA viruses, and provide an outlook for future research.

Closed, Circular Genome Sequence of Aureococcus anophagefferens Virus, a Lytic Virus of a Brown Tide-Forming Alga

Here, we report the genomic sequence of Aureococcus anophagefferens virus, assembled into one circular contig from both Nanopore and Illumina reads. The genome is 381,717 bp long with a GC content of 29.1%, which includes an additional 5-kb region between the previously predicted polar ends of the reference genome.

The Astounding World of Glycans from Giant Viruses

Viruses are a heterogeneous ensemble of entities, all sharing the need for a suitable host to replicate. They are extremely diverse, varying in morphology, size, nature, and complexity of their genomic content. Typically, viruses use host-encoded glycosyltransferases and glycosidases to add and remove sugar residues from their glycoproteins. Thus, the structure of the glycans on the viral proteins have, to date, typically been considered to mimick those of the host. However, the more recently discovered large and giant viruses differ from this paradigm. At least some of these viruses code for an (almost) autonomous glycosylation pathway. These viral genes include those that encode the production of activated sugars, glycosyltransferases, and other enzymes able to manipulate sugars at various levels. This review focuses on large and giant viruses that produce carbohydrate-processing enzymes. A brief description of those harboring these features at the genomic level will be discussed, followed by the achievements reached with regard to the elucidation of the glycan structures, the activity of the proteins able to manipulate sugars, and the organic synthesis of some of these virus-encoded glycans. During this progression, we will also comment on many of the challenging questions on this subject that remain to be addressed.

Metabolic arsenal of giant viruses: Host hijack or self-use?

Viruses generally are defined as lacking the fundamental properties of living organisms in that they do not harbor an energy metabolism system or protein synthesis machinery. However, the discovery of giant viruses of amoeba has fundamentally challenged this view because of their exceptional genome properties, particle sizes and encoding of the enzyme machinery for some steps of protein synthesis. Although giant viruses are not able to replicate autonomously and still require a host for their multiplication, numerous metabolic genes involved in energy production have been recently detected in giant virus genomes from many environments. These findings have further blurred the boundaries that separate viruses and living organisms. Herein, we summarize information concerning genes and proteins involved in cellular metabolic pathways and their orthologues that have, surprisingly, been discovered in giant viruses. The remarkable diversity of metabolic genes described in giant viruses include genes encoding enzymes involved in glycolysis, gluconeogenesis, tricarboxylic acid cycle, photosynthesis, and β-oxidation. These viral genes are thought to have been acquired from diverse biological sources through lateral gene transfer early in the evolution of Nucleo-Cytoplasmic Large DNA Viruses, or in some cases more recently. It was assumed that viruses are capable of hijacking host metabolic networks. But the giant virus auxiliary metabolic genes also may represent another form of host metabolism manipulation, by expanding the catalytic capabilities of the host cells especially in harsh environments, providing the infected host cells with a selective evolutionary advantage compared to non-infected cells and hence favoring the viral replication. However, the mechanism of these genes' functionality remains unclear to date.

Novel Viral Communities Potentially Assisting in Carbon, Nitrogen, and Sulfur Metabolism in the Upper Slope Sediments of Mariana Trench

Viruses are ubiquitous in the oceans. Even in the deep sediments of the Mariana Trench, viruses have high productivity. However, little is known about their species composition and survival strategies in that environment. Here, we uncovered novel viral communities (3,206 viral scaffolds) in the upper slope sediments of the Mariana Trench via metagenomic analysis of 15 sediment samples. Most (99%) of the viral scaffolds lack known viral homologs, and ca. 59% of the high-quality viral genomes (total of 111 with completeness of >90%) represent novel genera, including some Phycodnaviridae and jumbo phages. These viruses contain various auxiliary metabolic genes (AMGs) potentially involved in organic carbon degradation, inorganic carbon fixation, denitrification, and assimilatory sulfate reduction, etc. This study provides novel insight into the almost unknown benthic viral communities in the Mariana Trench. IMPORTANCE The Mariana Trench harbors a substantial number of infective viral particles. However, very little is known about the identity, survival strategy, and potential functions of viruses in the trench sediments. Here, through metagenomic analysis, unusual benthic viral communities with high diversity and novelty were discovered. Among them, 59% of the viruses with a genome completeness of >90% represent novel genera. Various auxiliary metabolic genes carried by these viruses reflect the potential adaptive characteristics of viruses in this extreme environment and the biogeochemical cycles that they may participate in. This study gives us a deeper understanding of the peculiarities of viral communities in deep-sea/hadal sediments.

Aureococcus anophagefferens (Pelagophyceae) genomes improve evaluation of nutrient acquisition strategies involved in brown tide dynamics

The pelagophyte Aureococcus anophagefferens causes harmful brown tide blooms in marine embayments on three continents. Aureococcus anophagefferens was the first harmful algal bloom species to have its genome sequenced, an advance that evidenced genes important for adaptation to environmental conditions that prevail during brown tides. To expand the genomic tools available for this species, genomes for four strains were assembled, including three newly sequenced strains and one assembled from publicly available data. These genomes ranged from 57.11 to 73.62 Mb, encoding 13,191-17,404 potential proteins. All strains shared ~90% of their encoded proteins as determined by homology searches and shared most functional orthologs as determined by KEGG, although each strain also possessed coding sequences with unique functions. Like the original reference genome, the genomes assembled in this study possessed genes hypothesized to be important in bloom proliferation, including genes involved in organic compound metabolism and growth at low light. Cross-strain informatics and culture experiments suggest that the utilization of purines is a potentially important source of organic nitrogen for brown tides. Analyses of metatranscriptomes from a brown tide event demonstrated that use of a single genome yielded a lower read mapping percentage (~30% of library reads) as compared to a database generated from all available genomes (~43%), suggesting novel information about bloom ecology can be gained from expanding genomic space. This work demonstrates the continued need to sequence ecologically relevant algae to understand the genomic potential and their ecology in the environment.

Year-round dynamics of amplicon sequence variant communities differ among eukaryotes, Imitervirales, and prokaryotes in a coastal ecosystem

Coastal microbial communities are affected by seasonal environmental change, biotic interactions, and fluctuating nutrient availability. We investigated the seasonal dynamics of communities of eukaryotes, a major group of double-stranded DNA viruses that infect eukaryotes (order Imitervirales; phylum Nucleocytoviricota), and prokaryotes in the Uranouchi Inlet, Kochi, Japan. We performed metabarcoding using ribosomal RNA genes and viral polB genes as markers in 43 seawater samples collected over 20 months. Eukaryotes, prokaryotes, and Imitervirales communities characterized by the compositions of amplicon sequence variants (ASVs) showed synchronic seasonal cycles. However, the community dynamics showed intriguing differences in several aspects, such as the recovery rate after a year. We also showed that the differences in community dynamics were at least partially explained by differences in recurrence/persistence levels of individual ASVs among eukaryotes, prokaryotes, and Imitervirales. Prokaryotic ASVs were the most persistent, followed by eukaryotic ASVs and Imitervirales ASVs, which were the least persistent. We argue that the differences in the specificity of interactions (virus-eukaryote vs. prokaryote-eukaryote) as well as the niche breadth of community members were at the origin of the distinct community dynamics among eukaryotes, their viruses, and prokaryotes.

A phylogenomic framework for charting the diversity and evolution of giant viruses

Large DNA viruses of the phylum Nucleocytoviricota have recently emerged as important members of ecosystems around the globe that challenge traditional views of viral complexity. Numerous members of this phylum that cannot be classified within established families have recently been reported, and there is presently a strong need for a robust phylogenomic and taxonomic framework for these viruses. Here, we report a comprehensive phylogenomic analysis of the Nucleocytoviricota, present a set of giant virus orthologous groups (GVOGs) together with a benchmarked reference phylogeny, and delineate a hierarchical taxonomy within this phylum. We show that the majority of Nucleocytoviricota diversity can be partitioned into 6 orders, 32 families, and 344 genera, substantially expanding the number of currently recognized taxonomic ranks for these viruses. We integrate our results within a taxonomy that has been adopted for all viruses to establish a unifying framework for the study of Nucleocytoviricota diversity, evolution, and environmental distribution.

Giant viruses have transformed our understanding of viral complexity, but we lack a framework for examining their diversity in the biosphere. This study presents a phylogenomic resource for charting the diversity, ecology, and evolution of giant viruses.

High Transcriptional Activity and Diverse Functional Repertoires of Hundreds of Giant Viruses in a Coastal Marine System

Viruses belonging to the Nucleocytoviricota phylum are globally distributed and include members with notably large genomes and complex functional repertoires. Recent studies have shown that these viruses are particularly diverse and abundant in marine systems, but the magnitude of actively replicating Nucleocytoviricota present in ocean habitats remains unclear. In this study, we compiled a curated database of 2,431 Nucleocytoviricota genomes and used it to examine the gene expression of these viruses in a 2.5-day metatranscriptomic time-series from surface waters of the California Current. We identified 145 viral genomes with high levels of gene expression, including 90 Imitervirales and 49 Algavirales viruses. In addition to recovering high expression of core genes involved in information processing that are commonly expressed during viral infection, we also identified transcripts of diverse viral metabolic genes from pathways such as glycolysis, the TCA cycle, and the pentose phosphate pathway, suggesting that virus-mediated reprogramming of central carbon metabolism is common in oceanic surface waters. Surprisingly, we also identified viral transcripts with homology to actin, myosin, and kinesin domains, suggesting that viruses may use these gene products to manipulate host cytoskeletal dynamics during infection. We performed phylogenetic analysis on the virus-encoded myosin and kinesin proteins, which demonstrated that most belong to deep-branching viral clades, but that others appear to have been acquired from eukaryotes more recently. Our results highlight a remarkable diversity of active Nucleocytoviricota in a coastal marine system and underscore the complex functional repertoires expressed by these viruses during infection. IMPORTANCE The discovery of giant viruses has transformed our understanding of viral complexity. Although viruses have traditionally been viewed as filterable infectious agents that lack metabolism, giant viruses can reach sizes rivalling cellular lineages and possess genomes encoding central metabolic processes. Recent studies have shown that giant viruses are widespread in aquatic systems, but the activity of these viruses and the extent to which they reprogram host physiology in situ remains unclear. Here, we show that numerous giant viruses consistently express central metabolic enzymes in a coastal marine system, including components of glycolysis, the TCA cycle, and other pathways involved in nutrient homeostasis. Moreover, we found expression of several viral-encoded actin, myosin, and kinesin genes, indicating viral manipulation of the host cytoskeleton during infection. Our study reveals a high activity of giant viruses in a coastal marine system and indicates they are a diverse and underappreciated component of microbial diversity in the ocean.

Metatranscriptome Library Preparation Influences Analyses of Viral Community Activity During a Brown Tide Bloom

There is growing interest in the use of metatranscriptomics to study virus community dynamics. We used RNA samples collected from harmful brown tides caused by the eukaryotic alga Aureococcus anophagefferens within New York (United States) estuaries and in the process observed how preprocessing of libraries by either selection for polyadenylation or reduction in ribosomal RNA (rRNA) influenced virus community analyses. As expected, more reads mapped to the A. anophagefferens genome in polyadenylation-selected libraries compared to the rRNA-reduced libraries, with reads mapped in each sample correlating to one another regardless of preprocessing of libraries. Yet, this trend was not seen for reads mapping to the Aureococcus anophagefferens Virus (AaV), where significantly more reads (approximately two orders of magnitude) were mapped to the AaV genome in the rRNA-reduced libraries. In the rRNA-reduced libraries, there was a strong and significant correlation between reads mappings to AaV and A. anophagefferens. Overall, polyadenylation-selected libraries produced fewer viral contigs, fewer reads mapped to viral contigs, and different proportions across viral realms and families, compared to their rRNA-reduced pairs. This study provides evidence that libraries generated by rRNA reduction and not selected for polyadenylation are more appropriate for quantitative characterization of viral communities in aquatic ecosystems by metatranscriptomics.

Quantitative Assessment of Nucleocytoplasmic Large DNA Virus and Host Interactions Predicted by Co-occurrence Analyses

Nucleocytoplasmic large DNA viruses (NCLDVs) are highly diverse and abundant in marine environments. However, the knowledge of their hosts is limited because only a few NCLDVs have been isolated so far. Taking advantage of the recent large-scale marine metagenomics census, in silico host prediction approaches are expected to fill the gap and further expand our knowledge of virus-host relationships for unknown NCLDVs. In this study, we built co-occurrence networks of NCLDVs and eukaryotic taxa to predict virus-host interactions using Tara Oceans sequencing data. Using the positive likelihood ratio to assess the performance of host prediction for NCLDVs, we benchmarked several co-occurrence approaches and demonstrated an increase in the odds ratio of predicting true positive relationships 4-fold compared to random host predictions. To further refine host predictions from high-dimensional co-occurrence networks, we developed a phylogeny-informed filtering method, Taxon Interaction Mapper, and showed it further improved the prediction performance by 12-fold. Finally, we inferred virophage-NCLDV networks to corroborate that co-occurrence approaches are effective for predicting interacting partners of NCLDVs in marine environments.IMPORTANCE NCLDVs can infect a wide range of eukaryotes, although their life cycle is less dependent on hosts compared to other viruses. However, our understanding of NCLDV-host systems is highly limited because few of these viruses have been isolated so far. Co-occurrence information has been assumed to be useful to predict virus-host interactions. In this study, we quantitatively show the effectiveness of co-occurrence inference for NCLDV host prediction. We also improve the prediction performance with a phylogeny-guided method, which leads to a concise list of candidate host lineages for three NCLDV families. Our results underpin the usage of co-occurrence approaches for the metagenomic exploration of the ecology of this diverse group of viruses.

Mitochondrial Genome Evolution in Pelagophyte Algae

The Pelagophyceae are marine stramenopile algae that include Aureoumbra lagunensis and Aureococcus anophagefferens, two microbial species notorious for causing harmful algal blooms. Despite their ecological significance, relatively few genomic studies of pelagophytes have been carried out. To improve understanding of the biology and evolution of pelagophyte algae, we sequenced complete mitochondrial genomes for A. lagunensis (CCMP1510), Pelagomonas calceolata (CCMP1756), and five strains of Aureoc. anophagefferens (CCMP1707, CCMP1708, CCMP1850, CCMP1984, and CCMP3368) using Nanopore long-read sequencing. All pelagophyte mitochondrial genomes assembled into single, circular mapping contigs between 39,376 bp (P. calceolata) and 55,968 bp (A. lagunensis) in size. Mitochondrial genomes for the five Aureoc. anophagefferens strains varied slightly in length (42,401–42,621 bp) and were 99.4–100.0% identical. Gene content and order were highly conserved between the Aureoc. anophagefferens and P. calceolata genomes, with the only major difference being a unique region in Aureoc. anophagefferens containingDNA adenine and cytosine methyltransferase (dam/dcm) genes that appear to be the product of lateral gene transfer from a prokaryotic or viral donor. Although the A. lagunensis mitochondrial genome shares seven distinct syntenic blocks with the other pelagophyte genomes, it has a tandem repeat expansion comprising ∼40% of its length, and lacks identifiable rps19 and glycine tRNA genes. Laterally acquired self-splicing introns were also found in the 23S rRNA (rnl) gene of P. calceolata and the coxI gene of the five Aureoc. anophagefferens genomes. Overall, these data provide baseline knowledge about the genetic diversity of bloom-forming pelagophytes relative to nonbloom-forming species.

Evolution of a major virion protein of the giant pandoraviruses from an inactivated bacterial glycoside hydrolase

The diverse viruses in the phylum Nucleocytoviricota (also known as NLCDVs, Nucleo-cytoplasmic Large DNA Viruses) typically possess large icosahedral virions. However, in several families of Nucleocytoviricota, the icosahedral capsid was replaced by irregular particle shapes, most notably, the amphora-like virions of pandoraviruses and pithoviruses, the largest known virus particles in the entire virosphere. Pandoraviruses appear to be the most highly derived viruses in this phylum because their evolution involved not only the change in the virion shape, but also, the actual loss of the gene encoding double-jelly roll major capsid protein (DJR MCP), the main building block of icosahedral capsids in this virus assemblage. Instead, pandoravirus virions are built of unrelated abundant proteins. Here we show that the second most abundant virion protein of pandoraviruses, major virion protein 2 (MVP2), evolved from an inactivated derivative of a bacterial glycoside hydrolase of the GH16 family. The ancestral form of MVP2 was apparently acquired early in the evolution of the Nucleocytoviricota, to become a minor virion protein. After a duplication in the common ancestor of pandoraviruses and molliviruses, one of the paralogs displaces DJR MCP in pandoraviruses, conceivably, opening the way for a major increase in the size of the virion and the genome. Exaptation of a carbohydrate-binding protein for the function of the MVP is a general trend in virus evolution and might underlie the transformation of the virion shape in other groups of the Nucleocytoviricota as well.

Fundamental Difficulties Prevent the Reconstruction of the Deep Phylogeny of Viruses

The extension of virology beyond its traditional medical, veterinary, or agricultural applications, now called environmental virology, has shown that viruses are both the most numerous and diverse biological entities on Earth. In particular, virus isolations from unicellular eukaryotic hosts (heterotrophic and photosynthetic protozoans) revealed numerous viral types previously unexpected in terms of virion structure, gene content, or mode of replication. Complemented by large-scale metagenomic analyses, these discoveries have rekindled interest in the enigma of the origin of viruses, for which a description encompassing all their diversity remains not available. Several laboratories have repeatedly tackled the deep reconstruction of the evolutionary history of viruses, using various methods of molecular phylogeny applied to the few shared "core" genes detected in certain virus groups (e.g., the Nucleocytoviricota). Beyond the practical difficulties of establishing reliable homology relationships from extremely divergent sequences, I present here conceptual arguments highlighting several fundamental limitations plaguing the reconstruction of the deep evolutionary history of viruses, and even more the identification of their unique or multiple origin(s). These arguments also underline the risk of establishing premature high level viral taxonomic classifications. Those limitations are direct consequences of the random mechanisms governing the reductive/retrogressive evolution of all obligate intracellular parasites.

Current capsid assembly models of icosahedral nucleocytoviricota viruses

Nucleocytoviricota viruses (NCVs) belong to a newly established phylum originally grouped as Nucleocytoplasmic large DNA viruses. NCVs are unique because of their large and complicated genomes that contain cellular genes with homologs from all kingdoms of life, raising intensive debates on their evolutional origins. Many NCVs pack their genomes inside massive icosahedral capsids assembled from thousands of proteins. Studying the assembly mechanism of such capsids has been challenging until breakthroughs from structural studies. Subsequently, several models of the capsid assembly were proposed, which provided some interesting insights on this elaborate process. In this review, we discuss three of the most recent assembly models as well as supporting experimental observations. Furthermore, we propose a new model that combines research developments from multiple sources. Investigation of the assembly process of these vast NCV capsids will facilitate future deciphering of the molecular mechanisms driving the formation of similar supramolecular complexes.

A Functional K+ Channel from Tetraselmis Virus 1, a Member of the Mimiviridae

Potassium ion (K⁺) channels have been observed in diverse viruses that infect eukaryotic marine and freshwater algae. However, experimental evidence for functional K⁺ channels among these alga-infecting viruses has thus far been restricted to members of the family Phycodnaviridae, which are large, double-stranded DNA viruses within the phylum Nucleocytoviricota. Recent sequencing projects revealed that alga-infecting members of Mimiviridae, another family within this phylum, may also contain genes encoding K⁺ channels. Here we examine the structural features and the functional properties of putative K⁺ channels from four cultivated members of Mimiviridae. While all four proteins contain variations of the conserved selectivity filter sequence of K⁺ channels, structural prediction algorithms suggest that only two of them have the required number and position of two transmembrane domains that are present in all K⁺ channels. After in vitro translation and reconstitution of the four proteins in planar lipid bilayers, we confirmed that one of them, a 79 amino acid protein from the virus Tetraselmis virus 1 (TetV-1), forms a functional ion channel with a distinct selectivity for K⁺ over Na⁺ and a sensitivity to Ba²⁺. Thus, virus-encoded K⁺ channels are not limited to Phycodnaviridae but also occur in the members of Mimiviridae. The large sequence diversity among the viral K⁺ channels implies multiple events of lateral gene transfer.

Structural and Proteomic Studies of the Aureococcus anophagefferens Virus Demonstrate a Global Distribution of Virus-Encoded Carbohydrate Processing

Viruses modulate the function(s) of environmentally relevant microbial populations, yet considerations of the metabolic capabilities of individual virus particles themselves are rare. We used shotgun proteomics to quantitatively identify 43 virus-encoded proteins packaged within purified Aureococcus anophagefferens Virus (AaV) particles, normalizing data to the per-virion level using a 9.5-Å-resolution molecular reconstruction of the 1900-Å (AaV) particle that we generated with cryogenic electron microscopy. This packaged proteome was used to determine similarities and differences between members of different giant virus families. We noted that proteins involved in sugar degradation and binding (e.g., carbohydrate lyases) were unique to AaV among characterized giant viruses. To determine the extent to which this virally encoded metabolic capability was ecologically relevant, we examined the TARA Oceans dataset and identified genes and transcripts of viral origin. Our analyses demonstrated that putative giant virus carbohydrate lyases represented up to 17% of the marine pool for this function. In total, our observations suggest that the AaV particle has potential prepackaged metabolic capabilities and that these may be found in other giant viruses that are widespread and abundant in global oceans.

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.

Internal Nitrogen Pools Shape the Infection of Aureococcus anophagefferens CCMP 1984 by a Giant Virus

The pelagophyte Aureococcus anophagefferens blooms annually in shallow bays around the world, where it is hypothesized to outcompete other phytoplankton in part by using alternative nitrogen sources. The high proportion of natural populations that are infected during the late stages of the bloom suggest viruses cause bloom collapse. We hypothesized that the Aureococcus anophagefferens Virus (AaV) infection cycle would be negatively influenced in cultures acclimated to decreasing external nitrogen conditions, but that the real-time external nitrogen concentration would not influence the infection cycle. Cultures acclimated in NO 3 - concentrations (0.0147 mM; N:P = 0.1225) that showed reduced end point cell abundances, forward scatter (a proxy for size) and red fluorescence (a proxy for chlorophyll a), also produced fewer viruses per cell at a slower rate. Decreasing the external concentration of nitrogen post infection did not alter burst size or time to lysis. These data suggest that the nitrogen used for new viral progeny is present within host cells at the time of infection. Flow cytometric data of an infection cycle showed a reduction in red fluorescence around twelve hours post infection, consistent with degradation of nitrogen-rich chloroplasts during the infection cycle. Using cell and virus quota estimates, we determined that A. anophagefferens cells had sufficient nitrogen and carbon for the lower ranges of burst sizes determined but did not contain enough phosphorous. Consistent with this observation, expression of nitrate and sugar transporters did not increase in the publicly available transcriptome data of the infection cycle, while several phosphorus transporters were. Our data demonstrate that dynamics of viruses infecting Aureococcus over the course of a bloom is dictated by the host cell state upon infection, which is set a priori by external nutrient supplies.

Dynamic genome evolution and complex virocell metabolism of globally-distributed giant viruses

The discovery of eukaryotic giant viruses has transformed our understanding of the limits of viral complexity, but the extent of their encoded metabolic diversity remains unclear. Here we generate 501 metagenome-assembled genomes of Nucleo-Cytoplasmic Large DNA Viruses (NCLDV) from environments around the globe, and analyze their encoded functional capacity. We report a remarkable diversity of metabolic genes in widespread giant viruses, including many involved in nutrient uptake, light harvesting, and nitrogen metabolism. Surprisingly, numerous NCLDV encode the components of glycolysis and the TCA cycle, suggesting that they can re-program fundamental aspects of their host's central carbon metabolism. Our phylogenetic analysis of NCLDV metabolic genes and their cellular homologs reveals distinct clustering of viral sequences into divergent clades, indicating that these genes are virus-specific and were acquired in the distant past. Overall our findings reveal that giant viruses encode complex metabolic capabilities with evolutionary histories largely independent of cellular life, strongly implicating them as important drivers of global biogeochemical cycles.

An Optimized Metabarcoding Method for Mimiviridae

Mimiviridae is a group of viruses with large genomes and virions. Ecological relevance of Mimiviridae in marine environments has been increasingly recognized through the discoveries of novel isolates and metagenomic studies. To facilitate ecological profiling of Mimiviridae, we previously proposed a meta-barcoding approach based on 82 degenerate primer pairs (i.e., MEGAPRIMER) targeting the DNA polymerase gene of Mimiviridae. The method detected a larger number of operational taxonomic units (OTUs) in environmental samples than previous methods. However, it required large quantities of DNA and was laborious due to the use of individual primer pairs. Here, we examined coastal seawater samples using varying PCR conditions and purification protocols to streamline the MEGAPRIMER method. Mixing primer pairs in "cocktails" reduced the required amount of environmental DNA by 90%, while reproducing the results obtained by the original protocol. We compared the results obtained by the meta-barcoding approach with quantifications using qPCR for selected OTUs. This revealed possible amplification biases among different OTUs, but the frequency profiles for individual OTUs across multiple samples were similar to those obtained by qPCR. We anticipate that the newly developed MEGAPRIMER protocols will be useful for ecological investigation of Mimiviridae in a larger set of environmental samples.

Influence of light on the infection of Aureococcus anophagefferens CCMP 1984 by a "giant virus"

The pelagophyte Aureococcus anophagefferens has caused recurrent brown tide blooms along the northeast coast of the United States since the mid-1980's, and more recently spread to other regions of the globe. These blooms, due to the high cell densities, are associated with severe light attenuation that destroys the sea grass beds which provide the basis for many fisheries. Data collected by transmission electron microscopy, PCR, and metatranscriptomic studies of the blooms, support the hypothesis that large dsDNA viruses play a role in bloom dynamics. While a large (~140 nm) icosahedral virus, with a 371 kbp genome, was first isolated more than a decade ago, the constraints imposed by environmental parameters on bloom infection dynamics by Aureococcus anophagefferens Virus, (AaV) remain unknown. To investigate the role light plays in infection by this virus, we acclimated A. anophagefferens to light intensities of 30 (low), 60 (medium) or 90 μmol photons m-2 s-1 (high) and infected cultures at these irradiance levels. Moreover, we completed light shift experiments where acclimated cultures were exposed to even lower light intensities (0, 5, and 15 μmol photons m-2 s-1) consistent with irradiance found during the peak of the bloom when cell concentrations are highest. The abundance of viruses produced per lytic event (burst size) was lower in the low irradiance acclimated cultures compared to the medium and high acclimated cultures. Transferring infected cultures to more-limiting light availabilities further decreased burst size and increased the length of time it took for cultures to lyse, regardless of acclimation irradiance level. A hypothetical mechanism for the reduced efficiency of the infection cycle in low light due to ribosome biogenesis was predicted from pre-existing transcriptomes. Overall, these studies provide a framework for understanding light effects on infection dynamics over the course of the summer months when A. anophagefferens blooms occur.

Targeted metagenomic recovery of four divergent viruses reveals shared and distinctive characteristics of giant viruses of marine eukaryotes

Giant viruses have remarkable genomic repertoires-blurring the line with cellular life-and act as top-down controls of eukaryotic plankton. However, to date only six cultured giant virus genomes are available from the pelagic ocean. We used at-sea flow cytometry with staining and sorting designed to target wild predatory eukaryotes, followed by DNA sequencing and assembly, to recover novel giant viruses from the Pacific Ocean. We retrieved four 'PacV' partial genomes that range from 421 to 1605 Kb, with 13 contigs on average, including the largest marine viral genomic assembly reported to date. Phylogenetic analyses indicate that three of the new viruses span a clade with deep-branching members of giant Mimiviridae, incorporating the Cafeteria roenbergensis virus, the uncultivated terrestrial Faunusvirus, one PacV from a choanoflagellate and two PacV with unclear hosts. The fourth virus, oPacV-421, is phylogenetically related to viruses that infect haptophyte algae. About half the predicted proteins in each PacV have no matches in NCBI nr (e-value < 10-5), totalling 1735 previously unknown proteins; the closest affiliations of the other proteins were evenly distributed across eukaryotes, prokaryotes and viruses of eukaryotes. The PacVs encode many translational proteins and two encode eukaryotic-like proteins from the Rh family of the ammonium transporter superfamily, likely influencing the uptake of nitrogen during infection. cPacV-1605 encodes a microbial viral rhodopsin (VirR) and the biosynthesis pathway for the required chromophore, the second finding of a choanoflagellate-associated virus that encodes these genes. In co-collected metatranscriptomes, 85% of cPacV-1605 genes were expressed, with capsids, heat shock proteins and proteases among the most highly expressed. Based on orthologue presence-absence patterns across the PacVs and other eukaryotic viruses, we posit the observed viral groupings are connected to host lifestyles as heterotrophs or phototrophs. This article is part of a discussion meeting issue 'Single cell ecology'.

Metabolic and biogeochemical consequences of viral infection in aquatic ecosystems

Ecosystems are controlled by 'bottom-up' (resources) and 'top-down' (predation) forces. Viral infection is now recognized as a ubiquitous top-down control of microbial growth across ecosystems but, at the same time, cell death by viral predation influences, and is influenced by, resource availability. In this Review, we discuss recent advances in understanding the biogeochemical impact of viruses, focusing on how metabolic reprogramming of host cells during lytic viral infection alters the flow of energy and nutrients in aquatic ecosystems. Our synthesis revealed several emerging themes. First, viral infection transforms host metabolism, in part through virus-encoded metabolic genes; the functions performed by these genes appear to alleviate energetic and biosynthetic bottlenecks to viral production. Second, viral infection depends on the physiological state of the host cell and on environmental conditions, which are challenging to replicate in the laboratory. Last, metabolic reprogramming of infected cells and viral lysis alter nutrient cycling and carbon export in the oceans, although the net impacts remain uncertain. This Review highlights the need for understanding viral infection dynamics in realistic physiological and environmental contexts to better predict their biogeochemical consequences.

A distinct lineage of giant viruses brings a rhodopsin photosystem to unicellular marine predators

Giant viruses are remarkable for their large genomes, often rivaling those of small bacteria, and for having genes thought exclusive to cellular life. Most isolated to date infect nonmarine protists, leaving their strategies and prevalence in marine environments largely unknown. Using eukaryotic single-cell metagenomics in the Pacific, we discovered a Mimiviridae lineage of giant viruses, which infects choanoflagellates, widespread protistan predators related to metazoans. The ChoanoVirus genomes are the largest yet from pelagic ecosystems, with 442 of 862 predicted proteins lacking known homologs. They are enriched in enzymes for modifying organic compounds, including degradation of chitin, an abundant polysaccharide in oceans, and they encode 3 divergent type-1 rhodopsins (VirR) with distinct evolutionary histories from those that capture sunlight in cellular organisms. One (VirR_DTS) is similar to the only other putative rhodopsin from a virus (PgV) with a known host (a marine alga). Unlike the algal virus, ChoanoViruses encode the entire pigment biosynthesis pathway and cleavage enzyme for producing the required chromophore, retinal. We demonstrate that the rhodopsin shared by ChoanoViruses and PgV binds retinal and pumps protons. Moreover, our 1.65-Å resolved VirR_DTS crystal structure and mutational analyses exposed differences from previously characterized type-1 rhodopsins, all of which come from cellular organisms. Multiple VirR types are present in metagenomes from across surface oceans, where they are correlated with and nearly as abundant as a canonical marker gene from Mimiviridae Our findings indicate that light-dependent energy transfer systems are likely common components of giant viruses of photosynthetic and phagotrophic unicellular marine eukaryotes.

Metagenomic Analysis of Virus Diversity and Relative Abundance in a Eutrophic Freshwater Harbour

Aquatic viruses have been extensively studied over the past decade, yet fundamental aspects of freshwater virus communities remain poorly described. Our goal was to characterize virus communities captured in the >0.22 µm size-fraction seasonally and spatially in a freshwater harbour. Community DNA was extracted from water samples and sequenced on an Illumina HiSeq platform. Assembled contigs were annotated as belonging to the virus groups (i.e., order or family) Caudovirales, Mimiviridae, Phycodnaviridae, and virophages (Lavidaviridae), or to other groups of undefined viruses. Virophages were often the most abundant group, and discrete virophage taxa were remarkably stable across sites and dates despite fluctuations in Mimiviridae community composition. Diverse Mimiviridae contigs were detected in the samples and the two sites contained distinct Mimiviridae communities, suggesting that Mimiviridae are important algal viruses in this system. Caudovirales and Phycodnaviridae were present at low abundances in most samples. Of the 18 environmental parameters tested, only chlorophyll a explained the variation in the data at the order or family level of classification. Overall, our findings provide insight into freshwater virus community assemblages by expanding the documented diversity of freshwater virus communities, highlighting the potential ecological importance of virophages, and revealing distinct communities over small spatial scales.

The Earth Is Small for "Leviathans": Long Distance Dispersal of Giant Viruses across Aquatic Environments

Giant viruses of ‘Megaviridae’ have the ability to widely disperse around the globe. We herein examined ‘Megaviridae’ communities in four distinct aquatic environments (coastal and offshore seawater, brackish water, and hot spring freshwater), which are distantly located from each other (between 74 and 1,765 km), using a meta-barcoding method. We identified between 593 and 3,627 OTUs in each sample. Some OTUs were detected in all five samples tested as well as in many of the Tara Oceans metagenomes, suggesting the existence of viruses of this family in a wide range of habitats and the ability to circulate on the planet.

Two classes of EF1-family translational GTPases encoded by giant viruses

Giant viruses have extraordinarily large dsDNA genomes, and exceptionally, they encode various components of the translation apparatus, including tRNAs, aminoacyl-tRNA synthetases and translation factors. Here, we focused on the elongation factor 1 (EF1) family of viral translational GTPases (trGTPases), using computational and functional approaches to shed light on their functions. Multiple sequence alignment indicated that these trGTPases clustered into two groups epitomized by members of Mimiviridae and Marseilleviridae, respectively. trGTPases in the first group were more closely related to GTP-binding protein 1 (GTPBP1), whereas trGTPases in the second group were closer to eEF1A, eRF3 and Hbs1. Functional characterization of representative GTPBP1-like trGTPases (encoded by Hirudovirus, Catovirus and Moumouvirus) using in vitro reconstitution revealed that they possess eEF1A-like activity and can deliver cognate aa-tRNAs to the ribosomal A site during translation elongation. By contrast, representative eEF1A/eRF3/Hbs1-like viral trGTPases, encoded by Marseillevirus and Lausannevirus, have eRF3-like termination activity and stimulate peptide release by eRF1. Our analysis identified specific aspects of the functioning of these viral trGTPases with eRF1 of human, amoebal and Marseillevirus origin.

Genome and Environmental Activity of a Chrysochromulina parva Virus and Its Virophages

Some giant viruses are ecological agents that are predicted to be involved in the top-down control of single-celled eukaryotic algae populations in aquatic ecosystems. Despite an increased interest in giant viruses since the discovery and characterization of Mimivirus and other viral giants, little is known about their physiology and ecology. In this study, we characterized the genome and functional potential of a giant virus that infects the freshwater haptophyte Chrysochromulina parva, originally isolated from Lake Ontario. This virus, CpV-BQ2, is a member of the nucleo-cytoplasmic large DNA virus (NCLDV) group and possesses a 437 kb genome encoding 503 ORFs with a GC content of 25%. Phylogenetic analyses of core NCLDV genes place CpV-BQ2 amongst the emerging group of algae-infecting Mimiviruses informally referred to as the “extended Mimiviridae,” making it the first virus of this group to be isolated from a freshwater ecosystem. During genome analyses, we also captured and described the genomes of three distinct virophages that co-occurred with CpV-BQ2 and likely exploit CpV for their own replication. These virophages belong to the polinton-like viruses (PLV) group and encompass 19–23 predicted genes, including all of the core PLV genes as well as several genes implicated in genome modifications. We used the CpV-BQ2 and virophage reference sequences to recruit reads from available environmental metatranscriptomic data to estimate their activity in fresh waters. We observed moderate recruitment of both virus and virophage transcripts in samples obtained during Microcystis aeruginosa blooms in Lake Erie and Lake Tai, China in 2013, with a spike in activity in one sample. Virophage transcript abundance for two of the three isolates strongly correlated with that of the CpV-BQ2. Together, the results highlight the importance of giant viruses in the environment and establish a foundation for future research on the physiology and ecology CpV-BQ2 as a model system for algal Mimivirus dynamics in freshwaters.

Diversity of Active Viral Infections within the Sphagnum Microbiome

Sphagnum-dominated peatlands play an important role in global carbon storage and represent significant sources of economic and ecological value. While recent efforts to describe microbial diversity and metabolic potential of the Sphagnum microbiome have demonstrated the importance of its microbial community, little is known about the viral constituents. We used metatranscriptomics to describe the diversity and activity of viruses infecting microbes within the Sphagnum peat bog. The vegetative portions of six Sphagnum plants were obtained from a peatland in northern Minnesota, and the total RNA was extracted and sequenced. Metatranscriptomes were assembled and contigs were screened for the presence of conserved virus marker genes. Using bacteriophage capsid protein gp23 as a marker for phage diversity, we identified 33 contigs representing undocumented phages that were active in the community at the time of sampling. Similarly, RNA-dependent RNA polymerase and the nucleocytoplasmic large DNA virus (NCLDV) major capsid protein were used as markers for single-stranded RNA (ssRNA) viruses and NCLDV, respectively. In total, 114 contigs were identified as originating from undescribed ssRNA viruses, 22 of which represent nearly complete genomes. An additional 64 contigs were identified as being from NCLDVs. Finally, 7 contigs were identified as putative virophage or polinton-like viruses. We developed co-occurrence networks with these markers in relation to the expression of potential-host housekeeping gene rpb1 to predict virus-host relationships, identifying 13 groups. Together, our approach offers new tools for the identification of virus diversity and interactions in understudied clades and suggests that viruses may play a considerable role in the ecology of the Sphagnum microbiome.IMPORTANCE Sphagnum-dominated peatlands play an important role in maintaining atmospheric carbon dioxide levels by modifying conditions in the surrounding soil to favor the growth of Sphagnum over that of other plant species. This lowers the rate of decomposition and facilitates the accumulation of fixed carbon in the form of partially decomposed biomass. The unique environment produced by Sphagnum enriches for the growth of a diverse microbial consortia that benefit from and support the moss's growth, while also maintaining the hostile soil conditions. While a growing body of research has begun to characterize the microbial groups that colonize Sphagnum, little is currently known about the ecological factors that constrain community structure and define ecosystem function. Top-down population control by viruses is almost completely undescribed. This study provides insight into the significant viral influence on the Sphagnum microbiome and identifies new potential model systems to study virus-host interactions in the peatland ecosystem.

Mimiviridae: An Expanding Family of Highly Diverse Large dsDNA Viruses Infecting a Wide Phylogenetic Range of Aquatic Eukaryotes

Since 1998, when Jim van Etten’s team initiated its characterization, Paramecium bursaria Chlorella virus 1 (PBCV-1) had been the largest known DNA virus, both in terms of particle size and genome complexity. In 2003, the Acanthamoeba-infecting Mimivirus unexpectedly superseded PBCV-1, opening the era of giant viruses, i.e., with virions large enough to be visible by light microscopy and genomes encoding more proteins than many bacteria. During the following 15 years, the isolation of many Mimivirus relatives has made Mimiviridae one of the largest and most diverse families of eukaryotic viruses, most of which have been isolated from aquatic environments. Metagenomic studies of various ecosystems (including soils) suggest that many more remain to be isolated. As Mimiviridae members are found to infect an increasing range of phytoplankton species, their taxonomic position compared to the traditional Phycodnaviridae (i.e., etymologically “algal viruses”) became a source of confusion in the literature. Following a quick historical review of the key discoveries that established the Mimiviridae family, we describe its current taxonomic structure and propose a set of operational criteria to help in the classification of future isolates.

Degenerate PCR Primers to Reveal the Diversity of Giant Viruses in Coastal Waters

“Megaviridae” is a proposed family of giant viruses infecting unicellular eukaryotes. These viruses are ubiquitous in the sea and have impact on marine microbial community structure and dynamics through their lytic infection cycle. However, their diversity and biogeography have been poorly characterized due to the scarce detection of Megaviridae sequences in metagenomes, as well as the limitation of reference sequences used to design specific primers for this viral group. Here, we propose a set of 82 degenerated primers (referred to as MEGAPRIMER), targeting DNA polymerase genes (polBs) of Megaviridae. MEGAPRIMER was designed based on 921 Megaviridae polBs from sequenced genomes and metagenomes. By applying this primer set to environmental DNA meta-barcoding of a coastal seawater sample, we report 5595 non-singleton operational taxonomic units (OTUs) of Megaviridae at 97% nucleotide sequence identity. The majority of the OTUs were found to form diverse clades, which were phylogenetically distantly related to known viruses such as Mimivirus. The Megaviridae OTUs detected in this study outnumber the giant virus OTUs identified in previous individual studies by more than an order of magnitude. Hence, MEGAPRIMER represents a useful tool to study the diversity of Megaviridae at the population level in natural environments.

Viruses of Eukaryotic Algae: Diversity, Methods for Detection, and Future Directions

The scope for ecological studies of eukaryotic algal viruses has greatly improved with the development of molecular and bioinformatic approaches that do not require algal cultures. Here, we review the history and perceived future opportunities for research on eukaryotic algal viruses. We begin with a summary of the 65 eukaryotic algal viruses that are presently in culture collections, with emphasis on shared evolutionary traits (e.g., conserved core genes) of each known viral type. We then describe how core genes have been used to enable molecular detection of viruses in the environment, ranging from PCR-based amplification to community scale "-omics" approaches. Special attention is given to recent studies that have employed network-analyses of -omics data to predict virus-host relationships, from which a general bioinformatics pipeline is described for this type of approach. Finally, we conclude with acknowledgement of how the field of aquatic virology is adapting to these advances, and highlight the need to properly characterize new virus-host systems that may be isolated using preliminary molecular surveys. Researchers can approach this work using lessons learned from the Chlorella virus system, which is not only the best characterized algal-virus system, but is also responsible for much of the foundation in the field of aquatic virology.

Transcriptional responses of the marine green alga Micromonas pusilla and an infecting prasinovirus under different phosphate conditions

Prasinophytes are widespread marine algae for which responses to nutrient limitation and viral infection are not well understood. We studied the picoprasinophyte, Micromonas pusilla, grown under phosphate-replete (0.65 ± 0.07 d^-1 ) and 10-fold lower (low)-phosphate (0.11 ± 0.04 d^-1 ) conditions, and infected by the phycodnavirus MpV-SP1. Expression of 17% of Micromonas genes in uninfected cells differed by >1.5-fold (q < 0.01) between nutrient conditions, with genes for P-metabolism and the uniquely-enriched Sel1-like repeat (SLR) family having higher relative transcript abundances, while phospholipid-synthesis genes were lower in low-P than P-replete. Approximately 70% (P-replete) and 30% (low-P) of cells were lysed 24 h post-infection, and expression of ≤5.8% of host genes changed relative to uninfected treatments. Host genes for CAZymes and glycolysis were activated by infection, supporting importance in viral production, which was significantly lower in slower growing (low-P) hosts. All MpV-SP1 genes were expressed, and our analyses suggest responses to differing host-phosphate backgrounds involve few viral genes, while the temporal program of infection involves many more, and is largely independent of host-phosphate background. Our study (i) identifies genes previously unassociated with nutrient acclimation or viral infection, (ii) provides insights into the temporal program of prasinovirus gene expression by hosts and (iii) establishes cell biological aspects of an ecologically important host-prasinovirus system that differ from other marine algal-virus systems.

Quantitative Infection Dynamics of Cafeteria Roenbergensis Virus

The discovery of giant viruses in unicellular eukaryotic hosts has raised new questions on the nature of viral life. Although many steps in the infection cycle of giant viruses have been identified, the quantitative life history traits associated with giant virus infection remain unknown or poorly constrained. In this study, we provide the first estimates of quantitative infection traits of a giant virus by tracking the infection dynamics of the bacterivorous protist Cafeteria roenbergensis and its lytic virus CroV. Leveraging mathematical models of infection, we quantitatively estimate the adsorption rate, onset of DNA replication, latency time, and burst size from time-series data. Additionally, by modulating the initial ratio of viruses to hosts, we also provide evidence of a potential MOI-dependence on adsorption and burst size. Our work provides a baseline characterization of giant virus infection dynamics relevant to ongoing efforts to understand the ecological role of giant viruses.

Taxon Richness of "Megaviridae" Exceeds those of Bacteria and Archaea in the Ocean

Since the discovery of the giant mimivirus, evolutionarily related viruses have been isolated or identified from various environments. Phylogenetic analyses of this group of viruses, tentatively referred to as the family “Megaviridae”, suggest that it has an ancient origin that may predate the emergence of major eukaryotic lineages. Environmental genomics has since revealed that Megaviridae represents one of the most abundant and diverse groups of viruses in the ocean. In the present study, we compared the taxon richness and phylogenetic diversity of Megaviridae, Bacteria, and Archaea using DNA-dependent RNA polymerase as a common marker gene. By leveraging existing microbial metagenomic data, we found higher richness and phylogenetic diversity in this single viral family than in the two prokaryotic domains. We also obtained results showing that the evolutionary rate alone cannot account for the observed high diversity of Megaviridae lineages. These results suggest that the Megaviridae family has a deep co-evolutionary history with diverse marine protists since the early “Big-Bang” radiation of the eukaryotic tree of life.

Microscopic and Molecular Evidence of the First Elasmobranch Adomavirus, the Cause of Skin Disease in a Giant Guitarfish, Rhynchobatus djiddensis

Only eight families of double-stranded DNA (dsDNA) viruses are known to infect vertebrate animals. During an investigation of papillomatous skin disease in an elasmobranch species, the giant guitarfish (Rhynchobatus djiddensis), a novel virus, distinct from all known viral families in regard to particle size, morphology, genome organization, and helicase phylogeny was discovered. Large inclusion bodies containing 75-nm icosahedral viral particles were present within epithelial cell nuclei in the proliferative skin lesions. Deep metagenomic sequencing revealed a 22-kb circular dsDNA viral genome, tentatively named guitarfish "adomavirus" (GAdoV), with only distant homology to two other fish viruses, Japanese eel endothelial cell-infecting virus (JEECV) and a recently reported marbled eel virus. Phylogenetic analysis of the helicase domain places the guitarfish virus in a novel clade that is equidistant between members of the Papillomaviridae and Polyomaviridae families. Specific PCR, quantitative PCR, and in situ hybridization were used to detect, quantify, and confirm that GAdoV DNA was localized to affected epithelial cell nuclei. Changes in the viral titer, as well as the presence of a hybridization signal, coincided with the progression and then final resolution of gross and microscopic lesions. The results indicate that GAdoV is the causative agent of the proliferative skin lesions.IMPORTANCE Cartilaginous fish, including the sharks and rays, evolved from ancestral fish species at least 400 million years ago. Even though they are the descendants of one of the most ancient vertebrate lineages, reports of viral diseases in these species are rare and poorly documented. Deep sequencing revealed a highly divergent virus, tentatively named guitarfish adomavirus, that is distantly related to known papillomaviruses and polyomaviruses. Out of the eight predicted viral genes, only the helicase could be identified as viral by sequence homology searches (BLAST), exemplifying the difficulties of discovering novel viruses within seas of unidentifiable "dark matter" associated with deep sequencing data. The novel adomavirus represents the first viral genome shown to cause clinical disease in a cartilaginous fish species, the giant guitarfish. Our findings demonstrate that emerging fish viruses are fertile ground to expand our understanding of viral evolution in vertebrates.

Infection by a Giant Virus (AaV) Induces Widespread Physiological Reprogramming in Aureococcus anophagefferens CCMP1984 - A Harmful Bloom Algae

While viruses with distinct phylogenetic origins and different nucleic acid types can infect and lyse eukaryotic phytoplankton, “giant” dsDNA viruses have been found to be associated with important ecological processes, including the collapse of algal blooms. However, the molecular aspects of giant virus–host interactions remain largely unknown. Aureococcus anophagefferens virus (AaV), a giant virus in the Mimiviridae clade, is known to play a critical role in regulating the fate of brown tide blooms caused by the pelagophyte Aureococcus anophagefferens. To understand the physiological response of A. anophagefferens CCMP1984 upon AaV infection, we studied the transcriptomic landscape of this host–virus pair over an entire infection cycle using a RNA-sequencing approach. A massive transcriptional response of the host was evident as early as 5 min post-infection, with modulation of specific processes likely related to both host defense mechanism(s) and viral takeover of the cell. Infected Aureococcus showed a relative suppression of host-cell transcripts associated with photosynthesis, cytoskeleton formation, fatty acid, and carbohydrate biosynthesis. In contrast, host cell processes related to protein synthesis, polyamine biosynthesis, cellular respiration, transcription, and RNA processing were overrepresented compared to the healthy cultures at different stages of the infection cycle. A large number of redox active host-selenoproteins were overexpressed, which suggested that viral replication and assembly progresses in a highly oxidative environment. The majority (99.2%) of annotated AaV genes were expressed at some point during the infection cycle and demonstrated a clear temporal–expression pattern and an increasing relative expression for the majority of the genes through the time course. We detected a putative early promoter motif for AaV, which was highly similar to the early promoter elements of two other Mimiviridae members, indicating some degree of evolutionary conservation of gene regulation within this clade. This large-scale transcriptome study provides insights into the Aureococcus cells infected by a giant virus and establishes a foundation to test hypotheses regarding metabolic and regulatory processes critical for AaV and other Mimiviridae members.

A giant virus infecting green algae encodes key fermentation genes

The family Mimiviridae contains uncommonly large viruses, many of which were isolated using a free-living amoeba as a host. Although the genomes of these and other mimivirids that infect marine heterokont and haptophyte protists have now been sequenced, there has yet to be a genomic investigation of a mimivirid that infects a member of the Viridiplantae lineage (green algae and land plants). Here we characterize the 668-kilobase complete genome of TetV-1, a mimivirid that infects the cosmopolitan green alga Tetraselmis (Chlorodendrophyceae). The analysis revealed genes not previously seen in viruses, such as the mannitol metabolism enzyme mannitol 1-phosphate dehydrogenase, the saccharide degradation enzyme alpha-galactosidase, and the key fermentation genes pyruvate formate-lyase and pyruvate formate-lyase activating enzyme. The TetV genome is the largest sequenced to date for a virus that infects a photosynthetic organism, and its genes reveal unprecedented mechanisms by which viruses manipulate their host's metabolism.

The kinetoplastid-infecting Bodo saltans virus (BsV), a window into the most abundant giant viruses in the sea

Giant viruses are ecologically important players in aquatic ecosystems that have challenged concepts of what constitutes a virus. Herein, we present the giant Bodo saltans virus (BsV), the first characterized representative of the most abundant group of giant viruses in ocean metagenomes, and the first isolate of a klosneuvirus, a subgroup of the Mimiviridae proposed from metagenomic data. BsV infects an ecologically important microzooplankton, the kinetoplastid Bodo saltans. Its 1.39 Mb genome encodes 1227 predicted ORFs, including a complex replication machinery. Yet, much of its translational apparatus has been lost, including all tRNAs. Essential genes are invaded by homing endonuclease-encoding self-splicing introns that may defend against competing viruses. Putative anti-host factors show extensive gene duplication via a genomic accordion indicating an ongoing evolutionary arms race and highlighting the rapid evolution and genomic plasticity that has led to genome gigantism and the enigma that is giant viruses.

Host Traits Drive Viral Life Histories across Phytoplankton Viruses

Viruses are integral to ecological and evolutionary processes, but we have a poor understanding of what drives variation in key traits across diverse viruses. For lytic viruses, burst size, latent period, and genome size are primary characteristics controlling host-virus dynamics. Here we synthesize data on these traits for 75 strains of phytoplankton viruses, which play an important role in global biogeochemistry. We find that primary traits of the host (genome size, growth rate) explain 40%-50% of variation in burst size and latent period. Specifically, burst size and latent period both exhibit saturating relationships versus the host∶virus genome size ratio, with both traits increasing at low genome size ratios while showing no relationship at high size ratios. In addition, latent period declines as host growth rate increases. We analyze a model of latent period evolution to explore mechanisms that could cause these patterns. The model predicts that burst size may often be set by the host genomic resources available for viral construction, while latent period evolves to permit this maximal burst size, modulated by host metabolic rate. These results suggest that general mechanisms may underlie the evolution of diverse viruses. Future extensions of this work could help explain viral regulation of host populations, viral influence on community structure and diversity, and viral roles in biogeochemical cycles.

Virus-host relationships of marine single-celled eukaryotes resolved from metatranscriptomics

Establishing virus–host relationships has historically relied on culture-dependent approaches. Here we report on the use of marine metatranscriptomics to probe virus–host relationships. Statistical co-occurrence analyses of dsDNA, ssRNA and dsRNA viral markers of polyadenylation-selected RNA sequences from microbial communities dominated by Aureococcus anophagefferens (Quantuck Bay, NY), and diatoms (Narragansett Bay, RI) show active infections by diverse giant viruses (NCLDVs) associated with algal and nonalgal hosts. Ongoing infections of A. anophagefferens by a known Mimiviridae (AaV) occur during bloom peak and decline. Bloom decline is also accompanied by increased activity of viruses other than AaV, including (+) ssRNA viruses. In Narragansett Bay, increased temporal resolution reveals active NCLDVs with both ‘boom-and-bust’ and ‘steady-state infection’-like ecologies that include known as well as novel virus–host interactions. Our approach offers a method for screening active viral infections and develops links between viruses and their potential hosts in situ. Our observations further demonstrate that previously unknown virus–host relationships in marine systems are abundant.

Viruses are partners in ecosystem ecology, yet their study has been primarily limited to laboratory models virus-host or derived from metagenomics. Here, Moniruzzaman et al. use metatranscriptomics to resolve interactions between giant viruses and single-celled eukaryotic hosts.

Comparative Genomics of Chrysochromulina Ericina Virus and Other Microalga-Infecting Large DNA Viruses Highlights Their Intricate Evolutionary Relationship with the Established Mimiviridae Family

Chrysochromulina ericina virus CeV-01B (CeV) was isolated from Norwegian coastal waters in 1998. Its icosahedral particle is 160 nm in diameter and encloses a 474-kb double-stranded DNA (dsDNA) genome. This virus, although infecting a microalga (the haptophyceae Haptolina ericina, formerly Chrysochromulina ericina), is phylogenetically related to members of the Mimiviridae family, initially established with the acanthamoeba-infecting mimivirus and megavirus as prototypes. This family was later split into two genera (Mimivirus and Cafeteriavirus) following the characterization of a virus infecting the heterotrophic stramenopile Cafeteria roenbergensis (CroV). CeV, as well as two of its close relatives, which infect the unicellular photosynthetic eukaryotes Phaeocystis globosa (Phaeocystis globosa virus [PgV]) and Aureococcus anophagefferens (Aureococcus anophagefferens virus [AaV]), are currently unclassified by the International Committee on Viral Taxonomy (ICTV). The detailed comparative analysis of the CeV genome presented here confirms the phylogenetic affinity of this emerging group of microalga-infecting viruses with the Mimiviridae but argues in favor of their classification inside a distinct clade within the family. Although CeV, PgV, and AaV share more common features among them than with the larger Mimiviridae, they also exhibit a large complement of unique genes, attesting to their complex evolutionary history. We identified several gene fusion events and cases of convergent evolution involving independent lateral gene acquisitions. Finally, CeV possesses an unusual number of inteins, some of which are closely related despite being inserted in nonhomologous genes. This appears to contradict the paradigm of allele-specific inteins and suggests that the Mimiviridae are especially efficient in spreading inteins while enlarging their repertoire of homing genes.IMPORTANCE Although it infects the microalga Chrysochromulina ericina, CeV is more closely related to acanthamoeba-infecting viruses of the Mimiviridae family than to any member of the Phycodnaviridae, the ICTV-approved family historically including all alga-infecting large dsDNA viruses. CeV, as well as its relatives that infect the microalgae Phaeocystic globosa (PgV) and Aureococcus anophagefferens (AaV), remains officially unclassified and a source of confusion in the literature. Our comparative analysis of the CeV genome in the context of this emerging group of alga-infecting viruses suggests that they belong to a distinct clade within the established Mimiviridae family. The presence of a large number of unique genes as well as specific gene fusion events, evolutionary convergences, and inteins integrated at unusual locations document the complex evolutionary history of the CeV lineage.

Genomic exploration of individual giant ocean viruses

Viruses are major pathogens in all biological systems. Virus propagation and downstream analysis remains a challenge, particularly in the ocean where the majority of their microbial hosts remain recalcitrant to current culturing techniques. We used a cultivation-independent approach to isolate and sequence individual viruses. The protocol uses high-speed fluorescence-activated virus sorting flow cytometry, multiple displacement amplification (MDA), and downstream genomic sequencing. We focused on 'giant viruses' that are readily distinguishable by flow cytometry. From a single-milliliter sample of seawater collected from off the dock at Boothbay Harbor, ME, USA, we sorted almost 700 single virus particles, and subsequently focused on a detailed genome analysis of 12. A wide diversity of viruses was identified that included Iridoviridae, extended Mimiviridae and even a taxonomically novel (unresolved) giant virus. We discovered a viral metacaspase homolog in one of our sorted virus particles and discussed its implications in rewiring host metabolism to enhance infection. In addition, we demonstrated that viral metacaspases are widespread in the ocean. We also discovered a virus that contains both a reverse transcriptase and a transposase; although highly speculative, we suggest such a genetic complement would potentially allow this virus to exploit a latency propagation mechanism. Application of single virus genomics provides a powerful opportunity to circumvent cultivation of viruses, moving directly to genomic investigation of naturally occurring viruses, with the assurance that the sequence data is virus-specific, non-chimeric and contains no cellular contamination.

Seasonal Dynamics of Haptophytes and dsDNA Algal Viruses Suggest Complex Virus-Host Relationship

Viruses influence the ecology and diversity of phytoplankton in the ocean. Most studies of phytoplankton host-virus interactions have focused on bloom-forming species like Emiliania huxleyi or Phaeocystis spp. The role of viruses infecting phytoplankton that do not form conspicuous blooms have received less attention. Here we explore the dynamics of phytoplankton and algal viruses over several sequential seasons, with a focus on the ubiquitous and diverse phytoplankton division Haptophyta, and their double-stranded DNA viruses, potentially with the capacity to infect the haptophytes. Viral and phytoplankton abundance and diversity showed recurrent seasonal changes, mainly explained by hydrographic conditions. By 454 tag-sequencing we revealed 93 unique haptophyte operational taxonomic units (OTUs), with seasonal changes in abundance. Sixty-one unique viral OTUs, representing Megaviridae and Phycodnaviridae, showed only distant relationship with currently isolated algal viruses. Haptophyte and virus community composition and diversity varied substantially throughout the year, but in an uncoordinated manner. A minority of the viral OTUs were highly abundant at specific time-points, indicating a boom-bust relationship with their host. Most of the viral OTUs were very persistent, which may represent viruses that coexist with their hosts, or able to exploit several host species.

A Student's Guide to Giant Viruses Infecting Small Eukaryotes: From Acanthamoeba to Zooxanthellae

The discovery of infectious particles that challenge conventional thoughts concerning "what is a virus" has led to the evolution a new field of study in the past decade. Here, we review knowledge and information concerning "giant viruses", with a focus not only on some of the best studied systems, but also provide an effort to illuminate systems yet to be better resolved. We conclude by demonstrating that there is an abundance of new host-virus systems that fall into this "giant" category, demonstrating that this field of inquiry presents great opportunities for future research.