The DNA methylation landscape of giant viruses

Refining the taxonomy of pithovirus-related giant DNA viruses within the order Pimascovirales

The first member of the family Pithoviridae (Pithovirus sibericum) was isolated from ancient Siberian permafrost and characterized in 2014. Since then, many relatives have been isolated, characterized, and classified as members of the genera Alphapithovirus, Alphacedratvirus, and Alphaorpheovirus. In addition, one complete circular genome sequence was assembled from metagenomic data (hydrivirus). All of these viruses form distinctive giant elongated ovoid particles, up to 2 µm in length, but they differ significantly in the size of their genome, their nucleotide composition, and their gene content. Based on their shared ovoid virion shape, common replication strategy, and core gene similarity, we recently proposed to update their taxonomic status by classifying them in three distinct families (Pithoviridae, Orpheoviridae, and Hydriviridae) within a new suborder, the Ocovirineae, to separate them clearly from the other more distant families (Marseilleviridae, Ascoviridae, Iridoviridae) of the order Pimascovirales. This new taxonomy, validated by the last ICTV Ratification vote held in March 2025, extends the previous partition from three clades to four (to include hydrivirus) while keeping the genera Alphacedratvirus and Alphapithovirus in the same family, Pithoviridae (but split into two subfamilies), due to their much greater similarity to each other than to orpheovirus and hydrivirus.

Host-encoded DNA methyltransferases modify the epigenome and host tropism of invading phages

Restriction modification (RM) systems are ubiquitous bacterial defense systems; however, some phages evade RM system and adapt to their bacterial hosts. In such cases, phages are thought to stochastically acquire DNA methylation from host-encoded DNA methyltransferases (MTases), facilitating host adaptation. However, no studies have directly compared the methylomes of host bacteria and their infecting phages. Here, we demonstrate the epigenetic landscape of adapted phages with diverse infection histories, focusing on the broad host-range phage KHP30T as its adapts to three Helicobacter pylori strains. Using a multistage infection system, we observed that the adapted phages displayed significantly high titers against the last infected H. pylori strain, suggesting an attendant change in host tropism. Single-molecule real-time sequencing revealed that methylated motifs were predominantly shared between the adapted phages and their most recent host. Our findings enhance our understanding of epigenetic phage-host interactions, which have significant implications for microbial ecology.

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.

Unravelling viral ecology and evolution over 20 years in a freshwater lake

As freshwater lakes undergo rapid anthropogenic change, long-term studies reveal key microbial dynamics, evolutionary shifts and biogeochemical interactions, yet the vital role of viruses remains overlooked. Here, leveraging a 20 year time series from Lake Mendota, WI, USA, we characterized 1.3 million viral genomes across time, seasonality and environmental factors. Double-stranded DNA phages from the class Caudoviricetes dominated the community. We identified 574 auxiliary metabolic gene families representing over 140,000 auxiliary metabolic genes, including important genes such as psbA (photosynthesis), pmoC (methane oxidation) and katG (hydrogen peroxide decomposition), which were consistently present and active across decades and seasons. Positive associations and niche differentiation between virus-host pairs, including keystone Cyanobacteria, methanotrophs and Nanopelagicales, emerged during seasonal changes. Inorganic carbon and ammonium influenced viral abundances, underscoring viral roles in both 'top-down' and 'bottom-up' interactions. Evolutionary processes favoured fitness genes, reduced genomic heterogeneity and dominant sub-populations. This study transforms understanding of viral ecology and evolution in Earth's microbiomes.

Metagenomic characterization of viruses and mobile genetic elements associated with the DPANN archaeal superphylum

The archaeal superphylum DPANN (an acronym formed from the initials of the first five phyla discovered: Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanohaloarchaeota and Nanoarchaeota) is a group of ultrasmall symbionts able to survive in extreme ecosystems. The diversity and dynamics between DPANN archaea and their virome remain largely unknown. Here we use a metagenomic clustered regularly interspaced short palindromic repeats (CRISPR) screening approach to identify 97 globally distributed, non-redundant viruses and unclassified mobile genetic elements predicted to infect hosts across 8 DPANN phyla, including 7 viral groups not previously characterized. Genomic analysis suggests a diversity of viral morphologies including head-tailed, tailless icosahedral and spindle-shaped viruses with the potential to establish lytic, chronic or lysogenic infections. We also find evidence of a virally encoded Cas12f1 protein (probably originating from uncultured DPANN archaea) and a mini-CRISPR array, which could play a role in modulating host metabolism. Many metagenomes have virus-to-host ratios >10, indicating that DPANN viruses play an important role in controlling host populations. Overall, our study illuminates the underexplored diversity, functional repertoires and host interactions of the DPANN virome.

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.

Long-read sequencing reveals extensive gut phageome structural variations driven by genetic exchange with bacterial hosts

Genetic variations are instrumental for unraveling phage evolution and deciphering their functional implications. Here, we explore the underlying fine-scale genetic variations in the gut phageome, especially structural variations (SVs). By using virome-enriched long-read metagenomic sequencing across 91 individuals, we identified a total of 14,438 nonredundant phage SVs and revealed their prevalence within the human gut phageome. These SVs are mainly enriched in genes involved in recombination, DNA methylation, and antibiotic resistance. Notably, a substantial fraction of phage SV sequences share close homology with bacterial fragments, with most SVs enriched for horizontal gene transfer (HGT) mechanism. Further investigations showed that these SV sequences were genetic exchanged between specific phage-bacteria pairs, particularly between phages and their respective bacterial hosts. Temperate phages exhibit a higher frequency of genetic exchange with bacterial chromosomes and then virulent phages. Collectively, our findings provide insights into the genetic landscape of the human gut phageome.

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.

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.

