Ecogenomics of virophages and their giant virus hosts assessed through time series metagenomics

Time-series metagenomics reveals changing protistan ecology of a temperate dimictic lake

BACKGROUND

Protists, single-celled eukaryotic organisms, are critical to food web ecology, contributing to primary productivity and connecting small bacteria and archaea to higher trophic levels. Lake Mendota is a large, eutrophic natural lake that is a Long-Term Ecological Research site and among the world's best-studied freshwater systems. Metagenomic samples have been collected and shotgun sequenced from Lake Mendota for the last 20 years. Here, we analyze this comprehensive time series to infer changes to the structure and function of the protistan community and to hypothesize about their interactions with bacteria.

RESULTS

Based on small subunit rRNA genes extracted from the metagenomes and metagenome-assembled genomes of microeukaryotes, we identify shifts in the eukaryotic phytoplankton community over time, which we predict to be a consequence of reduced zooplankton grazing pressures after the invasion of a invasive predator (the spiny water flea) to the lake. The metagenomic data also reveal the presence of the spiny water flea and the zebra mussel, a second invasive species to Lake Mendota, prior to their visual identification during routine monitoring. Furthermore, we use species co-occurrence and co-abundance analysis to connect the protistan community with bacterial taxa. Correlation analysis suggests that protists and bacteria may interact or respond similarly to environmental conditions. Cryptophytes declined in the second decade of the timeseries, while many alveolate groups (e.g., ciliates and dinoflagellates) and diatoms increased in abundance, changes that have implications for food web efficiency in Lake Mendota.

CONCLUSIONS

We demonstrate that metagenomic sequence-based community analysis can complement existing efforts to monitor protists in Lake Mendota based on microscopy-based count surveys. We observed patterns of seasonal abundance in microeukaryotes in Lake Mendota that corroborated expectations from other systems, including high abundance of cryptophytes in winter and diatoms in fall and spring, but with much higher resolution than previous surveys. Our study identified long-term changes in the abundance of eukaryotic microbes and provided context for the known establishment of an invasive species that catalyzes a trophic cascade involving protists. Our findings are important for decoding potential long-term consequences of human interventions, including invasive species introduction. Video Abstract.

Virophages, Satellite Viruses, Virophage Replication and Its Effects and Virophage Defence Mechanisms for Giant Virus Hosts and Giant Virus Defence Systems against Virophages

In this paper, the characteristics of 40 so far described virophages-parasites of giant viruses-are given, and the similarities and differences between virophages and satellite viruses, which also, like virophages, require helper viruses for replication, are described. The replication of virophages taking place at a specific site-the viral particle factory of giant viruses-and its consequences are presented, and the defence mechanisms of virophages for giant virus hosts, as a protective action for giant virus hosts-protozoa and algae-are approximated. The defence systems of giant viruses against virophages were also presented, which are similar to the CRISPR/Cas defence system found in bacteria and in Archea. These facts, and related to the very specific biological features of virophages (specific site of replication, specific mechanisms of their defensive effects for giant virus hosts, defence systems in giant viruses against virophages), indicate that virophages, and their host giant viruses, are biological objects, forming a 'novelty' in biology.

Widespread occurrence and diverse origins of polintoviruses influence lineage-specific genome dynamics in stony corals

Stony corals (Order: Scleractinia) are central to vital marine habitats known as coral reefs. Numerous stressors in the Anthropocene are contributing to the ongoing decline in coral reef health and coverage. While viruses are established modulators of marine microbial dynamics, their interactions within the coral holobiont and impact on coral health and physiology remain unclear. To address this key knowledge gap, we investigated diverse stony coral genomes for 'endogenous' viruses. Our study uncovered a remarkable number of integrated viral elements recognized as 'Polintoviruses' (Class Polintoviricetes) in thirty Scleractinia genomes; with several species harboring hundreds to thousands of polintoviruses. We reveal massive paralogous expansion of polintoviruses in stony coral genomes, alongside the presence of integrated elements closely related to Polinton-like viruses (PLVs), a group of viruses that exist as free virions. These results suggest multiple integrations of polintoviruses and PLV-relatives, along with paralogous expansions, shaped stony coral genomes. Re-analysis of existing gene expression data reveals all polintovirus structural and non-structural hallmark genes are expressed, providing support for free virion production from polintoviruses. Our results, revealing a significant diversity of polintovirus across the Scleractinia order, open a new research avenue into polintovirus and their possible roles in disease, genomic plasticity, and environmental adaptation in this key group of organisms.

Endogenous virophages are active and mitigate giant virus infection in the marine protist Cafeteria burkhardae

Endogenous viral elements (EVEs) are common genetic passengers in various protists. Some EVEs represent viral fossils, whereas others are still active. The marine heterotrophic flagellate Cafeteria burkhardae contains several EVE types related to the virophage mavirus, a small DNA virus that parasitizes the lytic giant virus CroV. We hypothesized that endogenous virophages may act as an antiviral defense system in protists, but no protective effect of virophages in wild host populations has been shown so far. Here, we tested the activity of virophage EVEs and studied their impact on giant virus replication. We found that endogenous mavirus-like elements (EMALEs) from globally distributed Cafeteria populations produced infectious virus particles specifically in response to CroV infection. However, reactivation was stochastic, often inefficient, and poorly reproducible. Interestingly, only one of eight EMALE types responded to CroV infection, implying that other EMALEs may be linked to different giant viruses. We isolated and cloned several reactivated virophages and characterized their particles, genomes, and infection dynamics. All tested virophages inhibited the production of CroV during coinfection, thereby preventing lysis of the host cultures in a dose-dependent manner. Comparative genomics of different C. burkhardae strains revealed that inducible EMALEs are common and are not linked to specific geographic locations. We demonstrate that naturally occurring virophage EVEs reactivate upon giant virus infection, thus providing a striking example that eukaryotic EVEs can become active under specific conditions. Moreover, our results support the hypothesis that virophages can act as an adaptive antiviral defense system in protists.

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

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

Community ecology and functional potential of bacteria, archaea, eukarya and viruses in Guerrero Negro microbial mat

In this study, the microbial ecology, potential environmental adaptive mechanisms, and the potential evolutionary interlinking of genes between bacterial, archaeal and viral lineages in Guerrero Negro (GN) microbial mat were investigated using metagenomic sequencing across a vertical transect at millimeter scale. The community composition based on unique genes comprised bacteria (98.01%), archaea (1.81%), eukarya (0.07%) and viruses (0.11%). A gene-focused analysis of bacteria archaea, eukarya and viruses showed a vertical partition of the community. The greatest coverages of genes of bacteria and eukarya were detected in first layers, while the highest coverages of genes of archaea and viruses were found in deeper layers. Many genes potentially related to adaptation to the local environment were detected, such as UV radiation, multidrug resistance, oxidative stress, heavy metals, salinity and desiccation. Those genes were found in bacterial, archaeal and viral lineages with 6477, 44, and 1 genes, respectively. The evolutionary histories of those genes were studied using phylogenetic analysis, showing an interlinking between domains in GN mat.

Virome of Australia's most endangered parrot in captivity evidenced of harboring hitherto unknown viruses

IMPORTANCE

The impact of circulating viruses on the critically endangered, orange-bellied parrot (OBP) population can be devastating. The OBP already faces numerous threats to its survival in the wild, including habitat loss, predation, and small population impacts. Conservation of the wild OBP population is heavily reliant on supplementation using OBPs from a managed captive breeding program. These birds may act as a source for introduction of a novel disease agent to the wild population that may affect survival and reproduction. It is, therefore, essential to monitor and assess the health of OBPs and take appropriate measures to prevent and control the spread of viral infections. This requires knowledge of the existing virome to identify novel and emerging viruses and support development of appropriate measures to manage associated risk. By monitoring and protecting these animals from emerging viral diseases, we can help ensure their ongoing survival and preserve the biodiversity of our planet.

