Single-cell genomics-based analysis of virus-host interactions in marine surface bacterioplankton

Phage puppet masters of the marine microbial realm

Viruses numerically dominate our oceans; however, we have only just begun to document the diversity, host range and infection dynamics of marine viruses, as well as the subsequent effects of infection on both host cell metabolism and oceanic biogeochemistry. Bacteriophages (that is, phages: viruses that infect bacteria) are highly abundant and are known to play critical roles in bacterial mortality, biogeochemical cycling and horizontal gene transfer. This Review Article summarizes current knowledge of marine viral ecology and highlights the importance of phage particles to the dissolved organic matter pool, as well as the complex interactions between phages and their bacterial hosts. We emphasize the newly recognized roles of phages as puppet masters of their bacterial hosts, where phages are capable of altering the metabolism of infected bacteria through the expression of auxiliary metabolic genes and the redirection of host gene expression patterns. Finally, we propose the 'royal family model' as a hypothesis to describe successional patterns of bacteria and phages over time in marine systems, where despite high richness and significant seasonal differences, only a small number of phages appear to continually dominate a given marine ecosystem. Although further testing is required, this model provides a framework for assessing the specificity and ecological consequences of phage-host dynamics.

A virus or more in (nearly) every cell: ubiquitous networks of virus-host interactions in extreme environments

The application of viral and cellular metagenomics to natural environments has expanded our understanding of the structure, functioning, and diversity of microbial and viral communities. The high diversity of many communities, e.g., soils, surface ocean waters, and animal-associated microbiomes, make it difficult to establish virus-host associations at the single cell (rather than population) level, assign cellular hosts, or determine the extent of viral host range from metagenomics studies alone. Here, we combine single-cell sequencing with environmental metagenomics to characterize the structure of virus-host associations in a Yellowstone National Park (YNP) hot spring microbial community. Leveraging the relatively low diversity of the YNP environment, we are able to overlay evidence at the single-cell level with contextualized viral and cellular community structure. Combining evidence from hexanucelotide analysis, single cell read mapping, network-based analytics, and CRISPR-based inference, we conservatively estimate that >60% of cells contain at least one virus type and a majority of these cells contain two or more virus types. Of the detected virus types, nearly 50% were found in more than 2 cellular clades, indicative of a broad host range. The new lens provided by the combination of metaviromics and single-cell genomics reveals a network of virus-host interactions in extreme environments, provides evidence that extensive virus-host associations are common, and further expands the unseen impact of viruses on cellular life.

High Diversity of Myocyanophage in Various Aquatic Environments Revealed by High-Throughput Sequencing of Major Capsid Protein Gene With a New Set of Primers

Myocyanophages, a group of viruses infecting cyanobacteria, are abundant and play important roles in elemental cycling. Here we investigated the particle-associated viral communities retained on 0.2 μm filters and in sediment samples (representing ancient cyanophage communities) from four ocean and three lake locations, using high-throughput sequencing and a newly designed primer pair targeting a gene fragment (∼145-bp in length) encoding the cyanophage gp23 major capsid protein (MCP). Diverse viral communities were detected in all samples. The fragments of 142-, 145-, and 148-bp in length were most abundant in the amplicons, and most sequences (>92%) belonged to cyanophages. Additionally, different sequencing depths resulted in different diversity estimates of the viral community. Operational taxonomic units obtained from deep sequencing of the MCP gene covered the majority of those obtained from shallow sequencing, suggesting that deep sequencing exhibited a more complete picture of cyanophage community than shallow sequencing. Our results also revealed a wide geographic distribution of marine myocyanophages, i.e., higher dissimilarities of the myocyanophage communities corresponded with the larger distances between the sampling sites. Collectively, this study suggests that the newly designed primer pair can be effectively used to study the community and diversity of myocyanophage from different environments, and the high-throughput sequencing represents a good method to understand viral diversity.

Single-cell metagenomics: challenges and applications

With the development of high throughput sequencing and single-cell genomics technologies, many uncultured bacterial communities have been dissected by combining these two techniques. Especially, by simultaneously leveraging of single-cell genomics and metagenomics, researchers can greatly improve the efficiency and accuracy of obtaining whole genome information from complex microbial communities, which not only allow us to identify microbes but also link function to species, identify subspecies variations, study host-virus interactions and etc. Here, we review recent developments and the challenges need to be addressed in single-cell metagenomics, including potential contamination, uneven sequence coverage, sequence chimera, genome assembly and annotation. With the development of sequencing and computational methods, single-cell metagenomics will undoubtedly broaden its application in various microbiome studies.

A comprehensive and quantitative exploration of thousands of viral genomes

The complete assembly of viral genomes from metagenomic datasets (short genomic sequences gathered from environmental samples) has proven to be challenging, so there are significant blind spots when we view viral genomes through the lens of metagenomics. One approach to overcoming this problem is to leverage the thousands of complete viral genomes that are publicly available. Here we describe our efforts to assemble a comprehensive resource that provides a quantitative snapshot of viral genomic trends – such as gene density, noncoding percentage, and abundances of functional gene categories – across thousands of viral genomes. We have also developed a coarse-grained method for visualizing viral genome organization for hundreds of genomes at once, and have explored the extent of the overlap between bacterial and bacteriophage gene pools. Existing viral classification systems were developed prior to the sequencing era, so we present our analysis in a way that allows us to assess the utility of the different classification systems for capturing genomic trends.

Major role of nitrite-oxidizing bacteria in dark ocean carbon fixation

Carbon fixation by chemoautotrophic microorganisms in the dark ocean has a major impact on global carbon cycling and ecological relationships in the ocean's interior, but the relevant taxa and energy sources remain enigmatic. We show evidence that nitrite-oxidizing bacteria affiliated with the Nitrospinae phylum are important in dark ocean chemoautotrophy. Single-cell genomics and community metagenomics revealed that Nitrospinae are the most abundant and globally distributed nitrite-oxidizing bacteria in the ocean. Metaproteomics and metatranscriptomics analyses suggest that nitrite oxidation is the main pathway of energy production in Nitrospinae. Microautoradiography, linked with catalyzed reporter deposition fluorescence in situ hybridization, indicated that Nitrospinae fix 15 to 45% of inorganic carbon in the mesopelagic western North Atlantic. Nitrite oxidation may have a greater impact on the carbon cycle than previously assumed.

Archaeal Viruses from High-Temperature Environments

Archaeal viruses are some of the most enigmatic viruses known, due to the small number that have been characterized to date. The number of known archaeal viruses lags behind known bacteriophages by over an order of magnitude. Despite this, the high levels of genetic and morphological diversity that archaeal viruses display has attracted researchers for over 45 years. Extreme natural environments, such as acidic hot springs, are almost exclusively populated by Archaea and their viruses, making these attractive environments for the discovery and characterization of new viruses. The archaeal viruses from these environments have provided insights into archaeal biology, gene function, and viral evolution. This review focuses on advances from over four decades of archaeal virology, with a particular focus on archaeal viruses from high temperature environments, the existing challenges in understanding archaeal virus gene function, and approaches being taken to overcome these limitations.

Rokubacteria: Genomic Giants among the Uncultured Bacterial Phyla

