Depth-stratified functional and taxonomic niche specialization in the 'core' and 'flexible' Pacific Ocean Virome

Viral Ecogenomics of Arctic Cryopeg Brine and Sea Ice

Arctic regions, which are changing rapidly as they warm 2 to 3 times faster than the global average, still retain microbial habitats that serve as natural laboratories for understanding mechanisms of microbial adaptation to extreme conditions. Seawater-derived brines within both sea ice (sea-ice brine) and ancient layers of permafrost (cryopeg brine) support diverse microbes adapted to subzero temperatures and high salinities, yet little is known about viruses in these extreme environments, which, if analogous to other systems, could play important evolutionary and ecosystem roles. Here, we characterized viral communities and their functions in samples of cryopeg brine, sea-ice brine, and melted sea ice. Viral abundance was high in cryopeg brine (1.2 × 108 ml-1) and much lower in sea-ice brine (1.3 × 105 to 2.1 × 105 ml-1), which roughly paralleled the differences in cell concentrations in these samples. Five low-input, quantitative viral metagenomes were sequenced to yield 476 viral populations (i.e., species level; ≥10 kb), only 12% of which could be assigned taxonomy by traditional database approaches, indicating a high degree of novelty. Additional analyses revealed that these viruses: (i) formed communities that differed between sample type and vertically with sea-ice depth; (ii) infected hosts that dominated these extreme ecosystems, including Marinobacter, Glaciecola, and Colwellia; and (iii) encoded fatty acid desaturase (FAD) genes that likely helped their hosts overcome cold and salt stress during infection, as well as mediated horizontal gene transfer of FAD genes between microbes. Together, these findings contribute to understanding viral abundances and communities and how viruses impact their microbial hosts in subzero brines and sea ice.IMPORTANCE This study explores viral community structure and function in remote and extreme Arctic environments, including subzero brines within marine layers of permafrost and sea ice, using a modern viral ecogenomics toolkit for the first time. In addition to providing foundational data sets for these climate-threatened habitats, we found evidence that the viruses had habitat specificity, infected dominant microbial hosts, encoded host-derived metabolic genes, and mediated horizontal gene transfer among hosts. These results advance our understanding of the virosphere and how viruses influence extreme ecosystems. More broadly, the evidence that virally mediated gene transfers may be limited by host range in these extreme habitats contributes to a mechanistic understanding of genetic exchange among microbes under stressful conditions in other systems.

Depth-related variability in viral communities in highly stratified sulfidic mine tailings

BACKGROUND

Recent studies have significantly expanded our knowledge of viral diversity and functions in the environment. Exploring the ecological relationships between viruses, hosts, and the environment is a crucial first step towards a deeper understanding of the complex and dynamic interplays among them.

RESULTS