Are Viruses Taxonomic Units? A Protein Domain and Loop-Centric Phylogenomic Assessment

Virus taxonomy uses a Linnaean-like subsumption hierarchy to classify viruses into taxonomic units at species and higher rank levels. Virus species are considered monophyletic groups of mobile genetic elements (MGEs) often delimited by the phylogenetic analysis of aligned genomic or metagenomic sequences. Taxonomic units are assumed to be independent organizational, functional and evolutionary units that follow a 'natural history' rationale. Here, I use phylogenomic and other arguments to show that viruses are not self-standing genetically-driven systems acting as evolutionary units. Instead, they are crucial components of holobionts, which are units of biological organization that dynamically integrate the genetics, epigenetic, physiological and functional properties of their co-evolving members. Remarkably, phylogenomic analyses show that viruses share protein domains and loops with cells throughout history via massive processes of reticulate evolution, helping spread evolutionary innovations across a wider taxonomic spectrum. Thus, viruses are not merely MGEs or microbes. Instead, their genomes and proteomes conduct cellularly integrated processes akin to those cataloged by the GO Consortium. This prompts the generation of compositional hierarchies that replace the 'is-a-kind-of' by a 'is-a-part-of' logic to better describe the mereology of integrated cellular and viral makeup. My analysis demands a new paradigm that integrates virus taxonomy into a modern evolutionarily centered taxonomy of organisms.

A taxonomic proposal for cedratviruses, orpheoviruses, and pithoviruses

Orpheoviruses, cedratviruses, and pithoviruses are large DNA viruses that cluster together taxonomically within the order Pimascovirales of the phylum Nucleocytoviricota. However, they were not classified previously by the International Committee on Taxonomy of Viruses (ICTV). Here, we present a comprehensive analysis of the gene content, morphology, and phylogenomics of these viruses, providing data that underpinned the recent proposal to establish new taxa for their initial classification. The new taxonomy, which has now been ratified by the ICTV, includes the family Orpheoviridae and genus Alphaorpheovirus, the family Pithoviridae and genus Alphapithovirus, and the family Cedratviridae and genus Alphacedratvirus, aiming to formally catalogue the isolates covered in this study. Additionally, as per the newly adopted rules, we applied standardized binomial names for the virus species created to classify isolates with complete genome sequences available in public databases at the time of the proposal. The specific epithet of each virus species was chosen as a reference to the location where the exemplar virus was isolated.

Dinucleotide biases in the genomes of prokaryotic and eukaryotic dsDNA viruses and their hosts

The genomes of cellular organisms display CpG and TpA dinucleotide composition biases. Such biases have been poorly investigated in dsDNA viruses. Here, we show that in dsDNA virus, bacterial, and eukaryotic genomes, the representation of TpA and CpG dinucleotides is strongly dependent on genomic G + C content. Thus, the classical observed/expected ratios do not fully capture dinucleotide biases across genomes. Because a larger portion of the variance in TpA frequency was explained by G + C content, we explored which additional factors drive the distribution of CpG dinucleotides. Using the residuals of the linear regressions as a measure of dinucleotide abundance and ancestral state reconstruction across eukaryotic and prokaryotic virus trees, we identified an important role for phylogeny in driving CpG representation. Nonetheless, phylogenetic ANOVA analyses showed that few host associations also account for significant variations. Among eukaryotic viruses, most significant differences were observed between arthropod-infecting viruses and viruses that infect vertebrates or unicellular organisms. However, an effect of viral DNA methylation status (either driven by the host or by viral-encoded methyltransferases) is also likely. Among prokaryotic viruses, cyanobacteria-infecting phages resulted to be significantly CpG-depleted, whereas phages that infect bacteria in the genera Burkolderia and Staphylococcus were CpG-rich. Comparison with bacterial genomes indicated that this effect is largely driven by the general tendency for phages to resemble the host's genomic CpG content. Notably, such tendency is stronger for temperate than for lytic phages. Our data shed light into the processes that shape virus genome composition and inform manipulation strategies for biotechnological applications.

Isolation of a widespread giant virus implicated in cryptophyte bloom collapse

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

Temperature modulates dominance of a superinfecting Arctic virus in its unicellular algal host

Complex virus-virus interactions can arise when multiple viruses coinfect the same host, impacting infection outcomes with broader ecological and evolutionary significance for viruses and host. Yet, our knowledge regarding virus competition is still limited, especially for single-celled eukaryotic host-virus systems. Here, we report on mutual interference of two dsDNA viruses, MpoV-45T and MpoV-46T, competing for their Arctic algal host Micromonas polaris. Both viruses affected each other's gene expression and displayed reduced genome replication during coinfection. MpoV-45T was the dominant virus, likely due to interference in the DNA replication of is competitor. Even when its coinfection was delayed, the dominant virus still prevailed while genome production of the other virus was strongly suppressed. This contrasts with typical superinfection exclusion, where the primary infection prevents secondary infection by other viruses. Higher temperature made the suppressed virus a stronger competitor, signifying that global warming is likely to alter virus-virus interactions in Arctic waters.