Recent advances in single-cell genomic and metagenomic techniques have facilitated the discovery of numerous previously unknown, deep branches of the tree of life that lack cultured representatives. Many of these candidate phyla are composed of microorganisms with minimalistic, streamlined genomes lacking some core metabolic pathways, which may contribute to their resistance to growth in pure culture. Here we analyzed single-cell genomes and metagenome bins to show that the "Candidate phylum Rokubacteria," formerly known as SPAM, represents an interesting exception, by having large genomes (6-8 Mbps), high GC content (66-71%), and the potential for a versatile, mixotrophic metabolism. We also observed an unusually high genomic heterogeneity among individual Rokubacteria cells in the studied samples. These features may have contributed to the limited recovery of sequences of this candidate phylum in prior cultivation and metagenomic studies. Our analyses suggest that Rokubacteria are distributed globally in diverse terrestrial ecosystems, including soils, the rhizosphere, volcanic mud, oil wells, aquifers, and the deep subsurface, with no reports from marine environments to date.

Viruses of archaea: Structural, functional, environmental and evolutionary genomics

Viruses of archaea represent one of the most enigmatic parts of the virosphere. Most of the characterized archaeal viruses infect extremophilic hosts and display remarkable diversity of virion morphotypes, many of which have never been observed among viruses of bacteria or eukaryotes. The uniqueness of the virion morphologies is matched by the distinctiveness of the genomes of these viruses, with ∼75% of genes encoding unique proteins, refractory to functional annotation based on sequence analyses. In this review, we summarize the state-of-the-art knowledge on various aspects of archaeal virus genomics. First, we outline how structural and functional genomics efforts provided valuable insights into the functions of viral proteins and revealed intricate details of the archaeal virus-host interactions. We then highlight recent metagenomics studies, which provided a glimpse at the diversity of uncultivated viruses associated with the ubiquitous archaea in the oceans, including Thaumarchaeota, Marine Group II Euryarchaeota, and others. These findings, combined with the recent discovery that archaeal viruses mediate a rapid turnover of thaumarchaea in the deep sea ecosystems, illuminate the prominent role of these viruses in the biosphere. Finally, we discuss the origins and evolution of archaeal viruses and emphasize the evolutionary relationships between viruses and non-viral mobile genetic elements. Further exploration of the archaeal virus diversity as well as functional studies on diverse virus-host systems are bound to uncover novel, unexpected facets of the archaeal virome.

The enigmatic archaeal virosphere

One of the most prominent features of archaea is the extraordinary diversity of their DNA viruses. Many archaeal viruses differ substantially in morphology from bacterial and eukaryotic viruses and represent unique virus families. The distinct nature of archaeal viruses also extends to the gene composition and architectures of their genomes and the properties of the proteins that they encode. Environmental research has revealed prominent roles of archaeal viruses in influencing microbial communities in ocean ecosystems, and recent metagenomic studies have uncovered new groups of archaeal viruses that infect extremophiles and mesophiles in diverse habitats. In this Review, we summarize recent advances in our understanding of the genomic and morphological diversity of archaeal viruses and the molecular biology of their life cycles and virus-host interactions, including interactions with archaeal CRISPR-Cas systems. We also examine the potential origins and evolution of archaeal viruses and discuss their place in the global virosphere.

Analysis of single-cell genome sequences of bacteria and archaea

Single-cell genome sequencing of individual archaeal and bacterial cells is a vital approach to decipher the genetic makeup of uncultured microorganisms. With this review, we describe single-cell genome analysis with a focus on the unique properties of single-cell sequence data and with emphasis on quality assessment and assurance.

Viruses in Soil Ecosystems: An Unknown Quantity Within an Unexplored Territory

Viral abundance in soils can range from below detection limits in hot deserts to over 1 billion per gram in wetlands. Abundance appears to be strongly influenced by water availability and temperature, but a lack of informational standards creates difficulties for cross-study analysis. Soil viral diversity is severely underestimated and undersampled, although current measures of viral richness are higher for soils than for aquatic ecosystems. Both morphometric and metagenomic analyses have raised questions about the prevalence of nontailed, ssDNA viruses in soils. Soil is complex and critically important to terrestrial biodiversity and human civilization, but impacts of viral activities on soil ecosystem services are poorly understood. While information from aquatic systems and medical microbiology suggests the potential for viral influences on nutrient cycles, food web interactions, gene transfer, and other key processes in soils, very few empirical data are available. To understand the soil virome, much work remains.

Genome diversity of marine phages recovered from Mediterranean metagenomes: Size matters

Marine viruses play a critical role not only in the global geochemical cycles but also in the biology and evolution of their hosts. Despite their importance, viral diversity remains underexplored mostly due to sampling and cultivation challenges. Direct sequencing approaches such as viromics has provided new insights into the marine viral world. As a complementary approach, we analysed 24 microbial metagenomes (>0.2 μm size range) obtained from six sites in the Mediterranean Sea that vary by depth, season and filter used to retrieve the fraction. Filter-size comparison showed a significant number of viral sequences that were retained on the larger-pore filters and were different from those found in the viral fraction from the same sample, indicating that some important viral information is missing using only assembly from viromes. Besides, we were able to describe 1,323 viral genomic fragments that were more than 10Kb in length, of which 36 represented complete viral genomes including some of them retrieved from a cross-assembly from different metagenomes. Host prediction based on sequence methods revealed new phage groups belonging to marine prokaryotes like SAR11, Cyanobacteria or SAR116. We also identified the first complete virophage from deep seawater and a new endemic clade of the recently discovered Marine group II Euryarchaeota virus. Furthermore, analysis of viral distribution using metagenomes and viromes indicated that most of the new phages were found exclusively in the Mediterranean Sea and some of them, mostly the ones recovered from deep metagenomes, do not recruit in any database probably indicating higher variability and endemicity in Mediterranean bathypelagic waters. Together these data provide the first detailed picture of genomic diversity, spatial and depth variations of viral communities within the Mediterranean Sea using metagenome assembly.

Marine archaea and archaeal viruses under global change

Global change is altering oceanic temperature, salinity, pH, and oxygen concentration, directly and indirectly influencing marine microbial food web structure and function. As microbes represent >90% of the ocean’s biomass and are major drivers of biogeochemical cycles, understanding their responses to such changes is fundamental for predicting the consequences of global change on ecosystem functioning. Recent findings indicate that marine archaea and archaeal viruses are active and relevant components of marine microbial assemblages, far more abundant and diverse than was previously thought. Further research is urgently needed to better understand the impacts of global change on virus–archaea dynamics and how archaea and their viruses can interactively influence the ocean’s feedbacks on global change.

Improved genome recovery and integrated cell-size analyses of individual uncultured microbial cells and viral particles

Microbial single-cell genomics can be used to provide insights into the metabolic potential, interactions, and evolution of uncultured microorganisms. Here we present WGA-X, a method based on multiple displacement amplification of DNA that utilizes a thermostable mutant of the phi29 polymerase. WGA-X enhances genome recovery from individual microbial cells and viral particles while maintaining ease of use and scalability. The greatest improvements are observed when amplifying high G+C content templates, such as those belonging to the predominant bacteria in agricultural soils. By integrating WGA-X with calibrated index-cell sorting and high-throughput genomic sequencing, we are able to analyze genomic sequences and cell sizes of hundreds of individual, uncultured bacteria, archaea, protists, and viral particles, obtained directly from marine and soil samples, in a single experiment. This approach may find diverse applications in microbiology and in biomedical and forensic studies of humans and other multicellular organisms.

Single-cell genomics can be used to study uncultured microorganisms. Here, Stepanauskas et al. present a method combining improved multiple displacement amplification and FACS, to obtain genomic sequences and cell size information from uncultivated microbial cells and viral particles in environmental samples.

