Mirusviruses link herpesviruses to giant viruses

Hot springs viruses at Yellowstone National Park have ancient origins and are adapted to thermophilic hosts

Geothermal springs house unicellular red algae in the class Cyanidiophyceae that dominate the microbial biomass at these sites. Little is known about host-virus interactions in these environments. We analyzed the virus community associated with red algal mats in three neighboring habitats (creek, endolithic, soil) at Lemonade Creek, Yellowstone National Park (YNP), USA. We find that despite proximity, each habitat houses a unique collection of viruses, with the giant viruses, Megaviricetes, dominant in all three. The early branching phylogenetic position of genes encoded on metagenome assembled virus genomes (vMAGs) suggests that the YNP lineages are of ancient origin and not due to multiple invasions from mesophilic habitats. The existence of genomic footprints of adaptation to thermophily in the vMAGs is consistent with this idea. The Cyanidiophyceae at geothermal sites originated ca. 1.5 Bya and are therefore relevant to understanding biotic interactions on the early Earth.

Symbiotic revolutions at the interface of genomics and microbiology

Symbiosis is an old idea with a contentious history. New genomic technologies and research paradigms are fueling a shift in some of its central tenets; we need to be humble and open-minded about what the data are telling us.

Xanthomonas Phage PBR31: Classifying the Unclassifiable

The ability of bacteriophages to destroy bacteria has made them the subject of extensive research. Interest in bacteriophages has recently increased due to the spread of drug-resistant bacteria, although genomic research has not kept pace with the growth of genomic data. Genomic analysis and, especially, the taxonomic description of bacteriophages are often difficult due to the peculiarities of the evolution of bacteriophages, which often includes the horizontal transfer of genes and genomic modules. The latter is particularly pronounced for temperate bacteriophages, which are capable of integration into the bacterial chromosome. Xanthomonas phage PBR31 is a temperate bacteriophage, which has been neither described nor classified previously, that infects the plant pathogen Xanthomonas campestris pv. campestris. Genomic analysis, including phylogenetic studies, indicated the separation of phage PBR31 from known classified bacteriophages, as well as its distant relationship with other temperate bacteriophages, including the Lederbervirus group. Bioinformatic analysis of proteins revealed distinctive features of PBR31, including the presence of a protein similar to the small subunit of D-family DNA polymerase and advanced lysis machinery. Taxonomic analysis showed the possibility of assigning phage PBR31 to a new taxon, although the complete taxonomic description of Xanthomonas phage PBR31 and other related bacteriophages is complicated by the complex evolutionary history of the formation of its genome. The general biological features of the PBR31 phage were analysed for the first time. Due to its presumably temperate lifestyle, there is doubt as to whether the PBR31 phage is appropriate for phage control purposes. Bioinformatics analysis, however, revealed the presence of cell wall-degrading enzymes that can be utilised for the treatment of bacterial infections.

Eukaryotic genomic data uncover an extensive host range of mirusviruses

A recent marine metagenomic study has revealed the existence of a novel group of viruses designated mirusviruses, which are proposed to form an evolutionary link between two realms of double-stranded DNA viruses, Varidnaviria and Duplodnaviria. Metagenomic data suggest that mirusviruses infect microeukaryotes in the photic layer of the ocean, but their host range remains largely unknown. In this study, we investigated the presence of mirusvirus marker genes in publicly available 1,901 eukaryotic genome assemblies, mainly derived from unicellular eukaryotes, to identify potential hosts of mirusviruses. Mirusvirus marker sequences were identified in 1,348 assemblies spanning 284 genera across eight supergroups of eukaryotes. The habitats of the putative mirusvirus hosts included not only marine but also other diverse environments. Among the major capsid protein (MCP) signals in the genome assemblies, we identified 85 sequences that showed high sequence and structural similarities to reference mirusvirus MCPs. A phylogenetic analysis of these sequences revealed their distant evolutionary relationships with the seven previously reported mirusvirus clades. Most of the scaffolds with these MCP sequences encoded multiple mirusvirus homologs, underscoring the impact of mirusviral infection on the evolution of the host genome. We also identified three circular mirusviral genomes within the genomic data of the oil producing thraustochytrid Schizochytrium sp. and the endolithic green alga Ostreobium quekettii. Overall, mirusviruses probably infect a wide spectrum of eukaryotes and are more diverse than previously reported.

