Giant virus diversity and host interactions through global metagenomics

Horizontal transmission of functionally diverse transposons is a major source of new introns

Since the discovery of spliceosomal introns in eukaryotic genomes, the proximate molecular and evolutionary processes that generate new introns have remained a critical mystery. Specialized transposable elements (TEs), introners, are thought to be one of the major drivers of intron gain in diverse eukaryotes. However, the molecular mechanism(s) and evolutionary processes driving introner propagation within and between lineages remain elusive. Here, we analyze 8,716 genomes, revealing 1,093 introner families in 201 species spanning 1.7 billion years of evolution. Introners are derived from functionally diverse TEs including families of terminal-inverted-repeat DNA TEs, retrotransposons, cryptons, and helitrons as well as mobile elements with unknown molecular mechanisms. We identify eight cases where introners recently transferred between divergent host species and show that giant viruses that integrate into genomes may facilitate introner transfer across lineages. We propose that ongoing intron gain is primarily a consequence of TE activity in eukaryotes, thereby resolving a key mystery of genome structure evolution.

Hybrid sequencing reveals the genome of a Chrysochromulina parva virus and highlight its distinct replication strategy

Chrysochromulina parva (C. parva) is a eukaryotic freshwater haptophyte algae found in lakes and rivers worldwide. It is known to be infected by viruses, yet knowledge of the diversity and activity of these viruses is still very limited. Based on sequences of PCR-amplified DNA polymerase B (polB) gene fragments, Chrysochromulina parva virus BQ1 (CpV-BQ1) was the first known lytic agent of C. parva, and was considered a member of the virus family Phycodnaviridae, order Algavirales. However, the genome of a different C. parva-infecting virus (CpV-BQ2, or Tethysvirus ontarioense) from another virus family, the Mesomimiviridae, order Imitervirales, was the first sequenced. Here, we report the complete genome sequence of the putative phycodnavirus CpV-BQ1, accession PQ783904. The complete CpV-BQ1 genome sequence is 165,454 bp with a GC content of 32.32% and it encodes 193 open reading frames. Phylogenetic analyses of several virus hallmark genes including the polB, the late gene transcription factor (VLTF-3), and the putative A32-like virion packaging ATPase (Viral A32) all demonstrate that CpV-BQ1 is most closely related to other viruses in the phylum Megaviricetes within the order Algavirales, family Phycodnaviridae.

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.

Ion-conducting and gating molecular mechanisms of channelrhodopsin revealed by true-atomic-resolution structures of open and closed states

Channelrhodopsins (ChRs) have emerged as major optogenetics tools, particularly in neuroscience. Despite their importance, the molecular mechanism of ChR opening remains elusive. Moreover, all reported structures of ChRs correspond to either a closed or an early intermediate state and lack the necessary level of detail owing to the limited resolution. Here we present the structures of the closed and open states of a cation-conducting ChR, OLPVR1, from Organic Lake phycodnavirus, belonging to the family of viral ChRs solved at 1.1- and 1.3-Å resolution at physiologically relevant pH conditions (pH 8.0). OLPVR1 was expressed in Escherichia coli and crystallized using an in meso approach, and the structures were solved by X-ray crystallography. We also present the structure of the OLPVR1 protonated state at acidic pH (pH 2.5) at 1.4-Å resolution. Together, these three structures elucidate the molecular mechanisms of the channel's opening and permeability in detail. Extensive functional studies support the proposed mechanisms. Channel opening is controlled by isomerization of the retinal cofactor, triggering protonation of proton acceptors and deprotonation of proton donors located in the three gates of the channel. The E51 residue in the core of the central gate (similar to E90 of ChR2 from Chlamydomonas reinhardtii) plays a key role in the opening of the channel. E51 flips out of the gate and towards the proton acceptor D200 (D253 in ChR2 in C. reinhardtii), establishing a hydrogen bond between them. Despite differences in subfamilies of ChRs, they share a common gate-cavity architecture, suggesting that they could have similar general gating mechanisms. These results enabled us to design viral rhodopsin with improved properties for optogenetic applications. The structural data and mechanisms might also be helpful for better understanding other ChRs and their engineering.

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.

Genomic and structural insights into Jyvaskylavirus, the first giant virus isolated from Finland

Giant viruses of protists are a diverse and likely ubiquitous group of organisms. Here, we describe Jyvaskylavirus, the first giant virus isolated from Finland. This clade B marseillevirus was found in Acanthamoeba castellanii from a composting soil sample in Jyväskylä, Central Finland. Its genome shares similarities with other marseilleviruses. Helium ion microscopy and electron microscopy of infected cells unraveled stages of the Jyvaskylavirus life cycle. We reconstructed the Jyvaskylavirus particle to 6.3 Å resolution using cryo-electron microscopy. The ~2500 Å diameter virion displays structural similarities to other Marseilleviridae giant viruses. The capsid comprises of 9240 copies of the major capsid protein, encoded by open reading frame (ORF) 184, which possesses a double jellyroll fold arranged in trimers forming pseudo-hexameric capsomers. Below the capsid shell, the internal membrane vesicle encloses the genome. Through cross-structural and -sequence comparisons with other Marseilleviridae using AI-based software in model building and prediction, we elucidated ORF142 as the penton protein, which plugs the 12 vertices of the capsid. Five additional ORFs were identified, with models predicted and fitted into densities that either cap the capsomers externally or stabilize them internally. The isolation of Jyvaskylavirus suggests that these viruses may be widespread in the boreal environment and provide structural insights extendable to other marseilleviruses.

VITAP: a high precision tool for DNA and RNA viral classification based on meta-omic data

The rapid growth in the number of newly identified DNA and RNA viral sequences underscores the need for an accurate and comprehensive classification system for all viral realms at different taxonomic levels. Here, we establish the Viral Taxonomic Assignment Pipeline (VITAP), which addresses classification challenges by integrating alignment-based techniques with graphs, offering high precision in classifying both DNA and RNA viral sequences and providing confidence level for each taxonomic unit. This tool automatically updates its database in sync with the latest references from the International Committee on Taxonomy of Viruses (ICTV), efficiently classifying viral sequences as short as 1,000 base pairs to genus level. VITAP possesses good generalization capabilities, maintaining accuracy comparable to other pipelines while achieving higher annotation rates across most DNA and RNA viral phyla. Its application in deep-sea viromes has led to significant taxonomic updates, providing comprehensive diversity information of viruses from deep-sea. VITAP is available at https://github.com/DrKaiyangZheng/VITAP .

vClean: assessing virus sequence contamination in viral genomes

