Biogeography of Viruses in the Sea

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.

Molecular characterization of plasma virome of hepatocellular carcinoma (HCC) patients

Hepatocellular carcinoma (HCC) stands as the most common cancer type, arising from various causes, and responsible for a substantial number of cancer-related fatalities. Recent advancements in viral metagenomics have empowered scientists to delve into the intricate diversity of the virosphere, viral evolution, interactions between viruses and their hosts, and the identification of viral causes behind disease outbreaks, the development of specific symptoms, and their potential role in altering the host's physiology. The present study had the objective of "Molecular Characterization of HBV, HCV, anelloviruses, CMV, SENV-D, SENV-H, HEV, and HPV viruses among individuals suffering from HCC." A total of 381 HCC patients contributed 10 cc of blood each for this study. The research encompassed the assessment of tumor markers, followed by molecular characterization of HBV, HCV, Anelloviruses (TTV, TTMV, and TTMDV), SENV-H and SENV-D viruses, HEV, CMV, and HPV, as well as histopathological examinations. The outcomes of this study revealed that majority of the HCC patients 72.4% (276/381) were male as compared to females. HCV infection, at 76.4% (291 out of 381), exhibited a significant association (p < 0.05) with HCC. Most patients displayed singular lesions in the liver, with Child Pugh Score Type B being the predominant finding in 45.2% of cases. Plasma virome analysis indicated the prevalence of TTMDV (75%), followed by TTMV (70%) and TTV (42.1%) among anelloviruses in HCC patients. Similarly, SENV-H (52%) was followed by SENV-D (20%), with co-infections at 15%. The presence of CMV and HEV among the HCC patients was recorded 5% each however 3.5% of the patients showed the presence of HPV. In conclusion, this study underscores that HCC patients serve as reservoirs for various pathogenic and non-pathogenic viruses, potentially contributing to the development, progression, and severity of the disease.

The structure and assembly mechanisms of T4-like cyanophages community in the South China Sea

Marine ecosystems contain an immense diversity of phages, many of which infect cyanobacteria (cyanophage) that are largely responsible for primary productivity. To characterize the genetic diversity and biogeographic distribution of the marine T4-like cyanophage community in the northern South China Sea, the T4-like cyanophage portal protein gene (g20) was amplified. Phylogenetic analysis revealed that marine T4-like cyanophages were highly diverse, with g20 operational taxonomic units being affiliated with five defined clades (Clusters I-V). Cluster II had a wide geographic distribution, Cluster IV was the most abundant in the open sea, and Cluster I was dominant in coastal shelf environments. Our results showed T4-like cyanophages (based on g20) community was generally shaped via heterogeneous selection. Highly variable environmental factors (such as salinity and temperature) can heterogeneously select different cyanophage communities. Nevertheless, the dominant drivers of the T4-like cyanophage community based on the g20 and g23 (T4-like phage major capsid protein gene) were different, probably due to different coverages by the primer sets. Furthermore, the community assembly processes of T4-like cyanophages were affected by host traits (abundance and distribution), viral traits (latent period, burst size, and host range), and environmental properties (temperature and salinity).IMPORTANCECyanophages are abundant and ubiquitous in the oceans, altering population structures and evolution of cyanobacteria, which account for a large portion of global carbon fixation, through host mortality, horizontal gene transfer, and the modulation of host metabolism. However, little is known about the biogeography and ecological drivers that shape the cyanophage community. Here, we use g20 and g23 genes to examine the biogeographic patterns and the assembly mechanisms of T4-like cyanophage community in the northern part of the South China Sea. The different coverages of primer sets might lead to the different dominant drivers of T4-like cyanophage community based on g20 and g23 genes. Our results showed that characteristics of viral traits (latent period, burst size, and host range) and host traits (abundance and distribution) were found to either limit or enhance the biogeographic distribution of T4-like cyanophages. Overall, both virus and host properties are critical to consider when determining rules of community assembly for viruses.

Environmental viromes reveal global virosphere of deep-sea sediment RNA viruses

INTRODUCTION

Viruses are the most abundant and diverse life forms on the earth. Both DNA viruses and RNA viruses play important roles in marine ecosystems via regulating biogeochemical cycles.

OBJECTIVES

However, the virome of marine RNA viruses has been rarely explored so far. In this study, therefore, the environmental viromes of deep-sea sediment RNA viruses were characterized on a global scale to reveal the global virosphere of deep-sea RNA viruses.

METHODS

The viral particles were purified from each of 133 deep-sea sediment samples and then characterized based on metagenomes of RNA viruses.

RESULTS

In this study, we established the global virome dataset of deep-sea RNA viruses purified from 133 sediment samples that were collected from typical deep-sea ecosystems of three oceans. A total of 85,059 viral operational taxonomic units (vOTUs) were identified, of which only 1.72% were hitherto known, indicating that the deep-sea sediment is a repository of novel RNA viruses. These vOTUs were classified into 20 viral families, including prokaryotic (7.09%) and eukaryotic (65.81%) RNA viruses. Furthermore, 1,463 deep-sea RNA viruses with complete genomes were obtained. The differentiation of RNA viral communities was driven by the deep-sea ecosystems as opposed to geographical region. Specifically, the virus-encoded metabolic genes took great effects on the differentiation of RNA viral communities by mediating the energy metabolism in the deep-sea ecosystems.

CONCLUSIONS

Therefore, our findings indicate that the deep sea is a vast reservoir of novel RNA viruses for the first time, and the differentiation of RNA viral communities is driven by the deep-sea ecosystems through energy metabolism.

CRISPR-Cas immunity is repressed by the LysR-type transcriptional regulator PigU

Bacteria protect themselves from infection by bacteriophages (phages) using different defence systems, such as CRISPR-Cas. Although CRISPR-Cas provides phage resistance, fitness costs are incurred, such as through autoimmunity. CRISPR-Cas regulation can optimise defence and minimise these costs. We recently developed a genome-wide functional genomics approach (SorTn-seq) for high-throughput discovery of regulators of bacterial gene expression. Here, we applied SorTn-seq to identify loci influencing expression of the two type III-A Serratia CRISPR arrays. Multiple genes affected CRISPR expression, including those involved in outer membrane and lipopolysaccharide synthesis. By comparing loci affecting type III CRISPR arrays and cas operon expression, we identified PigU (LrhA) as a repressor that co-ordinately controls both arrays and cas genes. By repressing type III-A CRISPR-Cas expression, PigU shuts off CRISPR-Cas interference against plasmids and phages. PigU also represses interference and CRISPR adaptation by the type I-F system, which is also present in Serratia. RNA sequencing demonstrated that PigU is a global regulator that controls secondary metabolite production and motility, in addition to CRISPR-Cas immunity. Increased PigU also resulted in elevated expression of three Serratia prophages, indicating their likely induction upon sensing PigU-induced cellular changes. In summary, PigU is a major regulator of CRISPR-Cas immunity in Serratia.

Megataxonomy and global ecology of the virosphere

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

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

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

Assessing the biogeography of marine giant viruses in four oceanic transects

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

Spatial and temporal metagenomics of river compartments reveals viral community dynamics in an urban impacted stream

Although river ecosystems comprise less than 1% of Earth's total non-glaciated area, they are critical modulators of microbially and virally orchestrated global biogeochemical cycles. However, most studies either use data that is not spatially resolved or is collected at timepoints that do not reflect the short life cycles of microorganisms. As a result, the relevance of microbiome interactions and the impacts they have over time on biogeochemical cycles are poorly understood. To assess how viral and microbial communities change over time, we sampled surface water and pore water compartments of the wastewater-impacted River Erpe in Germany every 3 hours over a 48-hour period resulting in 32 metagenomes paired to geochemical and metabolite measurements. We reconstructed 6,500 viral and 1,033 microbial genomes and found distinct communities associated with each river compartment. We show that 17% of our vMAGs clustered to viruses from other ecosystems like wastewater treatment plants and rivers. Our results also indicated that 70% of the viral community was persistent in surface waters, whereas only 13% were persistent in the pore waters taken from the hyporheic zone. Finally, we predicted linkages between 73 viral genomes and 38 microbial genomes. These putatively linked hosts included members of the Competibacteraceae, which we suggest are potential contributors to carbon and nitrogen cycling. Together, these findings demonstrate that microbial and viral communities in surface waters of this urban river can exist as stable communities along a flowing river; and raise important considerations for ecosystem models attempting to constrain dynamics of river biogeochemical cycles.

