vConTACT: an iVirus tool to classify double-stranded DNA viruses that infect Archaea and Bacteria

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.

Gut Virome Analysis of Cameroonians Reveals High Diversity of Enteric Viruses, Including Potential Interspecies Transmitted Viruses

Diarrhea remains one of the most common causes of deaths in children. A limited number of studies have investigated the prevalence of enteric pathogens in Cameroon, and as in many other African countries, the cause of many diarrheal episodes remains unexplained. A proportion of these unknown cases of diarrhea are likely caused by yet-unidentified viral agents, some of which could be the result of (recent) interspecies transmission from animal reservoirs, like bats. Using viral metagenomics, we screened fecal samples of 221 humans (almost all with gastroenteritis symptoms) between 0 and 89 years of age with different degrees of bat contact. We identified viruses belonging to families that are known to cause gastroenteritis such as Adenoviridae, Astroviridae, Caliciviridae, Picornaviridae, and Reoviridae Interestingly, a mammalian orthoreovirus, picobirnaviruses, a smacovirus, and a pecovirus were also found. Although there was no evidence of interspecies transmission of the most common human gastroenteritis-related viruses (Astroviridae, Caliciviridae, and Reoviridae), the phylogenies of the identified orthoreovirus, picobirnavirus, and smacovirus indicate a genetic relatedness of these viruses identified in stools of humans and those of bats and/or other animals. These findings points out the possibility of interspecies transmission or simply a shared host of these viruses (bacterial, fungal, parasitic, …) present in both animals (bats) and humans. Further screening of bat viruses in humans or vice versa will elucidate the epidemiological potential threats of animal viruses to human health. Furthermore, this study showed a huge diversity of highly divergent novel phages, thereby expanding the existing phageome considerably.IMPORTANCE Despite the availability of diagnostic tools for different enteric viral pathogens, a large fraction of human cases of gastroenteritis remains unexplained. This could be due to pathogens not tested for or novel divergent viruses of potential animal origin. Fecal virome analyses of Cameroonians showed a very diverse group of viruses, some of which are genetically related to those identified in animals. This is the first attempt to describe the gut virome of humans from Cameroon. Therefore, the data represent a baseline for future studies on enteric viral pathogens in this area and contribute to our knowledge of the world's virome. The studies also highlight the fact that more viruses may be associated with diarrhea than the typical known ones. Hence, it provides meaningful epidemiological information on diarrhea-related viruses in this area.

Minimum Information about an Uncultivated Virus Genome (MIUViG)

We present an extension of the Minimum Information about any (x) Sequence (MIxS) standard for reporting sequences of uncultivated virus genomes. Minimum Information about an Uncultivated Virus Genome (MIUViG) standards were developed within the Genomic Standards Consortium framework and include virus origin, genome quality, genome annotation, taxonomic classification, biogeographic distribution and in silico host prediction. Community-wide adoption of MIUViG standards, which complement the Minimum Information about a Single Amplified Genome (MISAG) and Metagenome-Assembled Genome (MIMAG) standards for uncultivated bacteria and archaea, will improve the reporting of uncultivated virus genomes in public databases. In turn, this should enable more robust comparative studies and a systematic exploration of the global virosphere.

Finer-Scale Phylosymbiosis: Insights from Insect Viromes

Viruses are the most abundant biological entity on the planet and interact with microbial communities with which they associate. The virome of animals is often dominated by bacterial viruses, known as bacteriophages or phages, which can (re)structure bacterial communities potentially vital to the animal host. Beta diversity relationships of animal-associated bacterial communities in laboratory and wild populations frequently parallel animal phylogenetic relationships, a pattern termed phylosymbiosis. However, little is known about whether viral communities also exhibit this eco-evolutionary pattern. Metagenomics of purified viruses from recently diverged species of Nasonia parasitoid wasps reared in the lab indicates for the first time that the community relationships of the virome can also exhibit complete phylosymbiosis. Therefore, viruses, particularly bacteriophages here, may also be influenced by animal evolutionary changes either directly or indirectly through the tripartite interactions among hosts, bacteria, and phage communities. Moreover, we report several new bacteriophage genomes from the common gut bacteria in Nasonia.

