Unveiling the role and life strategies of viruses from the surface to the dark ocean

DeepMicroClass sorts metagenomic contigs into prokaryotes, eukaryotes and viruses

Sequence classification facilitates a fundamental understanding of the structure of microbial communities. Binary metagenomic sequence classifiers are insufficient because environmental metagenomes are typically derived from multiple sequence sources. Here we introduce a deep-learning based sequence classifier, DeepMicroClass, that classifies metagenomic contigs into five sequence classes, i.e. viruses infecting prokaryotic or eukaryotic hosts, eukaryotic or prokaryotic chromosomes, and prokaryotic plasmids. DeepMicroClass achieved high performance for all sequence classes at various tested sequence lengths ranging from 500 bp to 100 kbps. By benchmarking on a synthetic dataset with variable sequence class composition, we showed that DeepMicroClass obtained better performance for eukaryotic, plasmid and viral contig classification than other state-of-the-art predictors. DeepMicroClass achieved comparable performance on viral sequence classification with geNomad and VirSorter2 when benchmarked on the CAMI II marine dataset. Using a coastal daily time-series metagenomic dataset as a case study, we showed that microbial eukaryotes and prokaryotic viruses are integral to microbial communities. By analyzing monthly metagenomes collected at HOT and BATS, we found relatively higher viral read proportions in the subsurface layer in late summer, consistent with the seasonal viral infection patterns prevalent in these areas. We expect DeepMicroClass will promote metagenomic studies of under-appreciated sequence types.

Giant viral signatures on the Greenland ice sheet

BACKGROUND

Dark pigmented snow and glacier ice algae on glaciers and ice sheets contribute to accelerating melt. The biological controls on these algae, particularly the role of viruses, remain poorly understood. Giant viruses, classified under the nucleocytoplasmic large DNA viruses (NCLDV) supergroup (phylum Nucleocytoviricota), are diverse and globally distributed. NCLDVs are known to infect eukaryotic cells in marine and freshwater environments, providing a biological control on the algal population in these ecosystems. However, there is very limited information on the diversity and ecosystem function of NCLDVs in terrestrial icy habitats.

RESULTS

In this study, we investigate for the first time giant viruses and their host connections on ice and snow habitats, such as cryoconite, dark ice, ice core, red and green snow, and genomic assemblies of five cultivated Chlorophyta snow algae. Giant virus marker genes were present in almost all samples; the highest abundances were recovered from red snow and the snow algae genomic assemblies, followed by green snow and dark ice. The variety of active algae and protists in these GrIS habitats containing NCLDV marker genes suggests that infection can occur on a range of eukaryotic hosts. Metagenomic data from red and green snow contained evidence of giant virus metagenome-assembled genomes from the orders Imitervirales, Asfuvirales, and Algavirales.

CONCLUSION

Our study highlights NCLDV family signatures in snow and ice samples from the Greenland ice sheet. Giant virus metagenome-assembled genomes (GVMAGs) were found in red snow samples, and related NCLDV marker genes were identified for the first time in snow algal culture genomic assemblies; implying a relationship between the NCLDVs and snow algae. Metatranscriptomic viral genes also aligned with metagenomic sequences, suggesting that NCLDVs are an active component of the microbial community and are potential "top-down" controls of the eukaryotic algal and protistan members. This study reveals the unprecedented presence of a diverse community of NCLDVs in a variety of glacial habitats dominated by algae.

Diversities and interactions of phages and bacteria in deep-sea sediments as revealed by metagenomics

Phages are found virtually everywhere, even in extreme environments, and are extremely diverse both in their virion structures and in their genomic content. They are thought to shape the taxonomic and functional composition of microbial communities as well as their stability. A number of studies on laboratory culture and viral metagenomic research provide deeper insights into the abundance, diversity, distribution, and interaction with hosts of phages across a wide range of ecosystems. Although most of these studies focus on easily accessible samples, such as soils, lakes, and shallow oceans, little is known about bathypelagic phages. In this study, through analyzing the 16S rRNA sequencing and viral metagenomic sequencing data of 25 samples collected from five different bathypelagic ecosystems, we detected a high diversity of bacteria and phages, particularly in the cold seep and hydrothermal vent ecosystems, which have stable chemical energy. The relative abundance of phages in these ecosystems was higher than in other three abyssal ecosystems. The low phage/host ratios obtained from host prediction were different from shallow ecosystems and indicated the prevalence of prophages, suggesting the complexity of phage-bacteria interactions in abyssal ecosystems. In the correlation analysis, we revealed several phages-bacteria interaction networks of potential ecological relevance. Our study contributes to a better understanding of the interactions between bathypelagic bacteria and their phages.

Phylogenomic analysis expands the known repertoire of single-stranded DNA viruses in benthic zones of the South Indian Ocean

Single-stranded (ss) DNA viruses are ubiquitous and constitute some of the most diverse entities on Earth. Most studies have focused on ssDNA viruses from terrestrial environments resulting in a significant deficit in benthic ecosystems including aphotic zones of the South Indian Ocean (SIO). Here, we assess the diversity and phylogeny of ssDNA in deep waters of the SIO using a combination of established viral taxonomy tools and a Hidden Markov Model based approach. Replication initiator protein-associated (Rep) phylogenetic reconstruction and sequence similarity networks were used to show that the SIO hosts divergent and as yet unknown circular Rep-encoding ssDNA viruses. Several sequences appear to represent entirely novel families, expanding the repertoire of known ssDNA viruses. Results suggest that a small proportion of these viruses may be circular genetic elements, which may strongly influence the diversity of both eukaryotes and prokaryotes in the SIO. Taken together, our data show that the SIO harbours a diverse assortment of previously unknown ssDNA viruses. Due to their potential to infect a variety of hosts, these viruses may be crucial for marine nutrient recycling through their influence of the biological carbon pump.

Viral enumeration using cost-effective wet-mount epifluorescence microscopy for aquatic ecosystems and modern microbialites

IMPORTANCE

Low-cost and robust viral enumeration is a critical first step toward understanding the global virome. Our method is a deep drive integration providing a window into viral dark matter within aquatic ecosystems. We enumerated the viruses within Green Lake and Great Salt Lake microbialites, EPS, and water column. The entire weight of all the viruses in Green Lake and Great Salt Lake are ~598 g and ~2.2 kg, respectively.

Viruses under the Antarctic Ice Shelf are active and potentially involved in global nutrient cycles

Viruses play an important role in the marine ecosystem. However, our comprehension of viruses inhabiting the dark ocean, and in particular, under the Antarctic Ice Shelves, remains limited. Here, we mine single-cell genomic, transcriptomic, and metagenomic data to uncover the viral diversity, biogeography, activity, and their role as metabolic facilitators of microbes beneath the Ross Ice Shelf. This is the largest Antarctic ice shelf with a major impact on global carbon cycle. The viral community found in the cavity under the ice shelf mainly comprises endemic viruses adapted to polar and mesopelagic environments. The low abundance of genes related to lysogenic lifestyle (<3%) does not support a predominance of the Piggyback-the-Winner hypothesis, consistent with a low-productivity habitat. Our results indicate a viral community actively infecting key ammonium and sulfur-oxidizing chemolithoautotrophs (e.g. Nitrosopumilus spp, Thioglobus spp.), supporting a "kill-the-winner" dynamic. Based on genome analysis, these viruses carry specific auxiliary metabolic genes potentially involved in nitrogen, sulfur, and phosphorus acquisition. Altogether, the viruses under Antarctic ice shelves are putatively involved in programming the metabolism of ecologically relevant microbes that maintain primary production in these chemosynthetically-driven ecosystems, which have a major role in global nutrient cycles.

Novel RNA viruses from the Atlantic Ocean: Ecogenomics, biogeography, and total virioplankton mass contribution from surface to the deep ocean

