Genomic exploration of individual giant ocean viruses

Single-cell RNA-seq of the rare virosphere reveals the native hosts of giant viruses in the marine environment

Giant viruses (phylum Nucleocytoviricota) are globally distributed in aquatic ecosystems. They play fundamental roles as evolutionary drivers of eukaryotic plankton and regulators of global biogeochemical cycles. However, we lack knowledge about their native hosts, hindering our understanding of their life cycle and ecological importance. In the present study, we applied a single-cell RNA sequencing (scRNA-seq) approach to samples collected during an induced algal bloom, which enabled pairing active giant viruses with their native protist hosts. We detected hundreds of single cells from multiple host lineages infected by diverse giant viruses. These host cells included members of the algal groups Chrysophycae and Prymnesiophycae, as well as heterotrophic flagellates in the class Katablepharidaceae. Katablepharids were infected with a rare Imitervirales-07 giant virus lineage expressing a large repertoire of cell-fate regulation genes. Analysis of the temporal dynamics of these host-virus interactions revealed an important role for the Imitervirales-07 in controlling the population size of the host Katablepharid population. Our results demonstrate that scRNA-seq can be used to identify previously undescribed host-virus interactions and study their ecological importance and impact.

New avenues for potentially seeking microbial responses to climate change beneath Antarctic ice shelves

The signs of climate change are undeniable, and the impact of these changes on ecosystem function heavily depends on the response of microbes that underpin the food web. Antarctic ice shelf is a massive mass of floating ice that extends from the continent into the ocean, exerting a profound influence on global carbon cycles. Beneath Antarctic ice shelves, marine ice stores valuable genetic information, where marine microbial communities before the industrial revolution are archived. Here, in this proof-of-concept, by employing a combination of single-cell technologiesand metagenomics, we have been able to sequence frozen microbial DNA (≈300 years old) stored in the marine ice core B15 collected from the Filchnner-Ronne Ice Shelf. Metagenomic data indicated that Proteobacteria and Thaumarchaeota (e.g., Nitrosopumilus spp.), followed by Actinobacteria (e.g., Actinomarinales), were abundant. Remarkably, our data allow us to "travel to the past" and calibrate genomic and genetic evolutionary changes for ecologically relevant microbes and functions, such as Nitrosopumilus spp., preserved in the marine ice (≈300 years old) with those collected recently in seawater under an ice shelf (year 2017). The evolutionary divergence for the ammonia monooxygenase gene amoA involved in chemolithoautotrophy was about 0.88 amino acid and 2.8 nucleotide substitution rate per 100 sites in a century, while the accumulated rate of genomic SNPs was 2,467 per 1 Mb of genome and 100 years. Whether these evolutionary changes remained constant over the last 300 years or accelerated during post-industrial periods remains an open question that will be further elucidated.

IMPORTANCE

Several efforts have been undertaken to predict the response of microbes under climate change, mainly based on short-term microcosm experiments under forced conditions. A common concern is that manipulative experiments cannot properly simulate the response of microbes to climate change, which is a long-term evolutionary process. In this proof-of-concept study with a limited sample size, we demonstrate a novel approach yet to be fully explored in science for accessing genetic information from putative past marine microbes preserved under Antarctic ice shelves before the industrial revolution. This potentially allows us estimating evolutionary changes as exemplified in our study. We advocate for gathering a more comprehensive Antarctic marine ice core data sets across various periods and sites. Such a data set would enable the establishment of a robust baseline, facilitating a better assessment of the potential effects of climate change on key genetic signatures of microbes.

Virologs, viral mimicry, and virocell metabolism: the expanding scale of cellular functions encoded in the complex genomes of giant viruses