Unveiling Prasinovirus diversity and host specificity through targeted enrichment in the South China Sea

Unicellular green picophytoplankton from the Mamiellales order are pervasive in marine ecosystems and susceptible to infections by prasinoviruses, large double-stranded DNA viruses within the Nucleocytoviricota phylum. We developed a double-stranded DNA virus enrichment and shotgun sequencing method, and successfully assembled 80 prasinovirus genomes from 43 samples in the South China Sea. Our research delivered the first direct estimation of 94% accuracy in correlating genome similarity to host range. Stirkingly, our analyses uncovered unexpected host-switching across diverse algal lineages, challenging the existing paradigms of host-virus co-speciation and revealing the dynamic nature of viral evolution. We also detected six instances of horizontal gene transfer between prasinoviruses and their hosts, including a novel alternative oxidase. Additionally, diversifying selection on a major capsid protein suggests an ongoing co-evolutionary arms race. These insights not only expand our understanding of prasinovirus genomic diversity but also highlight the intricate evolutionary mechanisms driving their ecological success and shaping broader virus-host interactions in marine environments.

A sequential one-pot approach for rapid and convenient characterization of putative restriction-modification systems

IMPORTANCE

The elucidation of the molecular basis of virus-host coevolutionary interactions is boosted with state-of-the-art sequencing technologies. However, the sequence-only information is often insufficient to output a conclusive argument without biochemical characterizations. We proposed a 1-day and one-pot approach to confirm the exact function of putative restriction-modification (R-M) genes that presumably mediate microbial coevolution. The experiments mainly focused on a series of putative R-M enzymes from a deep-sea virus and its host bacterium. The results quickly unveiled unambiguous substrate specificities, superior catalytic performance, and unique sequence preferences for two new restriction enzymes (capable of cleaving DNA) and two new methyltransferases (capable of modifying DNA with methyl groups). The reality of the functional R-M system reinforced a model of mutually beneficial interactions with the virus in the deep-sea microbial ecosystem. The cell culture-independent approach also holds great potential for exploring novel and biotechnologically significant R-M enzymes from microbial dark matter.

PacBio sequencing of human fecal samples uncovers the DNA methylation landscape of 22 673 gut phages

Gut phages have an important impact on human health. Methylation plays key roles in DNA recognition, gene expression regulation and replication for phages. However, the DNA methylation landscape of gut phages is largely unknown. Here, with PacBio sequencing (2120×, 4785 Gb), we detected gut phage methylation landscape based on 22 673 gut phage genomes, and presented diverse methylation motifs and methylation differences in genomic elements. Moreover, the methylation rate of phages was associated with taxonomy and host, and N6-methyladenine methylation rate was higher in temperate phages than in virulent phages, suggesting an important role for methylation in phage-host interaction. In particular, 3543 (15.63%) phage genomes contained restriction-modification system, which could aid in evading clearance by the host. This study revealed the DNA methylation landscape of gut phage and its potential roles, which will advance the understanding of gut phage survival and human health.

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

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

Pithoviruses Are Invaded by Repeats That Contribute to Their Evolution and Divergence from Cedratviruses

Pithoviridae are amoeba-infecting giant viruses possessing the largest viral particles known so far. Since the discovery of Pithovirus sibericum, recovered from a 30,000-yr-old permafrost sample, other pithoviruses, and related cedratviruses, were isolated from various terrestrial and aquatic samples. Here, we report the isolation and genome sequencing of 2 Pithoviridae from soil samples, in addition to 3 other recent isolates. Using the 12 available genome sequences, we conducted a thorough comparative genomic study of the Pithoviridae family to decipher the organization and evolution of their genomes. Our study reveals a nonuniform genome organization in 2 main regions: 1 concentrating core genes and another gene duplications. We also found that Pithoviridae genomes are more conservative than other families of giant viruses, with a low and stable proportion (5% to 7%) of genes originating from horizontal transfers. Genome size variation within the family is mainly due to variations in gene duplication rates (from 14% to 28%) and massive invasion by inverted repeats. While these repeated elements are absent from cedratviruses, repeat-rich regions cover as much as a quarter of the pithoviruses genomes. These regions, identified using a dedicated pipeline, are hotspots of mutations, gene capture events, and genomic rearrangements that contribute to their evolution.

A rapid genome-wide analysis of isolated giant viruses using MinION sequencing

Following the discovery of Acanthamoeba polyphaga mimivirus, diverse giant viruses have been isolated. However, only a small fraction of these isolates have been completely sequenced, limiting our understanding of the genomic diversity of giant viruses. MinION is a portable and low-cost long-read sequencer that can be readily used in a laboratory. Although MinION provides highly error-prone reads that require correction through additional short-read sequencing, recent studies assembled high-quality microbial genomes only using MinION sequencing. Here, we evaluated the accuracy of MinION-only genome assemblies for giant viruses by re-sequencing a prototype marseillevirus. Assembled genomes presented over 99.98% identity to the reference genome with a few gaps, demonstrating a high accuracy of the MinION-only assembly. As a proof of concept, we de novo assembled five newly isolated viruses. Average nucleotide identities to their closest known relatives suggest that the isolates represent new species of marseillevirus, pithovirus and mimivirus. The assembly of subsampled reads demonstrated that their taxonomy and genomic composition could be analysed at the 50× sequencing coverage. We also identified a pithovirus gene whose homologues were detected only in metagenome-derived relatives. Collectively, we propose that MinION-only assembly is an effective approach to rapidly perform a genome-wide analysis of isolated giant viruses.

