Quantification of diverse virus populations in the environment using the polony method

Hepcidin contributes to largemouth bass (Micropterus salmoides) against bacterial infections

The increasing emergence and dissemination of bacterial pathogens in largemouth bass culture accelerate the desire for new treatment measures. Antimicrobial peptides as the host's antimicrobial source dominate the preferred molecules for discovering antibacterial agents. Here, the potential of Hepcidin-1 from largemouth bass (Micropterus salmoides) (MsHep-1) against bacterial infection is demonstrated. MsHep-1 not only improved the survival rate in infection experiments involving Nocardia seriolae (12 %) and Aeromonas hydrophila (18 %) but also coped with iron overload conditions in vivo. Moreover, the antibacterial activity of MsHep-1 in vitro was identified against both gram-negative and gram-positive bacteria. Mechanistic studies show MsHep-1 leads to bacterial death by changing the bacterial membrane potential and disrupting the bacterial membrane structure. These findings demonstrate that MsHep-1 may play an important role in the host response to bacterial infection. It provides promising strategies in the application of immunosuppression prevention and control in fish. AMPs may be a promising and available reservoir for treating the current bacterial diseases.

Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. We find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

Biological interactions with Prochlorococcus: implications for the marine carbon cycle

The unicellular picocyanobacterium Prochlorococcus is the most abundant photoautotroph and contributes substantially to global CO2 fixation. In the vast euphotic zones of the open ocean, Prochlorococcus converts CO2 into organic compounds and supports diverse organisms, forming an intricate network of interactions that regulate the magnitude of carbon cycling and storage in the ocean. An understanding of the biological interactions with Prochlorococcus is critical for accurately estimating the contributions of Prochlorococcus and interacting organisms to the marine carbon cycle. This review synthesizes the primary production contributed by Prochlorococcus in the global ocean. We outline recent progress on the interactions of Prochlorococcus with heterotrophic bacteria, phages, and grazers that multifacetedly determine Prochlorococcus carbon production and fate. We discuss that climate change might affect the biological interactions with Prochlorococcus and thus the marine carbon cycle.

Increased cyanophage infection at the bottom of the euphotic zone, especially in the fall

Picocyanobacteria contribute greatly to offshore primary production with cells extending through the deep euphotic zone. Literature indicates high viral infection of cyanobacteria in ocean transition zones. We postulate that the bottom of the euphotic zone is a transition zone, where communities transition from phototrophic to aphotic processes. We use single-copy core genes to examine cyanophage to cyanobacteria ratios in cellular metagenomes in the subtropical North Atlantic and Pacific. Cyanophage to cyanobacteria terL/rpoB ratios generally increase to >10 in the deep euphotic zone. As light levels decrease in the fall, Prochlorococcus in the deep euphotic zone experience reduced light levels. We find clear differences between spring (Geotraces GA02) and fall (GA03) in the North Atlantic, with terL/rpoB ratios increasing to >40 in the fall. When examining 23 months of the North Pacific Hawaii Ocean Timeseries, the depth of elevated cyanophage to cyanobacteria ratios in cellular metagenomes negatively correlated with surface photosynthetic radiation (PAR), particularly with the change in PAR, which reflected the season. In fall, all picocyanobacteria ecotypes were found at depths enriched with viruses, while in summer, only low light ecotypes were affected. Thus, we find high cyanophage infection both in the deep euphotic zone and during seasonal transitions.

Daily turnover of active giant virus infection during algal blooms revealed by single-cell transcriptomics

Giant viruses infect many unicellular eukaryotes, including algae that form massive oceanic blooms. Despite the major impact of viruses on the marine ecosystem, the ability to quantify and assess active viral infection in nature remains a major challenge. We applied single-cell RNA sequencing, to profile virus and host transcriptomes of 12,000 single algal cells from a coccolithophore bloom. Viral infection was detected already at early exponential bloom phase, negatively correlating with the bloom intensity. A consistent percent of infected coccolithophores displayed the early phase of viral replication for several consecutive days, indicating a daily turnover and continuous virocell-associated metabolite production, potentially affecting the surrounding microbiome. Linking single-cell infection state to host physiology revealed that infected cells remained calcified even in the late infection stage. These findings stress the importance of studying host-virus dynamics in natural populations, at single-cell resolution, to better understand virus life cycle and its impact on microbial food webs.

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.

Development and validation of a real-time PCR assay protocol for the specific detection and quantification of pelagiphages in seawater samples

