Expanding the marine virosphere using metagenomics

The Evolutionary Success of the Marine Bacterium SAR11 Analyzed through a Metagenomic Perspective

As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

The SAR11 clade of Alphaproteobacteria is the most abundant group of planktonic cells in the near-surface epipelagic waters of the ocean, but the mechanisms underlying its exceptional success have not been fully elucidated. Here, we applied a metagenomic approach to explore microdiversity patterns by measuring the accumulation of synonymous and nonsynonymous mutations as well as homologous recombination in populations of SAR11 from different aquatic habitats (marine epipelagic, bathypelagic, and surface freshwater). The patterns of mutation accumulation and recombination were compared to those of other groups of representative marine microbes with multiple ecological strategies that share the same marine habitat, namely, Cyanobacteria (Prochlorococcus and Synechococcus), Archaea (“Candidatus Nitrosopelagicus” and Marine Group II Thalassoarchaea), and some heterotrophic marine bacteria (Alteromonas and Erythrobacter). SAR11 populations showed widespread recombination among distantly related members, preventing divergence leading to a genetically stable population. Moreover, their high intrapopulation sequence diversity with an enrichment in synonymous replacements supports the idea of a very ancient divergence and the coexistence of multiple different clones. However, other microbes analyzed seem to follow different evolutionary dynamics where processes of diversification driven by geographic and ecological instability produce a higher number of nonsynonymous replacements and lower intrapopulation sequence diversity. Together, these data shed light on some of the evolutionary and ecological processes that lead to the large genomic diversity in SAR11. Furthermore, this approach can be applied to other similar microbes that are difficult to culture in the laboratory, but abundant in nature, to investigate the underlying dynamics of their genomic evolution.

IMPORTANCE As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

Single-virus genomics and beyond

Viruses are extremely diverse and modulate important biological and ecological processes globally. However, much of viral diversity remains uncultured and yet to be discovered. Several powerful culture-independent tools, in particular metagenomics, have substantially advanced virus discovery. Among those tools is single-virus genomics, which yields sequenced reference genomes from individual sorted virus particles without the need for cultivation. This new method complements virus culturing and metagenomic approaches and its advantages include targeted investigation of specific virus groups and investigation of genomic microdiversity within viral populations. In this Review, we provide a brief history of single-virus genomics, outline how this emergent method has facilitated advances in virus ecology and discuss its current limitations and future potential. Finally, we address how this method may synergistically intersect with other single-virus and single-cell approaches.

Identification, Characterization, and Genomic Analysis of Novel Serratia Temperate Phages from a Gold Mine

Bacteria of the genus Serratia inhabit a variety of ecological niches like water, soil, and the bodies of animals, and have a wide range of lifestyles. Currently, the complete genome sequences of 25 Serratia phages are available in the NCBI database. All of them were isolated from nutrient-rich environments like sewage, with the use of clinical Serratia strains as hosts. In this study, we identified a novel Serratia myovirus named vB_SspM_BZS1. Both the phage and its host Serratia sp. OS31 were isolated from the same oligotrophic environment, namely, an abandoned gold mine (Zloty Stok, Poland). The BZS1 phage was thoroughly characterized here in terms of its genomics, morphology, and infection kinetics. We also demonstrated that Serratia sp. OS31 was lysogenized by mitomycin-inducible siphovirus vB_SspS_OS31. Comparative analyses revealed that vB_SspM_BZS1 and vB_SspS_OS31 were remote from the known Serratia phages. Moreover, vB_SspM_BZS1 was only distantly related to other viruses. However, we discovered similar prophage sequences in genomes of various bacteria here. Additionally, a protein-based similarity network showed a high diversity of Serratia phages in general, as they were scattered across nineteen different clusters. In summary, this work broadened our knowledge on the diverse relationships of Serratia phages.

Extraordinary diversity of viruses in deep-sea sediments as revealed by metagenomics without prior virion separation

Our current knowledge of the virosphere in deep-sea sediments remains rudimentary. Here we investigated viral diversity at both gene and genomic levels in deep-sea sediments of Southwest Indian Ocean. Analysis of 19 676 106 non-redundant genes from the metagenomic DNA sequences revealed a large number of unclassified viral groups in these samples. A total of 1106 high-confidence viral contigs were obtained after two runs of assemblies, and 217 of these contigs with sizes up to ~120 kb were shown to represent complete viral genomes. These contigs are clustered with no known viral genomes, and over 2/3 of the ORFs on the viral contigs encode no known functions. Furthermore, most of the complete viral contigs show limited similarity to known viral genomes in genome organization. Most of the classified viral contigs are derived from dsDNA viruses belonging to the order Caudovirales, including primarily members of the families Myoviridae, Podoviridae and Siphoviridae. Most of these viruses infect Proteobacteria and, less frequently, Planctomycetes, Firmicutes, Chloroflexi, etc. Auxiliary metabolic genes (AMGs), present in abundance on the viral contigs, appear to function in modulating the host ability to sense environmental gradients and community changes, and to uptake and metabolize nutrients.

Lysogenic host-virus interactions in SAR11 marine bacteria

Host-virus interactions structure microbial communities, drive biogeochemical cycles and enhance genetic diversity in nature1,2. Hypotheses proposed to explain the range of interactions that mediate these processes often invoke lysogeny3-6, a latent infection strategy used by temperate bacterial viruses to replicate in host cells until an induction event triggers the production and lytic release of free viruses. Most cultured bacteria harbour temperate viruses in their genomes (prophage)7. The absence of prophages in cultures of the dominant lineages of marine bacteria has contributed to an ongoing debate over the ecological significance of lysogeny and other viral life strategies in nature6,8-15. Here, we report the discovery of prophages in cultured SAR11, the ocean's most abundant clade of heterotrophic bacteria16,17. We show the concurrent production of cells and viruses, with enhanced virus production under carbon-limiting growth conditions. Evidence that related prophages are broadly distributed in the oceans suggests that similar interactions have contributed to the evolutionary success of SAR11 in nutrient-limited systems.

Comparative Genomics and Evolutionary Analysis of RNA-Binding Proteins of the CsrA Family in the Genus Pseudomonas

Gene expression is adjusted according to cellular needs through a combination of mechanisms acting at different layers of the flow of genetic information. At the posttranscriptional level, RNA-binding proteins are key factors controlling the fate of nascent and mature mRNAs. Among them, the members of the CsrA family are small dimeric proteins with heterogeneous distribution across the bacterial tree of life, that act as global regulators of gene expression because they recognize characteristic sequence/structural motifs (short hairpins with GGA triplets in the loop) present in hundreds of mRNAs. The regulatory output of CsrA binding to mRNAs is counteracted in most cases by molecular mimic, non-protein coding RNAs that titrate the CsrA dimers away from the target mRNAs. In γ-proteobacteria, the regulatory modules composed by CsrA homologs and the corresponding antagonistic sRNAs, are mastered by two-component systems of the GacS-GacA type, which control the transcription and the abundance of the sRNAs, thus constituting the rather linear cascade Gac-Rsm that responds to environmental or cellular signals to adjust and coordinate the expression of a set of target genes posttranscriptionally. Within the γ-proteobacteria, the genus Pseudomonas has been shown to contain species with different number of active CsrA (RsmA) homologs and of molecular mimic sRNAs. Here, with the help of the increasing availability of genomic data we provide a comprehensive state-of-the-art picture of the remarkable multiplicity of CsrA lineages, including novel yet uncharacterized paralogues, and discuss evolutionary aspects of the CsrA subfamilies of the genus Pseudomonas, and implications of the striking presence of csrA alleles in natural mobile genetic elements (phages and plasmids).

