Host Range and Coding Potential of Eukaryotic Giant Viruses

Expansion of the genomic and functional diversity of global ocean giant viruses

Giant viruses (GVs) play crucial roles in the global ocean microbial food web and biogeochemistry. Recent metagenomic advances have uncovered >1800 new GV genomes from the world's oceans. While this rapid increase in genomic information is impressive, it is nowhere close to the extensive genomic information available for other marine entities-e.g., prokaryotes and their "virome". We present 230 new high-quality GV genomes (genomes with 4 or more marker genes) and 398 partial genomes from nine global ocean datasets. Notably, we identified numerous GV genomes from the Baltic Sea, offering insights into their phylogenomics, metabolic potential, and environmental drivers in one of the largest brackish water ecosystems. We discovered new GV functions and identified a significant functional divide between the Imitervirales and Algavirales orders. Additionally, we evaluated factors affecting GV abundance through a case study on the Baltic Sea dataset. Our study significantly expands the marine GV genomic and functional diversity, broadening our understanding of their roles in the food web and biogeochemistry.

Haptophyte-infecting viruses change the genome condensing proteins of dinoflagellates

Giant viruses are extraordinary members of the virosphere due to their structural complexity and high diversity in gene content. Haptophytes are ecologically important primary producers in the ocean, and all known viruses that infect haptophytes are giant viruses. However, little is known about the specifics of their infection cycles and the responses they trigger in their host cells. Our in-depth electron microscopic, phylogenomic and virion proteomic analyses of two haptophyte-infecting giant viruses, Haptolina ericina virus RF02 (HeV RF02) and Prymnesium kappa virus RF02 (PkV RF02), unravel their large capacity for host manipulation and arsenals that function during the infection cycle from virus entry to release. The virus infection induces significant morphological changes in the host cell that is manipulated to build a virus proliferation factory. Both viruses' genomes encode a putative nucleoprotein (dinoflagellate/viral nucleoprotein; DVNP), which was also found in the virion proteome of PkV RF02. Phylogenetic analysis suggests that DVNPs are widespread in marine giant metaviromes. Furthermore, the analysis shows that the dinoflagellate homologues were possibly acquired from viruses of the order Imitervirales. These findings enhance our understanding of how viruses impact the biology of microalgae, providing insights into evolutionary biology, ecosystem dynamics, and nutrient cycling in the ocean.

Yaravirus brasiliense genomic structure analysis and its possible influence on the metabolism

Here we analyze the Yaravirus brasiliense, an amoeba-infecting 80-nm-sized virus with a 45-kbp dsDNA, using structural molecular modeling. Almost all of its 74 genes were previously identified as ORFans. Considering its unprecedented genetic content, we analyzed Yaravirus genome to understand its genetic organization, its proteome, and how it interacts with its host. We reported possible functions for all Yaravirus proteins. Our results suggest the first ever report of a fragment proteome, in which the proteins are separated in modules and joined together at a protein level. Given the structural resemblance between some Yaravirus proteins and proteins related to tricarboxylic acid cycle (TCA), glyoxylate cycle, and the respiratory complexes, our work also allows us to hypothesize that these viral proteins could be modulating cell metabolism by upregulation. The presence of these TCA cycle-related enzymes specifically could be trying to overcome the cycle's control points, since they are strategic proteins that maintain malate and oxaloacetate levels. Therefore, we propose that Yaravirus proteins are redirecting energy and resources towards viral production, and avoiding TCA cycle control points, "unlocking" the cycle. Altogether, our data helped understand a previously almost completely unknown virus, and a little bit more of the incredible diversity of viruses.

Host-Calibrated Time Tree Caps the Age of Giant Viruses

Viruses are widespread parasites with important impacts on public health, economy, and ecosystems. However, little is known about their origins, ages, and early evolutionary relationships with hosts. Here, we infer the maximum divergence times for eukaryotic giant DNA viruses (phylum Nucleocytoviricota) with dating analyses calibrated by host taxon ages of virus lineages with specific host ranges. The last common ancestor of Nucleocytoviricota existed after 1,000 million years ago, suggesting a much later origin than that of the eukaryotes. The early evolution of Nucleocytoviricota either coincided with or postdated a substantial increase in the oxygen levels on the Earth's surface during the Neoproterozoic Era. The lineage diversification of giant viruses was frequently associated with host shifts, including two major transitions from amoebozoan hosts to animal hosts that eventually led to the emergence of iridoviruses and African swine fever viruses within the last 450 million years. These results outline the evolutionary timescale of a major virus group and are pivotal for further understanding the virus-host interactions and their potential ecological roles in the Earth's history.

Antiviral Mx proteins have an ancient origin and widespread distribution among eukaryotes