Earth is inhabited by numerous adaptations of cellular forms shaped by the persistent scrutiny of natural selection. Thus, as natural selection has fixed beneficial adaptations of functional traits, cellular life has conquered almost all environmental niches on our planet. However, cellular life succumbs in number and genetic diversity to viruses. Among all viruses, phages are highly prevalent in diverse environments, and due to their vast genetic diversity and abundance, their relevant role as significant players in several ecological processes is now fully recognized. Pelagiphages, bacteriophages infecting bacteria of the SAR11 clade, are the most abundant viruses in the oceans. However, the ecological contribution of pelagiphages on populations of Pelagibacterales remains largely underestimated. An essential aspect of estimating the impact of bacteriophages is their absolute and precise quantification, which provides relevant information about the host-virus interactions and the structure of viral assemblages. Consequently, due to its abundance and claimed influence in the biogeochemical cycling of elements, the accurate quantification of pelagiphages results in an essential task. This study describes the development and validation of a sensitive, specific, accurate and reproducible qPCR platform targeting pelagiphages. Moreover, this method allowed the detection and quantification of pelagiphages in the Gulf of California for the first time.

Abundant and cosmopolitan lineage of cyanopodoviruses lacking a DNA polymerase gene

Cyanopodoviruses affect the mortality and population dynamics of the unicellular picocyanobacteria Prochlorococcus and Synechococcus, the dominant primary producers in the oceans. Known cyanopodoviruses all contain the DNA polymerase gene (DNA pol) that is important for phage DNA replication and widely used in field quantification and diversity studies. However, we isolated 18 cyanopodoviruses without identifiable DNA pol. They form a new MPP-C clade that was separated from the existing MPP-A, MPP-B, and P-RSP2 clades. The MPP-C phages have the smallest genomes (37.3-37.9 kb) among sequenced cyanophages, and show longer latent periods than the MPP-B phages. Metagenomic reads of both clades are highly abundant in surface waters, but the MPP-C phages show higher relative abundance in surface waters than in deeper waters, while MPP-B phages have higher relative abundance in deeper waters. Our study reveals that cyanophages with distinct genomic contents and infection kinetics can exhibit different depth profiles in the oceans.

Associations between picocyanobacterial ecotypes and cyanophage host genes across ocean basins and depth

Background

Cyanophages, viruses that infect cyanobacteria, are globally abundant in the ocean's euphotic zone and are a potentially important cause of mortality for marine picocyanobacteria. Viral host genes are thought to increase viral fitness by either increasing numbers of genes for synthesizing nucleotides for virus replication, or by mitigating direct stresses imposed by the environment. The encoding of host genes in viral genomes through horizontal gene transfer is a form of evolution that links viruses, hosts, and the environment. We previously examined depth profiles of the proportion of cyanophage containing various host genes in the Eastern Tropical North Pacific Oxygen Deficient Zone (ODZ) and at the subtropical North Atlantic (BATS). However, cyanophage host genes have not been previously examined in environmental depth profiles across the oceans.

Methodology

We examined geographical and depth distributions of picocyanobacterial ecotypes, cyanophage, and their viral-host genes across ocean basins including the North Atlantic, Mediterranean Sea, North Pacific, South Pacific, and Eastern Tropical North and South Pacific ODZs using phylogenetic metagenomic read placement. We determined the proportion of myo and podo-cyanophage containing a range of host genes by comparing to cyanophage single copy core gene terminase (terL). With this large dataset (22 stations), network analysis identified statistical links between 12 of the 14 cyanophage host genes examined here with their picocyanobacteria host ecotypes.

Results

Picyanobacterial ecotypes, and the composition and proportion of cyanophage host genes, shifted dramatically and predictably with depth. For most of the cyanophage host genes examined here, we found that the composition of host ecotypes predicted the proportion of viral host genes harbored by the cyanophage community. Terminase is too conserved to illuminate the myo-cyanophage community structure. Cyanophage cobS was present in almost all myo-cyanophage and did not vary in proportion with depth. We used the composition of cobS phylotypes to track changes in myo-cyanophage composition.

Conclusions

Picocyanobacteria ecotypes shift with changes in light, temperature, and oxygen and many common cyanophage host genes shift concomitantly. However, cyanophage phosphate transporter gene pstS appeared to instead vary with ocean basin and was most abundant in low phosphate regions. Abundances of cyanophage host genes related to nutrient acquisition may diverge from host ecotype constraints as the same host can live in varying nutrient concentrations. Myo-cyanophage community in the anoxic ODZ had reduced diversity. By comparison to the oxic ocean, we can see which cyanophage host genes are especially abundant (nirA, nirC, and purS) or not abundant (myo psbA) in ODZs, highlighting both the stability of conditions in the ODZ and the importance of nitrite as an N source to ODZ endemic LLV Prochlorococcus.

Viral infection in the ocean-A journey across scales