Marine viruses play a major role in the energy and nutrient cycle and affect the evolution of their hosts. Despite their importance, there is still little knowledge about RNA viruses. Here, we have explored the Atlantic Ocean, from surface to deep (4.296 m), and used viromics and quantitative methods to unveil the genomics, biogeography, and the mass contribution of RNA viruses to the total viroplankton. A total of 2481 putative RNA viral contigs (>500 bp) and 107 larger bona fide RNA viral genomes (>2.5 kb) were identified; 88 of them representing novel viruses belonging mostly to two clades: Yangshan assemblage (sister clade to the class Alsuviricetes) and Nodaviridae. These viruses were highly endemic and locally abundant, with little or no presence in other oceans since only ≈15% of them were found in at least one of the Tara sampling metatranscriptomes. Quantitative data indicated that the abundance of RNA viruses in the surface and deep chlorophyll maximum zone was within ≈106  VLP/mL representing a potential contribution of 5.2%-24.4% to the total viroplankton community (DNA and RNA viruses), with DNA viruses being the predominant members (≈107  VLP/mL). However, for the deep sample, the observed trend was the opposite, although as further discussed, several biases should be considered. Together these results contribute to our understanding of the diversity, abundance, and distribution of RNA viruses in the oceans and provide a basis for further investigation into their ecological roles and biogeography.

Depth-driven patterns in lytic viral diversity, auxiliary metabolic gene content, and productivity in offshore oligotrophic waters

Introduction

Marine viruses regulate microbial population dynamics and biogeochemical cycling in the oceans. The ability of viruses to manipulate hosts' metabolism through the expression of viral auxiliary metabolic genes (AMGs) was recently highlighted, having important implications in energy production and flow in various aquatic environments. Up to now, the presence and diversity of viral AMGs is studied using -omics data, and rarely using quantitative measures of viral activity alongside.

Methods

In the present study, four depth layers (5, 50, 75, and 1,000 m) with discrete hydrographic features were sampled in the Eastern Mediterranean Sea; we studied lytic viral community composition and AMG content through metagenomics, and lytic production rates through the viral reduction approach in the ultra-oligotrophic Levantine basin where knowledge regarding viral actions is rather limited.

Results and Discussion

Our results demonstrate depth-dependent patterns in viral diversity and AMG content, related to differences in temperature, nutrients availability, and host bacterial productivity and abundance. Although lytic viral production rates were similar along the water column, the virus-to-bacteria ratio was higher and the particular set of AMGs was more diverse in the bathypelagic (1,000 m) than the shallow epipelagic (5, 50, and 75 m) layers, revealing that the quantitative effect of viruses on their hosts may be the same along the water column through the intervention of different AMGs. In the resource- and energy-limited bathypelagic waters of the Eastern Mediterranean, the detected AMGs could divert hosts' metabolism toward energy production, through a boost in gluconeogenesis, fatty-acid and glycan biosynthesis and metabolism, and sulfur relay. Near the deep-chlorophyll maximum depth, an exceptionally high percentage of AMGs related to photosynthesis was noticed. Taken together our findings suggest that the roles of viruses in the deep sea might be even more important than previously thought as they seem to orchestrate energy acquisition and microbial community dynamics, and thus, biogeochemical turnover in the oceans.

Ubiquitous, B12-dependent virioplankton utilizing ribonucleotide-triphosphate reductase demonstrate interseasonal dynamics and associate with a diverse range of bacterial hosts in the pelagic ocean

Through infection and lysis of their coexisting bacterial hosts, viruses impact the biogeochemical cycles sustaining globally significant pelagic oceanic ecosystems. Currently, little is known of the ecological interactions between lytic viruses and their bacterial hosts underlying these biogeochemical impacts at ecosystem scales. This study focused on populations of lytic viruses carrying the B12-dependent Class II monomeric ribonucleotide reductase (RNR) gene, ribonucleotide-triphosphate reductase (Class II RTPR), documenting seasonal changes in pelagic virioplankton and bacterioplankton using amplicon sequences of Class II RTPR and the 16S rRNA gene, respectively. Amplicon sequence libraries were analyzed using compositional data analysis tools that account for the compositional nature of these data. Both virio- and bacterioplankton communities responded to environmental changes typically seen across seasonal cycles as well as shorter term upwelling-downwelling events. Defining Class II RTPR-carrying viral populations according to major phylogenetic clades proved a more robust means of exploring virioplankton ecology than operational taxonomic units defined by percent sequence homology. Virioplankton Class II RTPR populations showed positive associations with a broad phylogenetic diversity of bacterioplankton including dominant taxa within pelagic oceanic ecosystems such as Prochlorococcus and SAR11. Temporal changes in Class II RTPR virioplankton, occurring as both free viruses and within infected cells, indicated possible viral-host pairs undergoing sustained infection and lysis cycles throughout the seasonal study. Phylogenetic relationships inferred from Class II RTPR sequences mirrored ecological patterns in virio- and bacterioplankton populations demonstrating possible genome to phenome associations for an essential viral replication gene.

Seasonality of top-down control of bacterioplankton at two central Red Sea sites with different trophic status

The role of bottom-up (nutrient availability) and top-down (grazers and viruses mortality) controls on tropical bacterioplankton have been rarely investigated simultaneously from a seasonal perspective. We have assessed them through monthly samplings over 2 years in inshore and offshore waters of the central Red Sea differing in trophic status. Flow cytometric analysis allowed us to distinguish five groups of heterotrophic bacteria based on physiological properties (nucleic acid content, membrane integrity and active respiration), three groups of cyanobacteria (two populations of Synechococcus and Prochlorococcus), heterotrophic nanoflagellates (HNFs) and three groups of viruses based on nucleic acid content. The dynamics of bacterioplankton and their top-down controls varied with season and location, being more pronounced in inshore waters. HNFs abundances showed a strong preference for larger prey inshore (r = -0.62 to -0.59, p = 0.001-0.002). Positive relationships between viruses and heterotrophic bacterioplankton abundances were more marked inshore (r = 0.67, p < 0.001) than offshore (r = 0.44, p = 0.03). The negative correlation between HNFs and viruses abundances (r = -0.47, p = 0.02) in shallow waters indicates a persistent seasonal switch between protistan grazing and viral lysis that maintains the low bacterioplankton stocks in the central Red Sea area.

Virome reveals effect of Ulva prolifera green tide on the structural and functional profiles of virus communities in coastal environments

Viruses are widely distributed in marine environments, where they influence the transformation of matter and energy by modulating host metabolism. Driven by eutrophication, green tides are a rising concern in Chinese coastal areas, and are a serious ecological disaster that negatively affects coastal ecosystems and disrupts biogeochemical cycles. Although the composition of bacterial communities in green algae has been investigated, the diversity and roles of viruses in green algal blooms are largely unexplored. Therefore, the diversity, abundance, lifestyle, and metabolic potential of viruses in a natural bloom in Qingdao coastal area were investigated at three different stages (pre-bloom, during-bloom, and post-bloom) by metagenomics analysis. The dsDNA viruses, Siphoviridae, Myoviridae, Podoviridae, and Phycodnaviridae, were found to dominate the viral community. The viral dynamics exhibited distinct temporal patterns across different stages. The composition of the viral community varied during the bloom, especially in populations with low abundance. The lytic cycle was most predominant, and the abundance of lytic viruses increased slightly in the post-bloom stage. The diversity and richness of the viral communities varied distinctly during the green tide, and the post-bloom stage favored viral diversity and richness. The total organic carbon, dissolved oxygen, NO^3-, NO^2-, PO₄^3-, chlorophyll-a contents, and temperature variably co-influenced the viral communities. The primary hosts included bacteria, algae, and other microplankton. Network analysis revealed the closer links between the viral communities as the bloom progressed. Functional prediction revealed that the viruses possibly influenced the biodegradation of microbial hydrocarbons and carbon by metabolic augmentation via auxiliary metabolic genes. The composition, structure, metabolic potential, and interaction taxonomy of the viromes differed significantly across the different stages of the green tide. The study demonstrated that the ecological event shaped the viral communities during algal bloom, and the viral communities played a significant role in phycospheric microecology.

