Degenerate PCR Primers to Reveal the Diversity of Giant Viruses in Coastal Waters

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.

vAMPirus: A versatile amplicon processing and analysis program for studying viruses

Amplicon sequencing is an effective and increasingly applied method for studying viral communities in the environment. Here, we present vAMPirus, a user-friendly, comprehensive, and versatile DNA and RNA virus amplicon sequence analysis program, designed to support investigators in exploring virus amplicon sequencing data and running informed, reproducible analyses. vAMPirus intakes raw virus amplicon libraries and, by default, performs nucleotide- and amino acid-based analyses to produce results such as sequence abundance information, taxonomic classifications, phylogenies and community diversity metrics. The vAMPirus analytical framework leverages 16 different opensource tools and provides optional approaches that can increase the ratio of biological signal-to-noise and thereby reveal patterns that would have otherwise been masked. Here, we validate the vAMPirus analytical framework and illustrate its implementation as a general virus amplicon sequencing workflow by recapitulating findings from two previously published double-stranded DNA virus datasets. As a case study, we also apply the program to explore the diversity and distribution of a coral reef-associated RNA virus. vAMPirus is streamlined within Nextflow, offering straightforward scalability, standardization and communication of virus lineage-specific analyses. The vAMPirus framework is designed to be adaptable; community-driven analytical standards will continue to be incorporated as the field advances. vAMPirus supports researchers in revealing patterns of virus diversity and population dynamics in nature, while promoting study reproducibility and comparability.

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

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

From Mimivirus to Mirusvirus: The Quest for Hidden Giants

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

Virome reveals effect of Ulva prolifera green tide on the structural and functional profiles of virus communities in coastal environments

Viruses are widely distributed in marine environments, where they influence the transformation of matter and energy by modulating host metabolism. Driven by eutrophication, green tides are a rising concern in Chinese coastal areas, and are a serious ecological disaster that negatively affects coastal ecosystems and disrupts biogeochemical cycles. Although the composition of bacterial communities in green algae has been investigated, the diversity and roles of viruses in green algal blooms are largely unexplored. Therefore, the diversity, abundance, lifestyle, and metabolic potential of viruses in a natural bloom in Qingdao coastal area were investigated at three different stages (pre-bloom, during-bloom, and post-bloom) by metagenomics analysis. The dsDNA viruses, Siphoviridae, Myoviridae, Podoviridae, and Phycodnaviridae, were found to dominate the viral community. The viral dynamics exhibited distinct temporal patterns across different stages. The composition of the viral community varied during the bloom, especially in populations with low abundance. The lytic cycle was most predominant, and the abundance of lytic viruses increased slightly in the post-bloom stage. The diversity and richness of the viral communities varied distinctly during the green tide, and the post-bloom stage favored viral diversity and richness. The total organic carbon, dissolved oxygen, NO^3-, NO^2-, PO₄^3-, chlorophyll-a contents, and temperature variably co-influenced the viral communities. The primary hosts included bacteria, algae, and other microplankton. Network analysis revealed the closer links between the viral communities as the bloom progressed. Functional prediction revealed that the viruses possibly influenced the biodegradation of microbial hydrocarbons and carbon by metabolic augmentation via auxiliary metabolic genes. The composition, structure, metabolic potential, and interaction taxonomy of the viromes differed significantly across the different stages of the green tide. The study demonstrated that the ecological event shaped the viral communities during algal bloom, and the viral communities played a significant role in phycospheric microecology.

Covariance of Marine Nucleocytoplasmic Large DNA Viruses with Eukaryotic Plankton Communities in the Sub-Arctic Kongsfjorden Ecosystem: A Metagenomic Analysis of Marine Microbial Ecosystems