Viruses are the most abundant biological entity in the ocean and infect a wide range of microbial life across bacteria, archaea, and eukaryotes. In this essay, we take a journey across several orders of magnitude in the scales of biological organization, time, and space of host-virus interactions in the ocean, aiming to shed light on their ecological relevance. We start from viruses infecting microbial host cells by delivering their genetic material in seconds across nanometer-size membranes, which highjack their host's metabolism in a few minutes to hours, leading to a profound transcriptomic and metabolic rewiring. The outcome of lytic infection leads to a release of virions and signaling molecules that can reach neighboring cells a few millimeters away, resulting in a population whose heterogeneous infection level impacts the surrounding community for days. These population dynamics can leave unique metabolic and biogeochemical fingerprints across scales of kilometers and over several decades. One of the biggest challenges in marine microbiology is to assess the impact of viruses across these scales, from the single cell to the ecosystem level. Here, we argue that the advent of new methodologies and conceptual frameworks represents an exciting time to pursue these efforts and propose a set of important challenges for the field. A better understanding of host-virus interactions across scales will inform models of global ocean ecosystem function in different climate change scenarios.

Widespread yet persistent low abundance of TIM5-like cyanophages in the oceans

Ocean ecosystems are inhabited by a diverse set of viruses that impact microbial mortality and evolution. However, the distribution and abundances of specific viral lineages, particularly those from the large bank of rare viruses, remains largely unknown. Here, we assessed the diversity and abundance of the TIM5-like cyanophages. The sequencing of three new TIM5-like cyanophage genomes and environmental amplicons of a signature gene from the Red Sea revealed highly conserved gene content and sequence similarity. We adapted the polony method, a solid-phase polymerase chain reaction assay, to quantify TIM5-like cyanophages during three 2000 km expeditions in the Pacific Ocean and four annual cycles in the Red Sea. TIM5-like cyanophages were widespread, detected at all latitudes and seasons surveyed throughout the photic zone. Yet they were generally rare, ranging between <100 and 4000 viruses·ml^-1 . Occasional peaks in abundance of 10- to 100-fold were observed, reaching 71,000 viruses·ml^-1 . These peaks were ephemeral and seasonally variable in the Red Sea. Infection levels, quantified during one such peak, were very low. These characteristics of low diversity and abundance, as well as variable outbreaks, distinguishes the TIM5-like lineage from other major cyanophage lineages and illuminates that rare virus lineages can be persistent and widespread in the oceans.

Cyanophages from a less virulent clade dominate over their sister clade in global oceans

Environmental virus communities are highly diverse. However, the infection physiology underlying the evolution of diverse phage lineages and their ecological consequences are largely unknown. T7-like cyanophages are abundant in nature and infect the marine unicellular cyanobacteria, Synechococcus and Prochlorococcus, important primary producers in the oceans. Viruses belonging to this genus are divided into two distinct phylogenetic clades: clade A and clade B. These viruses have narrow host-ranges with clade A phages primarily infecting Synechococcus genotypes, while clade B phages are more diverse and can infect either Synechococcus or Prochlorococcus genotypes. Here we investigated infection properties (life history traits) and environmental abundances of these two clades of T7-like cyanophages. We show that clade A cyanophages have more rapid infection dynamics, larger burst sizes and greater virulence than clade B cyanophages. However, clade B cyanophages were at least 10-fold more abundant in all seasons, and infected more cyanobacteria, than clade A cyanophages in the Red Sea. Models predicted that steady-state cyanophage abundances, infection frequency, and virus-induced mortality, peak at intermediate virulence values. Our findings indicate that differences in infection properties are reflected in virus phylogeny at the clade level. They further indicate that infection properties, together with differences in subclade diversity and host repertoire, have important ecological consequences with the less aggressive, more diverse virus clade having greater ecological impacts.

Microbial Sources of Exocellular DNA in the Ocean

Exocellular DNA is operationally defined as the fraction of the total DNA pool that passes through a membrane filter (0.1 μm). It is composed of DNA-containing vesicles, viruses, and free DNA and is ubiquitous in all aquatic systems, although the sources, sinks, and ecological consequences are largely unknown. Using a method that provides separation of these three fractions, we compared open ocean depth profiles of DNA associated with each fraction. Pelagibacter-like DNA dominated the vesicle fractions for all samples examined over a depth range of 75 to 500 m. Viral DNA consisted predominantly of myovirus-like and podovirus-like DNA and contained the highest proportion of unannotated sequences. Euphotic zone free DNA (75 to 125 m) contained primarily bacterial and viral sequences, with bacteria dominating samples from the mesopelagic zone (500 to 1,000 m). A high proportion of mesopelagic zone free DNA sequences appeared to originate from surface waters, including a large amount of DNA contributed by high-light Prochlorococcus ecotypes. Throughout the water column, but especially in the mesopelagic zone, the composition of free DNA sequences was not always reflective of cooccurring microbial communities that inhabit the same sampling depth. These results reveal the composition of free DNA in different regions of the water column (euphotic and mesopelagic zones), with implications for dissolved organic matter cycling and export (by way of sinking particles and/or migratory zooplankton) as a delivery mechanism. IMPORTANCE With advances in metagenomic sequencing, the microbial composition of diverse environmental systems has been investigated, providing new perspectives on potential ecological dynamics and dimensions for experimental investigations. Here, we characterized exocellular free DNA via metagenomics, using a newly developed method that separates free DNA from cells, viruses, and vesicles, and facilitated the independent characterization of each fraction. The fate of this free DNA has both ecological consequences as a nutrient (N and P) source and potential evolutionary consequences as a source of genetic transformation. Here, we document different microbial sources of free DNA at the surface (0 to 200 m) versus depths of 250 to 1,000 m, suggesting that distinct free DNA production mechanisms may be present throughout the oligotrophic water column. Examining microbial processes through the lens of exocellular DNA provides insights into the production of labile dissolved organic matter (i.e., free DNA) at the surface (likely by viral lysis) and processes that influence the fate of sinking, surface-derived organic matter.