The phylum Nucleocytoviricota includes the largest and most complex viruses known. These "giant viruses" have a long evolutionary history that dates back to the early diversification of eukaryotes, and over time they have evolved elaborate strategies for manipulating the physiology of their hosts during infection. One of the most captivating of these mechanisms involves the use of genes acquired from the host-referred to here as viral homologs or "virologs"-as a means of promoting viral propagation. The best-known examples of these are involved in mimicry, in which viral machinery "imitates" immunomodulatory elements in the vertebrate defense system. But recent findings have highlighted a vast and rapidly expanding array of other virologs that include many genes not typically found in viruses, such as those involved in translation, central carbon metabolism, cytoskeletal structure, nutrient transport, vesicular trafficking, and light harvesting. Unraveling the roles of virologs during infection as well as the evolutionary pathways through which complex functional repertoires are acquired by viruses are important frontiers at the forefront of giant virus research.

Homing in on the rare virosphere reveals the native host of giant viruses

Giant viruses (phylum Nucleocytoviricota) are globally distributed in aquatic ecosystems1,2. They play major roles as evolutionary drivers of eukaryotic plankton3 and regulators of global biogeochemical cycles4. Recent metagenomic studies have significantly expanded the known diversity of marine giant viruses1,5-7, but we still lack fundamental knowledge about their native hosts, thereby hindering our understanding of their lifecycle and ecological importance. Here, we aim to discover the native hosts of giant viruses using a novel, sensitive single-cell metatranscriptomic approach. By applying this approach to natural plankton communities, we unraveled an active viral infection of several giant viruses, from multiple lineages, and identified their native hosts. We identify a rare lineage of giant virus (Imitervirales-07) infecting a minute population of protists (class Katablepharidaceae) and revealed the prevalence of highly expressed viral-encoded cell-fate regulation genes in infected cells. Further examination of this host-virus dynamics in a temporal resolution suggested this giant virus controls its host population demise. Our results demonstrate how single-cell metatranscriptomics is a sensitive approach for pairing viruses with their authentic hosts and studying their ecological significance in a culture-independent manner in the marine environment.

Genomics discovery of giant fungal viruses from subsurface oceanic crustal fluids

The oceanic igneous crust is a vast reservoir for microbial life, dominated by diverse and active bacteria, archaea, and fungi. Archaeal and bacterial viruses were previously detected in oceanic crustal fluids at the Juan de Fuca Ridge (JdFR). Here we report the discovery of two eukaryotic Nucleocytoviricota genomes from the same crustal fluids by sorting and sequencing single virions. Both genomes have a tRNA^Tyr gene with an intron (20 bps) at the canonical position between nucleotide 37 and 38, a common feature in eukaryotic and archaeal tRNA genes with short introns (<100 bps), and fungal genes acquired through horizontal gene transfer (HGT) events. The dominance of Ascomycota fungi as the main eukaryotes in crustal fluids and the evidence for HGT point to these fungi as the putative hosts, making these the first putative fungi-Nucleocytoviricota specific association. Our study suggests active host-viral dynamics for the only eukaryotic group found in the subsurface oceanic crust and raises important questions about the impact of viral infection on the productivity and biogeochemical cycling in this ecosystem.

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.

Antigen Presentation and Autophagy in Teleost Adaptive Immunity

Infectious diseases are a burden for aquaculture. Antigen processing and presentation (APP) to the immune effector cells that fight pathogens is key in the adaptive immune response. At the core of the adaptive immunity that appeared in lower vertebrates during evolution are the variable genes encoding the major histocompatibility complex (MHC). MHC class I molecules mainly present peptides processed in the cytosol by the proteasome and transported to the cell surface of all cells through secretory compartments. Professional antigen-presenting cells (pAPC) also express MHC class II molecules, which normally present peptides processed from exogenous antigens through lysosomal pathways. Autophagy is an intracellular self-degradation process that is conserved in all eukaryotes and is induced by starvation to contribute to cellular homeostasis. Self-digestion during autophagy mainly occurs by the fusion of autophagosomes, which engulf portions of cytosol and fuse with lysosomes (macroautophagy) or assisted by chaperones (chaperone-mediated autophagy, CMA) that deliver proteins to lysosomes. Thus, during self-degradation, antigens can be processed to be presented by the MHC to immune effector cells, thus, linking autophagy to APP. This review is focused on the essential components of the APP that are conserved in teleost fish and the increasing evidence related to the modulation of APP and autophagy during pathogen infection.