Virophages Found in Viromes from Lake Baikal

In this study, a previously little-studied group of viruses-virophages-was searched for and identified in the viromes of the ancient oligotrophic Lake Baikal. Virophages are small dsDNA viruses that parasitize giant viruses (e.g., Mimiviridae), which in turn affect unicellular eukaryotes. We analyzed eight viromes obtained from the deep-water areas of three basins of Lake Baikal and the shallow-water strait Maloye More in different seasons. The sequences of virophages were revealed in all viromes and were dominant after bacteriophages and algal viruses. Sixteen putative complete genomes of virophages were assembled, all of which contained four conserved genes encoding major capsid protein (MCP), minor capsid protein (mCP), maturation cysteine protease (PRO), and FtsK-HerA family DNA-packaging ATPase (ATPase). The MCP-based cluster analysis showed a sequence separation according to seasons, and a dependence on the geographical localization was not detected.

Viruses in deep-sea cold seep sediments harbor diverse survival mechanisms and remain genetically conserved within species

Deep sea cold seep sediments have been discovered to harbor novel, abundant, and diverse bacterial and archaeal viruses. However, little is known about viral genetic features and evolutionary patterns in these environments. Here, we examined the evolutionary ecology of viruses across active and extinct seep stages in the area of Haima cold seeps in the South China Sea. A total of 338 viral operational taxonomic units are identified and linked to 36 bacterial and archaeal phyla. The dynamics of host-virus interactions are informed by diverse antiviral defense systems across 43 families found in 487 microbial genomes. Cold seep viruses are predicted to harbor diverse adaptive strategies to persist in this environment, including counter-defense systems, auxiliary metabolic genes, reverse transcriptases, and alternative genetic code assignments. Extremely low nucleotide diversity is observed in cold seep viral populations, being influenced by factors including microbial host, sediment depth, and cold seep stage. Most cold seep viral genes are under strong purifying selection with trajectories that differ depending on whether cold seeps are active or extinct. This work sheds light on the understanding of environmental adaptation mechanisms and evolutionary patterns of viruses in the sub-seafloor biosphere.

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

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

Identifying local associations in biological time series: algorithms, statistical significance, and applications

Local associations refer to spatial-temporal correlations that emerge from the biological realm, such as time-dependent gene co-expression or seasonal interactions between microbes. One can reveal the intricate dynamics and inherent interactions of biological systems by examining the biological time series data for these associations. To accomplish this goal, local similarity analysis algorithms and statistical methods that facilitate the local alignment of time series and assess the significance of the resulting alignments have been developed. Although these algorithms were initially devised for gene expression analysis from microarrays, they have been adapted and accelerated for multi-omics next generation sequencing datasets, achieving high scientific impact. In this review, we present an overview of the historical developments and recent advances for local similarity analysis algorithms, their statistical properties, and real applications in analyzing biological time series data. The benchmark data and analysis scripts used in this review are freely available at http://github.com/labxscut/lsareview.

Astrovirology: how viruses enhance our understanding of life in the Universe

Viruses are the most numerically abundant biological entities on Earth. As ubiquitous replicators of molecular information and agents of community change, viruses have potent effects on the life on Earth, and may play a critical role in human spaceflight, for life-detection missions to other planetary bodies and planetary protection. However, major knowledge gaps constrain our understanding of the Earth's virosphere: (1) the role viruses play in biogeochemical cycles, (2) the origin(s) of viruses and (3) the involvement of viruses in the evolution, distribution and persistence of life. As viruses are the only replicators that span all known types of nucleic acids, an expanded experimental and theoretical toolbox built for Earth's viruses will be pivotal for detecting and understanding life on Earth and beyond. Only by filling in these knowledge and technical gaps we will obtain an inclusive assessment of how to distinguish and detect life on other planetary surfaces. Meanwhile, space exploration requires life-support systems for the needs of humans, plants and their microbial inhabitants. Viral effects on microbes and plants are essential for Earth's biosphere and human health, but virus-host interactions in spaceflight are poorly understood. Viral relationships with their hosts respond to environmental changes in complex ways which are difficult to predict by extrapolating from Earth-based proxies. These relationships should be studied in space to fully understand how spaceflight will modulate viral impacts on human health and life-support systems, including microbiomes. In this review, we address key questions that must be examined to incorporate viruses into Earth system models, life-support systems and life detection. Tackling these questions will benefit our efforts to develop planetary protection protocols and further our understanding of viruses in astrobiology.

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

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

Virophages-Known and Unknown Facts

The paper presents virophages, which, like their host, giant viruses, are "new" infectious agents whose role in nature, including mammalian health, is important. Virophages, along with their protozoan and algal hosts, are found in fresh inland waters and oceanic and marine waters, including thermal waters and deep-sea vents, as well as in soil, plants, and in humans and animals (ruminants). Representing "superparasitism", almost all of the 39 described virophages (except Zamilon) interact negatively with giant viruses by affecting their replication and morphogenesis and their "adaptive immunity". This causes them to become regulators and, at the same time, defenders of the host of giant viruses protozoa and algae, which are organisms that determine the homeostasis of the aquatic environment. They are classified in the family Lavidaviridae with two genus (Sputnikovirus, Mavirus). However, in 2023, a proposal was presented that they should form the class Maveriviricetes, with four orders and seven families. Their specific structure, including their microsatellite (SSR-Simple Sequence Repeats) and the CVV (cell-virus-virophage, or transpovirion) system described with them, as well as their function, makes them, together with the biological features of giant viruses, form the basis for discussing the existence of a fourth domain in addition to Bacteria, Archaea, and Eukaryota. The paper also presents the hypothetical possibility of using them as a vector for vaccine antigens.

The virome of the last eukaryotic common ancestor and eukaryogenesis

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

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.

Genomics discovery of giant fungal viruses from subsurface oceanic crustal fluids

The oceanic igneous crust is a vast reservoir for microbial life, dominated by diverse and active bacteria, archaea, and fungi. Archaeal and bacterial viruses were previously detected in oceanic crustal fluids at the Juan de Fuca Ridge (JdFR). Here we report the discovery of two eukaryotic Nucleocytoviricota genomes from the same crustal fluids by sorting and sequencing single virions. Both genomes have a tRNA^Tyr gene with an intron (20 bps) at the canonical position between nucleotide 37 and 38, a common feature in eukaryotic and archaeal tRNA genes with short introns (<100 bps), and fungal genes acquired through horizontal gene transfer (HGT) events. The dominance of Ascomycota fungi as the main eukaryotes in crustal fluids and the evidence for HGT point to these fungi as the putative hosts, making these the first putative fungi-Nucleocytoviricota specific association. Our study suggests active host-viral dynamics for the only eukaryotic group found in the subsurface oceanic crust and raises important questions about the impact of viral infection on the productivity and biogeochemical cycling in this ecosystem.

Updated Virophage Taxonomy and Distinction from Polinton-like Viruses