Nucleocytoplasmic large DNA viruses (NCLDVs) infect various marine eukaryotes. However, little is known about NCLDV diversity and their relationships with eukaryotic hosts in marine environments, the elucidation of which will advance the current understanding of marine ecosystems. This study characterizes the interplay between NCLDVs and the eukaryotic plankton community (EPC) in the sub-Arctic area using metagenomics and metabarcoding to investigate NCLDVs and EPC, respectively, in the Kongsfjorden ecosystem of Svalbard (Norway) in April and June 2018. Gyrodinium helveticum (Dinophyceae) is the most prevalent eukaryotic taxon in the EPC in April, during which time Mimiviridae (31.8%), Poxviridae (25.1%), Phycodnaviridae (14.7%) and Pandoraviridae (13.1%) predominate. However, in June, the predominant taxon is Aureococcus anophagefferens (Pelagophyceae), and the NCLDVs, Poxviridae (32.9%), Mimiviridae (29.1%), and Phycodnaviridae (18.5%) appear in higher proportions with an increase in Pelagophyceae, Bacillariophyceae, and Chlorophyta groups. Thus, differences in NCLDVs may be caused by changes in EPC composition in response to environmental changes, such as increases in water temperature and light intensity. Taken together, these findings are particularly relevant considering the anticipated impact of NCLDV-induced EPC control mechanisms on polar regions and, therefore, improve the understanding of the Sub-Arctic Kongsfjorden ecosystem.

Unravelling the ecological complexity of soil viromes: Challenges and opportunities

Viruses are extremely abundant and ubiquitous in soil, and significantly contribute to various terrestrial ecosystem processes such as biogeochemical nutrient cycling, microbiome regulation and community assembly, and host evolutionary dynamics. Despite their numerous dominance and functional importance, understanding soil viral ecology is a formidable challenge, because of the technological challenges to characterize the abundance, diversity and community compositions of viruses, and their interactions with other organisms in the complex soil environment. Viruses may engage in a myriad of biological interactions within soil food webs across a broad range of spatiotemporal scales and are exposed to various biotic and abiotic disturbances. Current studies on the soil viromes, however, often describe the complexity of their tremendous diversity, but lack of exploring their potential ecological roles. In this article, we summarized the major methods to decipher the ecology of soil viruses, discussed biotic and abiotic factors and global change factors that shape the diversity and composition of soil viromes, and the ecological roles of soil viruses. We also proposed a new framework to understand the ecological complexity of viruses from micro to macro ecosystem scales and to predict and unravel their activities in terrestrial ecosystems.

Degenerate PCR Targeting the Major Capsid Protein Gene of HcRNAV and Related Viruses

Heterocapsa circularisquama RNA virus (HcRNAV) is the only dinoflagellate-infecting RNA virus that has been isolated to date. We herein investigated the diversity of the major capsid protein gene of HcRNAV and related viruses using degenerate PCR and in silico ana-lyses. Diverse sequences related to HcRNAV were successfully amplified from marine sediments. Amplicons contained conserved and variable regions; the latter were predicted to be located on the outer surface of the capsid. Our approach provides insights into the diversity of viruses that are difficult to isolate in the environment and will enhance rapidly growing metagenome sequence repositories.

RNA Sequencing of Medusavirus Suggests Remodeling of the Host Nuclear Environment at an Early Infection Stage

Viruses of the phylum Nucleocytoviricota, or nucleo-cytoplasmic large DNA viruses (NCLDVs), undergo a cytoplasmic or nucleo-cytoplasmic cycle, the latter of which involves both nuclear and cytoplasmic compartments to proceed viral replication. Medusavirus, a recently isolated NCLDV, has a nucleo-cytoplasmic replication cycle in amoebas during which the host nuclear membrane apparently remains intact, a unique feature among amoeba-infecting NCLDVs. The medusavirus genome lacks most transcription genes but encodes a full set of histone genes. To investigate its infection strategy, we performed a time course RNA sequencing (RNA-seq) experiment. All viral genes were transcribed and classified into five temporal expression clusters. The immediate early genes (cluster 1, 42 genes) were mostly (83%) of unknown functions, frequently (95%) associated with a palindromic promoter-like motif, and often (45%) encoded putative nucleus-localized proteins. These results suggest massive reshaping of the host nuclear environment by viral proteins at an early stage of infection. Genes in other expression clusters (clusters 2 to 5) were assigned to various functional categories. The virally encoded core histone genes were in cluster 3, whereas the viral linker histone H1 gene was in cluster 1, suggesting they have distinct roles during the course of the virus infection. The transcriptional profile of the host Acanthamoeba castellanii genes was greatly altered postinfection. Several encystment-related host genes showed increased representation levels at 48 h postinfection, which is consistent with the previously reported amoeba encystment upon medusavirus infection. IMPORTANCE Medusavirus is an amoeba-infecting giant virus that was isolated from a hot spring in Japan. It belongs to the proposed family "Medusaviridae" in the phylum Nucleocytoviricota. Unlike other amoeba-infecting giant viruses, medusavirus initiates its DNA replication in the host nucleus without disrupting the nuclear membrane. Our RNA sequencing (RNA-seq) analysis of its infection course uncovered ordered viral gene expression profiles. We identified temporal expression clusters of viral genes and associated putative promoter motifs. The subcellular localization prediction showed a clear spatiotemporal correlation between gene expression timing and localization of the encoded proteins. Notably, the immediate early expression cluster was enriched in genes targeting the nucleus, suggesting the priority of remodeling the host intranuclear environment during infection. The transcriptional profile of amoeba genes was greatly altered postinfection.