VirFinder: a novel k-mer based tool for identifying viral sequences from assembled metagenomic data

BACKGROUND

Identifying viral sequences in mixed metagenomes containing both viral and host contigs is a critical first step in analyzing the viral component of samples. Current tools for distinguishing prokaryotic virus and host contigs primarily use gene-based similarity approaches. Such approaches can significantly limit results especially for short contigs that have few predicted proteins or lack proteins with similarity to previously known viruses.

METHODS

We have developed VirFinder, the first k-mer frequency based, machine learning method for virus contig identification that entirely avoids gene-based similarity searches. VirFinder instead identifies viral sequences based on our empirical observation that viruses and hosts have discernibly different k-mer signatures. VirFinder's performance in correctly identifying viral sequences was tested by training its machine learning model on sequences from host and viral genomes sequenced before 1 January 2014 and evaluating on sequences obtained after 1 January 2014.

RESULTS

VirFinder had significantly better rates of identifying true viral contigs (true positive rates (TPRs)) than VirSorter, the current state-of-the-art gene-based virus classification tool, when evaluated with either contigs subsampled from complete genomes or assembled from a simulated human gut metagenome. For example, for contigs subsampled from complete genomes, VirFinder had 78-, 2.4-, and 1.8-fold higher TPRs than VirSorter for 1, 3, and 5 kb contigs, respectively, at the same false positive rates as VirSorter (0, 0.003, and 0.006, respectively), thus VirFinder works considerably better for small contigs than VirSorter. VirFinder furthermore identified several recently sequenced virus genomes (after 1 January 2014) that VirSorter did not and that have no nucleotide similarity to previously sequenced viruses, demonstrating VirFinder's potential advantage in identifying novel viral sequences. Application of VirFinder to a set of human gut metagenomes from healthy and liver cirrhosis patients reveals higher viral diversity in healthy individuals than cirrhosis patients. We also identified contig bins containing crAssphage-like contigs with higher abundance in healthy patients and a putative Veillonella genus prophage associated with cirrhosis patients.

CONCLUSIONS

This innovative k-mer based tool complements gene-based approaches and will significantly improve prokaryotic viral sequence identification, especially for metagenomic-based studies of viral ecology.

VirFinder: a novel k-mer based tool for identifying viral sequences from assembled metagenomic data

BACKGROUND

Identifying viral sequences in mixed metagenomes containing both viral and host contigs is a critical first step in analyzing the viral component of samples. Current tools for distinguishing prokaryotic virus and host contigs primarily use gene-based similarity approaches. Such approaches can significantly limit results especially for short contigs that have few predicted proteins or lack proteins with similarity to previously known viruses.

METHODS

We have developed VirFinder, the first k-mer frequency based, machine learning method for virus contig identification that entirely avoids gene-based similarity searches. VirFinder instead identifies viral sequences based on our empirical observation that viruses and hosts have discernibly different k-mer signatures. VirFinder's performance in correctly identifying viral sequences was tested by training its machine learning model on sequences from host and viral genomes sequenced before 1 January 2014 and evaluating on sequences obtained after 1 January 2014.

RESULTS

VirFinder had significantly better rates of identifying true viral contigs (true positive rates (TPRs)) than VirSorter, the current state-of-the-art gene-based virus classification tool, when evaluated with either contigs subsampled from complete genomes or assembled from a simulated human gut metagenome. For example, for contigs subsampled from complete genomes, VirFinder had 78-, 2.4-, and 1.8-fold higher TPRs than VirSorter for 1, 3, and 5 kb contigs, respectively, at the same false positive rates as VirSorter (0, 0.003, and 0.006, respectively), thus VirFinder works considerably better for small contigs than VirSorter. VirFinder furthermore identified several recently sequenced virus genomes (after 1 January 2014) that VirSorter did not and that have no nucleotide similarity to previously sequenced viruses, demonstrating VirFinder's potential advantage in identifying novel viral sequences. Application of VirFinder to a set of human gut metagenomes from healthy and liver cirrhosis patients reveals higher viral diversity in healthy individuals than cirrhosis patients. We also identified contig bins containing crAssphage-like contigs with higher abundance in healthy patients and a putative Veillonella genus prophage associated with cirrhosis patients.

CONCLUSIONS

This innovative k-mer based tool complements gene-based approaches and will significantly improve prokaryotic viral sequence identification, especially for metagenomic-based studies of viral ecology.

Marine viruses discovered via metagenomics shed light on viral strategies throughout the oceans

Marine viruses are key drivers of host diversity, population dynamics and biogeochemical cycling and contribute to the daily flux of billions of tons of organic matter. Despite recent advancements in metagenomics, much of their biodiversity remains uncharacterized. Here we report a data set of 27,346 marine virome contigs that includes 44 complete genomes. These outnumber all currently known phage genomes in marine habitats and include members of previously uncharacterized lineages. We designed a new method for host prediction based on co-occurrence associations that reveals these viruses infect dominant members of the marine microbiome such as Prochlorococcus and Pelagibacter. A negative association between host abundance and the virus-to-host ratio supports the recently proposed Piggyback-the-Winner model of reduced phage lysis at higher host densities. An analysis of the abundance patterns of viruses throughout the oceans revealed how marine viral communities adapt to various seasonal, temperature and photic regimes according to targeted hosts and the diversity of auxiliary metabolic genes.

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.

Putative archaeal viruses from the mesopelagic ocean

Oceanic viruses that infect bacteria, or phages, are known to modulate host diversity, metabolisms, and biogeochemical cycling, while the viruses that infect marine Archaea remain understudied despite the critical ecosystem roles played by their hosts. Here we introduce "MArVD", for Metagenomic Archaeal Virus Detector, an annotation tool designed to identify putative archaeal virus contigs in metagenomic datasets. MArVD is made publicly available through the online iVirus analytical platform. Benchmarking analysis of MArVD showed it to be >99% accurate and 100% sensitive in identifying the 127 known archaeal viruses among the 12,499 viruses in the VirSorter curated dataset. Application of MArVD to 10 viral metagenomes from two depth profiles in the Eastern Tropical North Pacific (ETNP) oxygen minimum zone revealed 43 new putative archaeal virus genomes and large genome fragments ranging in size from 10 to 31 kb. Network-based classifications, which were consistent with marker gene phylogenies where available, suggested that these putative archaeal virus contigs represented six novel candidate genera. Ecological analyses, via fragment recruitment and ordination, revealed that the diversity and relative abundances of these putative archaeal viruses were correlated with oxygen concentration and temperature along two OMZ-spanning depth profiles, presumably due to structuring of the host Archaea community. Peak viral diversity and abundances were found in surface waters, where Thermoplasmata 16S rRNA genes are prevalent, suggesting these archaea as hosts in the surface habitats. Together these findings provide a baseline for identifying archaeal viruses in sequence datasets, and an initial picture of the ecology of such viruses in non-extreme environments.

Genome and epigenome of a novel marine Thaumarchaeota strain suggest viral infection, phosphorothioation DNA modification and multiple restriction systems