Recent advancements in viral metagenomics and single-virus genomics have improved our ability to obtain the draft genomes of environmental viruses. However, these methods can introduce virus sequence contaminations into viral genomes when short, fragmented partial sequences are present in the assembled contigs. These contaminations can lead to incorrect analyses; however, practical detection tools are lacking. In this study, we introduce vClean, a novel automated tool that detects contaminations in viral genomes. By applying machine learning to the nucleotide sequence features and gene patterns of the input viral genome, vClean could identify contaminations. Specifically, for tailed double-stranded DNA phages, we attempted accurate predictions by defining single-copy-like genes and counting their duplications. We evaluated the performance of vClean using simulated datasets derived from complete reference genomes, achieving a binary accuracy of 0.932. When vClean was applied to 4693 genomes of medium or higher quality derived from public ocean metagenomic data, 1604 genomes (34.2%) were identified as contaminated. We also demonstrated that vClean can detect contamination in single-virus genome data obtained from river water. vClean provides a new benchmark for quality control of environmental viral genomes and has the potential to become an essential tool for environmental viral genome analysis.

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.

RNA virus diversity highlights the potential biosecurity threat posed by Antarctic krill

Antarctic krill Euphausia superba, one of the most abundant species on the planet, is a keystone species of the Southern Ocean ecosystem. In the present study, we analyzed the RNA virome of Antarctic krill via metatranscription methods. The results showed that only 0.39% (49/12, 558) of the resultant unigenes could be assigned to known viral taxa, which were most similar to 17 known viruses, including nine invertebrate viruses, two vertebrate viruses, three protozoan viruses and three mycoviruses. However, most of the detected viruses possessed low amino acid similarity with counterparts in the viral databases. Penaeus vannamei picornavirus (PvPV; Family Picornaviridae) and covert mortality nodavirus (CMNV; Family Nodaviridae) were the two most abundant viruses in the Antarctic krill RNA virome. Notably, PvPV and CMNV are known pathogens to multiple aquatic animals according to epidemiological survey and exposure experiments, whereby PvPV positive krill caused clinical symptoms and histopathological lesions to P. vannamei and similarly, CMNV infection altered the swimming and feeding behavior of parent marine medaka Oryzias melastigma and caused tissue damage and even spinal curvature of the offspring. Results herein reveal, for the first time, the high abundance and taxonomic diversity of viruses in Antarctic krill while simultaneously highlighting the risk of an important virus reservoir to global aquaculture, and the potential impact on animals in the Antarctic ecosystem.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00270-w.

Co-occurrence patterns between Chlorophyta and nucleocytoplasmic large DNA virus in coastal ecosystem, South Korea

Nucleocytoplasmic large DNA viruses (NCLDVs) are known to infect phytoplankton and play a significant role in regulating their population dynamics. In this study, we aimed to investigate the co-occurrence patterns between phytoplankton and NCLDVs in the southern coastal ecosystem of South Korea. We collected seawater every month from March 2018 to December 2020 and analyzed the samples using Cytochrome c Oxidase subunit I metabarcoding and metagenomic analyses. Chlorophyta (36.08%) was the most abundant eukaryotic taxon, with Bathycoccaceae (58.25%) and Mamiellaceae (41.51%) being the most dominant families within Chlorophyta. Bathycoccaceae was dominant in winter, whereas Mamiellaceae was dominant in summer. In the NCLDV community, Phycodnaviridae (75.12%) was found to be the major family. The co-occurrence pattern of Phycodnaviridae showed a high correlation with Bathycoccaceae and Mamiellaceae, which is explained by the "boom-and-bust" concept. In particular, we predicted co-occurrence patterns between Bathycoccus prasinos and Prasnovirus, with known infectious relationships, and confirmed co-occurrence patterns between B. prasinos and Coccolithovirus and Micromonas pusilla and Prymnesiovirus, with unknown infectious relationships. These co-occurrence patterns between Chlorophyta and Phycodnaviridae provide valuable insights into the control of pico-sized primary production and the microbial loop of the coastal ecosystem.

The Detection of Mixed Infection with Canine Parvovirus, Canine Distemper Virus, and Rotavirus in Giant Pandas by Multiplex PCR

The well-being and subsistence of giant pandas, an endangered species with a limited distribution, are currently threatened by a number of viruses, including canine parvovirus (CPV-2), canine distemper virus (CDV), and giant panda rotavirus (GPRV). To allow for timely intervention upon viral infection, it is necessary to execute rapid and accurate diagnosis of potential mixed viral infections. In the present study, we developed and validated a multiplex PCR (mPCR) approach for the detection of CPV-2, CDV, and GPRV infections. The results indicate that the method could selectively amplify the three viruses with high sensitivity and specificity, which are necessary attributes in clinical settings. Utilizing the established method, (sub)clinical giant panda samples were examined, and CPV-2, CDV, and GPRV were found in 19.72% (43 out of 218), 7.34% (16 out of 218), and 6.42% (14 out of 218) of the samples, respectively. Noticeably, mixed infections of two or three viruses were common, and this was generally observed in CDV- or GPRV-positive samples. Meanwhile, mPCR results were further validated with sequencing and the phylogenetic analysis of full-length sequences of viral genes. Taken together, our study provides an approachable assay which enables the quick detection of the three viruses mentioned above, which will benefit clinical diagnosis and laboratory epidemiological-based investigations of the giant panda population.

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.

Novel High-Quality Amoeba Genomes Reveal Widespread Codon Usage Mismatch Between Giant Viruses and Their Hosts

The need for high-quality protist genomes has prevented in-depth computational and experimental studies of giant virus-host interactions. In addition, our current knowledge of host range is highly biased due to the few hosts used to isolate novel giant viruses. This study presents 6 high-quality amoeba genomes from known and potential giant virus hosts belonging to 2 distinct eukaryotic clades: Amoebozoa and Discoba. We employ their genomic data to investigate the predictability of giant virus host range. Using a combination of long- and short-read sequencing, we obtained highly contiguous and complete genomes of Acanthamoeba castellanii, Acanthamoeba griffini, Acanthamoeba terricola, Naegleria clarki, Vermamoeba vermiformis, and Willaertia magna, contributing to the collection of sequences for the eukaryotic tree of life. We found that the 6 amoebae have distinct codon usage patterns and that, contrary to other virus groups, giant viruses often have different and even opposite codon usage with their known hosts. Conversely, giant viruses with matching codon usage are frequently not known to infect or replicate in these hosts. Interestingly, analyses of integrated viral sequences in the amoeba host genomes reveal potential novel virus-host associations. Matching of codon usage preferences is often used to predict virus-host pairs. However, with the broad-scale analyses performed in this study, we demonstrate that codon usage alone appears to be a poor predictor of host range for giant viruses infecting amoeba. We discuss the potential strategies that giant viruses employ to ensure high viral fitness in nonmatching hosts. Moreover, this study emphasizes the need for more high-quality protist genomes. Finally, the amoeba genomes presented in this study set the stage for future experimental studies to better understand how giant viruses interact with different host species.