High Transcriptional Activity and Diverse Functional Repertoires of Hundreds of Giant Viruses in a Coastal Marine System

Viruses belonging to the Nucleocytoviricota phylum are globally distributed and include members with notably large genomes and complex functional repertoires. Recent studies have shown that these viruses are particularly diverse and abundant in marine systems, but the magnitude of actively replicating Nucleocytoviricota present in ocean habitats remains unclear. In this study, we compiled a curated database of 2,431 Nucleocytoviricota genomes and used it to examine the gene expression of these viruses in a 2.5-day metatranscriptomic time-series from surface waters of the California Current. We identified 145 viral genomes with high levels of gene expression, including 90 Imitervirales and 49 Algavirales viruses. In addition to recovering high expression of core genes involved in information processing that are commonly expressed during viral infection, we also identified transcripts of diverse viral metabolic genes from pathways such as glycolysis, the TCA cycle, and the pentose phosphate pathway, suggesting that virus-mediated reprogramming of central carbon metabolism is common in oceanic surface waters. Surprisingly, we also identified viral transcripts with homology to actin, myosin, and kinesin domains, suggesting that viruses may use these gene products to manipulate host cytoskeletal dynamics during infection. We performed phylogenetic analysis on the virus-encoded myosin and kinesin proteins, which demonstrated that most belong to deep-branching viral clades, but that others appear to have been acquired from eukaryotes more recently. Our results highlight a remarkable diversity of active Nucleocytoviricota in a coastal marine system and underscore the complex functional repertoires expressed by these viruses during infection. IMPORTANCE The discovery of giant viruses has transformed our understanding of viral complexity. Although viruses have traditionally been viewed as filterable infectious agents that lack metabolism, giant viruses can reach sizes rivalling cellular lineages and possess genomes encoding central metabolic processes. Recent studies have shown that giant viruses are widespread in aquatic systems, but the activity of these viruses and the extent to which they reprogram host physiology in situ remains unclear. Here, we show that numerous giant viruses consistently express central metabolic enzymes in a coastal marine system, including components of glycolysis, the TCA cycle, and other pathways involved in nutrient homeostasis. Moreover, we found expression of several viral-encoded actin, myosin, and kinesin genes, indicating viral manipulation of the host cytoskeleton during infection. Our study reveals a high activity of giant viruses in a coastal marine system and indicates they are a diverse and underappreciated component of microbial diversity in the ocean.

Quantitative Assessment of Nucleocytoplasmic Large DNA Virus and Host Interactions Predicted by Co-occurrence Analyses