Spatiotemporal Shift of T4-Like Phage Community Structure in the Three Largest Estuaries of China

Estuaries are one of the most highly productive and economically important ecosystems at the continent-ocean interface. Estuary productivity is largely determined by the microbial community structure and activity. Viruses are major agents of microbial mortality and are key drivers of global geochemical cycles. However, the taxonomic diversity of viral communities and their spatial-temporal distribution in estuarine ecosystems have been poorly studied. In this study, we investigated the T4-like viral community composition at three major Chinese estuaries in winter and in summer. Diverse T4-like viruses, which were divided into three main clusters (Clusters I to III), were revealed. The Marine Group of Cluster III, with seven identified subgroups, was the most dominant (averaging 76.5% of the total sequences) in the Chinese estuarine ecosystems. Significant variations of T4-like viral community composition were observed among estuaries and seasons, with higher diversity occurring in winter. Among various environmental variables, temperature was a main driver of the viral communities. This study demonstrates viral assemblage diversification and seasonality in Chinese estuarine ecosystems. IMPORTANCE Viruses are ubiquitous but largely uncharacterized members of aquatic environments that cause significant mortality in microbial communities. Recent large-scale oceanic projects have greatly advanced our understanding of viral ecology in marine environments, but those studies mostly focused on oceanic regions. There have yet to be spatiotemporal studies of viral communities in estuarine ecosystems, which are unique habitats that play a significant role in global ecology and biogeochemistry. This work is the first comprehensive study that provides a detailed picture of the spatial and seasonal variation of viral communities (specifically, T4-like viral communities) in three major estuarine ecosystems in China. These findings provide much-needed knowledge regarding estuarine viral ecosystems, which currently lags in oceanic ecosystem research.

High Level of Interaction between Phages and Bacteria in an Artisanal Raw Milk Cheese Microbial Community

Microbial starter cultures are used in the production of many cheeses around the world, such as Parmigiano-Reggiano, in Italy, Époisses, in France, and Canastra, in Brazil, providing many of the unique features of these cheeses. Bacteriophages (phages) are ubiquitous and well known to modulate the structure of bacterial communities, and recent data indicate that cheeses contain a high abundance of naturally occurring phages. Here, we analyze the viral and bacterial metagenomes of Canastra cheese: a traditional artisanal Brazilian cheese produced using an endogenous starter culture and raw milk. Over 1,200 viral operational taxonomic units were recovered using both isolated viral-like particles and complete metagenomic DNA. Common viral families identified included Siphoviridae and Myoviridae, with 40% of putative phage genomes unidentified at the family level of classification. We observed very high phage diversity, which varied greatly across different cheese producers, with 28% of phage genomes detected in only one producer. Several metagenome-assembled genomes were recovered for lactic acid-producing bacteria, as well as nonstarter bacterial species, and we identified several phage-bacterium interactions, at the strain level of resolution, varying across distinct cheese producers. We postulate that at least one bacterial strain detected could be endogenous and unique to the Canastra cheese-producing region in Brazil and that its growth seems to be modulated by autochthonous phages present in this artisanal production system. This phage-host relationship is likely to influence the fermentation dynamics and ultimately the sensorial profile of these cheeses, with implications for other similar cheese production systems around the world. IMPORTANCE Our work demonstrated a dynamic yet stable microbial ecosystem during cheese production using an endogenous starter culture. This was observed across several distinct producers and was marked by genomic evidence of continued phage-bacterium interactions, such as the presence of bacterial defense mechanisms. Furthermore, we provide evidence of unique microbial signatures for each individual cheese producer studied in the region, a fact that may have profound consequences on product traceability. This was the first effort to describe and understand the bacteriophage composition and ecological dynamics within the Brazilian Canastra cheese production system. The study of this prototypical backslopping production system provides a solid background for further mechanistic studies of the production of many cheeses around the world.

Diverse Marine T4-like Cyanophage Communities Are Primarily Comprised of Low-Abundance Species Including Species with Distinct Seasonal, Persistent, Occasional, or Sporadic Dynamics

Cyanophages exert important top-down controls on their cyanobacteria hosts; however, concurrent analysis of both phage and host populations is needed to better assess phage-host interaction models. We analyzed picocyanobacteria Prochlorococcus and Synechococcus and T4-like cyanophage communities in Pacific Ocean surface waters using five years of monthly viral and cellular fraction metagenomes. Cyanophage communities contained thousands of mostly low-abundance (<2% relative abundance) species with varying temporal dynamics, categorized as seasonally recurring or non-seasonal and occurring persistently, occasionally, or sporadically (detected in ≥85%, 15-85%, or <15% of samples, respectively). Viromes contained mostly seasonal and persistent phages (~40% each), while cellular fraction metagenomes had mostly sporadic species (~50%), reflecting that these sample sets capture different steps of the infection cycle-virions from prior infections or within currently infected cells, respectively. Two groups of seasonal phages correlated to Synechococcus or Prochlorococcus were abundant in spring/summer or fall/winter, respectively. Cyanophages likely have a strong influence on the host community structure, as their communities explained up to 32% of host community variation. These results support how both seasonally recurrent and apparent stochastic processes, likely determined by host availability and different host-range strategies among phages, are critical to phage-host interactions and dynamics, consistent with both the Kill-the-Winner and the Bank models.

Cross-Hemispheric Genetic Diversity and Spatial Genetic Structure of Callinectes sapidus Reovirus 1 (CsRV1)

The movement of viruses in aquatic systems is rarely studied over large geographic scales. Oceanic currents, host migration, latitude-based variation in climate, and resulting changes in host life history are all potential drivers of virus connectivity, adaptation, and genetic structure. To expand our understanding of the genetic diversity of Callinectes sapidus reovirus 1 (CsRV1) across a broad spatial and host life history range of its blue crab host (Callinectes sapidus), we obtained 22 complete and 96 partial genomic sequences for CsRV1 strains from the US Atlantic coast, Gulf of Mexico, Caribbean Sea, and the Atlantic coast of South America. Phylogenetic analyses of CsRV1 genomes revealed that virus genotypes were divided into four major genogroups consistent with their host geographic origins. However, some CsRV1 sequences from the US mid-Atlantic shared high genetic similarity with the Gulf of Mexico genotypes, suggesting potential human-mediated movement of CsRV1 between the US mid-Atlantic and Gulf coasts. This study advances our understanding of how climate, coastal geography, host life history, and human activity drive patterns of genetic structure and diversity of viruses in marine animals and contributes to the capacity to infer broadscale host population connectivity in marine ecosystems from virus population genetic data.

Abundant and cosmopolitan lineage of cyanopodoviruses lacking a DNA polymerase gene

Cyanopodoviruses affect the mortality and population dynamics of the unicellular picocyanobacteria Prochlorococcus and Synechococcus, the dominant primary producers in the oceans. Known cyanopodoviruses all contain the DNA polymerase gene (DNA pol) that is important for phage DNA replication and widely used in field quantification and diversity studies. However, we isolated 18 cyanopodoviruses without identifiable DNA pol. They form a new MPP-C clade that was separated from the existing MPP-A, MPP-B, and P-RSP2 clades. The MPP-C phages have the smallest genomes (37.3-37.9 kb) among sequenced cyanophages, and show longer latent periods than the MPP-B phages. Metagenomic reads of both clades are highly abundant in surface waters, but the MPP-C phages show higher relative abundance in surface waters than in deeper waters, while MPP-B phages have higher relative abundance in deeper waters. Our study reveals that cyanophages with distinct genomic contents and infection kinetics can exhibit different depth profiles in the oceans.