Viruses affect picocyanobacterial abundance and biogeography in the North Pacific Ocean

The photosynthetic picocyanobacteria Prochlorococcus and Synechococcus are models for dissecting how ecological niches are defined by environmental conditions, but how interactions with bacteriophages affect picocyanobacterial biogeography in open ocean biomes has rarely been assessed. We applied single-virus and single-cell infection approaches to quantify cyanophage abundance and infected picocyanobacteria in 87 surface water samples from five transects that traversed approximately 2,200 km in the North Pacific Ocean on three cruises, with a duration of 2-4 weeks, between 2015 and 2017. We detected a 550-km-wide hotspot of cyanophages and virus-infected picocyanobacteria in the transition zone between the North Pacific Subtropical and Subpolar gyres that was present in each transect. Notably, the hotspot occurred at a consistent temperature and displayed distinct cyanophage-lineage composition on all transects. On two of these transects, the levels of infection in the hotspot were estimated to be sufficient to substantially limit the geographical range of Prochlorococcus. Coincident with the detection of high levels of virally infected picocyanobacteria, we measured an increase of 10-100-fold in the Synechococcus populations in samples that are usually dominated by Prochlorococcus. We developed a multiple regression model of cyanophages, temperature and chlorophyll concentrations that inferred that the hotspot extended across the North Pacific Ocean, creating a biological boundary between gyres, with the potential to release organic matter comparable to that of the sevenfold-larger North Pacific Subtropical Gyre. Our results highlight the probable impact of viruses on large-scale phytoplankton biogeography and biogeochemistry in distinct regions of the oceans.

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

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

Isolation and characterization of a novel Lambda-like phage infecting the bloom-forming cyanobacteria Cylindrospermopsis raciborskii

Cylindrospermopsis raciborskii is a central bloom‐forming cyanobacteria. However, despite its ecological significance, little is known of its interactions with the phages that infect it. Currently, only a single sequenced genome of a Cylindrospermopsis‐infecting phage is publicly available. Here we describe the isolation and characterization of Cr‐LKS3, a second phage infecting Cylindrospermopsis. Cr‐LKS3 is a siphovirus with a higher genome similarity to prophages within heterotrophic bacteria genomes than to any other cyanophage/cyano‐prophage, suggesting that it represents a novel cyanophage group. The function, order and orientation of the 72 genes in the Cr‐LKS3 genome are highly similar to those of Escherichia virus Lambda (hereafter Lambda), despite the very low sequence similarity between these phages, showing high evolutionary convergence despite the substantial difference in host characteristics. Similarly to Lambda, the genome of Cr‐LKS3 contains various genes that are known to be central to lysogeny, suggesting it can enter a lysogenic cycle. Cr‐LKS3 has a unique ability to infect a host with a dramatically different GC content, without carrying any tRNA genes to compensate for this difference. This ability, together with its potential lysogenic lifestyle shed light on the complex interactions between C. raciborskii and its phages.

Adaptation of the polony technique to quantify Gokushovirinae, a diverse group of single-stranded DNA phage

Advances in metagenomics have revealed the ubiquity of single-stranded DNA (ssDNA) phage belonging to the subfamily Gokushovirinae in the oceans; however, the abundance and ecological roles of this group are unknown. Here, we quantify gokushoviruses through adaptation of the polony method, in which viral template DNA is immobilized in a gel, amplified by PCR, and subsequently detected by hybridization. Primers and probes for this assay were designed based on PCR amplicon diversity of gokushovirus major capsid protein gene sequences from a depth profile in the Gulf of Aqaba, Red Sea sampled in September 2015. At ≥95% identity, these 87 gokushovirus sequences formed 14 discrete clusters with the largest clades showing distinct depth distributions. The application of the polony method enabled the first quantification of gokushoviruses in any environment. The gokushoviruses were most abundant in the upper 40 m of the stratified water column, with a subsurface peak in abundance of 1.26 × 10⁵ viruses ml^-1 . These findings suggest that discrete gokushovirus genotypes infect bacterial hosts that differentially partition in the water column. Since the designed primers and probe are conserved across marine ecosystems, this polony method can be applied broadly for the quantification of gokushoviruses throughout the global oceans.

