Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses

Unravelling the ecological complexity of soil viromes: Challenges and opportunities

Viruses are extremely abundant and ubiquitous in soil, and significantly contribute to various terrestrial ecosystem processes such as biogeochemical nutrient cycling, microbiome regulation and community assembly, and host evolutionary dynamics. Despite their numerous dominance and functional importance, understanding soil viral ecology is a formidable challenge, because of the technological challenges to characterize the abundance, diversity and community compositions of viruses, and their interactions with other organisms in the complex soil environment. Viruses may engage in a myriad of biological interactions within soil food webs across a broad range of spatiotemporal scales and are exposed to various biotic and abiotic disturbances. Current studies on the soil viromes, however, often describe the complexity of their tremendous diversity, but lack of exploring their potential ecological roles. In this article, we summarized the major methods to decipher the ecology of soil viruses, discussed biotic and abiotic factors and global change factors that shape the diversity and composition of soil viromes, and the ecological roles of soil viruses. We also proposed a new framework to understand the ecological complexity of viruses from micro to macro ecosystem scales and to predict and unravel their activities in terrestrial ecosystems.

Microbial Diversity in the Indian Ocean Sediments: An Insight into the Distribution and Associated Factors

Indian Ocean is the third largest oceanic division of the world and shelter to a huge microbial diversity. These microbes play an important role in the metabolism of carbon, sulfur, nitrogen, and phosphorus in the ocean water. They are also major contributors of carbon fixing and sequestration, as much as terrestrial plants to achieve CO₂ emissions reduction. The prokaryotic community in the East Indian Ocean primarily comprises of heterotrophic bacteria like Alphaproteobacteria and Gammaproteobacteria, followed by Firmicutes and Actinobacteria. The Arabian Sea and the Bay of Bengal are typically characterized by presence of vast areas of oxygen minimum zones (OMZs) and have been witnessing a shift in the microbial diversity due to the changing conditions in the ocean water. Several canonical correspondence analyses reveal temperature, salinity, and phosphate levels as crucial environmental factors in propelling the distribution of diazotrophs. The viral consortia are dominated by the Caudovirales, an order of tailed bacteriophages. Due to the rapid change in the environmental factors such as topography, temperature, and sunlight contributing toward climate change, their role in sustaining the chemical composition of the ocean can be drastically affected especially with the evidence of several bacterial and fungal communities responding to latitudinal and temperature change. Therefore, we aim to critically review the status of microbial diversity in Indian Ocean to predict their response toward climate change as they are the sentinels of change in marine life and to understand the dynamics of microbial communities in the various locations of Indian Ocean.

Diversity of Pseudomonas aeruginosa Temperate Phages

Modern sequencing technologies have provided insight into the genetic diversity of numerous species, including the human pathogen Pseudomonas aeruginosa. Bacterial genomes often harbor bacteriophage genomes (prophages), which can account for upwards of 20% of the genome. Prior studies have found P. aeruginosa prophages that contribute to their host’s pathogenicity and fitness. These advantages come in many different forms, including the production of toxins, promotion of biofilm formation, and displacement of other P. aeruginosa strains. While several different genera and species of P. aeruginosa prophages have been studied, there has not been a comprehensive study of the overall diversity of P. aeruginosa-infecting prophages. Here, we present the results of just such an analysis. A total of 6,852 high-confidence prophages were identified from 5,383 P. aeruginosa genomes from strains isolated from the human body and other environments. In total, 3,201 unique prophage sequences were identified. While 53.1% of these prophage sequences displayed sequence similarity to publicly available phage genomes, novel and highly mosaic prophages were discovered. Among these prophages, there is extensive diversity, including diversity within the functionally conserved integrase and C repressor coding regions, two genes responsible for prophage entering and persisting through the lysogenic life cycle. Analysis of integrase, C repressor, and terminase coding regions revealed extensive reassortment among P. aeruginosa prophages. This catalog of P. aeruginosa prophages provides a resource for future studies into the evolution of the species.

IMPORTANCE Prophages play a critical role in the evolution of their host species and can also contribute to the virulence and fitness of pathogenic species. Here, we conducted a comprehensive investigation of prophage sequences from 5,383 publicly available Pseudomonas aeruginosa genomes from human as well as environmental isolates. We identified a diverse population of prophages, including tailed phages, inoviruses, and microviruses; 46.9% of the prophage sequences found share no significant sequence similarity with characterized phages, representing a vast array of novel P. aeruginosa-infecting phages. Our investigation into these prophages found substantial evidence of reassortment. In producing this, the first catalog of P. aeruginosa prophages, we uncovered both novel prophages as well as genetic content that have yet to be explored.

Novel Viral Communities Potentially Assisting in Carbon, Nitrogen, and Sulfur Metabolism in the Upper Slope Sediments of Mariana Trench

Viruses are ubiquitous in the oceans. Even in the deep sediments of the Mariana Trench, viruses have high productivity. However, little is known about their species composition and survival strategies in that environment. Here, we uncovered novel viral communities (3,206 viral scaffolds) in the upper slope sediments of the Mariana Trench via metagenomic analysis of 15 sediment samples. Most (99%) of the viral scaffolds lack known viral homologs, and ca. 59% of the high-quality viral genomes (total of 111 with completeness of >90%) represent novel genera, including some Phycodnaviridae and jumbo phages. These viruses contain various auxiliary metabolic genes (AMGs) potentially involved in organic carbon degradation, inorganic carbon fixation, denitrification, and assimilatory sulfate reduction, etc. This study provides novel insight into the almost unknown benthic viral communities in the Mariana Trench. IMPORTANCE The Mariana Trench harbors a substantial number of infective viral particles. However, very little is known about the identity, survival strategy, and potential functions of viruses in the trench sediments. Here, through metagenomic analysis, unusual benthic viral communities with high diversity and novelty were discovered. Among them, 59% of the viruses with a genome completeness of >90% represent novel genera. Various auxiliary metabolic genes carried by these viruses reflect the potential adaptive characteristics of viruses in this extreme environment and the biogeochemical cycles that they may participate in. This study gives us a deeper understanding of the peculiarities of viral communities in deep-sea/hadal sediments.

Distribution of soil viruses across China and their potential role in phosphorous metabolism

BACKGROUND

Viruses are the most abundant biological entities on the planet and drive biogeochemical cycling on a global scale. Our understanding of biogeography of soil viruses and their ecological functions lags significantly behind that of Bacteria and Fungi. Here, a viromic approach was used to investigate the distribution and ecological functions of viruses from 19 soils across China.

RESULTS

Soil viral community were clustered more significantly by geographical location than type of soil (agricultural and natural). Three clusters of viral communities were identified from North, Southeast and Southwest regions; these clusters differentiated using taxonomic composition and were mainly driven by geographic location and climate factors. A total of 972 viral populations (vOTUs) were detected spanning 23 viral families from the 19 viromes. Phylogenetic analyses of the phoH gene showed a remarkable diversity and the distribution of viral phoH genes was more dependent on the environment. Notably, five proteins involved in phosphorus (P) metabolism-related nucleotide synthesis functions, including dUTPase, MazG, PhoH, Thymidylate synthase complementing protein (Thy1), and Ribonucleoside reductase (RNR), were mainly identified in agricultural soils.

CONCLUSIONS

The present work revealed that soil viral communities were distributed across China according to geographical location and climate factors. In addition, P metabolism genes encoded by these viruses probably drive the synthesis of nucleotides for their own genomes inside bacterial hosts, thereby affecting P cycling in the soil ecosystems.

Highly diverse flavobacterial phages isolated from North Sea spring blooms

It is generally recognized that phages are a mortality factor for their bacterial hosts. This could be particularly true in spring phytoplankton blooms, which are known to be closely followed by a highly specialized bacterial community. We hypothesized that phages modulate these dense heterotrophic bacteria successions following phytoplankton blooms. In this study, we focused on Flavobacteriia, because they are main responders during these blooms and have an important role in the degradation of polysaccharides. A cultivation-based approach was used, obtaining 44 lytic flavobacterial phages (flavophages), representing twelve new species from two viral realms. Taxonomic analysis allowed us to delineate ten new phage genera and ten new families, from which nine and four, respectively, had no previously cultivated representatives. Genomic analysis predicted various life styles and genomic replication strategies. A likely eukaryote-associated host habitat was reflected in the gene content of some of the flavophages. Detection in cellular metagenomes and by direct-plating showed that part of these phages were actively replicating in the environment during the 2018 spring bloom. Furthermore, CRISPR/Cas spacers and re-isolation during two consecutive years suggested that, at least part of the new flavophages are stable components of the microbial community in the North Sea. Together, our results indicate that these diverse flavophages have the potential to modulate their respective host populations.

Host Cyanobacteria Killing by Novel Lytic Cyanophage YongM: A Protein Profiling Analysis