What Is Metagenomics Teaching Us, and What Is Missed?

Shotgun metagenomic sequencing has revolutionized our ability to detect and characterize the diversity and function of complex microbial communities. In this review, we highlight the benefits of using metagenomics as well as the breadth of conclusions that can be made using currently available analytical tools, such as greater resolution of species and strains across phyla and functional content, while highlighting challenges of metagenomic data analysis. Major challenges remain in annotating function, given the dearth of functional databases for environmental bacteria compared to model organisms, and the technical difficulties of metagenome assembly and phasing in heterogeneous environmental samples. In the future, improvements and innovation in technology and methodology will lead to lowered costs. Data integration using multiple technological platforms will lead to a better understanding of how to harness metagenomes. Subsequently, we will be able not only to characterize complex microbiomes but also to manipulate communities to achieve prosperous outcomes for health, agriculture, and environmental sustainability.

A network-based integrated framework for predicting virus-prokaryote interactions

Metagenomic sequencing has greatly enhanced the discovery of viral genomic sequences; however, it remains challenging to identify the host(s) of these new viruses. We developed VirHostMatcher-Net, a flexible, network-based, Markov random field framework for predicting virus–prokaryote interactions using multiple, integrated features: CRISPR sequences and alignment-free similarity measures (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$s_2^*$\end{document} and WIsH). Evaluation of this method on a benchmark set of 1462 known virus–prokaryote pairs yielded host prediction accuracy of 59% and 86% at the genus and phylum levels, representing 16–27% and 6–10% improvement, respectively, over previous single-feature prediction approaches. We applied our host prediction tool to crAssphage, a human gut phage, and two metagenomic virus datasets: marine viruses and viral contigs recovered from globally distributed, diverse habitats. Host predictions were frequently consistent with those of previous studies, but more importantly, this new tool made many more confident predictions than previous tools, up to nearly 3-fold more (n > 27 000), greatly expanding the diversity of known virus–host interactions.

Computational approaches in viral ecology

Dynamic virus-host interactions play a critical role in regulating microbial community structure and function. Yet for decades prior to the genomics era, viruses were largely overlooked in microbial ecology research, as only low-throughput culture-based methods of discovering viruses were available. With the advent of metagenomics, culture-independent techniques have provided exciting opportunities to discover and study new viruses. Here, we review recently developed computational methods for identifying viral sequences, exploring viral diversity in environmental samples, and predicting hosts from metagenomic sequence data. Methods to analyze viruses in silico utilize unconventional approaches to tackle challenges unique to viruses, such as vast diversity, mosaic viral genomes, and the lack of universal marker genes. As the field of viral ecology expands exponentially, computational advances have become increasingly important to gain insight into the role viruses in diverse habitats.

Alteromonas Myovirus V22 Represents a New Genus of Marine Bacteriophages Requiring a Tail Fiber Chaperone for Host Recognition

Marine phages play a variety of critical roles in regulating the microbial composition of our oceans. Despite constituting the majority of genetic diversity within these environments, there are relatively few isolates with complete genome sequences or in-depth analyses of their host interaction mechanisms, such as characterization of their receptor binding proteins (RBPs). Here, we present the 92,760-bp genome of the Alteromonas-targeting phage V22. Genomic and morphological analyses identify V22 as a myovirus; however, due to a lack of sequence similarity to any other known myoviruses, we propose that V22 be classified as the type phage of a new Myoalterovirus genus within the Myoviridae family. V22 shows gene homology and synteny with two different subfamilies of phages infecting enterobacteria, specifically within the structural region of its genome. To improve our understanding of the V22 adsorption process, we identified putative RBPs (gp23, gp24, and gp26) and tested their ability to decorate the V22 propagation strain, Alteromonas mediterranea PT11, as recombinant green fluorescent protein (GFP)-tagged constructs. Only GFP-gp26 was capable of bacterial recognition and identified as the V22 RBP. Interestingly, production of functional GFP-gp26 required coexpression with the downstream protein gp27. GFP-gp26 could be expressed alone but was incapable of host recognition. By combining size-exclusion chromatography with fluorescence microscopy, we reveal how gp27 is not a component of the final RBP complex but instead is identified as a new type of phage-encoded intermolecular chaperone that is essential for maturation of the gp26 RBP.IMPORTANCE Host recognition by phage-encoded receptor binding proteins (RBPs) constitutes the first step in all phage infections and the most critical determinant of host specificity. By characterizing new types of RBPs and identifying their essential chaperones, we hope to expand the repertoire of known phage-host recognition machineries. Due to their genetic plasticity, studying RBPs and their associated chaperones can shed new light onto viral evolution affecting phage-host interactions, which is essential for fields such as phage therapy or biotechnology. In addition, since marine phages constitute one of the most important reservoirs of noncharacterized genetic diversity on the planet, their genomic and functional characterization may be of paramount importance for the discovery of novel genes with potential applications.

Metatranscriptomic Analysis of Oil-Exposed Seawater Bacterial Communities Archived by an Environmental Sample Processor (ESP)

The use of natural marine bacteria as "oil sensors" for the detection of pollution events can be suggested as a novel way of monitoring oil occurrence at sea. Nucleic acid-based devices generically called genosensors are emerging as potentially promising tools for in situ detection of specific microbial marker genes suited for that purpose. Functional marker genes are particularly interesting as targets for oil-related genosensing but their identification remains a challenge. Here, seawater samples, collected in tanks with oil addition mimicking a realistic oil spill scenario, were filtered and archived by the Environmental Sample Processor (ESP), a fully robotized genosensor, and the samples were then used for post-retrieval metatranscriptomic analysis. After extraction, RNA from ESP-archived samples at start, Day 4 and Day 7 of the experiment was used for sequencing. Metatranscriptomics revealed that several KEGG pathways were significantly enriched in samples exposed to oil. However, these pathways were highly expressed also in the non-oil-exposed water samples, most likely as a result of the release of natural organic matter from decaying phytoplankton. Temporary peaks of aliphatic alcohol and aldehyde dehydrogenases and monoaromatic ring-degrading enzymes (e.g., ben, box, and dmp clusters) were observed on Day 4 in both control and oil-exposed and non-exposed tanks. Few alkane 1-monooxygenase genes were upregulated on oil, mostly transcribed by families Porticoccaceae and Rhodobacteraceae, together with aromatic ring-hydroxylating dioxygenases, mostly transcribed by Rhodobacteraceae. Few transcripts from obligate hydrocarbonoclastic genera of Alcanivorax, Oleispira and Cycloclasticus were significantly enriched in the oil-treated exposed tank in comparison to control the non-exposed tank, and these were mostly transporters and genes involved in nitrogen and phosphorous acquisition. This study highlights the importance of seasonality, i.e., phytoplankton occurrence and senescence leading to organic compound release which can be used preferentially by bacteria over oil compounds, delaying the latter process. As a result, such seasonal effect can reduce the sensitivity of genosensing tools employing bacterial functional genes to sense oil. A better understanding of the use of natural organic matter by bacteria involved in oil-biodegradation is needed to develop an array of functional markers enabling the rapid and specific in situ detection of anthropogenic pollution.

Diel cycling of the cosmopolitan abundant Pelagibacter virus 37-F6: one of the most abundant viruses on earth