Assessing the biogeography of marine giant viruses in four oceanic transects

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

Genetic Diversity and Characterization of Circular Replication (Rep)-Encoding Single-Stranded (CRESS) DNA Viruses

The CRESS-DNA viruses are the ubiquitous virus detected in almost all eukaryotic life trees and play an essential role in the maintaining ecosystem of the globe. Still, their genetic diversity is not fully understood. Here, we bring to light the genetic diversity of replication (Rep) and capsid (Cap) proteins of CRESS-DNA viruses. We divided the Rep protein of the CRESS-DNA virus into 10 clusters using CLANS and phylogenetic analyses. Also, most of the Rep protein in Rep cluster 1 (R1) and R2 (Circoviridae, Smacoviridae, Nanoviridae, and CRESSV1-5) contain the Viral_Rep superfamily and P-loop_NTPase superfamily domains, while the Rep protein of viruses in other clusters has no such characterized functional domain. The Circoviridae, Nanoviridae, and CRESSV1-3 viruses contain two domains, such as Viral_Rep and P-loop_NTPase; the CRESSV4 and CRESSV5 viruses have only the Viral_Rep domain; most of the sequences in the pCRESS-related group have only P-loop_NTPase; and Smacoviridae do not have these two domains. Further, we divided the Cap protein of the CRESS-DNA virus into 20 clusters using CLANS and phylogenetic analyses. The Rep and Cap proteins of Circoviridae and Smacoviridae are grouped into a specific cluster. Cap protein of CRESS-DNA viruses grouped with one cluster and Rep protein with another cluster. Further, our study reveals that selection pressure plays a significant role in the evolution of CRESS-DNA viruses' Rep and Cap genes rather than mutational pressure. We hope this study will help determine the genetic diversity of CRESS-DNA viruses as more sequences are discovered in the future. IMPORTANCE The genetic diversity of CRESS-DNA viruses is not fully understood. CRESS-DNA viruses are classified as CRESSV1 to CRESSV6 using only Rep protein. This study revealed that the Rep protein of the CRESS-DNA viruses is classified as CRESSV1 to CRESSV6 groups and the new Smacoviridae-related, CRESSV2-related, pCRESS-related, Circoviridae-related, and 1 to 4 outgroups, according to the Viral_Rep and P-loop_NTPase domain organization, CLANS, and phylogenetic analysis. Furthermore, for the first time in this study, the Cap protein of CRESS-DNA viruses was classified into 20 distinct clusters by CLANS and phylogenetic analysis. Through this classification, the genetic diversity of CRESS-DNA viruses clarifies the possibility of recombinations in Cap and Rep proteins. Finally, it has been shown that selection pressure plays a significant role in the evolution and genetic diversity of Cap and Rep proteins. This study explains the genetic diversity of CRESS-DNA viruses and hopes that it will help classify future detected viruses.

Phage family classification under Caudoviricetes: A review of current tools using the latest ICTV classification framework

Bacteriophages, which are viruses infecting bacteria, are the most ubiquitous and diverse entities in the biosphere. There is accumulating evidence revealing their important roles in shaping the structure of various microbiomes. Thanks to (viral) metagenomic sequencing, a large number of new bacteriophages have been discovered. However, lacking a standard and automatic virus classification pipeline, the taxonomic characterization of new viruses seriously lag behind the sequencing efforts. In particular, according to the latest version of ICTV, several large phage families in the previous classification system are removed. Therefore, a comprehensive review and comparison of taxonomic classification tools under the new standard are needed to establish the state-of-the-art. In this work, we retrained and tested four recently published tools on newly labeled databases. We demonstrated their utilities and tested them on multiple datasets, including the RefSeq, short contigs, simulated metagenomic datasets, and low-similarity datasets. This study provides a comprehensive review of phage family classification in different scenarios and a practical guidance for choosing appropriate taxonomic classification pipelines. To our best knowledge, this is the first review conducted under the new ICTV classification framework. The results show that the new family classification framework overall leads to better conserved groups and thus makes family-level classification more feasible.

Plasma Virome Reveals Blooms and Transmission of Anellovirus in Intravenous Drug Users with HIV-1, HCV, and/or HBV Infections

Intravenous drug users (IDUs) are a high-risk group for HIV-1, hepatitis C virus (HCV), and hepatitis B virus (HBV) infections, which are the leading causes of death in IDUs. However, the plasma virome of IDUs and how it is influenced by above viral infections remain unclear. Using viral metagenomics, we determined the plasma virome of IDUs and its association with HIV-1, HCV, and/or HBV infections. Compared with healthy individuals, IDUs especially those with major viral infections had higher viral abundance and diversity. Anelloviridae dominated plasma virome. Coinfections of multiple anelloviruses were common, and anelloviruses from the same genus tended to coexist together. In this study, 4,487 anellovirus ORF1 sequences were identified, including 1,620 (36.1%) with less than 69% identity to any known sequences, which tripled the current number. Compared with healthy controls (HC), more anellovirus sequences were observed in neg-IDUs, and HIV-1, HCV, and/or HBV infections further expanded the sequence number in IDUs, which was characterized by the emergence of novel divergent taxons and blooms of resident anelloviruses. Pegivirus was mainly identified in infected IDUs. Five main pegivirus transmission clusters (TCs) were identified by phylogenetic analysis, suggesting a transmission link. Similar anellovirus profiles were observed in IDUs within the same TC, suggesting transmission of anellome among IDUs. Our data suggested that IDUs suffered higher plasma viral burden especially anelloviruses, which was associated with HIV-1, HCV, and/or HBV infections. Blooms in abundance and unprecedented diversity of anellovirus highlighted active evolution and replication of this virus in blood circulation, and an uncharacterized role it may engage with the host.

IMPORTANCE Virome is associated with immune status and determines or influences disease progression through both pathogenic and resident viruses. Increased viral burden in IDUs especially those with major viral infections indicated the suboptimal immune status and high infection risks of these population. Blooms in abundance and unprecedented diversity of anellovirus highlighted its active evolution and replication in the blood circulation, and sensitive response to other viral infections. In addition, transmission cluster analysis revealed the transmission link of pegivirus among IDUs, and the individuals with transmission links shared similar anellome profiles. In-depth monitoring of the plasma virome in high-risk populations is not only needed for surveillance for emerging viruses and transmission networks of major and neglected bloodborne viruses, but also important for a better understanding of commensal viruses and their role it may engage with immune system.

Genomic and Biological Profile of a Novel Bacteriophage, Vibrio phage Virtus, Which Improves Survival of Sparus aurata Larvae Challenged with Vibrio harveyi

Due to the emergence of multidrug-resistant bacteria, commonly known as “superbugs”, phage therapy for the control of bacterial diseases rose in popularity. In this context, the use of phages for the management of many important bacterial diseases in the aquaculture environment is auspicious. Vibrio harveyi, a well-known and serious bacterial pathogen, is responsible for many disease outbreaks in aquaculture, resulting in huge economic and production losses. We isolated and fully characterized a novel bacteriophage, Vibrio phage Virtus, infecting V. harveyi strain VH2. Vibrio phage Virtus can infect a wide spectrum of Vibrio spp., including strains of V. harveyi, V. owensii, V. campbellii, V. parahaemolyticus, and V. mediterranei. It has a latent period of 40 min with an unusually high burst size of 3200 PFU/cell. Vibrio phage Virtus has a double-stranded DNA of 82,960 base pairs with 127 predicted open reading frames (ORFs). No virulence, antibiotic resistance, or integrase-encoding genes were detected. In vivo phage therapy trials in gilthead seabream, Sparus aurata, larvae demonstrated that Vibrio phage Virtus was able to significantly improve the survival of larvae for five days at a multiplicity of infection (MOI) of 10, which suggests that it can be an excellent candidate for phage therapy.

Cellular homologs of the double jelly-roll major capsid proteins clarify the origins of an ancient virus kingdom