Virophages are small dsDNA viruses that hijack the machinery of giant viruses during the co-infection of a protist (i.e., microeukaryotic) host and represent an exceptional case of "hyperparasitism" in the viral world. While only a handful of virophages have been isolated, a vast diversity of virophage-like sequences have been uncovered from diverse metagenomes. Their wide ecological distribution, idiosyncratic infection and replication strategy, ability to integrate into protist and giant virus genomes and potential role in antiviral defense have made virophages a topic of broad interest. However, one limitation for further studies is the lack of clarity regarding the nomenclature and taxonomy of this group of viruses. Specifically, virophages have been linked in the literature to other "virophage-like" mobile genetic elements and viruses, including polinton-like viruses (PLVs), but there are no formal demarcation criteria and proper nomenclature for either group, i.e., virophage or PLVs. Here, as part of the ICTV Virophage Study Group, we leverage a large set of genomes gathered from published datasets as well as newly generated protist genomes to propose delineation criteria and classification methods at multiple taxonomic ranks for virophages 'sensu stricto', i.e., genomes related to the prototype isolates Sputnik and mavirus. Based on a combination of comparative genomics and phylogenetic analyses, we show that this group of virophages forms a cohesive taxon that we propose to establish at the class level and suggest a subdivision into four orders and seven families with distinctive ecogenomic features. Finally, to illustrate how the proposed delineation criteria and classification method would be used, we apply these to two recently published datasets, which we show include both virophages and other virophage-related elements. Overall, we see this proposed classification as a necessary first step to provide a robust taxonomic framework in this area of the virosphere, which will need to be expanded in the future to cover other virophage-related viruses such as PLVs.

Unravelling Bile Viromes of Free-Range Laying Chickens Clinically Diagnosed with Spotty Liver Disease: Emergence of Many Novel Chaphamaparvoviruses into Multiple Lineages

Spotty liver disease (SLD) causes substantial egg production losses and chicken mortality; therefore, it is a disease that concerns Australian egg farmers. Over the last few decades, much research has been conducted to determine the etiologic agents of SLD and to develop potential therapeutics; however, SLD still remains a major issue for the chicken industries globally and remained without the elucidation of potentially multiple pathogens involved. To help fill this gap, this study was aimed at understanding the viral diversity of bile samples from which the SLD-causing bacterium, Campylobacter hepaticus, has been isolated and characterised. The collected samples were processed and sequenced using high-throughput next-generation sequencing. Remarkably, this study found 15 galliform chaphamaparvoviruses (GaChPVs), of which 14 are novel under the genus Chaphamaparvovirus. Among them, nine were complete genomes that showed between 41.7% and 78.3% genome-wide pairwise similarities to one another. Subsequent phylogenetic analysis using the NS1 gene exhibited a multiple incursion of chaphamaparvovirus lineages, including a novel lineage of unknown ancestral history in free-range laying chickens in Australia. This is the first evidence of circulating many parvoviruses in chickens in Australia, which has increased our knowledge of the pathogen diversity that may have an association with SLD in chickens.

Past and present giant viruses diversity explored through permafrost metagenomics

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

Although giant viruses are abundant in aquatic environments, less is known about giant viruses in soil. Here, the authors use permafrost metagenomics to reveal giant virus diversity and heterogeneity, as well as gene transfers between viruses from different families.

The Astounding World of Glycans from Giant Viruses

Viruses are a heterogeneous ensemble of entities, all sharing the need for a suitable host to replicate. They are extremely diverse, varying in morphology, size, nature, and complexity of their genomic content. Typically, viruses use host-encoded glycosyltransferases and glycosidases to add and remove sugar residues from their glycoproteins. Thus, the structure of the glycans on the viral proteins have, to date, typically been considered to mimick those of the host. However, the more recently discovered large and giant viruses differ from this paradigm. At least some of these viruses code for an (almost) autonomous glycosylation pathway. These viral genes include those that encode the production of activated sugars, glycosyltransferases, and other enzymes able to manipulate sugars at various levels. This review focuses on large and giant viruses that produce carbohydrate-processing enzymes. A brief description of those harboring these features at the genomic level will be discussed, followed by the achievements reached with regard to the elucidation of the glycan structures, the activity of the proteins able to manipulate sugars, and the organic synthesis of some of these virus-encoded glycans. During this progression, we will also comment on many of the challenging questions on this subject that remain to be addressed.

Simultaneous Giant Virus and Virophage Quantification Using Droplet Digital PCR

Viruses are an abundant component of aquatic systems, but their detection and quantification remain a challenge. Virophages co-replicate with giant viruses in the shared host cell, and can inhibit the production of new giant virus particles, thereby increasing the survival of the infected host population. Here, we present a protocol for Droplet Digital PCR (ddPCR) to quantify simultaneously giant virus and virophage in a mixed sample, enabling the rapid, culture-free and high throughput detection of virus and virophage. As virophage can be present as free virus particles or integrated into the virus host's genome as well as associated with organic particles, we developed a simple method that enables discrimination between free and particle-associated virophages. The latter include aggregated virophage particles as well as virophage integrated into the host genome. We used, for our experiments, a host-virus-virophage system consisting of Cafeteria burkhardae, CroV and mavirus. Our results show that ddPCR can be an efficient method to quantify virus and virophage, and we discuss potential applications of the method for studying ecological and evolutionary processes of virus and virophages.

Isolation and Identification of a Large Green Alga Virus (Chlorella Virus XW01) of Mimiviridae and Its Virophage (Chlorella Virus Virophage SW01) by Using Unicellular Green Algal Cultures

Virophages are a group of small double-stranded DNA viruses that infect protist hosts and parasitize the viral factory of host giant/large viruses to propagate. Here, we discover a novel cell-virus-virophage (CVv) tripartite interaction system by using unicellular micro-green algae (Chlorella sp.) as eukaryotic hosts for the first time. Viral particles, resembling known virophages and large alga viruses, are detected in culture supernatants and inside algal cells. Complete genomic sequences of the virophage (Chlorella virus virophage SW01 [CVv-SW01]; 24,744 bp) and large virus (Chlorella virus XW01 [CV-XW01]; 407,612 bp) are obtained from the cocultures. Both genomic and phylogenetic analyses show that CVv-SW01 is closely related to virophages previously found in Dishui Lake. CV-XW01 shares the greatest number of homologous genes (n = 82) with Cafeteria roenbergensis virus (CroV) and phylogenetically represents the closest relative to CroV. This is the first report of a large green alga virus being affiliated with a heterotrophic zooplankton-infecting Cafeteriavirus of the family Mimiviridae. Moreover, the codon usage preferences of CV-XW01 and CVv-SW01 are highly similar to those of CroV and its virophage Mavirus, respectively. The discovery of such a novel CVv system with the green alga Chlorella sp. as the single cellular eukaryotic host paves a way to further investigate the potential interaction mechanism of CVv and its significance in the ecology of green algae and the evolution of large/giant viruses and their parasitic viruses. IMPORTANCE Parasitic virophages are small unicellular eukaryotic dsDNA viruses that rely on the viral factories of coinfecting giant/large dsDNA viruses for propagation. Presently, the identified eukaryotic hosts of isolated virophages were restricted to a free-living amoeba, Acanthamoeba polyphaga, and a widespread marine heterotrophic flagellate, Cafeteria roenbergensis. In this study, we successfully discovered and identified a novel tripartite interaction system comprised of a micro-green alga (Chlorella sp.), Mimiviridae large green alga virus, and virophage at the coculture level, with Chlorella sp. as the eukaryotic host, based on combination analysis of infection, morphotype, genome, and phylogeny. The large green alga virus CV-XW01 represents the closest relative to the Mimiviridae giant virus Cafeteria roenbergensis virus, host virus of the virophage Mavirus, as well as a novel large virus of Mimiviridae that infects a non-protozoan protist host. The virophage CVv-SW01 highly resembles Mavirus in its codon usage frequency and preference, although they are phylogenetically distantly related. These findings give novel insights into the diversity of large/giant viruses and their virophages.