The spatiotemporal dynamics for marine viral populations has only recently been explored. However, nothing is known about temporal activities of the uncultured Pelagibacter virus vSAG 37-F6, which was discovered by single-virus genomics as potentially the most abundant marine virus. Here, we investigate the diel cycling of 37-F6 virus and the putative SAR11 host using coastal and oceanic transcriptomic and viromic time-series data from Osaka Bay and North Pacific Subtropical Gyre. Virus 37-F6 and relatives displayed diel cycling of transcriptional activities synchronized with its putative host. In both virus and host, the lowest transcription rates were observed at 14:00-15:00, coinciding roughly with maximum solar irradiance, while higher transcriptional rates were detected during the night/early morning and afternoon. Diel abundance of free viruses of 37-F6 in seawater roughly mirrored the transcriptional activities of both virus and host. In Osaka Bay, among viral relatives (genus level), virus 37-F6 specifically showed the highest ratio of transcriptional activity to virome abundance, a proxy for viral transcriptional activity relative to free viral particle abundance. This high ratio suggests high infection rate efficiencies in vSAG 37-F6 virus compared to viral relatives. Thus, time-series data revealed temporal transcript activities in one of the most abundant viruses in Earth.

Identification and Characterization of the First Virulent Phages, Including a Novel Jumbo Virus, Infecting Ochrobactrum spp

The Ochrobactrum genus consists of an extensive repertoire of biotechnologically valuable bacterial strains but also opportunistic pathogens. In our previous study, a novel strain, Ochrobactrum sp. POC9, which enhances biogas production in wastewater treatment plants (WWTPs) was identified and thoroughly characterized. Despite an insightful analysis of that bacterium, its susceptibility to bacteriophages present in WWTPs has not been evaluated. Using raw sewage sample from WWTP and applying the enrichment method, two virulent phages, vB_OspM_OC and vB_OspP_OH, which infect the POC9 strain, were isolated. These are the first virulent phages infecting Ochrobactrum spp. identified so far. Both phages were subjected to thorough functional and genomic analyses, which allowed classification of the vB_OspM_OC virus as a novel jumbo phage, with a genome size of over 227 kb. This phage encodes DNA methyltransferase, which mimics the specificity of cell cycle regulated CcrM methylase, a component of the epigenetic regulatory circuits in Alphaproteobacteria. In this study, an analysis of the overall diversity of Ochrobactrum-specific (pro)phages retrieved from databases and extracted in silico from bacterial genomes was also performed. Complex genome mining allowed us to build similarity networks to compare 281 Ochrobactrum-specific viruses. Analyses of the obtained networks revealed a high diversity of Ochrobactrum phages and their dissimilarity to the viruses infecting other bacteria.

Iroki: automatic customization and visualization of phylogenetic trees

Phylogenetic trees are an important analytical tool for evaluating community diversity and evolutionary history. In the case of microorganisms, the decreasing cost of sequencing has enabled researchers to generate ever-larger sequence datasets, which in turn have begun to fill gaps in the evolutionary history of microbial groups. However, phylogenetic analyses of these types of datasets create complex trees that can be challenging to interpret. Scientific inferences made by visual inspection of phylogenetic trees can be simplified and enhanced by customizing various parts of the tree. Yet, manual customization is time-consuming and error prone, and programs designed to assist in batch tree customization often require programming experience or complicated file formats for annotation. Iroki, a user-friendly web interface for tree visualization, addresses these issues by providing automatic customization of large trees based on metadata contained in tab-separated text files. Iroki’s utility for exploring biological and ecological trends in sequencing data was demonstrated through a variety of microbial ecology applications in which trees with hundreds to thousands of leaf nodes were customized according to extensive collections of metadata. The Iroki web application and documentation are available at https://www.iroki.net or through the VIROME portal http://virome.dbi.udel.edu. Iroki’s source code is released under the MIT license and is available at https://github.com/mooreryan/iroki.

Assembly-free single-molecule sequencing recovers complete virus genomes from natural microbial communities

Viruses are the most abundant biological entities on Earth and play key roles in host ecology, evolution, and horizontal gene transfer. Despite recent progress in viral metagenomics, the inherent genetic complexity of virus populations still poses technical difficulties for recovering complete virus genomes from natural assemblages. To address these challenges, we developed an assembly-free, single-molecule nanopore sequencing approach, enabling direct recovery of complete virus genome sequences from environmental samples. Our method yielded thousands of full-length, high-quality draft virus genome sequences that were not recovered using standard short-read assembly approaches. Additionally, our analyses discriminated between populations whose genomes had identical direct terminal repeats versus those with circularly permuted repeats at their termini, thus providing new insight into native virus reproduction and genome packaging. Novel DNA sequences were discovered, whose repeat structures, gene contents, and concatemer lengths suggest they are phage-inducible chromosomal islands, which are packaged as concatemers in phage particles, with lengths that match the size ranges of co-occurring phage genomes. Our new virus sequencing strategy can provide previously unavailable information about the genome structures, population biology, and ecology of naturally occurring viruses and viral parasites.

Metagenome Mining Reveals Hidden Genomic Diversity of Pelagimyophages in Aquatic Environments

SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

The SAR11 clade is one of the most abundant bacterioplankton groups in surface waters of most of the oceans and lakes. However, only 15 SAR11 phages have been isolated thus far, and only one of them belongs to the Myoviridae family (pelagimyophages). Here, we have analyzed 26 sequences of myophages that putatively infect the SAR11 clade. They have been retrieved by mining ca. 45 Gbp aquatic assembled cellular metagenomes and viromes. Most of the myophages were obtained from the cellular fraction (0.2 μm), indicating a bias against this type of virus in viromes. We have found the first myophages that putatively infect Candidatus Fonsibacter (freshwater SAR11) and another group putatively infecting bathypelagic SAR11 phylogroup Ic. The genomes have similar sizes and maintain overall synteny in spite of low average nucleotide identity values, revealing high similarity to marine cyanomyophages. Pelagimyophages recruited metagenomic reads widely from several locations but always much more from cellular metagenomes than from viromes, opposite to what happens with pelagipodophages. Comparing the genomes resulted in the identification of a hypervariable island that is related to host recognition. Interestingly, some genes in these islands could be related to host cell wall synthesis and coinfection avoidance. A cluster of curli-related proteins was widespread among the genomes, although its function is unclear.

IMPORTANCE SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

Double-stranded DNA virioplankton dynamics and reproductive strategies in the oligotrophic open ocean water column

Microbial communities are critical to ecosystem dynamics and biogeochemical cycling in the open oceans. Viruses are essential elements of these communities, influencing the productivity, diversity, and evolution of cellular hosts. To further explore the natural history and ecology of open-ocean viruses, we surveyed the spatiotemporal dynamics of double-stranded DNA (dsDNA) viruses in both virioplankton and bacterioplankton size fractions in the North Pacific Subtropical Gyre, one of the largest biomes on the planet. Assembly and clustering of viral genomes revealed a peak in virioplankton diversity at the base of the euphotic zone, where virus populations and host species richness both reached their maxima. Simultaneous characterization of both extracellular and intracellular viruses suggested depth-specific reproductive strategies. In particular, analyses indicated elevated lytic interactions in the mixed layer, more temporally variable temperate phage interactions at the base of the euphotic zone, and increased lysogeny in the mesopelagic ocean. Furthermore, the depth variability of auxiliary metabolic genes suggested habitat-specific strategies for viral influence on light-energy, nitrogen, and phosphorus acquisition during host infection. Most virus populations were temporally persistent over several years in this environment at the 95% nucleic acid identity level. In total, our analyses revealed variable distributional patterns and diverse reproductive and metabolic strategies of virus populations in the open-ocean water column.