Metagenomic Analyses Reveal the Influence of Depth Layers on Marine Biodiversity on Tropical and Subtropical Regions

The emergence of open ocean global-scale studies provided important information about the genomics of oceanic microbial communities. Metagenomic analyses shed light on the structure of marine habitats, unraveling the biodiversity of different water masses. Many biological and environmental factors can contribute to marine organism composition, such as depth. However, much remains unknown about microbial communities' taxonomic and functional features in different water layer depths. Here, we performed a metagenomic analysis of 76 publicly available samples from the Tara Ocean Project, distributed in 8 collection stations located in tropical or subtropical regions, and sampled from three layers of depth (surface water layer-SRF, deep chlorophyll maximum layer-DCM, and mesopelagic zone-MES). The SRF and DCM depth layers are similar in abundance and diversity, while the MES layer presents greater diversity than the other layers. Diversity clustering analysis shows differences regarding the taxonomic content of samples. At the domain level, bacteria prevail in most samples, and the MES layer presents the highest proportion of archaea among all samples. Taken together, our results indicate that the depth layer influences microbial sample composition and diversity.

Viral lysing can alleviate microbial nutrient limitations and accumulate recalcitrant dissolved organic matter components in soil

Viruses are critical for regulating microbial communities and biogeochemical processes affecting carbon/nutrient cycling. However, the role of soil phages in controlling microbial physiological traits and intrinsic dissolved organic matter (DOM) properties remains largely unknown. Herein, microcosm experiments with different soil phage concentrates (including no-added phages, inactive phages, and three dilutions of active phages) at two temperatures (15 °C and 25 °C) were conducted to disclose the nutrient and DOM dynamics associated with viral lysing. Results demonstrated three different phases of viral impacts on CO₂ emission at both temperatures, and phages played a role in maintaining Q₁₀ within bounds. At both temperatures, microbial nutrient limitations (especially P limitation) were alleviated by viral lysing as determined by extracellular enzyme activity (decreased V_angle with active phages). Additionally, the re-utilization of lysate-derived DOM by surviving microbes stimulated an increase of microbial metabolic efficiency and recalcitrant DOM components (e.g., SUV254, SUV260 and HIX). This research provides direct experimental evidence that the "viral shuttle" exists in soils, whereby soil phages increase recalcitrant DOM components. Our findings advance the understanding of viral controls on soil biogeochemical processes, and provide a new perspective for assessing whether soil phages provide a net "carbon sink" vs. "carbon source" in soils.

Water mass age structures the auxiliary metabolic gene content of free-living and particle-attached deep ocean viral communities

BACKGROUND

Viruses play important roles in the ocean's biogeochemical cycles. Yet, deep ocean viruses are one of the most under-explored fractions of the global biosphere. Little is known about the environmental factors that control the composition and functioning of their communities or how they interact with their free-living or particle-attached microbial hosts.

RESULTS

We analysed 58 viral communities associated with size-fractionated free-living (0.2-0.8 μm) and particle-attached (0.8-20 μm) cellular metagenomes from bathypelagic (2150-4018 m deep) microbiomes obtained during the Malaspina expedition. These metagenomes yielded 6631 viral sequences, 91% of which were novel, and 67 represented high-quality genomes. Taxonomic classification assigned 53% of the viral sequences to families of tailed viruses from the order Caudovirales. Computational host prediction associated 886 viral sequences to dominant members of the deep ocean microbiome, such as Alphaproteobacteria (284), Gammaproteobacteria (241), SAR324 (23), Marinisomatota (39), and Chloroflexota (61). Free-living and particle-attached viral communities had markedly distinct taxonomic composition, host prevalence, and auxiliary metabolic gene content, which led to the discovery of novel viral-encoded metabolic genes involved in the folate and nucleotide metabolisms. Water mass age emerged as an important factor driving viral community composition. We postulated this was due to changes in quality and concentration of dissolved organic matter acting on the host communities, leading to an increase of viral auxiliary metabolic genes associated with energy metabolism among older water masses.

CONCLUSIONS

These results shed light on the mechanisms by which environmental gradients of deep ocean ecosystems structure the composition and functioning of free-living and particle-attached viral communities. Video Abstract.

Characterization and genomic analysis of a novel Synechococcus phage S-H9-2 belonging to Bristolvirus genus isolated from the Yellow Sea

Cyanophages are known to influence the population dynamics and community structure of cyanobacteria and thus play an important role in biogeochemical cycles in aquatic ecosystems. In this study, a novel Synechococcus phage S-H9-2 infecting Synechococcus sp. WH 8102 was isolated from the coastal water of the Yellow Sea. Synechococcus phage S-H9-2 contains a 187,320 bp genome of double-stranded DNA with a G + C content of 40.3%, 202 potential open reading frames (ORFs), and 15 tRNAs. Phylogenetic analysis and nucleotide-based intergenomic similarity suggest that Synechococcus phage S-H9-2 belongs to the Bristolvirus genus under the family Kyanoviridae. Homologs of the S-H9-2 open reading frame can be found in a variety of marine environments, as shown by the results of mapping the genome sequence of S-H9-2 to the Global Ocean Viromes 2.0 dataset. The presence of auxiliary metabolic genes (AMGs) related to photosynthesis, carbon metabolism, and phosphorus assimilation, as well as phylogenetic relationships based on complete genome sequences, reflect the mechanism of phage-host interaction and host-specific strategies for adaptation to environmental conditions. This study enriches the current genomic database of cyanophage and contributed to our understanding of the virus-host interactions and their adaption to the environment.

The smallest in the deepest: the enigmatic role of viruses in the deep biosphere

It is commonly recognized that viruses control the composition, metabolism, and evolutionary trajectories of prokaryotic communities, with resulting vital feedback on ecosystem functioning and nutrient cycling in a wide range of ecosystems. Although the deep biosphere has been estimated to be the largest reservoir for viruses and their prokaryotic hosts, the biology and ecology of viruses therein remain poorly understood. The deep virosphere is an enigmatic field of study in which many critical questions are still to be answered. Is the deep virosphere simply a repository for deeply preserved, non-functioning virus particles? Or are deep viruses infectious agents that can readily infect suitable hosts and subsequently shape microbial populations and nutrient cycling? Can the cellular content released by viral lysis, and even the organic structures of virions themselves, serve as the source of bioavailable nutrients for microbial activity in the deep biosphere as in other ecosystems? In this review, we synthesize our current knowledge of viruses in the deep biosphere and seek to identify topics with the potential for substantial discoveries in the future.

Awaking the dormant virome in the rhizosphere

The rhizosphere is a vital soil compartment providing key plant-beneficial functions. However, little is known about the mechanisms driving viral diversity in the rhizosphere. Viruses can establish lytic or lysogenic interactions with their bacterial hosts. In the latter, they assume a dormant state integrated in the host genome and can be awakened by different perturbations that impact host cell physiology, triggering a viral bloom, which is potentially a fundamental mechanism driving soil viral diversity, as 22%-68% of soil bacteria are predicted to harbour dormant viruses. Here we assessed the viral bloom response in rhizospheric viromes by exposing them to three contrasting soil perturbation agents: earthworms, herbicide and antibiotic pollutant. The viromes were next screened for rhizosphere-relevant genes and also used as inoculant on microcosms incubations to test their impacts on pristine microbiomes. Our results show that while post-perturbation viromes diverged from control conditions, viral communities exposed to both herbicide and antibiotic pollutant were more similar to each other than those influenced by earthworms. The latter also favoured an increase in viral populations harbouring genes involved in plant-beneficial functions. Post-perturbation viromes inoculated on soil microcosms changed the diversity of pristine microbiomes, suggesting that viromes are important components of the soil ecological memory driving eco-evolutionary processes that determine future microbiome trajectories according to past events. Our findings demonstrate that viromes are active players in the rhizosphere and need to be considered in efforts to understand and control the microbial processes towards sustainable crop production.

