Ocean plankton. Patterns and ecological drivers of ocean viral communities

Spinal Cord Injury Changes the Structure and Functional Potential of Gut Bacterial and Viral Communities

Emerging data indicate that gut dysbiosis contributes to many human diseases, including several comorbidities that develop after traumatic spinal cord injury (SCI). To date, all analyses of SCI-induced gut dysbiosis have used 16S rRNA amplicon sequencing. This technique has several limitations, including being susceptible to taxonomic "blind spots," primer bias, and an inability to profile microbiota functions or identify viruses. Here, SCI-induced gut dysbiosis was assessed by applying genome- and gene-resolved metagenomic analysis of murine stool samples collected 21 days after an experimental SCI at the 4th thoracic spine (T4) or 10th thoracic spine (T10) spinal level. These distinct injuries partially (T10) or completely (T4) abolish sympathetic tone in the gut. Among bacteria, 105 medium- to high-quality metagenome-assembled genomes (MAGs) were recovered, with most (n = 96) representing new bacterial species. Read mapping revealed that after SCI, the relative abundance of beneficial commensals (Lactobacillus johnsonii and CAG-1031 spp.) decreased, while potentially pathogenic bacteria (Weissella cibaria, Lactococcus lactis _A, Bacteroides thetaiotaomicron) increased. Functionally, microbial genes encoding proteins for tryptophan, vitamin B6, and folate biosynthesis, essential pathways for central nervous system function, were reduced after SCI. Among viruses, 1,028 mostly novel viral populations were recovered, expanding known murine gut viral species sequence space ∼3-fold compared to that of public databases. Phages of beneficial commensal hosts (CAG-1031, Lactobacillus, and Turicibacter) decreased, while phages of pathogenic hosts (Weissella, Lactococcus, and class Clostridia) increased after SCI. Although the microbiomes and viromes were changed in all SCI mice, some of these changes varied as a function of spinal injury level, implicating loss of sympathetic tone as a mechanism underlying gut dysbiosis.IMPORTANCE To our knowledge, this is the first article to apply metagenomics to characterize changes in gut microbial population dynamics caused by a clinically relevant model of central nervous system (CNS) trauma. It also utilizes the most current approaches in genome-resolved metagenomics and viromics to maximize the biological inferences that can be made from these data. Overall, this article highlights the importance of autonomic nervous system regulation of a distal organ (gut) and its microbiome inhabitants after traumatic spinal cord injury (SCI). By providing information on taxonomy, function, and viruses, metagenomic data may better predict how SCI-induced gut dysbiosis influences systemic and neurological outcomes after SCI.

Role of Phylogenetic Structure in the Dynamics of Coastal Viral Assemblages

Marine microbes, including viruses, are an essential part of the marine ecosystem, forming the base of the food web and driving biogeochemical cycles. Within this system, the composition of viral assemblages changes markedly with time, and some of these changes are repeatable through time; however, the extent to which these dynamics are reflected within versus among evolutionarily related groups of viruses is largely unexplored. To examine these dynamics, changes in the composition of two groups of ecologically important viruses and communities of their potential hosts were sampled every 2 weeks for 13 months at a coastal site in British Columbia, Canada. We sequenced two marker genes for viruses-the gene encoding the major capsid protein of T4-like phages and their relatives (gp23) and the RNA-dependent RNA polymerase (RdRp) gene of marnavirus-like RNA viruses-as well as marker genes for their bacterial and eukaryotic host communities, the genes encoding 16S rRNA and 18S rRNA. There were strong lagged correlations between viral diversity and community similarity of putative hosts, implying that the viruses influenced the composition of the host communities. The results showed that for both viral assemblages, the dominant clusters of phylogenetically related viruses shifted over time, and this was correlated with environmental changes. Viral clusters contained many ephemeral taxa and few persistent taxa, but within a viral assemblage, the ephemeral and persistent taxa were closely related, implying ecological dynamics within these clusters. Furthermore, these dynamics occurred in both the RNA and DNA viral assemblages surveyed, implying that this structure is common in natural viral assemblages.IMPORTANCE Viruses are major agents of microbial mortality in marine systems, yet little is known about changes in the composition of viral assemblages in relation to those of the microbial communities that they infect. Here, we sampled coastal seawater every 2 weeks for 1 year and used high-throughput sequencing of marker genes to follow changes in the composition of two groups of ecologically important viruses, as well as the communities of bacteria and protists that serve as their respective hosts. Different subsets of genetically related viruses dominated at different times. These results demonstrate that although the genetic composition of viral assemblages is highly dynamic temporally, for the most part the shuffling of genotypes occurs within a few clusters of phylogenetically related viruses. Thus, it appears that even in temperate coastal waters with large seasonal changes, the highly dynamic shuffling of viral genotypes occurs largely within a few subsets of related individuals.

Analysis of Different Size Fractions Provides a More Complete Perspective of Viral Diversity in a Freshwater Embayment

Inspired by recent discoveries of the prevalence of large viruses in the environment, we reassessed the longstanding approach of filtering water through small-pore-size filters to separate viruses from cells before metagenomic analysis. We collected samples from three sites in Hamilton Harbour, an embayment of Lake Ontario, and studied 6 data sets derived from <0.45-μm- and >0.45-μm-size fractions to compare the diversity of viruses in these fractions. At the level of virus order/family, we observed highly diverse and distinct virus communities in the >0.45-μm-size fractions, whereas the <0.45-μm-size fractions were composed primarily of Caudovirales The relative abundances of Caudovirales for which hosts could be inferred varied widely between size fractions, with higher relative abundances of cyanophages in the >0.45-μm-size fractions, potentially indicating replication within cells during ongoing infections. Many viruses of eukaryotes, such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, were detected exclusively in the often-disregarded >0.45-μm-size fractions. In addition to observing unique virus communities associated with each size fraction from every site we examined, we detected viruses common to both fractions, suggesting that these are candidates for further exploration because they could be the product of ongoing or recent lytic events. Most importantly, our observations indicate that analysis of either fraction alone provides only a partial perspective of double-stranded DNA (dsDNA) viruses in the environment, highlighting the need for more comprehensive approaches for analyzing virus communities inferred from metagenomic sequencing.IMPORTANCE Most studies of aquatic virus communities analyze DNA sequences derived from the smaller-size "free-virus" fraction. Our study demonstrates that analysis of virus communities using only the smaller-size fraction can lead to erroneously low diversity estimates for many of the larger viruses such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, whereas analyzing only the larger->0.45-μm-size fraction can lead to underestimates of Caudovirales diversity and relative abundance. Similarly, our data show that examining only the smaller-size fraction can lead to underestimations of virophage and cyanophage relative abundances that could, in turn, cause researchers to assume their limited ecological importance. Given the considerable differences we observed in this study, we recommend cautious interpretations of environmental virus community assemblages and dynamics when based on metagenomic data derived from different size fractions.

BACPHLIP: predicting bacteriophage lifestyle from conserved protein domains