Marine Thaumarchaeota are abundant ammonia-oxidizers but have few representative laboratory-cultured strains. We report the cultivation of Candidatus Nitrosomarinus catalina SPOT01, a novel strain that is less warm-temperature tolerant than other cultivated Thaumarchaeota. Using metagenomic recruitment, strain SPOT01 comprises a major portion of Thaumarchaeota (4-54%) in temperate Pacific waters. Its complete 1.36 Mbp genome possesses several distinguishing features: putative phosphorothioation (PT) DNA modification genes; a region containing probable viral genes; and putative urea utilization genes. The PT modification genes and an adjacent putative restriction enzyme (RE) operon likely form a restriction modification (RM) system for defence from foreign DNA. PacBio sequencing showed >98% methylation at two motifs, and inferred PT guanine modification of 19% of possible TGCA sites. Metagenomic recruitment also reveals the putative virus region and PT modification and RE genes are present in 18-26%, 9-14% and <1.5% of natural populations at 150 m with ≥85% identity to strain SPOT01. The presence of multiple probable RM systems in a highly streamlined genome suggests a surprising importance for defence from foreign DNA for dilute populations that infrequently encounter viruses or other cells. This new strain provides new insights into the ecology, including viral interactions, of this important group of marine microbes.

Genomic exploration of individual giant ocean viruses

Viruses are major pathogens in all biological systems. Virus propagation and downstream analysis remains a challenge, particularly in the ocean where the majority of their microbial hosts remain recalcitrant to current culturing techniques. We used a cultivation-independent approach to isolate and sequence individual viruses. The protocol uses high-speed fluorescence-activated virus sorting flow cytometry, multiple displacement amplification (MDA), and downstream genomic sequencing. We focused on 'giant viruses' that are readily distinguishable by flow cytometry. From a single-milliliter sample of seawater collected from off the dock at Boothbay Harbor, ME, USA, we sorted almost 700 single virus particles, and subsequently focused on a detailed genome analysis of 12. A wide diversity of viruses was identified that included Iridoviridae, extended Mimiviridae and even a taxonomically novel (unresolved) giant virus. We discovered a viral metacaspase homolog in one of our sorted virus particles and discussed its implications in rewiring host metabolism to enhance infection. In addition, we demonstrated that viral metacaspases are widespread in the ocean. We also discovered a virus that contains both a reverse transcriptase and a transposase; although highly speculative, we suggest such a genetic complement would potentially allow this virus to exploit a latency propagation mechanism. Application of single virus genomics provides a powerful opportunity to circumvent cultivation of viruses, moving directly to genomic investigation of naturally occurring viruses, with the assurance that the sequence data is virus-specific, non-chimeric and contains no cellular contamination.

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

Taxonomic classification of archaeal and bacterial viruses is challenging, yet also fundamental for developing a predictive understanding of microbial ecosystems. Recent identification of hundreds of thousands of new viral genomes and genome fragments, whose hosts remain unknown, requires a paradigm shift away from traditional classification approaches and towards the use of genomes for taxonomy. Here we revisited the use of genomes and their protein content as a means for developing a viral taxonomy for bacterial and archaeal viruses. A network-based analytic was evaluated and benchmarked against authority-accepted taxonomic assignments and found to be largely concordant. Exceptions were manually examined and found to represent areas of viral genome 'sequence space' that are under-sampled or prone to excessive genetic exchange. While both cases are poorly resolved by genome-based taxonomic approaches, the former will improve as viral sequence space is better sampled and the latter are uncommon. Finally, given the largely robust taxonomic capabilities of this approach, we sought to enable researchers to easily and systematically classify new viruses. Thus, we established a tool, vConTACT, as an app at iVirus, where it operates as a fast, highly scalable, user-friendly app within the free and powerful CyVerse cyberinfrastructure.

Viral coinfection is shaped by host ecology and virus-virus interactions across diverse microbial taxa and environments

Infection of more than one virus in a host, coinfection, is common across taxa and environments. Viral coinfection can enable genetic exchange, alter the dynamics of infections, and change the course of viral evolution. Yet, a systematic test of the factors explaining variation in viral coinfection across different taxa and environments awaits completion. Here I employ three microbial data sets of virus–host interactions covering cross-infectivity, culture coinfection, and single-cell coinfection (total: 6,564 microbial hosts, 13,103 viruses) to provide a broad, comprehensive picture of the ecological and biological factors shaping viral coinfection. I found evidence that ecology and virus–virus interactions are recurrent factors shaping coinfection patterns. Host ecology was a consistent and strong predictor of coinfection across all three data sets: cross-infectivity, culture coinfection, and single-cell coinfection. Host phylogeny or taxonomy was a less consistent predictor, being weak or absent in the cross-infectivity and single-cell coinfection models, yet it was the strongest predictor in the culture coinfection model. Virus–virus interactions strongly affected coinfection. In the largest test of superinfection exclusion to date, prophage sequences reduced culture coinfection by other prophages, with a weaker effect on extrachromosomal virus coinfection. At the single-cell level, prophage sequences eliminated coinfection. Virus–virus interactions also increased culture coinfection with ssDNA–dsDNA coinfections >2× more likely than ssDNA-only coinfections. The presence of CRISPR spacers was associated with a ∼50% reduction in single-cell coinfection in a marine bacteria, despite the absence of exact spacer matches in any active infection. Collectively, these results suggest the environment bacteria inhabit and the interactions among surrounding viruses are two factors consistently shaping viral coinfection patterns. These findings highlight the role of virus–virus interactions in coinfection with implications for phage therapy, microbiome dynamics, and viral infection treatments.

Molecular Analysis of Arthrobacter Myovirus vB_ArtM-ArV1: We Blame It on the Tail

This is the first report on a myophage that infects Arthrobacter A novel virus, vB_ArtM-ArV1 (ArV1), was isolated from soil using Arthrobacter sp. strain 68b for phage propagation. Transmission electron microscopy showed its resemblance to members of the family Myoviridae: ArV1 has an isometric head (∼74 nm in diameter) and a contractile, nonflexible tail (∼192 nm). Phylogenetic and comparative sequence analyses, however, revealed that ArV1 has more genes in common with phages from the family Siphoviridae than it does with any myovirus characterized to date. The genome of ArV1 is a linear, circularly permuted, double-stranded DNA molecule (71,200 bp) with a GC content of 61.6%. The genome includes 101 open reading frames (ORFs) yet contains no tRNA genes. More than 50% of ArV1 genes encode unique proteins that either have no reliable identity to database entries or have homologues only in Arthrobacter phages, both sipho- and myoviruses. Using bioinformatics approaches, 13 ArV1 structural genes were identified, including those coding for head, tail, tail fiber, and baseplate proteins. A further 6 ArV1 ORFs were annotated as encoding putative structural proteins based on the results of proteomic analysis. Phylogenetic analysis based on the alignment of four conserved virion proteins revealed that Arthrobacter myophages form a discrete clade that seems to occupy a position somewhat intermediate between myo- and siphoviruses. Thus, the data presented here will help to advance our understanding of genetic diversity and evolution of phages that constitute the order CaudoviralesIMPORTANCE Bacteriophages, which likely originated in the early Precambrian Era, represent the most numerous population on the planet. Approximately 95% of known phages are tailed viruses that comprise three families: Podoviridae (with short tails), Siphoviridae (with long noncontractile tails), and Myoviridae (with contractile tails). Based on the current hypothesis, myophages, which may have evolved from siphophages, are thought to have first emerged among Gram-negative bacteria, whereas they emerged only later among Gram-positive bacteria. The results of the molecular characterization of myophage vB_ArtM-ArV1 presented here conform to the aforementioned hypothesis, since, at a glance, bacteriophage vB_ArtM-ArV1 appears to be a siphovirus that possesses a seemingly functional contractile tail. Our work demonstrates that such "chimeric" myophages are of cosmopolitan nature and are likely characteristic of the ecologically important soil bacterial genus Arthrobacter.