Pacmanvirus isolated from the Lost City hydrothermal field extends the concept of transpoviron beyond the family Mimiviridae

The microbial sampling of submarine hydrothermal vents remains challenging, with even fewer studies focused on viruses. Here we report what is to our knowledge the first isolation of a eukaryotic virus from the Lost City hydrothermal field, by co-culture with the laboratory host Acanthamoeba castellanii. This virus, named pacmanvirus lostcity, is closely related to previously isolated pacmanviruses (strains A23 and S19), clustering in a divergent clade within the long-established family Asfarviridae. The icosahedral particles of this virus are 200 nm in diameter, with an electron-dense core surrounded by an inner membrane. The viral genome of 395 708 bp (33% G + C) has been predicted to encode 473 proteins. However, besides these standard properties, pacmanvirus lostcity was found to be associated with a new type of selfish genetic element, 7 kb in length, whose architecture and gene content are reminiscent of those of transpovirons, hitherto specific to the family Mimiviridae. As in previously described transpovirons, this selfishg genetic element propagates as an episome within its host virus particles and exhibits partial recombination with its genome. In addition, an unrelated episome with a length of 2 kb was also found to be associated with pacmanvirus lostcity. Together, the transpoviron and the 2-kb episome might participate in exchanges between pacmanviruses and other DNA virus families. It remains to be elucidated if the presence of these mobile genetic elements is restricted to pacmanviruses or was simply overlooked in other members of the Asfarviridae.

How Viruses Shape Microbial Plankton Microdiversity

One major conundrum of modern microbiology is the large pangenome (gene pool) present in microbes, which is much larger than those found in complex organisms such as humans. Here, we argue that this diversity of gene pools carried by different strains is maintained largely due to the control exercised by viral predation. Viruses maintain a high strain diversity through time that we describe as constant-diversity equilibrium, preventing the hoarding of resources by specific clones. Thus, viruses facilitate the release and degradation of dissolved organic matter in the ocean, which may lead to better ecosystem functioning by linking top-down to bottom-up control. By maintaining this equilibrium, viruses act as a key element of the adaptation of marine microbes to their environment and likely behave as a single evolutionary unit.

Viral niche-partitioning: comparative genomics of giant viruses across environmental gradients in a high Arctic freshwater-saltwater lake

Giant viruses (GVs; Nucleocytoviricota) impact the biology and ecology of a wide range of eukaryotic hosts, with implications for global biogeochemical cycles. Here, we investigated GV niche separation in highly stratified Lake A at the northern coast of Ellesmere Island, Nunavut, Canada. This lake is composed of a layer of ice-covered freshwater that overlies saltwater derived from the ancient Arctic Ocean, and it therefore provides a broad gradient of environmental conditions and ecological habitats, each with a distinct protist community and rich assemblages of associated GVs. The upper layer (mixolimnion) had measurable light and oxygen, and contained diverse GVs linked to photosynthetic protists, indicating adaptation to surface biotic and abiotic conditions. In contrast, the saline lower layer (monimolimnion), lacking oxygen and light, hosted GVs associated with predicted heterotrophic protists, some of which are known for a predatory lifestyle, and with several viral genes suggesting adaptation to deep-water anaerobic conditions. Our observations underscore the coupling between physical and chemical gradients, microeukaryotes and their associated GVs in Lake A, and provide insight into the potential for GVs to directly and indirectly impact host metabolism. There were similarities between the genetic composition of GVs and the metabolic processes of their potential hosts, implying co-evolution and niche-adaptation within the lake habitats. Notably, we found a greater presence of viral rhodopsins in deeper water layers, suggesting an evolutionary relationship with potential hosts capable of supplementing their energetic needs to thrive in low energy, anoxic conditions.

Conservative taxonomy and quality assessment of giant virus genomes with GVClass

Large double-stranded DNA viruses of the phylum Nucleocytoviricota (Giant Viruses; GVs) are the largest known viruses, infecting various eukaryotic hosts, particularly protists and algae. These viruses impact biogeochemical cycles and host genome evolution but are challenging to identify and classify due to their complex genomes. We present GVClass, a tool for identifying giant viruses in sequence data, providing taxonomic assignments, and estimating genome completeness and contamination. GVClass employs optimized gene calling and a conservative approach using consensus single-protein phylogenies for robust taxonomic classification, relying on highly conserved orthologous groups. Benchmarking demonstrates over 90% accuracy at the genus-level and >99% at higher taxonomic ranks. GVClass addresses classification challenges and is available as a standalone tool and integrated into the Integrated Microbial Genomes/Virus database (IMG/VR).

Complete genomes of Asgard archaea reveal diverse integrated and mobile genetic elements

Asgard archaea are of great interest as the progenitors of Eukaryotes, but little is known about the mobile genetic elements (MGEs) that may shape their ongoing evolution. Here, we describe MGEs that replicate in Atabeyarchaeia, a wetland Asgard archaea lineage represented by two complete genomes. We used soil depth-resolved population metagenomic data sets to track 18 MGEs for which genome structures were defined and precise chromosome integration sites could be identified for confident host linkage. Additionally, we identified a complete 20.67 kbp circular plasmid and two family-level groups of viruses linked to Atabeyarchaeia, via CRISPR spacer targeting. Closely related 40 kbp viruses possess a hypervariable genomic region encoding combinations of specific genes for small cysteine-rich proteins structurally similar to restriction-homing endonucleases. One 10.9 kbp integrative conjugative element (ICE) integrates genomically into the Atabeyarchaeum deiterrae-1 chromosome and has a 2.5 kbp circularizable element integrated within it. The 10.9 kbp ICE encodes an expressed Type IIG restriction-modification system with a sequence specificity matching an active methylation motif identified by Pacific Biosciences (PacBio) high-accuracy long-read (HiFi) metagenomic sequencing. Restriction-modification of Atabeyarchaeia differs from that of another coexisting Asgard archaea, Freyarchaeia, which has few identified MGEs but possesses diverse defense mechanisms, including DISARM and Hachiman, not found in Atabeyarchaeia. Overall, defense systems and methylation mechanisms of Asgard archaea likely modulate their interactions with MGEs, and integration/excision and copy number variation of MGEs in turn enable host genetic versatility.

Seasonal dynamics and diversity of Antarctic marine viruses reveal a novel viral seascape