Visualizing active viral infection reveals diverse cell fates in synchronized algal bloom demise

Marine viruses are the most abundant biological entity in the ocean and are considered as major evolutionary drivers of microbial life [C. A. Suttle, Nat. Rev. Microbiol. 5, 801-812 (2007)]. Yet, we lack quantitative approaches to assess their impact on the marine ecosystem. Here, we provide quantification of active viral infection in the bloom forming single-celled phytoplankton Emiliania huxleyi infected by the large virus EhV, using high-throughput single-molecule messenger RNA in situ hybridization (smFISH) of both virus and host transcripts. In natural samples, viral infection reached only 25% of the population despite synchronized bloom demise exposing the coexistence of infected and noninfected subpopulations. We prove that photosynthetically active cells chronically release viral particles through nonlytic infection and that viral-induced cell lysis can occur without viral release, thus challenging major assumptions regarding the life cycle of giant viruses. We could also assess active infection in cell aggregates linking viral infection and carbon export to the deep ocean [C. P. Laber et al., Nat. Microbiol. 3, 537-547 (2018)] and suggest a potential host defense strategy by enrichment of infected cells in sinking aggregates. Our approach can be applied to diverse marine microbial systems, opening a mechanistic dimension to the study of biotic interactions in the ocean.

Revisiting the rules of life for viruses of microorganisms

Viruses that infect microbial hosts have traditionally been studied in laboratory settings with a focus on either obligate lysis or persistent lysogeny. In the environment, these infection archetypes are part of a continuum that spans antagonistic to beneficial modes. In this Review, we advance a framework to accommodate the context-dependent nature of virus-microorganism interactions in ecological communities by synthesizing knowledge from decades of virology research, eco-evolutionary theory and recent technological advances. We discuss that nuanced outcomes, rather than the extremes of the continuum, are particularly likely in natural communities given variability in abiotic factors, the availability of suboptimal hosts and the relevance of multitrophic partnerships. We revisit the 'rules of life' in terms of how long-term infections shape the fate of viruses and microbial cells, populations and ecosystems.

Integrating Viral Metagenomics into an Ecological Framework

Viral metagenomics has expanded our knowledge of the ecology of uncultured viruses, within both environmental (e.g., terrestrial and aquatic) and host-associated (e.g., plants and animals, including humans) contexts. Here, we emphasize the implementation of an ecological framework in viral metagenomic studies to address questions in virology rarely considered ecological, which can change our perception of viruses and how they interact with their surroundings. An ecological framework explicitly considers diverse variants of viruses in populations that make up communities of interacting viruses, with ecosystem-level effects. It provides a structure for the study of the diversity, distributions, dynamics, and interactions of viruses with one another, hosts, and the ecosystem, including interactions with abiotic factors. An ecological framework in viral metagenomics stands poised to broadly expand our knowledge in basic and applied virology. We highlight specific fundamental research needs to capitalize on its potential and advance the field. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

New technologies for developing phage-based tools to manipulate the human microbiome

Gut bacteria play an essential role in the human body by regulating multiple functions, producing essential metabolites, protecting against pathogen invasion, and much more. Conversely, changes in their community structure are linked to several gastrointestinal (GI) and non-GI conditions. Fortunately, these bacteria are amenable to external perturbations, but we need specific tools for their safe manipulation as nonspecific changes can cause unpredicted long-term consequences. Here, we mainly discuss recent advances in cultivation-independent technologies and argue their relevance to different key steps, that is, identifying the modulation targets and developing phage-based tools to precisely modulate gut bacteria and restore a sustainable microbiome in humans. We finally suggest multiple modulating strategies for different dysbiosis-associated diseases.

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.

A rigorous assessment and comparison of enumeration methods for environmental viruses