Megataxonomy and global ecology of the virosphere

Nearly all organisms are hosts to multiple viruses that collectively appear to be the most abundant biological entities in the biosphere. With recent advances in metagenomics and metatranscriptomics, the known diversity of viruses substantially expanded. Comparative analysis of these viruses using advanced computational methods culminated in the reconstruction of the evolution of major groups of viruses and enabled the construction of a virus megataxonomy, which has been formally adopted by the International Committee on Taxonomy of Viruses. This comprehensive taxonomy consists of six virus realms, which are aspired to be monophyletic and assembled based on the conservation of hallmark proteins involved in capsid structure formation or genome replication. The viruses in different major taxa substantially differ in host range and accordingly in ecological niches. In this review article, we outline the latest developments in virus megataxonomy and the recent discoveries that will likely lead to reassessment of some major taxa, in particular, split of three of the current six realms into two or more independent realms. We then discuss the correspondence between virus taxonomy and the distribution of viruses among hosts and ecological niches, as well as the abundance of viruses versus cells in different habitats. The distribution of viruses across environments appears to be primarily determined by the host ranges, i.e. the virome is shaped by the composition of the biome in a given habitat, which itself is affected by abiotic factors.

Phylogenomic analysis expands the known repertoire of single-stranded DNA viruses in benthic zones of the South Indian Ocean

Single-stranded (ss) DNA viruses are ubiquitous and constitute some of the most diverse entities on Earth. Most studies have focused on ssDNA viruses from terrestrial environments resulting in a significant deficit in benthic ecosystems including aphotic zones of the South Indian Ocean (SIO). Here, we assess the diversity and phylogeny of ssDNA in deep waters of the SIO using a combination of established viral taxonomy tools and a Hidden Markov Model based approach. Replication initiator protein-associated (Rep) phylogenetic reconstruction and sequence similarity networks were used to show that the SIO hosts divergent and as yet unknown circular Rep-encoding ssDNA viruses. Several sequences appear to represent entirely novel families, expanding the repertoire of known ssDNA viruses. Results suggest that a small proportion of these viruses may be circular genetic elements, which may strongly influence the diversity of both eukaryotes and prokaryotes in the SIO. Taken together, our data show that the SIO harbours a diverse assortment of previously unknown ssDNA viruses. Due to their potential to infect a variety of hosts, these viruses may be crucial for marine nutrient recycling through their influence of the biological carbon pump.

Microbiology: The curious case of the mysterious mirusvirus

Long-read sequencing of a marine stramenopile genome yields a trove of insights into protist genomics and solves a 50-year-old viral mystery.

The protist Aurantiochytrium has universal subtelomeric rDNAs and is a host for mirusviruses

Viruses are the most abundant biological entities in the world's oceans, where they play important ecological and biogeochemical roles. Metagenomics is revealing new groups of eukaryotic viruses, although disconnected from known hosts. Among these are the recently described mirusviruses, which share some similarities with herpesviruses.1 50 years ago, "herpes-type" viral particles2 were found in a thraustochytrid member of the labyrinthulomycetes, a diverse group of abundant and ecologically important marine eukaryotes,3,4 but could not be further characterized by methods then available. Long-read sequencing has allowed us to connect the biology of mirusviruses and thraustochytrids. We sequenced the genome of the genetically tractable model thraustochytrid Aurantiochytrium limacinum ATCC MYA-1381 and found that its 26 linear chromosomes have an extraordinary configuration. Subtelomeric ribosomal DNAs (rDNAs) found at all chromosome ends are interspersed with long repeated sequence elements denoted as long repeated-telomere and rDNA spacers (LORE-TEARS). We identified two genomic elements that are related to mirusvirus genomes. The first is a ∼300-kbp episome (circular element 1 [CE1]) present at a high copy number. Strikingly, the second, distinct, mirusvirus-like element is integrated between two sets of rDNAs and LORE-TEARS at the left end of chromosome 15 (LE-Chr15). Similar to metagenomically derived mirusviruses, these putative A. limacinum mirusviruses have a virion module related to that of herpesviruses along with an informational module related to nucleocytoplasmic large DNA viruses (NCLDVs). CE1 and LE-Chr15 bear striking similarities to episomal and endogenous latent forms of herpesviruses, respectively, and open new avenues of research into marine virus-host interactions.