Single-cell transcription mapping of murine and human mammary organoids responses to female hormones

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Long-Read Sequencing Reveals Extensive DNA Methylations in Human Gut Phagenome Contributed by Prevalently Phage-Encoded Methyltransferases

DNA methylation plays a crucial role in the survival of bacteriophages (phages), yet the understanding of their genome methylation remains limited. In this study, DNA methylation patterns are analyzed in 8848 metagenome-assembled high-quality phages from 104 fecal samples using single-molecule real-time sequencing. The results demonstrate that 97.60% of gut phages exhibit methylation, with certain factors correlating with methylation densities. Phages with higher methylation densities appear to have potential viability advantages. Strikingly, more than one-third of the phages possess their own DNA methyltransferases (MTases). Increased MTase copies are associated with higher genome methylation densities, specific methylation motifs, and elevated prevalence of certain phage groups. Notably, the majority of these MTases share close homology with those encoded by gut bacteria, suggesting their exchange during phage-bacterium interactions. Furthermore, these MTases can be employed to accurately predict phage-host relationships. Overall, the findings indicate the widespread utilization of DNA methylation by gut DNA phages as an evasion mechanism against host defense systems, with a substantial contribution from phage-encoded MTases.

Analysis of the Genomic Features and Evolutionary History of Pithovirus-Like Isolates Reveals Two Major Divergent Groups of Viruses

New representatives of the phylum Nucleocytoviricota have been rapidly described in the last decade. Despite this, not all viruses of this phylum are allocated to recognized taxonomic families, as is the case for orpheovirus, pithovirus, and cedratvirus, which form the proposed family Pithoviridae. In this study, we performed comprehensive comparative genomic analyses of 8 pithovirus-like isolates, aiming to understand their common traits and evolutionary history. Structural and functional genome annotation was performed de novo for all the viruses, which served as a reference for pangenome construction. The synteny analysis showed substantial differences in genome organization between these viruses, with very few and short syntenic blocks shared between orpheovirus and its relatives. It was possible to observe an open pangenome with a significant increase in the slope when orpheovirus was added, alongside a decrease in the core genome. Network analysis placed orpheovirus as a distant and major hub with a large fraction of unique clusters of orthologs, indicating a distant relationship between this virus and its relatives, with only a few shared genes. Additionally, phylogenetic analyses of strict core genes shared with other viruses of the phylum reinforced the divergence of orpheovirus from pithoviruses and cedratviruses. Altogether, our results indicate that although pithovirus-like isolates share common features, this group of ovoid-shaped giant viruses presents substantial differences in gene contents, genomic architectures, and the phylogenetic history of several core genes. Our data indicate that orpheovirus is an evolutionarily divergent viral entity, suggesting its allocation to a different viral family, Orpheoviridae. IMPORTANCE Giant viruses that infect amoebae form a monophyletic group named the phylum Nucleocytoviricota. Despite being genomically and morphologically very diverse, the taxonomic categories of some clades that form this phylum are not yet well established. With advances in isolation techniques, the speed at which new giant viruses are described has increased, escalating the need to establish criteria to define the emerging viral taxa. In this work, we performed a comparative genomic analysis of representatives of the putative family Pithoviridae. Based on the dissimilarity of orpheovirus from the other viruses of this putative family, we propose that orpheovirus be considered a member of an independent family, Orpheoviridae, and suggest criteria to demarcate families consisting of ovoid-shaped giant viruses.

A billion years arms-race between viruses, virophages, and eukaryotes

Bamfordviruses are arguably the most diverse group of viruses infecting eukaryotes. They include the Nucleocytoplasmic Large DNA viruses (NCLDVs), virophages, adenoviruses, Mavericks and Polinton-like viruses. Two main hypotheses for their origins have been proposed: the 'nuclear-escape' and 'virophage-first' hypotheses. The nuclear-escape hypothesis proposes an endogenous, Maverick-like ancestor which escaped from the nucleus and gave rise to adenoviruses and NCLDVs. In contrast, the virophage-first hypothesis proposes that NCLDVs coevolved with protovirophages; Mavericks then evolved from virophages that became endogenous, with adenoviruses escaping from the nucleus at a later stage. Here, we test the predictions made by both models and consider alternative evolutionary scenarios. We use a data set of the four core virion proteins sampled across the diversity of the lineage, together with Bayesian and maximum-likelihood hypothesis-testing methods, and estimate rooted phylogenies. We find strong evidence that adenoviruses and NCLDVs are not sister groups, and that Mavericks and Mavirus acquired the rve-integrase independently. We also found strong support for a monophyletic group of virophages (family Lavidaviridae) and a most likely root placed between virophages and the other lineages. Our observations support alternatives to the nuclear-escape scenario and a billion years evolutionary arms-race between virophages and NCLDVs.

Diagnosis and monitoring of virus-associated cancer using cell-free DNA