Morphological and Genomic Features of the New Klosneuvirinae Isolate Fadolivirus IHUMI-VV54

Since the discovery of Mimivirus, viruses with large genomes encoding components of the translation machinery and other cellular processes have been described as belonging to the nucleocytoplasmic large DNA viruses. Recently, genome-resolved metagenomics led to the discovery of more than 40 viruses that have been grouped together in a proposed viral subfamily named Klosneuvirinae. Members of this group had genomes of up to 2.4Mb in size and featured an expanded array of translation system genes. Yet, despite the large diversity of the Klosneuvirinae in metagenomic data, there are currently only two isolates available. Here, we report the isolation of a novel giant virus known as Fadolivirus from an Algerian sewage site and provide morphological data throughout its replication cycle in amoeba and a detailed genomic characterization. The Fadolivirus genome, which is more than 1.5Mb in size, encodes 1,452 predicted proteins and phylogenetic analyses place this viral isolate as a near relative of the metagenome assembled Klosneuvirus and Indivirus. The genome encodes for 66 tRNAs, 23 aminoacyl-tRNA synthetases and a wide range of transcription factors, surpassing Klosneuvirus and other giant viruses. The Fadolivirus genome also encodes putative vacuolar-type proton pumps with the domains D and A, potentially constituting a virus-derived system for energy generation. The successful isolation of Fadolivirus will enable future hypothesis-driven experimental studies providing deeper insights into the biology of the Klosneuvirinae.

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.

Exploring Viral Diversity in a Gypsum Karst Lake Ecosystem Using Targeted Single-Cell Genomics

Little is known about the diversity and distribution of viruses infecting green sulfur bacteria (GSB) thriving in euxinic (sulfuric and anoxic) habitats, including gypsum karst lake ecosystems. In this study, we used targeted cell sorting combined with single-cell sequencing to gain insights into the gene content and genomic potential of viruses infecting sulfur-oxidizing bacteria Chlorobium clathratiforme, obtained from water samples collected during summer stratification in gypsum karst Lake Kirkilai (Lithuania). In total, 82 viral contigs were bioinformatically identified in 62 single amplified genomes (SAGs) of C. clathratiforme. The majority of viral gene and protein sequences showed little to no similarity with phage sequences in public databases, uncovering the vast diversity of previously undescribed GSB viruses. We observed a high level of lysogenization in the C. clathratiforme population, as 87% SAGs contained intact prophages. Among the thirty identified auxiliary metabolic genes (AMGs), two, thiosulfate sulfurtransferase (TST) and thioredoxin-dependent phosphoadenosine phosphosulfate (PAPS) reductase (cysH), were found to be involved in the oxidation of inorganic sulfur compounds, suggesting that viruses can influence the metabolism and cycling of this essential element. Finally, the analysis of CRISPR spacers retrieved from the consensus C. clathratiforme genome imply persistent and active virus–host interactions for several putative phages prevalent among C. clathratiforme SAGs. Overall, this study provides a glimpse into the diversity of phages associated with naturally occurring and highly abundant sulfur-oxidizing bacteria.

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.

Benchmarking of single-virus genomics: a new tool for uncovering the virosphere

Metagenomics and single-cell genomics have enabled the discovery of relevant uncultured microbes. Recently, single-virus genomics (SVG), although still in an incipient stage, has opened new avenues in viral ecology by allowing the sequencing of one single virus at a time. The investigation of methodological alternatives and optimization of existing procedures for SVG is paramount to deliver high-quality genomic data. We report a sequencing dataset of viral single-amplified genomes (vSAGs) from cultured and uncultured viruses obtained by applying different conditions in each SVG step, from viral preservation and novel whole-genome amplification (WGA) to sequencing platforms and genome assembly. Sequencing data showed that cryopreservation and mild fixation were compatible with WGA, although fresh samples delivered better genome quality data. The novel TruPrime WGA, based on primase-polymerase features, and WGA-X employing a thermostable phi29 polymerase, were proven to be with sufficient sensitivity in SVG. The Oxford Nanopore (ON) sequencing platform did not provide a significant improvement of vSAG assembly compared to Illumina alone. Finally, the SPAdes assembler performed the best. Overall, our results represent a valuable genomic dataset that will help to standardized and advance new tools in viral ecology.