Mx proteins, first identified in mammals, encode potent antiviral activity against a wide range of viruses. Mx proteins arose within the Dynamin superfamily of proteins (DSP), which mediate critical cellular processes, such as endocytosis and mitochondrial, plastid, and peroxisomal dynamics. Despite their crucial role, the evolutionary origins of Mx proteins are poorly understood. Through comprehensive phylogenomic analyses with progressively expanded taxonomic sampling, we demonstrate that Mx proteins predate the interferon signaling system in vertebrates. Our analyses find an ancient monophyletic DSP lineage in eukaryotes that groups vertebrate and invertebrate Mx proteins with fungal MxF proteins, the largely uncharacterized plant and algal Dynamin 4A/4C proteins, and representatives from several other eukaryotic lineages, suggesting that Mx-like proteins date back close to the origin of Eukarya. Our phylogenetic analyses also find host-encoded and nucleocytoplasmic large DNA viruses-encoded DSPs interspersed in four distinct DSP lineages, indicating recurrent viral theft of host DSPs. Our analyses thus reveal an ancient history of viral and antiviral functions encoded by the Dynamin superfamily in eukaryotes.

Genome-resolved year-round dynamics reveal a broad range of giant virus microdiversity

Giant viruses are crucial for marine ecosystem dynamics because they regulate microeukaryotic community structure, accelerate carbon and nutrient cycles, and drive the evolution of their hosts through co-evolutionary processes. Previously reported long-term observations revealed that these viruses display seasonal fluctuations in abundance. However, the underlying genetic mechanisms driving such dynamics of these viruses remain largely unknown. In this study, we investigated the dynamics of giant viruses using time-series metagenomes from eutrophic coastal seawater samples collected over 20 months. A newly developed computational pipeline generated 1,065 high-quality genomes covering six major giant virus lineages. These genomic data revealed year-round recovery of the viral community structure at the study site and distinct dynamics of viral populations that were classified as persistent (n = 9), seasonal (n = 389), sporadic (n = 318), or others. By profiling the intra-species nucleotide-resolved microdiversity through read mapping, we also identified year-round recovery dynamics at subpopulation level for viruses classified as persistent or seasonal. Our results further indicated that giant viruses with broader niche breadth tended to exhibit higher levels of microdiversity. We argue that greater microdiversity of viruses likely enhances adaptability and thus survival under the virus-host arms race during prolonged interactions with their hosts.IMPORTANCERecent genome-resolved metagenomic surveys have uncovered the vast genomic diversity of giant viruses, which play significant roles in aquatic ecosystems by acting as bloom terminators and influencing biogeochemical cycles. However, the relationship between the ecological dynamics of giant viruses and underlying genetic structures of viral populations remains unresolved. In this study, we performed deep metagenomic sequencing of seawater samples collected across a time-series from a coastal area in Japan. The results revealed a significant positive correlation between microdiversity and temporal persistence of giant virus populations, suggesting that population structure is a crucial factor for adaptation and survival in the interactions with their hosts.

Protists as mediators of complex microbial and viral associations

Microbial eukaryotes (aka protists) are known for their important roles in nutrient cycling across different ecosystems. However, the composition and function of protist-associated microbiomes remains largely elusive. Here, we employ cultivation-independent single-cell isolation and genome-resolved metagenomics to provide detailed insights into underexplored microbiomes and viromes of over 100 currently uncultivable ciliates and amoebae isolated from diverse environments. Our findings reveal unique microbiome compositions and hint at an intricate network of complex interactions and associations with bacterial symbionts and viruses. We observed stark differences between ciliates and amoebae in terms of microbiome and virome compositions, highlighting the specificity of protist-microbe interactions. Over 115 of the recovered microbial genomes were affiliated with known endosymbionts of eukaryotes, including diverse members of the Holosporales, Rickettsiales, Legionellales, Chlamydiae, Dependentiae , and more than 250 were affiliated with possible host-associated bacteria of the phylum Patescibacteria. We also identified more than 80 giant viruses belonging to diverse viral lineages, of which some were actively expressing genes in single cell transcriptomes, suggesting a possible association with the sampled protists. We also revealed a wide range of other viruses that were predicted to infect eukaryotes or host-associated bacteria. Our results provide further evidence that protists serve as mediators of complex microbial and viral associations, playing a critical role in ecological networks. The frequent co-occurrence of giant viruses and diverse microbial symbionts in our samples suggests multipartite associations, particularly among amoebae. Our study provides a preliminary assessment of the microbial diversity associated with lesser-known protist lineages and paves the way for a deeper understanding of protist ecology and their roles in environmental and human health.

Viral dynamics in a high-rate algal pond reveals a burst of Phycodnaviridae diversity correlated with episodic algal mortality

This study explores virus-host dynamics in a unique environment: an industrial high-rate algal pond (HRAP). A wealth of novel DNA algal viruses are revealed, including members of Nucleocytoviricota "giant viruses" and the enigmatic Preplasmiviricota (e.g., virophages and polinton-like viruses). Several species of single-celled eukaryotic photosynthetic algae are identified (Chlorophyta) as putative hosts, with alternating dominant populations during the year of study. We specifically observe a surprising diversity of giant viruses from the family Phycodnaviridae (Nucleocytoviricota), including phylogenetically related but highly diversified genotypes appearing in the HRAP that we suggest are implicated in bloom collapse. We hypothesize that these related Phycodnaviridae lineages infect the same algal species of the genus Picochlorum that has been identified in the HRAP. This study establishes a baseline for comprehending the role viruses play in algal farming and emphasizes the necessity of controlling the viral load in future culture system development to optimize algal growth.