Ecological dynamics and co-occurrence among marine phytoplankton, bacteria and myoviruses shows microdiversity matters

Numerous ecological processes, such as bacteriophage infection and phytoplankton-bacterial interactions, often occur via strain-specific mechanisms. Therefore, studying the causes of microbial dynamics should benefit from highly resolving taxonomic characterizations. We sampled daily to weekly over 5 months following a phytoplankton bloom off Southern California and examined the extent of microdiversity, that is, significant variation within 99% sequence similarity clusters, operational taxonomic units (OTUs), of bacteria, archaea, phytoplankton chloroplasts (all via 16S or intergenic spacer (ITS) sequences) and T4-like-myoviruses (via g23 major capsid protein gene sequence). The extent of microdiversity varied between genes (ITS most, g23 least) and only temporally common taxa were highly microdiverse. Overall, 60% of taxa exhibited microdiversity; 59% of these had subtypes that changed significantly as a proportion of the parent taxon, indicating ecologically distinct taxa. Pairwise correlations between prokaryotes and myoviruses or phytoplankton (for example, highly microdiverse Chrysochromulina sp.) improved when using single-base variants. Correlations between myoviruses and SAR11 increased in number (172 vs 9, Spearman>0.65) and became stronger (0.61 vs 0.58, t-test: P<0.001) when using SAR11 ITS single-base variants vs OTUs. Whole-community correlation between SAR11 and myoviruses was much improved when using ITS single-base variants vs OTUs, with Mantel rho=0.49 vs 0.27; these results are consistent with strain-specific interactions. Mantel correlations suggested >1 μm (attached/large) prokaryotes are a major myovirus source. Consideration of microdiversity improved observation of apparent host and virus networks, and provided insights into the ecological and evolutionary factors influencing the success of lineages, with important implications to ecosystem resilience and microbial function.

Intriguing Interaction of Bacteriophage-Host Association: An Understanding in the Era of Omics

Innovations in next-generation sequencing technology have introduced new avenues in microbial studies through “omics” approaches. This technology has considerably augmented the knowledge of the microbial world without isolation prior to their identification. With an enormous volume of bacterial “omics” data, considerable attempts have been recently invested to improve an insight into virosphere. The interplay between bacteriophages and their host has created a significant influence on the biogeochemical cycles, microbial diversity, and bacterial population regulation. This review highlights various concepts such as genomics, transcriptomics, proteomics, and metabolomics to infer the phylogenetic affiliation and function of bacteriophages and their impact on diverse microbial communities. Omics technologies illuminate the role of bacteriophage in an environment, the influences of phage proteins on the bacterial host and provide information about the genes important for interaction with bacteria. These investigations will reveal some of bio-molecules and biomarkers of the novel phage which demand to be unveiled.

Discovery and Complete Genome Sequence of a Bacteriophage from an Obligate Intracellular Symbiont of a Cellulolytic Protist in the Termite Gut

Termites depend nutritionally on their gut microbes, and protistan, bacterial, and archaeal gut communities have been extensively studied. However, limited information is available on viruses in the termite gut. We herein report the complete genome sequence (99,517 bp) of a phage obtained during a genome analysis of “Candidatus Azobacteroides pseudotrichonymphae” phylotype ProJPt-1, which is an obligate intracellular symbiont of the cellulolytic protist Pseudotrichonympha sp. in the gut of the termite Prorhinotermes japonicus. The genome of the phage, designated ProJPt-Bp1, was circular or circularly permuted, and was not integrated into the two circular chromosomes or five circular plasmids composing the host ProJPt-1 genome. The phage was putatively affiliated with the order Caudovirales based on sequence similarities with several phage-related genes; however, most of the 52 protein-coding sequences had no significant homology to sequences in the databases. The phage genome contained a tRNA-Gln (CAG) gene, which showed the highest sequence similarity to the tRNA-Gln (CAA) gene of the host “Ca. A. pseudotrichonymphae” phylotype ProJPt-1. Since the host genome lacked a tRNA-Gln (CAG) gene, the phage tRNA gene may compensate for differences in codon usage bias between the phage and host genomes. The phage genome also contained a non-coding region with high nucleotide sequence similarity to a region in one of the host plasmids. No other phage-related sequences were found in the host ProJPt-1 genome. To the best of our knowledge, this is the first report of a phage from an obligate, mutualistic endosymbiont permanently associated with eukaryotic cells.

metaSPAdes: a new versatile metagenomic assembler

While metagenomics has emerged as a technology of choice for analyzing bacterial populations, the assembly of metagenomic data remains challenging, thus stifling biological discoveries. Moreover, recent studies revealed that complex bacterial populations may be composed from dozens of related strains, thus further amplifying the challenge of metagenomic assembly. metaSPAdes addresses various challenges of metagenomic assembly by capitalizing on computational ideas that proved to be useful in assemblies of single cells and highly polymorphic diploid genomes. We benchmark metaSPAdes against other state-of-the-art metagenome assemblers and demonstrate that it results in high-quality assemblies across diverse data sets.

Lysogeny in nature: mechanisms, impact and ecology of temperate phages

Viruses that infect bacteria (phages) can influence bacterial community dynamics, bacterial genome evolution and ecosystem biogeochemistry. These influences differ depending on whether phages establish lytic, chronic or lysogenic infections. Although the first two produce virion progeny, with lytic infections resulting in cell destruction, phages undergoing lysogenic infections replicate with cells without producing virions. The impacts of lysogeny are numerous and well-studied at the cellular level, but ecosystem-level consequences remain underexplored compared to those of lytic infections. Here, we review lysogeny from molecular mechanisms to ecological patterns to emerging approaches of investigation. Our goal is to highlight both its diversity and importance in complex communities. Altogether, using a combined viral ecology toolkit that is applied across broad model systems and environments will help us understand more of the diverse lifestyles and ecological impacts of lysogens in nature.

Environmental Viral Genomes Shed New Light on Virus-Host Interactions in the Ocean

Viruses are diverse and play significant ecological roles in marine ecosystems. However, our knowledge of genome-level diversity in viruses is biased toward those isolated from few culturable hosts. Here, we determined 1,352 nonredundant complete viral genomes from marine environments. Lifting the uncertainty that clouds short incomplete sequences, whole-genome-wide analysis suggests that these environmental genomes represent hundreds of putative novel viral genera. Predicted hosts include dominant groups of marine bacteria and archaea with no isolated viruses to date. Some of the viral genomes encode many functionally related enzymes, suggesting a strong selection pressure on these marine viruses to control cellular metabolisms by accumulating genes.

Metagenomics has revealed the existence of numerous uncharacterized viral lineages, which are referred to as viral “dark matter.” However, our knowledge regarding viral genomes is biased toward culturable viruses. In this study, we analyzed 1,600 (1,352 nonredundant) complete double-stranded DNA viral genomes (10 to 211 kb) assembled from 52 marine viromes. Together with 244 previously reported uncultured viral genomes, a genome-wide comparison delineated 617 genus-level operational taxonomic units (OTUs) for these environmental viral genomes (EVGs). Of these, 600 OTUs contained no representatives from known viruses, thus putatively corresponding to novel viral genera. Predicted hosts of the EVGs included major groups of marine prokaryotes, such as marine group II Euryarchaeota and SAR86, from which no viruses have been isolated to date, as well as Flavobacteriaceae and SAR116. Our analysis indicates that marine cyanophages are already well represented in genome databases and that one of the EVGs likely represents a new cyanophage lineage. Several EVGs encode many enzymes that appear to function for an efficient utilization of iron-sulfur clusters or to enhance host survival. This suggests that there is a selection pressure on these marine viruses to accumulate genes for specific viral propagation strategies. Finally, we revealed that EVGs contribute to a 4-fold increase in the recruitment of photic-zone viromes compared with the use of current reference viral genomes.