Host Range and Coding Potential of Eukaryotic Giant Viruses

Giant viruses are a group of eukaryotic double-stranded DNA viruses with large virion and genome size that challenged the traditional view of virus. Newly isolated strains and sequenced genomes in the last two decades have substantially advanced our knowledge of their host diversity, gene functions, and evolutionary history. Giant viruses are now known to infect hosts from all major supergroups in the eukaryotic tree of life, which predominantly comprises microbial organisms. The seven well-recognized viral clades (taxonomic families) have drastically different host range. Mimiviridae and Phycodnaviridae, both with notable intrafamilial genome variation and high abundance in environmental samples, have members that infect the most diverse eukaryotic lineages. Laboratory experiments and comparative genomics have shed light on the unprecedented functional potential of giant viruses, encoding proteins for genetic information flow, energy metabolism, synthesis of biomolecules, membrane transport, and sensing that allow for sophisticated control of intracellular conditions and cell-environment interactions. Evolutionary genomics can illuminate how current and past hosts shape viral gene repertoires, although it becomes more obscure with divergent sequences and deep phylogenies. Continued works to characterize giant viruses from marine and other environments will further contribute to our understanding of their host range, coding potential, and virus-host coevolution.

Shrinking of repeating unit length in leucine-rich repeats from double-stranded DNA viruses

Leucine-rich repeats (LRRs) are present in over 563,000 proteins from viruses to eukaryotes. LRRs repeat in tandem and have been classified into fifteen classes in which the repeat unit lengths range from 20 to 29 residues. Most LRR proteins are involved in protein-protein or ligand interactions. The amount of genome sequence data from viruses is increasing rapidly, and although viral LRR proteins have been identified, a comprehensive sequence analysis has not yet been done, and their structures, functions, and evolution are still unknown. In the present study, we characterized viral LRRs by sequence analysis and identified over 600 LRR proteins from 89 virus species. Most of these proteins were from double-stranded DNA (dsDNA) viruses, including nucleocytoplasmic large dsDNA viruses (NCLDVs). We found that the repeating unit lengths of 11 types are one to five residues shorter than those of the seven known corresponding LRR classes. The repeating units of six types are 19 residues long and are thus the shortest among all LRRs. In addition, two of the LRR types are unique and have not been observed in bacteria, archae or eukaryotes. Conserved strongly hydrophobic residues such as Leu, Val or Ile in the consensus sequences are replaced by Cys with high frequency. Phylogenetic analysis indicated that horizontal gene transfer of some viral LRR genes had occurred between the virus and its host. We suggest that the shortening might contribute to the survival strategy of viruses. The present findings provide a new perspective on the origin and evolution of LRRs.

Host-virus interactions and defense mechanisms for giant viruses

Giant viruses, with virions larger than 200 nm and genomes larger than 340 kilobase pairs, modified the now outdated perception of the virosphere. With virions now reported reaching up to 1.5 μm in size and genomes of up to 2.5 Mb encoding components shared with cellular life forms, giant viruses exhibit a complexity similar to microbes, such as bacteria and archaea. Here, we review interactions of giant viruses with their hosts and defense strategies of giant viruses against their hosts and coinfecting microorganisms or virophages. We also searched by comparative genomics for homologies with proteins described or suspected to be involved in defense mechanisms. Our search reveals that natural immunity and apoptosis seem to be crucial components of the host defense against giant virus infection. Conversely, giant viruses possess methods of hijacking host functions to counteract cellular antiviral responses. In addition, giant viruses may encode other unique and complex pathways to manipulate the host machinery and eliminate other competing microorganisms. Notably, giant viruses have evolved defense mechanisms against their virophages and they might trigger defense systems against other viruses through sequence integration. We anticipate that comparative genomics may help identifying genes involved in defense strategies of both giant viruses and their hosts.

Single-virus genomics and beyond