Phylosymbiosis was recently proposed to describe the eco-evolutionary pattern whereby the ecological relatedness (e.g., beta diversity relationships) of host-associated microbial communities parallels the phylogeny of the host species. Representing the most abundant biological entities on the planet and common members of the animal-associated microbiome, viruses can be influential members of host-associated microbial communities that may recapitulate, reinforce, or ablate phylosymbiosis. Here we sequence the metagenomes of purified viral communities from three different parasitic wasp Nasonia species, one cytonuclear introgression line of Nasonia, and the flour moth outgroup Ephestia kuehniella. Results demonstrate complete phylosymbiosis between the viral metagenome and insect phylogeny. Across all Nasonia contigs, 69% of the genes in the viral metagenomes are either new to the databases or uncharacterized, yet over 99% of the contigs have at least one gene with similarity to a known sequence. The core Nasonia virome spans 21% of the total contigs, and the majority of that core is likely derived from induced prophages residing in the genomes of common Nasonia-associated bacterial genera: Proteus, Providencia, and Morganella. We also assemble the first complete viral particle genomes from Nasonia-associated gut bacteria. Taken together, results reveal the first complete evidence for phylosymbiosis in viral metagenomes, new genome sequences of viral particles from Nasonia-associated gut bacteria, and a large set of novel or uncharacterized genes in the Nasonia virome. This work suggests that phylosymbiosis at the host-microbiome level will likely extend to the host-virome level in other systems as well.

IMPORTANCE Viruses are the most abundant biological entity on the planet and interact with microbial communities with which they associate. The virome of animals is often dominated by bacterial viruses, known as bacteriophages or phages, which can (re)structure bacterial communities potentially vital to the animal host. Beta diversity relationships of animal-associated bacterial communities in laboratory and wild populations frequently parallel animal phylogenetic relationships, a pattern termed phylosymbiosis. However, little is known about whether viral communities also exhibit this eco-evolutionary pattern. Metagenomics of purified viruses from recently diverged species of Nasonia parasitoid wasps reared in the lab indicates for the first time that the community relationships of the virome can also exhibit complete phylosymbiosis. Therefore, viruses, particularly bacteriophages here, may also be influenced by animal evolutionary changes either directly or indirectly through the tripartite interactions among hosts, bacteria, and phage communities. Moreover, we report several new bacteriophage genomes from the common gut bacteria in Nasonia.

High-Resolution Temporal and Spatial Patterns of Virome in Wastewater Treatment Systems

Wastewater treatment plants (WWTPs) are considered reservoirs of viruses, but the diversity and dynamic changes of viruses are not well understood. In this study, we recovered 8478 metagenomic viral contigs (mVCs; >5 kb) from two WWTPs (Shatin, 2806; Shek Wu Hui, 5672) in Hong Kong. Approximately 60% of the mVCs were poorly covered (<35% of genes in identified mVCs) by the current NCBI and IMG/VR viral databases. The temporal profile of the newly identified mVCs among 98 Shatin AS samples collected monthly (for approximately 9 years) revealed the presence of periodic dynamics at an interval of approximately one year (341 days). The spatial distribution pattern of the virome in the wastewater treatment systems showed that shared viral clusters (viral populations categorized based on shared gene content and network analysis) can be globally found among similar samples of wastewater treatment systems, indicating the presence of core viral communities among geographically isolated wastewater treatment systems. These results not only supplemented the current virome database of engineered systems but also, to some extent, expanded the understanding of long-term cyclical development and spatial distributions of viral communities in wastewater treatment systems.

Evaluation of the genomic diversity of viruses infecting bacteria, archaea and eukaryotes using a common bioinformatic platform: steps towards a unified taxonomy