Cyanobacteria are autotrophic prokaryotes that can proliferate robustly in eutrophic waters through photosynthesis. This can lead to outbreaks of lake “water blooms”, which result in water quality reduction and environmental pollution that seriously affect fisheries and aquaculture. The use of cyanophages to control the growth of cyanobacteria is an important strategy to tackle annual cyanobacterial blooms. YongM is a novel lytic cyanophage with a broad host spectrum and high efficiency in killing its host, cyanobacteria FACHB-596. However, changes in cyanophage protein profile during infestation and killing of the host remains unknown. To characterize the proteins and its regulation networks involved in the killing of host cyanobacteria by YongM and evaluate whether this strain YongM could be used as a chassis for further engineering to be a powerful tool in dealing with cyanobacterial blooms, we herein applied 4D label-free high-throughput quantitative proteomics to analyze differentially expressed proteins (DEPs) involved in cyanobacteria host response infected 1 and 8 h with YongM cyanophage. Metabolic pathways, such as photosynthesis, photosynthesis-antennal protein, oxidative phosphorylation, ribosome, carbon fixation, and glycolysis/glycol-isomerization were significantly altered in the infested host, whereas DEPs were associated with the metabolic processes of photosynthesis, precursor metabolites, energy production, and organic nitrogen compounds. Among these DEPs, key proteins involved in YongM-host interaction may be photosystem I P700 chlorophyll-a apolipoprotein, carbon dioxide concentration mechanism protein, cytochrome B, and some YongM infection lysis-related enzymes. Our results provide comprehensive information of protein profiles during the invasion and killing of host cyanobacteria by its cyanophage, which may shed light on future design and manipulation of artificial cyanophages against water blooms.

Hidden diversity of double-stranded DNA phages in symbiotic Rhizobium species

In this study, we addressed the extent of diversification of phages associated with nitrogen-fixing symbiotic Rhizobium species. Despite the ecological and economic importance of the Rhizobium genus, little is known about the diversity of the associated phages. A thorough assessment of viral diversity requires investigating both lytic phages and prophages harboured in diverse Rhizobium genomes. Protein-sharing networks identified 56 viral clusters (VCs) among a set of 425 isolated phages and predicted prophages. The VCs formed by phages had more proteins in common and a higher degree of synteny, and they group together in clades in the associated phylogenetic tree. By contrast, the VCs of prophages showed significant genetic variation and gene loss, with selective pressure on the remaining genes. Some VCs were found in various Rhizobium species and geographical locations, suggesting that they have wide host ranges. Our results indicate that the VCs represent distinct taxonomic units, probably representing taxa equivalent to genera or even species. The finding of previously undescribed phage taxa indicates the need for further exploration of the diversity of phages associated with Rhizobium species. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Microbial colonization and persistence in deep fractured shales is guided by metabolic exchanges and viral predation

BACKGROUND

Microbial colonization of subsurface shales following hydraulic fracturing offers the opportunity to study coupled biotic and abiotic factors that impact microbial persistence in engineered deep subsurface ecosystems. Shale formations underly much of the continental USA and display geographically distinct gradients in temperature and salinity. Complementing studies performed in eastern USA shales that contain brine-like fluids, here we coupled metagenomic and metabolomic approaches to develop the first genome-level insights into ecosystem colonization and microbial community interactions in a lower-salinity, but high-temperature western USA shale formation.

RESULTS

We collected materials used during the hydraulic fracturing process (i.e., chemicals, drill muds) paired with temporal sampling of water produced from three different hydraulically fractured wells in the STACK (Sooner Trend Anadarko Basin, Canadian and Kingfisher) shale play in OK, USA. Relative to other shale formations, our metagenomic and metabolomic analyses revealed an expanded taxonomic and metabolic diversity of microorganisms that colonize and persist in fractured shales. Importantly, temporal sampling across all three hydraulic fracturing wells traced the degradation of complex polymers from the hydraulic fracturing process to the production and consumption of organic acids that support sulfate- and thiosulfate-reducing bacteria. Furthermore, we identified 5587 viral genomes and linked many of these to the dominant, colonizing microorganisms, demonstrating the key role that viral predation plays in community dynamics within this closed, engineered system. Lastly, top-side audit sampling of different source materials enabled genome-resolved source tracking, revealing the likely sources of many key colonizing and persisting taxa in these ecosystems.

CONCLUSIONS

These findings highlight the importance of resource utilization and resistance to viral predation as key traits that enable specific microbial taxa to persist across fractured shale ecosystems. We also demonstrate the importance of materials used in the hydraulic fracturing process as both a source of persisting shale microorganisms and organic substrates that likely aid in sustaining the microbial community. Moreover, we showed that different physicochemical conditions (i.e., salinity, temperature) can influence the composition and functional potential of persisting microbial communities in shale ecosystems. Together, these results expand our knowledge of microbial life in deep subsurface shales and have important ramifications for management and treatment of microbial biomass in hydraulically fractured wells. Video Abstract.

Saline lakes on the Qinghai-Tibet Plateau harbor unique viral assemblages mediating microbial environmental adaption

The highest plateau on Earth, Qinghai-Tibet Plateau, contains thousands of lakes with broad salinity and diverse and unique microbial communities. However, little is known about their co-occurring viruses. Herein, we identify 4,560 viral Operational Taxonomic Units (vOTUs) from six viromes of three saline lakes on Qinghai-Tibet Plateau, with less than 1% that could be classified. Most of the predicted vOTUs were associated with the dominant bacterial and archaeal phyla. Virus-encoded auxiliary metabolic genes suggest that viruses influence microbial metabolisms of carbon, nitrogen, sulfur, and lipid; the antibiotic resistance mediation; and their salinity adaption. The six viromes clustered together with the ice core viromes and bathypelagic ocean viromes and might represent a new viral habitat. This study has revealed the unique characteristics and potential ecological roles of DNA viromes in the lakes of the highest plateau and established a foundation for the recognition of the viral roles in plateau lake ecosystems.

From Trees to Clouds: PhageClouds for Fast Comparison of ∼640,000 Phage Genomic Sequences and Host-Centric Visualization Using Genomic Network Graphs

Background: Fast and computationally efficient strategies are required to explore genomic relationships within an increasingly large and diverse phage sequence space. Here, we present PhageClouds, a novel approach using a graph database of phage genomic sequences and their intergenomic distances to explore the phage genomic sequence space. Methods: A total of 640,000 phage genomic sequences were retrieved from a variety of databases and public virome assemblies. Intergenomic distances were calculated with dashing, an alignment-free method suitable for handling massive data sets. These data were used to build a Neo4j^® graph database. Results: PhageClouds supported the search of related phages among all complete phage genomes from GenBank for a single query phage in just 10 s. Moreover, PhageClouds expanded the number of closely related phage sequences detected for both finished and draft phage genomes, in comparison with searches exclusively targeting phage entries from GenBank. Conclusions: PhageClouds is a novel resource that will facilitate the analysis of phage genomic sequences and the characterization of assembled phage genomes.

Multi-species host range of staphylococcal phages isolated from wastewater

The host range of bacteriophages defines their impact on bacterial communities and genome diversity. Here, we characterize 94 novel staphylococcal phages from wastewater and establish their host range on a diversified panel of 117 staphylococci from 29 species. Using this high-resolution phage-bacteria interaction matrix, we unveil a multi-species host range as a dominant trait of the isolated staphylococcal phages. Phage genome sequencing shows this pattern to prevail irrespective of taxonomy. Network analysis between phage-infected bacteria reveals that hosts from multiple species, ecosystems, and drug-resistance phenotypes share numerous phages. Lastly, we show that phages throughout this network can package foreign genetic material enclosing an antibiotic resistance marker at various frequencies. Our findings indicate a weak host specialism of the tested phages, and therefore their potential to promote horizontal gene transfer in this environment.

Minnesota peat viromes reveal terrestrial and aquatic niche partitioning for local and global viral populations

BACKGROUND

Peatlands are expected to experience sustained yet fluctuating higher temperatures due to climate change, leading to increased microbial activity and greenhouse gas emissions. Despite mounting evidence for viral contributions to these processes in peatlands underlain with permafrost, little is known about viruses in other peatlands. More generally, soil viral biogeography and its potential drivers are poorly understood at both local and global scales. Here, 87 metagenomes and five viral size-fraction metagenomes (viromes) from a boreal peatland in northern Minnesota (the SPRUCE whole-ecosystem warming experiment and surrounding bog) were analyzed for dsDNA viral community ecological patterns, and the recovered viral populations (vOTUs) were compared with our curated PIGEON database of 266,125 vOTUs from diverse ecosystems.

RESULTS