The Southern Ocean microbial ecosystem, with its pronounced seasonal shifts, is vulnerable to the impacts of climate change. Since viruses are key modulators of microbial abundance, diversity, and evolution, we need a better understanding of the effects of seasonality on the viruses in this region. Our comprehensive exploration of DNA viral diversity in the Southern Ocean reveals a unique and largely uncharted viral landscape, of which 75% was previously unidentified in other oceanic areas. We uncover novel viral taxa at high taxonomic ranks, expanding our understanding of crassphage, polinton-like virus, and virophage diversity. Nucleocytoviricota viruses represent an abundant and diverse group of Antarctic viruses, highlighting their potential as important regulators of phytoplankton population dynamics. Our temporal analysis reveals complex seasonal patterns in marine viral communities (bacteriophages, eukaryotic viruses) which underscores the apparent interactions with their microbial hosts, whilst deepening our understanding of their roles in the world's most sensitive and rapidly changing ecosystem.

Eukaryotic viruses encode the ribosomal protein eL40

Viruses in the phylum Nucleocytoviricota are large, complex and have an exceptionally diverse metabolic repertoire. Some encode hundreds of products involved in the translation of mRNA into protein, but none was known to encode any of the proteins in ribosomes, the central engines of translation. With the discovery of the eL40 gene in FloV-SA2, we report the first example of a eukaryotic virus encoding a ribosomal protein and show that this gene is also present and expressed in other uncultivated marine giant viruses. FloV-SA2 also encodes a "group II" viral rhodopsin, a viral light-activated protein of unknown function previously only reported in metagenomes. FloV-SA2 is thus a valuable model system for investigating new mechanisms by which viruses manipulate eukaryotic cell metabolism.

Robust optogenetic inhibition with red-light-sensitive anion-conducting channelrhodopsins

Channelrhodopsins (ChRs) are light-gated ion channels widely used to optically activate or silence selected electrogenic cells, such as individual brain neurons. Here, we describe identifying and characterizing a set of anion-conducting ChRs (ACRs) from diverse taxa and representing various branches of the ChR phylogenetic tree. The Mantoniella squamata ACR (MsACR1) showed high sensitivity to yellow-green light (λmax at 555 nm) and was further engineered for optogenetic applications. A single amino-acid substitution that mimicked red-light-sensitive rhodopsins like Chrimson shifted the photosensitivity 20 nm toward red light and accelerated photocurrent kinetics. Hence, it was named red and accelerated ACR, raACR. Both wild-type and mutant are capable optical silencers at low light intensities in mouse neurons in vitro and in vivo, while raACR offers a higher temporal resolution.

Complex transcriptional regulations of a hyperparasitic quadripartite system in giant viruses infecting protists

Hyperparasitism is a common pattern in nature that is not limited to cellular organisms. Giant viruses infecting protists can be hyperparasitized by smaller ones named virophages. In addition, both may carry episomal DNA molecules known as transpovirons in their particles. They all share transcriptional regulatory elements that dictate the expression of their genes within viral factories built by giant viruses in the host cytoplasm. This suggests the existence of interactions between their respective transcriptional networks. Here we investigated Acanthamoeba castellanii cells infected by a giant virus (megavirus chilensis), and coinfected with a virophage (zamilon vitis) and/or a transpoviron (megavirus vitis transpoviron). Infectious cycles were monitored through time-course RNA sequencing to decipher the transcriptional program of each partner and its impact on the gene expression of the others. We found highly diverse transcriptional responses. While the giant virus drastically reshaped the host cell transcriptome, the transpoviron had no effect on the gene expression of any of the players. In contrast, the virophage strongly modified the giant virus gene expression, albeit transiently, without altering the protein composition of mature viral particles. The virophage also induced the overexpression of transpoviron genes, likely through the indirect upregulation of giant virus-encoded transcription factors. Together, these analyses document the intricated transcriptionally regulated networks taking place in the infected cell.

Differentiation of the Anammox core microbiome: Unraveling the evolutionary impetus of scalable gene flow

Anaerobic ammonium oxidation bacteria (AAOB), distinguished by their unique autotrophic nitrogen metabolism, hold pivotal positions in the global nitrogen cycle and environmental biotechnologies. However, the ecophysiology and evolution of AAOB remain poorly understood, attributed to the absence of monocultures. Hence, a comprehensive elucidation of the AAOB-dominated core microbiome, anammox core, is imperative to further completing the theory of engineered nitrogen removal and ecological roles of anammox. Performing taxonomic and phylogenetic analyses on collected genome repertoires, we show here that Candidatus Brocadia and Candidatus Kuenenia possesses a more compact core than Candidatus Jettenia, which partly explains why the latter has a less common ecological presence. Evidence of gene flow is particularly striking in functions related to biosynthesis and oxygen detoxification, underscoring the evolutionary forces driving lineage and core differentiation. Furthermore, CRISPR spacer traceback of the AAOB metagenome-assembled genomes (MAGs) reveals a series of genetic traces for the concealed phages. By reconceptualizing the functional divergence of AAOB with the historical role of phages, we ultimately propose a coevolutionary framework to understand the evolutionary trajectory of anammox microecology. The discoveries provided in this study offer new insights into understanding the evolution of AAOB and the ecology of anammox.

Prophage-encoded antibiotic resistance genes are enriched in human-impacted environments

The spread of antibiotic resistance genes (ARGs) poses a substantial threat to human health. Phage-mediated transduction could exacerbate ARG transmission. While several case studies exist, it is yet unclear to what extent phages encode and mobilize ARGs at the global scale and whether human impacts play a role in this across different habitats. Here, we combine 38,605 bacterial genomes, 1432 metagenomes, and 1186 metatranscriptomes across 12 contrasting habitats to explore the distribution of prophages and their cargo ARGs in natural and human-impacted environments. Worldwide, we observe a significant increase in the abundance, diversity, and activity of prophage-encoded ARGs in human-impacted habitats linked with relatively higher risk of past antibiotic exposure. This effect was driven by phage-encoded cargo ARGs that could be mobilized to provide increased resistance in heterologous E. coli host for a subset of analyzed strains. Our findings suggest that human activities have altered bacteria-phage interactions, enriching ARGs in prophages and making ARGs more mobile across habitats globally.

Modern microbiology: Embracing complexity through integration across scales

Microbes were the only form of life on Earth for most of its history, and they still account for the vast majority of life's diversity. They convert rocks to soil, produce much of the oxygen we breathe, remediate our sewage, and sustain agriculture. Microbes are vital to planetary health as they maintain biogeochemical cycles that produce and consume major greenhouse gases and support large food webs. Modern microbiologists analyze nucleic acids, proteins, and metabolites; leverage sophisticated genetic tools, software, and bioinformatic algorithms; and process and integrate complex and heterogeneous datasets so that microbial systems may be harnessed to address contemporary challenges in health, the environment, and basic science. Here, we consider an inevitably incomplete list of emergent themes in our discipline and highlight those that we recognize as the archetypes of its modern era that aim to address the most pressing problems of the 21st century.