Metagenomic detection and characterisation of multiple viruses in apparently healthy Australian Neophema birds

Emerging viral pathogens are a significant concern, with potential consequences for human, animal and environmental health. Over the past several decades, many novel viruses have been found in animals, including birds, and often pose a significant threat to vulnerable species. However, despite enormous interest in virus research, little is known about virus communities (viromes) in Australian Neophema birds. Therefore, this study was designed to characterise the viromes of Neophema birds and track the evolutionary relationships of recently emerging psittacine siadenovirus F (PsSiAdV-F) circulating in the critically endangered, orange-bellied parrot (OBP, Neophema chrysogaster), using a viral metagenomic approach. This study identified 16 viruses belonging to the families Adenoviridae, Circoviridae, Endornaviridae, Picobirnaviridae and Picornaviridae. In addition, this study demonstrated a potential evolutionary relationship of a PsSiAdV-F sequenced previously from the critically endangered OBP. Strikingly, five adenoviral contigs identified in this study show the highest identities with human adenovirus 2 and human mastadenovirus C. This highlights an important and unexpected aspects of the avian virome and warrants further studies dedicated to this subject. Finally, the findings of this study emphasise the importance of testing birds used for trade or in experimental settings for potential pathogens to prevent the spread of infections.

Integrating Viral Metagenomics into an Ecological Framework

Viral metagenomics has expanded our knowledge of the ecology of uncultured viruses, within both environmental (e.g., terrestrial and aquatic) and host-associated (e.g., plants and animals, including humans) contexts. Here, we emphasize the implementation of an ecological framework in viral metagenomic studies to address questions in virology rarely considered ecological, which can change our perception of viruses and how they interact with their surroundings. An ecological framework explicitly considers diverse variants of viruses in populations that make up communities of interacting viruses, with ecosystem-level effects. It provides a structure for the study of the diversity, distributions, dynamics, and interactions of viruses with one another, hosts, and the ecosystem, including interactions with abiotic factors. An ecological framework in viral metagenomics stands poised to broadly expand our knowledge in basic and applied virology. We highlight specific fundamental research needs to capitalize on its potential and advance the field. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Ecology and molecular targets of hypermutation in the global microbiome

Changes in the sequence of an organism's genome, i.e., mutations, are the raw material of evolution. The frequency and location of mutations can be constrained by specific molecular mechanisms, such as diversity-generating retroelements (DGRs). DGRs have been characterized from cultivated bacteria and bacteriophages, and perform error-prone reverse transcription leading to mutations being introduced in specific target genes. DGR loci were also identified in several metagenomes, but the ecological roles and evolutionary drivers of these DGRs remain poorly understood. Here, we analyze a dataset of >30,000 DGRs from public metagenomes, establish six major lineages of DGRs including three primarily encoded by phages and seemingly used to diversify host attachment proteins, and demonstrate that DGRs are broadly active and responsible for >10% of all amino acid changes in some organisms. Overall, these results highlight the constraints under which DGRs evolve, and elucidate several distinct roles these elements play in natural communities.

Analysis of Different Size Fractions Provides a More Complete Perspective of Viral Diversity in a Freshwater Embayment

Inspired by recent discoveries of the prevalence of large viruses in the environment, we reassessed the longstanding approach of filtering water through small-pore-size filters to separate viruses from cells before metagenomic analysis. We collected samples from three sites in Hamilton Harbour, an embayment of Lake Ontario, and studied 6 data sets derived from <0.45-μm- and >0.45-μm-size fractions to compare the diversity of viruses in these fractions. At the level of virus order/family, we observed highly diverse and distinct virus communities in the >0.45-μm-size fractions, whereas the <0.45-μm-size fractions were composed primarily of Caudovirales The relative abundances of Caudovirales for which hosts could be inferred varied widely between size fractions, with higher relative abundances of cyanophages in the >0.45-μm-size fractions, potentially indicating replication within cells during ongoing infections. Many viruses of eukaryotes, such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, were detected exclusively in the often-disregarded >0.45-μm-size fractions. In addition to observing unique virus communities associated with each size fraction from every site we examined, we detected viruses common to both fractions, suggesting that these are candidates for further exploration because they could be the product of ongoing or recent lytic events. Most importantly, our observations indicate that analysis of either fraction alone provides only a partial perspective of double-stranded DNA (dsDNA) viruses in the environment, highlighting the need for more comprehensive approaches for analyzing virus communities inferred from metagenomic sequencing.IMPORTANCE Most studies of aquatic virus communities analyze DNA sequences derived from the smaller-size "free-virus" fraction. Our study demonstrates that analysis of virus communities using only the smaller-size fraction can lead to erroneously low diversity estimates for many of the larger viruses such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, whereas analyzing only the larger->0.45-μm-size fraction can lead to underestimates of Caudovirales diversity and relative abundance. Similarly, our data show that examining only the smaller-size fraction can lead to underestimations of virophage and cyanophage relative abundances that could, in turn, cause researchers to assume their limited ecological importance. Given the considerable differences we observed in this study, we recommend cautious interpretations of environmental virus community assemblages and dynamics when based on metagenomic data derived from different size fractions.

Seasonal dynamics of natural Ostreococcus viral infection at the single cell level using VirusFISH

Ostreococcus is a cosmopolitan marine genus of phytoplankton found in mesotrophic and oligotrophic waters, and the smallest free-living eukaryotes known to date, with a cell diameter close to 1 μm. Ostreococcus has been extensively studied as a model system to investigate viral-host dynamics in culture, yet the impact of viruses in naturally occurring populations is largely unknown. Here, we used Virus Fluorescence in situ Hybridization (VirusFISH) to visualize and quantify viral-host dynamics in natural populations of Ostreococcus during a seasonal cycle in the central Cantabrian Sea (Southern Bay of Biscay). Ostreococcus were predominantly found during summer and autumn at surface and 50 m depth, in coastal, mid-shelf and shelf waters, representing up to 21% of the picoeukaryotic communities. Viral infection was only detected in surface waters, and its impact was variable but highest from May to July and November to December, when up to half of the population was infected. Metatranscriptomic data available from the mid-shelf station unveiled that the Ostreococcus population was dominated by the species O. lucimarinus. This work represents a proof of concept that the VirusFISH technique can be used to quantify the impact of viruses on targeted populations of key microbes from complex natural communities.

Metagenomic sequencing reveals viral abundance and diversity in mosquitoes from the Shaanxi-Gansu-Ningxia region, China

BACKGROUND

Mosquitoes host and transmit numerous arthropod-borne viruses (arboviruses) that cause disease in both humans and animals. Effective surveillance of virome profiles in mosquitoes is vital to the prevention and control of mosquito-borne diseases in northwestern China, where epidemics occur frequently.

METHODS

Mosquitoes were collected in the Shaanxi-Gansu-Ningxia region (Shaanxi Province, Gansu Province, and Ningxia Hui Autonomous Region) of China from June to August 2019. Morphological methods were used for taxonomic identification of mosquito species. High-throughput sequencing and metagenomic analysis were used to characterize mosquito viromes.

RESULTS

A total of 22,959 mosquitoes were collected, including Culex pipiens (45.7%), Culex tritaeniorhynchus (40.6%), Anopheles sinensis (8.4%), Aedes (5.2%), and Armigeres subalbatus (0.1%). In total, 3,014,183 (0.95% of clean reads) viral sequences were identified and assigned to 116 viral species (including pathogens such as Japanese encephalitis virus and Getah virus) in 31 viral families, including Flaviviridae, Togaviridae, Phasmaviridae, Phenuiviridae, and some unclassified viruses. Mosquitoes collected in July (86 species in 26 families) showed greater viral diversity than those from June and August. Culex pipiens (69 species in 25 families) and Culex tritaeniorhynchus (73 species in 24 families) carried more viral species than Anopheles sinensis (50 species in 19 families) or Aedes (38 species in 20 families) mosquitoes.