Nucleocytoplasmic large DNA viruses (NCLDVs) are highly diverse and abundant in marine environments. However, the knowledge of their hosts is limited because only a few NCLDVs have been isolated so far. Taking advantage of the recent large-scale marine metagenomics census, in silico host prediction approaches are expected to fill the gap and further expand our knowledge of virus-host relationships for unknown NCLDVs. In this study, we built co-occurrence networks of NCLDVs and eukaryotic taxa to predict virus-host interactions using Tara Oceans sequencing data. Using the positive likelihood ratio to assess the performance of host prediction for NCLDVs, we benchmarked several co-occurrence approaches and demonstrated an increase in the odds ratio of predicting true positive relationships 4-fold compared to random host predictions. To further refine host predictions from high-dimensional co-occurrence networks, we developed a phylogeny-informed filtering method, Taxon Interaction Mapper, and showed it further improved the prediction performance by 12-fold. Finally, we inferred virophage-NCLDV networks to corroborate that co-occurrence approaches are effective for predicting interacting partners of NCLDVs in marine environments.IMPORTANCE NCLDVs can infect a wide range of eukaryotes, although their life cycle is less dependent on hosts compared to other viruses. However, our understanding of NCLDV-host systems is highly limited because few of these viruses have been isolated so far. Co-occurrence information has been assumed to be useful to predict virus-host interactions. In this study, we quantitatively show the effectiveness of co-occurrence inference for NCLDV host prediction. We also improve the prediction performance with a phylogeny-guided method, which leads to a concise list of candidate host lineages for three NCLDV families. Our results underpin the usage of co-occurrence approaches for the metagenomic exploration of the ecology of this diverse group of viruses.

Tara Oceans: towards global ocean ecosystems biology

A planetary-scale understanding of the ocean ecosystem, particularly in light of climate change, is crucial. Here, we review the work of Tara Oceans, an international, multidisciplinary project to assess the complexity of ocean life across comprehensive taxonomic and spatial scales. Using a modified sailing boat, the team sampled plankton at 210 globally distributed sites at depths down to 1,000 m. We describe publicly available resources of molecular, morphological and environmental data, and discuss how an ecosystems biology approach has expanded our understanding of plankton diversity and ecology in the ocean as a planetary, interconnected ecosystem. These efforts illustrate how global-scale concepts and data can help to integrate biological complexity into models and serve as a baseline for assessing ecosystem changes and the future habitability of our planet in the Anthropocene epoch.

Dynamic genome evolution and complex virocell metabolism of globally-distributed giant viruses

The discovery of eukaryotic giant viruses has transformed our understanding of the limits of viral complexity, but the extent of their encoded metabolic diversity remains unclear. Here we generate 501 metagenome-assembled genomes of Nucleo-Cytoplasmic Large DNA Viruses (NCLDV) from environments around the globe, and analyze their encoded functional capacity. We report a remarkable diversity of metabolic genes in widespread giant viruses, including many involved in nutrient uptake, light harvesting, and nitrogen metabolism. Surprisingly, numerous NCLDV encode the components of glycolysis and the TCA cycle, suggesting that they can re-program fundamental aspects of their host's central carbon metabolism. Our phylogenetic analysis of NCLDV metabolic genes and their cellular homologs reveals distinct clustering of viral sequences into divergent clades, indicating that these genes are virus-specific and were acquired in the distant past. Overall our findings reveal that giant viruses encode complex metabolic capabilities with evolutionary histories largely independent of cellular life, strongly implicating them as important drivers of global biogeochemical cycles.

An Optimized Metabarcoding Method for Mimiviridae

Mimiviridae is a group of viruses with large genomes and virions. Ecological relevance of Mimiviridae in marine environments has been increasingly recognized through the discoveries of novel isolates and metagenomic studies. To facilitate ecological profiling of Mimiviridae, we previously proposed a meta-barcoding approach based on 82 degenerate primer pairs (i.e., MEGAPRIMER) targeting the DNA polymerase gene of Mimiviridae. The method detected a larger number of operational taxonomic units (OTUs) in environmental samples than previous methods. However, it required large quantities of DNA and was laborious due to the use of individual primer pairs. Here, we examined coastal seawater samples using varying PCR conditions and purification protocols to streamline the MEGAPRIMER method. Mixing primer pairs in "cocktails" reduced the required amount of environmental DNA by 90%, while reproducing the results obtained by the original protocol. We compared the results obtained by the meta-barcoding approach with quantifications using qPCR for selected OTUs. This revealed possible amplification biases among different OTUs, but the frequency profiles for individual OTUs across multiple samples were similar to those obtained by qPCR. We anticipate that the newly developed MEGAPRIMER protocols will be useful for ecological investigation of Mimiviridae in a larger set of environmental samples.