Genomic and ecological attributes of marine bacteriophages encoding bacterial virulence genes

Background

Bacteriophages encode genes that modify bacterial functions during infection. The acquisition of phage-encoded virulence genes is a major mechanism for the rise of bacterial pathogens. In coral reefs, high bacterial density and lysogeny has been proposed to exacerbate reef decline through the transfer of phage-encoded virulence genes. However, the functions and distribution of these genes in phage virions on the reef remain unknown.

Results

Here, over 28,000 assembled viral genomes from the free viral community in Atlantic and Pacific Ocean coral reefs were queried against a curated database of virulence genes. The diversity of virulence genes encoded in the viral genomes was tested for relationships with host taxonomy and bacterial density in the environment. These analyses showed that bacterial density predicted the profile of virulence genes encoded by phages. The Shannon diversity of virulence-encoding phages was negatively related with bacterial density, leading to dominance of fewer genes at high bacterial abundances. A statistical learning analysis showed that reefs with high microbial density were enriched in viruses encoding genes enabling bacterial recognition and invasion of metazoan epithelium. Over 60% of phages could not have their hosts identified due to limitations of host prediction tools; for those which hosts were identified, host taxonomy was not an indicator of the presence of virulence genes.

Conclusions

This study described bacterial virulence factors encoded in the genomes of bacteriophages at the community level. The results showed that the increase in microbial densities that occurs during coral reef degradation is associated with a change in the genomic repertoire of bacteriophages, specifically in the diversity and distribution of bacterial virulence genes. This suggests that phages are implicated in the rise of pathogens in disturbed marine ecosystems.

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

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

Phage-specific metabolic reprogramming of virocells

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

Analysis of Spounaviruses as a Case Study for the Overdue Reclassification of Tailed Phages

Tailed bacteriophages are the most abundant and diverse viruses in the world, with genome sizes ranging from 10 kbp to over 500 kbp. Yet, due to historical reasons, all this diversity is confined to a single virus order-Caudovirales, composed of just four families: Myoviridae, Siphoviridae, Podoviridae, and the newly created Ackermannviridae family. In recent years, this morphology-based classification scheme has started to crumble under the constant flood of phage sequences, revealing that tailed phages are even more genetically diverse than once thought. This prompted us, the Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV), to consider overall reorganization of phage taxonomy. In this study, we used a wide range of complementary methods-including comparative genomics, core genome analysis, and marker gene phylogenetics-to show that the group of Bacillus phage SPO1-related viruses previously classified into the Spounavirinae subfamily, is clearly distinct from other members of the family Myoviridae and its diversity deserves the rank of an autonomous family. Thus, we removed this group from the Myoviridae family and created the family Herelleviridae-a new taxon of the same rank. In the process of the taxon evaluation, we explored the feasibility of different demarcation criteria and critically evaluated the usefulness of our methods for phage classification. The convergence of results, drawing a consistent and comprehensive picture of a new family with associated subfamilies, regardless of method, demonstrates that the tools applied here are particularly useful in phage taxonomy. We are convinced that creation of this novel family is a crucial milestone toward much-needed reclassification in the Caudovirales order.

Charting the diversity of uncultured viruses of Archaea and Bacteria

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Charting the diversity of uncultured viruses of Archaea and Bacteria

Background

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

Results

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

Conclusions

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

Prokaryotic Population Dynamics and Viral Predation in a Marine Succession Experiment Using Metagenomics

We performed an incubation experiment of seawater confined in plastic bottles with samples collected at three depths (15, 60, and 90 m) after retrieval from a single offshore location in the Mediterranean Sea, from a late summer stratified water column. Two samples representative of each depth were collected and stored in opaque bottles after two periods of 7 h. We took advantage of the "bottle effect" to investigate changes in the natural microbial communities (abundant and rare). We recovered 94 metagenome-assembled genomes (MAGs) and 1089 metagenomic viral contigs and examined their abundance using metagenomic recruitment. We detected a significant fast growth of copiotrophic bacteria such as Alteromonas or Erythrobacter throughout the entire water column with different dynamics that we assign to "clonal," "polyclonal," or "multispecies" depending on the recruitment pattern. Results also showed a marked ecotype succession in the phototropic picocyanobacteria that were able to grow at all the depths in the absence of light, highlighting the importance of their mixotrophic potential. In addition, "wall-chain-reaction" hypothesis based on the study of phage-host dynamics showed the higher impact of viral predation on archaea in deeper waters, evidencing their prominent role during incubations. Our results provide a step forward in understanding the mechanisms underlying dynamic patterns and ecology of the marine microbiome and the importance of processing the samples immediately after collection to avoid changes in the community structure.

Diverse, Abundant, and Novel Viruses Infecting the Marine Roseobacter RCA Lineage

Many major marine bacterial lineages such as SAR11, Prochlorococcus, SAR116, and several Roseobacter lineages have members that are abundant, relatively slow-growing, and genome streamlined. The isolation of phages that infect SAR11 and SAR116 have demonstrated the dominance of these phages in the marine virosphere. However, no phages have been isolated from bacteria in the Roseobacter RCA lineage, another abundant group of marine bacteria. In this study, seven RCA phages that infect three different RCA strains were isolated and characterized. All seven RCA phages belong to the Podoviridae family and have genome sizes ranging from 39.6 to 58.1 kb. Interestingly, three RCA phages (CRP-1, CRP-2, and CRP-3) show similar genomic content and architecture as SAR116 phage HMO-2011, which represents one of the most abundant known viral groups in the ocean. The high degree of homology among CRP-1, CRP-2, CRP-3, and HMO-2011 resulted in the contribution of RCA phages to the dominance of the HMO-2011-type group. CRP-4 and CRP-5 are similar to the Cobavirus group roseophages in terms of gene content and organization. The remaining two RCA phages, CRP-6 and CRP-7, show limited genomic similarity with known phages and represent two new phage groups. Metagenomic fragment recruitment analyses reveal that these RCA phage groups are much more abundant in the ocean than most existing marine roseophage groups. The characterization of these RCA phages has greatly expanded our understanding of the genomic diversity and evolution of marine roseophages and suggests the critical need for isolating phages from the abundant but "unculturable" bacteria.IMPORTANCE The RCA lineage of the marine Roseobacter group represents one of the slow-growing but dominant components of marine microbial communities. Although dozens of roseophages have been characterized, no phages infecting RCA strains have been reported. In this study, we reported on the first RCA phage genomes and investigated their distribution pattern and relative abundance in comparison with other important marine phage groups. Two of the four RCA phage groups were found closely related to previously reported SAR116 phage HMO-2011 and Cobavirus group roseophages, respectively. The remaining two groups are novel in the genome contents. Our study also revealed that RCA phages are widely distributed and exhibit high abundance in marine viromic data sets. Altogether, our findings have greatly broadened our understanding of RCA phages and emphasize the ecological and evolutionary importance of RCA phages in the marine virosphere.

Long-term stability and Red Queen-like strain dynamics in marine viruses