Viruses are the most abundant biological entities on Earth and ubiquitous parasites of cellular life forms. The general scenario for the origin of viruses involves evolution from nonviral replicators, such as plasmids and transposons, via recruitment of host proteins for virion formation. One of the most common virion core components, the double jelly-roll major capsid protein of a broad variety of viruses with double-stranded DNA genomes, so far has been thought to represent a virus innovation. However, we present evidence, obtained by protein structure comparison, that this type of virus capsid protein also evolved from a cellular ancestor, a distinct family of carbohydrate-active enzymes. These findings reinforce the chimeric scenario of virus origin.

Viruses are a distinct type of replicators that encode structural proteins encasing virus genomes in virions. For some of the widespread virus capsid proteins and other major components of virions, likely ancestors encoded by cellular life forms are identifiable. In particular, one of the most common capsid proteins, with the single jelly-roll (SJR) fold, appears to have evolved from a particular family of cellular carbohydrate-binding proteins. However, the double jelly-roll major capsid protein (DJR-MCP), the hallmark of the enormously diverse viruses of the kingdom Bamfordvirae within the realm Varidnaviria, which includes bacterial and archaeal icosahedral viruses as well as eukaryotic giant viruses, has been perceived as a virus innovation that evolved by duplication and fusion of the SJR capsid proteins. Here we employ protein structure comparison to show that the DJR fold is represented in several widespread families of cellular proteins, including several groups of carbohydrate-active enzymes. We show that DJR-MCPs share a common ancestry with a distinct family of bacterial DJR proteins (DUF2961) involved in carbohydrate metabolism. Based on this finding, we propose a scenario in which bamfordviruses evolved from nonviral replicators, in particular plasmids, by recruiting a host protein for capsid formation. This sequence of events appears to be the general route of virus origin. The results of this work indicate that virus kingdoms Bamfordvirae, with the DJR-MCPs, and Helvetiavirae that possess two SJR-MCPs, have distinct origins, suggesting a reappraisal of the realm Varidnaviria.

Communities of T4-like bacteriophages associated with bacteria in Lake Baikal: diversity and biogeography

Lake Baikal phage communities are important for lake ecosystem functioning. Here we describe the diversity of T4-bacteriophage associated with the bacterial fraction of filtered water samples collected from the pelagic zone, coastal zone and shallow bays. Although the study of the diversity of phages for the g23 gene has been carried out at Lake Baikal for more than ten years, shallow bays that comprise a significant part of the lake’s area have been neglected, and this gene has not previously been studied in the bacterial fraction. Phage communities were probed using amplicon sequencing methods targeting the gene of major capsid protein (g23) and compared phylogenetically across sample locations and with sequences previously retrieved from non-bacterial fractions (<0.2 um) and biofilms (non-fractionated). In this study, we examined six water samples, in which 24 to 74 viral OTUs were obtained. The sequences from shallow bays largely differed from those in the pelagic and coastal samples and formed individual subcluster in the UPGMA tree that was obtained from the comparison of phylogenetic distances of g23 sequence sets from various ecosystems, reflecting differences in viral communities depending on the productivity of various sites of Lake Baikal. According to the RefSeq database, from 58.3 to 73% of sequences of each sample had cultivated closest relatives belonging to cyanophages. In this study, for phylogenetic analysis, we chose the closest relatives not only from the RefSeq and GenBank NR databases but also from two marine and one freshwater viromes: eutrophic Osaka Bay (Japan), oligotrophic area of the Pacific Ocean (Station ALOHA) and mesotrophic and ancient Lake Biwa (Japan), which allowed us to more fully compare the diversity of marine and freshwater phages. The identity with marine sequences at the amino acid level ranged from 35 to 80%, and with the sequences from the viral fraction and bacterial one from Lake Biwa—from 35.3 to 98% and from 33.9 to 89.1%, respectively. Therefore, the sequences from marine viromes had a greater difference than those from freshwater viromes, which may indicate a close relationship between freshwater viruses and differences from marine viruses.

Saline lakes on the Qinghai-Tibet Plateau harbor unique viral assemblages mediating microbial environmental adaption

The highest plateau on Earth, Qinghai-Tibet Plateau, contains thousands of lakes with broad salinity and diverse and unique microbial communities. However, little is known about their co-occurring viruses. Herein, we identify 4,560 viral Operational Taxonomic Units (vOTUs) from six viromes of three saline lakes on Qinghai-Tibet Plateau, with less than 1% that could be classified. Most of the predicted vOTUs were associated with the dominant bacterial and archaeal phyla. Virus-encoded auxiliary metabolic genes suggest that viruses influence microbial metabolisms of carbon, nitrogen, sulfur, and lipid; the antibiotic resistance mediation; and their salinity adaption. The six viromes clustered together with the ice core viromes and bathypelagic ocean viromes and might represent a new viral habitat. This study has revealed the unique characteristics and potential ecological roles of DNA viromes in the lakes of the highest plateau and established a foundation for the recognition of the viral roles in plateau lake ecosystems.

Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes

The archaeal tailed viruses (arTV), evolutionarily related to tailed double-stranded DNA (dsDNA) bacteriophages of the class Caudoviricetes, represent the most common isolates infecting halophilic archaea. Only a handful of these viruses have been genomically characterized, limiting our appreciation of their ecological impacts and evolution. Here, we present 37 new genomes of haloarchaeal tailed virus isolates, more than doubling the current number of sequenced arTVs. Analysis of all 63 available complete genomes of arTVs, which we propose to classify into 14 new families and 3 orders, suggests ancient divergence of archaeal and bacterial tailed viruses and points to an extensive sharing of genes involved in DNA metabolism and counterdefense mechanisms, illuminating common strategies of virus-host interactions with tailed bacteriophages. Coupling of the comparative genomics with the host range analysis on a broad panel of haloarchaeal species uncovered 4 distinct groups of viral tail fiber adhesins controlling the host range expansion. The survey of metagenomes using viral hallmark genes suggests that the global architecture of the arTV community is shaped through recurrent transfers between different biomes, including hypersaline, marine, and anoxic environments.

Virus diversity in metagenomes of a lichen symbiosis (Umbilicaria phaea): complete viral genomes, putative hosts and elevational distributions

Viruses can play critical roles in symbioses by initiating horizontal gene transfer, affecting host phenotypes, or expanding their host's ecological niche. However, knowledge of viral diversity and distribution in symbiotic organisms remains elusive. Here we use deep-sequenced metagenomic DNA (PacBio Sequel II; two individuals), paired with a population genomics approach (Pool-seq; 11 populations, 550 individuals) to understand viral distributions in the lichen Umbilicaria phaea. We assess (i) viral diversity in lichen thalli, (ii) putative viral hosts (fungi, algae, bacteria) and (iii) viral distributions along two replicated elevation gradients. We identified five novel viruses, showing 28%-40% amino acid identity to known viruses. They tentatively belong to the families Caulimoviridae, Myoviridae, Podoviridae and Siphoviridae. Our analysis suggests that the Caulimovirus is associated with green algal photobionts (Trebouxia) of the lichen, and the remaining viruses with bacterial hosts. We did not detect viral sequences in the mycobiont. Caulimovirus abundance decreased with increasing elevation, a pattern reflected by a specific algal lineage hosting this virus. Bacteriophages showed population-specific patterns. Our work provides the first comprehensive insights into viruses associated with a lichen holobiont and suggests an interplay of viral hosts and environment in structuring viral distributions.

Spatial patterns in phage-Rhizobium coevolutionary interactions across regions of common bean domestication

Bacteriophages play significant roles in the composition, diversity, and evolution of bacterial communities. Despite their importance, it remains unclear how phage diversity and phage-host interactions are spatially structured. Local adaptation may play a key role. Nitrogen-fixing symbiotic bacteria, known as rhizobia, have been shown to locally adapt to domesticated common bean at its Mesoamerican and Andean sites of origin. This may affect phage-rhizobium interactions. However, knowledge about the diversity and coevolution of phages with their respective Rhizobium populations is lacking. Here, through the study of four phage-Rhizobium communities in Mexico and Argentina, we show that both phage and host diversity is spatially structured. Cross-infection experiments demonstrated that phage infection rates were higher overall in sympatric rhizobia than in allopatric rhizobia except for one Argentinean community, indicating phage local adaptation and host maladaptation. Phage-host interactions were shaped by the genetic identity and geographic origin of both the phage and the host. The phages ranged from specialists to generalists, revealing a nested network of interactions. Our results suggest a key role of local adaptation to resident host bacterial communities in shaping the phage genetic and phenotypic composition, following a similar spatial pattern of diversity and coevolution to that in the host.

