Viruses impact microbial systems through killing hosts, horizontal gene transfer, and altering cellular metabolism, consequently impacting nutrient cycles. A virus-infected cell, a “virocell,” is distinct from its uninfected sister cell as the virus commandeers cellular machinery to produce viruses rather than replicate cells. Problematically, virocell responses to the nutrient-limited conditions that abound in nature are poorly understood. Here we used a systems biology approach to investigate virocell metabolic reprogramming under nutrient limitation. Using transcriptomics, proteomics, lipidomics, and endo- and exo-metabolomics, we assessed how low phosphate (low-P) conditions impacted virocells of a marine Pseudoalteromonas host when independently infected by two unrelated phages (HP1 and HS2). With the combined stresses of infection and nutrient limitation, a set of nested responses were observed. First, low-P imposed common cellular responses on all cells (virocells and uninfected cells), including activating the canonical P-stress response, and decreasing transcription, translation, and extracellular organic matter consumption. Second, low-P imposed infection-specific responses (for both virocells), including enhancing nitrogen assimilation and fatty acid degradation, and decreasing extracellular lipid relative abundance. Third, low-P suggested virocell-specific strategies. Specifically, HS2-virocells regulated gene expression by increasing transcription and ribosomal protein production, whereas HP1-virocells accumulated host proteins, decreased extracellular peptide relative abundance, and invested in broader energy and resource acquisition. These results suggest that although environmental conditions shape metabolism in common ways regardless of infection, virocell-specific strategies exist to support viral replication during nutrient limitation, and a framework now exists for identifying metabolic strategies of nutrient-limited virocells in nature.

Summary: Microbes drive ecosystem functioning and their viruses modulate these impacts through mortality, gene transfer and metabolic reprogramming. Despite the importance of virus‐host interactions and likely variable infection efficiencies of individual phages across hosts, such variability is seldom quantified. Here, we quantify infection efficiencies of 38 phages against 19 host strains in aquatic Cellulophaga ( Bacteroidetes ) phage‐host model systems. Binary data revealed that some phages infected only one strain while others infected 17, whereas quantitative data revealed that efficiency of infection could vary 10 orders of magnitude, even among phages within one population. This provides a baseline for understanding and modeling intrapopulation host range variation. Genera specific host ranges were also informative. For example, the Cellulophaga Microviridae , showed a markedly broader intra‐species host range than previously observed in Escherichia coli systems. Further, one phage genus, Cba41, was examined to investigate nonheritable changes in plating efficiency and burst size that depended on which host strain it most recently infected. While consistent with host modification of phage DNA, no differences in nucleotide sequence or DNA modifications were detected, leaving the observation repeatable, but the mechanism unresolved. Overall, this study highlights the importance of quantitatively considering replication variations in studies of phage–host interactions.

Hemolytic–uremic syndrome is a life-threating disease most often associated with Shiga toxin-producing microorganisms like Escherichia coli (STEC), including E. coli O157:H7. Shiga toxin is encoded by resident prophages present within this bacterium, and both its production and release depend on the induction of Shiga toxin-encoding prophages. Consequently, treatment of STEC infections tend to be largely supportive rather than antibacterial, in part due to concerns about exacerbating such prophage induction. Here we explore STEC O157:H7 prophage induction in vitro as it pertains to phage therapy—the application of bacteriophages as antibacterial agents to treat bacterial infections—to curtail prophage induction events, while also reducing STEC O157:H7 presence. We observed that cultures treated with strictly lytic phages, despite being lysed, produce substantially fewer Shiga toxin-encoding temperate-phage virions than untreated STEC controls. We therefore suggest that phage therapy could have utility as a prophylactic treatment of individuals suspected of having been recently exposed to STEC, especially if prophage induction and by extension Shiga toxin production is not exacerbated.

Bacterial viruses (phages) are abundant and ecologically impactful, but laboratory-based experimental model systems vastly under-represent known phage diversity, particularly for ssDNA phages. Here, we characterize the genomes and infection properties of two unrelated marine flavophages—ssDNA generalist phage phi18:4 (6.5 Kbp) and dsDNA specialist phage phi18:1 (39.2 Kbp)—when infecting the same Cellulophaga baltica strain #18 (Cba18), of the class Flavobacteriia. Phage phi18:4 belongs to a new family of ssDNA phages, has an internal lipid membrane, and its genome encodes primarily structural proteins, as well as a DNA replication protein common to ssDNA phages and a unique lysis protein. Phage phi18:1 is a siphovirus that encodes several virulence genes, despite not having a known temperate lifestyle, a CAZy enzyme likely for regulatory purposes, and four DNA methyltransferases dispersed throughout the genome that suggest both host modulation and phage DNA protection against host restriction. Physiologically, ssDNA phage phi18:4 has a shorter latent period and smaller burst size than dsDNA phage phi18:1, and both phages efficiently infect this host. These results help augment the diversity of characterized environmental phage–host model systems by studying infections of genomically diverse phages (ssDNA vs. dsDNA) on the same host.

The world's ecosystems are shaped by microbiota. Their niches and their impacts depend on functional profiles influenced by gene gains and losses. While culture-based experiments demonstrate that mobile genetic elements (MGEs) can mediate gene flux, quantitative field data on the rates and impacts of MGE activity remains scarce. Here we leverage large-scale soil meta-omic data to develop and apply analytics for studying MGEs in complex natural systems. In our model permafrost-thaw ecosystem, Stordalen Mire, we identify ~2.1 million MGE recombinases across 89 microbial phyla to assess ecological distributions, affected functions, past mobility, and current activity. This revealed MGEs shaping natural genetic diversity via differential impacts on major phyla; affecting a wide range of functions, including diverse regulatory and metabolic genes affecting carbon flux and nutrient cycling; and moving at rates that should significantly influence the realized functional profiles of natural microbiomes. These findings and this systematic meta-omic framework open new avenues to better investigate MGE diversity, activity, mobility, and impacts in nature.