Viruses that infect microorganisms dominate marine microbial communities numerically, with impacts ranging from host evolution to global biogeochemical cycles^1,2. However, virus community dynamics, necessary for conceptual and mechanistic model development, remains difficult to assess. Here, we describe the long-term stability of a viral community by analysing the metagenomes of near-surface 0.02-0.2 μm samples from the San Pedro Ocean Time-series³ that were sampled monthly over 5 years. Of 19,907 assembled viral contigs (>5 kb, mean 15 kb), 97% were found in each sample (by >98% ID metagenomic read recruitment) to have relative abundances that ranged over seven orders of magnitude, with limited temporal reordering of rank abundances along with little change in richness. Seasonal variations in viral community composition were superimposed on the overall stability; maximum community similarity occurred at 12-month intervals. Despite the stability of viral genotypic clusters that had 98% sequence identity, viral sequences showed transient variations in single-nucleotide polymorphisms (SNPs) and constant turnover of minor population variants, each rising and falling over a few months, reminiscent of Red Queen dynamics⁴. The rise and fall of variants within populations, interpreted through the perspective of known virus-host interactions⁵, is consistent with the hypothesis that fluctuating selection acts on a microdiverse cloud of strains, and this succession is associated with ever-shifting virus-host defences and counterdefences. This results in long-term virus-host coexistence that is facilitated by perpetually changing minor variants.

Absolute quantification of infecting viral particles by chip-based digital polymerase chain reaction

In silico and empirical quantification of viruses is paramount for obtaining information on viral populations that have a major impact on biogeochemical cycles. The uncultured Pelagibacter virus vSAG 37-F6 discovered via single-virus genomics is one of the most abundant and cosmopolitan marine viruses; however, little is understood about its temporal variation. Here, we estimated the absolute number of infecting 37-F6 viruses in coastal bacterioplankton from the Mediterranean Sea by using a novel, feasible SYBR Green I chip-based digital PCR (SYBR dPCR) technique, not implemented before for enumerating (uncultured) microbes. Quantitative SYBR dPCR estimated 450-3480 genome copies of virus 37-F6 in cells/mL (i.e. infecting viruses) and a total of ≈10-400 putative infected cells/mL with a potential C release of 0.12-4.9 pg/ml in the analysed samples. Considering that virus 37-F6 is ubiquitous and abundant in all Tara samples, an enormous amount of C could be transformed by this virus through the 'viral shunt'. Thus, this SYBR dPCR technique has enabled the absolute quantification of an ecologically relevant uncultured virus in nature and the estimation of its potential contribution on biogeochemical cycles. Overall, our study also shows that this approach has a broad applicability for quantifying any other target loci in Microbiology and Virology.

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

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

Trends of Microdiversity Reveal Depth-Dependent Evolutionary Strategies of Viruses in the Mediterranean

Viruses are extremely abundant and diverse biological entities that contribute to the functioning of marine ecosystems. Despite their recognized importance, few studies have addressed trends of mutation accumulation in marine viral communities across depth gradients. By investigating these trends, we show that mutation frequencies differ among viral genes according to their molecular functions, with the highest microdiversity occurring among proteins related to host metabolism, followed by structural proteins and, lastly, genome replication proteins. This is in agreement with evolutionary theory that postulates that housekeeping genes are under strong purifying selection. We also observed a positive association between depth and microdiversity. One exception to this trend was the host recognition proteins from the deep chlorophyll maximum, which displayed strikingly high microdiversity, which we hypothesize to be associated with intraspecies competition for hosts. Finally, our data allowed us to propose a theoretical model for viral microdiversity across the depth gradient. These discoveries are of special relevance because many of the viral genomic sequences discovered here were predicted to infect some of the most abundant bacteria in marine ecosystems, such as “Candidatus Pelagibacter,” Puniceispirillum, and Prochlorococcus.

The evolutionary interactions between viruses and their prokaryotic hosts remain a little-known aspect of microbial evolution. Most studies on this topic were carried out in pure cultures that challenge one virus with one bacterial clone at a time, which is very removed from real-life situations. Few studies have addressed trends of microdiversity in marine viral communities throughout depth gradients. We analyzed metagenomes from both the cellular and viral fractions of Mediterranean seawater samples spanning the epipelagic to the bathypelagic zones at depths of 15, 45, 60, and 2,000 m during the summer stratification of the water column. We evaluated microdiversity patterns by measuring the accumulation of synonymous and nonsynonymous mutations in viral genes. Our results demonstrated clear depth-dependent trends in the frequency of polymorphic sites and nonsynonymous mutations among genes encoding metabolic, structural, and replication proteins. These differences were linked to changes in energy availability, host and viral densities, and the proportions of actively replicating viruses. We propose the hypothesis that in the energy-rich, high-host-density, euphotic depths, selection acts to favor diversity of the host recognition machinery to increase host range, while in energy-depleted aphotic waters, selection acts on viral replication fitness, enhancing diversity in auxiliary metabolic genes.

IMPORTANCE Viruses are extremely abundant and diverse biological entities that contribute to the functioning of marine ecosystems. Despite their recognized importance, few studies have addressed trends of mutation accumulation in marine viral communities across depth gradients. By investigating these trends, we show that mutation frequencies differ among viral genes according to their molecular functions, with the highest microdiversity occurring among proteins related to host metabolism, followed by structural proteins and, lastly, genome replication proteins. This is in agreement with evolutionary theory that postulates that housekeeping genes are under strong purifying selection. We also observed a positive association between depth and microdiversity. One exception to this trend was the host recognition proteins from the deep chlorophyll maximum, which displayed strikingly high microdiversity, which we hypothesize to be associated with intraspecies competition for hosts. Finally, our data allowed us to propose a theoretical model for viral microdiversity across the depth gradient. These discoveries are of special relevance because many of the viral genomic sequences discovered here were predicted to infect some of the most abundant bacteria in marine ecosystems, such as “Candidatus Pelagibacter,” Puniceispirillum, and Prochlorococcus.

Wide Distribution of Phage That Infect Freshwater SAR11 Bacteria

Fonsibacter (LD12 subclade) is among the most abundant bacterioplankton in freshwater ecosystems. These bacteria belong to the order Pelagibacterales (SAR11) and are related to Pelagibacter (marine SAR11), which dominates many marine habitats. Although a few Pelagibacter phage (Pelagiphage) have been described, no phage that infect Fonsibacter have been reported. In this study, we describe two groups of Podoviridae phage that infect Fonsibacter A complete Fonsibacter genome containing a prophage was reconstructed from metagenomic data. A circularized and complete genome related to the prophage, referred to as uv-Fonsiphage-EPL (lysogenic strategy), shows high similarity to marine Pelagiphage HTVC025P. Additionally, we reconstructed three complete genomes and one draft genome of phage related to marine Pelagiphage HTVC010P and predicted a lytic strategy. The similarity in codon usage and cooccurrence patterns of HTVC010P-related phage and Fonsibacter suggested that these phage infect Fonsibacter Similar phage were detected in Lake Mendota, Wisconsin, where Fonsibacter is also present. A search of related phage revealed the worldwide distribution of some genotypes in freshwater ecosystems, suggesting their substantial role in shaping indigenous microbial assemblages and influence on biogeochemical cycling. However, the uv-Fonsiphage-EPL and one group of HTVC010P-related phage have a more limited distribution in freshwater ecosystems. Overall, the findings provide insights into the genomic features of phage that infect Fonsibacter and expand understanding of the ecology and evolution of these important bacteria.IMPORTANCE Fonsibacter represents a significant microbial group of freshwater ecosystems. Although the genomic and metabolic features of these bacteria have been well studied, no phage infecting them has been reported. In this study, we reconstructed complete genomes of Fonsibacter and infecting phage and revealed their close relatedness to the phage infecting marine SAR11 members. Also, we illustrated that phage that infect Fonsibacter are widely distributed in freshwater habitats. In summary, the results contribute new insights into the ecology and evolution of Fonsibacter and phage.