Within the SPRUCE experiment, viral community composition was significantly correlated with peat depth, water content, and carbon chemistry, including CH₄ and CO₂ concentrations, but not with temperature during the first 2 years of warming treatments. Peat vOTUs with aquatic-like signatures (shared predicted protein content with marine and/or freshwater vOTUs) were significantly enriched in more waterlogged surface peat depths. Predicted host ranges for SPRUCE vOTUs were relatively narrow, generally within a single bacterial genus. Of the 4326 SPRUCE vOTUs, 164 were previously detected in other soils, mostly peatlands. None of the previously identified 202,371 marine and freshwater vOTUs in our PIGEON database were detected in SPRUCE peat, but 0.4% of 80,714 viral clusters (VCs, grouped by predicted protein content) were shared between soil and aquatic environments. On a per-sample basis, vOTU recovery was 32 times higher from viromes compared with total metagenomes.

CONCLUSIONS

Results suggest strong viral "species" boundaries between terrestrial and aquatic ecosystems and to some extent between peat and other soils, with differences less pronounced at higher taxonomic levels. The significant enrichment of aquatic-like vOTUs in more waterlogged peat suggests that viruses may also exhibit niche partitioning on more local scales. These patterns are presumably driven in part by host ecology, consistent with the predicted narrow host ranges. Although more samples and increased sequencing depth improved vOTU recovery from total metagenomes, the substantially higher per-sample vOTU recovery after viral particle enrichment highlights the utility of soil viromics. Video abstract The importance of Minnesota peat viromes in revealing terrestrial and aquatic niche partitioning for viral populations.

Perspective on taxonomic classification of uncultivated viruses

Historically, virus taxonomy has been limited to describing viruses that were readily cultivated in the laboratory or emerging in natural biomes. Metagenomic analyses, single-particle sequencing, and database mining efforts have yielded new sequence data on an astounding number of previously unknown viruses. As metagenomes are relatively free of biases, these data provide an unprecedented insight into the vastness of the virosphere, but to properly value the extent of this diversity it is critical that the viruses are taxonomically classified. Inclusion of uncultivated viruses has already improved the process as well as the understanding of the taxa, viruses, and their evolutionary relationships. The continuous development and testing of computational tools will be required to maintain a dynamic virus taxonomy that can accommodate the new discoveries.

Unexpected myriad of co-occurring viral strains and species in one of the most abundant and microdiverse viruses on Earth

Viral genetic microdiversity drives adaptation, pathogenicity, and speciation and has critical consequences for the viral-host arms race occurring at the strain and species levels, which ultimately impact microbial community structure and biogeochemical cycles. Despite the fact that most efforts have focused on viral macrodiversity, little is known about the microdiversity of ecologically important viruses on Earth. Recently, single-virus genomics discovered the putatively most abundant ocean virus in temperate and tropical waters: the uncultured dsDNA virus vSAG 37-F6 infecting Pelagibacter, the most abundant marine bacteria. In this study, we report the cooccurrence of up to ≈1,500 different viral strains (>95% nucleotide identity) and ≈30 related species (80-95% nucleotide identity) in a single oceanic sample. Viral microdiversity was maintained over space and time, and most alleles were the result of synonymous mutations without any apparent adaptive benefits to cope with host translation codon bias and efficiency. Gene flow analysis used to delimitate species according to the biological species concept (BSC) revealed the impact of recombination in shaping vSAG 37-F6 virus and Pelagibacter speciation. Data demonstrated that this large viral microdiversity somehow mirrors the host species diversity since ≈50% of the 926 analyzed Pelagibacter genomes were found to belong to independent BSC species that do not significantly engage in gene flow with one another. The host range of this evolutionarily successful virus revealed that a single viral species can infect multiple Pelagibacter BSC species, indicating that this virus crosses not only formal BSC barriers but also biomes since viral ancestors are found in freshwater.

Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes

The archaeal tailed viruses (arTV), evolutionarily related to tailed double-stranded DNA (dsDNA) bacteriophages of the class Caudoviricetes, represent the most common isolates infecting halophilic archaea. Only a handful of these viruses have been genomically characterized, limiting our appreciation of their ecological impacts and evolution. Here, we present 37 new genomes of haloarchaeal tailed virus isolates, more than doubling the current number of sequenced arTVs. Analysis of all 63 available complete genomes of arTVs, which we propose to classify into 14 new families and 3 orders, suggests ancient divergence of archaeal and bacterial tailed viruses and points to an extensive sharing of genes involved in DNA metabolism and counterdefense mechanisms, illuminating common strategies of virus-host interactions with tailed bacteriophages. Coupling of the comparative genomics with the host range analysis on a broad panel of haloarchaeal species uncovered 4 distinct groups of viral tail fiber adhesins controlling the host range expansion. The survey of metagenomes using viral hallmark genes suggests that the global architecture of the arTV community is shaped through recurrent transfers between different biomes, including hypersaline, marine, and anoxic environments.

Viruses infecting a warm water picoeukaryote shed light on spatial co-occurrence dynamics of marine viruses and their hosts

The marine picoeukaryote Bathycoccus prasinos has been considered a cosmopolitan alga, although recent studies indicate two ecotypes exist, Clade BI (B. prasinos) and Clade BII. Viruses that infect Bathycoccus Clade BI are known (BpVs), but not that infect BII. We isolated three dsDNA prasinoviruses from the Sargasso Sea against Clade BII isolate RCC716. The BII-Vs do not infect BI, and two (BII-V2 and BII-V3) have larger genomes (~210 kb) than BI-Viruses and BII-V1. BII-Vs share ~90% of their proteins, and between 65% to 83% of their proteins with sequenced BpVs. Phylogenomic reconstructions and PolB analyses establish close-relatedness of BII-V2 and BII-V3, yet BII-V2 has 10-fold higher infectivity and induces greater mortality on host isolate RCC716. BII-V1 is more distant, has a shorter latent period, and infects both available BII isolates, RCC716 and RCC715, while BII-V2 and BII-V3 do not exhibit productive infection of the latter in our experiments. Global metagenome analyses show Clade BI and BII algal relative abundances correlate positively with their respective viruses. The distributions delineate BI/BpVs as occupying lower temperature mesotrophic and coastal systems, whereas BII/BII-Vs occupy warmer temperature, higher salinity ecosystems. Accordingly, with molecular diagnostic support, we name Clade BII Bathycoccus calidus sp. nov. and propose that molecular diversity within this new species likely connects to the differentiated host-virus dynamics observed in our time course experiments. Overall, the tightly linked biogeography of Bathycoccus host and virus clades observed herein supports species-level host specificity, with strain-level variations in infection parameters.

Genomic Characterization of Two Novel RCA Phages Reveals New Insights into the Diversity and Evolution of Marine Viruses

Viruses are the most abundant living entities in marine ecosystems, playing critical roles in altering the structure and function of microbial communities and driving ocean biogeochemistry. Phages that infect Roseobacter clade-affiliated (RCA) cluster strains are an important component of marine viral communities. Here, we characterize the genome sequences of two new RCA phages, CRP-9 and CRP-13, which infect RCA strain FZCC0023. Genomic analysis reveals that CRP-9 and CRP-13 represent a novel evolutionary lineage of marine phages. They both have a DNA replication module most similar to those in Cobavirus group phages. In contrast, their morphogenesis and packaging modules are distinct from those in cobaviruses but homologous to those in HMO-2011-type phages. The genomic architecture of CRP-9 and CRP-13 suggests a genomic recombination event between distinct phage groups. Metagenomic data sets were examined for metagenome-assembled viral genomes (MAVGs) with similar recombinant genome architectures. Fifteen CRP-9-type MAVGs were identified from marine viromes. Additionally, 158 MAVGs were identified containing HMO-2011-type morphogenesis and packaging modules with other types of DNA replication genes, providing more evidence that recombination between different phage groups is a major driver of phage evolution. Altogether, this study significantly expands the understanding of diversity and evolution of marine roseophages. Meanwhile, the analysis of these novel RCA phages and MAVGs highlights the critical role of recombination in shaping phage diversity. These phage sequences are valuable resources for inferring the evolutionary connection of distinct phage groups. IMPORTANCE Diversity and evolution of phages that infect the relatively slow-growing but dominant Roseobacter lineages are largely unknown. In this study, RCA phages CRP-9 and CRP-13 have been isolated on a Roseobacter RCA strain and shown to have a unique genomic architecture, which appears to be the result of a recombination event. CRP-9 and CRP-13 have a DNA replication module most similar to those in Cobavirus group phages and morphogenesis and packaging modules most similar to those in HMO-2011-type phages. HMO-2011-type morphogenesis and packaging modules are found in combination with distinct types of DNA replication genes, suggesting compatibility with various DNA replication modules. Altogether, this study contributes toward a better understanding of marine viral diversity and evolution.

Viral Characteristics of the Warm Atlantic and Cold Arctic Water Masses in the Nordic Seas