Diverse Viruses in Deep-Sea Hydrothermal Vent Fluids Have Restricted Dispersal across Ocean Basins

In the ocean, viruses impact microbial mortality, regulate biogeochemical cycling, and alter the metabolic potential of microbial lineages. At deep-sea hydrothermal vents, abundant viruses infect a wide range of hosts among the archaea and bacteria that inhabit these dynamic habitats. However, little is known about viral diversity, host range, and biogeography across different vent ecosystems, which has important implications for how viruses manipulate microbial function and evolution. Here, we examined viral diversity, viral and host distribution, and virus-host interactions in microbial metagenomes generated from venting fluids from several vent sites within three different geochemically and geographically distinct hydrothermal systems: Piccard and Von Damm vent fields at the Mid-Cayman Rise in the Caribbean Sea, and at several vent sites within Axial Seamount in the Pacific Ocean. Analysis of viral sequences and clustered regularly interspaced short palindromic repeat (CRISPR) spacers revealed highly diverse viral assemblages and evidence of active infection. Network analysis revealed that viral host range was relatively narrow, with very few viruses infecting multiple microbial lineages. Viruses were largely endemic to individual vent sites, indicating restricted dispersal, and in some cases, viral assemblages persisted over time. Thus, we show that hydrothermal vent fluids are home to novel, diverse viral assemblages that are highly localized to specific regions and taxa. IMPORTANCE Viruses play important roles in manipulating microbial communities and their evolution in the ocean, yet not much is known about viruses in deep-sea hydrothermal vents. However, viral ecology and evolution are of particular interest in hydrothermal vent habitats because of their unique nature: previous studies have indicated that most viruses in hydrothermal vents are temperate rather than lytic, and it has been established that rates of horizontal gene transfer (HGT) are particularly high among thermophilic vent microbes, and viruses are common vectors for HGT. If viruses have broad host range or are widespread across vent sites, they have increased potential to act as gene-sharing "highways" between vent sites. By examining viral diversity, distribution, and infection networks across disparate vent sites, this study provides the opportunity to better characterize and constrain the viral impact on hydrothermal vent microbial communities. We show that viruses in hydrothermal vents are diverse and apparently active, but most have restricted host range and are not widely distributed among vent sites. Thus, the impacts of viral infection are likely to be highly localized and constrained to specific taxa in these habitats.

Reads Binning Improves the Assembly of Viral Genome Sequences From Metagenomic Samples

Metagenomes can be considered as mixtures of viral, bacterial, and other eukaryotic DNA sequences. Mining viral sequences from metagenomes could shed insight into virus-host relationships and expand viral databases. Current alignment-based methods are unsuitable for identifying viral sequences from metagenome sequences because most assembled metagenomic contigs are short and possess few or no predicted genes, and most metagenomic viral genes are dissimilar to known viral genes. In this study, I developed a Markov model-based method, VirMC, to identify viral sequences from metagenomic data. VirMC uses Markov chains to model sequence signatures and construct a scoring model using a likelihood test to distinguish viral and bacterial sequences. Compared with the other two state-of-the-art viral sequence-prediction methods, VirFinder and PPR-Meta, my proposed method outperformed VirFinder and had similar performance with PPR-Meta for short contigs with length less than 400 bp. VirMC outperformed VirFinder and PPR-Meta for identifying viral sequences in contaminated metagenomic samples with eukaryotic sequences. VirMC showed better performance in assembling viral-genome sequences from metagenomic data (based on filtering potential bacterial reads). Applying VirMC to human gut metagenomes from healthy subjects and patients with type-2 diabetes (T2D) revealed that viral contigs could help classify healthy and diseased statuses. This alignment-free method complements gene-based alignment approaches and will significantly improve the precision of viral sequence identification.

Deep viral blood metagenomics reveals extensive anellovirus diversity in healthy humans

Human blood metagenomics has revealed the presence of different types of viruses in apparently healthy subjects. By far, anelloviruses constitute the viral family that is more frequently found in human blood, although amplification biases and contaminations pose a major challenge in this field. To investigate this further, we subjected pooled plasma samples from 120 healthy donors in Spain to high-speed centrifugation, RNA and DNA extraction, random amplification, and massive parallel sequencing. Our results confirm the extensive presence of anelloviruses in such samples, which represented nearly 97% of the total viral sequence reads obtained. We assembled 114 different viral genomes belonging to this family, revealing remarkable diversity. Phylogenetic analysis of ORF1 suggested 28 potentially novel anellovirus species, 24 of which were validated by Sanger sequencing to discard artifacts. These findings underscore the importance of implementing more efficient purification procedures that enrich the viral fraction as an essential step in virome studies and question the suggested pathological role of anelloviruses.

The Unexplored Virome of Two Atlantic Coast Fish: Contribution of Next-Generation Sequencing to Fish Virology

Much of the knowledge on viruses is focused on those that can be propagated using cell-cultures or that can cause disease in humans or in economically important animals and plants. However, this only reflects a small portion of the virosphere. Therefore, in this study, we explore by targeted next-generation sequencing, how the virome varies between Atlantic horse mackerels and gilthead seabreams from fisheries and aquaculture from the center and south regions of Portugal. Viral genomes potentially pathogenic to fish and crustaceans, as well as to humans, were identified namelyese included Astroviridae, Nodaviridae, Hepadnaviridae, Birnaviridae, Caliciviridae, and Picornaviridae families. Also bacteriophages sequences were identified corresponding to the majority of sequencese detected, with Myoviridae, Podoviridae, and Siphoviridae, the most widespread families in both fish species. However, these findings can also be due to the presence of bacteria in fish tissues, or even to contamination. Overall, seabreams harbored viruses from a smaller number of families in comparison with mackerels. Therefore, the obtained data show that fish sold for consumption can harbor a high diversity of viruses, many of which are unknown, reflecting the overall uncharacterized virome of fish. While cross-species transmission of bonafide fish viruses to humans is unlikely, the finding of human pathogenic viruses in fish suggest that fish virome can be a potential threat regarding food safety.

Making sense of virus size and the tradeoffs shaping viral fitness

Viruses span an impressive size range, with genome length varying a thousandfold and virion volume nearly a millionfold. For cellular organisms the scaling of traits with size is a pervasive influence on ecological processes, but whether size plays a central role in viral ecology is unknown. Here, we focus on viruses of aquatic unicellular organisms, which exhibit the greatest known range of virus size. We outline hypotheses within a quantitative framework, and analyse data where available, to consider how size affects the primary components of viral fitness. We argue that larger viruses have fewer offspring per infection and slower contact rates with host cells, but a larger genome tends to increase infection efficiency, broaden host range, and potentially increase attachment success and decrease decay rate. These countervailing selective pressures may explain why a breadth of sizes exist and even coexist when infecting the same host populations. Oligotrophic ecosystems may be enriched in "giant" viruses, because environments with resource-limited phagotrophs at low concentrations may select for broader host range, better control of host metabolism, lower decay rate and a physical size that mimics bacterial prey. Finally, we describe where further research is needed to understand the ecology and evolution of viral size diversity.

T4-like myovirus community shaped by dispersal and deterministic processes in the South China Sea