Viral-associated cancers are a distinct group of malignancies with a unique pathogenesis and epidemiology. Liquid biopsy is a minimally invasive way to identify tumor-associated abnormalities in blood derivatives, such as plasma, to guide the diagnosis, prognosis, and treatment of patients with cancer. Liquid biopsy encompasses a multitude of circulating analytes with the most extensively studied being cell-free DNA (cfDNA). In recent decades, substantial advances have been made toward the study of circulating tumor DNA in nonviral-associated cancers. Many of these observations have been translated to the clinic to improve the outcomes of patients with cancer. The study of cfDNA in viral-associated cancers is rapidly evolving and reveals tremendous potential for clinical applications. This review provides an overview of the pathogenesis of viral-associated malignancies, the current state of cfDNA analysis in oncology, the current state of cfDNA analysis in viral-associated cancers, and perspectives for the future of liquid biopsies in viral-associated cancers.

A Unique m6A-Dependent Restriction Endonuclease from an Archaeal Virus

Prokaryotes possess numerous diverse defense systems to resist viral infections, while some viruses have also evolved antiviral defense systems to exclude other viruses in cases of multiple infections. Here, we report the first virus-derived modification-dependent restriction endonuclease (HHPV4I) from the archaeal virus HHPV4 (Haloarcula hispanica pleomorphic virus 4). HHPV4I contains an SRA domain, a winged helix (wH) domain, and an HNH domain; recognizes the Gm6ATC site; and specifically binds to Gm6ATC site-containing DNA. Both the wH domain and the HNH domain are responsible for DNA binding. Unlike the well-known m6A-specific restriction enzyme DpnI, HHPV4I only efficiently cleaves DNA with a fully methylated Gm6ATC site and cleaves DNA both upstream and downstream of the Gm6ATC sites on both DNA strands. Furthermore, HHPV4I preferentially cleaves DNA between VR bases (V = A/G/C, R = A/G) 4 to 20 nt away from the Gm6ATC site. Thus, the cleavage pattern of HHPV4I is distinct from those of all of the presently characterized restriction endonucleases. Mutations in the wH domain of HHPV4I do not alter m6A-dependent endonuclease activity, but they decrease recognition sequence specificity, thus expanding the cleaving capacity to more m6A-containing DNA sequences. The wH domain provides a target for searching, developing, and engineering novel m6A-dependent endonucleases. IMPORTANCE Many modification-dependent restriction endonucleases (MDREs) were identified in prokaryotes and recognized modified cytosine bases, such as 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC), and glucosyl-5-hydroxymethylcytosine (g5hmC). The first virus-derived MDRE (HHPV4I) from the archaeal virus HHPV4 was identified in this study. The viral MDRE suggested a new strategy employed by the virus to exclude other viruses in the case of multiple replications. HHPV4I is a novel N6-methyladenine (m6A)-dependent restriction endonuclease, while the cleavage pattern of HHPV4I is distinct from the well-known m6A-dependent restriction endonuclease DpnI. HHPV4I recognizes Gm6ATC sites and cleaves DNA both upstream and downstream of the Gm6ATC sites on both DNA strands. It preferentially cleaves DNA between VR bases (V = A/G/C, R = A/G) 4 to 20 nt away from the Gm6ATC sites. Furthermore, mutations in the HHPV4I wH domain can alter the sequence specificity without impeding the m6A-dependent DNA cleavage activity, providing a target for engineering more m6A-dependent endonucleases with different sequence specificities.

Restriction Endonuclease-Based Modification-Dependent Enrichment (REMoDE) of DNA for Metagenomic Sequencing

Metagenomic sequencing is a swift and powerful tool to ascertain the presence of an organism of interest in a sample. However, sequencing coverage of the organism of interest can be insufficient due to an inundation of reads from irrelevant organisms in the sample. Here, we report a nuclease-based approach to rapidly enrich for DNA from certain organisms, including enterobacteria, based on their differential endogenous modification patterns. We exploit the ability of taxon-specific methylated motifs to resist the action of cognate methylation-sensitive restriction endonucleases that thereby digest unwanted, unmethylated DNA. Subsequently, we use a distributive exonuclease or electrophoretic separation to deplete or exclude the digested fragments, thus enriching for undigested DNA from the organism of interest. As a proof of concept, we apply this method to enrich for the enterobacteria Escherichia coli and Salmonella enterica by 11- to 142-fold from mock metagenomic samples and validate this approach as a versatile means to enrich for genomes of interest in metagenomic samples. IMPORTANCE Pathogens that contaminate the food supply or spread through other means can cause outbreaks that bring devastating repercussions to the health of a populace. Investigations to trace the source of these outbreaks are initiated rapidly but can be drawn out due to the labored methods of pathogen isolation. Metagenomic sequencing can alleviate this hurdle but is often insufficiently sensitive. The approach and implementations detailed here provide a rapid means to enrich for many pathogens involved in foodborne outbreaks, thereby improving the utility of metagenomic sequencing as a tool in outbreak investigations. Additionally, this approach provides a means to broadly enrich for otherwise minute levels of modified DNA, which may escape unnoticed in metagenomic samples.

Virus diversity and interactions with hosts in deep-sea hydrothermal vents

BACKGROUND

