Lysogenic host-virus interactions in SAR11 marine bacteria

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.

A systematic analysis of marine lysogens and proviruses

Viruses are ubiquitous in the oceans, exhibiting high abundance and diversity. Here, we systematically analyze existing genomic sequences of marine prokaryotes to compile a Marine Prokaryotic Genome Dataset (MPGD, consisting of over 12,000 bacterial and archaeal genomes) and a Marine Temperate Viral Genome Dataset (MTVGD). At least 40% of the MPGD genomes contain one or more proviral sequences, indicating that they are lysogens. The MTVGD includes over 12,900 viral contigs or putative proviruses, clustered into 10,897 viral genera. We show that lysogens and proviruses are abundant in marine ecosystems, particularly in the deep sea, and marine lysogens differ from non-lysogens in multiple genomic features and growth properties. We reveal several virus-host interaction networks of potential ecological relevance, and identify proviruses that appear to be able to infect (or to be transferred between) different bacterial classes and phyla. Auxiliary metabolic genes in the MTVGD are enriched in functions related to carbohydrate metabolism. Finally, we experimentally demonstrate the impact of a prophage on the transcriptome of a representative marine Shewanella bacterium. Our work contributes to a better understanding of the ecology of marine prokaryotes and their viruses.

Virus-Host Interactions Drive Contrasting Bacterial Diel Dynamics in the Ocean

Marine organisms perform a sea of diel rhythmicity. Planktonic diel dynamics have been shown to be driven by light, energy resources, circadian rhythms, and the coordinated coupling of photoautotrophs and heterotrophic bacterioplankton. Here, we explore the diel fluctuation of viral production and decay and their impact on the total and active bacterial community in the coastal and open seawaters of the South China Sea. The results showed that the night-production diel pattern of lytic viral production was concurrent with the lower viral decay at night, contributing to the accumulation of the viral population size during the night for surface waters. The diel variations in bacterial activity, community composition, and diversity were found highly affected by viral dynamics. This was revealed by the finding that bacterial community diversity was positively correlated to lytic viral production in the euphotic zone of the open ocean but was negatively related to lysogenic viral production in the coastal ocean. Such distinct but contrasting correlations suggest that viral life strategies can not only contribute to diversifying bacterial community but also potentially piggyback their host to dominate bacterial community, suggesting the tightly synchronized depth-dependent and habitat-specific diel patterns of virus-host interactions. It further implies that viruses serve as an ecologically important driver of bacterial diel dynamics across the ocean, highlighting the viral roles in bacterial ecological and biogeochemical processes in the ocean.

In situ cell division and mortality rates of SAR11, SAR86, Bacteroidetes, and Aurantivirga during phytoplankton blooms reveal differences in population controls

Net growth of microbial populations, that is, changes in abundances over time, can be studied using 16S rRNA fluorescence in situ hybridization (FISH). However, this approach does not differentiate between mortality and cell division rates. We used FISH-based image cytometry in combination with dilution culture experiments to study net growth, cell division, and mortality rates of four bacterial taxa over two distinct phytoplankton blooms: the oligotrophs SAR11 and SAR86, and the copiotrophic phylum Bacteroidetes, and its genus Aurantivirga. Cell volumes, ribosome content, and frequency of dividing cells (FDC) co-varied over time. Among the three, FDC was the most suitable predictor to calculate cell division rates for the selected taxa. The FDC-derived cell division rates for SAR86 of up to 0.8/day and Aurantivirga of up to 1.9/day differed, as expected for oligotrophs and copiotrophs. Surprisingly, SAR11 also reached high cell division rates of up to 1.9/day, even before the onset of phytoplankton blooms. For all four taxonomic groups, the abundance-derived net growth (-0.6 to 0.5/day) was about an order of magnitude lower than the cell division rates. Consequently, mortality rates were comparably high to cell division rates, indicating that about 90% of bacterial production is recycled without apparent time lag within 1 day. Our study shows that determining taxon-specific cell division rates complements omics-based tools and provides unprecedented clues on individual bacterial growth strategies including bottom-up and top-down controls. IMPORTANCE The growth of a microbial population is often calculated from their numerical abundance over time. However, this does not take cell division and mortality rates into account, which are important for deriving ecological processes like bottom-up and top-down control. In this study, we determined growth by numerical abundance and calibrated microscopy-based methods to determine the frequency of dividing cells and subsequently calculate taxon-specific cell division rates in situ. The cell division and mortality rates of two oligotrophic (SAR11 and SAR86) and two copiotrophic (Bacteroidetes and Aurantivirga) taxa during two spring phytoplankton blooms showed a tight coupling for all four taxa throughout the blooms without any temporal offset. Unexpectedly, SAR11 showed high cell division rates days before the bloom while cell abundances remained constant, which is indicative of strong top-down control. Microscopy remains the method of choice to understand ecological processes like top-down and bottom-up control on a cellular level.

The Promise and Pitfalls of Prophages