As the most abundant and genetically diverse biological entities, viruses significantly influence ecological, biogeographical and evolutionary processes in the ocean. However, the biogeography of marine viruses and the drivers shaping viral community are unclear. Here, the biogeographic patterns of T4-like viruses and the relative impacts of deterministic (environmental selection) and dispersal (spatial distance) processes were investigated in the northern South China Sea. The dominant viral operational taxonomic units were affiliated with previously defined Marine, Estuary, Lake and Paddy Groups. A clear viral biogeographic pattern was observed along the environmental gradient from the estuary to open sea. Marine Groups I and IV had a wide geographical distribution, whereas Marine Groups II, III and V were abundant in lower-salinity continental or eutrophic environments. A significant distance-decay pattern was noted for the T4-like viral community, especially for those infecting cyanobacteria. Both deterministic and dispersal processes influenced viral community assembly, although environmental selection (e.g. temperature, salinity, bacterial abundance and community, etc.) had a greater impact than spatial distance. Network analysis confirmed the strong association between viral and bacterial community composition, and suggested a diverse ecological relationship (e.g. lysis, co-infection or mutualistic) between and within viruses and their potential bacterial hosts.

Lysogeny in the oceans: Lessons from cultivated model systems and a reanalysis of its prevalence

In the oceans, viruses that infect bacteria (phages) influence a variety of microbially mediated processes that drive global biogeochemical cycles. The nature of their influence is dependent upon infection mode, be it lytic or lysogenic. Temperate phages are predicted to be prevalent in marine systems where they are expected to execute both types of infection modes. Understanding the range and outcomes of temperate phage-host interactions is fundamental for evaluating their ecological impact. Here, we (i) review phage-mediated rewiring of host metabolism, with a focus on marine systems, (ii) consider the range and nature of temperate phage-host interactions, and (iii) draw on studies of cultivated model systems to examine the consequences of lysogeny among several dominant marine bacterial lineages. We also readdress the prevalence of lysogeny among marine bacteria by probing a collection of 1239 publicly available bacterial genomes, representing cultured and uncultivated strains, for evidence of complete prophages. Our conservative analysis, anticipated to underestimate true prevalence, predicts 18% of the genomes examined contain at least one prophage, the majority (97%) were found within genomes of cultured isolates. These results highlight the need for cultivation of additional model systems to better capture the diversity of temperate phage-host interactions in the oceans.

Unveiling Crucivirus Diversity by Mining Metagenomic Data

Viruses are the most abundant biological entities on Earth. In addition to their impact on animal and plant health, viruses have important roles in ecosystem dynamics as well as in the evolution of the biosphere. Circular Rep-encoding single-stranded (CRESS) DNA viruses are ubiquitous in nature, many are agriculturally important, and they appear to have multiple origins from prokaryotic plasmids. A subset of CRESS-DNA viruses, the cruciviruses, have homologues of capsid proteins encoded by RNA viruses. The genetic structure of cruciviruses attests to the transfer of capsid genes between disparate groups of viruses. However, the evolutionary history of cruciviruses is still unclear. By collecting and analyzing cruciviral sequence data, we provide a deeper insight into the evolutionary intricacies of cruciviruses. Our results reveal an unexpected diversity of this virus group, with frequent recombination as an important determinant of variability.

The discovery of cruciviruses revealed the most explicit example of a common protein homologue between DNA and RNA viruses to date. Cruciviruses are a novel group of circular Rep-encoding single-stranded DNA (ssDNA) (CRESS-DNA) viruses that encode capsid proteins that are most closely related to those encoded by RNA viruses in the family Tombusviridae. The apparent chimeric nature of the two core proteins encoded by crucivirus genomes suggests horizontal gene transfer of capsid genes between DNA and RNA viruses. Here, we identified and characterized 451 new crucivirus genomes and 10 capsid-encoding circular genetic elements through de novo assembly and mining of metagenomic data. These genomes are highly diverse, as demonstrated by sequence comparisons and phylogenetic analysis of subsets of the protein sequences they encode. Most of the variation is reflected in the replication-associated protein (Rep) sequences, and much of the sequence diversity appears to be due to recombination. Our results suggest that recombination tends to occur more frequently among groups of cruciviruses with relatively similar capsid proteins and that the exchange of Rep protein domains between cruciviruses is rarer than intergenic recombination. Additionally, we suggest members of the stramenopiles/alveolates/Rhizaria supergroup as possible crucivirus hosts. Altogether, we provide a comprehensive and descriptive characterization of cruciviruses.

Deep Roots and Splendid Boughs of the Global Plant Virome

Land plants host a vast and diverse virome that is dominated by RNA viruses, with major additional contributions from reverse-transcribing and single-stranded (ss) DNA viruses. Here, we introduce the recently adopted comprehensive taxonomy of viruses based on phylogenomic analyses, as applied to the plant virome. We further trace the evolutionary ancestry of distinct plant virus lineages to primordial genetic mobile elements. We discuss the growing evidence of the pivotal role of horizontal virus transfer from invertebrates to plants during the terrestrialization of these organisms, which was enabled by the evolution of close ecological associations between these diverse organisms. It is our hope that the emerging big picture of the formation and global architecture of the plant virome will be of broad interest to plant biologists and virologists alike and will stimulate ever deeper inquiry into the fascinating field of virus-plant coevolution.

Doubling of the known set of RNA viruses by metagenomic analysis of an aquatic virome

RNA viruses in aquatic environments remain poorly studied. Here, we analysed the RNA virome from approximately 10 l water from Yangshan Deep-Water Harbour near the Yangtze River estuary in China and identified more than 4,500 distinct RNA viruses, doubling the previously known set of viruses. Phylogenomic analysis identified several major lineages, roughly, at the taxonomic ranks of class, order and family. The 719-member-strong Yangshan virus assemblage is the sister clade to the expansive class Alsuviricetes and consists of viruses with simple genomes that typically encode only RNA-dependent RNA polymerase (RdRP), capping enzyme and capsid protein. Several clades within the Yangshan assemblage independently evolved domain permutation in the RdRP. Another previously unknown clade shares ancestry with Potyviridae, the largest known plant virus family. The ‘Aquatic picorna-like viruses/Marnaviridae’ clade was greatly expanded, with more than 800 added viruses. Several RdRP-linked protein domains not previously detected in any RNA viruses were identified, such as the small ubiquitin-like modifier (SUMO) domain, phospholipase A2 and PrsW-family protease domain. Multiple viruses utilize alternative genetic codes implying protist (especially ciliate) hosts. The results reveal a vast RNA virome that includes many previously unknown groups. However, phylogenetic analysis of the RdRPs supports the previously established five-branch structure of the RNA virus evolutionary tree, with no additional phyla.

Metagenomic analysis of a single RNA virome from the Yangshan Deep-Water Harbour in China enabled the recovery of more than 4,500 distinct RNA viruses, doubling the known set of RNA viruses to date, and provided insights into their biology.

Quantification of T4-Like and T7-Like Cyanophages Using the Polony Method Show They Are Significant Members of the Virioplankton in the North Pacific Subtropical Gyre

The North Pacific Subtropical Gyre (NPSG) is one of the largest biomes on Earth, with the cyanobacterium Prochlorococcus being the most abundant primary producer year-round. Viruses that infect cyanobacteria (cyanophages) influence cyanobacterial mortality, diversity and evolution. Two major cyanophage families are the T4-like cyanomyoviruses and T7-like cyanopodoviruses, yet their abundances and distribution patterns remain unknown due to difficulty in quantifying their populations. To address this limitation, we previously adapted the polony method (for PCR colony) to quantify T7-like cyanophages and applied it to spring populations in the Red Sea. Here, we further adapted the method for the quantification of T4-like cyanophages and analyzed the abundances of T4-like and T7-like cyanophage populations in the photic zone of the NPSG in summer 2015 and spring 2016. Combined, the peak abundances of these two cyanophage families reached 2.8 × 10⁶ and 1.1 × 10⁶ cyanophages ⋅ ml^–1 in the summer and spring, respectively. They constituted between 3 and 16% of total virus-like particles (VLPs), comprising a substantial component of the virioplankton in the NPSG. While both cyanophage families were highly abundant, the T4-like cyanophages were generally 1.3–4.4 fold more so. In summer, cyanophages had similar and reproducible distribution patterns with depth. Abundances were relatively low in the upper mixed layer and increased to form a pronounced subsurface peak at 100 m (1.9 × 10⁶ and 9.1 × 10⁵ phages ⋅ ml^–1 for the T4-like and T7-like cyanophages, respectively), coincident with the maximum in Prochlorococcus populations. Less vertical structure in cyanophage abundances was apparent in the spring profile, despite a subsurface peak in Prochlorococcus numbers. In the summer upper mixed layer, cyanophages constituted a smaller proportion of VLPs than below it and cyanophage to cyanobacteria ratios were considerably lower (1.3–2.8) than those of VLPs to bacteria (8.1–21.2). Differences in abundances between the two families and their contribution to VLPs with depth suggest differences in cyanophage production and/or decay processes relative to other members of the virioplankton in the upper mixed layer. These findings highlight the importance of quantifying distinct populations within the virioplankton to gain accurate understanding of their distribution patterns.