Over two decades of research on the marine RNA virosphere

RNA viruses (realm: Riboviria), including RNA phages and eukaryote-infecting RNA viruses, are essential components of marine ecosystems. A large number of marine RNA viruses have been discovered in the last two decades because of the rapid development of next-generation sequencing (NGS) technology. Indeed, the combination of NGS and state-of-the-art meta-omics methods (viromics, the study of all viruses in a specific environment) has led to a fundamental understanding of the taxonomy and genetic diversity of RNA viruses in the sea, suggesting the complex ecological roles played by RNA viruses in this complex ecosystem. Furthermore, comparisons of viromes in the context of highly variable marine niches reveal the biogeographic patterns and ecological impact of marine RNA viruses, whose role in global ecology is becoming increasingly clearer. In this review, we summarize the characteristics of the global marine RNA virosphere and outline the taxonomic hierarchy of RNA viruses with a specific focus on their ancient evolutionary history. We also review the development of methodology and the major progress resulting from its applications in RNA viromics. The aim of this review is not only to provide an in-depth understanding of multifaceted aspects of marine RNA viruses, but to offer future perspectives on developing a better methodology for discovery, and exploring the evolutionary origin and major ecological significance of marine RNA virosphere.

Microbial ecology of the Southern Ocean

The Southern Ocean (SO) distributes climate signals and nutrients worldwide, playing a pivotal role in global carbon sequestration. Microbial communities are essential mediators of primary productivity and carbon sequestration, yet we lack a comprehensive understanding of microbial diversity and functionality in the SO. Here, we examine contemporary studies in this unique polar system, focusing on prokaryotic communities and their relationships with other trophic levels (i.e. phytoplankton and viruses). Strong seasonal variations and the characteristic features of this ocean are directly linked to community composition and ecosystem functions. Specifically, we discuss characteristics of SO microbial communities and emphasise differences from the Arctic Ocean microbiome. We highlight the importance of abundant bacteria in recycling photosynthetically derived organic matter. These heterotrophs appear to control carbon flux to higher trophic levels when light and iron availability favour primary production in spring and summer. Conversely, during winter, evidence suggests that chemolithoautotrophs contribute to prokaryotic production in Antarctic waters. We conclude by reviewing the effects of climate change on marine microbiota in the SO.

Marine viral particles reveal an expansive repertoire of phage-parasitizing mobile elements

Phage satellites are mobile genetic elements that parasitize viruses, exerting profound biological and ecological impacts. Phage satellites are known to infect several gram-positive genera and a few gram-negative bacterial species, most associated with the human microbiome. Direct inspection of “wild” virus particles, however, revealed that marine phage satellites are widely distributed and abundant in the global oceans. Their genetic diversity, gene repertoires, and host ranges appear much greater than has been previously reported. Genetic analyses now provide clues about the parasitic life cycles, helper bacteriophage interactions, and reproductive strategies of these newly recognized marine phage satellites. Their properties, diversity, and environmental distributions suggest they may exert substantial influence on microbial ecology and evolution in the sea.

Phage satellites are mobile genetic elements that propagate by parasitizing bacteriophage replication. We report here the discovery of abundant and diverse phage satellites that were packaged as concatemeric repeats within naturally occurring bacteriophage particles in seawater. These same phage-parasitizing mobile elements were found integrated in the genomes of dominant co-occurring bacterioplankton species. Like known phage satellites, many marine phage satellites encoded genes for integration, DNA replication, phage interference, and capsid assembly. Many also contained distinctive gene suites indicative of unique virus hijacking, phage immunity, and mobilization mechanisms. Marine phage satellite sequences were widespread in local and global oceanic virioplankton populations, reflecting their ubiquity, abundance, and temporal persistence in marine planktonic communities worldwide. Their gene content and putative life cycles suggest they may impact host-cell phage immunity and defense, lateral gene transfer, bacteriophage-induced cell mortality and cellular host and virus productivity. Given that marine phage satellites cannot be distinguished from bona fide viral particles via commonly used microscopic techniques, their predicted numbers (∼3.2 × 10²⁶ in the ocean) may influence current estimates of virus densities, production, and virus-induced mortality. In total, the data suggest that marine phage satellites have potential to significantly impact the ecology and evolution of bacteria and their viruses throughout the oceans. We predict that any habitat that harbors bacteriophage will also harbor similar phage satellites, making them a ubiquitous feature of most microbiomes on Earth.

Viral shunt in tropical oligotrophic ocean

Viruses cause massive bacterial mortality and thus modulate bacteria-governed carbon transfer and nutrient recycling at global scale. The viral shunt hypothesis states the crucial role of viral lysis in retaining microbial carbon into food web processes, while its applicability to nature has not been well identified for over two decades. Here, we conducted nine diel surveys in the tropical South China Sea and suggested that the time scale adopted in sampling and system trophic status determine the "visibility" of the viral shunt in the field. Specifically, viral abundance (VA), bacterial biomass (BB), and bacterial specific growth rate (SGR) varied synchronously and presented the significant VA-BB and VA-SGR linkages at an hourly scale, which reveals direct interactions between viruses and their hosts. The differential responses of the viral shunt to temperature, i.e., looser VA-SGR coupling in warm and tighter VA-SGR coupling in cold environments, imply an altered carbon cycling in tropical oceans under climatic warming.

The Life Cycle Transitions of Temperate Phages: Regulating Factors and Potential Ecological Implications

Phages are viruses that infect bacteria. They affect various microbe-mediated processes that drive biogeochemical cycling on a global scale. Their influence depends on whether the infection is lysogenic or lytic. Temperate phages have the potential to execute both infection types and thus frequently switch their infection modes in nature, potentially causing substantial impacts on the host-phage community and relevant biogeochemical cycling. Understanding the regulating factors and outcomes of temperate phage life cycle transition is thus fundamental for evaluating their ecological impacts. This review thus systematically summarizes the effects of various factors affecting temperate phage life cycle decisions in both culturable phage-host systems and natural environments. The review further elucidates the ecological implications of the life cycle transition of temperate phages with an emphasis on phage/host fitness, host-phage dynamics, microbe diversity and evolution, and biogeochemical cycles.

Behavioral responses to riparian and anthropogenic edge effects in mantled howler monkeys (Alouatta palliata) in a disturbed riverine forest

Fragmented forests contain natural edges, including riparian zones, and anthropogenic edges. Edges generally have lower plant density and fewer large trees than forest interior. Riparian edges, however, contain gap-specialist trees yielding leaves with high protein content, providing primates with important resources. We examined mantled howler monkeys' behavioral responses to riparian and anthropogenic edges at La Suerte Biological Research Station (LSBRS), Costa Rica. We predicted the monkeys would spend more time resting and feeding and less time traveling, and be less spatially cohesive, in both anthropogenic and riparian edges compared to forest interior due to lower resource abundance in edges, and in anthropogenic compared to riparian edge due to higher leaf quality in riparian zones. From 2017 to 2020, we collected data across forest zones on activity and spatial cohesion patterns via focal sampling, recording data every 2 min. Howler monkeys were significantly more likely to rest and significantly less likely to travel in both anthropogenic and riparian edges compared to forest interior; however, there were no differences between these edge types. There were significantly more monkeys within a 5-m radius of focal subjects in both anthropogenic and riparian edges compared to forest interior, but no differences between these edge types. While prior research found no differences across zones when only anthropogenic edge and forest interior were compared, results of this study demonstrate that howler monkeys at LSBRS modify their activity patterns in anthropogenic and riparian edge zones compared to forest interior, highlighting the importance of focusing on both natural and anthropogenic edge zones to fully understand primates' behavioral responses in fragmented landscapes.

Virus-to-prokaryote ratio in the Salar de Huasco and different ecosystems of the Southern hemisphere and its relationship with physicochemical and biological parameters