Targeted metagenomic recovery of four divergent viruses reveals shared and distinctive characteristics of giant viruses of marine eukaryotes

Giant viruses have remarkable genomic repertoires-blurring the line with cellular life-and act as top-down controls of eukaryotic plankton. However, to date only six cultured giant virus genomes are available from the pelagic ocean. We used at-sea flow cytometry with staining and sorting designed to target wild predatory eukaryotes, followed by DNA sequencing and assembly, to recover novel giant viruses from the Pacific Ocean. We retrieved four 'PacV' partial genomes that range from 421 to 1605 Kb, with 13 contigs on average, including the largest marine viral genomic assembly reported to date. Phylogenetic analyses indicate that three of the new viruses span a clade with deep-branching members of giant Mimiviridae, incorporating the Cafeteria roenbergensis virus, the uncultivated terrestrial Faunusvirus, one PacV from a choanoflagellate and two PacV with unclear hosts. The fourth virus, oPacV-421, is phylogenetically related to viruses that infect haptophyte algae. About half the predicted proteins in each PacV have no matches in NCBI nr (e-value < 10-5), totalling 1735 previously unknown proteins; the closest affiliations of the other proteins were evenly distributed across eukaryotes, prokaryotes and viruses of eukaryotes. The PacVs encode many translational proteins and two encode eukaryotic-like proteins from the Rh family of the ammonium transporter superfamily, likely influencing the uptake of nitrogen during infection. cPacV-1605 encodes a microbial viral rhodopsin (VirR) and the biosynthesis pathway for the required chromophore, the second finding of a choanoflagellate-associated virus that encodes these genes. In co-collected metatranscriptomes, 85% of cPacV-1605 genes were expressed, with capsids, heat shock proteins and proteases among the most highly expressed. Based on orthologue presence-absence patterns across the PacVs and other eukaryotic viruses, we posit the observed viral groupings are connected to host lifestyles as heterotrophs or phototrophs. This article is part of a discussion meeting issue 'Single cell ecology'.

The Earth Is Small for "Leviathans": Long Distance Dispersal of Giant Viruses across Aquatic Environments

Giant viruses of ‘Megaviridae’ have the ability to widely disperse around the globe. We herein examined ‘Megaviridae’ communities in four distinct aquatic environments (coastal and offshore seawater, brackish water, and hot spring freshwater), which are distantly located from each other (between 74 and 1,765 km), using a meta-barcoding method. We identified between 593 and 3,627 OTUs in each sample. Some OTUs were detected in all five samples tested as well as in many of the Tara Oceans metagenomes, suggesting the existence of viruses of this family in a wide range of habitats and the ability to circulate on the planet.

Fifteen Marseilleviruses Newly Isolated From Three Water Samples in Japan Reveal Local Diversity of Marseilleviridae

The family Marseilleviridae, defined as a group of icosahedral double-stranded DNA viruses with particle size of approximately 250 nm and genome size of 350-380 kbp, belongs to the nucleo-cytoplasmic family of large DNA viruses. The family Marseilleviridae is currently classified into lineages A-E. In this study, we isolated 12 or 15 new members of the family Marseilleviridae from three sampling locations in Japan. Molecular phylogenetic analysis of the MCP genes showed that the new viruses could be further classified into three groups, hokutoviruses, kashiwazakiviruses, and kyotoviruses. Hokutoviruses were closely related to lineage B, kyotoviruses were related to lineage A, and kashiwazakiviruses were also classified into lineage B but a new putative subgroup of lineage B, revealing the diversity of this lineage. Interestingly, more than two viruses with slightly different MCP genes were isolated from a single water sample from a single location, i.e., two hokutoviruses and one kashiwazakivirus were isolated from a small reservoir, five kashiwazakiviruses from the mouth of a river, and five kyotoviruses from fresh water of a river, suggesting that several milliliters of water samples contain several types of giant viruses. Amoeba cells infected with hokutoviruses or kashiwazakiviruses exhibited a "bunch" formation consisting of normal and infected cells similarly to a tupanvirus, whereas cells infected with kyotoviruses or tokyovirus did not. These results suggest the previously unrecognized local diversity of the family Marseilleviridae in aquatic environments.