Low activity of lytic pelagiphages in coastal marine waters

Long-read powered viral metagenomics in the oligotrophic Sargasso Sea

Dominant microorganisms of the Sargasso Sea are key drivers of the global carbon cycle. However, associated viruses that shape microbial community structure and function are not well characterised. Here, we combined short and long read sequencing to survey Sargasso Sea phage communities in virus- and cellular fractions at viral maximum (80 m) and mesopelagic (200 m) depths. We identified 2,301 Sargasso Sea phage populations from 186 genera. Over half of the phage populations identified here lacked representation in global ocean viral metagenomes, whilst 177 of the 186 identified genera lacked representation in genomic databases of phage isolates. Viral fraction and cell-associated viral communities were decoupled, indicating viral turnover occurred across periods longer than the sampling period of three days. Inclusion of long-read data was critical for capturing the breadth of viral diversity. Phage isolates that infect the dominant bacterial taxa Prochlorococcus and Pelagibacter, usually regarded as cosmopolitan and abundant, were poorly represented.

Globally occurring pelagiphage infections create ribosome-deprived cells

Phages play an essential role in controlling bacterial populations. Those infecting Pelagibacterales (SAR11), the dominant bacteria in surface oceans, have been studied in silico and by cultivation attempts. However, little is known about the quantity of phage-infected cells in the environment. Using fluorescence in situ hybridization techniques, we here show pelagiphage-infected SAR11 cells across multiple global ecosystems and present evidence for tight community control of pelagiphages on the SAR11 hosts in a case study. Up to 19% of SAR11 cells were phage-infected during a phytoplankton bloom, coinciding with a ~90% reduction in SAR11 cell abundance within 5 days. Frequently, a fraction of the infected SAR11 cells were devoid of detectable ribosomes, which appear to be a yet undescribed possible stage during pelagiphage infection. We dubbed such cells zombies and propose, among other possible explanations, a mechanism in which ribosomal RNA is used as a resource for the synthesis of new phage genomes. On a global scale, we detected phage-infected SAR11 and zombie cells in the Atlantic, Pacific, and Southern Oceans. Our findings illuminate the important impact of pelagiphages on SAR11 populations and unveil the presence of ribosome-deprived zombie cells as part of the infection cycle.

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.

Marine viruses and climate change: Virioplankton, the carbon cycle, and our future ocean

Interactions between marine viruses and microbes are a critical part of the oceanic carbon cycle. The impacts of virus-host interactions range from short-term disruptions in the mobility of microbial biomass carbon to higher trophic levels through cell lysis (i.e., the viral shunt) to long-term reallocation of microbial biomass carbon to the deep sea through accelerating the biological pump (i.e., the viral shuttle). The biogeochemical backdrop of the ocean-the physical, chemical, and biological landscape-influences the likelihood of both virus-host interactions and particle formation, and the fate and flow of carbon. As climate change reshapes the oceanic landscape through large-scale shifts in temperature, circulation, stratification, and acidification, virus-mediated carbon flux is likely to shift in response. Dynamics in the directionality and magnitude of changes in how, where, and when viruses mediate the recycling or storage of microbial biomass carbon is largely unknown. Integrating viral infection dynamics data obtained from experimental models and field systems, with particle motion microphysics and global observations of oceanic biogeochemistry, into improved ecosystem models will enable viral oceanographers to better predict the role of viruses in marine carbon cycling in the future ocean.

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.

Dynamics of Baltic Sea phages driven by environmental changes

Phage predation constitutes a major mortality factor for bacteria in aquatic ecosystems, and thus, directly impacts nutrient cycling and microbial community dynamics. Yet, the population dynamics of specific phages across time scales from days to months remain largely unexplored, which limits our understanding of their influence on microbial succession. To investigate temporal changes in diversity and abundance of phages infecting particular host strains, we isolated 121 phage strains that infected three bacterial hosts during a Baltic Sea mesocosm experiment. Genome analysis revealed a novel Flavobacterium phage genus harboring gene sets putatively coding for synthesis of modified nucleotides and glycosylation of bacterial cell surface components. Another novel phage genus revealed a microdiversity of phage species that was largely maintained during the experiment and across mesocosms amended with different nutrients. In contrast to the newly described Flavobacterium phages, phages isolated from a Rheinheimera strain were highly similar to previously isolated genotypes, pointing to genomic consistency in this population. In the mesocosm experiment, the investigated phages were mainly detected after a phytoplankton bloom peak. This concurred with recurrent detection of the phages in the Baltic Proper during summer months, suggesting an influence on the succession of heterotrophic bacteria associated with phytoplankton blooms.