IMPORTANCE Viruses are diverse and play significant ecological roles in marine ecosystems. However, our knowledge of genome-level diversity in viruses is biased toward those isolated from few culturable hosts. Here, we determined 1,352 nonredundant complete viral genomes from marine environments. Lifting the uncertainty that clouds short incomplete sequences, whole-genome-wide analysis suggests that these environmental genomes represent hundreds of putative novel viral genera. Predicted hosts include dominant groups of marine bacteria and archaea with no isolated viruses to date. Some of the viral genomes encode many functionally related enzymes, suggesting a strong selection pressure on these marine viruses to control cellular metabolisms by accumulating genes.

Seasonal Dynamics and Metagenomic Characterization of Marine Viruses in Goseong Bay, Korea

Viruses are the most abundant biological entities in the oceans, and account for a significant amount of the genetic diversity of marine ecosystems. However, there is little detailed information about the biodiversity of viruses in marine environments. Rapid advances in metagenomics have enabled the identification of previously unknown marine viruses. We performed metagenomic profiling of seawater samples collected at 6 sites in Goseong Bay (South Sea, Korea) during the spring, summer, autumn, and winter of 2014. The results indicated the presence of highly diverse virus communities. The DNA libraries from samples collected during four seasons were sequenced using Illumina HiSeq 2000. The number of viral reads was 136,850 during March, 70,651 during June, 66,165 during September, and 111,778 during December. Species identification indicated that Pelagibacter phage HTVC010P, Ostreococcus lucimarinus OIV5 and OIV1, and Roseobacter phage SIO1 were the most common species in all samples. For viruses with at least 10 reads, there were 204 species during March, 189 during June, 170 during September, and 173 during December. Analysis of virus families indicated that the Myoviridae was the most common during all four seasons, and viruses in the Polyomaviridae were only present during March. Viruses in the Iridoviridae were only present during three seasons. Additionally, viruses in the Iridoviridae, Herpesviridae, and Poxviridae, which may affect fish and marine animals, appeared during different seasons. These results suggest that seasonal changes in temperature contribute to the dynamic structure of the viral community in the study area. The information presented here will be useful for comparative analyses with other marine viral communities.

PCR-activated cell sorting as a general, cultivation-free method for high-throughput identification and enrichment of virus hosts

Characterizing virus-host relationships is critical for understanding the impact of a virus on an ecosystem, but is challenging with existing techniques, particularly for uncultivable species. We present a general, cultivation-free approach for identifying phage-associated bacterial cells. Using PCR-activated cell sorting, we interrogate millions of individual bacteria for the presence of specific phage nucleic acids. If the nucleic acids are present, the bacteria are recovered via sorting and their genomes analyzed. This allows targeted recovery of all possible host species in a diverse population associated with a specific phage, and can be easily targeted to identify the hosts of different phages by modifying the PCR primers used for detection. Moreover, this technique allows quantification of free phage particles, as benchmarked against the "gold standard" of virus enumeration, the plaque assay.

Assembly of long error-prone reads using de Bruijn graphs

The recent breakthroughs in assembling long error-prone reads were based on the overlap-layout-consensus (OLC) approach and did not utilize the strengths of the alternative de Bruijn graph approach to genome assembly. Moreover, these studies often assume that applications of the de Bruijn graph approach are limited to short and accurate reads and that the OLC approach is the only practical paradigm for assembling long error-prone reads. We show how to generalize de Bruijn graphs for assembling long error-prone reads and describe the ABruijn assembler, which combines the de Bruijn graph and the OLC approaches and results in accurate genome reconstructions.

Virus-mediated archaeal hecatomb in the deep seafloor

Viruses are the most abundant biological entities in the world's oceans, and they play a crucial role in global biogeochemical cycles. In deep-sea ecosystems, archaea and bacteria drive major nutrient cycles, and viruses are largely responsible for their mortality, thereby exerting important controls on microbial dynamics. However, the relative impact of viruses on archaea compared to bacteria is unknown, limiting our understanding of the factors controlling the functioning of marine systems at a global scale. We evaluate the selectivity of viral infections by using several independent approaches, including an innovative molecular method based on the quantification of archaeal versus bacterial genes released by viral lysis. We provide evidence that, in all oceanic surface sediments (from 1000- to 10,000-m water depth), the impact of viral infection is higher on archaea than on bacteria. We also found that, within deep-sea benthic archaea, the impact of viruses was mainly directed at members of specific clades of Marine Group I Thaumarchaeota. Although archaea represent, on average, ~12% of the total cell abundance in the top 50 cm of sediment, virus-induced lysis of archaea accounts for up to one-third of the total microbial biomass killed, resulting in the release of ~0.3 to 0.5 gigatons of carbon per year globally. Our results indicate that viral infection represents a key mechanism controlling the turnover of archaea in surface deep-sea sediments. We conclude that interactions between archaea and their viruses might play a profound, previously underestimated role in the functioning of deep-sea ecosystems and in global biogeochemical cycles.

The genomic content and context of auxiliary metabolic genes in marine cyanomyoviruses

Viruses of marine cyanobacteria frequently contain auxiliary metabolic genes (AMGs) that augment host metabolism during infection, but little is known about their adaptive significance. We analyzed the distribution and genomic context of 33 AMGs across 60 cyanomyovirus genomes. Similarity in AMG content among cyanomyoviruses was only weakly correlated with phylogenetic relatedness; however, AMG content was generally conserved within the same operational taxonomic unit (OTU). A virus' AMG repertoire was also correlated with its isolation host and environment (coastal versus open ocean). A new analytical method based on shared co-linear blocks revealed that variation in the genomic location of an AMG was negatively correlated with its frequency across the genomes. We propose that rare AMGs are more frequently gained or lost as a result of fluctuating selection pressures, whereas common AMGs are associated with stable selection pressures. Finally, we describe a unique cyanomyovirus (S-CAM7) that lacks many AMGs including the photosynthesis gene psbA.

Superposition of Individual Activities: Urea-Mediated Suppression of Nitrate Uptake in the Dinoflagellate Prorocentrum minimum Revealed at the Population and Single-Cell Levels

Dinoflagellates readily use diverse inorganic and organic compounds as nitrogen sources, which is advantageous in eutrophied coastal areas exposed to high loads of anthropogenic nutrients, e.g., urea, one of the most abundant organic nitrogen substrates in seawater. Cell-to-cell variability in nutritional physiology can further enhance the diversity of metabolic strategies among dinoflagellates of the same species, but it has not been studied in free-living microalgae. We applied stable isotope tracers, isotope ratio mass spectrometry and nanoscale secondary ion mass spectrometry (NanoSIMS) to investigate the response of cultured nitrate-acclimated dinoflagellates Prorocentrum minimum to a sudden input of urea and the effect of urea on the concurrent nitrate uptake at the population and single-cell levels. We demonstrate that inputs of urea lead to suppression of nitrate uptake by P. minimum, and urea uptake exceeds the concurrent uptake of nitrate. Individual dinoflagellate cells within a population display significant heterogeneity in the rates of nutrient uptake and extent of the urea-mediated inhibition of the nitrate uptake, thus forming several groups characterized by different modes of nutrition. We conclude that urea originating from sporadic sources is rapidly utilized by dinoflagellates and can be used in biosynthesis or stored intracellularly depending on the nutrient status; therefore, sudden urea inputs can represent one of the factors triggering or supporting harmful algal blooms. Significant physiological heterogeneity revealed at the single-cell level is likely to play a role in alleviation of intra-population competition for resources and can affect the dynamics of phytoplankton populations and their maintenance in natural environments.