Determining exact viral titers in a given sample is essential for many environmental and clinical applications, e.g., for studying viral ecology or application of bacteriophages for food safety. However, virus quantification is not a simple task, especially for complex environmental samples. While clonal viral isolates can be quantified with relative high accuracy using virus-specific methods, i.e., plaque assay or quantitative real-time PCR, these methods are not valid for complex and diverse environmental samples. Moreover, it is not yet known how precisely laser-based methods, i.e., epifluorescence microscopy, flow cytometry, and nanoparticle tracking analysis, quantify environmental viruses. In the present study, we compared five state-of-the-art viral quantification methods by enumerating four model viral isolates of different genome and size characteristics as well as four different environmental water samples. Although Nanoparticle tracking analysis combined with gentle staining at 30 °C could be confirmed by this study to be a reliable quantification technique for tested environmental samples, environmental samples still lack an universally applicable and accurate quantification method. Special attention has to be put on optimal sample concentrations as well as optimized sample preparations, which are specific for each method. As our results show the inefficiency when enumerating small, or single-stranded DNA or RNA viruses, the global population of viruses is presumably higher than expected.

Cyanophage host-derived genes reflect contrasting selective pressures with depth in the oxic and anoxic water column of the Eastern Tropical North Pacific

Cyanophages encode host-derived genes that may increase their fitness. We examined the relative abundance of 18 host-derived cyanophages genes in metagenomes and viromes along depth profiles from the Eastern Tropical North Pacific Oxygen Deficient Zone (ETNP ODZ) where Prochlorococcus dominates a secondary chlorophyll maximum within the ODZ. Cyanophages at the oxic primary chlorophyll maximum encoded genes related to light and phosphate stress (psbA, psbD and pstS in T4-like and psbA in T7-like), but the proportion of cyanophage with these genes decreased with depth. The proportion of cyanophage with purine biosynthesis genes increased with depth in T4-like, but not T7-like cyanophages. No additional host-derived genes were found in deep T7-like cyanophages, suggesting that T4-like and T7-like cyanophages have different host-derived gene acquisition strategies, possibly linked to their different genome packaging mechanisms. In contrast to the ETNP, in the oxic North Atlantic T4-like cyanophages encoded psbA and pstS throughout the euphotic zone. Differences in pstS between the ETNP and the North Atlantic stations were consistent with differences in phosphate concentrations in those regimes. We suggest that the low proportion of cyanophage with psbA within the ODZ reflects the stably stratified low-light conditions occupied by their hosts, a Prochlorococcus ecotype endemic to ODZs.

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.

Epigenetic biosensors for bacteriophage detection and phage receptor discrimination

Environmental monitoring of bacteria using phage-based biosensors has been widely developed for many different species. However, there are only a few available methods to detect specific bacteriophages in raw environmental samples. In this work, we developed a simple and efficient assay to rapidly monitor the phage content of a given sample. The assay is based on the bistable expression of the Salmonella enterica opvAB operon. Under regular growth conditions, opvAB is only expressed by a small fraction of the bacterial subpopulation. In the OpvABON subpopulation, synthesis of the OpvA and OpvB products shortens the O-antigen and confers resistance to phages that use LPS as a receptor. As a consequence, the OpvABON subpopulation is selected in the presence of such phages. Using an opvAB::gfp fusion, we could monitor LPS-binding phages in various media, including raw water samples. To enlarge our phage-biosensor panoply, we also developed biosensors able to detect LPS, as well as protein-binding coliphages. Moreover, the combination of these tools allowed to identify the bacterial receptor triggering phage infection. The epigenetic opvAB::gfp biosensor thus comes in different flavours to detect a wide range of bacteriophages and identify the type of receptor they recognize.

Visualization of Viral Infection Dynamics in a Unicellular Eukaryote and Quantification of Viral Production Using Virus Fluorescence in situ Hybridization

One of the major challenges in viral ecology is to assess the impact of viruses in controlling the abundance of specific hosts in the environment. To this end, techniques that enable the detection and quantification of virus-host interactions at the single-cell level are essential. With this goal in mind, we implemented virus fluorescence in situ hybridization (VirusFISH) using as a model the marine picoeukaryote Ostreococcus tauri and its virus Ostreococcus tauri virus 5 (OtV5). VirusFISH allowed the visualization and quantification of the proportion of infected cells during an infection cycle in experimental conditions. We were also able to quantify the abundance of free viruses released during cell lysis, discriminating OtV5 from other mid-level fluorescence phages in our non-axenic infected culture that were not easily distinguishable with flow cytometry. Our results showed that although the major lysis of the culture occurred between 24 and 48 h after OtV5 inoculation, some new viruses were already produced between 8 and 24 h. With this work, we demonstrate that VirusFISH is a promising technique to study specific virus-host interactions in non-axenic cultures and establish a framework for its application in complex natural communities.

Quantification of T4-Like and T7-Like Cyanophages Using the Polony Method Show They Are Significant Members of the Virioplankton in the North Pacific Subtropical Gyre