Efficient dilution-to-extinction isolation of novel virus-host model systems for fastidious heterotrophic bacteria

Microbes and their associated viruses are key drivers of biogeochemical processes in marine and soil biomes. While viruses of phototrophic cyanobacteria are well-represented in model systems, challenges of isolating marine microbial heterotrophs and their viruses have hampered experimental approaches to quantify the importance of viruses in nutrient recycling. A resurgence in cultivation efforts has improved the availability of fastidious bacteria for hypothesis testing, but this has not been matched by similar efforts to cultivate their associated bacteriophages. Here, we describe a high-throughput method for isolating important virus-host systems for fastidious heterotrophic bacteria that couples advances in culturing of hosts with sequential enrichment and isolation of associated phages. Applied to six monthly samples from the Western English Channel, we first isolated one new member of the globally dominant bacterial SAR11 clade and three new members of the methylotrophic bacterial clade OM43. We used these as bait to isolate 117 new phages, including the first known siphophage-infecting SAR11, and the first isolated phage for OM43. Genomic analyses of 13 novel viruses revealed representatives of three new viral genera, and infection assays showed that the viruses infecting SAR11 have ecotype-specific host ranges. Similar to the abundant human-associated phage ɸCrAss001, infection dynamics within the majority of isolates suggested either prevalent lysogeny or chronic infection, despite a lack of associated genes, or host phenotypic bistability with lysis putatively maintained within a susceptible subpopulation. Broader representation of important virus-host systems in culture collections and genomic databases will improve both our understanding of virus-host interactions, and accuracy of computational approaches to evaluate ecological patterns from metagenomic data.

Seasonal dynamics of natural Ostreococcus viral infection at the single cell level using VirusFISH

Ostreococcus is a cosmopolitan marine genus of phytoplankton found in mesotrophic and oligotrophic waters, and the smallest free-living eukaryotes known to date, with a cell diameter close to 1 μm. Ostreococcus has been extensively studied as a model system to investigate viral-host dynamics in culture, yet the impact of viruses in naturally occurring populations is largely unknown. Here, we used Virus Fluorescence in situ Hybridization (VirusFISH) to visualize and quantify viral-host dynamics in natural populations of Ostreococcus during a seasonal cycle in the central Cantabrian Sea (Southern Bay of Biscay). Ostreococcus were predominantly found during summer and autumn at surface and 50 m depth, in coastal, mid-shelf and shelf waters, representing up to 21% of the picoeukaryotic communities. Viral infection was only detected in surface waters, and its impact was variable but highest from May to July and November to December, when up to half of the population was infected. Metatranscriptomic data available from the mid-shelf station unveiled that the Ostreococcus population was dominated by the species O. lucimarinus. This work represents a proof of concept that the VirusFISH technique can be used to quantify the impact of viruses on targeted populations of key microbes from complex natural communities.

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.

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

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

Single-virus genomics and beyond

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

Phages as a Cohesive Prophylactic and Therapeutic Approach in Aquaculture Systems

Facing antibiotic resistance has provoked a continuously growing focus on phage therapy. Although the greatest emphasis has always been placed on phage treatment in humans, behind phage application lies a complex approach that can be usefully adopted by the food industry, from hatcheries and croplands to ready-to-eat products. Such diverse businesses require an efficient method for combating highly pathogenic bacteria since antibiotic resistance concerns every aspect of human life. Despite the vast abundance of phages on Earth, the aquatic environment has been considered their most natural habitat. Water favors multidirectional Brownian motion and increases the possibility of contact between phage particles and their bacterial hosts. As the global production of aquatic organisms has rapidly grown over the past decades, phage treatment of bacterial infections seems to be an obvious and promising solution in this market sector. Pathogenic bacteria, such as Aeromonas and Vibrio, have already proved to be responsible for mass mortalities in aquatic systems, resulting in economic losses. The main objective of this work is to summarize, from a scientific and industry perspective, the recent data regarding phage application in the form of targeted probiotics and therapeutic agents in aquaculture niches.

Expanding the Diversity of Bacterioplankton Isolates and Modeling Isolation Efficacy with Large-Scale Dilution-to-Extinction Cultivation