IMPORTANCE

The virosphere is ubiquitous, but we have yet to characterize many environments where viruses exist. In an industrial polyculture of microalgae, a wealth of viruses persist, their diversity and dynamics changing over time and consequently give evidence of their evolution and ecological strategies. Several notable infectious agents of the culture's algae appear, including giant viruses, polinton-like viruses, and a virophage. As our reliance and interest in algal compound-based cosmetics, pharmaceuticals, and bio-plastics increases, so must our understanding of these systems, including the unique viruses that appear there.

Genomics and evolutionary analysis of Chlorella variabilis-infecting viruses demarcate criteria for defining species of giant viruses

Chloroviruses exhibit a close relationship with their hosts with the phenotypic aspect of their ability to form lytic plaques having primarily guided the taxonomy. However, with the isolation of viruses that are only able to complete their replication cycle in one strain of Chlorella variabilis, systematic challenges emerged. In this study, we described the genomic features of 53 new chlorovirus isolates and used them to elucidate part of the evolutionary history and taxonomy of this clade. Our analysis revealed new chloroviruses with the largest genomes to date (>400 kbp) and indicated that four genomic features are statistically different in the viruses that only infect the Syngen 2-3 strain of C. variabilis (OSy viruses). We found large regions of dissimilarity in the genomes of viruses PBCV-1 and OSy-NE5 when compared with the other genomes. These regions contained genes related to the interaction with the host cell machinery and viral capsid proteins, which provided insights into the evolution of the replicative and structural modules in these giant viruses. Phylogenetic analysis using hallmark genes of Nucleocytoviricota revealed that OSy-viruses evolved from the NC64A-viruses, possibly emerging as a result of the strict relationship with their hosts. Merging phylogenetics and nucleotide identity analyses, we propose strategies to demarcate viral species, resulting in seven new species of chloroviruses. Collectively, our results show how genomic data can be used as lines of evidence to demarcate viral species. Using the chloroviruses as a case study, we expect that similar initiatives will emerge using the basis exhibited here.IMPORTANCEChloroviruses are a group of giant viruses with long dsDNA genomes that infect different species of Chlorella-like green algae. They are host-specific, and some isolates can only replicate within a single strain of Chlorella variabilis. The genomics of these viruses is still poorly explored, and the characterization of new isolates provides important data on their genetic diversity and evolution. In this work, we describe 53 new chlorovirus genomes, including many isolated from alkaline lakes for the first time. Through comparative genomics and molecular phylogeny, we provide evidence of genomic gigantism in chloroviruses and show that a subset of viruses became highly specific for their hosts at a particular point in evolutionary history. We propose criteria to demarcate species of chloroviruses, paving the way for an update in the taxonomy of other groups of viruses. This study is a new and important piece in the complex puzzle of giant algal viruses.

Proteome-scale structural prediction of the giant Marseillevirus reveals conserved folds and putative homologs of the hypothetical proteins

A significant proportion of the highly divergent and novel proteins of giant viruses are termed "hypothetical" due to the absence of detectable homologous sequences in the existing databases. The quality of genome and proteome annotations often relies on the identification of signature sequences and motifs in order to assign putative functions to the gene products. These annotations serve as the first set of information for researchers to develop workable hypotheses for further experimental research. The structure-function relationship of proteins suggests that proteins with similar functions may also exhibit similar folding patterns. Here, we report the first proteome-wide structure prediction of the giant Marseillevirus. We use AlphaFold-predicted structures and their comparative analysis with the experimental structures in the PDB database to preliminarily annotate the viral proteins. Our work highlights the conservation of structural folds in proteins with highly divergent sequences and reveals potentially paralogous relationships among them. We also provide evidence for gene duplication and fusion as contributing factors to giant viral genome expansion and evolution. With the easily accessible AlphaFold and other advanced bioinformatics tools for high-confidence de novo structure prediction, we propose a combined sequence and predicted-structure-based proteome annotation approach for the initial characterization of novel and complex organisms or viruses.

Antiviral Mx proteins have an ancient origin and widespread distribution among eukaryotes