Genome Relationship Applied to Virus Taxonomy (GRAViTy) is a genetics-based tool that computes sequence relatedness between viruses. Composite generalized Jaccard (CGJ) distances combine measures of homology between encoded viral genes and similarities in genome organizational features (gene orders and orientations). This scoring framework effectively recapitulates the current, largely morphology and phenotypic-based, family-level classification of eukaryotic viruses. Eukaryotic virus families typically formed monophyletic groups with consistent CGJ distance cut-off dividing between and within family divergence ranges. In the current study, a parallel analysis of prokaryotic virus families revealed quite different sequence relationships, particularly those of tailed phage families (Siphoviridae, Myoviridae and Podoviridae), where members of the same family were generally far more divergent and often not detectably homologous to each other. Analysis of the 20 currently classified prokaryotic virus families indeed split them into 70 separate clusters of tailed phages genetically equivalent to family-level assignments of eukaryotic viruses. It further divided several bacterial (Sphaerolipoviridae, Tectiviridae) and archaeal (Lipothrixviridae) families. We also found that the subfamily-level groupings of tailed phages were generally more consistent with the family assignments of eukaryotic viruses, and this supports ongoing reclassifications, including Spounavirinae and Vi1virus taxa as new virus families. The current study applied a common benchmark with which to compare taxonomies of eukaryotic and prokaryotic viruses. The findings support the planned shift away from traditional morphology-based classifications of prokaryotic viruses towards a genome-based taxonomy. They demonstrate the feasibility of a unified taxonomy of viruses into which the vast body of metagenomic viral sequences may be consistently assigned.

Soil Viruses Are Underexplored Players in Ecosystem Carbon Processing