Viruses are extremely diverse and modulate important biological and ecological processes globally. However, much of viral diversity remains uncultured and yet to be discovered. Several powerful culture-independent tools, in particular metagenomics, have substantially advanced virus discovery. Among those tools is single-virus genomics, which yields sequenced reference genomes from individual sorted virus particles without the need for cultivation. This new method complements virus culturing and metagenomic approaches and its advantages include targeted investigation of specific virus groups and investigation of genomic microdiversity within viral populations. In this Review, we provide a brief history of single-virus genomics, outline how this emergent method has facilitated advances in virus ecology and discuss its current limitations and future potential. Finally, we address how this method may synergistically intersect with other single-virus and single-cell approaches.

Giant virus diversity and host interactions through global metagenomics

Our current knowledge about nucleocytoplasmic large DNA viruses (NCLDVs) is largely derived from viral isolates that are co-cultivated with protists and algae. Here we reconstructed 2,074 NCLDV genomes from sampling sites across the globe by building on the rapidly increasing amount of publicly available metagenome data. This led to an 11-fold increase in phylogenetic diversity and a parallel 10-fold expansion in functional diversity. Analysis of 58,023 major capsid proteins from large and giant viruses using metagenomic data revealed the global distribution patterns and cosmopolitan nature of these viruses. The discovered viral genomes encoded a wide range of proteins with putative roles in photosynthesis and diverse substrate transport processes, indicating that host reprogramming is probably a common strategy in the NCLDVs. Furthermore, inferences of horizontal gene transfer connected viral lineages to diverse eukaryotic hosts. We anticipate that the global diversity of NCLDVs that we describe here will establish giant viruses-which are associated with most major eukaryotic lineages-as important players in ecosystems across Earth's biomes.

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

On the occurrence of cytochrome P450 in viruses

Genes encoding cytochrome P450 (CYP; P450) enzymes occur widely in the Archaea, Bacteria, and Eukarya, where they play important roles in metabolism of endogenous regulatory molecules and exogenous chemicals. We now report that genes for multiple and unique P450s occur commonly in giant viruses in the Mimiviridae, Pandoraviridae, and other families in the proposed order Megavirales. P450 genes were also identified in a herpesvirus (Ranid herpesvirus 3) and a phage (Mycobacterium phage Adler). The Adler phage P450 was classified as CYP102L1, and the crystal structure of the open form was solved at 2.5 Å. Genes encoding known redox partners for P450s (cytochrome P450 reductase, ferredoxin and ferredoxin reductase, and flavodoxin and flavodoxin reductase) were not found in any viral genome so far described, implying that host redox partners may drive viral P450 activities. Giant virus P450 proteins share no more than 25% identity with the P450 gene products we identified in Acanthamoeba castellanii, an amoeba host for many giant viruses. Thus, the origin of the unique P450 genes in giant viruses remains unknown. If giant virus P450 genes were acquired from a host, we suggest it could have been from an as yet unknown and possibly ancient host. These studies expand the horizon in the evolution and diversity of the enormously important P450 superfamily. Determining the origin and function of P450s in giant viruses may help to discern the origin of the giant viruses themselves.

Medusavirus, a Novel Large DNA Virus Discovered from Hot Spring Water