Bacteriophages are broadly classified into two distinct lifestyles: temperate and virulent. Temperate phages are capable of a latent phase of infection within a host cell (lysogenic cycle), whereas virulent phages directly replicate and lyse host cells upon infection (lytic cycle). Accurate lifestyle identification is critical for determining the role of individual phage species within ecosystems and their effect on host evolution. Here, we present BACPHLIP, a BACterioPHage LIfestyle Predictor. BACPHLIP detects the presence of a set of conserved protein domains within an input genome and uses this data to predict lifestyle via a Random Forest classifier that was trained on a dataset of 634 phage genomes. On an independent test set of 423 phages, BACPHLIP has an accuracy of 98% greatly exceeding that of the previously existing tools (79%). BACPHLIP is freely available on GitHub (https://github.com/adamhockenberry/bacphlip) and the code used to build and test the classifier is provided in a separate repository (https://github.com/adamhockenberry/bacphlip-model-dev) for users wishing to interrogate and re-train the underlying classification model.

Comparison of ultrafiltration and iron chloride flocculation in the preparation of aquatic viromes from contrasting sample types

Viral metagenomes (viromes) are a valuable untargeted tool for studying viral diversity and the central roles viruses play in host disease, ecology, and evolution. Establishing effective methods to concentrate and purify viral genomes prior to sequencing is essential for high quality viromes. Using virus spike-and-recovery experiments, we stepwise compared two common approaches for virus concentration, ultrafiltration and iron chloride flocculation, across diverse matrices: wastewater influent, wastewater secondary effluent, river water, and seawater. Viral DNA was purified by removing cellular DNA via chloroform cell lysis, filtration, and enzymatic degradation of extra-viral DNA. We found that viral genomes were concentrated 1-2 orders of magnitude more with ultrafiltration than iron chloride flocculation for all matrices and resulted in higher quality DNA suitable for amplification-free and long-read sequencing. Given its widespread use and utility as an inexpensive field method for virome sampling, we nonetheless sought to optimize iron flocculation. We found viruses were best concentrated in seawater with five-fold higher iron concentrations than the standard used, inhibition of DNase activity reduced purification effectiveness, and five-fold more iron was needed to flocculate viruses from freshwater than seawater—critical knowledge for those seeking to apply this broadly used method to freshwater virome samples. Overall, our results demonstrated that ultrafiltration and purification performed better than iron chloride flocculation and purification in the tested matrices. Given that the method performance depended on the solids content and salinity of the samples, we suggest spike-and-recovery experiments be applied when concentrating and purifying sample types that diverge from those tested here.

Functional Restoration of Bacteriomes and Viromes by Fecal Microbiota Transplantation

BACKGROUND & AIMS

Fecal microbiota transplantation (FMT) is an effective therapy for recurrent Clostridioides difficile infection (rCDI). However, the overall mechanisms underlying FMT success await comprehensive elucidation, and the safety of FMT has recently become a serious concern because of the occurrence of drug-resistant bacteremia transmitted by FMT. We investigated whether functional restoration of the bacteriomes and viromes by FMT could be an indicator of successful FMT.

METHODS

The human intestinal bacteriomes and viromes from 9 patients with rCDI who had undergone successful FMT and their donors were analyzed. Prophage-based and CRISPR spacer-based host bacteria-phage associations in samples from recipients before and after FMT and in donor samples were examined. The gene functions of intestinal microorganisms affected by FMT were evaluated.

RESULTS

Metagenomic sequencing of both the viromes and bacteriomes revealed that FMT does change the characteristics of intestinal bacteriomes and viromes in recipients after FMT compared with those before FMT. In particular, many Proteobacteria, the fecal abundance of which was high before FMT, were eliminated, and the proportion of Microviridae increased in recipients. Most temperate phages also behaved in parallel with the host bacteria that were altered by FMT. Furthermore, the identification of bacterial and viral gene functions before and after FMT revealed that some distinctive pathways, including fluorobenzoate degradation and secondary bile acid biosynthesis, were significantly represented.

CONCLUSIONS

The coordinated action of phages and their host bacteria restored the recipients' intestinal flora. These findings show that the restoration of intestinal microflora functions reflects the success of FMT.

Interaction dynamics and virus-host range for estuarine actinophages captured by epicPCR

Viruses impact microbial diversity, gene flow and function through virus-host interactions. Although metagenomics surveys are rapidly cataloguing viral diversity, methods are needed to capture specific virus-host interactions in situ. Here, we leveraged metagenomics and repurposed emulsion paired isolation-concatenation PCR (epicPCR) to investigate viral diversity and virus-host interactions in situ over time in an estuarine environment. The method fuses a phage marker, the ribonucleotide reductase gene, with the host 16S rRNA gene of infected bacterial cells within emulsion droplets providing single-cell resolution for dozens of samples. EpicPCR captured in situ virus-host interactions for viral clades with no closely related database representatives. Abundant freshwater Actinobacteria lineages, in particular Rhodoluna sp., were the most common hosts for these poorly characterized viruses, with interactions correlated with environmental factors. Multiple methods used to identify virus-host interactions, including epicPCR, identified different and largely non-overlapping interactions within the vast virus-host interaction space. Tracking virus-host interaction dynamics also revealed that multi-host viruses had significantly longer periods with observed virus-host interactions, whereas single-host viruses were observed interacting with hosts at lower minimum abundances, suggesting more efficient interactions. Capturing in situ interactions with epicPCR revealed environmental and ecological factors shaping virus-host interactions, highlighting epicPCR as a valuable technique in viral ecology.

Cyanophage Diversity and Community Structure in Dead Zone Sediments

Cyanophages are viruses that target cyanobacteria and directly control their abundance via viral lysis. Cyanobacteria are known to cause large blooms in water bodies, substantially contributing to oxygen depletion in bottom waters resulting in areas called dead zones.

Up to 20% of prokaryotic organisms in the oceans are estimated to die every day due to viral infection and lysis. Viruses can therefore alter microbial diversity, community structure, and biogeochemical processes driven by these organisms. Cyanophages are viruses that infect and lyse cyanobacterial cells, adding bioavailable carbon and nutrients into the environment. Cyanobacteria are photosynthesizing bacteria, with some species capable of N₂ fixation, which are known to form large blooms as well as resistant resting cells known as akinetes. Here, we investigated cyanophage diversity and community structure plus cyanobacteria in dead zone sediments. We sampled surface sediments and sequenced DNA and RNA, along an oxygen gradient—representing oxic, hypoxic, and anoxic conditions—in one of the world’s largest dead zones located in the Baltic Sea. Cyanophages were detected at all stations and, based on partial genome contigs, had a higher alpha diversity and different beta diversity in the hypoxic-anoxic sediments, suggesting that cyanobacteria in dead zone sediments and/or environmental conditions select for specific cyanophages. Some of these cyanophages can infect cyanobacteria with potential consequences for gene expression related to their photosystem and phosphate regulation. Top cyanobacterial genera detected in the anoxic sediment included Dolichospermum/Anabaena, Synechococcus, and Cyanobium. RNA transcripts classified to cyanobacteria were associated with numerous pathways, including anaerobic carbon metabolism and N₂ fixation. Cyanobacterial blooms are known to fuel oxygen-depleted ecosystems with phosphorus (so-called internal loading), and our cyanophage data indicate the potential for viral lysis of cyanobacteria which might explain the high nutrient turnover in these environments.

IMPORTANCE Cyanophages are viruses that target cyanobacteria and directly control their abundance via viral lysis. Cyanobacteria are known to cause large blooms in water bodies, substantially contributing to oxygen depletion in bottom waters resulting in areas called dead zones. Our knowledge of cyanophages in dead zones is very scarce, and so far, no studies have assembled partial cyanophage genomes and investigated their associated cyanobacteria in these dark and anoxic sediments. Here, we present the first study using DNA and RNA sequencing to investigate in situ diversity of cyanophages and cyanobacteria in dead zones. Our study shows that dead zone sediments contain different cyanophages compared to oxic sediments and suggest that these viruses are able to affect cyanobacterial photosystem and phosphate regulation. Furthermore, cyanophage-controlled lysis of cyanobacteria might also increase the turnover of carbon, phosphorus, and nitrogen in these oxygen-free environments at the bottom of the sea.

Potential virus-mediated nitrogen cycling in oxygen-depleted oceanic waters

Viruses play an important role in the ecology and biogeochemistry of marine ecosystems. Beyond mortality and gene transfer, viruses can reprogram microbial metabolism during infection by expressing auxiliary metabolic genes (AMGs) involved in photosynthesis, central carbon metabolism, and nutrient cycling. While previous studies have focused on AMG diversity in the sunlit and dark ocean, less is known about the role of viruses in shaping metabolic networks along redox gradients associated with marine oxygen minimum zones (OMZs). Here, we analyzed relatively quantitative viral metagenomic datasets that profiled the oxygen gradient across Eastern Tropical South Pacific (ETSP) OMZ waters, assessing whether OMZ viruses might impact nitrogen (N) cycling via AMGs. Identified viral genomes encoded six N-cycle AMGs associated with denitrification, nitrification, assimilatory nitrate reduction, and nitrite transport. The majority of these AMGs (80%) were identified in T4-like Myoviridae phages, predicted to infect Cyanobacteria and Proteobacteria, or in unclassified archaeal viruses predicted to infect Thaumarchaeota. Four AMGs were exclusive to anoxic waters and had distributions that paralleled homologous microbial genes. Together, these findings suggest viruses modulate N-cycling processes within the ETSP OMZ and may contribute to nitrogen loss throughout the global oceans thus providing a baseline for their inclusion in the ecosystem and geochemical models.

VirION2: a short- and long-read sequencing and informatics workflow to study the genomic diversity of viruses in nature

Microbes play fundamental roles in shaping natural ecosystem properties and functions, but do so under constraints imposed by their viral predators. However, studying viruses in nature can be challenging due to low biomass and the lack of universal gene markers. Though metagenomic short-read sequencing has greatly improved our virus ecology toolkit—and revealed many critical ecosystem roles for viruses—microdiverse populations and fine-scale genomic traits are missed. Some of these microdiverse populations are abundant and the missed regions may be of interest for identifying selection pressures that underpin evolutionary constraints associated with hosts and environments. Though long-read sequencing promises complete virus genomes on single reads, it currently suffers from high DNA requirements and sequencing errors that limit accurate gene prediction. Here we introduce VirION2, an integrated short- and long-read metagenomic wet-lab and informatics pipeline that updates our previous method (VirION) to further enhance the utility of long-read viral metagenomics. Using a viral mock community, we first optimized laboratory protocols (polymerase choice, DNA shearing size, PCR cycling) to enable 76% longer reads (now median length of 6,965 bp) from 100-fold less input DNA (now 1 nanogram). Using a virome from a natural seawater sample, we compared viromes generated with VirION2 against other library preparation options (unamplified, original VirION, and short-read), and optimized downstream informatics for improved long-read error correction and assembly. VirION2 assemblies combined with short-read based data (‘enhanced’ viromes), provided significant improvements over VirION libraries in the recovery of longer and more complete viral genomes, and our optimized error-correction strategy using long- and short-read data achieved 99.97% accuracy. In the seawater virome, VirION2 assemblies captured 5,161 viral populations (including all of the virus populations observed in the other assemblies), 30% of which were uniquely assembled through inclusion of long-reads, and 22% of the top 10% most abundant virus populations derived from assembly of long-reads. Viral populations unique to VirION2 assemblies had significantly higher microdiversity means, which may explain why short-read virome approaches failed to capture them. These findings suggest the VirION2 sample prep and workflow can help researchers better investigate the virosphere, even from challenging low-biomass samples. Our new protocols are available to the research community on protocols.io as a ‘living document’ to facilitate dissemination of updates to keep pace with the rapid evolution of long-read sequencing technology.

Genomic Characterization and Distribution Pattern of a Novel Marine OM43 Phage

Bacteriophages have a significant impact on the structure and function of marine microbial communities. Phages of some major bacterial lineages have recently been shown to dominate the marine viral communities. However, phages that infect many important bacterial clades still remained unexplored. Members of the marine OM43 clade are methylotrophs that play important roles in C1 metabolism. OM43 phages (phages that infect the OM43 bacteria) represent an understudied viral group with only one known isolate. In this study, we describe the genomic characterization and biogeography of an OM43 phage that infects the strain HTCC2181, designated MEP301. MEP301 has a genome size of 34,774 bp. We found that MEP301 is genetically distinct from other known phage isolates and only displays significant sequence similarity with some metagenomic viral genomes (MVGs). A total of 12 MEP301-type MVGs were identified from metagenomic datasets. Comparative genomic and phylogenetic analyses revealed that MEP301-type phages can be separated into two subgroups (subgroup I and subgroup II). We also performed a metagenomic recruitment analysis to determine the relative abundance of reads mapped to these MEP301-type phages, which suggested that subgroup I MEP301-type phages are present predominantly in the cold upper waters with lower salinity. Notably, subgroup II phages have an inverse different distribution pattern, implying that they may infect hosts from a distinct OM43 subcluster. Our study has expanded the knowledge about the genomic diversity of marine OM43 phages and identified a new phage group that is widespread in the ocean.

Hybrid assembly of an agricultural slurry virome reveals a diverse and stable community with the potential to alter the metabolism and virulence of veterinary pathogens

BACKGROUND

Viruses are the most abundant biological entities on Earth, known to be crucial components of microbial ecosystems. However, there is little information on the viral community within agricultural waste. There are currently ~ 2.7 million dairy cattle in the UK producing 7-8% of their own bodyweight in manure daily, and 28 million tonnes annually. To avoid pollution of UK freshwaters, manure must be stored and spread in accordance with guidelines set by DEFRA. Manures are used as fertiliser, and widely spread over crop fields, yet little is known about their microbial composition. We analysed the virome of agricultural slurry over a 5-month period using short and long-read sequencing.

RESULTS

Hybrid sequencing uncovered more high-quality viral genomes than long or short-reads alone; yielding 7682 vOTUs, 174 of which were complete viral genomes. The slurry virome was highly diverse and dominated by lytic bacteriophage, the majority of which represent novel genera (~ 98%). Despite constant influx and efflux of slurry, the composition and diversity of the slurry virome was extremely stable over time, with 55% of vOTUs detected in all samples over a 5-month period. Functional annotation revealed a diverse and abundant range of auxiliary metabolic genes and novel features present in the community, including the agriculturally relevant virulence factor VapE, which was widely distributed across different phage genera that were predicted to infect several hosts. Furthermore, we identified an abundance of phage-encoded diversity-generating retroelements, which were previously thought to be rare on lytic viral genomes. Additionally, we identified a group of crAssphages, including lineages that were previously thought only to be found in the human gut.

CONCLUSIONS

The cattle slurry virome is complex, diverse and dominated by novel genera, many of which are not recovered using long or short-reads alone. Phages were found to encode a wide range of AMGs that are not constrained to particular groups or predicted hosts, including virulence determinants and putative ARGs. The application of agricultural slurry to land may therefore be a driver of bacterial virulence and antimicrobial resistance in the environment. Video abstract.

Uncultivated Viral Populations Dominate Estuarine Viromes on the Spatiotemporal Scale

Viruses are ubiquitous and abundant in the oceans, and viral metagenomes (viromes) have been investigated extensively via several large-scale ocean sequencing projects. However, there have not been any systematic viromic studies in estuaries. Here, we investigated the viromes of the Delaware Bay and Chesapeake Bay, two Mid-Atlantic estuaries. Deep sequencing generated a total of 48,190 assembled viral sequences (>5 kb) and 26,487 viral populations (9,204 virus clusters and 17,845 singletons), including 319 circular viral contigs between 7.5 kb and 161.8 kb. Unknown viruses represented the vast majority of the dominant populations, while the composition of known viruses, such as pelagiphage and cyanophage, appeared to be relatively consistent across a wide range of salinity gradients and in different seasons. A difference between estuarine and ocean viromes was reflected by the proportions of Myoviridae, Podoviridae, Siphoviridae, Phycodnaviridae, and a few well-studied virus representatives. The difference in viral community between the Delaware Bay and Chesapeake Bay is significantly more pronounced than the difference caused by temperature or salinity, indicating strong local profiles caused by the unique ecology of each estuary. Interestingly, a viral contig similar to phages infecting Acinetobacter baumannii (“Iraqibacter”) was found to be highly abundant in the Delaware Bay but not in the Chesapeake Bay, the source of which is yet to be identified. Highly abundant viruses in both estuaries have close hits to viral sequences derived from the marine single-cell genomes or long-read single-molecule sequencing, suggesting that important viruses are still waiting to be discovered in the estuarine environment.

IMPORTANCE This is the first systematic study about spatial and temporal variation of virioplankton communities in estuaries using deep metagenomics sequencing. It is among the highest-quality viromic data sets to date, showing remarkably consistent sequencing depth and quality across samples. Our results indicate that there exists a large pool of abundant and diverse viruses in estuaries that have not yet been cultivated, their genomes only available thanks to single-cell genomics or single-molecule sequencing, demonstrating the importance of these methods for viral discovery. The spatiotemporal pattern of these abundant uncultivated viruses is more variable than that of cultured viruses. Despite strong environmental gradients, season and location had surprisingly little impact on the viral community within an estuary, but we saw a significant distinction between the two estuaries and also between estuarine and open ocean viromes.

The communities and functional profiles of virioplankton along a salinity gradient in a subtropical estuary

Viruses are the major drivers shaping microorganismal communities, and impact marine biogeochemical cycling. They are affected by various environmental parameters, such as salinity. Although the spatiotemporal distribution and dynamics of virioplankton have been extensively studied in saline environments, few detailed studies of community structure and function of viruses along salinity gradients have been conducted. Here, we used the 16S and 18S rRNA gene amplicon and metagenomic sequencing from a subtropical estuary (Pearl River Estuary, PRE; located in Shenzhen, Guangdong Province, China) to explore how viral community composition and function vary along a salinity gradient. Results showed that the detected viruses were mainly bacteriophages. The double-stranded DNA viruses were the most abundant (especially Siphoviridae, Myoviridae, Mimiviridae, Phycodnaviridae, and Podoviridae), followed by a small number of single-stranded DNA (Circoviridae) and RNA (Retroviridae) viruses. Viral biodiversity significantly declined and community structure varied greatly along the salinity gradient. The salinity, ammonium and dissolved oxygen were dominated factors influencing the community composition of viruses. Association network analysis showed that viruses had a negative effect on multiple host taxa (prokaryotic and eukaryotic species). Metagenomic data revealed that the main viral functional potential was involved in organic matter metabolism by carbohydrate-active enzymes (CAZymes). Deeper comparative functional analyses showed that viruses in the low-salinity environment had more carbohydrate-binding module and glycosidase hydrolases activities than those under high-salinity conditions. However, an opposite pattern was observed for carbohydrate esterases. These results suggest that virus-encoded CAZyme genes may alter the bacterial metabolism in estuaries. Overall, our results demonstrate that there is a spatial heterogeneity in the composition and function of virioplankton along a salinity gradient. This study enhances our understanding of viral distribution and their contribution to regulating carbon degradation throughout environments with varying salinities in subtropical estuaries.

The Fennoscandian Shield deep terrestrial virosphere suggests slow motion 'boom and burst' cycles

The deep biosphere contains members from all three domains of life along with viruses. Here we investigate the deep terrestrial virosphere by sequencing community nucleic acids from three groundwaters of contrasting chemistries, origins, and ages. These viromes constitute a highly unique community compared to other environmental viromes and sequenced viral isolates. Viral host prediction suggests that many of the viruses are associated with Firmicutes and Patescibacteria, a superphylum lacking previously described active viruses. RNA transcript-based activity implies viral predation in the shallower marine water-fed groundwater, while the deeper and more oligotrophic waters appear to be in ‘metabolic standby’. Viral encoded antibiotic production and resistance systems suggest competition and antagonistic interactions. The data demonstrate a viral community with a wide range of predicted hosts that mediates nutrient recycling to support a higher microbial turnover than previously anticipated. This suggests the presence of ‘kill-the-winner’ oscillations creating slow motion ‘boom and burst’ cycles.

Karin Holmfeldt et al. sequence metagenomes and metatranscriptomes of viruses in deep groundwaters down to 448 m below the surface. The results reveal ecological dynamics of viruses including slow motion ‘boom and burst’ cycles and a ‘kill the winner’ model potentially driven by viral predation.

Metaviromics coupled with phage-host identification to open the viral 'black box'

Viruses are found in almost all biomes on Earth, with bacteriophages (phages) accounting for the majority of viral particles in most ecosystems. Phages have been isolated from natural environments using the plaque assay and liquid medium-based dilution culturing. However, phage cultivation is restricted by the current limitations in the number of culturable bacterial strains. Unlike prokaryotes, which possess universally conserved 16S rRNA genes, phages lack universal marker genes for viral taxonomy, thus restricting cultureindependent analyses of viral diversity. To circumvent these limitations, shotgun viral metagenome sequencing (i.e., metaviromics) has been developed to enable the extensive sequencing of a variety of viral particles present in the environment and is now widely used. Using metaviromics, numerous studies on viral communities have been conducted in oceans, lakes, rivers, and soils, resulting in many novel phage sequences. Furthermore, auxiliary metabolic genes such as ammonic monooxygenase C and β-lactamase have been discovered in viral contigs assembled from viral metagenomes. Current attempts to identify putative bacterial hosts of viral metagenome sequences based on sequence homology have been limited due to viral sequence variations. Therefore, culture-independent approaches have been developed to predict bacterial hosts using single-cell genomics and fluorescentlabeling. This review focuses on recent viral metagenome studies conducted in natural environments, especially in aquatic ecosystems, and their contributions to phage ecology. Here, we concluded that although metaviromics is a key tool for the study of viral ecology, this approach must be supplemented with phage-host identification, which in turn requires the cultivation of phage-bacteria systems.

Shotgun metagenomic sequencing reveals the full taxonomic, trophic, and functional diversity of a coral reef benthic cyanobacterial mat from Bonaire, Caribbean Netherlands

Anthropogenic forcing is spurring cyanobacterial proliferation in aquatic ecosystems worldwide. While planktonic cyanobacterial blooms have received substantial research attention, benthic blooms of mat-forming cyanobacteria have received considerably less attention, especially benthic mat blooms on coral reefs. Resultingly, numerous aspects of coral reef benthic cyanobacterial bloom ecology remain unknown, including underlying biodiversity in the mat communities. Most previous characterizations of coral reef cyanobacterial mat composition have only considered the cyanobacterial component. Without an unbiased characterization of full community diversity, we cannot predict whole-community response to anthropogenic inputs or effectively determine appropriate mitigation strategies. Here, we advocate for the implementation of shotgun sequencing techniques to study coral reef cyanobacterial mats worldwide, utilizing a case study of a coral reef benthic cyanobacterial mat sampled from the island of Bonaire, Caribbean Netherlands. Read-based taxonomic profiling revealed that Cyanobacteria was present at only 47.57% relative abundance in a coral reef cyanobacterial mat, with non-cyanobacterial members of the sampled mat community, including diatoms (0.78%), fungi (0.25%), Archaea (0.34%), viruses (0.08%), and other bacteria (45.78%), co-dominating the community. We found numerous gene families for regulatory systems and for functional pathways (both aerobic and anaerobic). These gene families were involved in community coordination; photosynthesis; nutrient scavenging; and the cycling of sulfur, nitrogen, phosphorous, and iron. We also report bacteriophage (including prophage) sequences associated with this subtidal coral reef cyanobacterial mat, which could contribute to intra-mat nutrient cycling and bloom dynamics. Overall, our results suggest that Cyanobacteria-focused analysis of coral reef cyanobacterial mats underestimates mat diversity and fails to capture community members possessing broad metabolic potential for intra-mat nutrient scavenging, recycling, and retention that likely contribute to the contemporary success of cyanobacterial mats on reefs. We advocate for increased collaboration between microbiologists and coral reef ecologists to unite insights from each discipline and improve efforts to understand mat ecology.

VirSorter2: a multi-classifier, expert-guided approach to detect diverse DNA and RNA viruses

BACKGROUND

Viruses are a significant player in many biosphere and human ecosystems, but most signals remain "hidden" in metagenomic/metatranscriptomic sequence datasets due to the lack of universal gene markers, database representatives, and insufficiently advanced identification tools.

RESULTS

Here, we introduce VirSorter2, a DNA and RNA virus identification tool that leverages genome-informed database advances across a collection of customized automatic classifiers to improve the accuracy and range of virus sequence detection. When benchmarked against genomes from both isolated and uncultivated viruses, VirSorter2 uniquely performed consistently with high accuracy (F1-score > 0.8) across viral diversity, while all other tools under-detected viruses outside of the group most represented in reference databases (i.e., those in the order Caudovirales). Among the tools evaluated, VirSorter2 was also uniquely able to minimize errors associated with atypical cellular sequences including eukaryotic genomes and plasmids. Finally, as the virosphere exploration unravels novel viral sequences, VirSorter2's modular design makes it inherently able to expand to new types of viruses via the design of new classifiers to maintain maximal sensitivity and specificity.

CONCLUSION

With multi-classifier and modular design, VirSorter2 demonstrates higher overall accuracy across major viral groups and will advance our knowledge of virus evolution, diversity, and virus-microbe interaction in various ecosystems. Source code of VirSorter2 is freely available ( https://bitbucket.org/MAVERICLab/virsorter2 ), and VirSorter2 is also available both on bioconda and as an iVirus app on CyVerse ( https://de.cyverse.org/de ). Video abstract.

Eukaryotic virus composition can predict the efficiency of carbon export in the global ocean

The biological carbon pump, in which carbon fixed by photosynthesis is exported to the deep ocean through sinking, is a major process in Earth's carbon cycle. The proportion of primary production that is exported is termed the carbon export efficiency (CEE). Based on in-lab or regional scale observations, viruses were previously suggested to affect the CEE (i.e., viral “shunt” and “shuttle”). In this study, we tested associations between viral community composition and CEE measured at a global scale. A regression model based on relative abundance of viral marker genes explained 67% of the variation in CEE. Viruses with high importance in the model were predicted to infect ecologically important hosts. These results are consistent with the view that the viral shunt and shuttle functions at a large scale and further imply that viruses likely act in this process in a way dependent on their hosts and ecosystem dynamics.

•

Eukaryotic virus community composition is shown to predict carbon export efficiency

•

Tens of viruses are highly important in the prediction of the efficiency

•

These viruses are inferred to infect ecologically important hosts

Classifying the Lifestyle of Metagenomically-Derived Phages Sequences Using Alignment-Free Methods

Phages are viruses that infect bacteria. The phages can be classified into two different categories based on their lifestyles: temperate and lytic. Now, the metavirome can generate a large number of fragments from the viral genomic sequences of entire environmental community, which makes it impossible to determine their lifestyles through experiments. Thus, there is a need to development computational methods for annotating phage contigs and making prediction of their lifestyles. Alignment-based methods for classifying phage lifestyle are limited by incomplete assembled genomes and nucleotide databases. Alignment-free methods based on the frequencies of k-mers were widely used for genome and metagenome comparison which did not rely on the completeness of genome or nucleotide databases. To mimic fragmented metagenomic sequences, the temperate and lytic phages genomic sequences were split into non-overlapping fragments with different lengths, then, I comprehensively compared nine alignment-free dissimilarity measures with a wide range of choices of k-mer length and Markov orders for predicting the lifestyles of these phage contigs. The dissimilarity measure, d2S, performed better than other dissimilarity measures for classifying the lifestyles of phages. Thus, I propose that the alignment-free method, d2S, can be used for predicting the lifestyles of phages which derived from the metagenomic data.

Gut Microbiota in Dholes During Estrus

The co-evolution of gut microbes and the host plays a vital role in the survival and reproduction of the host. The dhole (Cuon alpinus) has been listed as endangered species by the International Union for Conservation of Nature; therefore, conservation and effective breeding of dholes are essential. Effective estrus can promote reproduction. However, little is known about the relative contribution of estrus in shaping the structure and the functions of fecal microbiota. Here, we investigated the potential association between estrus and the fecal microbiota in dholes using shotgun metagenomic sequencing. We found that the estrus stages in dholes vary significantly in terms of gut bacterial composition and microbiome metabolism and function. Compared with that of non-estrus, adult dholes, the microbiome of estrus adult dholes had a significantly higher abundance of Bacillus faecalis and Veillonella, which play a key role in the synthesis of sex hormones and nucleic acids, energy production, and reproductive cell division. The insulin and energy metabolism-related pathways are significantly enhanced in the gut microbes and the related gluconeogenic enzymes are significantly enriched during estrus. These findings suggest that the structure and metagenome of the fecal microbiome during estrus have a significant effect in promoting estrus in dholes, thus providing a new perspective for dhole conservation.

Historical Discoveries on Viruses in the Environment and Their Impact on Public Health

BACKGROUND

Transmission of many viruses occurs by direct transmission during a close contact between two hosts, or by an indirect transmission through the environment. Several and often interconnected factors, both abiotic and biotic, determine the persistence of these viruses released in the environment, which can last from a few seconds to several years. Moreover, viruses in the environment are able to travel short to very long distances, especially in the air or in water.

SUMMARY

Although well described now, the role of these environments as intermediaries or as reservoirs in virus transmission has been extensively studied and debated in the last century. The majority of these discoveries, such as the pioneer work on bacteria transmission, the progressive discoveries of viruses, as well as the persistence of the influenza virus in the air varying along with droplet sizes, or the role of water in the transmission of poliovirus, have contributed to the improvement of public health. Recent outbreaks of human coronavirus, influenza virus, and Ebola virus have also demonstrated the contemporaneity of these research studies and the need to study virus persistence in the environment. Key Messages: In this review, we discuss historical discoveries that contributed to describe biotic and abiotic factors determining viral persistence in the environment.

Genome-resolved viral ecology in a marine oxygen minimum zone

Oxygen minimum zones (OMZs) are critical to marine nitrogen cycling and global climate change. While OMZ microbial communities are relatively well-studied, little is known about their viruses. Here, we assess the viral community ecology of 22 deeply sequenced viral metagenomes along a gradient of oxygenated to anoxic waters (<0.02 μmol/l O₂ ) in the Eastern Tropical South Pacific (ETSP) OMZ. We identified 46 127 viral populations (≥5 kb), which augments the known viruses from ETSP by 10-fold. Viral communities clustered into six groups that correspond to oceanographic features. Oxygen concentration was the predominant environmental feature driving viral community structure. Alpha and beta diversity of viral communities in the anoxic zone were lower than in surface waters, which parallels the low microbial diversity seen in other studies. ETSP viruses were largely endemic, with the majority of shared viruses (87%) also present in other OMZ samples. We detected 543 putative viral-encoded auxiliary metabolic genes (AMGs), of which some have a distribution that reflects physico-chemical characteristics across depth. Together these findings provide an ecological baseline for viral community structure, drivers and population variability in OMZs that will help future studies assess the role of viruses in these climate-critical environments.

Analysis of the vaginal microbiome of giant pandas using metagenomics sequencing

In this study, a total of 14 vaginal samples (GPV1‐14) from giant pandas were analyzed. These vaginal samples were divided into two groups as per the region and age of giant pandas. All the vaginal samples were analyzed using metagenomic sequencing. As per the outcomes of metagenomic analysis, Proteobacteria (39.04%), Firmicutes (5.27%), Actinobacteria (2.94%), and Basidiomycota (2.77%) were found to be the dominant phyla in the microbiome of the vaginal samples. At the genus level, Pseudomonas (21.90%) was found to be the most dominant genus, followed by Streptococcus (3.47%), Psychrobacter (1.89%), and Proteus (1.38%). Metastats analysis of the microbial species in the vaginal samples of giant pandas from Wolong Nature Reserve, Dujiangyan and Ningbo Youngor Zoo, and Ya'an Bifengxia Nature Reserve was found to be significantly different (p < 0.05). Age groups, that is, AGE1 (5‐10 years old) and AGE2 (11‐16 years old), also demonstrated significantly different inter‐group microbial species (p < 0.05). For the first time, Chlamydia and Neisseria gonorrhoeae were detected in giant pandas’ reproductive tract. GPV3 vaginal sample (2.63%) showed highest Chlamydia content followed by GPV14 (0.91%), and GPV7 (0.62%). GPV5 vaginal sample (7.17%) showed the highest Neisseria gonorrhoeae content, followed by GPV14 (7.02%), and GPV8 (6.50%). Furthermore, we employed eggNOG, CAZy, KEGG, and NCBI databases to investigate the functional significance of giant panda's vaginal microbial community. The outcomes indicated that giant panda's vaginal microbes were involved in biological processes. The data from this study will help in improving the reproductive health of giant pandas.

New viral biogeochemical roles revealed through metagenomic analysis of Lake Baikal

BACKGROUND

Lake Baikal is the largest body of liquid freshwater on Earth. Previous studies have described the microbial composition of this habitat, but the viral communities from this ecosystem have not been characterized in detail.

RESULTS

Here, we describe the viral diversity of this habitat across depth and seasonal gradients. We discovered 19,475 bona fide viral sequences, which are derived from viruses predicted to infect abundant and ecologically important taxa that reside in Lake Baikal, such as Nitrospirota, Methylophilaceae, and Crenarchaeota. Diversity analysis revealed significant changes in viral community composition between epipelagic and bathypelagic zones. Analysis of the gene content of individual viral populations allowed us to describe one of the first bacteriophages that infect Nitrospirota, and their extensive repertoire of auxiliary metabolic genes that might enhance carbon fixation through the reductive TCA cycle. We also described bacteriophages of methylotrophic bacteria with the potential to enhance methanol oxidation and the S-adenosyl-L-methionine cycle.

CONCLUSIONS

These findings unraveled new ways by which viruses influence the carbon cycle in freshwater ecosystems, namely, by using auxiliary metabolic genes that act upon metabolisms of dark carbon fixation and methylotrophy. Therefore, our results shed light on the processes through which viruses can impact biogeochemical cycles of major ecological relevance. Video Abstract.

The Gut Virome Database Reveals Age-Dependent Patterns of Virome Diversity in the Human Gut

The gut microbiome profoundly affects human health and disease, and their infecting viruses are likely as important, but often missed because of reference database limitations. Here, we (1) built a human Gut Virome Database (GVD) from 2,697 viral particle or microbial metagenomes from 1,986 individuals representing 16 countries, (2) assess its effectiveness, and (3) report a meta-analysis that reveals age-dependent patterns across healthy Westerners. The GVD contains 33,242 unique viral populations (approximately species-level taxa) and improves average viral detection rates over viral RefSeq and IMG/VR nearly 182-fold and 2.6-fold, respectively. GVD meta-analyses show highly personalized viromes, reveal that inter-study variability from technical artifacts is larger than any "disease" effect at the population level, and document how viral diversity changes from human infancy into senescence. Together, this compact foundational resource, these standardization guidelines, and these meta-analysis findings provide a systematic toolkit to help maximize our understanding of viral roles in health and disease.

Seasonal and diel patterns of abundance and activity of viruses in the Red Sea

Virus-microbe interactions have been studied in great molecular details for many years in cultured model systems, yielding a plethora of knowledge on how viruses use and manipulate host machinery. Since the advent of molecular techniques and high-throughput sequencing, methods such as cooccurrence, nucleotide composition, and other statistical frameworks have been widely used to infer virus-microbe interactions, overcoming the limitations of culturing methods. However, their accuracy and relevance is still debatable as cooccurrence does not necessarily mean interaction. Here we introduce an ecological perspective of marine viral communities and potential interaction with their hosts, using analyses that make no prior assumptions on specific virus-host pairs. By size fractionating water samples into free viruses and microbes (i.e., also viruses inside or attached to their hosts) and looking at how viral group abundance changes over time along both fractions, we show that the viral community is undergoing a change in rank abundance across seasons, suggesting a seasonal succession of viruses in the Red Sea. We use abundance patterns in the different size fractions to classify viral clusters, indicating potential diverse interactions with their hosts and potential differences in life history traits between major viral groups. Finally, we show hourly resolved variations of intracellular abundance of similar viral groups, which might indicate differences in their infection cycles or metabolic capacities.

VIRIDIC-A Novel Tool to Calculate the Intergenomic Similarities of Prokaryote-Infecting Viruses

Nucleotide-based intergenomic similarities are useful to understand how viruses are related with each other and to classify them. Here we have developed VIRIDIC, which implements the traditional algorithm used by the International Committee on Taxonomy of Viruses (ICTV), Bacterial and Archaeal Viruses Subcommittee, to calculate virus intergenomic similarities. When compared with other software, VIRIDIC gave the best agreement with the traditional algorithm, which is based on the percent identity between two genomes determined by BLASTN. Furthermore, VIRIDIC proved best at estimating the relatedness between more distantly-related phages, relatedness that other tools can significantly overestimate. In addition to the intergenomic similarities, VIRIDIC also calculates three indicators of the alignment ability to capture the relatedness between viruses: the aligned fractions for each genome in a pair and the length ratio between the two genomes. The main output of VIRIDIC is a heatmap integrating the intergenomic similarity values with information regarding the genome lengths and the aligned genome fraction. Additionally, VIRIDIC can group viruses into clusters, based on user-defined intergenomic similarity thresholds. The sensitivity of VIRIDIC is given by the BLASTN. Thus, it is able to capture relationships between viruses having in common even short genomic regions, with as low as 65% similarity. Below this similarity level, protein-based analyses should be used, as they are the best suited to capture distant relationships. VIRIDIC is available at viridic.icbm.de, both as a web-service and a stand-alone tool. It allows fast analysis of large phage genome datasets, especially in the stand-alone version, which can be run on the user's own servers and can be integrated in bioinformatics pipelines. VIRIDIC was developed having viruses of Bacteria and Archaea in mind; however, it could potentially be used for eukaryotic viruses as well, as long as they are monopartite.

Exploration into the origins and mobilization of di-hydrofolate reductase genes and the emergence of clinical resistance to trimethoprim

Trimethoprim is a synthetic antibacterial agent that targets folate biosynthesis by competitively binding to the di-hydrofolate reductase enzyme (DHFR). Trimethoprim is often administered synergistically with sulfonamide, another chemotherapeutic agent targeting the di-hydropteroate synthase (DHPS) enzyme in the same pathway. Clinical resistance to both drugs is widespread and mediated by enzyme variants capable of performing their biological function without binding to these drugs. These mutant enzymes were assumed to have arisen after the discovery of these synthetic drugs, but recent work has shown that genes conferring resistance to sulfonamide were present in the bacterial pangenome millions of years ago. Here, we apply phylogenetics and comparative genomics methods to study the largest family of mobile trimethoprim-resistance genes (dfrA). We show that most of the dfrA genes identified to date map to two large clades that likely arose from independent mobilization events. In contrast to sulfonamide resistance (sul) genes, we find evidence of recurrent mobilization in dfrA genes. Phylogenetic evidence allows us to identify novel dfrA genes in the emerging pathogen Acinetobacter baumannii, and we confirm their resistance phenotype in vitro. We also identify a cluster of dfrA homologues in cryptic plasmid and phage genomes, but we show that these enzymes do not confer resistance to trimethoprim. Our methods also allow us to pinpoint the chromosomal origin of previously reported dfrA genes, and we show that many of these ancient chromosomal genes also confer resistance to trimethoprim. Our work reveals that trimethoprim resistance predated the clinical use of this chemotherapeutic agent, but that novel mutations have likely also arisen and become mobilized following its widespread use within and outside the clinic. Hence, this work confirms that resistance to novel drugs may already be present in the bacterial pangenome, and stresses the importance of rapid mobilization as a fundamental element in the emergence and global spread of resistance determinants.

Metatranscriptomic Identification of Diverse and Divergent RNA Viruses in Green and Chlorarachniophyte Algae Cultures

Our knowledge of the diversity and evolution of the virosphere will likely increase dramatically with the study of microbial eukaryotes, including the microalgae within which few RNA viruses have been documented. By combining total RNA sequencing with sequence and structural-based homology detection, we identified 18 novel RNA viruses in cultured samples from two major groups of microbial algae: the chlorophytes and the chlorarachniophytes. Most of the RNA viruses identified in the green algae class Ulvophyceae were related to the Tombusviridae and Amalgaviridae viral families commonly associated with land plants. This suggests that the evolutionary history of these viruses extends to divergence events between algae and land plants. Seven Ostreobium sp-associated viruses exhibited sequence similarity to the mitoviruses most commonly found in fungi, compatible with horizontal virus transfer between algae and fungi. We also document, for the first time, RNA viruses associated with chlorarachniophytes, including the first negative-sense (bunya-like) RNA virus in microalgae, as well as a distant homolog of the plant virus Virgaviridae, potentially signifying viral inheritance from the secondary chloroplast endosymbiosis that marked the origin of the chlorarachniophytes. More broadly, these data suggest that the scarcity of RNA viruses in algae results from limited investigation rather than their absence.

Culturing novel and abundant pelagiphages in the ocean

Viruses play a key role in biogeochemical cycling and host mortality, metabolism, physiology and evolution in the ocean. Viruses that infect the globally abundant SAR11 bacteria (pelagiphages) were reported to be an important component of the marine viral communities. Our current knowledge of pelagiphages is based on a few studies and therefore is limited. In this study, 10 new pelagiphages were isolated and genomically characterized. These pelagiphages represent the first cultivated representatives of four viral lineages only found in metagenomic sequencing datasets previously. Many abundant environmental viral sequences, i.e., single-virus vSAG 37-F6 and several Global Ocean Viromes (GOV) viral populations, are now further confirmed with these pelagiphages. Viromic read mapping reveals that these new pelagiphages are globally distributed in the ocean and can be detected throughout the water column. Remarkably, isolation of these pelagiphages contributed up to 12% of all viromic reads annotated in the analysed viromes. Altogether, this study has greatly broadened our understanding of pelagiphages regarding their morphology, genetic diversity, infection strategies, and distribution pattern. The availability of these newly isolated pelagiphages and their genome sequences will allow us to further explore their infectivities and ecological strategies.

Metatranscriptomic Identification of Diverse and Divergent RNA Viruses in Green and Chlorarachniophyte Algae Cultures

Our knowledge of the diversity and evolution of the virosphere will likely increase dramatically with the study of microbial eukaryotes, including the microalgae within which few RNA viruses have been documented. By combining total RNA sequencing with sequence and structural-based homology detection, we identified 18 novel RNA viruses in cultured samples from two major groups of microbial algae: the chlorophytes and the chlorarachniophytes. Most of the RNA viruses identified in the green algae class Ulvophyceae were related to the Tombusviridae and Amalgaviridae viral families commonly associated with land plants. This suggests that the evolutionary history of these viruses extends to divergence events between algae and land plants. Seven Ostreobium sp-associated viruses exhibited sequence similarity to the mitoviruses most commonly found in fungi, compatible with horizontal virus transfer between algae and fungi. We also document, for the first time, RNA viruses associated with chlorarachniophytes, including the first negative-sense (bunya-like) RNA virus in microalgae, as well as a distant homolog of the plant virus Virgaviridae, potentially signifying viral inheritance from the secondary chloroplast endosymbiosis that marked the origin of the chlorarachniophytes. More broadly, these data suggest that the scarcity of RNA viruses in algae results from limited investigation rather than their absence.

Single-virus genomics and beyond

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

Viral ecogenomics across the Porifera

BACKGROUND

Viruses directly affect the most important biological processes in the ocean via their regulation of prokaryotic and eukaryotic populations. Marine sponges form stable symbiotic partnerships with a wide diversity of microorganisms and this high symbiont complexity makes them an ideal model for studying viral ecology. Here, we used morphological and molecular approaches to illuminate the diversity and function of viruses inhabiting nine sponge species from the Great Barrier Reef and seven from the Red Sea.

RESULTS

Viromic sequencing revealed host-specific and site-specific patterns in the viral assemblages, with all sponge species dominated by the bacteriophage order Caudovirales but also containing variable representation from the nucleocytoplasmic large DNA virus families Mimiviridae, Marseilleviridae, Phycodnaviridae, Ascoviridae, Iridoviridae, Asfarviridae and Poxviridae. Whilst core viral functions related to replication, infection and structure were largely consistent across the sponge viromes, functional profiles varied significantly between species and sites largely due to differential representation of putative auxiliary metabolic genes (AMGs) and accessory genes, including those associated with herbicide resistance, heavy metal resistance and nylon degradation. Furthermore, putative AMGs varied with the composition and abundance of the sponge-associated microbiome. For instance, genes associated with antimicrobial activity were enriched in low microbial abundance sponges, genes associated with nitrogen metabolism were enriched in high microbial abundance sponges and genes related to cellulose biosynthesis were enriched in species that host photosynthetic symbionts.

CONCLUSIONS

Our results highlight the diverse functional roles that viruses can play in marine sponges and are consistent with our current understanding of sponge ecology. Differential representation of putative viral AMGs and accessory genes across sponge species illustrate the diverse suite of beneficial roles viruses can play in the functional ecology of these complex reef holobionts. Video Abstract.

Lysogeny in the oceans: Lessons from cultivated model systems and a reanalysis of its prevalence

In the oceans, viruses that infect bacteria (phages) influence a variety of microbially mediated processes that drive global biogeochemical cycles. The nature of their influence is dependent upon infection mode, be it lytic or lysogenic. Temperate phages are predicted to be prevalent in marine systems where they are expected to execute both types of infection modes. Understanding the range and outcomes of temperate phage-host interactions is fundamental for evaluating their ecological impact. Here, we (i) review phage-mediated rewiring of host metabolism, with a focus on marine systems, (ii) consider the range and nature of temperate phage-host interactions, and (iii) draw on studies of cultivated model systems to examine the consequences of lysogeny among several dominant marine bacterial lineages. We also readdress the prevalence of lysogeny among marine bacteria by probing a collection of 1239 publicly available bacterial genomes, representing cultured and uncultivated strains, for evidence of complete prophages. Our conservative analysis, anticipated to underestimate true prevalence, predicts 18% of the genomes examined contain at least one prophage, the majority (97%) were found within genomes of cultured isolates. These results highlight the need for cultivation of additional model systems to better capture the diversity of temperate phage-host interactions in the oceans.

Partitiviruses Infecting Drosophila melanogaster and Aedes aegypti Exhibit Efficient Biparental Vertical Transmission

Partitiviruses are segmented, multipartite double-stranded RNA (dsRNA) viruses that until recently were only known to infect fungi, plants, and protozoans. Metagenomic surveys have revealed that partitivirus-like sequences are also commonly associated with arthropods. One arthropod-associated partitivirus, galbut virus, is common in wild populations of Drosophila melanogaster To begin to understand the processes that underlie this virus's high global prevalence, we established colonies of wild-caught infected flies. Infection remained at stably high levels over 3 years, with between 63 and 100% of individual flies infected. Galbut virus infects fly cells and replicates in tissues throughout infected adults, including reproductive tissues and the gut epithelium. We detected no evidence of horizontal transmission via ingestion, but vertical transmission from either infected females or infected males was ∼100% efficient. Vertical transmission of a related partitivirus, verdadero virus, that we discovered in a laboratory colony of Aedes aegypti mosquitoes was similarly efficient. This suggests that efficient biparental vertical transmission may be a feature of at least a subset of insect-infecting partitiviruses. To study the impact of galbut virus infection free from the confounding effect of other viruses, we generated an inbred line of flies with galbut virus as the only detectable virus infection. We were able to transmit infection experimentally via microinjection of homogenate from these galbut-only flies. This sets the stage for experiments to understand the biological impact and possible utility of partitiviruses infecting model organisms and disease vectors.IMPORTANCE Galbut virus is a recently discovered partitivirus that is extraordinarily common in wild populations of the model organism Drosophila melanogaster Like for most viruses discovered through metagenomics, most of the basic biological questions about this virus remain unanswered. We found that galbut virus, along with a closely related partitivirus found in Aedes aegypti mosquitoes, is transmitted from infected females or males to offspring with ∼100% efficiency and can be maintained in laboratory colonies over years. This efficient transmission mechanism likely underlies the successful spread of these viruses through insect populations. We created Drosophila lines that contained galbut virus as the only virus infection and showed that these flies can be used as a source for experimental infections. This provides insight into how arthropod-infecting partitiviruses may be maintained in nature and sets the stage for exploration of their biology and potential utility.

Harnessing the strategy of metagenomics for exploring the intestinal microecology of sable (Martes zibellina), the national first-level protected animal

Sable (Martes zibellina), a member of family Mustelidae, order Carnivora, is primarily distributed in the cold northern zone of Eurasia. The purpose of this study was to explore the intestinal flora of the sable by metagenomic library-based techniques. Libraries were sequenced on an Illumina HiSeq 4000 instrument. The effective sequencing data of each sample was above 6000 M, and the ratio of clean reads to raw reads was over 98%. The total ORF length was approximately 603,031, equivalent to 347.36 Mbp. We investigated gene functions with the KEGG database and identified 7140 KEGG ortholog (KO) groups comprising 129,788 genes across all of the samples. We selected a subset of genes with the highest abundances to construct cluster heat maps. From the results of the KEGG metabolic pathway annotations, we acquired information on gene functions, as represented by the categories of metabolism, environmental information processing, genetic information processing, cellular processes and organismal systems. We then investigated gene function with the CAZy database and identified functional carbohydrate hydrolases corresponding to genes in the intestinal microorganisms of sable. This finding is consistent with the fact that the sable is adapted to cold environments and requires a large amount of energy to maintain its metabolic activity. We also investigated gene functions with the eggNOG database; the main functions of genes included gene duplication, recombination and repair, transport and metabolism of amino acids, and transport and metabolism of carbohydrates. In this study, we attempted to identify the complex structure of the microbial population of sable based on metagenomic sequencing methods, which use whole metagenomic data, and to map the obtained sequences to known genes or pathways in existing databases, such as CAZy, KEGG, and eggNOG. We then explored the genetic composition and functional diversity of the microbial community based on the mapped functional categories.

Identification and Distribution of Novel Cressdnaviruses and Circular molecules in Four Penguin Species in South Georgia and the Antarctic Peninsula

There is growing interest in uncovering the viral diversity present in wild animal species. The remote Antarctic region is home to a wealth of uncovered microbial diversity, some of which is associated with its megafauna, including penguin species, the dominant avian biota. Penguins interface with a number of other biota in their roles as marine mesopredators and several species overlap in their ranges and habitats. To characterize the circular single-stranded viruses related to those in the phylum Cressdnaviricota from these environmental sentinel species, cloacal swabs (n = 95) were obtained from King Penguins in South Georgia, and congeneric Adélie Penguins, Chinstrap Penguins, and Gentoo Penguins across the South Shetland Islands and Antarctic Peninsula. Using a combination of high-throughput sequencing, abutting primers-based PCR recovery of circular genomic elements, cloning, and Sanger sequencing, we detected 97 novel sequences comprising 40 ssDNA viral genomes and 57 viral-like circular molecules from 45 individual penguins. We present their detection patterns, with Chinstrap Penguins harboring the highest number of new sequences. The novel Antarctic viruses identified appear to be host-specific, while one circular molecule was shared between sympatric Chinstrap and Gentoo Penguins. We also report viral genotype sharing between three adult-chick pairs, one in each Pygoscelid species. Sequence similarity network approaches coupled with Maximum likelihood phylogenies of the clusters indicate the 40 novel viral genomes do not fall within any known viral families and likely fall within the recently established phylum Cressdnaviricota based on their replication-associated protein sequences. Similarly, 83 capsid protein sequences encoded by the viruses or viral-like circular molecules identified in this study do not cluster with any of those encoded by classified viral groups. Further research is warranted to expand knowledge of the Antarctic virome and would help elucidate the importance of viral-like molecules in vertebrate host evolution.

A single-cell polony method reveals low levels of infected Prochlorococcus in oligotrophic waters despite high cyanophage abundances

Long-term stability of picocyanobacteria in the open oceans is maintained by a balance between synchronous division and death on daily timescales. Viruses are considered a major source of microbial mortality, however, current methods to measure infection have significant methodological limitations. Here we describe a method that pairs flow-cytometric sorting with a PCR-based polony technique to simultaneously screen thousands of taxonomically resolved individual cells for intracellular virus DNA, enabling sensitive, high-throughput, and direct quantification of infection by different virus lineages. Under controlled conditions with picocyanobacteria-cyanophage models, the method detected infection throughout the lytic cycle and discriminated between varying infection levels. In North Pacific subtropical surface waters, the method revealed that only a small percentage of Prochlorococcus (0.35–1.6%) were infected, predominantly by T4-like cyanophages, and that infection oscillated 2-fold in phase with the diel cycle. This corresponds to 0.35–4.8% of Prochlorococcus mortality daily. Cyanophages were 2–4-fold more abundant than Prochlorococcus, indicating that most encounters did not result in infection and suggesting infection is mitigated via host resistance, reduced phage infectivity and inefficient adsorption. This method will enable quantification of infection for key microbial taxa across oceanic regimes and will help determine the extent that viruses shape microbial communities and ecosystem level processes.

Viral elements and their potential influence on microbial processes along the permanently stratified Cariaco Basin redoxcline

Little is known about viruses in oxygen-deficient water columns (ODWCs). In surface ocean waters, viruses are known to act as gene vectors among susceptible hosts. Some of these genes may have metabolic functions and are thus termed auxiliary metabolic genes (AMGs). AMGs introduced to new hosts by viruses can enhance viral replication and/or potentially affect biogeochemical cycles by modulating key microbial pathways. Here we identify 748 viral populations that cluster into 94 genera along a vertical geochemical gradient in the Cariaco Basin, a permanently stratified and euxinic ocean basin. The viral communities in this ODWC appear to be relatively novel as 80 of these viral genera contained no reference viral sequences, likely due to the isolation and unique features of this system. We identify viral elements that encode AMGs implicated in distinctive processes, such as sulfur cycling, acetate fermentation, signal transduction, [Fe–S] formation, and N-glycosylation. These AMG-encoding viruses include two putative Mu-like viruses, and viral-like regions that may constitute degraded prophages that have been modified by transposable elements. Our results provide an insight into the ecological and biogeochemical impact of viruses oxygen-depleted and euxinic habitats.

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

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

Viral Metagenomic Content Reflects Seawater Ecological Quality in the Coastal Zone

Viruses interfere with their host’s metabolism through the expression of auxiliary metabolic genes (AMGs) that, until now, are mostly studied under large physicochemical gradients. Here, we focus on coastal marine ecosystems and we sequence the viral metagenome (virome) of samples with discrete levels of human-driven disturbances. We aim to describe the relevance of viromics with respect to ecological quality status, defined by the classic seawater trophic index (TRIX). Neither viral (family level) nor bacterial (family level, based on 16S rRNA sequencing) community structure correlated with TRIX. AMGs involved in the Calvin and tricarboxylic acid cycles were found at stations with poor ecological quality, supporting viral lysis by modifying the host’s energy supply. AMGs involved in “non-traditional” energy-production pathways (3HP, sulfur oxidation) were found irrespective of ecological quality, highlighting the importance of recognizing the prevalent metabolic paths and their intermediate byproducts. Various AMGs explained the variability between stations with poor vs. good ecological quality. Our study confirms the pivotal role of the virome content in ecosystem functioning, acting as a “pool” of available functions that may be transferred to the hosts. Further, it suggests that AMGs could be used as an ultra-sensitive metric of energy-production pathways with relevance in the vulnerable coastal zone and its ecological quality.

Metagenome Data on Intestinal Phage-Bacteria Associations Aids the Development of Phage Therapy against Pathobionts

The application of bacteriophages (phages) is proposed as a highly specific therapy for intestinal pathobiont elimination. However, the infectious associations between phages and bacteria in the human intestine, which is essential information for the development of phage therapies, have yet to be fully elucidated. Here, we report the intestinal viral microbiomes (viromes), together with bacterial microbiomes (bacteriomes), in 101 healthy Japanese individuals. Based on the genomic sequences of bacteriomes and viromes from the same fecal samples, the host bacteria-phage associations are illustrated for both temperate and virulent phages. To verify the usefulness of the comprehensive host bacteria-phage information, we screened Clostridioides difficile-specific phages and identified antibacterial enzymes whose activity is confirmed both in vitro and in vivo. These comprehensive metagenome analyses reveal not only host bacteria-phage associations in the human intestine but also provide vital information for the development of phage therapies against intestinal pathobionts.

Assessing Viral Abundance and Community Composition in Four Contrasting Regions of the Southern Ocean

We explored how changes of viral abundance and community composition among four contrasting regions in the Southern Ocean relied on physicochemical and microbiological traits. During January-February 2015, we visited areas north and south of the South Orkney Islands (NSO and SSO) characterized by low temperature and salinity and high inorganic nutrient concentration, north of South Georgia Island (NSG) and west of Anvers Island (WA), which have relatively higher temperatures and lower inorganic nutrient concentrations. Surface viral abundance (VA) was highest in NSG (21.50 ± 10.70 × 10⁶ viruses mL^-1) and lowest in SSO (2.96 ± 1.48 × 10⁶ viruses mL^-1). VA was positively correlated with temperature, prokaryote abundance and prokaryotic heterotrophic production, chlorophyll a, diatoms, haptophytes, fluorescent organic matter, and isoprene concentration, and was negatively correlated with inorganic nutrients (NO_3-, SiO₄^2-, PO₄^3-), and dimethyl sulfide (DMS) concentrations. Viral communities determined by randomly amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR) were grouped according to the sampling location, being more similar within them than among regions. The first two axes of a canonical correspondence analysis, including physicochemical (temperature, salinity, inorganic nutrients-NO_3-, SiO₄^2-, and dimethyl sulfoniopropionate -DMSP- and isoprene concentrations) and microbiological (chlorophyll a, haptophytes and diatom, and prokaryote abundance and prokaryotic heterotrophic production) factors accounted for 62.9% of the variance. The first axis, temperature-related, accounted for 33.8%; the second one, salinity-related, accounted for 29.1%. Thus, different environmental situations likely select different hosts for viruses, leading to distinct viral communities.

Diverse and unique viruses discovered in the surface water of the East China Sea

Background

Viruses are the most abundant biological entities on earth and play import roles in marine biogeochemical cycles. Here, viral communities in the surface water of the East China Sea (ECS) were collected from three representative regions of Yangshan Harbor (YSH), Gouqi Island (GQI), and the Yangtze River Estuary (YRE) and explored primarily through epifluorescence microscopy (EM), transmission electron microscopy (TEM), and metagenomics analysis.

Results

The virus-like particles (VLPs) in the surface water of the ECS were measured to be 10⁶ to 10⁷ VLPs/ml. Most of the isolated viral particles possessed a head-and-tail structure, but VLPs with unique morphotypes that had never before been observed in the realm of viruses were also found. The sequences related to known viruses in GenBank accounted for 21.1–22.8% of the viromic datasets from YSH, GQI, and YRE. In total, 1029 viral species were identified in the surface waters of the ECS. Among them, tailed phages turn out to make up the majority of viral communities, however a small number of Phycodnaviridae or Mimiviridae related sequences were also detected. The diversity of viruses did not appear to be a big difference among these three aquatic environments but their relative abundance was geographically variable. For example, the Pelagibacter phage HTVC010P accounted for 50.4% of the identified viral species in GQI, but only 9.1% in YSH and 11.7% in YRE. Sequences, almost identical to those of uncultured marine thaumarchaeal dsDNA viruses and magroviruses that infect Marine Group II Euryarchaeota, were confidently detected in the ECS viromes. The predominant classes of virome ORFs with functional annotations that were found were those involved in viral biogenesis. Virus-host connections, inferred from CRISPR spacer-protospacer mapping, implied newly discovered infection relationships in response to arms race between them.

Conclusions

Together, both identified viruses and unknown viral assemblages observed in this study were indicative of the complex viral community composition found in the ECS. This finding fills a major gap in the dark world of oceanic viruses of China and additionally contributes to the better understanding of global marine viral diversity, composition, and distribution.

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

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

Initial Mapping of the New York City Wastewater Virome

Wastewater is a rich source of microbial life and contains bacteria, viruses, and other microbes found in human waste as well as environmental runoff sources. As part of an effort to characterize the New York City wastewater metagenome, we profiled the viral community of sewage samples across all five boroughs of NYC and found that local sampling sites have unique sets of viruses. We focused on bacteriophages, or viruses of bacteria, to understand how they may influence the microbial ecology of this system. We identified several new clusters of phages and successfully associated them with bacterial hosts, providing insight into virus-host interactions in urban wastewater. This study provides a first look into the viral communities present across the wastewater system in NYC and points to their functional importance in this environment.

Bacteriophages are abundant members of all microbiomes studied to date, influencing microbial communities through interactions with their bacterial hosts. Despite their functional importance and ubiquity, phages have been underexplored in urban environments compared to their bacterial counterparts. We profiled the viral communities in New York City (NYC) wastewater using metagenomic data collected in November 2014 from 14 wastewater treatment plants. We show that phages accounted for the largest viral component of the sewage samples and that specific virus communities were associated with local environmental conditions within boroughs. The vast majority of the virus sequences had no homology matches in public databases, forming an average of 1,700 unique virus clusters (putative genera). These new clusters contribute to elucidating the overwhelming proportion of data that frequently goes unidentified in viral metagenomic studies. We assigned potential hosts to these phages, which appear to infect a wide range of bacterial genera, often outside their presumed host. We determined that infection networks form a modular-nested pattern, indicating that phages include a range of host specificities, from generalists to specialists, with most interactions organized into distinct groups. We identified genes in viral contigs involved in carbon and sulfur cycling, suggesting functional importance of viruses in circulating pathways and gene functions in the wastewater environment. In addition, we identified virophage genes as well as a nearly complete novel virophage genome. These findings provide an understanding of phage abundance and diversity in NYC wastewater, previously uncharacterized, and further examine geographic patterns of phage-host association in urban environments.

IMPORTANCE Wastewater is a rich source of microbial life and contains bacteria, viruses, and other microbes found in human waste as well as environmental runoff sources. As part of an effort to characterize the New York City wastewater metagenome, we profiled the viral community of sewage samples across all five boroughs of NYC and found that local sampling sites have unique sets of viruses. We focused on bacteriophages, or viruses of bacteria, to understand how they may influence the microbial ecology of this system. We identified several new clusters of phages and successfully associated them with bacterial hosts, providing insight into virus-host interactions in urban wastewater. This study provides a first look into the viral communities present across the wastewater system in NYC and points to their functional importance in this environment.

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.