Nordic Seas are the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean with complex water masses, experiencing an abrupt climate change. Though knowledge of the marine virosphere has expanded rapidly, the diversity of viruses and their relationships with host cells and water masses in the Nordic Seas remain to be fully revealed. Here, we establish the Nordic Sea DNA virome (NSV) data set of 55,315 viral contigs including 1,478 unique viral populations from seven stations influenced by both the warm Atlantic and cold Arctic water masses. Caudovirales dominated in the seven NSVs, especially in the warm Atlantic waters. The major giant nucleocytoplasmic large DNA viruses (NCLDVs) contributed a significant proportion of the classified viral contigs in the NSVs (32.2%), especially in the cold Arctic waters (44.9%). The distribution patterns of Caudovirales and NCLDVs were a reflection of the community structure of their hosts in the corresponding water masses and currents. Latitude, pH, and flow speed were found to be key factors influencing the microbial communities and coinfluencing the variation of viral communities. Network analysis illustrated the tight coupling between the variation of viral communities and microbial communities in the Nordic Seas. This study suggests a probable linkage between viromes, host cells, and surface water masses from both the cool Arctic and warm Atlantic Oceans. IMPORTANCE This is a systematic study of Nordic Sea viromes using metagenomic analysis. The viral diversity, community structure, and their relationships with host cells and the complex water masses from both the cool Arctic and the warm Atlantic oceans were illustrated. The NCLDVs and Caudovirales are proposed as the viral characteristics of the cold Arctic and warm Atlantic waters, respectively. This study provides an important background for the viromes in the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean and sheds light on their responses to abrupt climate change in the future.

Imaging Techniques for Detecting Prokaryotic Viruses in Environmental Samples

Viruses are the most abundant biological entities on Earth with an estimate of 1031 viral particles across all ecosystems. Prokaryotic viruses-bacteriophages and archaeal viruses-influence global biogeochemical cycles by shaping microbial communities through predation, through the effect of horizontal gene transfer on the host genome evolution, and through manipulating the host cellular metabolism. Imaging techniques have played an important role in understanding the biology and lifestyle of prokaryotic viruses. Specifically, structure-resolving microscopy methods, for example, transmission electron microscopy, are commonly used for understanding viral morphology, ultrastructure, and host interaction. These methods have been applied mostly to cultivated phage-host pairs. However, recent advances in environmental genomics have demonstrated that the majority of viruses remain uncultivated, and thus microscopically uncharacterized. Although light- and structure-resolving microscopy of viruses from environmental samples is possible, quite often the link between the visualization and the genomic information of uncultivated prokaryotic viruses is missing. In this minireview, we summarize the current state of the art of imaging techniques available for characterizing viruses in environmental samples and discuss potential links between viral imaging and environmental genomics for shedding light on the morphology of uncultivated viruses and their lifestyles in Earth's ecosystems.

Bottom-up controls, ecological revolutions and diversification in the oceans through time

Animals originated in the oceans and evolved there for hundreds of millions of years before adapting to terrestrial environments. Today, oceans cover more than two-thirds of Earth and generate as much primary production as land. The path from the first macrobiota to modern marine biodiversity involved parallel increases in terrestrial nutrient input, marine primary production, species' abundance, metabolic rates, ecotypic diversity and taxonomic diversity. Bottom-up theories of ecosystem cascades arrange these changes in a causal sequence. At the base of marine food webs, nutrient fluxes and atmosphere-ocean chemistry interact with phytoplankton to regulate production. First-order consumers (e.g., zooplankton) might propagate changes in quantity and quality of phytoplankton to changes in abundance and diversity of larger predators (e.g., nekton). However, many uncertainties remain about the mechanisms and effect size of bottom-up control, particularly in oceans across the entire history of animal life. Here, we review modern and fossil evidence for hypothesized bottom-up pathways, and we assess the ramifications of these processes for four key intervals in marine ecosystems: the Ediacaran-Cambrian (635-485 million years ago), the Ordovician (485-444 million years ago), the Devonian (419-359 million years ago) and the Mesozoic (252-66 million years ago). We advocate for a clear articulation of bottom-up hypotheses to better understand causal relationships and proposed effects, combined with additional ecological experiments, paleontological documentation, isotope geochemistry and geophysical reconstructions. How small-scale ecological change transitions into large-scale evolutionary change remains an outstanding question for empirical and theoretical research.

Taxonomy-aware, sequence similarity ranking reliably predicts phage-host relationships

Background

Characterizing phage–host interactions is critical to understanding the ecological role of both partners and effective isolation of phage therapeuticals. Unfortunately, experimental methods for studying these interactions are markedly slow, low-throughput, and unsuitable for phages or hosts difficult to maintain in laboratory conditions. Therefore, a number of in silico methods emerged to predict prokaryotic hosts based on viral sequences. One of the leading approaches is the application of the BLAST tool that searches for local similarities between viral and microbial genomes. However, this prediction method has three major limitations: (i) top-scoring sequences do not always point to the actual host; (ii) mosaic virus genomes may match to many, typically related, bacteria; and (iii) viral and host sequences may diverge beyond the point where their relationship can be detected by a BLAST alignment.

Results

We created an extension to BLAST, named Phirbo, that improves host prediction quality beyond what is obtainable from standard BLAST searches. The tool harnesses information concerning sequence similarity and bacteria relatedness to predict phage–host interactions. Phirbo was evaluated on three benchmark sets of known virus–host pairs, and it improved precision and recall by 11–40 percentage points over currently available, state-of-the-art, alignment-based, alignment-free, and machine-learning host prediction tools. Moreover, the discriminatory power of Phirbo for the recognition of virus–host relationships surpassed the results of other tools by at least 10 percentage points (area under the curve = 0.95), yielding a mean host prediction accuracy of 57% and 68% at the genus and family levels, respectively, and drops by 12 percentage points when using only a fraction of viral genome sequences (3 kb). Finally, we provide insights into a repertoire of protein and ncRNA genes that are shared between phages and hosts and may be prone to horizontal transfer during infection.

Conclusions

Our results suggest that Phirbo is a simple and effective tool for predicting phage–host relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01146-6.

Diversity and distribution of viruses inhabiting the deepest ocean on Earth

As the most abundant biological entities on the planet, viruses significantly influence the overall functioning of marine ecosystems. The abundance, distribution, and biodiversity of viral communities in the upper ocean have been relatively well studied, but our understanding of viruses in the hadal biosphere remains poor. Here, we established the oceanic trench viral genome dataset (OTVGD) by analysing 19 microbial metagenomes derived from seawater and sediment samples of the Mariana, Yap, and Kermadec Trenches. The trench viral communities harbored remarkably high novelty, and they were predicted to infect ecologically important microbial clades, including Thaumarchaeota and Oleibacter. Significant inter-trench and intra-trench exchange of viral communities was proposed. Moreover, viral communities in different habitats (seawater/sediment and depth-stratified ocean zones) exhibited distinct niche-dependent distribution patterns and genomic properties. Notably, microbes and viruses in the hadopelagic seawater seemed to preferably adopt lysogenic lifestyles compared to those in the upper ocean. Furthermore, niche-specific auxiliary metabolic genes were identified in the hadal viral genomes, and a novel viral D-amino acid oxidase was functionally and phylogenetically characterized, suggesting the contribution of these genes in the utilization of refractory organic matter. Together, these findings highlight the genomic novelty, dynamic movement, and environment-driven diversification of viral communities in oceanic trenches, and suggest that viruses may influence the hadal ecosystem by reprogramming the metabolism of their hosts and modulating the community of keystone microbes.

Prophage Tracer: precisely tracing prophages in prokaryotic genomes using overlapping split-read alignment

The life cycle of temperate phages includes a lysogenic cycle stage when the phage integrates into the host genome and becomes a prophage. However, the identification of prophages that are highly divergent from known phages remains challenging. In this study, by taking advantage of the lysis-lysogeny switch of temperate phages, we designed Prophage Tracer, a tool for recognizing active prophages in prokaryotic genomes using short-read sequencing data, independent of phage gene similarity searching. Prophage Tracer uses the criterion of overlapping split-read alignment to recognize discriminative reads that contain bacterial (attB) and phage (attP) att sites representing prophage excision signals. Performance testing showed that Prophage Tracer could predict known prophages with precise boundaries, as well as novel prophages. Two novel prophages, dsDNA and ssDNA, encoding highly divergent major capsid proteins, were identified in coral-associated bacteria. Prophage Tracer is a reliable data mining tool for the identification of novel temperate phages and mobile genetic elements. The code for the Prophage Tracer is publicly available at https://github.com/WangLab-SCSIO/Prophage_Tracer.

High-throughput identification of viral termini and packaging mechanisms in virome datasets using PhageTermVirome