The virus-to-prokaryote ratio (VPR) has been used in many ecosystems to study the relationship between viruses and their hosts. While high VPR values indicate a high rate of prokaryotes' cell lysis, low values are interpreted as a decrease in or absence of viral activity. Salar de Huasco is a high-altitude wetland characterized by a rich microbial diversity associated with aquatic sites like springs, ponds, streams and a lagoon with variable physicochemical conditions. Samples from two ponds, Poza Rosada (PR) and Poza Verde (PV), were analyzed by epifluorescence microscopy to determine variability of viral and prokaryotic abundance and to calculate the VPR in a dry season. In addition, to put Salar de Huasco results into perspective, a compilation of research articles on viral and prokaryotic abundance, VPR, and metadata from various Southern hemisphere ecosystems was revised. The ecosystems were grouped into six categories: high-altitude wetlands, Pacific, Atlantic, Indian, and Southern Oceans and Antarctic lakes. Salar de Huasco ponds recorded similar VPR values (an average of 7.4 and 1.7 at PR and PV, respectively), ranging from 3.22 to 15.99 in PR. The VPR variability was associated with VA and chlorophyll a, when considering all data available for this ecosystem. In general, high-altitude wetlands recorded the highest VPR average (53.22 ± 95.09), followed by the Oceans, Southern (21.91 ± 25.72), Atlantic (19.57 ± 15.77) and Indian (13.43 ± 16.12), then Antarctic lakes (11.37 ± 15.82) and the Pacific Ocean (6.34 ± 3.79). Physicochemical variables, i.e., temperature, conductivity, nutrients (nitrate, ammonium, and phosphate) and chlorophyll a as a biological variable, were found to drive the VPR in the ecosystems analyzed. Thus, the viral activity in the Wetland followed similar trends of previous reports based on larger sets of metadata analyses. In total, this study highlights the importance of including viruses as a biological variable to study microbial temporal dynamics in wetlands considering their crucial role in the carbon budgets of these understudied ecosystems in the southern hemisphere.

Virioplankton assemblages from challenger deep, the deepest place in the oceans

Hadal ocean biosphere, that is, the deepest part of the world’s oceans, harbors a unique microbial community, suggesting a potential uncovered co-occurring virioplankton assemblage. Herein, we reveal the unique virioplankton assemblages of the Challenger Deep, comprising 95,813 non-redundant viral contigs from the surface to the hadal zone. Almost all of the dominant viral contigs in the hadal zone were unclassified, potentially related to Alteromonadales and Oceanospirillales. 2,586 viral auxiliary metabolic genes from 132 different KEGG orthologous groups were mainly related to the carbon, nitrogen, sulfur, and arsenic metabolism. Lysogenic viral production and integrase genes were augmented in the hadal zone, suggesting the prevalence of viral lysogenic life strategy. Abundant rve genes in the hadal zone, which function as transposase in the caudoviruses, further suggest the prevalence of viral-mediated horizontal gene transfer. This study provides fundamental insights into the virioplankton assemblages of the hadal zone, reinforcing the necessity of incorporating virioplankton into the hadal biogeochemical cycles.

•

The unique virioplankton assemblages of the Challenger Deep were revealed

•

Virus encoded auxiliary metabolic genes relating to the biogeochemical cycling

•

Viruses in deep and hadal zone tend to be lysogenic, and potentially mediate the horizontal gene transfer

Into the Dark: Exploring the Deep Ocean with Single-Virus Genomics

Single-virus genomics (SVGs) has been successfully applied to ocean surface samples allowing the discovery of widespread dominant viruses overlooked for years by metagenomics, such as the uncultured virus vSAG 37-F6 infecting the ubiquitous Pelagibacter spp. In SVGs, one uncultured virus at a time is sorted from the environmental sample, whole-genome amplified, and sequenced. Here, we have applied SVGs to deep-ocean samples (200–4000 m depth) from global Malaspina and MEDIMAX expeditions, demonstrating the feasibility of this method in deep-ocean samples. A total of 1328 virus-like particles were sorted from the North Atlantic Ocean, the deep Mediterranean Sea, and the Pacific Ocean oxygen minimum zone (OMZ). For this proof of concept, sixty single viruses were selected at random for sequencing. Genome annotation identified 27 of these genomes as bona fide viruses, and detected three auxiliary metabolic genes involved in nucleotide biosynthesis and sugar metabolism. Massive protein profile analysis confirmed that these viruses represented novel viral groups not present in databases. Although they were not previously assembled by viromics, global fragment recruitment analysis showed a conserved profile of relative abundance of these viruses in all analyzed samples spanning different oceans. Altogether, these results reveal the feasibility in using SVGs in this vast environment to unveil the genomes of relevant viruses.

Reduced bacterial mortality and enhanced viral productivity during sinking in the ocean

Particle sinking is an important process in the ocean, influencing the biogeochemical cycle and driving the long-term preservation of carbon into the deep sea via the biological pump. However, as an important component of marine ecosystems, the role of viruses during sinking is still poorly understood. Therefore, we performed a series of transplantation experiments in the South China Sea to simulate environmental changes during sinking and investigate their effects on viral eco-dynamics and life strategy. Our study demonstrated increased viral production but decreased virus-mediated bacterial mortality after transplantation. A larger burst size and switch from the lysogenic to lytic strategy were shown to contribute to enhanced viral productivity. We provide experimental evidence that surface viral ecological characteristics changed dramatically after transplantation into deep-sea waters, indicating a potential importance of viruses during vertical sinking in the ocean. This effect probably provides positive feedback on the efficiency of the biological pump.

Deciphering the Virus Signal Within the Marine Dissolved Organic Matter Pool

Viruses are ubiquitously distributed in the marine environment, influencing microbial population dynamics and biogeochemical cycles on a large scale. Due to their small size, they fall into the oceanographic size-class definition of dissolved organic matter (DOM; <0.7 μm). The purpose of our study was to investigate if there is a detectable imprint of virus particles in natural DOM following standard sample preparation and molecular analysis routines using ultrahigh-resolution mass spectrometry (FT-ICR-MS). Therefore, we tested if a molecular signature deriving from virus particles can be detected in the DOM fingerprint of a bacterial culture upon prophage induction and of seawater containing the natural microbial community. Interestingly, the virus-mediated lysate of the infected bacterial culture differed from the cell material of a physically disrupted control culture in its molecular composition. Overall, a small subset of DOM compounds correlated significantly with virus abundances in the bacterial culture setup, accounting for <1% of the detected molecular formulae and <2% of the total signal intensity of the DOM dataset. These were phosphorus- and nitrogen-containing compounds and they were partially also detected in DOM samples from other studies that included high virus abundances. While some of these formulae matched with typical biomolecules that are constituents of viruses, others matched with bacterial cell wall components. Thus, the identified DOM molecular formulae were probably not solely derived from virus particles but were partially also derived from processes such as the virus-mediated bacterial cell lysis. Our results indicate that a virus-derived DOM signature is part of the natural DOM and barely detectable within the analytical window of ultrahigh-resolution mass spectrometry when a high natural background is present.

Heterogeneous viral contribution to dissolved organic matter processing in a long-term macrocosm experiment

Viruses saturate environments throughout the world and play key roles in microbial food webs, yet how viral activities affect dissolved organic matter (DOM) processing in natural environments remains elusive. We established a large-scale long-term macrocosm experiment to explore viral dynamics and their potential impacts on microbial mortality and DOM quantity and quality in starved and stratified ecosystems. High viral infection dynamics and the virus-induced cell lysis (6.23-64.68% d^-1) was found in the starved seawater macrocosm, which contributed to a significant transformation of microbial biomass into DOM (0.72-5.32 μg L^-1 d^-1). In the stratified macrocosm, a substantial amount of viral lysate DOM (2.43-17.87 μg L^-1 d^-1) was released into the upper riverine water, and viral lysis and DOM release (0.35-5.75 μg L^-1 d^-1) were reduced in the mixed water layer between riverine water and seawater. Viral lysis was stimulated at the bottom of stratified macrocosm, potentially fueled by the sinking of particulate organic carbon. Significant positive and negative associations between lytic viral production and different fluorescent DOM components were found in the starved and stratified macrocosm, indicating the potentially complex viral impacts on the production and utilization of DOM. Results also revealed the significant viral contribution to pools of both relatively higher molecular weight labile DOM and lower molecular weight recalcitrant DOM. Our study suggests that viruses have heterogeneous impact on the cycling and fate of DOM in aquatic environments.