First identified in mammals, Mx proteins are potent antivirals against a broad swathe of viruses. Mx proteins arose within the Dynamin superfamily of proteins (DSP), mediating critical cellular processes, such as endocytosis and mitochondrial, plastid, and peroxisomal dynamics. And yet, the evolutionary origins of Mx proteins are poorly understood. Using a series of phylogenomic analyses with stepwise increments in taxonomic coverage, we show that Mx proteins predate the interferon signaling system in vertebrates. Our analyses find an ancient monophyletic DSP lineage in eukaryotes that groups vertebrate and invertebrate Mx proteins with previously undescribed fungal MxF proteins, the relatively uncharacterized plant and algal Dynamin 4A/4C proteins, and representatives from several early-branching eukaryotic lineages. Thus, Mx-like proteins date back close to the origin of Eukarya. Our phylogenetic analyses also reveal that host-encoded and NCLDV (nucleocytoplasmic large DNA viruses)-encoded DSPs are interspersed in four distinct DSP lineages, indicating recurrent viral theft of host DSPs. Our analyses thus reveal an ancient history of viral and antiviral functions encoded by the Dynamin superfamily in eukaryotes.

Gene Transfer Among Viruses Substantially Contributes to Gene Gain of Giant Viruses

The phylum Nucleocytoviricota comprises a diverse group of double-stranded DNA viruses that display a wide range of gene repertoires. Although these gene repertoires determine the characteristics of individual viruses, the evolutionary processes that have shaped the gene repertoires of extant viruses since their common ancestor are poorly characterized. In this study, we aimed to address this gap in knowledge by using amalgamated likelihood estimation, a probabilistic tree reconciliation method that infers evolutionary scenarios by distinguishing origination, gene duplications, virus-to-virus horizontal gene transfer (vHGT), and gene losses. We analyzed over 4,700 gene families from 195 genomes spanning all known viral orders. The evolutionary reconstruction suggests a history of extensive gene gains and losses during the evolution of these viruses, notably with vHGT contributing to gene gains at a comparable level to duplications and originations. The vHGT frequently occurred between phylogenetically closely related viruses, as well as between distantly related viruses with an overlapping host range. We observed a pattern of massive gene duplications that followed vHGTs for gene families that was potentially related to host range control and virus-host arms race. These results suggest that vHGT represents a previously overlooked, yet important, evolutionary force that integrates the evolutionary paths of multiple viruses and affects shaping of Nucleocytoviricota virus gene repertoires.

Iron‑sulfur clusters in viral proteins: Exploring their elusive nature, roles and new avenues for targeting infections

Viruses have evolved complex mechanisms to exploit host factors for replication and assembly. In response, host cells have developed strategies to block viruses, engaging in a continuous co-evolutionary battle. This dynamic interaction often revolves around the competition for essential resources necessary for both host cell and virus replication. Notably, iron, required for the biosynthesis of several cofactors, including iron‑sulfur (FeS) clusters, represents a critical element in the ongoing competition for resources between infectious agents and host. Although several recent studies have identified FeS cofactors at the core of virus replication machineries, our understanding of their specific roles and the cellular processes responsible for their incorporation into viral proteins remains limited. This review aims to consolidate our current knowledge of viral components that have been characterized as FeS proteins and elucidate how viruses harness these versatile cofactors to their benefit. Its objective is also to propose that viruses may depend on incorporation of FeS cofactors more extensively than is currently known. This has the potential to revolutionize our understanding of viral replication, thereby carrying significant implications for the development of strategies to target infections.

A giant virus infecting the amoeboflagellate Naegleria

Giant viruses (Nucleocytoviricota) are significant lethality agents of various eukaryotic hosts. Although metagenomics indicates their ubiquitous distribution, available giant virus isolates are restricted to a very small number of protist and algal hosts. Here we report on the first viral isolate that replicates in the amoeboflagellate Naegleria. This genus comprises the notorious human pathogen Naegleria fowleri, the causative agent of the rare but fatal primary amoebic meningoencephalitis. We have elucidated the structure and infection cycle of this giant virus, Catovirus naegleriensis (a.k.a. Naegleriavirus, NiV), and show its unique adaptations to its Naegleria host using fluorescence in situ hybridization, electron microscopy, genomics, and proteomics. Naegleriavirus is only the fourth isolate of the highly diverse subfamily Klosneuvirinae, and like its relatives the NiV genome contains a large number of translation genes, but lacks transfer RNAs (tRNAs). NiV has acquired genes from its Naegleria host, which code for heat shock proteins and apoptosis inhibiting factors, presumably for host interactions. Notably, NiV infection was lethal to all Naegleria species tested, including the human pathogen N. fowleri. This study expands our experimental framework for investigating giant viruses and may help to better understand the basic biology of the human pathogen N. fowleri.

Adaptation strategies of giant viruses to low-temperature marine ecosystems