Here, we obtained extensive 16S rRNA gene amplicon, metagenomics sequencing, and geochemical datasets from different depths of two highly stratified sulfidic mine tailings cores with steep geochemical gradients especially pH, and explored how variations in viral community composition and functions were coupled to the co-existing prokaryotic assemblages and the varying environmental conditions. Our data showed that many viruses in the mine tailings represented novel genera, based on gene-sharing networks. Siphoviridae, Podoviridae, and Myoviridae dominated the classified viruses in the surface tailings and deeper layers. Both viral richness and normalized coverage increased with depth in the tailings cores and were significantly correlated with geochemical properties, for example, pH. Viral richness was also coupled to prokaryotic richness (Pearson's r = 0.65, P = 0.032). The enrichment of prophages in the surface mine tailings suggested a preference of lysogenic viral lifestyle in more acidic conditions. Community-wide comparative analyses clearly showed that viruses in the surface tailings encoded genes mostly with unknown functions while viruses in the deeper layers contained genes mainly annotated as conventional functions related to metabolism and structure. Notably, significantly abundant assimilatory sulfate reduction genes were identified from the deeper tailings layers and they were widespread in viruses predicted to infect diverse bacterial phyla.

CONCLUSIONS

Overall, our results revealed a depth-related distribution of viral populations in the extreme and heterogeneous tailings system. The viruses may interact with diverse hosts and dynamic environmental conditions and likely play a role in the functioning of microbial community and modulate sulfur cycles in situ. Video Abstract.

Tara Oceans: towards global ocean ecosystems biology

A planetary-scale understanding of the ocean ecosystem, particularly in light of climate change, is crucial. Here, we review the work of Tara Oceans, an international, multidisciplinary project to assess the complexity of ocean life across comprehensive taxonomic and spatial scales. Using a modified sailing boat, the team sampled plankton at 210 globally distributed sites at depths down to 1,000 m. We describe publicly available resources of molecular, morphological and environmental data, and discuss how an ecosystems biology approach has expanded our understanding of plankton diversity and ecology in the ocean as a planetary, interconnected ecosystem. These efforts illustrate how global-scale concepts and data can help to integrate biological complexity into models and serve as a baseline for assessing ecosystem changes and the future habitability of our planet in the Anthropocene epoch.

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.

Investigating the Diversity of Marine Bacteriophage in Contrasting Water Masses Associated with the East Australian Current (EAC) System

Virus- and bacteriophage-induced mortality can have a significant impact on marine productivity and alter the flux of nutrients in marine microbial food-webs. Viral mediated horizontal gene transfer can also influence host fitness and community composition. However, there are very few studies of marine viral diversity in the Southern Hemisphere, which hampers our ability to fully understand the complex interplay of biotic and abiotic factors that shape microbial communities. We carried out the first genetic study of bacteriophage communities within a dynamic western boundary current (WBC) system, the east Australian current (EAC). Virus DNA sequences were extracted from 63 assembled metagenomes and six metaviromes obtained from various depths at 24 different locations. More than 1700 bacteriophage genomic fragments (>9 kbps) were recovered from the assembled sequences. Bacteriophage diversity displayed distinct depth and regional patterns. There were clear differences in the bacteriophage populations associated with the EAC and Tasman Sea euphotic zones, at both the taxonomic and functional level. In contrast, bathypelagic phages were similar across the two oceanic regions. These data provide the first characterisation of viral diversity across a dynamic western boundary current, which is an emerging model for studying the response of microbial communities to climate change.

Phage diversity, genomics and phylogeny

Recent advances in viral metagenomics have enabled the rapid discovery of an unprecedented catalogue of phages in numerous environments, from the human gut to the deep ocean. Although these advances have expanded our understanding of phage genomic diversity, they also revealed that we have only scratched the surface in the discovery of novel viruses. Yet, despite the remarkable diversity of phages at the nucleotide sequence level, the structural proteins that form viral particles show strong similarities and conservation. Phages are uniquely interconnected from an evolutionary perspective and undergo multiple events of genetic exchange in response to the selective pressure of their hosts, which drives their diversity. In this Review, we explore phage diversity at the structural, genomic and community levels as well as the complex evolutionary relationships between phages, moulded by the mosaicity of their genomes.

Gut Bacteriophage: Current Understanding and Challenges

The gut microbiome is widely accepted to have a significant impact on human health yet, despite years of research on this complex ecosystem, the contributions of different forces driving microbial population structure remain to be fully elucidated. The viral component of the human gut microbiome is dominated by bacteriophage, which are known to play crucial roles in shaping microbial composition, driving bacterial diversity, and facilitating horizontal gene transfer. Bacteriophage are also one of the most poorly understood components of the human gut microbiome, with the vast majority of viral sequences sharing little to no homology to reference databases. If we are to understand the dynamics of bacteriophage populations, their interaction with the human microbiome and ultimately their influence on human health, we will depend heavily on sequence based approaches and in silico tools. This is complicated by the fact that, as with any research field in its infancy, methods of analyses vary and this can impede our ability to compare the outputs of different studies. Here, we discuss the major findings to date regarding the human virome and reflect on our current understanding of how gut bacteriophage shape the microbiome. We consider whether or not the virome field is built on unstable foundations and if so, how can we provide a solid basis for future experimentation. The virome is a challenging yet crucial piece of the human microbiome puzzle. In order to develop our understanding, we will discuss the need to underpin future studies with robust research methods and suggest some solutions to existing challenges.

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.

Fungal Community Composition and Potential Depth-Related Driving Factors Impacting Distribution Pattern and Trophic Modes from Epi- to Abyssopelagic Zones of the Western Pacific Ocean

Fungi play an important role in cycling organic matter and nutrients in marine ecosystems. However, the distribution of fungal communities in the ocean, especially the vertical distribution along depth in the water column, remained poorly understood. Here, we assess the depth-related distribution pattern of fungal communities along the water column from epi- to abyssopelagic zones of the Western Pacific Ocean using internal transcribed spacer 2 (ITS2) metabarcoding. Majority of the assigned OTUs were affiliated to Ascomycota, followed by three other minor phyla (Basidiomycota, Chytridiomycota, and Mucoromycota). The epipelagic zone harbored a higher OTU richness with distinct fungal communities as compared with meso-, bathy-, and abyssopelagic zones. Across the whole water column, depth appears as a key parameter for both fungal α- and β-diversity. However, when the dataset was split into the upper (5-500 m) and deeper (below 500 m) layers, no significant correlation was observed between depth and community compositions. In the upper layer, temperature and dissolved oxygen were recognized as the primary environmental factors shaping fungal α- and β- diversity. By parsing fungal OTUs into ecological categories, multi-trophic mode of nutrition was found to be more prevalent with increasing depth, suggesting a potential adaptation to the extreme conditions of the deep sea. This study provides new and meaningful information on the depth-stratified fungal diversity, community structure, and putative ecological roles in the open sea.

Potential Interactions between Clade SUP05 Sulfur-Oxidizing Bacteria and Phages in Hydrothermal Vent Sponges

In deep-sea hydrothermal vent environments, sulfur-oxidizing bacteria belonging to the clade SUP05 are crucial symbionts of invertebrate animals. Marine viruses, as the most abundant biological entities in the ocean, play essential roles in regulating the sulfur metabolism of the SUP05 bacteria. To date, vent sponge-associated SUP05 and their phages have not been well documented. The current study analyzed microbiomes of Haplosclerida sponges from hydrothermal vents in the Okinawa Trough and recovered the dominant SUP05 genome, designated VS-SUP05. Phylogenetic analysis showed that VS-SUP05 was closely related to endosymbiotic SUP05 strains from mussels living in deep-sea hydrothermal vent fields. Homology and metabolic pathway comparisons against free-living and symbiotic SUP05 strains revealed that the VS-SUP05 genome shared many features with the deep-sea mussel symbionts. Supporting a potentially symbiotic lifestyle, the VS-SUP05 genome contained genes involved in the synthesis of essential amino acids and cofactors that are desired by the host. Analysis of sponge-associated viral sequences revealed putative VS-SUP05 phages, all of which were double-stranded viruses belonging to the families Myoviridae, Siphoviridae, Podoviridae, and Microviridae Among the phage sequences, one contig contained metabolic genes (iscR, iscS, and iscU) involved in iron-sulfur cluster formation. Interestingly, genome sequence comparison revealed horizontal transfer of the iscS gene among phages, VS-SUP05, and other symbiotic SUP05 strains, indicating an interaction between marine phages and SUP05 symbionts. Overall, our findings confirm the presence of SUP05 bacteria and their phages in sponges from deep-sea vents and imply a beneficial interaction that allows adaptation of the host sponge to the hydrothermal vent environment.IMPORTANCE Chemosynthetic SUP05 bacteria dominate the microbial communities of deep-sea hydrothermal vents around the world, SUP05 bacteria utilize reduced chemical compounds in vent fluids and commonly form symbioses with invertebrate organisms. This symbiotic relationship could be key to adapting to such unique and extreme environments. Viruses are the most abundant biological entities on the planet and have been identified in hydrothermal vent environments. However, their interactions with the symbiotic microbes of the SUP05 clade, along with their role in the symbiotic system, remain unclear. Here, using metagenomic sequence-based analyses, we determined that bacteriophages may support metabolism in SUP05 bacteria and play a role in the sponge-associated symbiosis system in hydrothermal vent environments.

A Phage Protein Aids Bacterial Symbionts in Eukaryote Immune Evasion

Phages are increasingly recognized as important members of host-associated microbiomes, with a vast genomic diversity. The new frontier is to understand how phages may affect higher order processes, such as in the context of host-microbe interactions. Here, we use marine sponges as a model to investigate the interplay between phages, bacterial symbionts, and eukaryotic hosts. Using viral metagenomics, we find that sponges, although massively filtering seawater, harbor species-specific and even individually unique viral signatures that are taxonomically distinct from other environments. We further discover a symbiont phage-encoded ankyrin-domain-containing protein, which is widely spread in phages of many host-associated contexts including human. We confirm in macrophage infection assays that the ankyrin protein (ANKp) modulates the eukaryotic host immune response against bacteria. We predict that the role of ANKp in nature is to facilitate coexistence in the tripartite interplay between phages, symbionts, and sponges and possibly many other host-microbe associations.

Functional Diversity: An Epistemic Roadmap

Functional diversity holds the promise of understanding ecosystems in ways unattainable by taxonomic diversity studies. Underlying this promise is the intuition that investigating the diversity of what organisms actually do (i.e., their functional traits) within ecosystems will generate more reliable insights into the ways these ecosystems behave, compared to considering only species diversity. But this promise also rests on several conceptual and methodological (i.e., epistemic) assumptions that cut across various theories and domains of ecology. These assumptions should be clearly addressed, notably for the sake of an effective comparison and integration across domains, and for assessing whether or not to use functional diversity approaches for developing ecological management strategies. The objective of this contribution is to identify and critically analyze the most salient of these assumptions. To this aim, we provide an epistemic roadmap that pinpoints these assumptions along a set of historical, conceptual, empirical, theoretical, and normative dimensions.

Viral Communities in the Global Deep Ocean Conveyor Belt Assessed by Targeted Viromics

Viruses are an abundant, diverse and dynamic component of marine and terrestrial ecosystems. In the ocean, viruses play a key role in the biogeochemical cycles and controlling microbial abundance, diversity and evolution. Recent metagenomics studies assessed the structure of the viral community in the upper ocean. However, little is known about the compositional changes in viral communities along the deep ocean conveyor belt. To assess potential changes in the viral community in the global deep-water circulation system, water samples were collected in the core of the North Atlantic Deep Water (NADW) (∼2,500 m) and Pacific Antarctic Bottom Water (∼4,000 m). Microbial and viral abundance were evaluated by flow cytometry. Subsequently, flow cytometry was used to sort virus-like particles and next generation sequencing was applied to build DNA libraries from the sorted virus populations. The viral communities were highly diverse across different oceanic regions with high dissimilarity between samples. Only 18% of the viral protein clusters were shared between the NADW and the Pacific Antarctic Bottom Water. Few viral groups, mainly associated with uncultured environmental and uncultured Mediterranean viruses were ubiquitously distributed along the global deep-water circulation system. Thus, our results point to a few groups of widely distributed abundant viruses in addition to the presence of rare and diverse types of viruses at a local scale.

Diversities and potential biogeochemical impacts of mangrove soil viruses

BACKGROUND

Mangroves are ecologically and economically important forests of the tropics. As one of the most carbon-rich biomes, mangroves account for 11% of the total input of terrestrial carbon into oceans. Although viruses are considered to significantly influence local and global biogeochemical cycles, little information is available regarding the community structure, genetic diversity and ecological roles of viruses in mangrove ecosystems.

METHODS

Here, we utilised viral metagenomics sequencing and virome-specific bioinformatics tools to study viral communities in six mangrove soil samples collected from different mangrove habitats in Southern China.

RESULTS

Mangrove soil viruses were found to be largely uncharacterised. Phylogenetic analyses of the major viral groups demonstrated extensive diversity and previously unknown viral clades and suggested that global mangrove viral communities possibly comprise evolutionarily close genotypes. Comparative analysis of viral genotypes revealed that mangrove soil viromes are mainly affected by marine waters, with less influence coming from freshwaters. Notably, we identified abundant auxiliary carbohydrate-active enzyme (CAZyme) genes from mangrove viruses, most of which participate in biolysis of complex polysaccharides, which are abundant in mangrove soils and organism debris. Host prediction results showed that viral CAZyme genes are diverse and probably widespread in mangrove soil phages infecting diverse bacteria of different phyla.

CONCLUSIONS

Our results showed that mangrove viruses are diverse and probably directly manipulate carbon cycling by participating in biomass recycling of complex polysaccharides, providing the knowledge essential in revealing the ecological roles of viruses in mangrove ecosystems.

Numerous cultivated and uncultivated viruses encode ribosomal proteins

Viruses modulate ecosystems by directly altering host metabolisms through auxiliary metabolic genes. However, viral genomes are not known to encode the core components of translation machinery, such as ribosomal proteins (RPs). Here, using reference genomes and global-scale viral metagenomic datasets, we identify 14 different RPs across viral genomes arising from cultivated viral isolates and metagenome-assembled viruses. Viruses tend to encode dynamic RPs, easily exchangeable between ribosomes, suggesting these proteins can replace cellular versions in host ribosomes. Functional assays confirm that the two most common virus-encoded RPs, bS21 and bL12, are incorporated into 70S ribosomes when expressed in Escherichia coli. Ecological distribution of virus-encoded RPs suggests some level of ecosystem adaptations as aquatic viruses and viruses of animal-associated bacteria are enriched for different subsets of RPs. Finally, RP genes are under purifying selection and thus likely retained an important function after being horizontally transferred into virus genomes.

Bacterial 'Grounded' Prophages: Hotspots for Genetic Renovation and Innovation

Bacterial genomes are highly plastic allowing the generation of variants through mutations and acquisition of genetic information. The fittest variants are then selected by the econiche thereby allowing the bacterial adaptation and colonization of the habitat. Larger genomes, however, may impose metabolic burden and hence bacterial genomes are optimized by the loss of frivolous genetic information. The activity of temperate bacteriophages has acute consequences on the bacterial population as well as the bacterial genome through lytic and lysogenic cycles. Lysogeny is a selective advantage as the prophage provides immunity to the lysogen against secondary phage attack. Since the non-lysogens are eliminated by the lytic phages, lysogens multiply and colonize the habitat. Nevertheless, all lysogens have an imminent risk of lytic cycle activation and cell lysis. However, a mutation in the attachment sites or in the genes that encode the specific recombinase responsible for prophage excision could result in 'grounding' of the prophage. Since the lysogens with grounded prophage are immune to respective phage infection as well as dodge the induction of lytic cycle, we hypothesize that the selection of these mutant lysogens is favored relative to their normal lysogenic counterparts. These grounded prophages offer several advantages to the bacterial genome evolution through propensity for genetic variations including inversions, deletions, and insertions via horizontal gene transfer. We propose that the grounded prophages expedite bacterial genome evolution by acting as 'genetic buffer zones' thereby increasing the frequency as well as the diversity of variations on which natural selection favors the beneficial variants. The grounded prophages are also hotspots for horizontal gene transfer wherein several ecologically significant genes such as those involved in stress tolerance, antimicrobial resistance, and novel metabolic pathways, are integrated. Moreover, the high frequency of genetic changes within prophages also allows proportionate probability for the de novo genesis of genetic information. Through sequence analyses of well-characterized E. coli prophages we exemplify various roles of grounded prophages in E. coli ecology and evolution. Therefore, the temperate prophages are one of the most significant drivers of bacterial genome evolution and sites of biogenesis of genetic information.

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

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

Host-hijacking and planktonic piracy: how phages command the microbial high seas

Microbial communities living in the oceans are major drivers of global biogeochemical cycles. With nutrients limited across vast swathes of the ocean, marine microbes eke out a living under constant assault from predatory viruses. Viral concentrations exceed those of their bacterial prey by an order of magnitude in surface water, making these obligate parasites the most abundant biological entities in the ocean. Like the pirates of the 17th and 18th centuries that hounded ships plying major trade and exploration routes, viruses have evolved mechanisms to hijack microbial cells and repurpose their cargo and indeed the vessels themselves to maximise viral propagation. Phenotypic reconfiguration of the host is often achieved through Auxiliary Metabolic Genes - genes originally derived from host genomes but maintained and adapted in viral genomes to redirect energy and substrates towards viral synthesis. In this review, we critically evaluate the literature describing the mechanisms used by bacteriophages to reconfigure host metabolism and to plunder intracellular resources to optimise viral production. We also highlight the mechanisms used when, in challenging environments, a 'batten down the hatches' strategy supersedes that of 'plunder and pillage'. Here, the infecting virus increases host fitness through phenotypic augmentation in order to ride out the metaphorical storm, with a concomitant impact on host substrate uptake and metabolism, and ultimately, their interactions with their wider microbial community. Thus, the traditional view of the virus-host relationship as predator and prey does not fully characterise the variety or significance of the interactions observed. Recent advances in viral metagenomics have provided a tantalising glimpse of novel mechanisms of viral metabolic reprogramming in global oceans. Incorporation of these new findings into global biogeochemical models requires experimental evidence from model systems and major improvements in our ability to accurately predict protein function from sequence data.

Libra: scalable k-mer-based tool for massive all-vs-all metagenome comparisons

Background

Shotgun metagenomics provides powerful insights into microbial community biodiversity and function. Yet, inferences from metagenomic studies are often limited by dataset size and complexity and are restricted by the availability and completeness of existing databases. De novo comparative metagenomics enables the comparison of metagenomes based on their total genetic content.

Results

We developed a tool called Libra that performs an all-vs-all comparison of metagenomes for precise clustering based on their k-mer content. Libra uses a scalable Hadoop framework for massive metagenome comparisons, Cosine Similarity for calculating the distance using sequence composition and abundance while normalizing for sequencing depth, and a web-based implementation in iMicrobe (http://imicrobe.us) that uses the CyVerse advanced cyberinfrastructure to promote broad use of the tool by the scientific community.

Conclusions

A comparison of Libra to equivalent tools using both simulated and real metagenomic datasets, ranging from 80 million to 4.2 billion reads, reveals that methods commonly implemented to reduce compute time for large datasets, such as data reduction, read count normalization, and presence/absence distance metrics, greatly diminish the resolution of large-scale comparative analyses. In contrast, Libra uses all of the reads to calculate k-mer abundance in a Hadoop architecture that can scale to any size dataset to enable global-scale analyses and link microbial signatures to biological processes.

Interspecies cross-feeding orchestrates carbon degradation in the rumen ecosystem

Because of their agricultural value, there is a great body of research dedicated to understanding the microorganisms responsible for rumen carbon degradation. However, we lack a holistic view of the microbial food web responsible for carbon processing in this ecosystem. Here, we sampled rumen-fistulated moose, allowing access to rumen microbial communities actively degrading woody plant biomass in real time. We resolved 1,193 viral contigs and 77 unique, near-complete microbial metagenome-assembled genomes, many of which lacked previous metabolic insights. Plant-derived metabolites were measured with NMR and carbohydrate microarrays to quantify the carbon nutrient landscape. Network analyses directly linked measured metabolites to expressed proteins from these unique metagenome-assembled genomes, revealing a genome-resolved three-tiered carbohydrate-fuelled trophic system. This provided a glimpse into microbial specialization into functional guilds defined by specific metabolites. To validate our proteomic inferences, the catalytic activity of a polysaccharide utilization locus from a highly connected metabolic hub genome was confirmed using heterologous gene expression. Viral detected proteins and linkages to microbial hosts demonstrated that phage are active controllers of rumen ecosystem function. Our findings elucidate the microbial and viral members, as well as their metabolic interdependencies, that support in situ carbon degradation in the rumen ecosystem.

A combination of proteomics, metagenome-assembled genomes and heterologous gene expression experiments reveals a trophic system for carbon utilization in the moose rumen microbiome and provides insights into phage dynamics in this ecosystem.

Soil Viruses Are Underexplored Players in Ecosystem Carbon Processing

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

Genetic Diversity and Cooccurrence Patterns of Marine Cyanopodoviruses and Picocyanobacteria

Picocyanobacteria Prochlorococcus and Synechococcus are abundant in the global oceans and subject to active viral infection. In this study, the genetic diversity of picocyanobacteria and the genetic diversity of cyanopodoviruses were synchronously investigated along water columns in the equatorial Indian Ocean and over a seasonal time course in the coastal Sanya Bay, South China Sea. Using the 16S-23S rRNA internal transcribed spacer (ITS)-based clone library and quantitative PCR (qPCR) analyses, the picocyanobacterial community composition and abundance were determined. Sanya Bay was dominated by clade II Synechococcus during all the seasons, and a typical population shift from high-light-adapted Prochlorococcus to low-light-adapted Prochlorococcus was found along the vertical profiles. Strikingly, the DNA polymerase gene sequences of cyanopodoviruses revealed a much greater genetic diversity than we expected. Nearly one-third of the phylogenetic groups were newly described here. No apparent seasonal pattern was observed for the Sanya Bay picocyanobacterial or cyanopodoviral communities. Different dominant cyanopodovirus lineages were identified for the coastal area, upper euphotic zone, and middle-to-lower euphotic zone of the open ocean. Diversity indices of both picocyanobacteria and cyanopodoviruses were highest in the middle euphotic zone and both were lower in the upper euphotic zone, reflecting a host-virus interaction. Cyanopodoviral communities differed significantly between the upper euphotic zone and the middle-to-lower euphotic zone, showing a vertical pattern similar to that of picocyanobacteria. However, in the surface waters of the open ocean, cyanopodoviruses exhibited no apparent biogeographic pattern, differing from picocyanobacteria. This study demonstrates correlated distribution patterns of picocyanobacteria and cyanopodoviruses, as well as the complex biogeography of cyanopodoviruses.IMPORTANCE Picocyanobacteria are highly diverse and abundant in the ocean and display remarkable global biogeography and a vertical distribution pattern. However, how the diversity and distribution of picocyanobacteria affect those of the viruses that infect them remains largely unknown. Here we synchronously analyzed the community structures of cyanopodoviruses and picocyanobacteria at spatial and temporal scales. Both spatial and temporal variations of cyanopodoviral communities can be linked to those of picocyanobacteria. The coastal area, upper euphotic zone, and middle-to-lower euphotic zone of the open ocean have distinct cyanopodoviral communities, showing horizontal and vertical variation patterns closely related to those of picocyanobacteria. These findings emphasize the driving force of host community in shaping the biogeographic structure of viruses. Our work provides important information for future assessments of the ecological roles of viruses and hosts for each other.

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.

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.

Water Masses and Depth Structure Prokaryotic and T4-Like Viral Communities Around Hydrothermal Systems of the Nordic Seas

The oceanographic features of the Nordic Seas, situated between Iceland and Svalbard, have been extensively studied over the last decades. As well, the Nordic Seas hydrothermal systems situated on the Arctic Mid-Ocean Ridge System have received an increasing interest. However, there is very little knowledge on the microbial communities inhabiting the water column of the Nordic Seas, and nothing is known about the influence of the different water masses and hydrothermal plumes on the microbial community structures. In this study, we aimed at characterizing the impact of hydrothermal plumes on prokaryotic and T4-like viral communities around the island of Jan Mayen. To this end, we used 16S rRNA-gene and g23-gene profiling as well as flow cytometry counts to examine prokaryotic and viral communities in 27 samples obtained from different water masses in this area. While Thaumarchaeota and Marine group II Archaea dominated the waters deeper than 500 m, members of Flavobacteria generally dominated the shallower waters. Furthermore, extensive chemical and physical characteristics of all samples were obtained, including temperature measurements and concentrations of major ions and gases. The effect of these physiochemical variables on the communities was measured by using constrained and unconstrained multivariate analyzes, Mantel tests, network analyzes, phylogenetic analyzes, taxonomic analyzes and temperature-salinity (Θ-S) plots. Our results suggest that hydrothermal activity has little effect on pelagic microbial communities in hydrothermal plumes of the Nordic Seas. However, we provide evidences that observed differences in prokaryotic community structure can largely be attributed to which water mass each sample was taken from. In contrast, depth was the major factor structuring the T4-like viral communities. Our results also show that it is crucial to include water masses when studying the influence of hydrothermal plumes on microbial communities, as it could prevent to falsely associate a change in community structure with the presence of a plume.

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.

The 'Neglected' Soil Virome - Potential Role and Impact

Bacteriophages are among the most abundant and diverse biological units in the biosphere. They have contributed to our understanding of the central dogma of biology and have been instrumental in the evolutionary success of bacterial pathogens. In contrast to our current understanding of marine viral communities, the soil virome and its function in terrestrial ecosystems has remained relatively understudied. Here, we examine, in a comparative fashion, the knowledge gathered from studies performed in soil versus marine settings. We address the information with respect to the abundance, diversity, ecological significance, and effects of, in particular, bacteriophages on their host's evolutionary trajectories. We also identify the main challenges that soil virology faces and the studies that are required to accompany the current developments in marine settings.

Iron Chloride Flocculation of Bacteriophages from Seawater

Viruses influence ecosystem dynamics by modulating microbial host population dynamics, evolutionary trajectories and metabolic outputs. While they are ecologically important across diverse ecosystems, viruses are challenging to study due to minimal biomass often obtained when sampling natural communities. Here we describe a technique using chemical flocculation, filtration and resuspension to recover bacteriophages from seawater and other natural waters. The method uses iron to precipitate viruses which are recovered by filtration onto large-pore size membranes and then resuspended using a buffer containing magnesium and a reductant (ascorbic acid or oxalic acid) at slightly acid pH (6-6.5). The recovery of bacteriophages using iron flocculation is efficient (>90%), inexpensive and reliable, resulting in preparations that are amenable to downstream analysis by next generation DNA sequencing, proteomics and, in some cases, can be used to study virus-host interactions.

Dietary energy drives the dynamic response of bovine rumen viral communities

BACKGROUND

Rumen microbes play a greater role in host energy acquisition than that of gut-associated microbes in monogastric animals. Although genome-enabled advancements are providing access to the vast diversity of uncultivated microbes, our understanding of variables shaping rumen microbial communities is in its infancy. Viruses have been shown to impact microbial populations through a myriad of processes, including cell lysis and reprogramming of host metabolism. However, little is known about the processes shaping the distribution of rumen viruses or how viruses may modulate microbial-driven processes in the rumen. To this end, we investigated how rumen bacterial and viral community structure and function responded in five steers fed four randomized dietary treatments in a crossover design.

RESULTS

Total digestible nutrients (TDN), a measure of dietary energy, best explained the variation in bacterial and viral communities. Additional ecological drivers of viral communities included dietary zinc content and microbial functional diversity. Using partial least squares regression, we demonstrate significant associations between the abundances of 267 viral populations and variables driving the variation in rumen viral communities. While rumen viruses were dynamic, 14 near ubiquitous viral populations were identified, suggesting the presence of a core rumen virome largely comprised of novel viruses. Moreover, analysis of virally encoded auxiliary metabolic genes (AMGs) indicates rumen viruses have glycosidic hydrolases to potentially augment the breakdown of complex carbohydrates to increase energy production. Other AMGs identified have a role in redirecting carbon to the pentose phosphate pathway and one carbon pools by folate to boost viral replication.

CONCLUSIONS

We demonstrate that rumen bacteria and viruses have differing responses and ecological drivers to dietary perturbation. Our results show that rumen viruses have implications for understanding the structuring of the previously identified core rumen microbiota and impacting microbial metabolism through a vast array of AMGs. AMGs in the rumen appear to have consequences for microbial metabolism that are largely in congruence with the current paradigm established in marine systems. This study provides a foundation for future hypotheses regarding the dynamics of viral-mediated processes in the rumen.

Phage hunters: Computational strategies for finding phages in large-scale 'omics datasets

A plethora of tools exist for identifying phage sequences in bacterial genomes, single cell amplified genomes, and host-associated and environmental metagenomes. Yet because the genetics of phages and their hosts are closely intertwined, distinguishing viral from bacterial signal remains an ongoing challenge. Further the size, quantity and fragmentary nature of modern 'omics datasets ushers in a new set of computational challenges. Here, we detail the promises and pitfalls of using currently available gene-centric or k-mer based tools for identifying prophage sequences in genomes and prophage and viral contigs in metagenomes. Each of these methods offers a unique piece of the puzzle to elucidating the intriguing signatures of phage-host coevolution.

Virioplankton Assemblage Structure in the Lower River and Ocean Continuum of the Amazon

The Amazon River forms a vast plume in the Atlantic Ocean that can extend for more than 1,000 km. Microbial communities promote a globally relevant carbon sink system in the plume. Despite the importance of viruses for the global carbon cycle, the diversity and the possible roles of viruses in the Amazonia are poorly understood. The present work assesses, for the first time, the abundance and diversity of viruses simultaneously in the river and ocean in order to elucidate their possible roles. DNA sequence assembly yielded 29,358 scaffolds, encoding 82,546 viral proteins, with 15 new complete viral genomes from the 12 river and ocean locations. Viral diversity was clearly distinguished by river and ocean. Bacteriophages were the most abundant and occurred throughout the continuum. Viruses that infect eukaryotes were more abundant in the river, whereas phages appeared to have strong control over the host prokaryotic populations in the plume.

The Amazon River watershed and its associated plume comprise a vast continental and oceanic area. The microbial activities along this continuum contribute substantially to global carbon and nutrient cycling, and yet there is a dearth of information on the diversity, abundance, and possible roles of viruses in this globally important river. The aim of this study was to elucidate the diversity and structure of virus assemblages of the Amazon River-ocean continuum. Environmental viral DNA sequences were obtained for 12 locations along the river’s lower reach (n = 5) and plume (n = 7). Sequence assembly yielded 29,358 scaffolds, encoding 82,546 viral proteins, with 15 new complete viral genomes. Despite the spatial connectivity mediated by the river, virome analyses and physical-chemical water parameters clearly distinguished river and plume ecosystems. Bacteriophages were ubiquitous in the continuum and were more abundant in the transition region. Eukaryotic viruses occurred mostly in the river, while the plume had more viruses of autotrophic organisms (Prochlorococcus, Synechococcus) and heterotrophic bacteria (Pelagibacter). The viral families Microviridae and Myoviridae were the most abundant and occurred throughout the continuum. The major functions of the genes in the continuum involved viral structures and life cycles, and viruses from plume locations and Tapajós River showed the highest levels of functional diversity. The distribution patterns of the viral assemblages were defined not only by the occurrence of possible hosts but also by water physical and chemical parameters, especially salinity. The findings presented here help to improve understanding of the possible roles of viruses in the organic matter cycle along the river-ocean continuum.

IMPORTANCE The Amazon River forms a vast plume in the Atlantic Ocean that can extend for more than 1,000 km. Microbial communities promote a globally relevant carbon sink system in the plume. Despite the importance of viruses for the global carbon cycle, the diversity and the possible roles of viruses in the Amazon are poorly understood. The present work assesses, for the first time, the abundance and diversity of viruses simultaneously in the river and ocean in order to elucidate their possible roles. DNA sequence assembly yielded 29,358 scaffolds, encoding 82,546 viral proteins, with 15 new complete viral genomes from the 12 river and ocean locations. Viral diversity was clearly distinguished by river and ocean. Bacteriophages were the most abundant and occurred throughout the continuum. Viruses that infect eukaryotes were more abundant in the river, whereas phages appeared to have strong control over the host prokaryotic populations in the plume.

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.

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

Background

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

Results

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

Conclusions

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

Genome analysis of three novel lytic Vibrio coralliilyticus phages isolated from seawater, Okinawa, Japan

Three novel Vibrio phages were isolated from seawater in Okinawa. The Vibrio phage RYC infected Vibrio coralliilyticus SWA 07, while Vibrio phages CKB-S1 and CKB-S2 infected the coral pathogen V. coralliilyticus P1 (LMG 23696). The Vibrio phages CKB-S1 and CKB-S2 displayed head-tail structures whereas the Vibrio phage RYC showed a tailless non-enveloped capsid. All these Vibrio phages contained linear and double-stranded DNA. The whole genome sequencing revealed that Vibrio phage RYC has a larger genome size compared to Vibrio phages CKB-S1 and CKB-S2, and six tRNAs genes were found only in Vibrio phage RYC. Genome-wide comparison showed that Vibrio phage CKB-S1 was closely related, but was not identical, to Vibrio parahaemolyticus phages VP16T and VP16C. Meanwhile, the Vibrio phages RYC and CKB-S2 did not show high genome-wide similarity to any phages. These results suggest that the Vibrio phages CKB-S1, CKB-S2 and RYC are novel phages, which need further exploration, especially for their potential applications in phage therapy.

Comparative Omics and Trait Analyses of Marine Pseudoalteromonas Phages Advance the Phage OTU Concept

Viruses influence the ecology and evolutionary trajectory of microbial communities. Yet our understanding of their roles in ecosystems is limited by the paucity of model systems available for hypothesis generation and testing. Further, virology is limited by the lack of a broadly accepted conceptual framework to classify viral diversity into evolutionary and ecologically cohesive units. Here, we introduce genomes, structural proteomes, and quantitative host range data for eight Pseudoalteromonas phages isolated from Helgoland (North Sea, Germany) and use these data to advance a genome-based viral operational taxonomic unit (OTU) definition. These viruses represent five new genera and inform 498 unaffiliated or unannotated protein clusters (PCs) from global virus metagenomes. In a comparison of previously sequenced Pseudoalteromonas phage isolates (n = 7) and predicted prophages (n = 31), the eight phages are unique. They share a genus with only one other isolate, Pseudoalteromonas podophage RIO-1 (East Sea, South Korea) and two Pseudoalteromonas prophages. Mass-spectrometry of purified viral particles identified 12–20 structural proteins per phage. When combined with 3-D structural predictions, these data led to the functional characterization of five previously unidentified major capsid proteins. Protein functional predictions revealed mechanisms for hijacking host metabolism and resources. Further, they uncovered a hybrid sipho-myovirus that encodes genes for Mu-like infection rarely described in ocean systems. Finally, we used these data to evaluate a recently introduced definition for virus populations that requires members of the same population to have >95% average nucleotide identity across at least 80% of their genes. Using physiological traits and genomics, we proposed a conceptual model for a viral OTU definition that captures evolutionarily cohesive and ecologically distinct units. In this trait-based framework, sensitive hosts are considered viral niches, while host ranges and infection efficiencies are tracked as viral traits. Quantitative host range assays revealed conserved traits within virus OTUs that break down between OTUs, suggesting the defined units capture niche and fitness differentiation. Together these analyses provide a foundation for model system-based hypothesis testing that will improve our understanding of marine copiotrophs, as well as phage–host interactions on the ocean particles and aggregates where Pseudoalteromonas thrive.

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.

Contribution of Next-Generation Sequencing to Aquatic and Fish Virology

The recent technological advances in nucleic acid sequencing, called next-generation sequencing (NGS), have revolutionized the field of genomics and have also influenced viral research. Aquatic viruses, and especially those infecting fish, have also greatly benefited from NGS technologies, which provide a huge amount of molecular information at a low cost in a relatively short period of time. Here, we review the use of the current high-throughput sequencing platforms with a special focus on the associated challenges (regarding sample preparation and bioinformatics) in their applications to the field of aquatic virology, especially for: (i) discovering novel viruses that may be associated with fish mortalities, (ii) elucidating the mechanisms of pathogenesis, and finally (iii) studying the molecular epidemiology of these pathogens.

Metagenomic Approaches to Assess Bacteriophages in Various Environmental Niches

Bacteriophages are ubiquitous and numerous parasites of bacteria and play a critical evolutionary role in virtually every ecosystem, yet our understanding of the extent of the diversity and role of phages remains inadequate for many ecological niches, particularly in cases in which the host is unculturable. During the past 15 years, the emergence of the field of viral metagenomics has drastically enhanced our ability to analyse the so-called viral ‘dark matter’ of the biosphere. Here, we review the evolution of viral metagenomic methodologies, as well as providing an overview of some of the most significant applications and findings in this field of research.

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.

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.

Visualizing Adsorption of Cyanophage P-SSP7 onto Marine Prochlorococcus

Marine cyanobacteria perform roughly a quarter of global carbon fixation, and cyanophages that infect them liberate some of this carbon during infection and cell lysis. Studies of the cyanobacterium Prochlorococcus MED4 and its associated cyanophage P-SSP7 have revealed complex gene expression dynamics once infection has begun, but the initial cyanophage-host interactions remain poorly understood. Here, we used single particle cryo-electron tomography (cryo-ET) to investigate cyanophage-host interactions in this model system, based on 170 cyanophage-to-host adsorption events. Subtomogram classification and averaging revealed three main conformations characterized by different angles between the phage tail and the cell surface. Namely, phage tails were (i) parallel to, (ii) ~45 degrees to, or (iii) perpendicular to the cell surface. Furthermore, different conformations of phage tail fibers correlated with the aforementioned orientations of the tails. We also observed density beyond the tail tip in vertically-oriented phages that had penetrated the cell wall, capturing the final stage of adsorption. Together, our data provide a quantitative characterization of the orientation of phages as they adsorb onto cells, and suggest that cyanophages that abut their cellular targets are only transiently in the “perpendicular” orientation required for successful infection.

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.

Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses

Ocean microbes drive biogeochemical cycling on a global scale. However, this cycling is constrained by viruses that affect community composition, metabolic activity, and evolutionary trajectories. Owing to challenges with the sampling and cultivation of viruses, genome-level viral diversity remains poorly described and grossly understudied, with less than 1% of observed surface-ocean viruses known. Here we assemble complete genomes and large genomic fragments from both surface- and deep-ocean viruses sampled during the Tara Oceans and Malaspina research expeditions, and analyse the resulting 'global ocean virome' dataset to present a global map of abundant, double-stranded DNA viruses complete with genomic and ecological contexts. A total of 15,222 epipelagic and mesopelagic viral populations were identified, comprising 867 viral clusters (defined as approximately genus-level groups). This roughly triples the number of known ocean viral populations and doubles the number of candidate bacterial and archaeal virus genera, providing a near-complete sampling of epipelagic communities at both the population and viral-cluster level. We found that 38 of the 867 viral clusters were locally or globally abundant, together accounting for nearly half of the viral populations in any global ocean virome sample. While two-thirds of these clusters represent newly described viruses lacking any cultivated representative, most could be computationally linked to dominant, ecologically relevant microbial hosts. Moreover, we identified 243 viral-encoded auxiliary metabolic genes, of which only 95 were previously known. Deeper analyses of four of these auxiliary metabolic genes (dsrC, soxYZ, P-II (also known as glnB) and amoC) revealed that abundant viruses may directly manipulate sulfur and nitrogen cycling throughout the epipelagic ocean. This viral catalog and functional analyses provide a necessary foundation for the meaningful integration of viruses into ecosystem models where they act as key players in nutrient cycling and trophic networks.

Genomes of Abundant and Widespread Viruses from the Deep Ocean

The deep sea is a massive, largely oligotrophic ecosystem, stretched over nearly 65% of the planet’s surface. Deep-sea planktonic communities are almost completely dependent upon organic carbon sinking from the productive surface, forming a vital component of global biogeochemical cycles. However, despite their importance, viruses from the deep ocean remain largely unknown. Here, we describe the first complete genomes of deep-sea viruses assembled from metagenomic fosmid libraries. “Candidatus Pelagibacter” (SAR11) phage HTVC010P and Puniceispirillum phage HMO-2011 are considered the most abundant cultured marine viruses known to date. Remarkably, some of the viruses described here recruited as many reads from deep waters as these viruses do in the photic zone, and, considering the gigantic scale of the bathypelagic habitat, these genomes provide information about what could be some of the most abundant viruses in the world at large. Their role in the viral shunt in the global ocean could be very significant. Despite the challenges encountered in inferring the identity of their hosts, we identified one virus predicted to infect members of the globally distributed SAR11 cluster. We also identified a number of putative proviruses from diverse taxa, including deltaproteobacteria, bacteroidetes, SAR11, and gammaproteobacteria. Moreover, our findings also indicate that lysogeny is the preferred mode of existence for deep-sea viruses inhabiting an energy-limited environment, in sharp contrast to the predominantly lytic lifestyle of their photic-zone counterparts. Some of the viruses show a widespread distribution, supporting the tenet “everything is everywhere” for the deep-ocean virome.

The deep sea is among the largest known habitats and a critical cog in biogeochemical cycling but remains underexplored in its microbiology. Even more than is the case for its prokaryotic community, our knowledge of its viral component has remained limited by the paucity of information provided by studies dependent upon short sequence fragments. In this work, we attempt to fill this existing gap by using a combination of classical fosmid libraries with next-generation sequencing and assembly to recover long viral genomic fragments. We have sequenced ca. 6,000 fosmids from two metagenomics libraries made from prokaryotic biomass from the deep Mediterranean Sea and recovered twenty-eight complete viral genomes, all of them novel and quite distinct from all previously described viral genomes. They are preferentially found in deeper waters and are widely distributed all over the oceans. To our knowledge, this is the first report on complete and cosmopolitan viral genomes from the bathypelagic habitat.

iVirus: facilitating new insights in viral ecology with software and community data sets imbedded in a cyberinfrastructure

Microbes affect nutrient and energy transformations throughout the world's ecosystems, yet they do so under viral constraints. In complex communities, viral metagenome (virome) sequencing is transforming our ability to quantify viral diversity and impacts. Although some bottlenecks, for example, few reference genomes and nonquantitative viromics, have been overcome, the void of centralized data sets and specialized tools now prevents viromics from being broadly applied to answer fundamental ecological questions. Here we present iVirus, a community resource that leverages the CyVerse cyberinfrastructure to provide access to viromic tools and data sets. The iVirus Data Commons contains both raw and processed data from 1866 samples and 73 projects derived from global ocean expeditions, as well as existing and legacy public repositories. Through the CyVerse Discovery Environment, users can interrogate these data sets using existing analytical tools (software applications known as ‘Apps') for assembly, open reading frame prediction and annotation, as well as several new Apps specifically developed for analyzing viromes. Because Apps are web based and powered by CyVerse supercomputing resources, they enable scalable analyses for a broad user base. Finally, a use-case scenario documents how to apply these advances toward new data. This growing iVirus resource should help researchers utilize viromics as yet another tool to elucidate viral roles in nature.