Time of day of infection shapes development of a eukaryotic algal-Nucleocytoviricota virocell

Aureococcus anophagefferens forms a model host-virus system with the "giant virus" Kratosvirus quantuckense. Studies to define its ribocell (uninfected) and virocell (virus-infected) forms are needed as these states co-occur during algal blooms. Previously, a link between light-derived energy, virus particle production, and virocell formation was noted. We explored how the time of day (morning, midday, or late day) of virus-host contact shaped virocell ontogeny. In parallel, we explored the dependence on light-derived energy in this mixotrophic plankter by inhibiting photosystem II, testing the role of heterotrophic energy in infection dynamics. Using flow cytometry and photochemical assessments, we examined the physiology of infected cells and controls, and estimated virus particle production. We observed differences between ribocell and virocell response to treatments, including reductions in virus particle production during reduced light duration) and PSII inhibition (i.e. "forced heterotrophy"). This work demonstrates the importance of light in shaping the fate of infected cells and provides insight into factors that constrain in situ blooms. Most significantly, we show that time of the solar day when a virus and host come into contact influences viral particle production, and therefore bloom dynamics; a factor that needs to be considered in bloom modeling work.

Giant viruses as reservoirs of antibiotic resistance genes

Nucleocytoplasmic large DNA viruses (NCLDVs; also called giant viruses), constituting the phylum Nucleocytoviricota, can infect a wide range of eukaryotes and exchange genetic material with not only their hosts but also prokaryotes and phages. A few NCLDVs were reported to encode genes conferring resistance to beta‑lactam, trimethoprim, or pyrimethamine, suggesting that they are potential vehicles for the transmission of antibiotic resistance genes (ARGs) in the biome. However, the incidence of ARGs across the phylum Nucleocytoviricota, their evolutionary characteristics, their dissemination potential, and their association with virulence factors remain unexplored. Here, we systematically investigated ARGs of 1416 NCLDV genomes including those of almost all currently available cultured isolates and high-quality metagenome-assembled genomes from diverse habitats across the globe. We reveal that 39.5% of them carry ARGs, which is approximately 37 times higher than that for phage genomes. A total of 12 ARG types are encoded by NCLDVs. Phylogenies of the three most abundant NCLDV-encoded ARGs hint that NCLDVs acquire ARGs from not only eukaryotes but also prokaryotes and phages. Two NCLDV-encoded trimethoprim resistance genes are demonstrated to confer trimethoprim resistance in Escherichia coli. The presence of ARGs in NCLDV genomes is significantly correlated with mobile genetic elements and virulence factors.

Prokaryotic-virus-encoded auxiliary metabolic genes throughout the global oceans

BACKGROUND

Prokaryotic microbes have impacted marine biogeochemical cycles for billions of years. Viruses also impact these cycles, through lysis, horizontal gene transfer, and encoding and expressing genes that contribute to metabolic reprogramming of prokaryotic cells. While this impact is difficult to quantify in nature, we hypothesized that it can be examined by surveying virus-encoded auxiliary metabolic genes (AMGs) and assessing their ecological context.

RESULTS

We systematically developed a global ocean AMG catalog by integrating previously described and newly identified AMGs and then placed this catalog into ecological and metabolic contexts relevant to ocean biogeochemistry. From 7.6 terabases of Tara Oceans paired prokaryote- and virus-enriched metagenomic sequence data, we increased known ocean virus populations to 579,904 (up 16%). From these virus populations, we then conservatively identified 86,913 AMGs that grouped into 22,779 sequence-based gene clusters, 7248 (~ 32%) of which were not previously reported. Using our catalog and modeled data from mock communities, we estimate that ~ 19% of ocean virus populations carry at least one AMG. To understand AMGs in their metabolic context, we identified 340 metabolic pathways encoded by ocean microbes and showed that AMGs map to 128 of them. Furthermore, we identified metabolic "hot spots" targeted by virus AMGs, including nine pathways where most steps (≥ 0.75) were AMG-targeted (involved in carbohydrate, amino acid, fatty acid, and nucleotide metabolism), as well as other pathways where virus-encoded AMGs outnumbered cellular homologs (involved in lipid A phosphates, phosphatidylethanolamine, creatine biosynthesis, phosphoribosylamine-glycine ligase, and carbamoyl-phosphate synthase pathways).

CONCLUSIONS

Together, this systematically curated, global ocean AMG catalog and analyses provide a valuable resource and foundational observations to understand the role of viruses in modulating global ocean metabolisms and their biogeochemical implications. Video Abstract.

Phylogenetic proximity drives temporal succession of marine giant viruses in a five-year metagenomic time-series

Nucleocytoplasmic Large DNA Viruses (NCLDVs, also called giant viruses) are widespread in marine systems and infect a broad range of microbial eukaryotes (protists). Recent biogeographic work has provided global snapshots of NCLDV diversity and community composition across the world's oceans, yet little information exists about the guiding 'rules' underpinning their community dynamics over time. We leveraged a five-year monthly metagenomic time-series to quantify the community composition of NCLDVs off the coast of Southern California and characterize these populations' temporal dynamics. NCLDVs were dominated by Algavirales (Phycodnaviruses, 59%) and Imitervirales (Mimiviruses, 36%). We identified clusters of NCLDVs with distinct classes of seasonal and non-seasonal temporal dynamics. Overall, NCLDV population abundances were often highly dynamic with a strong seasonal signal. The Imitervirales group had highest relative abundance in the more oligotrophic late summer and fall, while Algavirales did so in winter. Generally, closely related strains had similar temporal dynamics, suggesting that evolutionary history is a key driver of the temporal niche of marine NCLDVs. However, a few closely-related strains had drastically different seasonal dynamics, suggesting that while phylogenetic proximity often indicates ecological similarity, occasionally phenology can shift rapidly, possibly due to host-switching. Finally, we identified distinct functional content and possible host interactions of two major NCLDV orders-including connections of Imitervirales with primary producers like the diatom Chaetoceros and widespread marine grazers like Paraphysomonas and Spirotrichea ciliates. Together, our results reveal key insights on season-specific effect of phylogenetically distinct giant virus communities on marine protist metabolism, biogeochemical fluxes and carbon cycling.

Unveiling the unknown viral world in groundwater