Wide Distribution of Phage That Infect Freshwater SAR11 Bacteria

Fonsibacter represents a significant microbial group of freshwater ecosystems. Although the genomic and metabolic features of these bacteria have been well studied, no phage infecting them has been reported. In this study, we reconstructed complete genomes of Fonsibacter and infecting phage and revealed their close relatedness to the phage infecting marine SAR11 members. Also, we illustrated that phage that infect Fonsibacter are widely distributed in freshwater habitats. In summary, the results contribute new insights into the ecology and evolution of Fonsibacter and phage.

Fonsibacter (LD12 subclade) is among the most abundant bacterioplankton in freshwater ecosystems. These bacteria belong to the order Pelagibacterales (SAR11) and are related to Pelagibacter (marine SAR11), which dominates many marine habitats. Although a few Pelagibacter phage (Pelagiphage) have been described, no phage that infect Fonsibacter have been reported. In this study, we describe two groups of Podoviridae phage that infect Fonsibacter. A complete Fonsibacter genome containing a prophage was reconstructed from metagenomic data. A circularized and complete genome related to the prophage, referred to as uv-Fonsiphage-EPL (lysogenic strategy), shows high similarity to marine Pelagiphage HTVC025P. Additionally, we reconstructed three complete genomes and one draft genome of phage related to marine Pelagiphage HTVC010P and predicted a lytic strategy. The similarity in codon usage and cooccurrence patterns of HTVC010P-related phage and Fonsibacter suggested that these phage infect Fonsibacter. Similar phage were detected in Lake Mendota, Wisconsin, where Fonsibacter is also present. A search of related phage revealed the worldwide distribution of some genotypes in freshwater ecosystems, suggesting their substantial role in shaping indigenous microbial assemblages and influence on biogeochemical cycling. However, the uv-Fonsiphage-EPL and one group of HTVC010P-related phage have a more limited distribution in freshwater ecosystems. Overall, the findings provide insights into the genomic features of phage that infect Fonsibacter and expand understanding of the ecology and evolution of these important bacteria.

IMPORTANCEFonsibacter represents a significant microbial group of freshwater ecosystems. Although the genomic and metabolic features of these bacteria have been well studied, no phage infecting them has been reported. In this study, we reconstructed complete genomes of Fonsibacter and infecting phage and revealed their close relatedness to the phage infecting marine SAR11 members. Also, we illustrated that phage that infect Fonsibacter are widely distributed in freshwater habitats. In summary, the results contribute new insights into the ecology and evolution of Fonsibacter and phage.

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

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

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

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

A Review on Viral Metagenomics in Extreme Environments

Viruses are the most abundant biological entities in the biosphere, and have the ability to infect Bacteria, Archaea, and Eukaryotes. The virome is estimated to be at least ten times more abundant than the microbiome with 10⁷ viruses per milliliter and 10⁹ viral particles per gram in marine waters and sediments or soils, respectively. Viruses represent a largely unexplored genetic diversity, having an important role in the genomic plasticity of their hosts. Moreover, they also play a significant role in the dynamics of microbial populations. In recent years, metagenomic approaches have gained increasing popularity in the study of environmental viromes, offering the possibility of extending our knowledge related to both virus diversity and their functional characterization. Extreme environments represent an interesting source of both microbiota and their virome due to their particular physicochemical conditions, such as very high or very low temperatures and >1 atm hydrostatic pressures, among others. Despite the fact that some progress has been made in our understanding of the ecology of the microbiota in these habitats, few metagenomic studies have described the viromes present in extreme ecosystems. Thus, limited advances have been made in our understanding of the virus community structure in extremophilic ecosystems, as well as in their biotechnological potential. In this review, we critically analyze recent progress in metagenomic based approaches to explore the viromes in extreme environments and we discuss the potential for new discoveries, as well as methodological challenges and perspectives.

Cooccurrence of Broad- and Narrow-Host-Range Viruses Infecting the Bloom-Forming Toxic Cyanobacterium Microcystis aeruginosa

Viruses play important roles in regulating the abundance and composition of bacterial populations in aquatic ecosystems. The bloom-forming toxic cyanobacterium Microcystis aeruginosa is predicted to interact with diverse cyanoviruses, resulting in Microcystis population diversification. However, current knowledge of the genomes from these viruses and their infection programs is limited to those of Microcystis virus Ma-LMM01. Here, we performed a time series sampling at a small pond in Japan during a Microcystis bloom and then investigated the genomic information and transcriptional dynamics of Microcystis-interacting viruses using metagenomic and metatranscriptomic approaches. We identified 15 viral genomic fragments classified into three groups, groups I (including Ma-LMM01), II (high abundance and transcriptional activity), and III (new lineages). According to the phylogenetic distribution of Microcystis strains possessing spacers against each viral group, the group II-original viruses interacted with all three phylogenetically distinct Microcystis population types (phylotypes), whereas the groups I and III-original viruses interacted with only one or two phylotypes, indicating the cooccurrence of broad- (group II) and narrow (groups I and III)-host-range viruses in the bloom. These viral fragments showed the highest transcriptional levels during daytime regardless of their genomic differences. Interestingly, M. aeruginosa expressed antiviral defense genes in the environment, unlike what was seen with an Ma-LMM01 infection in a previous culture experiment. Given that broad-host-range viruses often induce antiviral responses within alternative hosts, our findings suggest that such antiviral responses might inhibit viral multiplication, mainly that of broad-host-range viruses like those in group II.IMPORTANCE The bloom-forming toxic cyanobacterium Microcystis aeruginosa is thought to have diversified its population through the interactions between host and viruses in antiviral defense systems. However, current knowledge of viral genomes and infection programs is limited to those of Microcystis virus Ma-LMM01, which was a narrow host range in which it can escape from the highly abundant host defense systems. Our metagenomic approaches unveiled the cooccurrence of narrow- and broad-host-range Microcystis viruses, which included fifteen viral genomic fragments from Microcystis blooms that were classified into three groups. Interestingly, Microcystis antiviral defense genes were expressed against viral infection in the environment, unlike what was seen in a culture experiment with Ma-LMM01. Given that viruses with a broad host range often induce antiviral responses within alternative hosts, our findings suggest that antiviral responses inhibit viral reproduction, especially that of broad-range viruses like those in group II. This paper augments our understanding of the interactions between M. aeruginosa and its viruses and fills an important knowledge gap.

Long-read metagenomic exploration of extrachromosomal mobile genetic elements in the human gut

BACKGROUND

Elucidating the ecological and biological identity of extrachromosomal mobile genetic elements (eMGEs), such as plasmids and bacteriophages, in the human gut remains challenging due to their high complexity and diversity.

RESULTS