Microbes in marine ecosystems have evolved their gene content to thrive successfully in the cold. Although this process has been reasonably well studied in bacteria and selected eukaryotes, less is known about the impact of cold environments on the genomes of viruses that infect eukaryotes. Here, we analyzed cold adaptations in giant viruses (Nucleocytoviricota and Mirusviricota) from austral marine environments and compared them with their Arctic and temperate counterparts. We recovered giant virus metagenome-assembled genomes (98 Nucleocytoviricota and 12 Mirusviricota MAGs) from 61 newly sequenced metagenomes and metaviromes from sub-Antarctic Patagonian fjords and Antarctic seawater samples. When analyzing our data set alongside Antarctic and Arctic giant viruses MAGs already deposited in the Global Ocean Eukaryotic Viral database, we found that Antarctic and Arctic giant viruses predominantly inhabit sub-10°C environments, featuring a high proportion of unique phylotypes in each ecosystem. In contrast, giant viruses in Patagonian fjords were subject to broader temperature ranges and showed a lower degree of endemicity. However, despite differences in their distribution, giant viruses inhabiting low-temperature marine ecosystems evolved genomic cold-adaptation strategies that led to changes in genetic functions and amino acid frequencies that ultimately affect both gene content and protein structure. Such changes seem to be absent in their mesophilic counterparts. The uniqueness of these cold-adapted marine giant viruses may now be threatened by climate change, leading to a potential reduction in their biodiversity.

Isolation of a widespread giant virus implicated in cryptophyte bloom collapse

Photosynthetic cryptophytes are ubiquitous protists that are major participants in the freshwater phytoplankton bloom at the onset of spring. Mortality due to change in environmental conditions and grazing have been recognized as key factors contributing to bloom collapse. In contrast, the role of viral outbreaks as factors terminating phytoplankton blooms remains unknown from freshwaters. Here, we isolated and characterized a cryptophyte virus contributing to the annual collapse of a natural cryptophyte spring bloom population. This viral isolate is also representative for a clade of abundant giant viruses (phylum Nucleocytoviricota) found in freshwaters all over the world.

Hecatomb: an integrated software platform for viral metagenomics

BACKGROUND

Modern sequencing technologies offer extraordinary opportunities for virus discovery and virome analysis. Annotation of viral sequences from metagenomic data requires a complex series of steps to ensure accurate annotation of individual reads and assembled contigs. In addition, varying study designs will require project-specific statistical analyses.

FINDINGS

Here we introduce Hecatomb, a bioinformatic platform coordinating commonly used tasks required for virome analysis. Hecatomb means "a great sacrifice." In this setting, Hecatomb is "sacrificing" false-positive viral annotations using extensive quality control and tiered-database searches. Hecatomb processes metagenomic data obtained from both short- and long-read sequencing technologies, providing annotations to individual sequences and assembled contigs. Results are provided in commonly used data formats useful for downstream analysis. Here we demonstrate the functionality of Hecatomb through the reanalysis of a primate enteric and a novel coral reef virome.

CONCLUSION

Hecatomb provides an integrated platform to manage many commonly used steps for virome characterization, including rigorous quality control, host removal, and both read- and contig-based analysis. Each step is managed using the Snakemake workflow manager with dependency management using Conda. Hecatomb outputs several tables properly formatted for immediate use within popular data analysis and visualization tools, enabling effective data interpretation for a variety of study designs. Hecatomb is hosted on GitHub (github.com/shandley/hecatomb) and is available for installation from Bioconda and PyPI.

From Mimivirus to Mirusvirus: The Quest for Hidden Giants

Our perception of viruses has been drastically evolving since the inception of the field of virology over a century ago. In particular, the discovery of giant viruses from the Nucleocytoviricota phylum marked a pivotal moment. Their previously concealed diversity and abundance unearthed an unprecedented complexity in the virus world, a complexity that called for new definitions and concepts. These giant viruses underscore the intricate interactions that unfold over time between viruses and their hosts, and are themselves suspected to have played a significant role as a driving force in the evolution of eukaryotes since the dawn of this cellular domain. Whether they possess exceptional relationships with their hosts or whether they unveil the actual depths of evolutionary connections between viruses and cells otherwise hidden in smaller viruses, the attraction giant viruses exert on the scientific community and beyond continues to grow. Yet, they still hold surprises. Indeed, the recent identification of mirusviruses connects giant viruses to herpesviruses, each belonging to distinct viral realms. This discovery substantially broadens the evolutionary landscape of Nucleocytoviricota. Undoubtedly, the years to come will reveal their share of surprises.

Widespread Distribution and Evolution of Poxviral Entry-Fusion Complex Proteins in Giant Viruses