Viruses as the prevailing biological entities are poorly understood in underground realms. Here, we establish the first metagenomic Groundwater Virome Catalogue (GWVC) comprising 280,420 viral species ( ≥ 5 kb) detected from 607 monitored wells in seven geo-environmental zones throughout China. In expanding ~10-fold the global portfolio of known groundwater viruses, we uncover over 99% novel viruses and about 95% novel viral clusters. By linking viruses to hosts from 119 prokaryotic phyla, we double the number of microbial phyla known to be virus-infected in groundwater. As keystone ultrasmall symbionts in aquifers, CPR bacteria and DPANN archaea are susceptible to virulent viruses. Certain complete CPR viruses even likely infect non-CPR bacteria, while partial CPR/DPANN viruses harbor cell-surface modification genes that assist symbiont cell adhesion to free-living microbes. This study reveals the unknown viral world and auxiliary metabolism associated with methane, nitrogen, sulfur, and phosphorus cycling in groundwater, and highlights the importance of subsurface virosphere in viral ecology.

Complex Genomes of Early Nucleocytoviruses Revealed by Ancient Origins of Viral Aminoacyl-tRNA Synthetases

Aminoacyl-tRNA synthetases (aaRSs), also known as tRNA ligases, are essential enzymes in translation. Owing to their functional essentiality, these enzymes are conserved in all domains of life and used as informative markers to trace the evolutionary history of cellular organisms. Unlike cellular organisms, viruses generally lack aaRSs because of their obligate parasitic nature, but several large and giant DNA viruses in the phylum Nucleocytoviricota encode aaRSs in their genomes. The discovery of viral aaRSs led to the idea that the phylogenetic analysis of aaRSs can shed light on ancient viral evolution. However, conflicting results have been reported from previous phylogenetic studies: one posited that nucleocytoviruses recently acquired their aaRSs from their host eukaryotes, while another hypothesized that the viral aaRSs have ancient origins. Here, we investigated 4,168 nucleocytovirus genomes, including metagenome-assembled genomes (MAGs) derived from large-scale metagenomic studies. In total, we identified 780 viral aaRS sequences in 273 viral genomes. We generated and examined phylogenetic trees of these aaRSs with a large set of cellular sequences to trace evolutionary relationships between viral and cellular aaRSs. The analyses suggest that the origins of some viral aaRSs predate the last common eukaryotic ancestor. Inside viral aaRS clades, we identify intricate evolutionary trajectories of viral aaRSs with horizontal transfers, losses, and displacements. Overall, these results suggest that ancestral nucleocytoviruses already developed complex genomes with an expanded set of aaRSs in the proto-eukaryotic era.

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.

Identification of mobile genetic elements with geNomad

Identifying and characterizing mobile genetic elements in sequencing data is essential for understanding their diversity, ecology, biotechnological applications and impact on public health. Here we introduce geNomad, a classification and annotation framework that combines information from gene content and a deep neural network to identify sequences of plasmids and viruses. geNomad uses a dataset of more than 200,000 marker protein profiles to provide functional gene annotation and taxonomic assignment of viral genomes. Using a conditional random field model, geNomad also detects proviruses integrated into host genomes with high precision. In benchmarks, geNomad achieved high classification performance for diverse plasmids and viruses (Matthews correlation coefficient of 77.8% and 95.3%, respectively), substantially outperforming other tools. Leveraging geNomad's speed and scalability, we processed over 2.7 trillion base pairs of sequencing data, leading to the discovery of millions of viruses and plasmids that are available through the IMG/VR and IMG/PR databases. geNomad is available at https://portal.nersc.gov/genomad .

Pathogen dynamics and discovery of novel viruses and enzymes by deep nucleic acid sequencing of wastewater

Wastewater contains an extensive reservoir of genetic information, yet largely unexplored. Here, we analyzed by high-throughput sequencing total nucleic acids extracted from wastewater samples collected during a 17 month-period in Berlin, Germany. By integrating global wastewater datasets and applying a novel computational approach to accurately identify viral strains within sewage RNA-sequencing data, we demonstrated the emergence and global dissemination of a specific astrovirus strain. Astrovirus abundance and sequence variation mirrored temporal and spatial patterns of infection, potentially serving as footprints of specific timeframes and geographical locations. Additionally, we revealed more than 100,000 sequence contigs likely originating from novel viral species, exhibiting distinct profiles in total RNA and DNA datasets and including undescribed bunyaviruses and parvoviruses. Finally, we identified thousands of new CRISPR-associated protein sequences, including Transposase B (TnpB), a class of compact, RNA-guided DNA editing enzymes. Collectively, our findings underscore the potential of high-throughput sequencing of total nucleic acids derived from wastewater for a broad range of applications.

Novel Z-DNA binding domains in giant viruses

Z-nucleic acid structures play vital roles in cellular processes and have implications in innate immunity due to their recognition by Zα domains containing proteins (Z-DNA/Z-RNA binding proteins, ZBPs). Although Zα domains have been identified in six proteins, including viral E3L, ORF112, and I73R, as well as, cellular ADAR1, ZBP1, and PKZ, their prevalence across living organisms remains largely unexplored. In this study, we introduce a computational approach to predict Zα domains, leading to the revelation of previously unidentified Zα domain-containing proteins in eukaryotic organisms, including non-metazoan species. Our findings encompass the discovery of new ZBPs in previously unexplored giant viruses, members of the Nucleocytoviricota phylum. Through experimental validation, we confirm the Zα functionality of select proteins, establishing their capability to induce the B-to-Z conversion. Additionally, we identify Zα-like domains within bacterial proteins. While these domains share certain features with Zα domains, they lack the ability to bind to Z-nucleic acids or facilitate the B-to-Z DNA conversion. Our findings significantly expand the ZBP family across a wide spectrum of organisms and raise intriguing questions about the evolutionary origins of Zα-containing proteins. Moreover, our study offers fresh perspectives on the functional significance of Zα domains in virus sensing and innate immunity and opens avenues for exploring hitherto undiscovered functions of ZBPs.

Unraveling the habitat preferences, ecological drivers, potential hosts, and auxiliary metabolism of soil giant viruses across China

BACKGROUND

Soil giant viruses are increasingly believed to have profound effects on ecological functioning by infecting diverse eukaryotes. However, their biogeography and ecology remain poorly understood.

RESULTS