Viruses that infect bacteria (phages) are increasingly recognized for their importance in diverse ecosystems but identifying and annotating them in large-scale sequence datasets is still challenging. Although efficient scalable virus identification tools are emerging, defining the exact ends (termini) of phage genomes is still particularly difficult. The proper identification of termini is crucial, as it helps in characterizing the packaging mechanism of bacteriophages and provides information on various aspects of phage biology. Here, we introduce PhageTermVirome (PTV) as a tool for the easy and rapid high-throughput determination of phage termini and packaging mechanisms using modern large-scale metagenomics datasets. We successfully tested the PTV algorithm on a mock virome dataset and then used it on two real virome datasets to achieve the rapid identification of more than 100 phage termini and packaging mechanisms, with just a few hours of computing time. Because PTV allows the identification of free fully formed viral particles (by recognition of termini present only in encapsidated DNA), it can also complement other virus identification softwares to predict the true viral origin of contigs in viral metagenomics datasets. PTV is a novel and unique tool for high-throughput characterization of phage genomes, including phage termini identification and characterization of genome packaging mechanisms. This software should help researchers better visualize, map and study the virosphere. PTV is freely available for downloading and installation at https://gitlab.pasteur.fr/vlegrand/ptv.

Tracing the fate of wastewater viruses reveals catchment-scale virome diversity and connectivity

The discharge of wastewater-derived viruses in aquatic environments impacts catchment-scale virome composition. To explore this, we used viromic analysis of RNA and DNA virus-like particles to holistically track virus communities entering and leaving wastewater treatment plants and the connecting river catchment system and estuary. We reconstructed >40 000 partial viral genomes into 10 149 species-level groups, dominated by dsDNA and (+)ssRNA bacteriophages (Caudoviricetes and Leviviricetes) and a small number of genomes that could pose a risk to human health. We found substantial viral diversity and geographically distinct virus communities associated with different wastewater treatment plants. River and estuarine water bodies harboured more diverse viral communities in downstream locations, influenced by tidal movement and proximity to wastewater treatment plants. Shellfish and beach sand were enriched in viral communities when compared with the surrounding water, acting as entrapment matrices for virus particles. Extensive phylogenetic analyses of environmental-derived and reference sequences showed the presence of human-associated sapovirus GII in all sample types, multiple rotavirus A strains in wastewater and a diverse set of picorna-like viruses associated with shellfish. We conclude that wastewater-derived viral genetic material is commonly deposited in the environment and can be traced throughout the freshwater-marine continuum of the river catchment, where it is influenced by local geography, weather events and tidal effects. Our data illustrate the utility of viromic analyses for wastewater- and environment-based ecology and epidemiology, and we present a conceptual model for the circulation of all types of viruses in a freshwater catchment.

Finding functional associations between prokaryotic virus orthologous groups: a proof of concept

Background

The field of viromics has greatly benefited from recent developments in metagenomics, with significant efforts focusing on viral discovery. However, functional annotation of the increasing number of viral genomes is lagging behind. This is highlighted by the degree of annotation of the protein clusters in the prokaryotic Virus Orthologous Groups (pVOGs) database, with 83% of its current 9518 pVOGs having an unknown function.

Results

In this study we describe a machine learning approach to explore potential functional associations between pVOGs. We measure seven genomic features and use them as input to a Random Forest classifier to predict protein–protein interactions between pairs of pVOGs. After systematic evaluation of the model’s performance on 10 different datasets, we obtained a predictor with a mean accuracy of 0.77 and Area Under Receiving Operation Characteristic (AUROC) score of 0.83. Its application to a set of 2,133,027 pVOG-pVOG interactions allowed us to predict 267,265 putative interactions with a reported probability greater than 0.65. At an expected false discovery rate of 0.27, we placed 95.6% of the previously unannotated pVOGs in a functional context, by predicting their interaction with a pVOG that is functionally annotated.

Conclusions

We believe that this proof-of-concept methodology, wrapped in a reproducible and automated workflow, can represent a significant step towards obtaining a more complete picture of bacteriophage biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04343-w.

Prokaryotic viruses impact functional microorganisms in nutrient removal and carbon cycle in wastewater treatment plants

As one of the largest biotechnological applications, activated sludge (AS) systems in wastewater treatment plants (WWTPs) harbor enormous viruses, with 10-1,000-fold higher concentrations than in natural environments. However, the compositional variation and host-connections of AS viruses remain poorly explored. Here, we report a catalogue of ~50,000 prokaryotic viruses from six WWTPs, increasing the number of described viral species of AS by 23-fold, and showing the very high viral diversity which is largely unknown (98.4-99.6% of total viral contigs). Most viral genera are represented in more than one AS system with 53 identified across all. Viral infection widely spans 8 archaeal and 58 bacterial phyla, linking viruses with aerobic/anaerobic heterotrophs, and other functional microorganisms controlling nitrogen/phosphorous removal. Notably, Mycobacterium, notorious for causing AS foaming, is associated with 402 viral genera. Our findings expand the current AS virus catalogue and provide reference for the phage treatment to control undesired microorganisms in WWTPs.

Circoviruses and cycloviruses identified in Weddell seal fecal samples from McMurdo Sound, Antarctica

Circoviridae is a family of circular single-stranded DNA viruses whose members infect a wide variety of hosts. While well characterized in avian and mammalian hosts, little is known about circoviruses associated with Antarctic animals. From 48 Weddell seal (Leptonychotes weddellii) fecal samples collected on the sea ice in McMurdo between Nov 2014 and Dec 2014, we identified and determined the genomes of novel viruses that fall within two genera of the family Circoviridae, i.e. Circovirus (n = 7) and Cyclovirus (n = 45). We named these viruses as werosea circovirus (WerCV) and werosea cyclovirus (WerCyV). The genomes of WerCV and WerCyV share ~63-64% genome-wide pairwise identity with classified circoviruses and cycloviruses, respectively. Based on the species demarcation threshold of 80% for members of the Circoviridae, the genomes of WerCV and WerCyV represent new species in their respective genera. Evidence indicated recombination in five of the 45 WerCyV genomes identified in this study. These are the first circoviruses found associated with Antarctic pinnipeds, adding to those recently identified associated with Adélie (Pygoscelis adeliae) and chinstrap penguins (P. antarcticus).

Arms race in a cell: genomic, transcriptomic, and proteomic insights into intracellular phage-bacteria interplay in deep-sea snail holobionts

BACKGROUND

Deep-sea animals in hydrothermal vents often form endosymbioses with chemosynthetic bacteria. Endosymbionts serve essential biochemical and ecological functions, but the prokaryotic viruses (phages) that determine their fate are unknown.

RESULTS

We conducted metagenomic analysis of a deep-sea vent snail. We assembled four genome bins for Caudovirales phages that had developed dual endosymbiosis with sulphur-oxidising bacteria (SOB) and methane-oxidising bacteria (MOB). Clustered regularly interspaced short palindromic repeat (CRISPR) spacer mapping, genome comparison, and transcriptomic profiling revealed that phages Bin1, Bin2, and Bin4 infected SOB and MOB. The observation of prophages in the snail endosymbionts and expression of the phage integrase gene suggested the presence of lysogenic infection, and the expression of phage structural protein and lysozyme genes indicated active lytic infection. Furthermore, SOB and MOB appear to employ adaptive CRISPR-Cas systems to target phage DNA. Additional expressed defence systems, such as innate restriction-modification systems and dormancy-inducing toxin-antitoxin systems, may co-function and form multiple lines for anti-viral defence. To counter host defence, phages Bin1, Bin2, and Bin3 appear to have evolved anti-restriction mechanisms and expressed methyltransferase genes that potentially counterbalance host restriction activity. In addition, the high-level expression of the auxiliary metabolic genes narGH, which encode nitrate reductase subunits, may promote ATP production, thereby benefiting phage DNA packaging for replication.

CONCLUSIONS

This study provides new insights into phage-bacteria interplay in intracellular environments of a deep-sea vent snail. Video Abstract.

Ecogenomics and Adaptation Strategies of Southern Ocean Viral Communities

The Southern Ocean (SO) represents up to one-fifth of the total carbon drawdown worldwide. Intense selective pressures (low temperature, high UV radiation, and strong seasonality) and physical isolation characterize the SO, serving as a "natural" laboratory for the study of ecogenomics and unique adaptations of endemic viral populations. Here, we report 2,416 novel viral genomes from the SO, obtained from newly sequenced viral metagenomes in combination with mining of publicly available data sets, which represents a 25% increase in the SO viral genomes reported to date. They comprised 567 viral clusters (defined as approximately genus-level groups), with 186 genera endemic to the SO, demonstrating that the SO viral community is predominantly constituted by a large pool of genetically divergent viral species from widespread viral families. The predicted proteome from SO viruses revealed that several protein clusters related to cold-shock-event responses and quorum-sensing mechanisms involved in the lysogenic-lytic cycle shift decision were under positive selection, which is ultimately important for fine adaptation of viral populations in response to the strong selective pressures of the SO. Finally, changes in the hydrophobicity patterns and amino acid frequencies suggested marked temperature-driven genetic selection of the SO viral proteome. Our data provide valuable insights into how viruses adapt and remain successful in this extreme polar marine environment. IMPORTANCE Viruses are the most abundant biologic entities in marine systems and strongly influence the microbial community composition and diversity. However, little is known about viral communities' adaptation and diversification in the ocean. In this work, we take advantage of the geographical isolation and the intense selective pressures of the SO, to which viruses are exposed, to identify potential viral adaptations due to positive environmental selection and dispersal limitation. To that end, we recovered more than two thousand novel viral genomes, revealing a high degree of divergence in these SO endemic communities. Furthermore, we describe remarkable viral adaptations in amino acid frequencies and accessory proteins related to cold shock response and quorum sensing that allow them to thrive at lower temperatures. Consequently, our work greatly expands the understanding of the diversification of the viral communities of the SO and their particular adaptations to low temperatures.