Phages dominate every ecosystem on the planet. While virulent phages sculpt the microbiome by killing their bacterial hosts, temperate phages provide unique growth advantages to their hosts through lysogenic conversion. Many prophages benefit their host, and prophages are responsible for genotypic and phenotypic differences that separate individual microbial strains. However, the microbes also endure a cost to maintain those phages: additional DNA to replicate and proteins to transcribe and translate. We have never quantified those benefits and costs. Here, we analysed over two and a half million prophages from over half a million bacterial genome assemblies. Analysis of the whole dataset and a representative subset of taxonomically diverse bacterial genomes demonstrated that the normalised prophage density was uniform across all bacterial genomes above 2 Mbp. We identified a constant carrying capacity of phage DNA per bacterial DNA. We estimated that each prophage provides cellular services equivalent to approximately 2.4 % of the cell's energy or 0.9 ATP per bp per hour. We demonstrate analytical, taxonomic, geographic, and temporal disparities in identifying prophages in bacterial genomes that provide novel targets for identifying new phages. We anticipate that the benefits bacteria accrue from the presence of prophages balance the energetics involved in supporting prophages. Furthermore, our data will provide a new framework for identifying phages in environmental datasets, diverse bacterial phyla, and from different locations.

Prevalence of Viral Frequency-Dependent Infection in Coastal Marine Prokaryotes Revealed Using Monthly Time Series Virome Analysis

Viruses infecting marine prokaryotes have a large impact on the diversity and dynamics of their hosts. Model systems suggest that viral infection is frequency dependent and constrained by the virus-host encounter rate. However, it is unclear whether frequency-dependent infection is pervasive among the abundant prokaryotic populations with different temporal dynamics. To address this question, we performed a comparison of prokaryotic and viral communities using 16S rRNA amplicon and virome sequencing based on samples collected monthly for 2 years at a Japanese coastal site, Osaka Bay. Concurrent seasonal shifts observed in prokaryotic and viral community dynamics indicated that the abundance of viruses correlated with that of their predicted host phyla (or classes). Cooccurrence network analysis between abundant prokaryotes and viruses revealed 6,423 cooccurring pairs, suggesting a tight coupling of host and viral abundances and their "one-to-many" correspondence. Although stable dominant species, such as SAR11, showed few cooccurring viruses, a fast succession of their viruses suggests that viruses infecting these populations changed continuously. Our results suggest that frequency-dependent viral infection prevails in coastal marine prokaryotes regardless of host taxa and temporal dynamics. IMPORTANCE There is little room for doubt that viral infection is prevalent among abundant marine prokaryotes regardless of their taxa or growth strategy. However, comprehensive evaluations of viral infections in natural prokaryotic communities are still technically difficult. In this study, we examined viral infection in abundant prokaryotes by monitoring the monthly dynamics of prokaryotic and viral communities at a eutrophic coastal site, Osaka Bay. We compared the community dynamics of viruses with those of their putative hosts based on genome-based in silico host prediction. We observed frequent cooccurrence among the predicted virus-host pairs, suggesting that viral infection is prevalent in abundant prokaryotes regardless of their taxa or temporal dynamics. This likely indicates that frequent lysis of the abundant prokaryotes via viral infection has a considerable contribution to the biogeochemical cycling and maintenance of prokaryotic community diversity.

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

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 × 1026 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.

Marine DNA methylation patterns are associated with microbial community composition and inform virus-host dynamics

BACKGROUND

DNA methylation in prokaryotes is involved in many different cellular processes including cell cycle regulation and defense against viruses. To date, most prokaryotic methylation systems have been studied in culturable microorganisms, resulting in a limited understanding of DNA methylation from a microbial ecology perspective. Here, we analyze the distribution patterns of several microbial epigenetics marks in the ocean microbiome through genome-centric metagenomics across all domains of life.

RESULTS

We reconstructed 15,056 viral, 252 prokaryotic, 56 giant viral, and 6 eukaryotic metagenome-assembled genomes from northwest Pacific Ocean seawater samples using short- and long-read sequencing approaches. These metagenome-derived genomes mostly represented novel taxa, and recruited a majority of reads. Thanks to single-molecule real-time (SMRT) sequencing technology, base modification could also be detected for these genomes. This showed that DNA methylation can readily be detected across dominant oceanic bacterial, archaeal, and viral populations, and microbial epigenetic changes correlate with population differentiation. Furthermore, our genome-wide epigenetic analysis of Pelagibacter suggests that GANTC, a DNA methyltransferase target motif, is related to the cell cycle and is affected by environmental conditions. Yet, the presence of this motif also partitions the phylogeny of the Pelagibacter phages, possibly hinting at a competitive co-evolutionary history and multiple effects of a single methylation mark.

CONCLUSIONS

Overall, this study elucidates that DNA methylation patterns are associated with ecological changes and virus-host dynamics in the ocean microbiome. Video Abstract.

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.

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.

Genetic engineering of marine cyanophages reveals integration but not lysogeny in T7-like cyanophages