In this study, we analyzed 333 soil metagenomes from 5 habitat types (farmland, forest, grassland, Gobi desert, and mine wasteland) across China and identified 533 distinct giant virus phylotypes affiliated with nine families, thereby greatly expanding the diversity of soil giant viruses. Among the nine families, Pithoviridae were the most diverse. The majority of phylotypes exhibited a heterogeneous distribution among habitat types, with a remarkably high proportion of unique phylotypes in mine wasteland. The abundances of phylotypes were negatively correlated with their environmental ranges. A total of 76 phylotypes recovered in this study were detectable in a published global topsoil metagenome dataset. Among climatic, geographical, edaphic, and biotic characteristics, soil eukaryotes were identified as the most important driver of beta-diversity of giant viral communities across habitat types. Moreover, co-occurrence network analysis revealed some pairings between giant viral phylotypes and eukaryotes (protozoa, fungi, and algae). Analysis of 44 medium- to high-quality giant virus genomes recovered from our metagenomes uncovered not only their highly shared functions but also their novel auxiliary metabolic genes related to carbon, sulfur, and phosphorus cycling.

CONCLUSIONS

These findings extend our knowledge of diversity, habitat preferences, ecological drivers, potential hosts, and auxiliary metabolism of soil giant viruses. Video Abstract.

DNA methylation enables recurrent endogenization of giant viruses in an animal relative

5-Methylcytosine (5mC) is a widespread silencing mechanism that controls genomic parasites. In eukaryotes, 5mC has gained complex roles in gene regulation beyond parasite control, yet 5mC has also been lost in many lineages. The causes for 5mC retention and its genomic consequences are still poorly understood. Here, we show that the protist closely related to animals Amoebidium appalachense features both transposon and gene body methylation, a pattern reminiscent of invertebrates and plants. Unexpectedly, hypermethylated genomic regions in Amoebidium derive from viral insertions, including hundreds of endogenized giant viruses, contributing 14% of the proteome. Using a combination of inhibitors and genomic assays, we demonstrate that 5mC silences these giant virus insertions. Moreover, alternative Amoebidium isolates show polymorphic giant virus insertions, highlighting a dynamic process of infection, endogenization, and purging. Our results indicate that 5mC is critical for the controlled coexistence of newly acquired viral DNA into eukaryotic genomes, making Amoebidium a unique model to understand the hybrid origins of eukaryotic DNA.

Characterization of CRESS-DNA viruses in human vaginal secretions: An exploratory metagenomic investigation

The Phylum Cressdnaviricota consists of a large number of circular Rep-encoding single-stranded (CRESS)-DNA viruses. Recently, metagenomic analyzes revealed their ubiquitous distribution in a diverse range of eukaryotes. Data relating to CRESS-DNA viruses in humans remains scarce. Our study investigated the presence and genetic diversity of CRESS-DNA viruses in human vaginal secretions. Vaginal swabs were collected from 28 women between 29 and 43 years old attending a fertility clinic in New York City. An exploratory metagenomic analysis was performed and detection of CRESS-DNA viruses was confirmed through analysis of near full-length sequences of the viral isolates. A phylogenetic tree was based on the REP open reading frame sequences of the CRESS-DNA virus genome. Eleven nearly complete CRESS-DNA viral genomes were identified in 16 (57.1%) women. There were no associations between the presence of these viruses and any demographic or clinical parameters. Phylogenetic analysis indicated that one of the sequences belonged to the genus Gemycircularvirus within the Genomoviridae family, while ten sequences represented previously unclassified species of CRESS-DNA viruses. Novel species of CRESS-DNA viruses are present in the vaginal tract of adult women. Although they be transient commensal agents, the potential clinical implications for their presence at this site cannot be dismissed.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

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.

Giant viral signatures on the Greenland ice sheet

BACKGROUND

Dark pigmented snow and glacier ice algae on glaciers and ice sheets contribute to accelerating melt. The biological controls on these algae, particularly the role of viruses, remain poorly understood. Giant viruses, classified under the nucleocytoplasmic large DNA viruses (NCLDV) supergroup (phylum Nucleocytoviricota), are diverse and globally distributed. NCLDVs are known to infect eukaryotic cells in marine and freshwater environments, providing a biological control on the algal population in these ecosystems. However, there is very limited information on the diversity and ecosystem function of NCLDVs in terrestrial icy habitats.

RESULTS

In this study, we investigate for the first time giant viruses and their host connections on ice and snow habitats, such as cryoconite, dark ice, ice core, red and green snow, and genomic assemblies of five cultivated Chlorophyta snow algae. Giant virus marker genes were present in almost all samples; the highest abundances were recovered from red snow and the snow algae genomic assemblies, followed by green snow and dark ice. The variety of active algae and protists in these GrIS habitats containing NCLDV marker genes suggests that infection can occur on a range of eukaryotic hosts. Metagenomic data from red and green snow contained evidence of giant virus metagenome-assembled genomes from the orders Imitervirales, Asfuvirales, and Algavirales.

CONCLUSION

Our study highlights NCLDV family signatures in snow and ice samples from the Greenland ice sheet. Giant virus metagenome-assembled genomes (GVMAGs) were found in red snow samples, and related NCLDV marker genes were identified for the first time in snow algal culture genomic assemblies; implying a relationship between the NCLDVs and snow algae. Metatranscriptomic viral genes also aligned with metagenomic sequences, suggesting that NCLDVs are an active component of the microbial community and are potential "top-down" controls of the eukaryotic algal and protistan members. This study reveals the unprecedented presence of a diverse community of NCLDVs in a variety of glacial habitats dominated by algae.

Continuous year-round isolation of giant viruses from brackish shoreline soils

Giant viruses, categorized under Nucleocytoviricota, are believed to exist ubiquitously in natural environments. However, comprehensive reports on isolated giant viruses remain scarce, with limited information available on unrecoverable strains, viral proliferation sites, and natural hosts. Previously, the author highlighted Pandoravirus hades, Pandoravirus persephone, and Mimivirus sp. styx, isolated from brackish water soil, as potential hotspots for giant virus multiplication. This study presents findings from nearly a year of monthly sampling within the same brackish water region after isolating the three aforementioned strains. This report details the recurrent isolation of a wide range of giant viruses. Each month, four soil samples were randomly collected from an approximately 5 × 10 m plot, comprising three soil samples and one water sample containing sediment from the riverbed. Acanthamoeba castellanii was used as a host for virus isolation. These efforts consistently yielded at least one viral species per month, culminating in a total of 55 giant virus isolates. The most frequently isolated species was Mimiviridae (24 isolates), followed by Marseilleviridae (23 isolates), Pandoravirus (6 isolates), and singular isolates of Pithovirus and Cedratvirus. Notably, viruses were not consistently isolated from any of the four samples every month, with certain sites yielding no viruses. Cluster analysis based on isolate numbers revealed that soil samples from May and water and sediment samples from January produced the highest number of viral strains. These findings underscore brackish coastal soil as a significant site for isolating numerous giant viruses, highlighting the non-uniform distribution along coastlines.

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.