The North Pacific Subtropical Gyre (NPSG) is one of the largest biomes on Earth, with the cyanobacterium Prochlorococcus being the most abundant primary producer year-round. Viruses that infect cyanobacteria (cyanophages) influence cyanobacterial mortality, diversity and evolution. Two major cyanophage families are the T4-like cyanomyoviruses and T7-like cyanopodoviruses, yet their abundances and distribution patterns remain unknown due to difficulty in quantifying their populations. To address this limitation, we previously adapted the polony method (for PCR colony) to quantify T7-like cyanophages and applied it to spring populations in the Red Sea. Here, we further adapted the method for the quantification of T4-like cyanophages and analyzed the abundances of T4-like and T7-like cyanophage populations in the photic zone of the NPSG in summer 2015 and spring 2016. Combined, the peak abundances of these two cyanophage families reached 2.8 × 10⁶ and 1.1 × 10⁶ cyanophages ⋅ ml^–1 in the summer and spring, respectively. They constituted between 3 and 16% of total virus-like particles (VLPs), comprising a substantial component of the virioplankton in the NPSG. While both cyanophage families were highly abundant, the T4-like cyanophages were generally 1.3–4.4 fold more so. In summer, cyanophages had similar and reproducible distribution patterns with depth. Abundances were relatively low in the upper mixed layer and increased to form a pronounced subsurface peak at 100 m (1.9 × 10⁶ and 9.1 × 10⁵ phages ⋅ ml^–1 for the T4-like and T7-like cyanophages, respectively), coincident with the maximum in Prochlorococcus populations. Less vertical structure in cyanophage abundances was apparent in the spring profile, despite a subsurface peak in Prochlorococcus numbers. In the summer upper mixed layer, cyanophages constituted a smaller proportion of VLPs than below it and cyanophage to cyanobacteria ratios were considerably lower (1.3–2.8) than those of VLPs to bacteria (8.1–21.2). Differences in abundances between the two families and their contribution to VLPs with depth suggest differences in cyanophage production and/or decay processes relative to other members of the virioplankton in the upper mixed layer. These findings highlight the importance of quantifying distinct populations within the virioplankton to gain accurate understanding of their distribution patterns.

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

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

Host Range of Bacteriophages Against a World-Wide Collection of Erwinia amylovora Determined Using a Quantitative PCR Assay

Erwinia amylovora is a globally devastating pathogen of apple, pear, and other Rosaceous plants. The use of lytic bacteriophages for disease management continues to garner attention as a possible supplement or alternative to antibiotics. A quantitative productive host range was established for 10 Erwinia phages using 106 wild type global isolates of E. amylovora, and the closely related Erwinia pyrifoliae, to investigate the potential regional efficacy of these phages within a biopesticide. Each host was individually infected with each of the 10 Erwinia phages and phage production after 8 h incubation was measured using quantitative real time PCR (qPCR) in conjunction with a standardized plasmid. PCR amplicons for all phages used in the study were incorporated into a single plasmid, allowing standardized quantification of the phage genome copy number after the infection process. Nine of the tested phages exhibited a broad host range, replicating their genomes by at least one log in over 88% of tested hosts. Also, every Amygdaloideae infecting E. amylovora host was able to increase at least one phage by three logs. Bacterial hosts isolated in western North America were less susceptible to most phages, as the mean genomic titre produced dropped by nearly two logs, and this phenomenon was strongly correlated to the amount of exopolysaccharide produced by the host. This method of host range analysis is faster and requires less effort than traditional plaque assay techniques, and the resulting quantitative data highlight subtle differences in phage host preference not observable with typical plaque-based host range assays. These quantitative host range data will be useful to determine which phages should be incorporated into a phage-mediated biocontrol formulation to be tested for regional and universal control of E. amylovora.

Resistance in marine cyanobacteria differs against specialist and generalist cyanophages

Long-term coexistence between unicellular cyanobacteria and their lytic viruses (cyanophages) in the oceans is thought to be due to the presence of sensitive cells in which cyanophages reproduce, ultimately killing the cell, while other cyanobacteria survive due to resistance to infection. Here, we investigated resistance in marine cyanobacteria from the genera Synechococcus and Prochlorococcus and compared modes of resistance against specialist and generalist cyanophages belonging to the T7-like and T4-like cyanophage families. Resistance was extracellular in most interactions against specialist cyanophages irrespective of the phage family, preventing entry into the cell. In contrast, resistance was intracellular in practically all interactions against generalist T4-like cyanophages. The stage of intracellular arrest was interaction-specific, halting at various stages of the infection cycle. Incomplete infection cycles proceeded to various degrees of phage genome transcription and translation as well as phage genome replication in numerous interactions. In a particularly intriguing case, intracellular capsid assembly was observed, but the phage genome was not packaged. The cyanobacteria survived the encounter despite late-stage infection and partial genome degradation. We hypothesize that this is tolerated due to genome polyploidy, which we found for certain strains of both Synechococcus and Prochlorococcus Our findings unveil a heavy cost of promiscuous entry of generalist phages into nonhost cells that is rarely paid by specialist phages and suggests the presence of unknown mechanisms of intracellular resistance in the marine unicellular cyanobacteria. Furthermore, these findings indicate that the range for virus-mediated horizontal gene transfer extends beyond hosts to nonhost cyanobacterial cells.