Large-scale maps of variable infection efficiencies in aquatic Bacteroidetes phage-host model systems

Microbes drive ecosystem functioning and their viruses modulate these impacts through mortality, gene transfer and metabolic reprogramming. Despite the importance of virus-host interactions and likely variable infection efficiencies of individual phages across hosts, such variability is seldom quantified. Here, we quantify infection efficiencies of 38 phages against 19 host strains in aquatic Cellulophaga (Bacteroidetes) phage-host model systems. Binary data revealed that some phages infected only one strain while others infected 17, whereas quantitative data revealed that efficiency of infection could vary 10 orders of magnitude, even among phages within one population. This provides a baseline for understanding and modeling intrapopulation host range variation. Genera specific host ranges were also informative. For example, the Cellulophaga Microviridae, showed a markedly broader intra-species host range than previously observed in Escherichia coli systems. Further, one phage genus, Cba41, was examined to investigate nonheritable changes in plating efficiency and burst size that depended on which host strain it most recently infected. While consistent with host modification of phage DNA, no differences in nucleotide sequence or DNA modifications were detected, leaving the observation repeatable, but the mechanism unresolved. Overall, this study highlights the importance of quantitatively considering replication variations in studies of phage-host interactions.

Assessing viral taxonomic composition in benthic marine ecosystems: reliability and efficiency of different bioinformatic tools for viral metagenomic analyses

In benthic deep-sea ecosystems, which represent the largest biome on Earth, viruses have a recognised key ecological role, but their diversity is still largely unknown. Identifying the taxonomic composition of viruses is crucial for understanding virus-host interactions, their role in food web functioning and evolutionary processes. Here, we compared the performance of various bioinformatic tools (BLAST, MG-RAST, NBC, VMGAP, MetaVir, VIROME) for analysing the viral taxonomic composition in simulated viromes and viral metagenomes from different benthic deep-sea ecosystems. The analyses of simulated viromes indicate that all the BLAST tools, followed by MetaVir and VMGAP, are more reliable in the affiliation of viral sequences and strains. When analysing the environmental viromes, tBLASTx, MetaVir, VMGAP and VIROME showed a similar efficiency of sequence annotation; however, MetaVir and tBLASTx identified a higher number of viral strains. These latter tools also identified a wider range of viral families than the others, providing a wider view of viral taxonomic diversity in benthic deep-sea ecosystems. Our findings highlight strengths and weaknesses of available bioinformatic tools for investigating the taxonomic diversity of viruses in benthic ecosystems in order to improve our comprehension of viral diversity in the oceans and its relationships with host diversity and ecosystem functioning.

Optimization of viral resuspension methods for carbon-rich soils along a permafrost thaw gradient

Permafrost stores approximately 50% of global soil carbon (C) in a frozen form; it is thawing rapidly under climate change, and little is known about viral communities in these soils or their roles in C cycling. In permafrost soils, microorganisms contribute significantly to C cycling, and characterizing them has recently been shown to improve prediction of ecosystem function. In other ecosystems, viruses have broad ecosystem and community impacts ranging from host cell mortality and organic matter cycling to horizontal gene transfer and reprogramming of core microbial metabolisms. Here we developed an optimized protocol to extract viruses from three types of high organic-matter peatland soils across a permafrost thaw gradient (palsa, moss-dominated bog, and sedge-dominated fen). Three separate experiments were used to evaluate the impact of chemical buffers, physical dispersion, storage conditions, and concentration and purification methods on viral yields. The most successful protocol, amended potassium citrate buffer with bead-beating or vortexing and BSA, yielded on average as much as 2-fold more virus-like particles (VLPs) g⁻¹ of soil than other methods tested. All method combinations yielded VLPs g⁻¹ of soil on the 10⁸ order of magnitude across all three soil types. The different storage and concentration methods did not yield significantly more VLPs g⁻¹ of soil among the soil types. This research provides much-needed guidelines for resuspending viruses from soils, specifically carbon-rich soils, paving the way for incorporating viruses into soil ecology studies.

Illuminating structural proteins in viral "dark matter" with metaproteomics

Viruses are ecologically important, yet environmental virology is limited by dominance of unannotated genomic sequences representing taxonomic and functional "viral dark matter." Although recent analytical advances are rapidly improving taxonomic annotations, identifying functional dark matter remains problematic. Here, we apply paired metaproteomics and dsDNA-targeted metagenomics to identify 1,875 virion-associated proteins from the ocean. Over one-half of these proteins were newly functionally annotated and represent abundant and widespread viral metagenome-derived protein clusters (PCs). One primarily unannotated PC dominated the dataset, but structural modeling and genomic context identified this PC as a previously unidentified capsid protein from multiple uncultivated tailed virus families. Furthermore, four of the five most abundant PCs in the metaproteome represent capsid proteins containing the HK97-like protein fold previously found in many viruses that infect all three domains of life. The dominance of these proteins within our dataset, as well as their global distribution throughout the world's oceans and seas, supports prior hypotheses that this HK97-like protein fold is the most abundant biological structure on Earth. Together, these culture-independent analyses improve virion-associated protein annotations, facilitate the investigation of proteins within natural viral communities, and offer a high-throughput means of illuminating functional viral dark matter.

Understanding the ecology and evolution of host-parasite interactions across scales

Predicting the emergence, spread and evolution of parasites within and among host populations requires insight to both the spatial and temporal scales of adaptation, including an understanding of within-host up through community-level dynamics. Although there are very few pathosystems for which such extensive data exist, there has been a recent push to integrate studies performed over multiple scales or to simultaneously test for dynamics occurring across scales. Drawing on examples from the literature, with primary emphasis on three diverse host-parasite case studies, we first examine current understanding of the spatial structure of host and parasite populations, including patterns of local adaptation and spatial variation in host resistance and parasite infectivity. We then explore the ways to measure temporal variation and dynamics in host-parasite interactions and discuss the need to examine change over both ecological and evolutionary timescales. Finally, we highlight new approaches and syntheses that allow for simultaneous analysis of dynamics across scales. We argue that there is great value in examining interplay among scales in studies of host-parasite interactions.

Computational approaches to predict bacteriophage-host relationships

Metagenomics has changed the face of virus discovery by enabling the accurate identification of viral genome sequences without requiring isolation of the viruses. As a result, metagenomic virus discovery leaves the first and most fundamental question about any novel virus unanswered: What host does the virus infect? The diversity of the global virosphere and the volumes of data obtained in metagenomic sequencing projects demand computational tools for virus–host prediction. We focus on bacteriophages (phages, viruses that infect bacteria), the most abundant and diverse group of viruses found in environmental metagenomes. By analyzing 820 phages with annotated hosts, we review and assess the predictive power of in silico phage–host signals. Sequence homology approaches are the most effective at identifying known phage–host pairs. Compositional and abundance-based methods contain significant signal for phage–host classification, providing opportunities for analyzing the unknowns in viral metagenomes. Together, these computational approaches further our knowledge of the interactions between phages and their hosts. Importantly, we find that all reviewed signals significantly link phages to their hosts, illustrating how current knowledge and insights about the interaction mechanisms and ecology of coevolving phages and bacteria can be exploited to predict phage–host relationships, with potential relevance for medical and industrial applications.