Predicting the capsid architecture of phages from metagenomic data

Tailed phages are viruses that infect bacteria and are the most abundant biological entities on Earth. Their ecological, evolutionary, and biogeochemical roles in the planet stem from their genomic diversity. Known tailed phage genomes range from 10 to 735 kilobase pairs thanks to the size variability of the protective protein capsids that store them. However, the role of tailed phage capsids’ diversity in ecosystems is unclear. A fundamental gap is the difficulty of associating genomic information with viral capsids in the environment. To address this problem, here, we introduce a computational approach to predict the capsid architecture (T-number) of tailed phages using the sequence of a single gene—the major capsid protein. This approach relies on an allometric model that relates the genome length and capsid architecture of tailed phages. This allometric model was applied to isolated phage genomes to generate a library that associated major capsid proteins and putative capsid architectures. This library was used to train machine learning methods, and the most computationally scalable model investigated (random forest) was applied to human gut metagenomes. Compared to isolated phages, the analysis of gut data reveals a large abundance of mid-sized (T = 7) capsids, as expected, followed by a relatively large frequency of jumbo-like tailed phage capsids (T ≥ 25) and small capsids (T = 4) that have been under-sampled. We discussed how to increase the method’s accuracy and how to extend the approach to other viruses. The computational pipeline introduced here opens the doors to monitor the ongoing evolution and selection of viral capsids across ecosystems.

Viral Production in Seawater Filtered Through 0.2-μm Pore-Size Filters: A Hidden Biogeochemical Cycle in a Neglected Realm

Viral production is a key parameter for assessing virus-mediated biogeochemical cycles. One widely used method for the determination of viral production, called the virus reduction assay, reduces viral abundance, while maintaining bacterial abundance, using 0.2-μm pore-size filters. Viral production is estimated from the increase of viral abundance during incubation. We hypothesized that small-cell-sized bacterial communities can pass through 0.2-μm filters and drive viral production, representing a missing fraction of viral production that is missed by the virus reduction assay. Coastal seawater was filtered through 0.2-μm filters and diluted with virus-free seawater. Viral production in the <0.2-μm filtrate was estimated from changes in viral abundance determined through flow cytometry. We found that viruses were produced in the <0.2-μm communities, which were strongly enriched with low nucleic acid content bacteria. Estimated viral production in the <0.2-μm filtrates accounted for up to 43% of total viral production and 10% of dissolved organic carbon production mediated by viral lysis of bacterial cells. By not considering viral production in these <0.2-μm communities, the virus reduction assay may underestimate viral production. Virus–bacteria interactions in <0.2-μm communities may represent a significant and overlooked role of viruses in marine food webs and carbon fluxes.

The global ocean size spectrum from bacteria to whales

It has long been hypothesized that aquatic biomass is evenly distributed among logarithmic body mass size classes. Although this community structure has been observed regionally, mostly among plankton groups, its generality has never been formally tested across all marine life over the global ocean, nor have the impacts of humans on it been globally assessed. Here, we bring together data at the global scale to test the hypothesis from bacteria to whales. We find that biomass within most order of magnitude size classes is indeed remarkably constant, near 1 gigatonne (Gt) wet weight (10¹⁵ g), but bacteria and large marine mammals are markedly above and below this value, respectively. Furthermore, human impacts appear to have significantly truncated the upper one-third of the spectrum. This dramatic alteration to what is possibly life’s largest-scale regularity underscores the global extent of human activities.

Diversity and distribution of viruses inhabiting the deepest ocean on Earth

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

The landscape of lysogeny across microbial community density, diversity and energetics

Lysogens are common at high bacterial densities, an observation that contrasts with the prevailing view of lysogeny as a low-density refugium strategy. Here, we review the mechanisms regulating lysogeny in complex communities and show that the additive effects of coinfections, diversity and host energic status yield a bimodal distribution of lysogeny as a function of microbial densities. At high cell densities (above 10⁶  cells ml^-1 or g^-1 ) and low diversity, coinfections by two or more phages are frequent and excess energy availability stimulates inefficient metabolism. Both mechanisms favour phage integration and characterize the Piggyback-the-Winner dynamic. At low densities (below 10⁵  cells ml^-1 or g^-1 ), starvation represses lytic genes and extends the time window for lysogenic commitment, resulting in a higher frequency of coinfections that cause integration. This pattern follows the predictions of the refugium hypothesis. At intermediary densities (between 10⁵ and 10⁶  cells ml^-1 or g^-1 ), encounter rates and efficient energy metabolism favour lysis. This may involve Kill-the-Winner lytic dynamics and induction. Based on these three regimes, we propose a framework wherein phage integration occurs more frequently at both ends of the host density gradient, with distinct underlying molecular mechanisms (coinfections and host metabolism) dominating at each extreme.

Expanding standards in viromics: in silico evaluation of dsDNA viral genome identification, classification, and auxiliary metabolic gene curation

BACKGROUND

Viruses influence global patterns of microbial diversity and nutrient cycles. Though viral metagenomics (viromics), specifically targeting dsDNA viruses, has been critical for revealing viral roles across diverse ecosystems, its analyses differ in many ways from those used for microbes. To date, viromics benchmarking has covered read pre-processing, assembly, relative abundance, read mapping thresholds and diversity estimation, but other steps would benefit from benchmarking and standardization. Here we use in silico-generated datasets and an extensive literature survey to evaluate and highlight how dataset composition (i.e., viromes vs bulk metagenomes) and assembly fragmentation impact (i) viral contig identification tool, (ii) virus taxonomic classification, and (iii) identification and curation of auxiliary metabolic genes (AMGs).

RESULTS

The in silico benchmarking of five commonly used virus identification tools show that gene-content-based tools consistently performed well for long (≥3 kbp) contigs, while k-mer- and blast-based tools were uniquely able to detect viruses from short (≤3 kbp) contigs. Notably, however, the performance increase of k-mer- and blast-based tools for short contigs was obtained at the cost of increased false positives (sometimes up to ∼5% for virome and ∼75% bulk samples), particularly when eukaryotic or mobile genetic element sequences were included in the test datasets. For viral classification, variously sized genome fragments were assessed using gene-sharing network analytics to quantify drop-offs in taxonomic assignments, which revealed correct assignations ranging from ∼95% (whole genomes) down to ∼80% (3 kbp sized genome fragments). A similar trend was also observed for other viral classification tools such as VPF-class, ViPTree and VIRIDIC, suggesting that caution is warranted when classifying short genome fragments and not full genomes. Finally, we highlight how fragmented assemblies can lead to erroneous identification of AMGs and outline a best-practices workflow to curate candidate AMGs in viral genomes assembled from metagenomes.

CONCLUSION

Together, these benchmarking experiments and annotation guidelines should aid researchers seeking to best detect, classify, and characterize the myriad viruses 'hidden' in diverse sequence datasets.

Moderate Seasonal Dynamics Indicate an Important Role for Lysogeny in the Red Sea

Viruses are the most abundant microorganisms in marine environments and viral infections can be either lytic (virulent) or lysogenic (temperate phage) within the host cell. The aim of this study was to quantify viral dynamics (abundance and infection) in the coastal Red Sea, a narrow oligotrophic basin with high surface water temperatures (22–32 °C degrees), high salinity (37.5–41) and continuous high insolation, thus making it a stable and relatively unexplored environment. We quantified viral and environmental changes in the Red Sea (two years) and the occurrence of lysogenic bacteria (induced by mitomycin C) on the second year. Water temperatures ranged from 24.0 to 32.5 °C, and total viral and bacterial abundances ranged from 1.5 to 8.7 × 10⁶ viruses mL⁻¹ and 1.9 to 3.2 × 10⁵ bacteria mL⁻¹, respectively. On average, 12.24% ± 4.8 (SE) of the prophage bacteria could be induced by mitomycin C, with the highest percentage of 55.8% observed in January 2018 when bacterial abundances were low; whereas no induction was measurable in spring when bacterial abundances were highest. Thus, despite the fact that the Red Sea might be perceived as stable, warm and saline, relatively modest changes in seasonal conditions were associated with large swings in the prevalence of lysogeny.