Recent discoveries of new large DNA viruses reveal high diversity in their morphologies, genetic repertoires, and replication strategies. Here, we report the novel features of medusavirus, a large DNA virus newly isolated from hot spring water in Japan. Medusavirus, with a diameter of 260 nm, shows a T=277 icosahedral capsid with unique spherical-headed spikes on its surface. It has a 381-kb genome encoding 461 putative proteins, 86 of which have their closest homologs in Acanthamoeba, whereas 279 (61%) are orphan genes. The virus lacks the genes encoding DNA topoisomerase II and RNA polymerase, showing that DNA replication takes place in the host nucleus, whereas the progeny virions are assembled in the cytoplasm. Furthermore, the medusavirus genome harbored genes for all five types of histones (H1, H2A, H2B, H3, and H4) and one DNA polymerase, which are phylogenetically placed at the root of the eukaryotic clades. In contrast, the host amoeba encoded many medusavirus homologs, including the major capsid protein. These facts strongly suggested that amoebae are indeed the most promising natural hosts of medusavirus, and that lateral gene transfers have taken place repeatedly and bidirectionally between the virus and its host since the early stage of their coevolution. Medusavirus reflects the traces of direct evolutionary interactions between the virus and eukaryotic hosts, which may be caused by sharing the DNA replication compartment and by evolutionarily long lasting virus-host relationships. Based on its unique morphological characteristics and phylogenomic relationships with other known large DNA viruses, we propose that medusavirus represents a new family, Medusaviridae IMPORTANCE We have isolated a new nucleocytoplasmic large DNA virus (NCLDV) from hot spring water in Japan, named medusavirus. This new NCLDV is phylogenetically placed at the root of the eukaryotic clades based on the phylogenies of several key genes, including that encoding DNA polymerase, and its genome surprisingly encodes the full set of histone homologs. Furthermore, its laboratory host, Acanthamoeba castellanii, encodes many medusavirus homologs in its genome, including the major capsid protein, suggesting that the amoeba is the genuine natural host from ancient times of this newly described virus and that lateral gene transfers have repeatedly occurred between the virus and amoeba. These results suggest that medusavirus is a unique NCLDV preserving ancient footprints of evolutionary interactions with its hosts, thus providing clues to elucidate the evolution of NCLDVs, eukaryotes, and virus-host interaction. Based on the dissimilarities with other known NCLDVs, we propose that medusavirus represents a new viral family, Medusaviridae.

Cryopreservation of Paramecium bursaria Chlorella Virus-1 during an active infection cycle of its host

Best practices in laboratory culture management often include cryopreservation of microbiota, but this can be challenging with some virus particles. By preserving viral isolates researchers can mitigate genetic drift and laboratory-induced selection, thereby maintaining genetically consistent strains between experiments. To this end, we developed a method to cryopreserve the model, green-alga infecting virus, Paramecium bursaria Chlorella virus 1 (PBCV-1). We explored cryotolerance of the infectivity of this virus particle, whereby freezing without cryoprotectants was found to maintain the highest infectivity (~2.5%). We then assessed the cryopreservation potential of PBCV-1 during an active infection cycle in its Chlorella variabilis NC64A host, and found that virus survivorship was highest (69.5 ± 16.5%) when the infected host is cryopreserved during mid-late stages of infection (i.e., coinciding with virion assembly). The most optimal condition for cryopreservation was observed at 240 minutes post-infection. Overall, utilizing the cell as a vehicle for viral cryopreservation resulted in 24.9–30.1 fold increases in PBCV-1 survival based on 95% confidence intervals of frozen virus particles and virus cryopreserved at 240 minutes post-infection. Given that cryoprotectants are often naturally produced by psychrophilic organisms, we suspect that cryopreservation of infected hosts may be a reliable mechanism for virus persistence in non-growth permitting circumstances in the environment, such as ancient permafrosts.

Bipartite Network Analysis of Gene Sharings in the Microbial World

Extensive microbial gene flows affect how we understand virology, microbiology, medical sciences, genetic modification, and evolutionary biology. Phylogenies only provide a narrow view of these gene flows: plasmids and viruses, lacking core genes, cannot be attached to cellular life on phylogenetic trees. Yet viruses and plasmids have a major impact on cellular evolution, affecting both the gene content and the dynamics of microbial communities. Using bipartite graphs that connect up to 149,000 clusters of homologous genes with 8,217 related and unrelated genomes, we can in particular show patterns of gene sharing that do not map neatly with the organismal phylogeny. Homologous genes are recycled by lateral gene transfer, and multiple copies of homologous genes are carried by otherwise completely unrelated (and possibly nested) genomes, that is, viruses, plasmids and prokaryotes. When a homologous gene is present on at least one plasmid or virus and at least one chromosome, a process of "gene externalization," affected by a postprocessed selected functional bias, takes place, especially in Bacteria. Bipartite graphs give us a view of vertical and horizontal gene flow beyond classic taxonomy on a single very large, analytically tractable, graph that goes beyond the cellular Web of Life.

Minimum Information about an Uncultivated Virus Genome (MIUViG)