Rapidly thawing permafrost harbors ∼30 to 50% of global soil carbon, and the fate of this carbon remains unknown. Microorganisms will play a central role in its fate, and their viruses could modulate that impact via induced mortality and metabolic controls. Because of the challenges of recovering viruses from soils, little is known about soil viruses or their role(s) in microbial biogeochemical cycling. Here, we describe 53 viral populations (viral operational taxonomic units [vOTUs]) recovered from seven quantitatively derived (i.e., not multiple-displacement-amplified) viral-particle metagenomes (viromes) along a permafrost thaw gradient at the Stordalen Mire field site in northern Sweden. Only 15% of these vOTUs had genetic similarity to publicly available viruses in the RefSeq database, and ∼30% of the genes could be annotated, supporting the concept of soils as reservoirs of substantial undescribed viral genetic diversity. The vOTUs exhibited distinct ecology, with different distributions along the thaw gradient habitats, and a shift from soil-virus-like assemblages in the dry palsas to aquatic-virus-like assemblages in the inundated fen. Seventeen vOTUs were linked to microbial hosts (in silico), implicating viruses in infecting abundant microbial lineages from Acidobacteria, Verrucomicrobia, and Deltaproteobacteria, including those encoding key biogeochemical functions such as organic matter degradation. Thirty auxiliary metabolic genes (AMGs) were identified and suggested virus-mediated modulation of central carbon metabolism, soil organic matter degradation, polysaccharide binding, and regulation of sporulation. Together, these findings suggest that these soil viruses have distinct ecology, impact host-mediated biogeochemistry, and likely impact ecosystem function in the rapidly changing Arctic. IMPORTANCE This work is part of a 10-year project to examine thawing permafrost peatlands and is the first virome-particle-based approach to characterize viruses in these systems. This method yielded >2-fold-more viral populations (vOTUs) per gigabase of metagenome than vOTUs derived from bulk-soil metagenomes from the same site (J. B. Emerson, S. Roux, J. R. Brum, B. Bolduc, et al., Nat Microbiol 3:870-880, 2018, https://doi.org/10.1038/s41564-018-0190-y). We compared the ecology of the recovered vOTUs along a permafrost thaw gradient and found (i) habitat specificity, (ii) a shift in viral community identity from soil-like to aquatic-like viruses, (iii) infection of dominant microbial hosts, and (iv) carriage of host metabolic genes. These vOTUs can impact ecosystem carbon processing via top-down (inferred from lysing dominant microbial hosts) and bottom-up (inferred from carriage of auxiliary metabolic genes) controls. This work serves as a foundation which future studies can build upon to increase our understanding of the soil virosphere and how viruses affect soil ecosystem services.

Biogeography and taxonomic overview of terrestrial hot spring thermophilic phages

Bacterial viruses ("phages") play important roles in the regulation and evolution of microbial communities in most ecosystems. Terrestrial hot springs typically contain thermophilic bacterial communities, but the diversity and impacts of its associated viruses ("thermophilic phages") are largely unexplored. Here, we provide a taxonomic overview of phages that have been isolated strictly from terrestrial hot springs around the world. In addition, we placed 17 thermophilic phage genomes in a global phylogenomic context to detect evolutionary patterns. Thermophilic phages have diverse morphologies (e.g., tailed, filamentous), unique virion structures (e.g., extremely long tailed siphoviruses), and span five taxonomic families encompassing strictly thermophilic phage genera. Within the phage proteomic tree, six thermophilic phage-related clades were identified, with evident genomic relatedness between thermophilic phages and archaeal viruses. Moreover, whole proteome analyses showed clustering between phages that infect distinct host phyla, such as Firmicutes and Deinococcus-Thermus. The potential for discovery of novel phage-host systems in terrestrial hot springs remain mostly untapped, thus additional emphasis on thermophilic phages in ecological prospecting is encouraged to gain insights into the microbial population dynamics of these environments.

Viral and metabolic controls on high rates of microbial sulfur and carbon cycling in wetland ecosystems

BACKGROUND

Microorganisms drive high rates of methanogenesis and carbon mineralization in wetland ecosystems. These signals are especially pronounced in the Prairie Pothole Region of North America, the tenth largest wetland ecosystem in the world. Sulfate reduction rates up to 22 μmol cm-3 day-1 have been measured in these wetland sediments, as well as methane fluxes up to 160 mg m-2 h-1-some of the highest emissions ever measured in North American wetlands. While pore waters from PPR wetlands are characterized by high concentrations of sulfur species and dissolved organic carbon, the constraints on microbial activity are poorly understood. Here, we utilized metagenomics to investigate candidate sulfate reducers and methanogens in this ecosystem and identify metabolic and viral controls on microbial activity.

RESULTS

We recovered 162 dsrA and 206 dsrD sequences from 18 sediment metagenomes and reconstructed 24 candidate sulfate reducer genomes assigned to seven phyla. These genomes encoded the potential for utilizing a wide variety of electron donors, such as methanol and other alcohols, methylamines, and glycine betaine. We also identified 37 mcrA sequences spanning five orders and recovered two putative methanogen genomes representing the most abundant taxa-Methanosaeta and Methanoregulaceae. However, given the abundance of Methanofollis-affiliated mcrA sequences, the detection of F420-dependent alcohol dehydrogenases, and millimolar concentrations of ethanol and 2-propanol in sediment pore fluids, we hypothesize that these alcohols may drive a significant fraction of methanogenesis in this ecosystem. Finally, extensive viral novelty was detected, with approximately 80% of viral populations being unclassified at any known taxonomic levels and absent from publicly available databases. Many of these viral populations were predicted to target dominant sulfate reducers and methanogens.

CONCLUSIONS

Our results indicate that diversity is likely key to extremely high rates of methanogenesis and sulfate reduction observed in these wetlands. The inferred genomic diversity and metabolic versatility could result from dynamic environmental conditions, viral infections, and niche differentiation in the heterogeneous sediment matrix. These processes likely play an important role in modulating carbon and sulfur cycling in this ecosystem.

Virus classification - where do you draw the line?

High-throughput sequencing (HTS) and its use in recovering and assembling novel virus sequences from environmental, human clinical, veterinary and plant samples has unearthed a vast new catalogue of viruses. Their classification, known by their sequences alone, sets a major challenge to traditional virus taxonomy, especially at the family and species levels, which have been historically based largely on descriptive taxon definitions. These typically entail some knowledge of their phenotypic properties, including replication strategies, virion structure and clinical and epidemiological features, such as host range, geographical distribution and disease outcomes. Little to no information on these attributes is available, however, for viruses identified in metagenomic datasets. If such viruses are to be included in virus taxonomy, their assignments will have to be guided largely or entirely by metrics of genetic relatedness. The immediate problem here is that the International Committee on Taxonomy of Viruses (ICTV), an organisation that authorises the taxonomic classification of viruses, provides little or no guidance on how similar or how divergent viruses must be in order to be considered members of new species or new families. We have recently developed a method for scoring genomic (dis)similarity between viruses (Genome Relationships Applied to Virus Taxonomy - GRAViTy) among the eukaryotic and prokaryotic viruses currently classified by the ICTV. At the family and genus levels, we found large-scale consistency between genetic relationships and their taxonomic assignments for eukaryotic viruses of all genome configurations and genome sizes. Family assignments of prokaryotic viruses have, however, been made at a quite different genetic level, and groupings currently classified as sub-families are a much better match to the eukaryotic virus family level. These findings support the ongoing reorganisation of bacteriophage taxonomy by the ICTV Phage Study Group. A rapid and objective means to explore metagenomic viral diversity and make evidence-based assignments for such viruses at each taxonomic layer is essential. Analysis of sequences by GRAViTy provides evidence that family (and genus) assignments of currently classified viruses are largely underpinned by genomic relatedness, and these features could serve as a guide towards an evidence-based classification of metagenomic viruses in the future.

Host-linked soil viral ecology along a permafrost thaw gradient

Climate change threatens to release abundant carbon that is sequestered at high latitudes, but the constraints on microbial metabolisms that mediate the release of methane and carbon dioxide are poorly understood^1–7. The role of viruses, which are known to affect microbial dynamics, metabolism and biogeochemistry in the oceans^8–10, remains largely unexplored in soil. Here, we aimed to investigate how viruses influence microbial ecology and carbon metabolism in peatland soils along a permafrost thaw gradient in Sweden. We recovered 1,907 viral populations (genomes and large genome fragments) from 197 bulk soil and size-fractionated metagenomes, 58% of which were detected in metatranscriptomes and presumed to be active. In silico predictions linked 35% of the viruses to microbial host populations, highlighting likely viral predators of key carbon-cycling microorganisms, including methanogens and methanotrophs. Lineage-specific virus/host ratios varied, suggesting that viral infection dynamics may differentially impact microbial responses to a changing climate. Virus-encoded glycoside hydrolases, including an endomannanase with confirmed functional activity, indicated that viruses influence complex carbon degradation and that viral abundances were significant predictors of methane dynamics. These findings suggest that viruses may impact ecosystem function in climate-critical, terrestrial habitats and identify multiple potential viral contributions to soil carbon cycling.

The recovery of viral populations from peatland soils across a permafrost thaw gradient provides insights into soil viral diversity, their hosts and the potential impacts on carbon cycling in this environment.

Multiple mechanisms drive phage infection efficiency in nearly identical hosts

Phage-host interactions are critical to ecology, evolution, and biotechnology. Central to those is infection efficiency, which remains poorly understood, particularly in nature. Here we apply genome-wide transcriptomics and proteomics to investigate infection efficiency in nature's own experiment: two nearly identical (genetically and physiologically) Bacteroidetes bacterial strains (host18 and host38) that are genetically intractable, but environmentally important, where phage infection efficiency varies. On host18, specialist phage phi18:3 infects efficiently, whereas generalist phi38:1 infects inefficiently. On host38, only phi38:1 infects, and efficiently. Overall, phi18:3 globally repressed host18's transcriptome and proteome, expressed genes that likely evaded host restriction/modification (R/M) defenses and controlled its metabolism, and synchronized phage transcription with translation. In contrast, phi38:1 failed to repress host18's transcriptome and proteome, did not evade host R/M defenses or express genes for metabolism control, did not synchronize transcripts with proteins and its protein abundances were likely targeted by host proteases. However, on host38, phi38:1 globally repressed host transcriptome and proteome, synchronized phage transcription with translation, and infected host38 efficiently. Together these findings reveal multiple infection inefficiencies. While this contrasts the single mechanisms often revealed in laboratory mutant studies, it likely better reflects the phage-host interaction dynamics that occur in nature.

Biological species in the viral world

Due to their dependence on cellular organisms for metabolism and replication, viruses are typically named and assigned to species according to their genome structure and the original host that they infect. But because viruses often infect multiple hosts and the numbers of distinct lineages within a host can be vast, their delineation into species is often dictated by arbitrary sequence thresholds, which are highly inconsistent across lineages. Here we apply an approach to determine the boundaries of viral species based on the detection of gene flow within populations, thereby defining viral species according to the biological species concept (BSC). Despite the potential for gene transfer between highly divergent genomes, viruses, like the cellular organisms they infect, assort into reproductively isolated groups and can be organized into biological species. This approach revealed that BSC-defined viral species are often congruent with the taxonomic partitioning based on shared gene contents and host tropism, and that bacteriophages can similarly be classified in biological species. These results open the possibility to use a single, universal definition of species that is applicable across cellular and acellular lifeforms.

Bacteriophage research - What we have learnt and what still needs to be addressed

Research on bacteriophages has significantly enhanced our understanding of molecular biology, the genomes of prokaryotic cells, and viral ecology. Phages and lysins offer a viable alternative to the declining utility of antibiotics in this post-antibiotic era. They also provide ideal teaching tools for genomics and bioinformatics. This article touches on the first 100 years of phage research with the author commenting on what he thinks are the highlights, and what needs to be addressed.

Lysogeny is prevalent and widely distributed in the murine gut microbiota

Bacteriophages are central members and potential modulators of the gut microbiome; however, the ecological and evolutionary relationships of gut bacteria and phages are poorly understood. Here we investigated the abundance and diversity of lysogenic bacteria (lysogens) in the bacterial community of C57BL/6J mice by detecting integrated prophages in genomes reconstructed from the metagenome of commensal bacteria. For the activities of lysogens and prophages, we compared the prophage genomes with the metagenome of free phages. The majority of commensal bacteria in different taxa were identified as lysogens. More lysogens were found among Firmicutes and Proteobacteria, than among Bacteroidetes and Actinobacteria. The prophage genomes shared high sequence similarity with the metagenome of free phages, indicating that most lysogens appeared to be active, and that prophages are spontaneously induced as active phages; dietary interventions changed the composition of the induced prophages. By contrast, CRISPR-Cas systems were present in few commensal bacteria, and were rarely active against gut phages. The structure of the bacteria-phage infection networks was "nested-modular", with modularity emerging across taxonomic scales, indicating that temperate phage features have developed over a long phylogenetic timescale. We concluded that phage generalists contribute to the prevalence of lysogeny in the gut ecosystem.

The lytic Myoviridae of Enterobacteriaceae form tight recombining assemblages separated by discontinuities in genome average nucleotide identity and lateral gene flow

In Bacteria, a working consensus of species circumscription may have been reached and one of the most prominent criteria is high average nucleotide identity (ANI). ANI in effect groups strains that may recombine more or less frequently, depending on their biology, as opposed to rare interspecies gene transfer. For bacteriophages, which show various lifestyles, the nature of the fundamental natural unit, if it exists in a biological sense, is not well understood and defined. The approaches based on dot-plots are useful to group similar bacteriophages, yet are not quantitative and use arbitrarily set cut-offs. Here, we focus on lytic Myoviridae and test the ANI metric for group delineation. We show that ANI-based groups are in agreement with the International Committee on Taxonomy of Viruses (ICTV) classification and already established dot-plot groups, which are occasionally further refined owing to higher resolution of ANI analysis. Furthermore, these groups are separated among themselves by clear ANI discontinuities. Their members readily exchange core genes with each other while they do not with bacteriophages of other ANI groups, not even with the most similar. Thus, ANI-delineated groups may represent the natural units in lytic Myoviridae evolution with features that resemble those encountered in bacterial species.

Deciphering the Human Virome with Single-Virus Genomics and Metagenomics

Single-cell genomics has unveiled the metabolic potential of dominant microbes inhabiting different environments, including the human body. The lack of genomic information for predominant microbes of the human body, such as bacteriophages, hinders our ability to answer fundamental questions about our viral communities. Here, we applied single-virus genomics (SVGs) to natural human salivary samples in combination with viral metagenomics to gain some insights into the viral community structure of the oral cavity. Saliva samples were processed for viral metagenomics (n = 15) and SVGs (n = 3). A total of 1328 uncultured single viruses were sorted by fluorescence-activated virus sorting followed by whole genome amplification. Sequencing of 24 viral single amplified genomes (vSAGs) showed that half of the vSAGs contained viral hallmark genes. Among those bona fide viruses, the uncultured single virus 92-C13 putatively infecting oral Streptococcus-like species was within the top ≈10 most abundant viruses in the oral virome. Viral gene network and viral metagenomics analyses of 439 oral viruses from cultures, metagenomics, and SVGs revealed that salivary viruses were tentatively structured into ≈200 major viral clusters, corresponding to approximately genus-level groupings. Data showed that none of the publicly available viral isolates, excepting an Actinomyces phage, were significantly abundant in the oral viromes. In addition, none of the obtained viral contigs and vSAGs from this study were present in all viromes. Overall, the data demonstrates that most viral isolates are not naturally abundant in saliva, and furthermore, the predominant viruses in the oral cavity are yet uncharacterized. Results suggest a variable, complex, and interpersonal viral profile. Finally, we demonstrated the power of SVGs in combination with viral metagenomics to unveil the genetic information of the uncultured viruses of the human virome.

The genomic underpinnings of eukaryotic virus taxonomy: creating a sequence-based framework for family-level virus classification

BACKGROUND

The International Committee on Taxonomy of Viruses (ICTV) classifies viruses into families, genera and species and provides a regulated system for their nomenclature that is universally used in virus descriptions. Virus taxonomic assignments have traditionally been based upon virus phenotypic properties such as host range, virion morphology and replication mechanisms, particularly at family level. However, gene sequence comparisons provide a clearer guide to their evolutionary relationships and provide the only information that may guide the incorporation of viruses detected in environmental (metagenomic) studies that lack any phenotypic data.

RESULTS

The current study sought to determine whether the existing virus taxonomy could be reproduced by examination of genetic relationships through the extraction of protein-coding gene signatures and genome organisational features. We found large-scale consistency between genetic relationships and taxonomic assignments for viruses of all genome configurations and genome sizes. The analysis pipeline that we have called 'Genome Relationships Applied to Virus Taxonomy' (GRAViTy) was highly effective at reproducing the current assignments of viruses at family level as well as inter-family groupings into orders. Its ability to correctly differentiate assigned viruses from unassigned viruses, and classify them into the correct taxonomic group, was evaluated by threefold cross-validation technique. This predicted family membership of eukaryotic viruses with close to 100% accuracy and specificity potentially enabling the algorithm to predict assignments for the vast corpus of metagenomic sequences consistently with ICTV taxonomy rules. In an evaluation run of GRAViTy, over one half (460/921) of (near)-complete genome sequences from several large published metagenomic eukaryotic virus datasets were assigned to 127 novel family-level groupings. If corroborated by other analysis methods, these would potentially more than double the number of eukaryotic virus families in the ICTV taxonomy.

CONCLUSIONS

A rapid and objective means to explore metagenomic viral diversity and make informed recommendations for their assignments at each taxonomic layer is essential. GRAViTy provides one means to make rule-based assignments at family and order levels in a manner that preserves the integrity and underlying organisational principles of the current ICTV taxonomy framework. Such methods are increasingly required as the vast virosphere is explored.

New insight into the RNA aquatic virosphere via viromics

RNA viruses that infect microbes are now recognized as an active, persistent and important component of the aquatic microbial community. While some information about the diversity and dynamics of the RNA virioplankton has been derived from culture-based and single gene approaches, research based on viromic and metatransciptomic methods has generated unprecedented insight into this relatively understudied class of microbes. Here, the relevant literature is summarized and discussed, including viromic studies of extracellular aquatic RNA viral assemblages, and transcriptomic studies of active and associated RNA viruses from aquatic environments followed by commentary on the present challenges and future directions of this field of research.

Benchmarking viromics: an in silico evaluation of metagenome-enabled estimates of viral community composition and diversity

Background

Viral metagenomics (viromics) is increasingly used to obtain uncultivated viral genomes, evaluate community diversity, and assess ecological hypotheses. While viromic experimental methods are relatively mature and widely accepted by the research community, robust bioinformatics standards remain to be established. Here we used in silico mock viral communities to evaluate the viromic sequence-to-ecological-inference pipeline, including (i) read pre-processing and metagenome assembly, (ii) thresholds applied to estimate viral relative abundances based on read mapping to assembled contigs, and (iii) normalization methods applied to the matrix of viral relative abundances for alpha and beta diversity estimates.

Results

Tools specifically designed for metagenomes, specifically metaSPAdes, MEGAHIT, and IDBA-UD, were the most effective at assembling viromes. Read pre-processing, such as partitioning, had virtually no impact on assembly output, but may be useful when hardware is limited. Viral populations with 2–5 × coverage typically assembled well, whereas lesser coverage led to fragmented assembly. Strain heterogeneity within populations hampered assembly, especially when strains were closely related (average nucleotide identity, or ANI ≥97%) and when the most abundant strain represented <50% of the population. Viral community composition assessments based on read recruitment were generally accurate when the following thresholds for detection were applied: (i) ≥10 kb contig lengths to define populations, (ii) coverage defined from reads mapping at ≥90% identity, and (iii) ≥75% of contig length with ≥1 × coverage. Finally, although data are limited to the most abundant viruses in a community, alpha and beta diversity patterns were robustly estimated (±10%) when comparing samples of similar sequencing depth, but more divergent (up to 80%) when sequencing depth was uneven across the dataset. In the latter cases, the use of normalization methods specifically developed for metagenomes provided the best estimates.

Conclusions

These simulations provide benchmarks for selecting analysis cut-offs and establish that an optimized sample-to-ecological-inference viromics pipeline is robust for making ecological inferences from natural viral communities. Continued development to better accessing RNA, rare, and/or diverse viral populations and improved reference viral genome availability will alleviate many of viromics remaining limitations.

Single-virus genomics reveals hidden cosmopolitan and abundant viruses

Microbes drive ecosystems under constraints imposed by viruses. However, a lack of virus genome information hinders our ability to answer fundamental, biological questions concerning microbial communities. Here we apply single-virus genomics (SVGs) to assess whether portions of marine viral communities are missed by current techniques. The majority of the here-identified 44 viral single-amplified genomes (vSAGs) are more abundant in global ocean virome data sets than published metagenome-assembled viral genomes or isolates. This indicates that vSAGs likely best represent the dsDNA viral populations dominating the oceans. Species-specific recruitment patterns and virome simulation data suggest that vSAGs are highly microdiverse and that microdiversity hinders the metagenomic assembly, which could explain why their genomes have not been identified before. Altogether, SVGs enable the discovery of some of the likely most abundant and ecologically relevant marine viral species, such as vSAG 37-F6, which were overlooked by other methodologies.