Droplet Digital PCR for Estimating Absolute Abundances of Widespread Pelagibacter Viruses

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

Cyanophage MazG is a pyrophosphohydrolase but unable to hydrolyse magic spot nucleotides

Bacteriophage possess a variety of auxiliary metabolic genes of bacterial origin. These proteins enable them to maximize infection efficiency, subverting bacterial metabolic processes for the purpose of viral genome replication and synthesis of the next generation of virion progeny. Here, we examined the enzymatic activity of a cyanophage MazG protein - a putative pyrophosphohydrolase previously implicated in regulation of the stringent response via reducing levels of the central alarmone molecule (p)ppGpp. We demonstrate, however, that the purified viral MazG shows no binding or hydrolysis activity against (p)ppGpp. Instead, dGTP and dCTP appear to be the preferred substrates of this protein, consistent with a role preferentially hydrolysing deoxyribonucleotides from the high GC content host Synechococcus genome. This showcases a new example of the fine-tuned nature of viral metabolic processes.

Dynamic marine viral infections and major contribution to photosynthetic processes shown by spatiotemporal picoplankton metatranscriptomes

Viruses provide top-down control on microbial communities, yet their direct study in natural environments was hindered by culture limitations. The advance of bioinformatics enables cultivation-independent study of viruses. Many studies assemble new viral genomes and study viral diversity using marker genes from free viruses. Here we use cellular metatranscriptomics to study active community-wide viral infections. Recruitment to viral contigs allows tracking infection dynamics over time and space. Our assemblies represent viral populations, but appear biased towards low diversity viral taxa. Tracking relatives of published T4-like cyanophages and pelagiphages reveals high genomic continuity. We determine potential hosts by matching dynamics of infection with abundance of particular microbial taxa. Finally, we quantify the relative contribution of cyanobacteria and viruses to photosystem-II psbA (reaction center) expression in our study sites. We show sometimes >50% of all cyanobacterial+viral psbA expression is of viral origin, highlighting the contribution of viruses to photosynthesis and oxygen production.

Here, Sieradzki et al. use metatranscriptomics to study active community-wide viral infections at three coastal California sites throughout a year, identify potential viral hosts, and show that viruses can contribute a substantial amount to photosystem-II psbA expression.

Viruses of Eukaryotic Algae: Diversity, Methods for Detection, and Future Directions

The scope for ecological studies of eukaryotic algal viruses has greatly improved with the development of molecular and bioinformatic approaches that do not require algal cultures. Here, we review the history and perceived future opportunities for research on eukaryotic algal viruses. We begin with a summary of the 65 eukaryotic algal viruses that are presently in culture collections, with emphasis on shared evolutionary traits (e.g., conserved core genes) of each known viral type. We then describe how core genes have been used to enable molecular detection of viruses in the environment, ranging from PCR-based amplification to community scale "-omics" approaches. Special attention is given to recent studies that have employed network-analyses of -omics data to predict virus-host relationships, from which a general bioinformatics pipeline is described for this type of approach. Finally, we conclude with acknowledgement of how the field of aquatic virology is adapting to these advances, and highlight the need to properly characterize new virus-host systems that may be isolated using preliminary molecular surveys. Researchers can approach this work using lessons learned from the Chlorella virus system, which is not only the best characterized algal-virus system, but is also responsible for much of the foundation in the field of aquatic virology.

Phage puppet masters of the marine microbial realm

Viruses numerically dominate our oceans; however, we have only just begun to document the diversity, host range and infection dynamics of marine viruses, as well as the subsequent effects of infection on both host cell metabolism and oceanic biogeochemistry. Bacteriophages (that is, phages: viruses that infect bacteria) are highly abundant and are known to play critical roles in bacterial mortality, biogeochemical cycling and horizontal gene transfer. This Review Article summarizes current knowledge of marine viral ecology and highlights the importance of phage particles to the dissolved organic matter pool, as well as the complex interactions between phages and their bacterial hosts. We emphasize the newly recognized roles of phages as puppet masters of their bacterial hosts, where phages are capable of altering the metabolism of infected bacteria through the expression of auxiliary metabolic genes and the redirection of host gene expression patterns. Finally, we propose the 'royal family model' as a hypothesis to describe successional patterns of bacteria and phages over time in marine systems, where despite high richness and significant seasonal differences, only a small number of phages appear to continually dominate a given marine ecosystem. Although further testing is required, this model provides a framework for assessing the specificity and ecological consequences of phage-host dynamics.