New viruses infecting bacteria are increasingly being discovered in many environments through sequence-based explorations. To understand their role in microbial ecosystems, computational tools are indispensable to prioritize and guide experimental efforts. This review assesses and discusses a range of bioinformatic approaches to predict bacteriophage–host relationships when all that is known is their genome sequence.

hybridSPAdes: an algorithm for hybrid assembly of short and long reads

MOTIVATION

Recent advances in single molecule real-time (SMRT) and nanopore sequencing technologies have enabled high-quality assemblies from long and inaccurate reads. However, these approaches require high coverage by long reads and remain expensive. On the other hand, the inexpensive short reads technologies produce accurate but fragmented assemblies. Thus, a hybrid approach that assembles long reads (with low coverage) and short reads has a potential to generate high-quality assemblies at reduced cost.

RESULTS

We describe hybridSPAdes algorithm for assembling short and long reads and benchmark it on a variety of bacterial assembly projects. Our results demonstrate that hybridSPAdes generates accurate assemblies (even in projects with relatively low coverage by long reads) thus reducing the overall cost of genome sequencing. We further present the first complete assembly of a genome from single cells using SMRT reads.

AVAILABILITY AND IMPLEMENTATION

hybridSPAdes is implemented in C++ as a part of SPAdes genome assembler and is publicly available at http://bioinf.spbau.ru/en/spades

CONTACT

d.antipov@spbu.ru

SUPPLEMENTARY INFORMATION

supplementary data are available at Bioinformatics online.

Biogeography of Viruses in the Sea

Viral ecology is a rapidly progressing area of research, as molecular methods have improved significantly for targeted research on specific populations and whole communities. To interpret and synthesize global viral diversity and distribution, it is feasible to assess whether macroecology concepts can apply to marine viruses. We review how viral and host life history and physical properties can influence viral distribution in light of biogeography and metacommunity ecology paradigms. We highlight analytical approaches that can be applied to emerging global data sets and meta-analyses to identify individual taxa with global influence and drivers of emergent properties that influence microbial community structure by drawing on examples across the spectrum of viral taxa, from RNA to ssDNA and dsDNA viruses.

Viral dark matter and virus-host interactions resolved from publicly available microbial genomes

The ecological importance of viruses is now widely recognized, yet our limited knowledge of viral sequence space and virus–host interactions precludes accurate prediction of their roles and impacts. In this study, we mined publicly available bacterial and archaeal genomic data sets to identify 12,498 high-confidence viral genomes linked to their microbial hosts. These data augment public data sets 10-fold, provide first viral sequences for 13 new bacterial phyla including ecologically abundant phyla, and help taxonomically identify 7–38% of ‘unknown’ sequence space in viromes. Genome- and network-based classification was largely consistent with accepted viral taxonomy and suggested that (i) 264 new viral genera were identified (doubling known genera) and (ii) cross-taxon genomic recombination is limited. Further analyses provided empirical data on extrachromosomal prophages and coinfection prevalences, as well as evaluation of in silico virus–host linkage predictions. Together these findings illustrate the value of mining viral signal from microbial genomes.

DOI:http://dx.doi.org/10.7554/eLife.08490.001

Viruses are infectious particles that can only multiply inside the cells of microbes and other organisms. Little is known about the genetic differences between virus particles (so-called ‘genetic diversity’), especially compared to what we know about the diversity of bacteria, archaea, and other single-celled microbes. This lack of knowledge hampers our understanding of the role viruses play in the evolution of microbial communities and their associated ecosystems.

Studying the genetics of the viruses in these communities is challenging. There is no single ‘marker’ gene that can be used to identify all viruses in environmental samples. Also, many of the fragments of viral genomes that have been identified have not yet been linked to their host microbes. Many viruses integrate their genome into the DNA of their host cell, and there are computational tools available that exploit this ability to identify viruses and link them to their host. However, other viruses can live and multiply inside cells without integrating their genome into the host's DNA.

Earlier in 2015, researchers developed a new computational tool called VirSorter that can predict virus genome sequences within the DNA extracted from microbes. VirSorter identifies viral genome sequences based on the presence of ‘hallmark’ genes that encode for components found in many virus particles, together with a reference database of genomes from many viruses.

Now, Roux et al.—including some of the researchers from the earlier work—use VirSorter to predict viral DNA from publicly available bacteria and archaea genome data. The study identifies over 12,000 viral genomes and links them to their microbial hosts. These data increase the number of viral genome sequences that are publically available by a factor of ten and identify the first viruses associated with 13 new types of bacteria, which include species that are abundant in particular environments.

It is possible for several different viruses to infect a single cell at the same time. Some viruses are known to be able to exchange DNA, and if this happens frequently in other viruses, it could have a big impact on how viruses evolve. Roux et al.'s findings suggest that although it is common for several different viruses to infect the same cell, it is relatively rare for these viruses to exchange genetic material.

Roux et al.'s findings demonstrate the value of searching publicly available microbial genome data for fragments of viral genomes. These new viral genomes will serve as a useful resource for researchers as they explore the communities of viruses and microbes in natural environments, the human body and in industrial processes.

DOI:http://dx.doi.org/10.7554/eLife.08490.002

VirSorter: mining viral signal from microbial genomic data

Viruses of microbes impact all ecosystems where microbes drive key energy and substrate transformations including the oceans, humans and industrial fermenters. However, despite this recognized importance, our understanding of viral diversity and impacts remains limited by too few model systems and reference genomes. One way to fill these gaps in our knowledge of viral diversity is through the detection of viral signal in microbial genomic data. While multiple approaches have been developed and applied for the detection of prophages (viral genomes integrated in a microbial genome), new types of microbial genomic data are emerging that are more fragmented and larger scale, such as Single-cell Amplified Genomes (SAGs) of uncultivated organisms or genomic fragments assembled from metagenomic sequencing. Here, we present VirSorter, a tool designed to detect viral signal in these different types of microbial sequence data in both a reference-dependent and reference-independent manner, leveraging probabilistic models and extensive virome data to maximize detection of novel viruses. Performance testing shows that VirSorter’s prophage prediction capability compares to that of available prophage predictors for complete genomes, but is superior in predicting viral sequences outside of a host genome (i.e., from extrachromosomal prophages, lytic infections, or partially assembled prophages). Furthermore, VirSorter outperforms existing tools for fragmented genomic and metagenomic datasets, and can identify viral signal in assembled sequence (contigs) as short as 3kb, while providing near-perfect identification (>95% Recall and 100% Precision) on contigs of at least 10kb. Because VirSorter scales to large datasets, it can also be used in “reverse” to more confidently identify viral sequence in viral metagenomes by sorting away cellular DNA whether derived from gene transfer agents, generalized transduction or contamination. Finally, VirSorter is made available through the iPlant Cyberinfrastructure that provides a web-based user interface interconnected with the required computing resources. VirSorter thus complements existing prophage prediction softwares to better leverage fragmented, SAG and metagenomic datasets in a way that will scale to modern sequencing. Given these features, VirSorter should enable the discovery of new viruses in microbial datasets, and further our understanding of uncultivated viral communities across diverse ecosystems.