Here, we show efficient identification of eMGEs as complete circular or linear contigs from PacBio long-read metagenomic data. De novo assembly of PacBio long reads from 12 faecal samples generated 82 eMGE contigs (2.5~666.7-kb), which were classified as 71 plasmids and 11 bacteriophages, including 58 novel plasmids and six bacteriophages, and complete genomes of five diverse crAssphages with terminal direct repeats. In a dataset of 413 gut metagenomes from five countries, many of the identified plasmids were highly abundant and prevalent. The ratio of gut plasmids by our plasmid data is more than twice that in the public database. Plasmids outnumbered bacterial chromosomes three to one on average in this metagenomic dataset. Host prediction suggested that Bacteroidetes-associated plasmids predominated, regardless of microbial abundance. The analysis found several plasmid-enriched functions, such as inorganic ion transport, while antibiotic resistance genes were harboured mostly in low-abundance Proteobacteria-associated plasmids.

CONCLUSIONS

Overall, long-read metagenomics provided an efficient approach for unravelling the complete structure of human gut eMGEs, particularly plasmids.

A novel uncultured marine cyanophage lineage with lysogenic potential linked to a putative marine Synechococcus 'relic' prophage

Marine cyanobacteria are important contributors to primary production in the ocean and their viruses (cyanophages) affect the ocean microbial communities. Despite reports of lysogeny in marine cyanobacteria, a genome sequence of such temperate cyanophages remains unknown although genomic analysis indicate potential for lysogeny in certain marine cyanophages. Using assemblies from Red Sea and Tara Oceans metagenomes, we recovered genomes of a novel uncultured marine cyanophage lineage, which contain, in addition to common cyanophage genes, a phycobilisome degradation protein NblA, an integrase and a split DNA polymerase. The DNA polymerase forms a monophyletic clade with a DNA polymerase from a genomic island in Synechococcus WH8016. The island contains a relic prophage that does not resemble any previously reported cyanophage but shares several genes with the newly identified cyanophages reported here. Metagenomic recruitment indicates that the novel cyanophages are widespread, albeit at low abundance. Here, we describe a novel potentially lysogenic cyanophage family, their abundance and distribution in the marine environment.

Droplet Digital PCR for Estimating Absolute Abundances of Widespread Pelagibacter Viruses

Absolute abundances of prokaryotes are typically determined by FISH. Due to the lack of a universal conserved gene among all viruses, metagenomic fragment recruitment is commonly used to estimate the relative viral abundance. However, the paucity of absolute virus abundance data hinders our ability to fully understand how viruses drive global microbial populations. The cosmopolitan marine Pelagibacter ubique is host for the highly widespread HTVC010P pelagiphage isolate and the extremely abundant uncultured virus vSAG 37-F6 recently discovered by single-virus genomics. Here we applied droplet digital PCR (ddPCR) to calculate the absolute abundance of these pelagiphage genotypes in the Mediterranean Sea and the Gulf of Maine. Abundances were between 360 and 8,510 virus mL-1 and 1,270–14,400 virus mL-1 for vSAG 37-F6 and HTVC010P, respectively. Illumina PCR-amplicon sequencing corroborated the absence of ddPCR non-specific amplifications for vSAG 37-F6, but showed an overestimation of 6% for HTVC010P from off-targets, genetically unrelated viruses. Absolute abundances of both pelagiphages, two of the most abundance marine viruses, suggest a large viral pelagiphage diversity in marine environments, and show the efficiency and power of ddPCR to disentangle the structure of marine viral communities. Results also highlight the need for a standardized workflow to obtain accurate quantification that allows cross data comparison.

TAR-VIR: a pipeline for TARgeted VIRal strain reconstruction from metagenomic data

Background

Strain-level RNA virus characterization is essential for developing prevention and treatment strategies. Viral metagenomic data, which can contain sequences of both known and novel viruses, provide new opportunities for characterizing RNA viruses. Although there are a number of pipelines for analyzing viruses in metagenomic data, they have different limitations. First, viruses that lack closely related reference genomes cannot be detected with high sensitivity. Second, strain-level analysis is usually missing.

Results

In this study, we developed a hybrid pipeline named TAR-VIR that reconstructs viral strains without relying on complete or high-quality reference genomes. It is optimized for identifying RNA viruses from metagenomic data by combining an effective read classification method and our in-house strain-level de novo assembly tool. TAR-VIR was tested on both simulated and real viral metagenomic data sets. The results demonstrated that TAR-VIR competes favorably with other tested tools.

Conclusion

TAR-VIR can be used standalone for viral strain reconstruction from metagenomic data. Or, its read recruiting stage can be used with other de novo assembly tools for superior viral functional and taxonomic analyses. The source code and the documentation of TAR-VIR are available at https://github.com/chjiao/TAR-VIR.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2878-2) contains supplementary material, which is available to authorized users.

Long-read viral metagenomics captures abundant and microdiverse viral populations and their niche-defining genomic islands

Marine viruses impact global biogeochemical cycles via their influence on host community structure and function, yet our understanding of viral ecology is constrained by limitations in host culturing and a lack of reference genomes and 'universal' gene markers to facilitate community surveys. Short-read viral metagenomic studies have provided clues to viral function and first estimates of global viral gene abundance and distribution, but their assemblies are confounded by populations with high levels of strain evenness and nucleotide diversity (microdiversity), limiting assembly of some of the most abundant viruses on Earth. Such features also challenge assembly across genomic islands containing niche-defining genes that drive ecological speciation. These populations and features may be successfully captured by single-virus genomics and fosmid-based approaches, at least in abundant taxa, but at considerable cost and technical expertise. Here we established a low-cost, low-input, high throughput alternative sequencing and informatics workflow to improve viral metagenomic assemblies using short-read and long-read technology. The 'VirION' (Viral, long-read metagenomics via MinION sequencing) approach was first validated using mock communities where it was found to be as relatively quantitative as short-read methods and provided significant improvements in recovery of viral genomes. We then then applied VirION to the first metagenome from a natural viral community from the Western English Channel. In comparison to a short-read only approach, VirION: (i) increased number and completeness of assembled viral genomes; (ii) captured abundant, highly microdiverse virus populations, and (iii) captured more and longer genomic islands. Together, these findings suggest that VirION provides a high throughput and cost-effective alternative to fosmid and single-virus genomic approaches to more comprehensively explore viral communities in nature.

Host Taxon Predictor - A Tool for Predicting Taxon of the Host of a Newly Discovered Virus

Recent advances in metagenomics provided a valuable alternative to culture-based approaches for better sampling viral diversity. However, some of newly identified viruses lack sequence similarity to any of previously sequenced ones, and cannot be easily assigned to their hosts. Here we present a bioinformatic approach to this problem. We developed classifiers capable of distinguishing eukaryotic viruses from the phages achieving almost 95% prediction accuracy. The classifiers are wrapped in Host Taxon Predictor (HTP) software written in Python which is freely available at https://github.com/wojciech-galan/viruses_classifier. HTP’s performance was later demonstrated on a collection of newly identified viral genomes and genome fragments. In summary, HTP is a culture- and alignment-free approach for distinction between phages and eukaryotic viruses. We have also shown that it is possible to further extend our method to go up the evolutionary tree and predict whether a virus can infect narrower taxa.

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.

Family A DNA Polymerase Phylogeny Uncovers Diversity and Replication Gene Organization in the Virioplankton