Novel virocell metabolic potential revealed in agricultural soils by virus-enriched soil metagenome analysis

Viruses are now recognized as important players in microbial dynamics and biogeochemical cycles in the oceans. Yet, compared with aquatic ecosystems, virus discovery in terrestrial ecosystems has been challenging partly due to the inherent complexity of soils. To expand our understanding of soil viruses and their putative contributions to soil microbial processes, we analysed metagenomes of community-level virus-enriched suspensions by tangential flow filtration obtained from two French agricultural soils. We found viral sequences representing a total of 239 viral operational taxonomic units that corresponded to 29.5% of the mapping reads in the metagenomic datasets. The analysis of their genomic sequences revealed novel virocell metabolic potential with implications to virus-host interactions, carbon cycling, plant-beneficial functions in the rhizosphere, horizontal gene transfer and other relevant microbial strategies applied to survive in soils.

Virome composition in marine fish revealed by meta-transcriptomics

Revealing the determinants of virome composition is central to placing disease emergence in a broader evolutionary context. Fish are the most species-rich group of vertebrates and so provide an ideal model system to study the factors that shape virome compositions and their evolution. We characterized the viromes of nineteen wild-caught species of marine fish using total RNA sequencing (meta-transcriptomics) combined with analyses of sequence and protein structural homology to identify divergent viruses that often evade characterization. From this, we identified twenty-five new vertebrate-associated viruses and a further twenty-two viruses likely associated with fish diet or their microbiomes. The vertebrate-associated viruses identified here included the first fish virus in the Matonaviridae (single-strand, negative-sense RNA virus). Other viruses fell within the Astroviridae, Picornaviridae, Arenaviridae, Reoviridae, Hepadnaviridae, Paramyxoviridae, Rhabdoviridae, Hantaviridae, Filoviridae, and Flaviviridae, and were sometimes phylogenetically distinct from known fish viruses. We also show how key metrics of virome composition-viral richness, abundance, and diversity-can be analysed along with host ecological and biological factors as a means to understand virus ecology. Accordingly, these data suggest that that the vertebrate-associated viromes of the fish sampled here are predominantly shaped by the phylogenetic history (i.e. taxonomic order) of their hosts, along with several biological factors including water temperature, habitat depth, community diversity and swimming behaviour. No such correlations were found for viruses associated with porifera, molluscs, arthropods, fungi, and algae, that are unlikely to replicate in fish hosts. Overall, these data indicate that fish harbour particularly large and complex viromes and the vast majority of fish viromes are undescribed.

Elevated Contribution of Low Nucleic Acid Prokaryotes and Viral Lysis to the Prokaryotic Community Along the Nutrient Gradient From an Estuary to Open Ocean Transect

Prokaryotes represent the largest living biomass reservoir in aquatic environments and play a crucial role in the global ocean. However, the factors that shape the abundance and potential growth rate of the ecologically distinct prokaryotic subgroups [i.e., high nucleic acid (HNA) and low nucleic acid (LNA) cells] along varying trophic conditions in the ocean remain poorly understood. This study conducted a series of modified dilution experiments to investigate how the abundance and potential growth rate of HNA and LNA prokaryotes and their regulating factors (i.e., protozoan grazing and viral lysis) change along a cross-shore nutrient gradient in the northern South China Sea. The results showed that the abundance of both HNA and LNA cells was significantly positively correlated with the abundance of heterotrophic nanoflagellates and viruses, whereas only HNA abundance exhibited a significant positive correlation with nutrient level. With a decreasing nutrient concentration, the potential growth rate of the HNA subgroup declined significantly, while that of the LNA subgroup was significantly enhanced, leading to an elevated relative potential growth rate of the LNA to HNA subgroup under decreasing nutrient levels. Furthermore, our data revealed different regulatory roles of protozoan grazing and viral lysis on the HNA and LNA subgroups, with HNA suffering higher mortality pressure from grazing than from lysis in contrast to LNA, which experienced equivalent pressures. As the nutrient levels declined, the relative contribution of lysis to the mortality of the HNA subgroup increased significantly, in contrast to the insignificant change in that of the LNA subgroup. Our results indicated the elevated role of LNA cells in the prokaryotic community and the enhanced viral lysis pressure on the total prokaryotes under oligotrophic conditions. This implies a weakened efficiency of carbon cycling within the microbial loop and enhanced viral lysis to shunt more carbon and energy flow in the future ocean, in which oligotrophication will be strengthened due to global warming.

Weekly variations of viruses and heterotrophic nanoflagellates and their potential impact on bacterioplankton in shallow waters of the central Red Sea

Bacterioplankton play a pivotal role in marine ecosystems. However, their temporal dynamics and underlying control mechanisms are poorly understood in tropical regions such as the Red Sea. Here, we assessed the impact of bottom-up (resource availability) and top-down (viruses and heterotrophic nanoflagellates) controls on bacterioplankton abundances by weekly sampling a coastal central Red Sea site in 2017. We monitored microbial abundances by flow cytometry together with a set of environmental variables including temperature, salinity, dissolved organic and inorganic nutrients and chlorophyll a. We distinguished five groups of heterotrophic bacteria depending on their physiological properties relative nucleic acid content, membrane integrity and cell-specific respiratory activity, two groups of Synechococcus cyanobacteria and three groups of viruses. Viruses controlled heterotrophic bacteria for most of the year, as supported by a negative correlation between their respective abundances and a positive one between bacterial mortality rates and mean viral abundances. On the contrary, heterotrophic nanoflagellates abundance covaried with that of heterotrophic bacteria. Heterotrophic nanoflagellates showed preference for larger bacteria from both the high and low nucleic acid content groups. Our results demonstrate that top-down control is fundamental in keeping heterotrophic bacterioplankton abundances low (< 5 × 10 ⁵ cells mL⁻¹) in Red Sea coastal waters.

Significance of Viral Activity for Regulating Heterotrophic Prokaryote Community Dynamics along a Meridional Gradient of Stratification in the Northeast Atlantic Ocean

How microbial populations interact influences the availability and flux of organic carbon in the ocean. Understanding how these interactions vary over broad spatial scales is therefore a fundamental aim of microbial oceanography. In this study, we assessed variations in the abundances, production, virus and grazing induced mortality of heterotrophic prokaryotes during summer along a meridional gradient in stratification in the North Atlantic Ocean. Heterotrophic prokaryote abundance and activity varied with phytoplankton biomass, while the relative distribution of prokaryotic subpopulations (ratio of high nucleic acid fluorescent (HNA) and low nucleic acid fluorescent (LNA) cells) was significantly correlated to phytoplankton mortality mode (i.e., viral lysis to grazing rate ratio). Virus-mediate morality was the primary loss process regulating the heterotrophic prokaryotic communities (average 55% of the total mortality), which may be attributed to the strong top-down regulation of the bacterivorous protozoans. Host availability, encounter rate, and HNA:LNA were important factors regulating viral dynamics. Conversely, the abundance and activity of bacterivorous protozoans were largely regulated by temperature and turbulence. The ratio of total microbial mediated mortality to total available prokaryote carbon reveals that over the latitudinal gradient the heterotrophic prokaryote community gradually moved from a near steady state system regulated by high turnover in subtropical region to net heterotrophic production in the temperate region.

Seasonal dynamics in taxonomy and function within bacterial and viral metagenomic assemblages recovered from a freshwater agricultural pond

BACKGROUND