Structure-guided metagenome mining to tap microbial functional diversity

Scientists now have access to millions of accurate three-dimensional (3D) models of protein structures. How do we leverage 3D structural models to learn about microbial functions encoded in metagenomes? Here, we review recent developments using protein structural features to mine metagenomes from diverse environments ranging from the human gut to soil and ocean viromes. We compare 3D protein structural methods to characterize antibiotic resistance phenotypes, nutrient cycling, and host-drug-microbe interactions. Broadly, we encourage the scientific community to look beyond global sequence and structure alignments by considering fine-grained descriptors such as distance to ligand, active site, and tertiary interactions between amino acid residues scaling to microbiomes. Finally, we highlight structure-inspired approaches to chart new areas of microbial protein-coding sequence space.

Giant variations in giant virus genome packaging

Giant viruses (Nucleocytoviricota) have a largely conserved lifecycle, yet how they cram their large genomes into viral capsids is mostly unknown. The major capsid protein and the packaging ATPase (pATPase) comprise a highly conserved morphogenesis module in giant viruses, yet some giant viruses dispense with an icosahedral capsid, and others encode multiple versions of pATPases, including conjoined ATPase doublets, or encode none. Some giant viruses have acquired DNA-condensing proteins to compact their genomes, including sheath-like structures encasing folded DNA or densely packed viral nucleosomes that show a resemblance to eukaryotic nucleosomes at the telomeres. Here, we review what is known and unknown about these ATPases and condensing proteins, and place these variations in the context of viral lifecycles.

Automated classification of giant virus genomes using a random forest model built on trademark protein families

Viruses of the phylum Nucleocytoviricota , often referred to as "giant viruses," are prevalent in various environments around the globe and play significant roles in shaping eukaryotic diversity and activities in global ecosystems. Given the extensive phylogenetic diversity within this viral group and the highly complex composition of their genomes, taxonomic classification of giant viruses, particularly incomplete metagenome-assembled genomes (MAGs) can present a considerable challenge. Here we developed TIGTOG ( T axonomic Information of G iant viruses using T rademark O rthologous G roups), a machine learning-based approach to predict the taxonomic classification of novel giant virus MAGs based on profiles of protein family content. We applied a random forest algorithm to a training set of 1,531 quality-checked, phylogenetically diverse Nucleocytoviricota genomes using pre-selected sets of giant virus orthologous groups (GVOGs). The classification models were predictive of viral taxonomic assignments with a cross-validation accuracy of 99.6% to the order level and 97.3% to the family level. We found that no individual GVOGs or genome features significantly influenced the algorithm's performance or the models' predictions, indicating that classification predictions were based on a comprehensive genomic signature, which reduced the necessity of a fixed set of marker genes for taxonomic assigning purposes. Our classification models were validated with an independent test set of 823 giant virus genomes with varied genomic completeness and taxonomy and demonstrated an accuracy of 98.6% and 95.9% to the order and family level, respectively. Our results indicate that protein family profiles can be used to accurately classify large DNA viruses at different taxonomic levels and provide a fast and accurate method for the classification of giant viruses. This approach could easily be adapted to other viral groups.

Analysis of Pseudomonas aeruginosa Isolates from Patients with Cystic Fibrosis Revealed Novel Groups of Filamentous Bacteriophages