The Influence of Habitat on Viral Diversity in Neotropical Rodent Hosts

Rodents are important reservoirs of numerous viruses, some of which have significant impacts on public health. Ecosystem disturbances and decreased host species richness have been associated with the emergence of zoonotic diseases. In this study, we aimed at (a) characterizing the viral diversity in seven neotropical rodent species living in four types of habitats and (b) exploring how the extent of environmental disturbance influences this diversity. Through a metagenomic approach, we identified 77,767 viral sequences from spleen, kidney, and serum samples. These viral sequences were attributed to 27 viral families known to infect vertebrates, invertebrates, plants, and amoeba. Viral diversities were greater in pristine habitats compared with disturbed ones, and lowest in peri-urban areas. High viral richness was observed in savannah areas. Differences in these diversities were explained by rare viruses that were generally more frequent in pristine forest and savannah habitats. Moreover, changes in the ecology and behavior of rodent hosts, in a given habitat, such as modifications to the diet in disturbed vs. pristine forests, are major determinants of viral composition. Lastly, the phylogenetic relationships of four vertebrate-related viral families (Polyomaviridae, Flaviviridae, Togaviridae, and Phenuiviridae) highlighted the wide diversity of these viral families, and in some cases, a potential risk of transmission to humans. All these findings provide significant insights into the diversity of rodent viruses in Amazonia, and emphasize that habitats and the host’s dietary ecology may drive viral diversity. Linking viral richness and abundance to the ecology of their hosts and their responses to habitat disturbance could be the starting point for a better understanding of viral emergence and for future management of ecosystems.

Interactions between bacterial and phage communities in natural environments

We commonly acknowledge that bacterial viruses (phages) shape the composition and evolution of bacterial communities in nature and therefore have important roles in ecosystem functioning. This view stems from studies in the 1990s to the first decade of the twenty-first century that revealed high viral abundance, high viral diversity and virus-induced microbial death in aquatic ecosystems as well as an association between collapses in bacterial density and peaks in phage abundance. The recent surge in metagenomic analyses has provided deeper insight into the abundance, genomic diversity and spatio-temporal dynamics of phages in a wide variety of ecosystems, ranging from deep oceans to soil and the mammalian digestive tract. However, the causes and consequences of variations in phage community compositions remain poorly understood. In this Review, we explore current knowledge of the composition and evolution of phage communities, as well as their roles in controlling the population and evolutionary dynamics of bacterial communities. We discuss the need for greater ecological realism in laboratory studies to capture the complexity of microbial communities that thrive in natural environments.

Hostile Takeover: How Viruses Reprogram Prokaryotic Metabolism

To reproduce, prokaryotic viruses must hijack the cellular machinery of their hosts and redirect it toward the production of viral particles. While takeover of the host replication and protein synthesis apparatus has long been considered an essential feature of infection, recent studies indicate that extensive reprogramming of host primary metabolism is a widespread phenomenon among prokaryotic viruses that is required to fulfill the biosynthetic needs of virion production. In this review we provide an overview of the most significant recent findings regarding virus-induced reprogramming of prokaryotic metabolism and suggest how quantitative systems biology approaches may be used to provide a holistic understanding of metabolic remodeling during lytic viral infection. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Deep sea sediments associated with cold seeps are a subsurface reservoir of viral diversity

In marine ecosystems, viruses exert control on the composition and metabolism of microbial communities, influencing overall biogeochemical cycling. Deep sea sediments associated with cold seeps are known to host taxonomically diverse microbial communities, but little is known about viruses infecting these microorganisms. Here, we probed metagenomes from seven geographically diverse cold seeps across global oceans to assess viral diversity, virus-host interaction, and virus-encoded auxiliary metabolic genes (AMGs). Gene-sharing network comparisons with viruses inhabiting other ecosystems reveal that cold seep sediments harbour considerable unexplored viral diversity. Most cold seep viruses display high degrees of endemism with seep fluid flux being one of the main drivers of viral community composition. In silico predictions linked 14.2% of the viruses to microbial host populations with many belonging to poorly understood candidate bacterial and archaeal phyla. Lysis was predicted to be a predominant viral lifestyle based on lineage-specific virus/host abundance ratios. Metabolic predictions of prokaryotic host genomes and viral AMGs suggest that viruses influence microbial hydrocarbon biodegradation at cold seeps, as well as other carbon, sulfur and nitrogen cycling via virus-induced mortality and/or metabolic augmentation. Overall, these findings reveal the global diversity and biogeography of cold seep viruses and indicate how viruses may manipulate seep microbial ecology and biogeochemistry.

Virus-associated organosulfur metabolism in human and environmental systems

Viruses influence the fate of nutrients and human health by killing microorganisms and altering metabolic processes. Organosulfur metabolism and biologically derived hydrogen sulfide play dynamic roles in manifestation of diseases, infrastructure degradation, and essential biological processes. Although microbial organosulfur metabolism is well studied, the role of viruses in organosulfur metabolism is unknown. Here, we report the discovery of 39 gene families involved in organosulfur metabolism encoded by 3,749 viruses from diverse ecosystems, including human microbiomes. The viruses infect organisms from all three domains of life. Six gene families encode for enzymes that degrade organosulfur compounds into sulfide, whereas others manipulate organosulfur compounds and may influence sulfide production. We show that viral metabolic genes encode key enzymatic domains, are translated into protein, and are maintained after recombination, and sulfide provides a fitness advantage to viruses. Our results reveal viruses as drivers of organosulfur metabolism with important implications for human and environmental health.

Glacier ice archives nearly 15,000-year-old microbes and phages

BACKGROUND

Glacier ice archives information, including microbiology, that helps reveal paleoclimate histories and predict future climate change. Though glacier-ice microbes are studied using culture or amplicon approaches, more challenging metagenomic approaches, which provide access to functional, genome-resolved information and viruses, are under-utilized, partly due to low biomass and potential contamination.

RESULTS

We expand existing clean sampling procedures using controlled artificial ice-core experiments and adapted previously established low-biomass metagenomic approaches to study glacier-ice viruses. Controlled sampling experiments drastically reduced mock contaminants including bacteria, viruses, and free DNA to background levels. Amplicon sequencing from eight depths of two Tibetan Plateau ice cores revealed common glacier-ice lineages including Janthinobacterium, Polaromonas, Herminiimonas, Flavobacterium, Sphingomonas, and Methylobacterium as the dominant genera, while microbial communities were significantly different between two ice cores, associating with different climate conditions during deposition. Separately, ~355- and ~14,400-year-old ice were subject to viral enrichment and low-input quantitative sequencing, yielding genomic sequences for 33 vOTUs. These were virtually all unique to this study, representing 28 novel genera and not a single species shared with 225 environmentally diverse viromes. Further, 42.4% of the vOTUs were identifiable temperate, which is significantly higher than that in gut, soil, and marine viromes, and indicates that temperate phages are possibly favored in glacier-ice environments before being frozen. In silico host predictions linked 18 vOTUs to co-occurring abundant bacteria (Methylobacterium, Sphingomonas, and Janthinobacterium), indicating that these phages infected ice-abundant bacterial groups before being archived. Functional genome annotation revealed four virus-encoded auxiliary metabolic genes, particularly two motility genes suggest viruses potentially facilitate nutrient acquisition for their hosts. Finally, given their possible importance to methane cycling in ice, we focused on Methylobacterium viruses by contextualizing our ice-observed viruses against 123 viromes and prophages extracted from 131 Methylobacterium genomes, revealing that the archived viruses might originate from soil or plants.

CONCLUSIONS

Together, these efforts further microbial and viral sampling procedures for glacier ice and provide a first window into viral communities and functions in ancient glacier environments. Such methods and datasets can potentially enable researchers to contextualize new discoveries and begin to incorporate glacier-ice microbes and their viruses relative to past and present climate change in geographically diverse regions globally. Video Abstract.

Revealing the full biosphere structure and versatile metabolic functions in the deepest ocean sediment of the Challenger Deep

BACKGROUND