Poxviruses are known to encode a set of proteins that form an entry-fusion complex (EFC) to mediate virus entry. However, the diversity, evolution, and origin of these EFC proteins remain poorly understood. Here, we identify the EFC protein homologs in poxviruses and other giant viruses of the phylum Nucleocytoviricota. The 11 EFC genes are present in almost all poxviruses, with the two smallest, G3 and O3, being absent in Entomopoxvirinae and basal lineages of Chordopoxvirinae. Five of the EFC genes are further grouped into two families, A16/G9/J5 and F9/L1, which are widely distributed across other major lineages of Nucleocytoviricota, including metagenome-assembled genomes, but are generally absent in viruses infecting algae or nonamoebozoan heterotrophic protists. The A16/G9/J5 and F9/L1 families cooccur, mostly as single copies, in 93% of the non-Poxviridae giant viruses that have at least one of them. Distribution and phylogenetic patterns suggest that both families originated in the ancestor of Nucleocytoviricota. In addition to the Poxviridae genes, homologs from each of the other Nucleocytoviricota families are largely clustered together, suggesting their ancient presence and vertical inheritance. Despite deep sequence divergences, we observed noticeable conservation of cysteine residues and predicted structures between EFC proteins of Poxviridae and other families. Overall, our study reveals widespread distribution of these EFC protein homologs beyond poxviruses, implies the existence of a conserved membrane fusion mechanism, and sheds light on host range and ancient evolution of Nucleocytoviricota. IMPORTANCE Fusion between virus and host membranes is critical for viruses to release genetic materials and to initiate infection. Whereas most viruses use a single protein for membrane fusion, poxviruses employ a multiprotein entry-fusion complex (EFC). We report that two major families of the EFC proteins are widely distributed within the virus phylum Nucleocytoviricota, which includes poxviruses and other double-stranded (dsDNA) giant viruses that infect animals, amoebozoans, algae, and various microbial eukaryotes. Each of these two protein families is structurally conserved, traces its origin to the root of Nucleocytoviricota, was passed down to the major subclades of Nucleocytoviricota, and is retained in most giant viruses known to infect animals and amoebozoans. The EFC proteins therefore represent a potential mechanism for virus entry in diverse giant viruses. We hypothesize that they may have facilitated the infection of an animal/amoebozoan-like host by the last Nucleocytoviricota common ancestor.

Elimination of Foreign Sequences in Eukaryotic Viral Reference Genomes Improves the Accuracy of Virome Analysis

Widespread in public databases, foreign contaminant sequences pose a substantial obstacle in genomic analyses. Such contamination in viral genome databases is also notorious but more complicated and often causes questionable results in various applications, particularly in virome-based virus detection. Here, we conducted comprehensive screening and identification of the foreign sequences hidden in the largest eukaryotic viral genome collections of GenBank and UniProt using a scrutiny pipeline, which enables us to rigorously detect those problematic viral sequences (PVSs) with origins in hosts, vectors, and laboratory components. As a result, a total of 766 nucleotide PVSs and 276 amino acid PVSs with lengths up to 6,605 bp were determined, which were widely distributed in 39 families with many involving highly public health-concerning viruses, such as hepatitis C virus, Crimean-Congo hemorrhagic fever virus, and filovirus. The majority of these PVSs are genomic fragments of hosts including humans and bacteria. However, they cannot simply be regarded as foreign contaminants, since parts of them are results of natural occurrence or artificial engineering of viruses. Nevertheless, they severely disturb such sequence-based analyses as genome annotation, taxonomic assignment, and virome profiling. Therefore, we provide a clean version of the eukaryotic viral reference data set by the removal of these PVSs, which allows more accurate virome analysis with less time consumed than with other comprehensive databases. IMPORTANCE High-throughput sequencing-based viromics highly depends on reference databases, but foreign contamination is widespread in public databases and often leads to confusing and even wrong conclusions in genomic analysis and viromic profiling. To address this issue, we systematically detected and identified the contamination in the largest viral sequence collections of GenBank and UniProt based on a stringent scrutiny pipeline. We found hundreds of PVSs that are related to hosts, vectors, and laboratory components. By the removal of them, the resulting data set greatly improves the accuracy and efficiency of eukaryotic virome profiling. These results refresh our knowledge of the type and origin of PVSs and also have warning implications for viromic analysis. Viromic practitioners should be aware of these problems caused by PVSs and need to realize that a careful review of bioinformatic results is necessary for a reliable conclusion.

Giant virus biology and diversity in the era of genome-resolved metagenomics

The discovery of giant viruses, with capsids as large as some bacteria, megabase-range genomes and a variety of traits typically found only in cellular organisms, was one of the most remarkable breakthroughs in biology. Until recently, most of our knowledge of giant viruses came from ~100 species-level isolates for which genome sequences were available. However, these isolates were primarily derived from laboratory-based co-cultivation with few cultured protists and algae and, thus, did not reflect the true diversity of giant viruses. Although virus co-cultures enabled valuable insights into giant virus biology, many questions regarding their origin, evolution and ecological importance remain unanswered. With advances in sequencing technologies and bioinformatics, our understanding of giant viruses has drastically expanded. In this Review, we summarize our understanding of giant virus diversity and biology based on viral isolates as laboratory cultivation has enabled extensive insights into viral morphology and infection strategies. We then explore how cultivation-independent approaches have heightened our understanding of the coding potential and diversity of the Nucleocytoviricota. We discuss how metagenomics has revolutionized our perspective of giant viruses by revealing their distribution across our planet's biomes, where they impact the biology and ecology of a wide range of eukaryotic hosts and ultimately affect global nutrient cycles.

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.

Complementary Effects of Virus Population Are Required for Efficient Virus Infection