Global Organization and Proposed Megataxonomy of the Virus World

Viruses and mobile genetic elements are molecular parasites or symbionts that coevolve with nearly all forms of cellular life. The route of virus replication and protein expression is determined by the viral genome type. Comparison of these routes led to the classification of viruses into seven "Baltimore classes" (BCs) that define the major features of virus reproduction. However, recent phylogenomic studies identified multiple evolutionary connections among viruses within each of the BCs as well as between different classes. Due to the modular organization of virus genomes, these relationships defy simple representation as lines of descent but rather form complex networks. Phylogenetic analyses of virus hallmark genes combined with analyses of gene-sharing networks show that replication modules of five BCs (three classes of RNA viruses and two classes of reverse-transcribing viruses) evolved from a common ancestor that encoded an RNA-directed RNA polymerase or a reverse transcriptase. Bona fide viruses evolved from this ancestor on multiple, independent occasions via the recruitment of distinct cellular proteins as capsid subunits and other structural components of virions. The single-stranded DNA (ssDNA) viruses are a polyphyletic class, with different groups evolving by recombination between rolling-circle-replicating plasmids, which contributed the replication protein, and positive-sense RNA viruses, which contributed the capsid protein. The double-stranded DNA (dsDNA) viruses are distributed among several large monophyletic groups and arose via the combination of distinct structural modules with equally diverse replication modules. Phylogenomic analyses reveal the finer structure of evolutionary connections among RNA viruses and reverse-transcribing viruses, ssDNA viruses, and large subsets of dsDNA viruses. Taken together, these analyses allow us to outline the global organization of the virus world. Here, we describe the key aspects of this organization and propose a comprehensive hierarchical taxonomy of viruses.

Extreme Viral Partitioning in a Marine-Derived High Arctic Lake

The Arctic is warming at an accelerating pace, and the rise in temperature has increasing impacts on the Arctic biome. Lakes are integrators of their surroundings and thus excellent sentinels of environmental change. Despite their importance in the regulation of key microbial processes, viruses remain largely uncharacterized in Arctic lacustrine environments. We sampled a highly stratified meromictic lake near the northern limit of the Canadian High Arctic, a region in rapid transition due to climate change. We found that the different layers of the lake harbored viral communities that were strikingly dissimilar and highly divergent from known viruses. Viruses were more abundant in the deepest part of the lake containing ancient Arctic Ocean seawater that was trapped during glacial retreat and were genomically unlike any viruses previously described. This research demonstrates the complexity and novelty of viral communities in an environment that is vulnerable to ongoing perturbation.

High-latitude, perennially stratified (meromictic) lakes are likely to be especially vulnerable to climate warming because of the importance of ice in maintaining their water column structure and associated distribution of microbial communities. This study aimed to characterize viral abundance, diversity, and distribution in a meromictic lake of marine origin on the far northern coast of Ellesmere Island, in the Canadian High Arctic. We collected triplicate samples for double-stranded DNA (dsDNA) viromics from five depths that encompassed the major features of the lake, as determined by limnological profiling of the water column. Viral abundance and virus-to-prokaryote ratios were highest at greater depths, while bacterial and cyanobacterial counts were greatest in the surface waters. The viral communities from each zone of the lake defined by salinity, temperature, and dissolved oxygen concentrations were markedly distinct, suggesting that there was little exchange of viral types among lake strata. Ten viral assembled genomes were obtained from our libraries, and these also segregated with depth. This well-defined structure of viral communities was consistent with that of potential hosts. Viruses from the monimolimnion, a deep layer of ancient Arctic Ocean seawater, were more diverse and relatively abundant, with few similarities to available viral sequences. The Lake A viral communities also differed from published records from the Arctic Ocean and meromictic Ace Lake in Antarctica. This first characterization of viral diversity from this sentinel environment underscores the microbial richness and complexity of an ecosystem type that is increasingly exposed to major perturbations in the fast-changing Arctic.

IMPORTANCE The Arctic is warming at an accelerating pace, and the rise in temperature has increasing impacts on the Arctic biome. Lakes are integrators of their surroundings and thus excellent sentinels of environmental change. Despite their importance in the regulation of key microbial processes, viruses remain largely uncharacterized in Arctic lacustrine environments. We sampled a highly stratified meromictic lake near the northern limit of the Canadian High Arctic, a region in rapid transition due to climate change. We found that the different layers of the lake harbored viral communities that were strikingly dissimilar and highly divergent from known viruses. Viruses were more abundant in the deepest part of the lake containing ancient Arctic Ocean seawater that was trapped during glacial retreat and were genomically unlike any viruses previously described. This research demonstrates the complexity and novelty of viral communities in an environment that is vulnerable to ongoing perturbation.

Origin of viruses: primordial replicators recruiting capsids from hosts

Viruses are ubiquitous parasites of cellular life and the most abundant biological entities on Earth. It is widely accepted that viruses are polyphyletic, but a consensus scenario for their ultimate origin is still lacking. Traditionally, three scenarios for the origin of viruses have been considered: descent from primordial, precellular genetic elements, reductive evolution from cellular ancestors and escape of genes from cellular hosts, achieving partial replicative autonomy and becoming parasitic genetic elements. These classical scenarios give different timelines for the origin(s) of viruses and do not explain the provenance of the two key functional modules that are responsible, respectively, for viral genome replication and virion morphogenesis. Here, we outline a 'chimeric' scenario under which different types of primordial, selfish replicons gave rise to viruses by recruiting host proteins for virion formation. We also propose that new groups of viruses have repeatedly emerged at all stages of the evolution of life, often through the displacement of ancestral structural and genome replication genes.

High-throughput sequencing (HTS) for the analysis of viral populations

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

Gene Transfer Agents in Symbiotic Microbes

Prokaryotes commonly undergo genome reduction, particularly in the case of symbiotic bacteria. Genome reductions tend toward the energetically favorable removal of unnecessary, redundant, or nonfunctional genes. However, without mechanisms to compensate for these losses, deleterious mutation and genetic drift might otherwise overwhelm a population. Among the mechanisms employed to counter gene loss and share evolutionary success within a population, gene transfer agents (GTAs) are increasingly becoming recognized as important contributors. Although viral in origin, GTA particles package fragments of their "host" genome for distribution within a population of cells, often in a synchronized manner, rather than selfishly packaging genes necessary for their spread. Microbes as diverse as archaea and alpha-proteobacteria have been known to produce GTA particles, which are capable of transferring selective advantages such as virulence factors and antibiotic resistance. In this review, we discuss the various types of GTAs identified thus far, focusing on a defined set of symbiotic alpha-proteobacteria known to carry them. Drawing attention to the predicted presence of these genes, we discuss their potential within the selective marine and terrestrial environments occupied by mutualistic, parasitic, and endosymbiotic microbes.

An Ecosystems Perspective on Virus Evolution and Emergence