The full biosphere structure and functional exploration of the microbial communities of the Challenger Deep of the Mariana Trench, the deepest known hadal zone on Earth, lag far behind that of other marine realms.

RESULTS

We adopt a deep metagenomics approach to investigate the microbiome in the sediment of Challenger Deep, Mariana Trench. We construct 178 metagenome-assembled genomes (MAGs) representing 26 phyla, 16 of which are reported from hadal sediment for the first time. Based on the MAGs, we find the microbial community functions are marked by enrichment and prevalence of mixotrophy and facultative anaerobic metabolism. The microeukaryotic community is found to be dominated by six fungal groups that are characterized for the first time in hadal sediment to possess the assimilatory and dissimilatory nitrate/sulfate reduction, and hydrogen sulfide oxidation pathways. By metaviromic analysis, we reveal novel hadal Caudovirales clades, distinctive virus-host interactions, and specialized auxiliary metabolic genes for modulating hosts' nitrogen/sulfur metabolism. The hadal microbiome is further investigated by large-scale cultivation that cataloged 1070 bacterial and 19 fungal isolates from the Challenger Deep sediment, many of which are found to be new species specialized in the hadal habitat.

CONCLUSION

Our hadal MAGs and isolates increase the diversity of the Challenger Deep sediment microbial genomes and isolates present in the public. The deep metagenomics approach fills the knowledge gaps in structure and diversity of the hadal microbiome, and provides novel insight into the ecology and metabolism of eukaryotic and viral components in the deepest biosphere on earth.

The Epistemic Revolution Induced by Microbiome Studies: An Interdisciplinary View

Simple Summary

This interdisciplinary study, conducted by experts in evolutionary biology, ecology, ecosystem studies, arts, medicine, forensic analyses, agriculture, law, and philosophy of science describe how microbiome studies are convergently affecting the concepts and practices of diverse fields and practices, that now consider microbiomes within their legitimate scope. Consequently, it describes what seems to be an ongoing pluridisciplinary epistemic revolution, with the potential to fundamentally change how we understand the world through an ecologization of pre-existing concepts, a greater focus on interactions, the use of multi-scalar interaction networks as explanatory frameworks, the reconceptualization of the usual definitions of individuals, and a de-anthropocentrification of our perception of phenomena.

Abstract

Many separate fields and practices nowadays consider microbes as part of their legitimate focus. Therefore, microbiome studies may act as unexpected unifying forces across very different disciplines. Here, we summarize how microbiomes appear as novel major biological players, offer new artistic frontiers, new uses from medicine to laws, and inspire novel ontologies. We identify several convergent emerging themes across ecosystem studies, microbial and evolutionary ecology, arts, medicine, forensic analyses, law and philosophy of science, as well as some outstanding issues raised by microbiome studies across these disciplines and practices. An ‘epistemic revolution induced by microbiome studies’ seems to be ongoing, characterized by four features: (i) an ecologization of pre-existing concepts within disciplines, (ii) a growing interest in systemic analyses of the investigated or represented phenomena and a greater focus on interactions as their root causes, (iii) the intent to use openly multi-scalar interaction networks as an explanatory framework to investigate phenomena to acknowledge the causal effects of microbiomes, (iv) a reconceptualization of the usual definitions of which individuals are worth considering as an explanans or as an explanandum by a given field, which result in a fifth strong trend, namely (v) a de-anthropocentrification of our perception of the world.

Genomic diversity, life strategies and ecology of marine HTVC010P-type pelagiphages

SAR11 bacteria dominate ocean surface bacterioplankton communities, and play an important role in marine carbon and nutrient cycling. The biology and ecology of SAR11 are impacted by SAR11 phages (pelagiphages) that are highly diverse and abundant in the ocean. Among the currently known pelagiphages, HTVC010P represents an extremely abundant but under-studied phage group in the ocean. In this study, we have isolated seven new HTVC010P-type pelagiphages, and recovered 77 nearly full-length HTVC010P-type metagenomic viral genomes from marine metagenomes. Comparative genomic and phylogenomic analyses showed that HTVC010P-type pelagiphages display genome synteny and can be clustered into two major subgroups, with subgroup I consisting of strictly lytic phages and subgroup II mostly consisting of phages with potential lysogenic life cycles. All but one member of the subgroup II contain an integrase gene. Site-specific integration of subgroup II HTVC010P-type pelagiphage was either verified experimentally or identified by in silico genomic sequence analyses, which revealed that various SAR11 tRNA genes can serve as the integration sites of HTVC010P-type pelagiphages. Moreover, HTVC010P-type pelagiphage integration was confirmed by the detection of several Global Ocean Survey (GOS) fragments that contain hybrid phage–host integration sites. Metagenomic recruitment analysis revealed that these HTVC010P-type phages were globally distributed and most lytic subgroup I members exhibited higher relative abundance. Altogether, this study significantly expands our knowledge about the genetic diversity, life strategies and ecology of HTVC010P-type pelagiphages.

Spatial patterns in phage-Rhizobium coevolutionary interactions across regions of common bean domestication

Bacteriophages play significant roles in the composition, diversity, and evolution of bacterial communities. Despite their importance, it remains unclear how phage diversity and phage-host interactions are spatially structured. Local adaptation may play a key role. Nitrogen-fixing symbiotic bacteria, known as rhizobia, have been shown to locally adapt to domesticated common bean at its Mesoamerican and Andean sites of origin. This may affect phage-rhizobium interactions. However, knowledge about the diversity and coevolution of phages with their respective Rhizobium populations is lacking. Here, through the study of four phage-Rhizobium communities in Mexico and Argentina, we show that both phage and host diversity is spatially structured. Cross-infection experiments demonstrated that phage infection rates were higher overall in sympatric rhizobia than in allopatric rhizobia except for one Argentinean community, indicating phage local adaptation and host maladaptation. Phage-host interactions were shaped by the genetic identity and geographic origin of both the phage and the host. The phages ranged from specialists to generalists, revealing a nested network of interactions. Our results suggest a key role of local adaptation to resident host bacterial communities in shaping the phage genetic and phenotypic composition, following a similar spatial pattern of diversity and coevolution to that in the host.

Revealing the Viral Community in the Hadal Sediment of the New Britain Trench

Marine viruses are widely distributed and influence matter and energy transformation in ecosystems by modulating hosts’ metabolism. The hadal trenches represent the deepest marine habitat on Earth, for which the viral communities and related biogeochemical functions are least explored and poorly understood. Here, using the sediment samples (8720 m below sea level) collected from the New Britain Trench (NBT), we investigated the viral community, diversity, and genetic potentials in the hadal sediment habitat for the first time by deep shotgun metagenomic sequencing. We found the NBT sediment viral community was dominated by Siphoviridae, Myoviridae, Podoviridae, Mimiviridae, and Phycodnaviridae, which belong to the dsDNA viruses. However, the large majority of them remained uncharacterized. We found the hadal sediment virome had some common components by comparing the hadal sediment viruses with those of hadal aquatic habitats and those of bathypelagic and terrestrial habitats. It was also distinctive in community structure and had many novel viral clusters not associated with the other habitual virome included in our analyses. Further phylogenetic analysis on its Caudovirales showed novel diversities, including new clades specially evolved in the hadal sediment habitat. Annotation of the NBT sediment viruses indicated the viruses might influence microbial hydrocarbon biodegradation and carbon and sulfur cycling via metabolic augmentation through auxiliary metabolic genes (AMGs). Our study filled in the knowledge gaps on the virome of the hadal sediment habitats and provided insight into the evolution and the potential metabolic functions of the hadal sediment virome.

Ecogenomics of Groundwater Phages Suggests Niche Differentiation Linked to Specific Environmental Tolerance

Viruses are ubiquitous microbiome components, shaping ecosystems via strain-specific predation, horizontal gene transfer and redistribution of nutrients through host lysis. Viral impacts are important in groundwater ecosystems, where microbes drive many nutrient fluxes and metabolic processes; however, little is known about the diversity of viruses in these environments. We analyzed four groundwater plasmidomes (the entire plasmid content of an environment) and identified 200 viral sequences, which clustered into 41 genus-level viral clusters (approximately equivalent to viral genera) including 9 known and 32 putative new genera. We used publicly available bacterial whole-genome sequences (WGS) and WGS from 261 bacterial isolates from this groundwater environment to identify potential viral hosts. We linked 76 of the 200 viral sequences to a range of bacterial phyla, the majority associated with Proteobacteria, followed by Firmicutes, Bacteroidetes, and Actinobacteria. The publicly available WGS enabled mapping bacterial hosts to several viral sequences. The WGS of groundwater isolates increased the depth of host prediction by allowing host identification at the strain level. The latter included 4 viruses that were almost entirely (>99% query coverage, >99% identity) identified as integrated in the genomes of Pseudomonas, Acidovorax, and Castellaniella strains, resulting in high-confidence host assignments. Lastly, 21 of these viruses carried putative auxiliary metabolite genes for metal and antibiotic resistance, which might drive their infection cycles and/or provide selective advantage to infected hosts. Exploring the groundwater virome provides a necessary foundation for integration of viruses into ecosystem models where they are key players in microbial adaption to environmental stress. IMPORTANCE To our knowledge, this is the first study to identify the bacteriophage distribution in a groundwater ecosystem shedding light on their prevalence and distribution across metal-contaminated and background sites. Our study is uniquely based on selective sequencing of solely the extrachromosomal elements of a microbiome followed by analysis for viral signatures, thus establishing a more focused approach for phage identifications. Using this method, we detected several novel phage genera along with those previously established. Our approach of using the whole-genome sequences of hundreds of bacterial isolates from the same site enabled us to make host assignments with high confidence, several at strain levels. Certain phage genes suggest that they provide an environment-specific selective advantage to their bacterial hosts. Our study lays the foundation for future research on directed phage isolations using specific bacterial host strains to further characterize groundwater phages, their life cycles, and their effects on groundwater microbiome and biogeochemistry.