It is believed that the virions of a virus infecting a host may share the identical viral genome and characteristics. However, the role of genomic heterogeneity of the virions of a virus in virus infection has not been extensively explored. To address this issue, white spot syndrome virus (WSSV), a DNA virus infecting crustaceans, was characterized in the current study. In WSSV, differences in two nucleotides of the viral genome generated two types of WSSV, forming a virus population that consisted of Type A WSSV (encoding WSSV lncRNA-24) and Type B WSSV (encoding the wsv195 gene) at a ratio of 1:3. The virus populations in all virus-infected cells and tissues of different hosts exhibited a stable 1:3 structure. WSSV lncRNA-24 in Type A WSSV promoted virus infection by binding to shrimp and WSSV miRNAs, while the wsv195 gene in Type B WSSV played an essential role in virus infection. Loss of Type A WSSV or Type B WSSV in the WSSV population led to a 100-fold decrease in viral copy number in shrimp. Simultaneous loss of both types of WSSV prevented virus infection. These results indicated that the virus infection process was completed by two types of WSSV encoding different functional genes, revealing the complementary effects of WSSV population. Therefore, our study highlights the importance of the complementarity of virus population components in virus infection.

Systematic evaluation of horizontal gene transfer between eukaryotes and viruses

Gene exchange between viruses and their hosts acts as a key facilitator of horizontal gene transfer and is hypothesized to be a major driver of evolutionary change. Our understanding of this process comes primarily from bacteria and phage co-evolution, but the mode and functional importance of gene transfers between eukaryotes and their viruses remain anecdotal. Here we systematically characterized viral-eukaryotic gene exchange across eukaryotic and viral diversity, identifying thousands of transfers and revealing their frequency, taxonomic distribution and projected functions. Eukaryote-derived viral genes, abundant in the Nucleocytoviricota, highlighted common strategies for viral host-manipulation, including metabolic reprogramming, proteolytic degradation and extracellular modification. Furthermore, viral-derived eukaryotic genes implicate genetic exchange in the early evolution and diversification of eukaryotes, particularly through viral-derived glycosyltransferases, which have impacted structures as diverse as algal cell walls, trypanosome mitochondria and animal tissues. These findings illuminate the nature of viral-eukaryotic gene exchange and its impact on the evolution of viruses and their eukaryotic hosts.

Unraveling gene content variation across eukaryotic giant viruses based on network analyses and host associations

The nucleocytoplasmic large DNA viruses (NCLDVs, phylum Nucleocytoviricota) infect vertebrates, invertebrates, algae, amoebae, and other unicellular organisms across supergroups of eukaryotes and in various ecosystems. The expanding collection of their genome sequences has revolutionized our view of virus genome size and coding capacity. Phylogenetic trees based on a few core genes are commonly used as a model to understand their evolution. However, the tree topology can differ between analyses, and the vast majority of encoded genes might not share a common evolutionary history. To explore the whole-genome variation and evolution of NCLDVs, we dissected their gene contents using clustering, network, and comparative analyses. Our updated core-gene tree served as a framework to classify NCLDVs into families and intrafamilial lineages, but networks of individual genomes and family pangenomes showed patterns of gene sharing that contradict with the tree topology, in particular at higher taxonomic levels. Clustering of NCLDV genomes revealed variable granularity and degrees of gene sharing within each family, which cannot be inferred from the tree. At the level of NCLDV families, a correlation exists between gene content variation, but not core-gene sequence divergence, and host supergroup diversity. In addition, there is significantly higher gene sharing between divergent viruses that infect similar host types. The identified shared genes would be a useful resource for further functional analyses of NCLDV–host interactions. Overall this study provides a comprehensive view of gene repertoire variation in NCLDVs at different taxonomic levels, as well as a novel approach to studying the extremely diverse giant virus genomes.

The genomes of nucleocytoplasmic large DNA viruses: viral evolution writ large

The nucleocytoplasmic large DNA viruses (NCLDVs) are a diverse group that currently contain the largest known virions and genomes, also called giant viruses. The first giant virus was isolated and described nearly 20 years ago. Their genome sizes were larger than for any other known virus at the time and it contained a number of genes that had not been previously described in any virus. The origin and evolution of these unusually complex viruses has been puzzling, and various mechanisms have been put forward to explain how some NCLDVs could have reached genome sizes and coding capacity overlapping with those of cellular microbes. Here we critically discuss the evidence and arguments on this topic. We have also updated and systematically reanalysed protein families of the NCLDVs to further study their origin and evolution. Our analyses further highlight the small number of widely shared genes and extreme genomic plasticity among NCLDVs that are shaped via combinations of gene duplications, deletions, lateral gene transfers and de novo creation of protein-coding genes. The dramatic expansions of the genome size and protein-coding gene capacity characteristic of some NCLDVs is now increasingly understood to be driven by environmental factors rather than reflecting relationships to an ancient common ancestor among a hypothetical cellular lineage. Thus, the evolution of NCLDVs is writ large viral, and their origin, like all other viral lineages, remains unknown.

Giant Virus-Eukaryote Interactions as Ecological and Evolutionary Driving Forces