Isolation of viruses, including mollivirus, with the potential to infect Acanthamoeba from a Japanese warm temperate zone

Acanthamoeba castellanii is infected with diverse nucleocytoplasmic large DNA viruses. Here, we report the co-isolation of 12 viral strains from marine sediments in Uranouchi Inlet, Kochi, Japan. Based on the morphological features revealed by electron microscopy, these isolates were classified into four viral groups including Megamimiviridae, Molliviridae, Pandoraviridae, and Pithoviridae. Genomic analyses indicated that these isolates showed high similarities to the known viral genomes with which they are taxonomically clustered, and their phylogenetic relationships were also supported by core gene similarities. It is noteworthy that Molliviridae was isolated from the marine sediments in the Japanese warm temperate zone because other strains have only been found in the subarctic region. Furthermore, this strain has 19 and 4 strain-specific genes found in Mollivirus sibericum and Mollivirus kamchatka, respectively. This study extends our knowledge about the habitat and genomic diversity of Molliviridae.

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 (Taxonomic Information of Giant viruses using Trademark Orthologous Groups), 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 1531 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% at the order level and 97.3% at 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% at 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.

A novel and diverse family of filamentous DNA viruses associated with parasitic wasps

Large dsDNA viruses from the Naldaviricetes class are currently composed of four viral families infecting insects and/or crustaceans. Since the 1970s, particles described as filamentous viruses (FVs) have been observed by electronic microscopy in several species of Hymenoptera parasitoids but until recently, no genomic data was available. This study provides the first comparative morphological and genomic analysis of these FVs. We analyzed the genomes of seven FVs, six of which were newly obtained, to gain a better understanding of their evolutionary history. We show that these FVs share all genomic features of the Naldaviricetes while encoding five specific core genes that distinguish them from their closest relatives, the Hytrosaviruses. By mining public databases, we show that FVs preferentially infect Hymenoptera with parasitoid lifestyle and that these viruses have been repeatedly integrated into the genome of many insects, particularly Hymenoptera parasitoids, overall suggesting a long-standing specialization of these viruses to parasitic wasps. Finally, we propose a taxonomical revision of the class Naldaviricetes in which FVs related to the Leptopilina boulardi FV constitute a fifth family. We propose to name this new family, Filamentoviridae.

COBRA improves the completeness and contiguity of viral genomes assembled from metagenomes

Viruses are often studied using metagenome-assembled sequences, but genome incompleteness hampers comprehensive and accurate analyses. Contig Overlap Based Re-Assembly (COBRA) resolves assembly breakpoints based on the de Bruijn graph and joins contigs. Here we benchmarked COBRA using ocean and soil viral datasets. COBRA accurately joined the assembled sequences and achieved notably higher genome accuracy than binning tools. From 231 published freshwater metagenomes, we obtained 7,334 bacteriophage clusters, ~83% of which represent new phage species. Notably, ~70% of these were circular, compared with 34% before COBRA analyses. We expanded sampling of huge phages (≥200 kbp), the largest of which was curated to completion (717 kbp). Improved phage genomes from Rotsee Lake provided context for metatranscriptomic data and indicated the in situ activity of huge phages, whiB-encoding phages and cysC- and cysH-encoding phages. COBRA improves viral genome assembly contiguity and completeness, thus the accuracy and reliability of analyses of gene content, diversity and evolution.

A broad wastewater screening and clinical data surveillance for virus-related diseases in the metropolitan Detroit area in Michigan

BACKGROUND

Periodic bioinformatics-based screening of wastewater for assessing the diversity of potential human viral pathogens circulating in a given community may help to identify novel or potentially emerging infectious diseases. Any identified contigs related to novel or emerging viruses should be confirmed with targeted wastewater and clinical testing.

RESULTS

During the COVID-19 pandemic, untreated wastewater samples were collected for a 1-year period from the Great Lakes Water Authority Wastewater Treatment Facility in Detroit, MI, USA, and viral population diversity from both centralized interceptor sites and localized neighborhood sewersheds was investigated. Clinical cases of the diseases caused by human viruses were tabulated and compared with data from viral wastewater monitoring. In addition to Betacoronavirus, comparison using assembled contigs against a custom Swiss-Prot human virus database indicated the potential prevalence of other pathogenic virus genera, including: Orthopoxvirus, Rhadinovirus, Parapoxvirus, Varicellovirus, Hepatovirus, Simplexvirus, Bocaparvovirus, Molluscipoxvirus, Parechovirus, Roseolovirus, Lymphocryptovirus, Alphavirus, Spumavirus, Lentivirus, Deltaretrovirus, Enterovirus, Kobuvirus, Gammaretrovirus, Cardiovirus, Erythroparvovirus, Salivirus, Rubivirus, Orthohepevirus, Cytomegalovirus, Norovirus, and Mamastrovirus. Four nearly complete genomes were recovered from the Astrovirus, Enterovirus, Norovirus and Betapolyomavirus genera and viral species were identified.

CONCLUSIONS

The presented findings in wastewater samples are primarily at the genus level and can serve as a preliminary "screening" tool that may serve as indication to initiate further testing for the confirmation of the presence of species that may be associated with human disease. Integrating innovative environmental microbiology technologies like metagenomic sequencing with viral epidemiology offers a significant opportunity to improve the monitoring of, and predictive intelligence for, pathogenic viruses, using wastewater.

Ecological and evolutionary dynamics of cell-virus-virophage systems

Microbial eukaryotes, giant viruses and virophages form a unique hyperparasitic system. Virophages are parasites of the virus transcription machinery and can interfere with virus replication, resulting in a benefit to the eukaryotic host population. Surprisingly, virophages can integrate into the genomes of their cell or virus hosts, and have been shown to reactivate during coinfection. This raises questions about the role of integration in the dynamics of cell-virus-virophage systems. We use mathematical models and computational simulations to understand the effect of virophage integration on populations of cells and viruses. We also investigate multicellularity and programmed cell-death (PCD) as potential antiviral defence strategies used by cells. We found that virophages which enter the cell independently of the host virus, such as Mavirus, are expected to integrate commonly into the genomes of their cell hosts. Our models suggest that integrations from virophages without an independent mode of entry like Sputnik, are less likely to become fixed in the cell host population. Alternatively, we found that Sputnik virophages can stably persist integrated in the virus population, as long as they do not completely inhibit virus replication. We also show that increasing virophage inhibition can stabilise oscillatory dynamics, which may explain the long-term persistence of viruses and virophages in the environment. Our results demonstrate that inhibition by virophages and multicellularity are effective antiviral strategies that may act in synergy against viral infection in microbial species.