Understanding the emergence of pathogenic viruses has dominated studies of virus evolution. However, new metagenomic studies imply that relatively few of an immense number of viruses may lead to overt disease. This suggests a change in emphasis, from viruses as habitual pathogens to integral components of ecosystems. Here we show how viruses alter interactions between host individuals, populations, and ecosystems, impacting ecosystem health, resilience, and function, and how host ecology in turn impacts viral abundance and diversity. Moving to an ecosystems perspective will put virus evolution and disease emergence in its true context, and enhance our understanding of ecological processes.

Virus-borne mini-CRISPR arrays are involved in interviral conflicts

CRISPR-Cas immunity is at the forefront of antivirus defense in bacteria and archaea and specifically targets viruses carrying protospacers matching the spacers catalogued in the CRISPR arrays. Here, we perform deep sequencing of the CRISPRome-all spacers contained in a microbiome-associated with hyperthermophilic archaea of the order Sulfolobales recovered directly from an environmental sample and from enrichment cultures established in the laboratory. The 25 million CRISPR spacers sequenced from a single sampling site dwarf the diversity of spacers from all available Sulfolobales isolates and display complex temporal dynamics. Comparison of closely related virus strains shows that CRISPR targeting drives virus genome evolution. Furthermore, we show that some archaeal viruses carry mini-CRISPR arrays with 1-2 spacers and preceded by leader sequences but devoid of cas genes. Closely related viruses present in the same population carry spacers against each other. Targeting by these virus-borne spacers represents a distinct mechanism of heterotypic superinfection exclusion and appears to promote archaeal virus speciation.

Demographic and environmental drivers of metagenomic viral diversity in vampire bats

Viruses infect all forms of life and play critical roles as agents of disease, drivers of biochemical cycles and sources of genetic diversity for their hosts. Our understanding of viral diversity derives primarily from comparisons among host species, precluding insight into how intraspecific variation in host ecology affects viral communities or how predictable viral communities are across populations. Here we test spatial, demographic and environmental hypotheses explaining viral richness and community composition across populations of common vampire bats, which occur in diverse habitats of North, Central and South America. We demonstrate marked variation in viral communities that was not consistently predicted by a null model of declining community similarity with increasing spatial or genetic distances separating populations. We also find no evidence that larger bat colonies host greater viral diversity. Instead, viral diversity follows an elevational gradient, is enriched by juvenile-biased age structure, and declines with local anthropogenic food resources as measured by livestock density. Our results establish the value of linking the modern influx of metagenomic sequence data with comparative ecology, reveal that snapshot views of viral diversity are unlikely to be representative at the species level, and affirm existing ecological theories that link host ecology not only to single pathogen dynamics but also to viral communities.

Multiple origins of prokaryotic and eukaryotic single-stranded DNA viruses from bacterial and archaeal plasmids

Single-stranded (ss) DNA viruses are a major component of the earth virome. In particular, the circular, Rep-encoding ssDNA (CRESS-DNA) viruses show high diversity and abundance in various habitats. By combining sequence similarity network and phylogenetic analyses of the replication proteins (Rep) belonging to the HUH endonuclease superfamily, we show that the replication machinery of the CRESS-DNA viruses evolved, on three independent occasions, from the Reps of bacterial rolling circle-replicating plasmids. The CRESS-DNA viruses emerged via recombination between such plasmids and cDNA copies of capsid genes of eukaryotic positive-sense RNA viruses. Similarly, the rep genes of prokaryotic DNA viruses appear to have evolved from HUH endonuclease genes of various bacterial and archaeal plasmids. Our findings also suggest that eukaryotic polyomaviruses and papillomaviruses with dsDNA genomes have evolved via parvoviruses from CRESS-DNA viruses. Collectively, our results shed light on the complex evolutionary history of a major class of viruses revealing its polyphyletic origins.

Most single-stranded DNA viruses have small genomes replicated by rolling circle mechanism which is initiated by the Rep protein. Here, using sequence similarity network and phylogenetic analyses, Kazlauskas et al. show that viral Reps evolved from Reps of bacterial and archaeal plasmids on multiple independent occasions.

Marine RNA Virus Quasispecies Are Distributed throughout the Oceans

Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution.

RNA viruses, particularly genetically diverse members of the Picornavirales, are widespread and abundant in the ocean. Gene surveys suggest that there are spatial and temporal patterns in the composition of RNA virus assemblages, but data on their diversity and genetic variability in different oceanographic settings are limited. Here, we show that specific RNA virus genomes have widespread geographic distributions and that the dominant genotypes are under purifying selection. Genomes from three previously unknown picorna-like viruses (BC-1, -2, and -3) assembled from a coastal site in British Columbia, Canada, as well as marine RNA viruses JP-A, JP-B, and Heterosigma akashiwo RNA virus exhibited different biogeographical patterns. Thus, biotic factors such as host specificity and viral life cycle, and not just abiotic processes such as dispersal, affect marine RNA virus distribution. Sequence differences relative to reference genomes imply that virus quasispecies are under purifying selection, with synonymous single-nucleotide variations dominating in genomes from geographically distinct regions resulting in conservation of amino acid sequences. Conversely, sequences from coastal South Africa that mapped to marine RNA virus JP-A exhibited more nonsynonymous mutations, probably representing amino acid changes that accumulated over a longer separation. This biogeographical analysis of marine RNA viruses demonstrates that purifying selection is occurring across oceanographic provinces. These data add to the spectrum of known marine RNA virus genomes, show the importance of dispersal and purifying selection for these viruses, and indicate that closely related RNA viruses are pathogens of eukaryotic microbes across oceans.

IMPORTANCE Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution.

Classifying the Unclassified: A Phage Classification Method

This work reports the method ClassiPhage to classify phage genomes using sequence derived taxonomic features. ClassiPhage uses a set of phage specific Hidden Markov Models (HMMs) generated from clusters of related proteins. The method was validated on all publicly available genomes of phages that are known to infect Vibrionaceae. The phages belong to the well-described phage families of Myoviridae, Podoviridae, Siphoviridae, and Inoviridae. The achieved classification is consistent with the assignments of the International Committee on Taxonomy of Viruses (ICTV), all tested phages were assigned to the corresponding group of the ICTV-database. In addition, 44 out of 58 genomes of Vibrio phages not yet classified could be assigned to a phage family. The remaining 14 genomes may represent phages of new families or subfamilies. Comparative genomics indicates that the ability of the approach to identify and classify phages is correlated to the conserved genomic organization. ClassiPhage classifies phages exclusively based on genome sequence data and can be applied on distinct phage genomes as well as on prophage regions within host genomes. Possible applications include (a) classifying phages from assembled metagenomes; and (b) the identification and classification of integrated prophages and the splitting of phage families into subfamilies.

The Environment and Cyanophage Diversity: Insights From Environmental Sequencing of DNA Polymerase

Globally distributed and abundant cyanophages in the family Myoviridae have dsDNA genomes with variable gene content, including host-derived auxiliary metabolic genes (AMGs) that potentially can facilitate viral replication. However, it is not well understood how this variation in gene content interacts with environmental variables to shape cyanomyovirus communities. This project correlated the genetic repertoire of cyanomyoviruses with their phyologeny, and investigated cyanomyovirus ecotype distribution as a function of environmental conditions across locations and seasons. Reference cyanomyovirus genomes were compared for their overlap in gene content to infer phyologenetic distances, and these distances were compared to distances calculated based on DNA polymerase (gp43) gene sequences. In turn, gp43 partial gene sequences amplified from natural cyanophage communities were used to describe cyanomyovirus community composition and to assess the relationship between environmental variables. The results showed the following: (1) DNA polymerase gene phylogeny generally correlated with the similarity in gene content among reference cyanomyoviruses, and thus can be used to describe environmental cyanomyovirus communities; (2) spatial and seasonal patterns in cyanomyovirus communities were related to environmental variables; (3) salinity and temperature, combined with nutrient concentration were predictors of cyanomyovirus richness, diversity and community composition. This study shows that environmental variables shape viral communities by drawing on a diverse seed bank of viral genotypes. From these results it is evident that that viral ecotypes with their corresponding genetic repertoires underlie selection pressures. However, the mechanisms involved in selecting for specific viral genotypes remain to be fully understood.