Giant DNA viruses of eukaryotes are notable for their extraordinary genome size and coding capacity. Once thought to be oddities in the virus world, these elusive microbes have turned out to be widely occurring in marine, freshwater, and terrestrial ecosystems and are commonly associated with diverse hosts, in particular microbial eukaryotes. This commentary discusses how new sequencing techniques and information can inform us about the interactions between giant viruses and eukaryotic hosts during the viral replication cycle and their implications for ecological and evolutionary processes across different spatiotemporal scales.

Discovery of Viral Myosin Genes With Complex Evolutionary History Within Plankton

Nucleocytoplasmic large DNA viruses (NCLDVs) infect diverse eukaryotes and form a group of viruses with capsids encapsulating large genomes. Recent studies are increasingly revealing a spectacular array of functions encoded in their genomes, including genes for energy metabolisms, nutrient uptake, as well as cytoskeleton. Here, we report the discovery of genes homologous to myosins, the major eukaryotic motor proteins previously unrecognized in the virosphere, in environmental genomes of NCLDVs from the surface of the oceans. Phylogenetic analyses indicate that most viral myosins (named "virmyosins") belong to the Imitervirales order, except for one belonging to the Phycodnaviridae family. On the one hand, the phylogenetic positions of virmyosin-encoding Imitervirales are scattered within the Imitervirales. On the other hand, Imitervirales virmyosin genes form a monophyletic group in the phylogeny of diverse myosin sequences. Furthermore, phylogenetic trends for the virmyosin genes and viruses containing them were incongruent. Based on these results, we argue that multiple transfers of myosin homologs have occurred not only from eukaryotes to viruses but also between viruses, supposedly during co-infections of the same host. Like other viruses that use host motor proteins for their intracellular transport or motility, these viruses may use the virally encoded myosins for the intracellular trafficking of giant viral particles.

Characterization of an aminotransferase from Acanthamoeba polyphaga Mimivirus

Acanthamoeba polyphaga Mimivirus, a complex virus that infects amoeba, was first reported in 2003. It is now known that its DNA genome encodes for nearly 1,000 proteins including enzymes that are required for the biosynthesis of the unusual sugar 4-amino-4,6-dideoxy-d-glucose, also known as d-viosamine. As observed in some bacteria, the pathway for the production of this sugar initiates with a nucleotide-linked sugar, which in the Mimivirus is thought to be UDP-d-glucose. The enzyme required for the installment of the amino group at the C-4' position of the pyranosyl moiety is encoded in the Mimivirus by the L136 gene. Here, we describe a structural and functional analysis of this pyridoxal 5'-phosphate-dependent enzyme, referred to as L136. For this analysis, three high-resolution X-ray structures were determined: the wildtype enzyme/pyridoxamine 5'-phosphate/dTDP complex and the site-directed mutant variant K185A in the presence of either UDP-4-amino-4,6-dideoxy-d-glucose or dTDP-4-amino-4,6-dideoxy-d-glucose. Additionally, the kinetic parameters of the enzyme utilizing either UDP-d-glucose or dTDP-d-glucose were measured and demonstrated that L136 is efficient with both substrates. This is in sharp contrast to the structurally related DesI from Streptomyces venezuelae, whose three-dimensional architecture was previously reported by this laboratory. As determined in this investigation, DesI shows a profound preference in its catalytic efficiency for the dTDP-linked sugar substrate. This difference can be explained in part by a hydrophobic patch in DesI that is missing in L136. Notably, the structure of L136 reported here represents the first three-dimensional model for a virally encoded PLP-dependent enzyme and thus provides new information on sugar aminotransferases in general.

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.

ViralRecall-A Flexible Command-Line Tool for the Detection of Giant Virus Signatures in 'Omic Data

Giant viruses are widespread in the biosphere and play important roles in biogeochemical cycling and host genome evolution. Also known as nucleo-cytoplasmic large DNA viruses (NCLDVs), these eukaryotic viruses harbor the largest and most complex viral genomes known. Studies have shown that NCLDVs are frequently abundant in metagenomic datasets, and that sequences derived from these viruses can also be found endogenized in diverse eukaryotic genomes. The accurate detection of sequences derived from NCLDVs is therefore of great importance, but this task is challenging owing to both the high level of sequence divergence between NCLDV families and the extraordinarily high diversity of genes encoded in their genomes, including some encoding for metabolic or translation-related functions that are typically found only in cellular lineages. Here, we present ViralRecall, a bioinformatic tool for the identification of NCLDV signatures in ‘omic data. This tool leverages a library of giant virus orthologous groups (GVOGs) to identify sequences that bear signatures of NCLDVs. We demonstrate that this tool can effectively identify NCLDV sequences with high sensitivity and specificity. Moreover, we show that it can be useful both for removing contaminating sequences in metagenome-assembled viral genomes as well as the identification of eukaryotic genomic loci that derived from NCLDV. ViralRecall is written in Python 3.5 and is freely available on GitHub: https://github.com/faylward/viralrecall.