Even before the coining of the term “great plate count anomaly” in the 1980s, scientists had noted the discrepancy between the number of microorganisms observed under the microscope and the number of colonies that grew on traditional agar media. New cultivation approaches have reduced this disparity, resulting in the isolation of some of the “most wanted” bacterial lineages. Nevertheless, the vast majority of microorganisms remain uncultured, hampering progress toward answering fundamental biological questions about many important microorganisms. Furthermore, few studies have evaluated the underlying factors influencing cultivation success, limiting our ability to improve cultivation efficacy. Our work details the use of dilution-to-extinction (DTE) cultivation to expand the phylogenetic and geographic diversity of available axenic cultures. We also provide a new model of the DTE approach that uses cultivation results and natural abundance information to predict taxon-specific viability and iteratively constrain DTE experimental design to improve cultivation success.

Cultivated bacterioplankton representatives from diverse lineages and locations are essential for microbiology, but the large majority of taxa either remain uncultivated or lack isolates from diverse geographic locales. We paired large-scale dilution-to-extinction (DTE) cultivation with microbial community analysis and modeling to expand the phylogenetic and geographic diversity of cultivated bacterioplankton and to evaluate DTE cultivation success. Here, we report results from 17 DTE experiments totaling 7,820 individual incubations over 3 years, yielding 328 repeatably transferable isolates. Comparison of isolates to microbial community data for source waters indicated that we successfully isolated 5% of the observed bacterioplankton community throughout the study; 43% and 26% of our isolates matched operational taxonomic units and amplicon single-nucleotide variants, respectively, within the top 50 most abundant taxa. Isolates included those from previously uncultivated clades such as SAR11 LD12 and Actinobacteria acIV, as well as geographically novel members from other ecologically important groups like SAR11 subclade IIIa, SAR116, and others, providing isolates in eight putatively new genera and seven putatively new species. Using a newly developed DTE cultivation model, we evaluated taxon viability by comparing relative abundance with cultivation success. The model (i) revealed the minimum attempts required for successful isolation of taxa amenable to growth on our media and (ii) identified possible subpopulation viability variation in abundant taxa such as SAR11 that likely impacts cultivation success. By incorporating viability in experimental design, we can now statistically constrain the effort necessary for successful cultivation of specific taxa on a defined medium.

IMPORTANCE Even before the coining of the term “great plate count anomaly” in the 1980s, scientists had noted the discrepancy between the number of microorganisms observed under the microscope and the number of colonies that grew on traditional agar media. New cultivation approaches have reduced this disparity, resulting in the isolation of some of the “most wanted” bacterial lineages. Nevertheless, the vast majority of microorganisms remain uncultured, hampering progress toward answering fundamental biological questions about many important microorganisms. Furthermore, few studies have evaluated the underlying factors influencing cultivation success, limiting our ability to improve cultivation efficacy. Our work details the use of dilution-to-extinction (DTE) cultivation to expand the phylogenetic and geographic diversity of available axenic cultures. We also provide a new model of the DTE approach that uses cultivation results and natural abundance information to predict taxon-specific viability and iteratively constrain DTE experimental design to improve cultivation success.

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.

Transcriptional Patterns of Biogeochemically Relevant Marker Genes by Temperate Marine Bacteria

Environmental microbial gene expression patterns remain largely unexplored, particularly at interannual time scales. We analyzed the variability in the expression of marker genes involved in ecologically relevant biogeochemical processes at a temperate Atlantic site over two consecutive years. Most of nifH transcripts, involved in nitrogen (N) fixation, were affiliated with the symbiotic cyanobacterium Candidatus Atelocyanobacterium thalassa, suggesting a key role as N providers in this system. The expression of nifH and amoA (i.e., marker for ammonia oxidation) showed consistent maxima in summer and autumn, respectively, suggesting a temporal succession of these important N cycling processes. The patterns of expression of genes related to the oxidation of carbon monoxide (coxL) and reduced sulfur (soxB) were different from that of amoA, indicating alternate timings for these energy conservation strategies. We detected expression of alkaline phosphatases, induced under phosphorus limitation, in agreement with the reported co-limitation by this nutrient at the study site. In contrast, low-affinity phosphate membrane transporters (pit) typically expressed under phosphorus luxury conditions, were mainly detected in post-bloom conditions. Rhodobacteraceae dominated the expression of soxB, coxL and ureases, while Pelagibacteraceae dominated the expression of proteorhodopsins. Bacteroidetes and Gammaproteobacteria were major contributors to the uptake of inorganic nutrients (pit and amt transporters). Yet, in autumn, Thauma- and Euryarchaeota unexpectedly contributed importantly to the uptake of ammonia and phosphate, respectively. We provide new hints on the active players and potential dynamics of ecologically relevant functions in situ, highlighting the potential of metatranscriptomics to provide significant input to future omics-driven marine ecosystem assessment.

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.

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.