Pseudomonas aeruginosa is an opportunistic pathogen that can cause infections in humans, especially in hospital patients with compromised host defence mechanisms, including patients with cystic fibrosis. Filamentous bacteriophages represent a group of single-stranded DNA viruses infecting different bacteria, including P. aeruginosa and other human and animal pathogens; many of them can replicate when integrated into the bacterial chromosome. Filamentous bacteriophages can contribute to the virulence of P. aeruginosa and influence the course of the disease. There are just a few isolated and officially classified filamentous bacteriophages infecting P. aeruginosa, but genomic studies indicated the frequent occurrence of integrated prophages in many P. aeruginosa genomes. An analysis of sequenced genomes of P. aeruginosa isolated from upper respiratory tract (throat and nasal swabs) and sputum specimens collected from Russian patients with cystic fibrosis indicated a higher diversity of filamentous bacteriophages than first thought. A detailed analysis of predicted bacterial proteins revealed prophage regions representing the filamentous phages known to be quite distantly related to known phages. Genomic comparisons and phylogenetic studies enabled the proposal of several new taxonomic groups of filamentous bacteriophages.

Long-Read-Based Genome Assembly Reveals Numerous Endogenous Viral Elements in the Green Algal Bacterivore Cymbomonas tetramitiformis

The marine tetraflagellate Cymbomonas tetramitiformis has drawn attention as an early diverging green alga that uses a phago-mixotrophic mode of nutrition (i.e., the ability to derive nourishment from both photosynthesis and bacterial prey). The Cymbomonas nuclear genome was sequenced previously, but due to the exclusive use of short-read (Illumina) data, the assembly suffered from missing a large proportion of the genome's repeat regions. For this study, we generated Oxford Nanopore long-read and additional short-read Illumina data and performed a hybrid assembly that significantly improved the total assembly size and contiguity. Numerous endogenous viral elements were identified in the repeat regions of the new assembly. These include the complete genome of a giant Algavirales virus along with many genomes of integrated Polinton-like viruses (PLVs) from two groups: Gezel-like PLVs and a novel group of prasinophyte-specific PLVs. The integrated ∼400 kb genome of the giant Algavirales virus is the first account of the association of the uncultured viral family AG_03 with green algae. The complete PLV genomes from C. tetramitiformis ranged between 15 and 25 kb in length and showed a diverse gene content. In addition, heliorhodopsin gene-containing repeat elements of putative mirusvirus origin were identified. These results illustrate past (and possibly ongoing) multiple alga-virus interactions that accompanied the genome evolution of C. tetramitiformis.

A 1.5-Mb continuous endogenous viral region in the arbuscular mycorrhizal fungus Rhizophagus irregularis

Most fungal viruses are RNA viruses, and no double-stranded DNA virus that infects fungi is known to date. A recent study detected DNA polymerase genes that originated from large dsDNA viruses in the genomes of basal fungi, suggestive of the existence of dsDNA viruses capable of infecting fungi. In this study, we searched for viral infection signatures in chromosome-level genome assemblies of the arbuscular mycorrhizal fungus Rhizophagus irregularis. We identified a continuous 1.5-Mb putative viral region on a chromosome in R. irregularis strain 4401. Phylogenetic analyses revealed that the viral region is related to viruses in the family Asfarviridae of the phylum Nucleocytoviricota. This viral region was absent in the genomes of four other R. irregularis strains and had fewer signals of fungal transposable elements than the other genomic regions, suggesting a recent and single insertion of a large dsDNA viral genome in the genome of this fungal strain. We also incidentally identified viral-like sequences in the genome assembly of the sea slug Elysia marginata that are evolutionally close to the 1.5-Mb putative viral region. In conclusion, our findings provide strong evidence of the recent infection of the fungus by a dsDNA virus.

Genomic adaptation of giant viruses in polar oceans

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

Metagenome-derived virus-microbe ratios across ecosystems