The deep sea harbors many viruses, yet their diversity and interactions with hosts in hydrothermal ecosystems are largely unknown. Here, we analyzed the viral composition, distribution, host preference, and metabolic potential in different habitats of global hydrothermal vents, including vent plumes, background seawater, diffuse fluids, and sediments.

RESULTS

From 34 samples collected at eight vent sites, a total of 4662 viral populations (vOTUs) were recovered from the metagenome assemblies, encompassing diverse phylogenetic groups and defining many novel lineages. Apart from the abundant unclassified viruses, tailed phages are most predominant across the global hydrothermal vents, while single-stranded DNA viruses, including Microviridae and small eukaryotic viruses, also constitute a significant part of the viromes. As revealed by protein-sharing network analysis, hydrothermal vent viruses formed many novel genus-level viral clusters and are highly endemic to specific vent sites and habitat types. Only 11% of the vOTUs can be linked to hosts, which are the key microbial taxa of hydrothermal habitats, such as Gammaproteobacteria and Campylobacterota. Intriguingly, vent viromes share some common metabolic features in that they encode auxiliary genes that are extensively involved in the metabolism of carbohydrates, amino acids, cofactors, and vitamins. Specifically, in plume viruses, various auxiliary genes related to methane, nitrogen, and sulfur metabolism were observed, indicating their contribution to host energy conservation. Moreover, the prevalence of sulfur-relay pathway genes indicated the significant role of vent viruses in stabilizing the tRNA structure, which promotes host adaptation to steep environmental gradients.

CONCLUSIONS

The deep-sea hydrothermal systems hold untapped viral diversity with novelty. They may affect both vent prokaryotic and eukaryotic communities and modulate host metabolism related to vent adaptability. More explorations are needed to depict global vent virus diversity and its roles in this unique ecosystem. Video Abstract.

Giant Viruses as a Source of Novel Enzymes for Biotechnological Application

The global demand for industrial enzymes has been increasing in recent years, and the search for new sources of these biological products is intense, especially in microorganisms. Most known viruses have limited genetic machinery and, thus, have been overlooked by the enzyme industry for years. However, a peculiar group of viruses breaks this paradigm. Giant viruses of the phylum Nucleocytoviricota infect protists (i.e., algae and amoebae) and have complex genomes, reaching up to 2.7 Mb in length and encoding hundreds of genes. Different giant viruses have robust metabolic machinery, especially those in the Phycodnaviridae and Mimiviridae families. In this review, we present some peculiarities of giant viruses that infect protists and discuss why they should be seen as an outstanding source of new enzymes. We revisited the genomes of representatives of different groups of giant viruses and put together information about their enzymatic machinery, highlighting several genes to be explored in biotechnology involved in carbohydrate metabolism, DNA replication, and RNA processing, among others. Finally, we present additional evidence based on structural biology using chitinase as a model to reinforce the role of giant viruses as a source of novel enzymes for biotechnological application.

MSNet-4mC: learning effective multi-scale representations for identifying DNA N4-methylcytosine sites

MOTIVATION

N4-methylcytosine (4mC) is an essential kind of epigenetic modification that regulates a wide range of biological processes. However, experimental methods for detecting 4mC sites are time-consuming and labor-intensive. As an alternative, computational methods that are capable of automatically identifying 4mC with data analysis techniques become a reasonable option. A major challenge is how to develop effective methods to fully exploit the complex interactions within the DNA sequences to improve the predictive capability.

RESULTS

In this work, we propose MSNet-4mC, a lightweight neural network building upon convolutional operations with multi-scale receptive fields to perceive cross-element relationships over both short and long ranges of given DNA sequences. With strong imbalances in the number of candidates in different species in mind, we compute and apply class weights in the cross-entropy loss to balance the training process. Extensive benchmarking experiments show that our method achieves a significant performance improvement and outperforms other state-of-the-art methods.

AVAILABILITY AND IMPLEMENTATION

The source code and models are freely available for download at https://github.com/LIU-CT/MSNet-4mC, implemented in Python and supported on Linux and Windows.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

SARS-CoV-2-induced hypomethylation of the ferritin heavy chain (FTH1) gene underlies serum hyperferritinemia in severe COVID-19 patients

High serum ferritin (hyperferritinemia), a reliable hallmark of severe COVID-19 often associates with a moderate decrease in serum iron (hypoferremia) and a moderate increase in serum hepcidin. This suggests that hyperferritinemia in severe COVID-19 is reflective of inflammation rather than iron overload. To test this possibility, the expression status of ferritin heavy chain (FTH1), transferrin receptor 1 (TFRC), hepcidin (HAMP), and ferroportin (SLC40A1) genes and promoter methylation status of FTH1 and TFRC genes were examined in blood samples obtained from COVID-19 patients showing no, mild or severe symptoms and in healthy-donor monocytes stimulated with SARS-CoV-2-derived peptides. Severe COVID-19 samples showed a significant increase in FTH1 expression and hypomethylation relative to mild or asymptomatic COVID-19 samples. S-peptide treated monocytes also showed a significant increase in FTH1 expression and hypomethylation relative to that in controls; treatment with ECD or NP did not change FTH1 expression nor its methylation status. In silico and in vitro analysis showed a significant increase in the expression of the TET3 demethylase in S peptide-treated monocytes. Findings presented here suggest that S peptide-driven hypomethylation of the FTH1 gene promoter underlies hyperferritinemia in severe COVID-19 disease.