Application of machine learning in bacteriophage research

Phages are one of the key components in the structure, dynamics, and interactions of microbial communities in different bins. It has a clear impact on human health and the food industry. Bacteriophage characterization using in vitro approaches are time/cost consuming and laborious tasks. On the other hand, with the advent of new high-throughput sequencing technology, the development of a powerful computational framework to characterize the newly identified bacteriophages is inevitable for future research. Machine learning includes powerful techniques that enable the analysis of complex datasets for knowledge discovery and pattern recognition. In this study, we have conducted a comprehensive review of machine learning methods application using different types of features were applied in various aspects of bacteriophage research including, automated curation, identification, classification, host species recognition, virion protein identification, and life cycle prediction. Moreover, potential limitations and advantages of the developed frameworks were discussed.

Prospects for viruses infecting eukaryotic microalgae in biotechnology

Besides being considered pathogens, viruses are important drivers of evolution and they can shape large ecological and biogeochemical processes, by influencing host fitness, population dynamics, and community structures. Moreover, they are simple systems that can be used and manipulated to be beneficial and useful for biotechnological applications. In this context, microalgae biotechnology is a growing field of research, which investigated the usage of photosynthetic microorganisms for the sustainable production of food, fuel, chemical, and pharmaceutical sectors. Viruses infecting microalgae have become important subject of ecological studies related to marine and aquatic environments only four decades ago when virus-like-particles associated with bloom-forming algae were discovered. These first findings have opened new questions on evolution and identity. To date, 63 viruses that infect eukaryotic microalgae have been isolated and cultured. In this short review we briefly summarize what is known about viruses infecting eukaryotic microalgae, and how acknowledging their importance can shape future research focussed not only on marine ecology and evolutionary biology but also on biotechnological applications related to microalgae cell factories.

Diverse Viruses in Deep-Sea Hydrothermal Vent Fluids Have Restricted Dispersal across Ocean Basins

In the ocean, viruses impact microbial mortality, regulate biogeochemical cycling, and alter the metabolic potential of microbial lineages. At deep-sea hydrothermal vents, abundant viruses infect a wide range of hosts among the archaea and bacteria that inhabit these dynamic habitats. However, little is known about viral diversity, host range, and biogeography across different vent ecosystems, which has important implications for how viruses manipulate microbial function and evolution. Here, we examined viral diversity, viral and host distribution, and virus-host interactions in microbial metagenomes generated from venting fluids from several vent sites within three different geochemically and geographically distinct hydrothermal systems: Piccard and Von Damm vent fields at the Mid-Cayman Rise in the Caribbean Sea, and at several vent sites within Axial Seamount in the Pacific Ocean. Analysis of viral sequences and clustered regularly interspaced short palindromic repeat (CRISPR) spacers revealed highly diverse viral assemblages and evidence of active infection. Network analysis revealed that viral host range was relatively narrow, with very few viruses infecting multiple microbial lineages. Viruses were largely endemic to individual vent sites, indicating restricted dispersal, and in some cases, viral assemblages persisted over time. Thus, we show that hydrothermal vent fluids are home to novel, diverse viral assemblages that are highly localized to specific regions and taxa. IMPORTANCE Viruses play important roles in manipulating microbial communities and their evolution in the ocean, yet not much is known about viruses in deep-sea hydrothermal vents. However, viral ecology and evolution are of particular interest in hydrothermal vent habitats because of their unique nature: previous studies have indicated that most viruses in hydrothermal vents are temperate rather than lytic, and it has been established that rates of horizontal gene transfer (HGT) are particularly high among thermophilic vent microbes, and viruses are common vectors for HGT. If viruses have broad host range or are widespread across vent sites, they have increased potential to act as gene-sharing "highways" between vent sites. By examining viral diversity, distribution, and infection networks across disparate vent sites, this study provides the opportunity to better characterize and constrain the viral impact on hydrothermal vent microbial communities. We show that viruses in hydrothermal vents are diverse and apparently active, but most have restricted host range and are not widely distributed among vent sites. Thus, the impacts of viral infection are likely to be highly localized and constrained to specific taxa in these habitats.

Viral footprints across Gulfs of Kathiawar Peninsula and Arabian Sea: Unraveled from pelagic sediment metagenomic data

Marine biosphere is one of the largest, diverse and dynamic system hosting numerous of microorganisms. Viruses being the most abundant under explored lifeforms in ocean, represent a reservoir of great genetic diversity. We report the metagenomic insights on the viral communities in the deep sediments of the two Gulfs of Gujarat i.e. Gulf of Khambhat and Gulf of Kutch, with one sample from Arabian Sea, treated as open sea control. The viral reads were filtered from the whole dataset, assembled and studied for viral diversity, which was visualized by Pavian. The sequences were checked for the viral abundance, diversity and functionality. The resulting viral taxonomic classification contained 6 orders, 8 families and 47 genera. The results revealed that the phages infecting Cyanobacterium, Bacillus and Vibrio dominated the sediments. Further, it was observed that majority of viral sequences belonged to double-stranded DNA phages. The present study attempts to provide a primary insight of the viral signals and potential genetic content in the Gulfs of Kathiawar.

Viromes in marine ecosystems reveal remarkable invertebrate RNA virus diversity

Little is known about ocean viromes and the ecological drivers of the evolution of aquatic RNA viruses. This study employed a meta-transcriptomic approach to characterize the viromes of 58 marine invertebrate species across three seas. This revealed the presence of 315 newly identified RNA viruses in nine viral families or orders (Durnavirales, Totiviridae, Bunyavirales, Hantaviridae, Picornavirales, Flaviviridae, Hepelivirales, Solemoviridae, and Tombusviridae), with most of them being sufficiently divergent to the already documented viruses. Notably, this study revealed three marine invertebrate hantaviruses that are rooted to vertebrate hantaviruses, further supporting that hantaviruses may have a marine origin. We have also found evidence for possible host sharing and switch events during virus evolution. Overall, we have revealed the hidden diversity of marine invertebrate RNA viruses.

efam: an expanded, metaproteome-supported HMM profile database of viral protein families

Motivation

Viruses infect, reprogram and kill microbes, leading to profound ecosystem consequences, from elemental cycling in oceans and soils to microbiome-modulated diseases in plants and animals. Although metagenomic datasets are increasingly available, identifying viruses in them is challenging due to poor representation and annotation of viral sequences in databases.

Results

Here, we establish efam, an expanded collection of Hidden Markov Model (HMM) profiles that represent viral protein families conservatively identified from the Global Ocean Virome 2.0 dataset. This resulted in 240 311 HMM profiles, each with at least 2 protein sequences, making efam >7-fold larger than the next largest, pan-ecosystem viral HMM profile database. Adjusting the criteria for viral contig confidence from ‘conservative’ to ‘eXtremely Conservative’ resulted in 37 841 HMM profiles in our efam-XC database. To assess the value of this resource, we integrated efam-XC into VirSorter viral discovery software to discover viruses from less-studied, ecologically distinct oxygen minimum zone (OMZ) marine habitats. This expanded database led to an increase in viruses recovered from every tested OMZ virome by ∼24% on average (up to ∼42%) and especially improved the recovery of often-missed shorter contigs (<5 kb). Additionally, to help elucidate lesser-known viral protein functions, we annotated the profiles using multiple databases from the DRAM pipeline and virion-associated metaproteomic data, which doubled the number of annotations obtainable by standard, single-database annotation approaches. Together, these marine resources (efam and efam-XC) are provided as searchable, compressed HMM databases that will be updated bi-annually to help maximize viral sequence discovery and study from any ecosystem.

Availability and implementation

The resources are available on the iVirus platform at (doi.org/10.25739/9vze-4143).

Supplementary information

Supplementary data are available at Bioinformatics online.