Marine cyanobacteria of the genera Synechococcus and Prochlorococcus are the most abundant photosynthetic organisms on earth, spanning vast regions of the oceans and contributing significantly to global primary production. Their viruses (cyanophages) greatly influence cyanobacterial ecology and evolution. Although many cyanophage genomes have been sequenced, insight into the functional role of cyanophage genes is limited by the lack of a cyanophage genetic engineering system. Here, we describe a simple, generalizable method for genetic engineering of cyanophages from multiple families, that we named REEP for REcombination, Enrichment and PCR screening. This method enables direct investigation of key cyanophage genes, and its simplicity makes it adaptable to other ecologically relevant host-virus systems. T7-like cyanophages often carry integrase genes and attachment sites, yet exhibit lytic infection dynamics. Here, using REEP, we investigated their ability to integrate and maintain a lysogenic life cycle. We found that these cyanophages integrate into the host genome and that the integrase and attachment site are required for integration. However, stable lysogens did not form. The frequency of integration was found to be low in both lab cultures and the oceans. These findings suggest that T7-like cyanophage integration is transient and is not part of a classical lysogenic cycle.

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.

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

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

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

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

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.

Ocean Acidification Induces Changes in Virus-Host Relationships in Mediterranean Benthic Ecosystems

Acidified marine systems represent "natural laboratories", which provide opportunities to investigate the impacts of ocean acidification on different living components, including microbes. Here, we compared the benthic microbial response in four naturally acidified sites within the Southern Tyrrhenian Sea characterized by different acidification sources (i.e., CO₂ emissions at Ischia, mixed gases at Panarea and Basiluzzo and acidified freshwater from karst rocks at Presidiana) and pH values. We investigated prokaryotic abundance, activity and biodiversity, viral abundance and prokaryotic infections, along with the biochemical composition of the sediment organic matter. We found that, despite differences in local environmental dynamics, viral life strategies change in acidified conditions from mainly lytic to temperate lifestyles (e.g., chronic infection), also resulting in a lowered impact on prokaryotic communities, which shift towards (chemo)autotrophic assemblages, with lower organic matter consumption. Taken together, these results suggest that ocean acidification exerts a deep control on microbial benthic assemblages, with important feedbacks on ecosystem functioning.

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.

Culturing novel and abundant pelagiphages in the ocean

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

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

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

Quantification of Lysogeny Caused by Phage Coinfections in Microbial Communities from Biophysical Principles

The association of temperate phages and bacterial hosts during lysogeny manipulates microbial dynamics from the oceans to the human gut. Lysogeny is well studied in laboratory models, but its environmental drivers remain unclear. Here, we quantified the probability of lysogenization caused by phage coinfections, a well-known trigger of lysogeny, in marine and gut microbial environments. Coinfections were quantified by developing a biophysical model that incorporated the traits of viral and bacterial communities. Lysogenization via coinfection was more frequent in highly productive environments like the gut, due to higher microbial densities and higher phage adsorption rates. At low cell densities, lysogenization occurred in bacteria with long duplication times. These results bridge the molecular understanding of lysogeny with the ecology of complex microbial communities.

Temperate phages can associate with their bacterial host to form a lysogen, often modifying the phenotype of the host. Lysogens are dominant in the microbially dense environment of the mammalian gut. This observation contrasts with the long-standing hypothesis of lysogeny being favored at low microbial densities, such as in oligotrophic marine environments. Here, we hypothesized that phage coinfections—a well-understood molecular mechanism of lysogenization—increase at high microbial abundances. To test this hypothesis, we developed a biophysical model of coinfection for marine and gut microbiomes. The model stochastically sampled ranges of phage and bacterial concentrations, adsorption rates, lysogenic commitment times, and community diversity from each environment. In 90% of the sampled marine communities, less than 10% of the bacteria were predicted to be lysogenized via coinfection. In contrast, 25% of the sampled gut communities displayed more than 25% of lysogenization. The probability of lysogenization in the gut was a consequence of the higher densities and higher adsorption rates. These results suggest that, on average, coinfections can form two trillion lysogens in the human gut every day. In marine microbiomes, which were characterized by lower densities and phage adsorption rates, lysogeny via coinfection was still possible for communities with long lysogenic commitment times. Our study indicates that different physical factors causing coinfections can reconcile the traditional view of lysogeny at poor host growth (long commitment times) and the recent Piggyback-the-Winner framework proposing that lysogeny is favored in rich environments (high densities and adsorption rates).

IMPORTANCE The association of temperate phages and bacterial hosts during lysogeny manipulates microbial dynamics from the oceans to the human gut. Lysogeny is well studied in laboratory models, but its environmental drivers remain unclear. Here, we quantified the probability of lysogenization caused by phage coinfections, a well-known trigger of lysogeny, in marine and gut microbial environments. Coinfections were quantified by developing a biophysical model that incorporated the traits of viral and bacterial communities. Lysogenization via coinfection was more frequent in highly productive environments like the gut, due to higher microbial densities and higher phage adsorption rates. At low cell densities, lysogenization occurred in bacteria with long duplication times. These results bridge the molecular understanding of lysogeny with the ecology of complex microbial communities.

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.