Multiple waves of viral invasions in Symbiodiniaceae algal genomes

Dinoflagellates from the family Symbiodiniaceae are phototrophic marine protists that engage in symbiosis with diverse hosts. Their large and distinct genomes are characterized by pervasive gene duplication and large-scale retroposition events. However, little is known about the role and scale of horizontal gene transfer (HGT) in the evolution of this algal family. In other dinoflagellates, high levels of HGTs have been observed, linked to major genomic transitions, such as the appearance of a viral-acquired nucleoprotein that originated via HGT from a large DNA algal virus. Previous work showed that Symbiodiniaceae from different hosts are actively infected by viral groups, such as giant DNA viruses and ssRNA viruses, that may play an important role in coral health. Latent viral infections may also occur, whereby viruses could persist in the cytoplasm or integrate into the host genome as a provirus. This hypothesis received experimental support; however, the cellular localization of putative latent viruses and their taxonomic affiliation are still unknown. In addition, despite the finding of viral sequences in some genomes of Symbiodiniaceae, viral origin, taxonomic breadth, and metabolic potential have not been explored. To address these questions, we searched for putative viral-derived proteins in thirteen Symbiodiniaceae genomes. We found fifty-nine candidate viral-derived HGTs that gave rise to twelve phylogenies across ten genomes. We also describe the taxonomic affiliation of these virus-related sequences, their structure, and their genomic context. These results lead us to propose a model to explain the origin and fate of Symbiodiniaceae viral acquisitions.

Past and present giant viruses diversity explored through permafrost metagenomics

Giant viruses are abundant in aquatic environments and ecologically important through the metabolic reprogramming of their hosts. Less is known about giant viruses from soil even though two of them, belonging to two different viral families, were reactivated from 30,000-y-old permafrost samples. This suggests an untapped diversity of Nucleocytoviricota in this environment. Through permafrost metagenomics we reveal a unique diversity pattern and a high heterogeneity in the abundance of giant viruses, representing up to 12% of the sum of sequence coverage in one sample. Pithoviridae and Orpheoviridae-like viruses were the most important contributors. A complete 1.6 Mb Pithoviridae-like circular genome was also assembled from a 42,000-y-old sample. The annotation of the permafrost viral sequences revealed a patchwork of predicted functions amidst a larger reservoir of genes of unknown functions. Finally, the phylogenetic reconstructions not only revealed gene transfers between cells and viruses, but also between viruses from different families.

Diverse DNA modification in marine prokaryotic and viral communities

DNA chemical modifications, including methylation, are widespread and play important roles in prokaryotes and viruses. However, current knowledge of these modification systems is severely biased towards a limited number of culturable prokaryotes, despite the fact that a vast majority of microorganisms have not yet been cultured. Here, using single-molecule real-time sequencing, we conducted culture-independent 'metaepigenomic' analyses (an integrated analysis of metagenomics and epigenomics) of marine microbial communities. A total of 233 and 163 metagenomic-assembled genomes (MAGs) were constructed from diverse prokaryotes and viruses, respectively, and 220 modified motifs and 276 DNA methyltransferases (MTases) were identified. Most of the MTase genes were not genetically linked with the endonuclease genes predicted to be involved in defense mechanisms against extracellular DNA. The MTase-motif correspondence found in the MAGs revealed 10 novel pairs, 5 of which showed novel specificities and experimentally confirmed the catalytic specificities of the MTases. We revealed novel alternative specificities in MTases that are highly conserved in Alphaproteobacteria, which may enhance our understanding of the co-evolutionary history of the methylation systems and the genomes. Our findings highlight diverse unexplored DNA modifications that potentially affect the ecology and evolution of prokaryotes and viruses in nature.

A Brief History of Giant Viruses’ Studies in Brazilian Biomes

Almost two decades after the isolation of the first amoebal giant viruses, indubitably the discovery of these entities has deeply affected the current scientific knowledge on the virosphere. Much has been uncovered since then: viruses can now acknowledge complex genomes and huge particle sizes, integrating remarkable evolutionary relationships that date as early as the emergence of life on the planet. This year, a decade has passed since the first studies on giant viruses in the Brazilian territory, and since then biomes of rare beauty and biodiversity (Amazon, Atlantic forest, Pantanal wetlands, Cerrado savannas) have been explored in the search for giant viruses. From those unique biomes, novel viral entities were found, revealing never before seen genomes and virion structures. To celebrate this, here we bring together the context, inspirations, and the major contributions of independent Brazilian research groups to summarize the accumulated knowledge about the diversity and the exceptionality of some of the giant viruses found in Brazil.

OUP accepted manuscript