CONCLUSION

Viral diversity and abundance were affected by mosquito species and collection time. The present study elucidates the virome compositions of various mosquito species in northwestern China, improving the understanding of virus transmission dynamics for comparison with those of disease outbreaks.

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.

Metaviromics coupled with phage-host identification to open the viral 'black box'

Viruses are found in almost all biomes on Earth, with bacteriophages (phages) accounting for the majority of viral particles in most ecosystems. Phages have been isolated from natural environments using the plaque assay and liquid medium-based dilution culturing. However, phage cultivation is restricted by the current limitations in the number of culturable bacterial strains. Unlike prokaryotes, which possess universally conserved 16S rRNA genes, phages lack universal marker genes for viral taxonomy, thus restricting cultureindependent analyses of viral diversity. To circumvent these limitations, shotgun viral metagenome sequencing (i.e., metaviromics) has been developed to enable the extensive sequencing of a variety of viral particles present in the environment and is now widely used. Using metaviromics, numerous studies on viral communities have been conducted in oceans, lakes, rivers, and soils, resulting in many novel phage sequences. Furthermore, auxiliary metabolic genes such as ammonic monooxygenase C and β-lactamase have been discovered in viral contigs assembled from viral metagenomes. Current attempts to identify putative bacterial hosts of viral metagenome sequences based on sequence homology have been limited due to viral sequence variations. Therefore, culture-independent approaches have been developed to predict bacterial hosts using single-cell genomics and fluorescentlabeling. This review focuses on recent viral metagenome studies conducted in natural environments, especially in aquatic ecosystems, and their contributions to phage ecology. Here, we concluded that although metaviromics is a key tool for the study of viral ecology, this approach must be supplemented with phage-host identification, which in turn requires the cultivation of phage-bacteria systems.

Virome composition in marine fish revealed by meta-transcriptomics

Revealing the determinants of virome composition is central to placing disease emergence in a broader evolutionary context. Fish are the most species-rich group of vertebrates and so provide an ideal model system to study the factors that shape virome compositions and their evolution. We characterized the viromes of nineteen wild-caught species of marine fish using total RNA sequencing (meta-transcriptomics) combined with analyses of sequence and protein structural homology to identify divergent viruses that often evade characterization. From this, we identified twenty-five new vertebrate-associated viruses and a further twenty-two viruses likely associated with fish diet or their microbiomes. The vertebrate-associated viruses identified here included the first fish virus in the Matonaviridae (single-strand, negative-sense RNA virus). Other viruses fell within the Astroviridae, Picornaviridae, Arenaviridae, Reoviridae, Hepadnaviridae, Paramyxoviridae, Rhabdoviridae, Hantaviridae, Filoviridae, and Flaviviridae, and were sometimes phylogenetically distinct from known fish viruses. We also show how key metrics of virome composition-viral richness, abundance, and diversity-can be analysed along with host ecological and biological factors as a means to understand virus ecology. Accordingly, these data suggest that that the vertebrate-associated viromes of the fish sampled here are predominantly shaped by the phylogenetic history (i.e. taxonomic order) of their hosts, along with several biological factors including water temperature, habitat depth, community diversity and swimming behaviour. No such correlations were found for viruses associated with porifera, molluscs, arthropods, fungi, and algae, that are unlikely to replicate in fish hosts. Overall, these data indicate that fish harbour particularly large and complex viromes and the vast majority of fish viromes are undescribed.

VirSorter2: a multi-classifier, expert-guided approach to detect diverse DNA and RNA viruses

BACKGROUND

Viruses are a significant player in many biosphere and human ecosystems, but most signals remain "hidden" in metagenomic/metatranscriptomic sequence datasets due to the lack of universal gene markers, database representatives, and insufficiently advanced identification tools.

RESULTS

Here, we introduce VirSorter2, a DNA and RNA virus identification tool that leverages genome-informed database advances across a collection of customized automatic classifiers to improve the accuracy and range of virus sequence detection. When benchmarked against genomes from both isolated and uncultivated viruses, VirSorter2 uniquely performed consistently with high accuracy (F1-score > 0.8) across viral diversity, while all other tools under-detected viruses outside of the group most represented in reference databases (i.e., those in the order Caudovirales). Among the tools evaluated, VirSorter2 was also uniquely able to minimize errors associated with atypical cellular sequences including eukaryotic genomes and plasmids. Finally, as the virosphere exploration unravels novel viral sequences, VirSorter2's modular design makes it inherently able to expand to new types of viruses via the design of new classifiers to maintain maximal sensitivity and specificity.

CONCLUSION