Shotgun metagenomics, which allows for broad sampling of viral diversity, has uncovered genes that are widely distributed among virioplankton populations and show linkages to important biological features of unknown viruses. Over 25% of known dsDNA phage carry the DNA polymerase I (polA) gene, making it one of the most widely distributed phage genes. Because of its pivotal role in DNA replication, this enzyme is linked to phage lifecycle characteristics. Previous research has suggested that a single amino acid substitution might be predictive of viral lifestyle. In this study Chesapeake Bay virioplankton were sampled by shotgun metagenomic sequencing (using long and short read technologies). More polA sequences were predicted from this single viral metagenome (virome) than from 86 globally distributed virome libraries (ca. 2,100, and 1,200, respectively). The PolA peptides predicted from the Chesapeake Bay virome clustered with 69% of PolA peptides from global viromes; thus, remarkably the Chesapeake Bay virome captured the majority of known PolA peptide diversity in viruses. This deeply sequenced virome also expanded the diversity of PolA sequences, increasing the number of PolA clusters by 44%. Contigs containing polA sequences were also used to examine relationships between phylogenetic clades of PolA and other genes within unknown viral populations. Phylogenic analysis revealed five distinct groups of phages distinguished by the amino acids at their 762 (Escherichia coli IAI39 numbering) positions and replication genes. DNA polymerase I sequences from Tyr762 and Phe762 groups were most often neighbored by ring-shaped superfamily IV helicases and ribonucleotide reductases (RNRs). The Leu762 groups had non-ring shaped helicases from superfamily II and were further distinguished by an additional helicase gene from superfamily I and the lack of any identifiable RNR genes. Moreover, we found that the inclusion of ribonucleotide reductase associated with PolA helped to further differentiate phage diversity, chiefly within lytic podovirus populations. Altogether, these data show that DNA Polymerase I is a useful marker for observing the diversity and composition of the virioplankton and may be a driving factor in the divergence of phage replication components.

The North Sea goes viral: Occurrence and distribution of North Sea bacteriophages

Marine viruses are dominated by phages and have an enormous influence on microbial population dynamics, due to lysis and horizontal gene transfer. The aim of this study is to analyze the occurrence and diversity of phages in the North Sea, considering the virus-host interactions and biogeographic factors. The virus community of four sampling stations were described using virus metagenomics (viromes). The results show that the virus community was not evenly distributed throughout the North Sea. The dominant phage members were identified as unclassified phage group, followed by Caudovirales order. Myoviridae was the dominant phage family in the North Sea, which occurrence decreased from the coast to the open sea. In contrast, the occurrence of Podoviridae increased and the occurrence of Siphoviridae was low throughout the North Sea. The occurrence of other groups such as Phycodnaviridae decreased from the coast to the open sea. The coastal virus community was genetically more diverse than the open sea community. The influence of riverine inflow and currents, for instance the English Channel flow affects the genetic virus diversity with the community carrying genes from a variety of metabolic pathways and other functions. The present study offers the first insights in the virus community in the North Sea using viromes and shows the variation in virus diversity and the genetic information moved from coastal to open sea areas.

A New Freshwater Cyanosiphovirus Harboring Integrase

Pelagic cyanobacteria are key players in the functioning of aquatic ecosystems, and their viruses (cyanophages) potentially affect the abundance and composition of cyanobacterial communities. Yet, there are few well-described freshwater cyanophages relative to their marine counterparts, and in general, few cyanosiphoviruses (family Siphoviridae) have been characterized, limiting our understanding of the biology and the ecology of this prominent group of viruses. Here, we characterize S-LBS1, a freshwater siphovirus lytic to a phycoerythrin-rich Synechococcus isolate (Strain TCC793). S-LBS1 has a narrow host range, a burst size of ∼400 and a relatively long infecting step before cell lysis occurs. It has a dsDNA 34,641 bp genome with putative genes for structure, DNA packing, lysis, replication, host interactions, DNA repair and metabolism. S-LBS1 is similar in genome size, genome architecture, and gene content, to previously described marine siphoviruses also infecting PE-rich Synechococcus, e.g., S-CBS1 and S-CBS3. However, unlike other Synechococcus phages, S-LBS1 encodes an integrase, suggesting its ability to establish lysogenic relationships with its host. Sequence recruitment from viral metagenomic data showed that S-LBS1-like viruses are diversely present in a wide range of aquatic environments, emphasizing their potential importance in controlling and structuring Synechococcus populations. A comparative analysis with 16 available sequenced cyanosiphoviruses reveals the absence of core genes within the genomes, suggesting high degree of genetic variability in siphoviruses infecting cyanobacteria. It is likely that cyanosiphoviruses have evolved as distinct evolutionary lineages and that adaptive co-evolution occurred between these viruses and their hosts (i.e., Synechococcus, Prochlorococcus, Nodularia, and Acaryochloris), constituting an important driving force for such phage diversification.

Single-cell genomics uncover Pelagibacter as the putative host of the extremely abundant uncultured 37-F6 viral population in the ocean

The identification of relevant virus-host pairs that globally account for a large pool of carbon and nutrients in the ocean is paramount to build accurate ecological models. A previous work using single-virus genomics led to the discovery of the uncultured single-virus vSAG 37-F6, originally sorted from the Mediterranean Sea (Blanes Bay Microbial Observatory), that represents one of the most abundant dsDNA viral population in the marine surface virosphere. Here, from same sampling site, we report that a Pelagibacter single-cell contained a viral member of vSAG 37-F6 population, by means of PCR screening of sorted, genome-amplified single cells with vSAG 37-F6-specific primers and whole-genome sequencing. Furthermore, viruses from this population were also found in three other Pelagibacter single cells from the South Pacific and Atlantic oceans. These new uncultured pelagiphages were genetically different from the previously characterized pelagiphage isolates. Data showed that the uncultured vSAG 37-F6 population represents the Pelagibacter phages that inhabit the sunlit ocean better, and contains a vast unrecognized microdiversity.

Metagenomics Sheds Light on the Ecology of Marine Microbes and Their Viruses

Advances brought about by omics-based approaches have revolutionized our understanding of the diversity and ecological processes involving marine archaea, bacteria, and their viruses. This broad review discusses recent examples of how genomics, metagenomics, and ecogenomics have been applied to reveal the ecology of these biological entities. Three major topics are covered in this revision: (i) the novel roles of microorganisms in ecosystem processes; (ii) virus-host associations; and (iii) ecological associations of microeukaryotes and other microbes. We also briefly comment on the discovery of novel taxa from marine ecosystems; development of a robust taxonomic framework for prokaryotes; breakthroughs on the diversity and ecology of cyanobacteria; and advances on ecological modelling. We conclude by discussing limitations of the field and suggesting directions for future research.

Pseudomonas PB1-Like Phages: Whole Genomes from Metagenomes Offer Insight into an Abundant Group of Bacteriophages

Despite the abundance, ubiquity and impact of environmental viruses, their inherent genomic plasticity and extreme diversity pose significant challenges for the examination of bacteriophages on Earth. Viral metagenomic studies have offered insight into broader aspects of phage ecology and repeatedly uncover genes to which we are currently unable to assign function. A combined effort of phage isolation and metagenomic survey of Chicago&rsquo;s nearshore waters of Lake Michigan revealed the presence of Pbunaviruses, relatives of the Pseudomonas phage PB1. This prompted our expansive investigation of PB1-like phages. Genomic signatures of PB1-like phages and Pbunaviruses were identified, permitting the unambiguous distinction between the presence/absence of these phages in soils, freshwater and wastewater samples, as well as publicly available viral metagenomic datasets. This bioinformatic analysis led to the de novo assembly of nine novel PB1-like phage genomes from a metagenomic survey of samples collected from Lake Michigan. While this study finds that Pbunaviruses are abundant in various environments of Northern Illinois, genomic variation also exists to a considerable extent within individual communities.