Potyviridae, the largest family of known RNA viruses (realm Riboviria), belongs to the picorna-like supergroup and has important agricultural and ecological impacts. Potyvirid genomes are translated into polyproteins, which are in turn hydrolyzed to release mature products. Recent sequencing efforts revealed an unprecedented number of potyvirids with a rich variability in gene content and genomic layouts. Here, we review the heterogeneity of non-core modules that expand the structural and functional diversity of the potyvirid proteomes. We provide a family-wide classification of P1 proteinases into the functional Types A and B, and discuss pretty interesting sweet potato potyviral ORF (PISPO), putative zinc fingers, and alkylation B (AlkB)—non-core modules found within P1 cistrons. The atypical inosine triphosphate pyrophosphatase (ITPase/HAM1), as well as the pseudo tobacco mosaic virus-like coat protein (TMV-like CP) are discussed alongside homologs of unrelated virus taxa. Family-wide abundance of the multitasking helper component proteinase (HC-pro) is revised. Functional connections between non-core modules are highlighted to support host niche adaptation and immune evasion as main drivers of the Potyviridae evolutionary radiation. Potential biotechnological and synthetic biology applications of potyvirid leader proteinases and non-core modules are finally explored.

Clandestinovirus: A Giant Virus With Chromatin Proteins and a Potential to Manipulate the Cell Cycle of Its Host Vermamoeba vermiformis

For several decades, the vast world of DNA viruses has been expanding constantly. Various discoveries in this field have broadened our knowledge and revealed that DNA viruses encode many functional features, which were once thought to be exclusive to cellular life. Here, we report the isolation of a giant virus named "clandestinovirus," grown on the amoebal host Vermamoeba vermiformis. This virus was discovered in a mixed co-culture associated with another giant virus, Faustovirus ST1. Clandestinovirus possesses a linear dsDNA genome of 581,987 base pairs containing 617 genes. Phylogenetically, clandestinovirus is most closely related to Acanthamoeba castellanii medusavirus and was considered a member of the proposed Medusaviridae family. However, clandestinovirus genome is 65% larger than that of medusavirus, emphasizing the considerable genome size variation within this virus family. Functional annotation of the clandestinovirus genes suggests that the virus encodes four core histones. Furthermore, clandestinovirus appears to orchestrate the cell cycle and mitochondrial activities of the infected host by virtue of encoding a panel of protein kinases and phosphatases, and a suite of functionally diverse mitochondrial protein homologs, respectively. Collectively, these observations illuminate a strategy employed by clandestinovirus to optimize the intracellular environment for efficient virus propagation.

Marseilleviruses: An Update in 2021

The family Marseilleviridae was the second family of giant viruses that was described in 2013, after the family Mimiviridae. Marseillevirus marseillevirus, isolated in 2007 by coculture on Acanthamoeba polyphaga, is the prototype member of this family. Afterward, the worldwide distribution of marseilleviruses was revealed through their isolation from samples of various types and sources. Thus, 62 were isolated from environmental water, one from soil, one from a dipteran, one from mussels, and two from asymptomatic humans, which led to the description of 67 marseillevirus isolates, including 21 by the IHU Méditerranée Infection in France. Recently, five marseillevirus genomes were assembled from deep sea sediment in Norway. Isolated marseilleviruses have ≈250 nm long icosahedral capsids and 348–404 kilobase long mosaic genomes that encode 386–545 predicted proteins. Comparative genomic analyses indicate that the family Marseilleviridae includes five lineages and possesses a pangenome composed of 3,082 clusters of genes. The detection of marseilleviruses in both symptomatic and asymptomatic humans in stool, blood, and lymph nodes, and an up-to-30-day persistence of marseillevirus in rats and mice, raise questions concerning their possible clinical significance that are still under investigation.

How cyanophage S-2L rejects adenine and incorporates 2-aminoadenine to saturate hydrogen bonding in its DNA

Bacteriophages have long been known to use modified bases in their DNA to prevent cleavage by the host's restriction endonucleases. Among them, cyanophage S-2L is unique because its genome has all its adenines (A) systematically replaced by 2-aminoadenines (Z). Here, we identify a member of the PrimPol family as the sole possible polymerase of S-2L and we find it can incorporate both A and Z in front of a T. Its crystal structure at 1.5 Å resolution confirms that there is no structural element in the active site that could lead to the rejection of A in front of T. To resolve this contradiction, we show that a nearby gene is a triphosphohydolase specific of dATP (DatZ), that leaves intact all other dNTPs, including dZTP. This explains the absence of A in S-2L genome. Crystal structures of DatZ with various ligands, including one at sub-angstrom resolution, allow to describe its mechanism as a typical two-metal-ion mechanism and to set the stage for its engineering.

Comparative Analysis of the Circular and Highly Asymmetrical Marseilleviridae Genomes

Marseilleviridae members are large dsDNA viruses with icosahedral particles 250 nm in diameter infecting Acanthamoeba. Their 340 to 390 kb genomes encode 450 to 550 protein-coding genes. Since the discovery of marseillevirus (the prototype of the family) in 2009, several strains were isolated from various locations, among which 13 are now fully sequenced. This allows the organization of their genomes to be deciphered through comparative genomics. Here, we first experimentally demonstrate that the Marseilleviridae genomes are circular. We then acknowledge a strong bias in sequence conservation, revealing two distinct genomic regions. One gathers most Marseilleviridae paralogs and has undergone genomic rearrangements, while the other, enriched in core genes, exhibits the opposite pattern. Most of the genes whose protein products compose the viral particles are located in the conserved region. They are also strongly biased toward a late gene expression pattern. We finally discuss the potential advantages of Marseilleviridae having a circular genome, and the possible link between the biased distribution of their genes and the transcription as well as DNA replication mechanisms that remain to be characterized.

Fundamental Difficulties Prevent the Reconstruction of the Deep Phylogeny of Viruses

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