We present an extension of the Minimum Information about any (x) Sequence (MIxS) standard for reporting sequences of uncultivated virus genomes. Minimum Information about an Uncultivated Virus Genome (MIUViG) standards were developed within the Genomic Standards Consortium framework and include virus origin, genome quality, genome annotation, taxonomic classification, biogeographic distribution and in silico host prediction. Community-wide adoption of MIUViG standards, which complement the Minimum Information about a Single Amplified Genome (MISAG) and Metagenome-Assembled Genome (MIMAG) standards for uncultivated bacteria and archaea, will improve the reporting of uncultivated virus genomes in public databases. In turn, this should enable more robust comparative studies and a systematic exploration of the global virosphere.

Seasonality Drives Microbial Community Structure, Shaping both Eukaryotic and Prokaryotic Host⁻Viral Relationships in an Arctic Marine Ecosystem

The Arctic marine environment experiences dramatic seasonal changes in light and nutrient availability. To investigate the influence of seasonality on Arctic marine virus communities, five research cruises to the west and north of Svalbard were conducted across one calendar year, collecting water from the surface to 1000 m in depth. We employed metabarcoding analysis of major capsid protein g23 and mcp genes in order to investigate T4-like myoviruses and large dsDNA viruses infecting prokaryotic and eukaryotic picophytoplankton, respectively. Microbial abundances were assessed using flow cytometry. Metabarcoding results demonstrated that seasonality was the key mediator shaping virus communities, whereas depth exerted a diversifying effect within seasonal virus assemblages. Viral diversity and virus-to-prokaryote ratios (VPRs) dropped sharply at the commencement of the spring bloom but increased across the season, ultimately achieving the highest levels during the winter season. These findings suggest that viral lysis may be an important process during the polar winter, when productivity is low. Furthermore, winter viral communities consisted of Operational Taxonomic Units (OTUs) distinct from those present during the spring-summer season. Our data provided a first insight into the diversity of viruses in a hitherto undescribed marine habitat characterized by extremes in light and productivity.

Hidden diversity of soil giant viruses

Known giant virus diversity is currently skewed towards viruses isolated from aquatic environments and cultivated in the laboratory. Here, we employ cultivation-independent metagenomics and mini-metagenomics on soils from the Harvard Forest, leading to the discovery of 16 novel giant viruses, chiefly recovered by mini-metagenomics. The candidate viruses greatly expand phylogenetic diversity of known giant viruses and either represented novel lineages or are affiliated with klosneuviruses, Cafeteria roenbergensis virus or tupanviruses. One assembled genome with a size of 2.4 Mb represents the largest currently known viral genome in the Mimiviridae, and others encode up to 80% orphan genes. In addition, we find more than 240 major capsid proteins encoded on unbinned metagenome fragments, further indicating that giant viruses are underexplored in soil ecosystems. The fact that most of these novel viruses evaded detection in bulk metagenomes suggests that mini-metagenomics could be a valuable approach to unearth viral giants.

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.

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.

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.

Deciphering the Human Virome with Single-Virus Genomics and Metagenomics

Single-cell genomics has unveiled the metabolic potential of dominant microbes inhabiting different environments, including the human body. The lack of genomic information for predominant microbes of the human body, such as bacteriophages, hinders our ability to answer fundamental questions about our viral communities. Here, we applied single-virus genomics (SVGs) to natural human salivary samples in combination with viral metagenomics to gain some insights into the viral community structure of the oral cavity. Saliva samples were processed for viral metagenomics (n = 15) and SVGs (n = 3). A total of 1328 uncultured single viruses were sorted by fluorescence-activated virus sorting followed by whole genome amplification. Sequencing of 24 viral single amplified genomes (vSAGs) showed that half of the vSAGs contained viral hallmark genes. Among those bona fide viruses, the uncultured single virus 92-C13 putatively infecting oral Streptococcus-like species was within the top ≈10 most abundant viruses in the oral virome. Viral gene network and viral metagenomics analyses of 439 oral viruses from cultures, metagenomics, and SVGs revealed that salivary viruses were tentatively structured into ≈200 major viral clusters, corresponding to approximately genus-level groupings. Data showed that none of the publicly available viral isolates, excepting an Actinomyces phage, were significantly abundant in the oral viromes. In addition, none of the obtained viral contigs and vSAGs from this study were present in all viromes. Overall, the data demonstrates that most viral isolates are not naturally abundant in saliva, and furthermore, the predominant viruses in the oral cavity are yet uncharacterized. Results suggest a variable, complex, and interpersonal viral profile. Finally, we demonstrated the power of SVGs in combination with viral metagenomics to unveil the genetic information of the uncultured viruses of the human virome.