Ponds are important freshwater habitats that support both human and environmental activities. However, relative to their larger counterparts (e.g. rivers, lakes), ponds are understudied, especially with regard to their microbial communities. Our study aimed to fill this knowledge gap by using culture-independent, high-throughput sequencing to assess the dynamics, taxonomy, and functionality of bacterial and viral communities in a freshwater agricultural pond.

RESULTS

Water samples (n = 14) were collected from a Mid-Atlantic agricultural pond between June 2017 and May 2018 and filtered sequentially through 1 and 0.2 μm filter membranes. Total DNA was then extracted from each filter, pooled, and subjected to 16S rRNA gene and shotgun sequencing on the Illumina HiSeq 2500 platform. Additionally, on eight occasions water filtrates were processed for viral metagenomes (viromes) using chemical concentration and then shotgun sequenced. A ubiquitous freshwater phylum, Proteobacteria was abundant at all sampling dates throughout the year. However, environmental characteristics appeared to drive the structure of the community. For instance, the abundance of Cyanobacteria (e.g. Nostoc) increased with rising water temperatures, while a storm event appeared to trigger an increase in overall bacterial diversity, as well as the relative abundance of Bacteroidetes. This event was also associated with an increase in the number of antibiotic resistance genes. The viral fractions were dominated by dsDNA of the order Caudovirales, namely Siphoviridae and Myovirdae.

CONCLUSIONS

Overall, this study provides one of the largest datasets on pond water microbial ecology to date, revealing seasonal trends in the microbial taxonomic composition and functional potential.

Low Host Abundance and High Temperature Determine Switching from Lytic to Lysogenic Cycles in Planktonic Microbial Communities in a Tropical Sea (Red Sea)

The lytic and lysogenic life cycles of marine phages are influenced by environmental conditions such as solar radiation, temperature, and host abundance. Temperature can regulate phage infection, but its role is difficult to discern in oligotrophic waters where there is typically low host abundance and high temperatures. Here, we study the temporal variability of viral dynamics and the occurrence of lysogeny using mitomycin C in a eutrophic coastal lagoon in the oligotrophic Red Sea, which showed strong seasonality in terms of temperature (22.1–33.3 °C) and large phytoplankton blooms. Viral abundances ranged from 2.2 × 10⁶ to 1.5 × 10⁷ viruses mL⁻¹ and were closely related to chlorophyll a (chl a) concentration. Observed high virus-to-bacterium ratio (VBR) (4–79; 16 ± 4 (SE)) suggests that phages exerted a tight control of their hosts as indicated by the significant decrease in bacterial abundance with increasing virus concentration. Heterotrophic bacterial abundance also showed a significant decrease with increasing temperature. However, viral abundance was not related to temperature changes and the interaction of water temperature, suggesting an indirect effect of temperature on decreased host abundance, which was observed at the end of the summertime. From the estimated burst size (BS), we observed lysogeny (undetectable to 29.1%) at low percentages of 5.0% ± 1.2 (SE) in half of the incubations with mitomycin C, while it increased to 23.9% ± 2.8 (SE) when the host abundance decreased. The results suggest that lytic phages predominate, switching to a moderate proportion of temperate phages when the host abundance reduces.

Between a Rock and a Soft Place: The Role of Viruses in Lithification of Modern Microbial Mats

Stromatolites are geobiological systems formed by complex microbial communities, and fossilized stromatolites provide a record of some of the oldest life on Earth. Microbial mats are precursors of extant stromatolites; however, the mechanisms of transition from mat to stromatolite are controversial and are still not well understood. To fully recognize the profound impact that these ecosystems have had on the evolution of the biosphere requires an understanding of modern lithification mechanisms and how they relate to the geological record. We propose here viral mechanisms in carbonate precipitation, leading to stromatolite formation, whereby viruses directly or indirectly impact microbial metabolisms that govern the transition from microbial mat to stromatolite. Finding a tangible link between host-virus interactions and changes in biogeochemical processes will provide tools to interpret mineral biosignatures through geologic time, including those on Earth and beyond.

Prokaryotic Capability to Use Organic Substrates Across the Global Tropical and Subtropical Ocean

Prokaryotes play a fundamental role in decomposing organic matter in the ocean, but little is known about how microbial metabolic capabilities vary at the global ocean scale and what are the drivers causing this variation. We aimed at obtaining the first global exploration of the functional capabilities of prokaryotes in the ocean, with emphasis on the under-sampled meso- and bathypelagic layers. We explored the potential utilization of 95 carbon sources with Biolog GN2 plates^® in 441 prokaryotic communities sampled from surface to bathypelagic waters (down to 4,000 m) at 111 stations distributed across the tropical and subtropical Atlantic, Indian, and Pacific oceans. The resulting metabolic profiles were compared with biological and physico-chemical properties such as fluorescent dissolved organic matter (DOM) or temperature. The relative use of the individual substrates was remarkably consistent across oceanic regions and layers, and only the Equatorial Pacific Ocean showed a different metabolic structure. When grouping substrates by categories, we observed some vertical variations, such as an increased relative utilization of polymers in bathypelagic layers or a higher relative use of P-compounds or amino acids in the surface ocean. The increased relative use of polymers with depth, together with the increases in humic DOM, suggest that deep ocean communities have the capability to process complex DOM. Overall, the main identified driver of the metabolic structure of ocean prokaryotic communities was temperature. Our results represent the first global depiction of the potential use of a variety of carbon sources by prokaryotic communities across the tropical and the subtropical ocean and show that acetic acid clearly emerges as one of the most widely potentially used carbon sources in the ocean.

Picocyanobacteria Community and Cyanophage Infection Responses to Nutrient Enrichment in a Mesocosms Experiment in Oligotrophic Waters

Prochlorococcus and Synechococcus are pico-sized cyanobacteria that play a fundamental role in oceanic primary production, being particularly important in warm, nutrient-poor waters. Their potential response to nutrient enrichment is expected to be contrasting and to differ from larger phytoplankton species. Here, we used a metagenomic approach to characterize the responses to nutrient enrichment in the community of picocyanobacteria and to analyze the cyanophage response during a mesocosms experiment in the oligotrophic Red Sea. Natural picoplankton community was dominated by Synechococcus clade II, with marginal presence of Prochlorococcus (0.3% bacterial reads). Increased nutrient input triggered a fast Synechococcus bloom, with clade II being the dominant, with no response of Prochlorococcus growth. The largest bloom developed in the mesocosms receiving a single initial input of nutrients, instead of daily additions. The relative abundances of cyanophage sequences in cellular metagenomes increased during the experiment from 12.6% of total virus reads up to 40% in the treatment with the largest Synechococcus bloom. The subsequent collapse of the bloom pointed to a cyanophage infection on Synechococcus that reduced its competitive capacity, and was then followed by a diatom bloom. The cyanophage attack appears to have preferentially affected the most abundant Synechococcus clade II, increasing the evenness within the host population. Our results highlight the relevance of host-phage interactions on determining population dynamics and diversity of Synechococcus populations.

Viral control of biomass and diversity of bacterioplankton in the deep sea

Viral abundance in deep-sea environments is high. However, the biological, ecological and biogeochemical roles of viruses in the deep sea are under debate. In the present study, microcosm incubations of deep-sea bacterioplankton (2,000 m deep) with normal and reduced pressure of viral lysis were conducted in the western Pacific Ocean. We observed a negative effect of viruses on prokaryotic abundance, indicating the top-down control of bacterioplankton by virioplankton in the deep-sea. The decreased bacterial diversity and a different bacterial community structure with diluted viruses indicate that viruses are sustaining a diverse microbial community in deep-sea environments. Network analysis showed that relieving viral pressure decreased the complexity and clustering coefficients but increased the proportion of positive correlations for the potentially active bacterial community, which suggests that viruses impact deep-sea bacterioplankton interactions. Our study provides experimental evidences of the crucial role of viruses in microbial ecology and biogeochemistry in deep-sea ecosystems.

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

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