It is generally assumed that viruses outnumber cells on Earth by at least tenfold. Virus-to-microbe ratios (VMR) are largely based on counts of fluorescently labelled virus-like particles. However, these exclude intracellular viruses and potentially include false positives (DNA-containing vesicles, gene-transfer agents, unspecifically stained inert particles). Here, we develop a metagenome-based VMR estimate (mVRM) that accounts for DNA viruses across all stages of their replication cycles (virion, intracellular lytic and lysogenic) by using normalised RPKM (reads per kilobase of gene sequence per million of mapped metagenome reads) counts of the major capsid protein (MCP) genes and cellular universal single-copy genes (USCGs) as proxies for virus and cell counts, respectively. After benchmarking this strategy using mock metagenomes with increasing VMR, we inferred mVMR across different biomes. To properly estimate mVMR in aquatic ecosystems, we generated metagenomes from co-occurring cellular and viral fractions (>50 kDa-200 µm size-range) in freshwater, seawater and solar saltern ponds (10 metagenomes, 2 control metaviromes). Viruses outnumbered cells in freshwater by ~13 fold and in plankton from marine and saline waters by ~2-4 fold. However, across an additional set of 121 diverse non-aquatic metagenomes including microbial mats, microbialites, soils, freshwater and marine sediments and metazoan-associated microbiomes, viruses, on average, outnumbered cells by barely two-fold. Although viruses likely are the most diverse biological entities on Earth, their global numbers might be closer to those of cells than previously estimated.

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.

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.

A critical analysis of the current state of virus taxonomy

Taxonomical classification has preceded evolutionary understanding. For that reason, taxonomy has become a battleground fueled by knowledge gaps, technical limitations, and a priorism. Here we assess the current state of the challenging field, focusing on fallacies that are common in viral classification. We emphasize that viruses are crucial contributors to the genomic and functional makeup of holobionts, organismal communities that behave as units of biological organization. Consequently, viruses cannot be considered taxonomic units because they challenge crucial concepts of organismality and individuality. Instead, they should be considered processes that integrate virions and their hosts into life cycles. Viruses harbor phylogenetic signatures of genetic transfer that compromise monophyly and the validity of deep taxonomic ranks. A focus on building phylogenetic networks using alignment-free methodologies and molecular structure can help mitigate the impasse, at least in part. Finally, structural phylogenomic analysis challenges the polyphyletic scenario of multiple viral origins adopted by virus taxonomy, defeating a polyphyletic origin and supporting instead an ancient cellular origin of viruses. We therefore, prompt abandoning deep ranks and urgently reevaluating the validity of taxonomic units and principles of virus classification.

The virome of the last eukaryotic common ancestor and eukaryogenesis

All extant eukaryotes descend from the last eukaryotic common ancestor (LECA), which is thought to have featured complex cellular organization. To gain insight into LECA biology and eukaryogenesis-the origin of the eukaryotic cell, which remains poorly understood-we reconstructed the LECA virus repertoire. We compiled an inventory of eukaryotic hosts of all major virus taxa and reconstructed the LECA virome by inferring the origins of these groups of viruses. The origin of the LECA virome can be traced back to a small set of bacterial-not archaeal-viruses. This provenance of the LECA virome is probably due to the bacterial origin of eukaryotic membranes, which is most compatible with two endosymbiosis events in a syntrophic model of eukaryogenesis. In the first endosymbiosis, a bacterial host engulfed an Asgard archaeon, preventing archaeal viruses from entry owing to a lack of archaeal virus receptors on the external membranes.

Assessing the biogeography of marine giant viruses in four oceanic transects

Viruses of the phylum Nucleocytoviricota are ubiquitous in ocean waters and play important roles in shaping the dynamics of marine ecosystems. In this study, we leveraged the bioGEOTRACES metagenomic dataset collected across the Atlantic and Pacific Oceans to investigate the biogeography of these viruses in marine environments. We identified 330 viral genomes, including 212 in the order Imitervirales and 54 in the order Algavirales. We found that most viruses appeared to be prevalent in shallow waters (<150 m), and that viruses of the Mesomimiviridae (Imitervirales) and Prasinoviridae (Algavirales) are by far the most abundant and diverse groups in our survey. Five mesomimiviruses and one prasinovirus are particularly widespread in oligotrophic waters; annotation of these genomes revealed common stress response systems, photosynthesis-associated genes, and oxidative stress modulation genes that may be key to their broad distribution in the pelagic ocean. We identified a latitudinal pattern in viral diversity in one cruise that traversed the North and South Atlantic Ocean, with viral diversity peaking at high latitudes of the northern hemisphere. Community analyses revealed three distinct Nucleocytoviricota communities across latitudes, categorized by latitudinal distance towards the equator. Our results contribute to the understanding of the biogeography of these viruses in marine systems.