RNase H As Gene Modifier, Driver of Evolution and Antiviral Defense

Retroviral infections are 'mini-symbiotic' events supplying recipient cells with sequences for viral replication, including the reverse transcriptase (RT) and ribonuclease H (RNase H). These proteins and other viral or cellular sequences can provide novel cellular functions including immune defense mechanisms. Their high error rate renders RT-RNases H drivers of evolutionary innovation. Integrated retroviruses and the related transposable elements (TEs) have existed for at least 150 million years, constitute up to 80% of eukaryotic genomes and are also present in prokaryotes. Endogenous retroviruses regulate host genes, have provided novel genes including the syncytins that mediate maternal-fetal immune tolerance and can be experimentally rendered infectious again. The RT and the RNase H are among the most ancient and abundant protein folds. RNases H may have evolved from ribozymes, related to viroids, early in the RNA world, forming ribosomes, RNA replicases and polymerases. Basic RNA-binding peptides enhance ribozyme catalysis. RT and ribozymes or RNases H are present today in bacterial group II introns, the precedents of TEs. Thousands of unique RTs and RNases H are present in eukaryotes, bacteria, and viruses. These enzymes mediate viral and cellular replication and antiviral defense in eukaryotes and prokaryotes, splicing, R-loop resolvation, DNA repair. RNase H-like activities are also required for the activity of small regulatory RNAs. The retroviral replication components share striking similarities with the RNA-induced silencing complex (RISC), the prokaryotic CRISPR-Cas machinery, eukaryotic V(D)J recombination and interferon systems. Viruses supply antiviral defense tools to cellular organisms. TEs are the evolutionary origin of siRNA and miRNA genes that, through RISC, counteract detrimental activities of TEs and chromosomal instability. Moreover, piRNAs, implicated in transgenerational inheritance, suppress TEs in germ cells. Thus, virtually all known immune defense mechanisms against viruses, phages, TEs, and extracellular pathogens require RNase H-like enzymes. Analogous to the prokaryotic CRISPR-Cas anti-phage defense possibly originating from TEs termed casposons, endogenized retroviruses ERVs and amplified TEs can be regarded as related forms of inheritable immunity in eukaryotes. This survey suggests that RNase H-like activities of retroviruses, TEs, and phages, have built up innate and adaptive immune systems throughout all domains of life.

A novel inducible prophage from the mycosphere inhabitant Paraburkholderia terrae BS437

Bacteriophages constitute key gene transfer agents in many bacteria. Specifically, they may confer gene mobility to Paraburkholderia spp. that dwells in soil and the mycosphere. In this study, we first screened mycosphere and bulk soils for phages able to produce plaques, however found these to be below detection. Then, prophage identification methods were applied to the genome sequences of the mycosphere-derived Paraburkholderia terrae strains BS001, BS007, BS110 and BS437, next to P. phytofirmans strains BS455, BIFAS53, J1U5 and PsJN. These analyses revealed all bacterial genomes to contain considerable amounts [up to 13.3%] of prophage-like sequences. One sequence predicted to encode a complete phage was found in the genome of P. terrae BS437. Using the inducing agent mitomycin C, we produced high-titered phage suspensions. These indeed encompassed the progeny of the identified prophage (denoted ɸ437), as evidenced using phage major capsid gene molecular detection. We obtained the full sequence of phage ɸ437, which, remarkably, had undergone a reshuffling of two large gene blocks. One predicted moron gene was found, and it is currently analyzed to understand the extent of its ecological significance for the host.