With multi-classifier and modular design, VirSorter2 demonstrates higher overall accuracy across major viral groups and will advance our knowledge of virus evolution, diversity, and virus-microbe interaction in various ecosystems. Source code of VirSorter2 is freely available ( https://bitbucket.org/MAVERICLab/virsorter2 ), and VirSorter2 is also available both on bioconda and as an iVirus app on CyVerse ( https://de.cyverse.org/de ). Video abstract.

Polinton-like viruses are abundant in aquatic ecosystems

BACKGROUND

Polintons are large mobile genetic elements found in the genomes of eukaryotic organisms that are considered the ancient ancestors of most eukaryotic dsDNA viruses. Originally considered as transposons, they have been found to encode virus capsid genes, suggesting they may actually be integrated viruses; however, an extracellular form has yet to be detected. Recently, circa 25 Polinton-like viruses have been discovered in environmental metagenomes and algal genomes, which shared distantly related genes to both Polintons and virophages (Lavidaviridae). These entities could be the first members of a major class of ancient eukaryotic viruses; however, owing to the lack of available genomes for analysis, information on their global diversity, evolutionary relationships, eukaryotic hosts, and status as free virus particles is limited.

RESULTS

Here, we analysed the metaviromes of an alpine lake to show that Polinton-like virus genome sequences are abundant in the water column. We identify major capsid protein genes belonging to 82 new Polinton-like viruses and use these to interrogate publicly available metagenomic datasets, identifying 543 genomes and a further 16 integrated into eukaryotic genomes. Using an analysis of shared gene content and major capsid protein phylogeny, we define large groups of Polinton-like viruses and link them to diverse eukaryotic hosts, including a new group of viruses, which possess all the core genes of virophages and infect oomycetes and Chrysophyceae.

CONCLUSIONS

Our study increased the number of known Polinton-like viruses by 25-fold, identifying five major new groups of eukaryotic viruses, which until now have been hidden in metagenomic datasets. The large enrichment (> 100-fold) of Polinton-like virus sequences in the virus-sized fraction of this alpine lake and the fact that their viral major capsid proteins are found in eukaryotic host transcriptomes support the hypothesis that Polintons in unicellular eukaryotes are viruses. In summary, our data reveals a diverse assemblage of globally distributed viruses, associated with a wide range of unicellular eukaryotic hosts. We anticipate that the methods we have developed for Polinton-like virus detection and the database of over 20,000 genes we present will allow for continued discovery and analysis of these new viral groups. Video abstract.

Characterization of a Novel Mitovirus of the Sand Fly Lutzomyia longipalpis Using Genomic and Virus-Host Interaction Signatures

Hematophagous insects act as the major reservoirs of infectious agents due to their intimate contact with a large variety of vertebrate hosts. Lutzomyia longipalpis is the main vector of Leishmania chagasi in the New World, but its role as a host of viruses is poorly understood. In this work, Lu. longipalpis RNA libraries were subjected to progressive assembly using viral profile HMMs as seeds. A sequence phylogenetically related to fungal viruses of the genus Mitovirus was identified and this novel virus was named Lul-MV-1. The 2697-base genome presents a single gene coding for an RNA-directed RNA polymerase with an organellar genetic code. To determine the possible host of Lul-MV-1, we analyzed the molecular characteristics of the viral genome. Dinucleotide composition and codon usage showed profiles similar to mitochondrial DNA of invertebrate hosts. Also, the virus-derived small RNA profile was consistent with the activation of the siRNA pathway, with size distribution and 5′ base enrichment analogous to those observed in viruses of sand flies, reinforcing Lu. longipalpis as a putative host. Finally, RT-PCR of different insect pools and sequences of public Lu. longipalpis RNA libraries confirmed the high prevalence of Lul-MV-1. This is the first report of a mitovirus infecting an insect host.

Seasonal Regime Shift in the Viral Communities of a Permafrost Thaw Lake

Permafrost thaw lakes including thermokarst lakes and ponds are ubiquitous features of Subarctic and Arctic landscapes and are hotspots of microbial activity. Input of terrestrial organic matter into the planktonic microbial loop of these lakes may greatly amplify global greenhouse gas emissions. This microbial loop, dominated in the summer by aerobic microorganisms including phototrophs, is radically different in the winter, when metabolic processes shift to the anaerobic degradation of organic matter. Little is known about the viruses that infect these microbes, despite evidence that viruses can control microbial populations and influence biogeochemical cycling in other systems. Here, we present the results of a metagenomics-based study of viruses in the larger than 0.22 µm fraction across two seasons (summer and winter) in a permafrost thaw lake in Subarctic Canada. We uncovered 351 viral populations (vOTUs) in the surface waters of this lake, with diversity significantly greater during the summer. We also identified and characterized several phage genomes and prophages, which were mostly present in the summer. Finally, we compared the viral community of this waterbody to other habitats and found unexpected similarities with distant bog lakes in North America.

Advantages and Limits of Metagenomic Assembly and Binning of a Giant Virus

Giant viruses have large genomes, often within the size range of cellular organisms. This distinguishes them from most other viruses and demands additional effort for the successful recovery of their genomes from environmental sequence data. Here, we tested the performance of genome-resolved metagenomics on a recently isolated giant virus, Fadolivirus, by spiking it into an environmental sample from which two other giant viruses were isolated. At high spike-in levels, metagenome assembly and binning led to the successful genomic recovery of Fadolivirus from the sample. A complementary survey of the major capsid protein indicated the presence of other giant viruses in the sample matrix but did not detect the two isolated from this sample. Our results indicate that genome-resolved metagenomics is a valid approach for the recovery of near-complete giant virus genomes given that sufficient clonal particles are present. However, our data also underline that a vast majority of giant viruses remain currently undetected, even in an era of terabase-scale metagenomics.IMPORTANCE The discovery of large and giant nucleocytoplasmic large DNA viruses (NCLDV) with genomes in the megabase range and equipped with a wide variety of features typically associated with cellular organisms was one of the most unexpected, intriguing, and spectacular breakthroughs in virology. Recent studies suggest that these viruses are highly abundant in the oceans, freshwater, and soil, impact the biology and ecology of their eukaryotic hosts, and ultimately affect global nutrient cycles. Genome-resolved metagenomics is becoming an increasingly popular tool to assess the diversity and coding potential of giant viruses, but this approach is currently lacking validation.

Initial Mapping of the New York City Wastewater Virome

Wastewater is a rich source of microbial life and contains bacteria, viruses, and other microbes found in human waste as well as environmental runoff sources. As part of an effort to characterize the New York City wastewater metagenome, we profiled the viral community of sewage samples across all five boroughs of NYC and found that local sampling sites have unique sets of viruses. We focused on bacteriophages, or viruses of bacteria, to understand how they may influence the microbial ecology of this system. We identified several new clusters of phages and successfully associated them with bacterial hosts, providing insight into virus-host interactions in urban wastewater. This study provides a first look into the viral communities present across the wastewater system in NYC and points to their functional importance in this environment.

Bacteriophages are abundant members of all microbiomes studied to date, influencing microbial communities through interactions with their bacterial hosts. Despite their functional importance and ubiquity, phages have been underexplored in urban environments compared to their bacterial counterparts. We profiled the viral communities in New York City (NYC) wastewater using metagenomic data collected in November 2014 from 14 wastewater treatment plants. We show that phages accounted for the largest viral component of the sewage samples and that specific virus communities were associated with local environmental conditions within boroughs. The vast majority of the virus sequences had no homology matches in public databases, forming an average of 1,700 unique virus clusters (putative genera). These new clusters contribute to elucidating the overwhelming proportion of data that frequently goes unidentified in viral metagenomic studies. We assigned potential hosts to these phages, which appear to infect a wide range of bacterial genera, often outside their presumed host. We determined that infection networks form a modular-nested pattern, indicating that phages include a range of host specificities, from generalists to specialists, with most interactions organized into distinct groups. We identified genes in viral contigs involved in carbon and sulfur cycling, suggesting functional importance of viruses in circulating pathways and gene functions in the wastewater environment. In addition, we identified virophage genes as well as a nearly complete novel virophage genome. These findings provide an understanding of phage abundance and diversity in NYC wastewater, previously uncharacterized, and further examine geographic patterns of phage-host association in urban environments.

IMPORTANCE Wastewater is a rich source of microbial life and contains bacteria, viruses, and other microbes found in human waste as well as environmental runoff sources. As part of an effort to characterize the New York City wastewater metagenome, we profiled the viral community of sewage samples across all five boroughs of NYC and found that local sampling sites have unique sets of viruses. We focused on bacteriophages, or viruses of bacteria, to understand how they may influence the microbial ecology of this system. We identified several new clusters of phages and successfully associated them with bacterial hosts, providing insight into virus-host interactions in urban wastewater. This study provides a first look into the viral communities present across the wastewater system in NYC and points to their functional importance in this environment.

Novel Cell-Virus-Virophage Tripartite Infection Systems Discovered in the Freshwater Lake Dishui Lake in Shanghai, China

Virophages are small parasitic double-stranded DNA (dsDNA) viruses of giant dsDNA viruses infecting unicellular eukaryotes. Except for a few isolated virophages characterized by parasitization mechanisms, features of virophages discovered in metagenomic data sets remain largely unknown. Here, the complete genomes of seven virophages (26.6 to 31.5 kbp) and four large DNA viruses (190.4 to 392.5 kbp) that coexist in the freshwater lake Dishui Lake, Shanghai, China, have been identified based on environmental metagenomic investigation. Both genomic and phylogenetic analyses indicate that Dishui Lake virophages (DSLVs) are closely related to each other and to other lake virophages, and Dishui Lake large DNA viruses are affiliated with the micro-green alga-infecting Prasinovirus of the Phycodnaviridae (named Dishui Lake phycodnaviruses [DSLPVs]) and protist (protozoan and alga)-infecting Mimiviridae (named Dishui Lake large alga virus [DSLLAV]). The DSLVs possess more genes with closer homology to that of large alga viruses than to that of giant protozoan viruses. Furthermore, the DSLVs are strongly associated with large green alga viruses, including DSLPV4 and DSLLAV1, based on codon usage as well as oligonucleotide frequency and correlation analyses. Surprisingly, a nonhomologous CRISPR-Cas like system is found in DSLLAV1, which appears to protect DSLLAV1 from the parasitization of DSLV5 and DSLV8. These results suggest that novel cell-virus-virophage (CVv) tripartite infection systems of green algae, large green alga virus (Phycodnaviridae- and Mimiviridae-related), and virophage exist in Dishui Lake, which will contribute to further deep investigations of the evolutionary interaction of virophages and large alga viruses as well as of the essential roles that the CVv plays in the ecology of algae.IMPORTANCE Virophages are small parasitizing viruses of large/giant viruses. To our knowledge, the few isolated virophages all parasitize giant protozoan viruses (Mimiviridae) for propagation and form a tripartite infection system with hosts, here named the cell-virus-virophage (CVv) system. However, the CVv system remains largely unknown in environmental metagenomic data sets. In this study, we systematically investigated the metagenomic data set from the freshwater lake Dishui Lake, Shanghai, China. Consequently, four novel large alga viruses and seven virophages were discovered to coexist in Dishui Lake. Surprisingly, a novel CVv tripartite infection system comprising green algae, large green alga viruses (Phycodnaviridae- and Mimiviridae-related), and virophages was identified based on genetic link, genomic signature, and CRISPR system analyses. Meanwhile, a nonhomologous CRISPR-like system was found in Dishui Lake large alga viruses, which appears to protect the virus host from the infection of Dishui Lake virophages (DSLVs). These findings are critical to give insight into the potential significance of CVv in global evolution and ecology.

Phage-specific metabolic reprogramming of virocells

Ocean viruses are abundant and infect 20–40% of surface microbes. Infected cells, termed virocells, are thus a predominant microbial state. Yet, virocells and their ecosystem impacts are understudied, thus precluding their incorporation into ecosystem models. Here we investigated how unrelated bacterial viruses (phages) reprogram one host into contrasting virocells with different potential ecosystem footprints. We independently infected the marine Pseudoalteromonas bacterium with siphovirus PSA-HS2 and podovirus PSA-HP1. Time-resolved multi-omics unveiled drastically different metabolic reprogramming and resource requirements by each virocell, which were related to phage–host genomic complementarity and viral fitness. Namely, HS2 was more complementary to the host in nucleotides and amino acids, and fitter during infection than HP1. Functionally, HS2 virocells hardly differed from uninfected cells, with minimal host metabolism impacts. HS2 virocells repressed energy-consuming metabolisms, including motility and translation. Contrastingly, HP1 virocells substantially differed from uninfected cells. They repressed host transcription, responded to infection continuously, and drastically reprogrammed resource acquisition, central carbon and energy metabolisms. Ecologically, this work suggests that one cell, infected versus uninfected, can have immensely different metabolisms that affect the ecosystem differently. Finally, we relate phage–host genome complementarity, virocell metabolic reprogramming, and viral fitness in a conceptual model to guide incorporating viruses into ecosystem models.

Charting the diversity of uncultured viruses of Archaea and Bacteria

BACKGROUND

Viruses of Archaea and Bacteria are among the most abundant and diverse biological entities on Earth. Unraveling their biodiversity has been challenging due to methodological limitations. Recent advances in culture-independent techniques, such as metagenomics, shed light on the unknown viral diversity, revealing thousands of new viral nucleotide sequences at an unprecedented scale. However, these novel sequences have not been properly classified and the evolutionary associations between them were not resolved.

RESULTS

Here, we performed phylogenomic analysis of nearly 200,000 viral nucleotide sequences to establish GL-UVAB: Genomic Lineages of Uncultured Viruses of Archaea and Bacteria. The pan-genome content of the identified lineages shed light on some of their infection strategies, potential to modulate host physiology, and mechanisms to escape host resistance systems. Furthermore, using GL-UVAB as a reference database for annotating metagenomes revealed elusive habitat distribution patterns of viral lineages and environmental drivers of community composition.

CONCLUSIONS

These findings provide insights about the genomic diversity and ecology of viruses of prokaryotes. The source code used in these analyses is freely available at https://sourceforge.net/projects/gluvab/.

Charting the diversity of uncultured viruses of Archaea and Bacteria

Background

Viruses of Archaea and Bacteria are among the most abundant and diverse biological entities on Earth. Unraveling their biodiversity has been challenging due to methodological limitations. Recent advances in culture-independent techniques, such as metagenomics, shed light on the unknown viral diversity, revealing thousands of new viral nucleotide sequences at an unprecedented scale. However, these novel sequences have not been properly classified and the evolutionary associations between them were not resolved.

Results

Here, we performed phylogenomic analysis of nearly 200,000 viral nucleotide sequences to establish GL-UVAB: Genomic Lineages of Uncultured Viruses of Archaea and Bacteria. The pan-genome content of the identified lineages shed light on some of their infection strategies, potential to modulate host physiology, and mechanisms to escape host resistance systems. Furthermore, using GL-UVAB as a reference database for annotating metagenomes revealed elusive habitat distribution patterns of viral lineages and environmental drivers of community composition.

Conclusions

These findings provide insights about the genomic diversity and ecology of viruses of prokaryotes. The source code used in these analyses is freely available at https://sourceforge.net/projects/gluvab/.

Genome-Resolved Metagenomics Extends the Environmental Distribution of the Verrucomicrobia Phylum to the Deep Terrestrial Subsurface

The Verrucomicrobia phylum of bacteria is widespread in many different ecosystems; however, its role in microbial communities remains poorly understood. Verrucomicrobia are often low-abundance community members, yet previous research suggests they play a major role in organic carbon degradation. While Verrucomicrobia remain poorly represented in culture collections, numerous genomes have been reconstructed from metagenomic data sets in recent years. The study of genomes from across the phylum allows for an extensive assessment of their potential ecosystem roles. The significance of this work is (i) the recovery of a novel genus of Verrucomicrobia from 2.3 km in the subsurface with the ability to withstand the extreme conditions that characterize this environment, and (ii) the most extensive assessment of ecophysiological traits encoded by Verrucomicrobia genomes to date. We show that members of this phylum are specialist organic polymer degraders that can withstand a wider range of environmental conditions than previously thought.

Bacteria of the phylum Verrucomicrobia are prevalent and are particularly common in soil and freshwater environments. Their cosmopolitan distribution and reported capacity for polysaccharide degradation suggests members of Verrucomicrobia are important contributors to carbon cycling across Earth’s ecosystems. Despite their prevalence, the Verrucomicrobia are underrepresented in isolate collections and genome databases; consequently, their ecophysiological roles may not be fully realized. Here, we expand genomic sampling of the Verrucomicrobia phylum by describing a novel genus, “Candidatus Marcellius,” belonging to the order Opitutales. “Ca. Marcellius” was recovered from a shale-derived produced fluid metagenome collected 313 days after hydraulic fracturing, the deepest environment from which a member of the Verrucomicrobia has been recovered to date. We uncover genomic attributes that may explain the capacity of this organism to inhabit a shale gas well, including the potential for utilization of organic polymers common in hydraulic fracturing fluids, nitrogen fixation, adaptation to high salinities, and adaptive immunity via CRISPR-Cas. To illuminate the phylogenetic and environmental distribution of these metabolic and adaptive traits across the Verrucomicrobia phylum, we performed a comparative genomic analysis of 31 publicly available, nearly complete Verrucomicrobia genomes. Our genomic findings extend the environmental distribution of the Verrucomicrobia 2.3 kilometers into the terrestrial subsurface. Moreover, we reveal traits widely encoded across members of the Verrucomicrobia, including the capacity to degrade hemicellulose and to adapt to physical and biological environmental perturbations, thereby contributing to the expansive habitat range reported for this phylum.

IMPORTANCE The Verrucomicrobia phylum of bacteria is widespread in many different ecosystems; however, its role in microbial communities remains poorly understood. Verrucomicrobia are often low-abundance community members, yet previous research suggests they play a major role in organic carbon degradation. While Verrucomicrobia remain poorly represented in culture collections, numerous genomes have been reconstructed from metagenomic data sets in recent years. The study of genomes from across the phylum allows for an extensive assessment of their potential ecosystem roles. The significance of this work is (i) the recovery of a novel genus of Verrucomicrobia from 2.3 km in the subsurface with the ability to withstand the extreme conditions that characterize this environment, and (ii) the most extensive assessment of ecophysiological traits encoded by Verrucomicrobia genomes to date. We show that members of this phylum are specialist organic polymer degraders that can withstand a wider range of environmental conditions than previously thought.

Diversity, evolution, and classification of virophages uncovered through global metagenomics

BACKGROUND

Virophages are small viruses with double-stranded DNA genomes that replicate along with giant viruses and co-infect eukaryotic cells. Due to the paucity of virophage reference genomes, a collective understanding of the global virophage diversity, distribution, and evolution is lacking.

RESULTS

Here we screened a public collection of over 14,000 metagenomes using the virophage-specific major capsid protein (MCP) as "bait." We identified 44,221 assembled virophage sequences, of which 328 represent high-quality (complete or near-complete) genomes from diverse habitats including the human gut, plant rhizosphere, and terrestrial subsurface. Comparative genomic analysis confirmed the presence of four core genes in a conserved block. We used these genes to establish a revised virophage classification including 27 clades with consistent genome length, gene content, and habitat distribution. Moreover, for eight high-quality virophage genomes, we computationally predicted putative eukaryotic virus hosts.

CONCLUSION

Overall, our approach has increased the number of known virophage genomes by 10-fold and revealed patterns of genome evolution and global virophage distribution. We anticipate that the expanded diversity presented here will provide the backbone for further virophage studies.

Genome-resolved viral and cellular metagenomes revealed potential key virus-host interactions in a deep freshwater lake

Metagenomics has dramatically expanded the known virosphere, but freshwater viral diversity and their ecological interaction with hosts remain poorly understood. Here, we conducted a metagenomic exploration of planktonic dsDNA prokaryotic viruses by sequencing both virion (<0.22 μm) and cellular (0.22-5.0 μm) fractions collected spatiotemporally from a deep freshwater lake (Lake Biwa, Japan). This simultaneously reconstructed 183 complete (i.e., circular) viral genomes and 57 bacterioplankton metagenome-assembled genomes. Analysis of metagenomic read coverage revealed vertical partitioning of the viral community analogous to the vertically stratified bacterioplankton community. The hypolimnetic community was generally stable during stratification, but occasionally shifted abruptly, presumably due to lysogenic induction. Genes involved in assimilatory sulfate reduction were encoded in 20 (10.9%) viral genomes, including those of dominant viruses, and may aid viral propagation in sulfur-limited freshwater systems. Hosts were predicted for 40 (21.9%) viral genomes, encompassing 10 phyla (or classes of Proteobacteria) including ubiquitous freshwater bacterioplankton lineages (e.g., Ca. Fonsibacter and Ca. Nitrosoarchaeum). Comparison with viral genomes derived from published metagenomes revealed viral phylogeographic connectivity in geographically isolated habitats. Notably, analogous to their hosts, actinobacterial viruses were among the most diverse, ubiquitous and abundant viral groups in freshwater systems, with potential high lytic activity in surface waters.

Phage-centric ecological interactions in aquatic ecosystems revealed through ultra-deep metagenomics

The persistent inertia in the ability to culture environmentally abundant microbes from aquatic ecosystems represents an obstacle in disentangling the complex web of ecological interactions spun by a diverse assortment of participants (pro- and eukaryotes and their viruses). In aquatic microbial communities, the numerically most abundant actors, the viruses, remain the most elusive, and especially in freshwaters their identities and ecology remain unknown. Here, using ultra-deep metagenomic sequencing from pelagic freshwater habitats, we recovered complete genomes of > 2000 phages, including small "miniphages" and large "megaphages" infecting iconic freshwater prokaryotic lineages. For instance, abundant freshwater Actinobacteria support infection by a very broad size range of phages (13-200 Kb). We describe many phages encoding genes that likely afford protection to their host from reactive oxygen species (ROS) in the aquatic environment and in the oxidative burst in protist phagolysosomes (phage-mediated ROS defense). Spatiotemporal abundance analyses of phage genomes revealed evanescence as the primary dynamic in upper water layers, where they displayed short-lived existences. In contrast, persistence was characteristic for the deeper layers where many identical phage genomes were recovered repeatedly. Phage and host abundances corresponded closely, with distinct populations displaying preferential distributions in different seasons and depths, closely mimicking overall stratification and mixis.

Phage-centric ecological interactions in aquatic ecosystems revealed through ultra-deep metagenomics

The persistent inertia in the ability to culture environmentally abundant microbes from aquatic ecosystems represents an obstacle in disentangling the complex web of ecological interactions spun by a diverse assortment of participants (pro- and eukaryotes and their viruses). In aquatic microbial communities, the numerically most abundant actors, the viruses, remain the most elusive, and especially in freshwaters their identities and ecology remain unknown. Here, using ultra-deep metagenomic sequencing from pelagic freshwater habitats, we recovered complete genomes of > 2000 phages, including small "miniphages" and large "megaphages" infecting iconic freshwater prokaryotic lineages. For instance, abundant freshwater Actinobacteria support infection by a very broad size range of phages (13-200 Kb). We describe many phages encoding genes that likely afford protection to their host from reactive oxygen species (ROS) in the aquatic environment and in the oxidative burst in protist phagolysosomes (phage-mediated ROS defense). Spatiotemporal abundance analyses of phage genomes revealed evanescence as the primary dynamic in upper water layers, where they displayed short-lived existences. In contrast, persistence was characteristic for the deeper layers where many identical phage genomes were recovered repeatedly. Phage and host abundances corresponded closely, with distinct populations displaying preferential distributions in different seasons and depths, closely mimicking overall stratification and mixis.

Metagenomic Analysis of Virus Diversity and Relative Abundance in a Eutrophic Freshwater Harbour

Aquatic viruses have been extensively studied over the past decade, yet fundamental aspects of freshwater virus communities remain poorly described. Our goal was to characterize virus communities captured in the >0.22 µm size-fraction seasonally and spatially in a freshwater harbour. Community DNA was extracted from water samples and sequenced on an Illumina HiSeq platform. Assembled contigs were annotated as belonging to the virus groups (i.e., order or family) Caudovirales, Mimiviridae, Phycodnaviridae, and virophages (Lavidaviridae), or to other groups of undefined viruses. Virophages were often the most abundant group, and discrete virophage taxa were remarkably stable across sites and dates despite fluctuations in Mimiviridae community composition. Diverse Mimiviridae contigs were detected in the samples and the two sites contained distinct Mimiviridae communities, suggesting that Mimiviridae are important algal viruses in this system. Caudovirales and Phycodnaviridae were present at low abundances in most samples. Of the 18 environmental parameters tested, only chlorophyll a explained the variation in the data at the order or family level of classification. Overall, our findings provide insight into freshwater virus community assemblages by expanding the documented diversity of freshwater virus communities, highlighting the potential ecological importance of virophages, and revealing distinct communities over small spatial scales.