The genomic content and context of auxiliary metabolic genes in marine cyanomyoviruses

Genomic analysis of Synechococcus phage S-B43 and its adaption to the coastal environment

Synechococcus dominate picocyanobacterial communities in coastal environments. However, only a few Synechococcus phages have been described from the coastal seas of the Northwest Pacific Ocean. Here a new Synechococcus phage, S-B43 was isolated from the Bohai Sea, a semi-closed coastal sea of the Northwest Pacific Ocean. S-B43 is a member of Myoviridae, containing 275 predicted open reading frames. Fourteen auxiliary metabolic genes (AMG) were identified from the genome of S-B43, including five photosynthetic associated genes and several AMGs related to its adaption to the high turbidity and eutrophic coastal environment with a low ratio of phosphorus to nitrogen (HNLP). The occurrences of 31 AMGs among 34 cyanophage genomes indicates that AMGs zwf, gnd, speD, petF and those coding for FECH and thioredoxin were more common in coastal areas than in the open ocean and AMGs pebS and ho1 were more prevalent in the open ocean. The occurrence of cyanophage AMGs in different environments might be a reflection of the environmental adaption of their hosts. This study contributes to our understanding of the interactions between cyanobacteria and cyanophages and their environmental adaption to the coastal environment.

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.

Exploring the Remarkable Diversity of Culturable Escherichia coli Phages in the Danish Wastewater Environment

Phages drive bacterial diversity, profoundly influencing microbial communities, from microbiomes to the drivers of global biogeochemical cycling. Aiming to broaden our understanding of Escherichiacoli (MG1655, K-12) phages, we screened 188 Danish wastewater samples and isolated 136 phages. Ninety-two of these have genomic sequences with less than 95% similarity to known phages, while most map to existing genera several represent novel lineages. The isolated phages are highly diverse, estimated to represent roughly one-third of the true diversity of culturable virulent dsDNA Escherichia phages in Danish wastewater, yet almost half (40%) are not represented in metagenomic databases, emphasising the importance of isolating phages to uncover diversity. Seven viral families, Myoviridae, Siphoviridae, Podoviridae,Drexlerviridae,Chaseviridae,Autographviridae, and Microviridae, are represented in the dataset. Their genomes vary drastically in length from 5.3 kb to 170.8 kb, with a guanine and cytosine (GC) content ranging from 35.3% to 60.0%. Hence, even for a model host bacterium, substantial diversity remains to be uncovered. These results expand and underline the range of coliphage diversity and demonstrate how far we are from fully disclosing phage diversity and ecology.

Genome Analysis of Two Novel Synechococcus Phages That Lack Common Auxiliary Metabolic Genes: Possible Reasons and Ecological Insights by Comparative Analysis of Cyanomyoviruses

The abundant and widespread unicellular cyanobacteria Synechococcus plays an important role in contributing to global phytoplankton primary production. In the present study, two novel cyanomyoviruses, S-N03 and S-H34 that infected Synechococcus MW02, were isolated from the coastal waters of the Yellow Sea. S-N03 contained a 167,069-bp genome comprising double-stranded DNA with a G + C content of 50.1%, 247 potential open reading frames and 1 tRNA; S-H34 contained a 167,040-bp genome with a G + C content of 50.1%, 246 potential open reading frames and 5 tRNAs. These two cyanophages contain fewer auxiliary metabolic genes (AMGs) than other previously isolated cyanophages. S-H34 in particular, is currently the only known cyanomyovirus that does not contain any AMGs related to photosynthesis. The absence of such common AMGs in S-N03 and S-H34, their distinct evolutionary history and ecological features imply that the energy for phage production might be obtained from other sources rather than being strictly dependent on the maintenance of photochemical ATP under high light. Phylogenetic analysis showed that the two isolated cyanophages clustered together and had a close relationship with two other cyanophages of low AMG content. Comparative genomic analysis, habitats and hosts across 81 representative cyanomyovirus showed that cyanomyovirus with less AMGs content all belonged to Synechococcus phages isolated from eutrophic waters. The relatively small genome size and high G + C content may also relate to the lower AMG content, as suggested by the significant correlation between the number of AMGs and G + C%. Therefore, the lower content of AMG in S-N03 and S-H34 might be a result of viral evolution that was likely shaped by habitat, host, and their genomic context. The genomic content of AMGs in cyanophages may have adaptive significance and provide clues to their evolution.

Temporal transcriptional patterns of cyanophage genes suggest synchronized infection of cyanobacteria in the oceans

BACKGROUND

Based on the peak expression times during infection, early, middle, and late genes have been characterized in viruses (cyanophages) that infect the unicellular cyanobacterium Prochlorococcus. Laboratory experiments show that some cyanophages can only replicate in the light and thus exhibit diurnal infection rhythms under light-dark cycles. Field evidence also suggests synchronized infection of Prochlorococcus by cyanophages in the oceans, which should result in progressive expression of cyanophage early, middle, and late genes. However, distinct temporal expression patterns have not been observed in cyanophage field populations.

RESULTS

In this study, we reanalyzed a previous metatranscriptomic dataset collected in the North Pacific Subtropical Gyre. In this dataset, it was previously shown that aggregate transcripts from cyanophage scaffolds display diurnal transcriptional rhythms with transcript abundances decreasing at night. By mapping metatranscriptomic reads to individual viral genes, we identified periodically expressed genes from putative viruses infecting the cyanobacteria Prochlorococcus and Synechococcus, heterotrophic bacteria, and algae. Of the 41 cyanophage genes, 35 were from cyanomyoviruses. We grouped the periodically expressed cyanomyovirus genes into early, middle, and late genes based on the conserved temporal expression patterns of their orthologs in cyanomyovirus laboratory cultures. We found that the peak expression times of late genes in cyanophage field populations were significantly later than those of early and middle genes, which were similar to the temporal expression patterns of synchronized cyanophage laboratory cultures.

CONCLUSIONS

The significantly later peak expression times of late genes in cyanomyovirus field populations suggest that cyanophage infection of Prochlorococcus is synchronized in the North Pacific Subtropical Gyre. The night-time peak expression of late genes also suggests synchronized lysis of Prochlorococcus at night, which might result in synchronized release of dissolved organic matter to the marine food web. Video abstract.

Isolation and complete genome sequence of a novel cyanophage, S-B05, infecting an estuarine Synechococcus strain: insights into environmental adaptation

A new cyanophage, S-B05, infecting a phycoerythrin-enriched (PE-type) Synechococcus strain was isolated by the liquid infection method, and its morphology and genetic features were examined. Phylogenetic analysis and morphological observation confirmed that S-B05 belongs to the family Myoviridae of the order Caudovirales. Its genome was fully sequenced, and found to be 208,857 bp in length with a G + C content of 39.9%. It contained 280 potential open reading frames and 123 conserved domains. Ninety-eight functional genes responsible for cyanophage structuring and packaging, DNA replication and regulation, and photosynthesis were identified, as well as genes encoding 172 hypothetical proteins. The genome of S-B05 is most similar to that of Prochlorococcus phage P-TIM68. Homologues of open reading frames of S-B05 can be found in various marine environments, as revealed by comparison of the S-B05 genome sequence to sequences in marine viral metagenomic databases. The presence of auxiliary metabolic genes (AMGs) related to photosynthesis, carbon metabolism, and phosphorus assimilation, as well as the phylogenetic relationships based on AMGs and the complete genome sequence, reflect the phage-host interaction mechanism or the specific adaptation strategy of the host to environmental conditions. The genome sequence information reported here will provide an important basis for further study of the adaptive evolution and ecological role of cyanophages and their hosts in the marine environment.

Complete Genome Sequences for Two Myoviridae Strains Infecting Cyanobacteria in a Subarctic Lake

We isolated two closely related strains that belong to the Myoviridae family and infect cyanobacteria in a shallow subarctic rock basin lake. Their host was identified as a member of the Synechococcus-Cyanobium complex. Sequenced genomes of the two phages were 244,930 bp and 243,633 bp. We describe their annotation and highlight some noteworthy features.

We isolated two closely related strains that belong to the Myoviridae family and infect cyanobacteria in a shallow subarctic rock basin lake. Their host was identified as a member of the Synechococcus-Cyanobium complex. Sequenced genomes of the two phages were 244,930 bp and 243,633 bp. We describe their annotation and highlight some noteworthy features.

Metagenome Mining Reveals Hidden Genomic Diversity of Pelagimyophages in Aquatic Environments

SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

The SAR11 clade is one of the most abundant bacterioplankton groups in surface waters of most of the oceans and lakes. However, only 15 SAR11 phages have been isolated thus far, and only one of them belongs to the Myoviridae family (pelagimyophages). Here, we have analyzed 26 sequences of myophages that putatively infect the SAR11 clade. They have been retrieved by mining ca. 45 Gbp aquatic assembled cellular metagenomes and viromes. Most of the myophages were obtained from the cellular fraction (0.2 μm), indicating a bias against this type of virus in viromes. We have found the first myophages that putatively infect Candidatus Fonsibacter (freshwater SAR11) and another group putatively infecting bathypelagic SAR11 phylogroup Ic. The genomes have similar sizes and maintain overall synteny in spite of low average nucleotide identity values, revealing high similarity to marine cyanomyophages. Pelagimyophages recruited metagenomic reads widely from several locations but always much more from cellular metagenomes than from viromes, opposite to what happens with pelagipodophages. Comparing the genomes resulted in the identification of a hypervariable island that is related to host recognition. Interestingly, some genes in these islands could be related to host cell wall synthesis and coinfection avoidance. A cluster of curli-related proteins was widespread among the genomes, although its function is unclear.

IMPORTANCE SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

A novel jumbo Tenacibaculum maritimum lytic phage with head-fiber-like appendages

A novel jumbo bacteriophage (myovirus) is described. The lytic phage of Tenacibaculum maritimum, which is the etiological agent of tenacibaculosis in a variety of farmed marine fish worldwide, was plaque-isolated from seawater around a fish aquaculture field in Japan. The phage had an isometric head 110-120 nm in diameter, from which several 50- to 100-nm-long flexible fiber-like appendages emanate, and a 150-nm-long rigid contractile tail. The full genomes of the two representative phages (PTm1 and PTm5) were 224,680 and 226,876 bp long, respectively, both with 29.7% GC content, and the number of predicted open reading frames (ORFs) was 308 and 306, respectively. The average nucleotide sequence identity between PTm1 and PTm5 was 99.95%, indicating they are quite similar to each other. A genetic relationship was found in 15.0-16.6% of the predicted ORFs among the T. maritimum phages PTm1 and PTm5, the Tenacibaculum spp. phage pT24, and the Sphingomonas paucimobilis phage PAU. Phylogenetic analysis based on the terminase large subunit genes revealed that these four phages (PTm1, PTm5, pT24 and PAU) are more closely related than the other 10 jumbo myoviruses that have similar genome sizes. Transmission electron microscopy observations suggest that the head fibers of the T. maritimum phage function as tentacles to search and recognize the host cell surface to facilitate infection.

Metabolic and biogeochemical consequences of viral infection in aquatic ecosystems

Ecosystems are controlled by 'bottom-up' (resources) and 'top-down' (predation) forces. Viral infection is now recognized as a ubiquitous top-down control of microbial growth across ecosystems but, at the same time, cell death by viral predation influences, and is influenced by, resource availability. In this Review, we discuss recent advances in understanding the biogeochemical impact of viruses, focusing on how metabolic reprogramming of host cells during lytic viral infection alters the flow of energy and nutrients in aquatic ecosystems. Our synthesis revealed several emerging themes. First, viral infection transforms host metabolism, in part through virus-encoded metabolic genes; the functions performed by these genes appear to alleviate energetic and biosynthetic bottlenecks to viral production. Second, viral infection depends on the physiological state of the host cell and on environmental conditions, which are challenging to replicate in the laboratory. Last, metabolic reprogramming of infected cells and viral lysis alter nutrient cycling and carbon export in the oceans, although the net impacts remain uncertain. This Review highlights the need for understanding viral infection dynamics in realistic physiological and environmental contexts to better predict their biogeochemical consequences.

Metabolic Modeling Elucidates the Transactions in the Rumen Microbiome and the Shifts Upon Virome Interactions

The complex microbial ecosystem within the bovine rumen plays a crucial role in host nutrition, health, and environmental impact. However, little is known about the interactions between the functional entities within the system, which dictates the community structure and functional dynamics and host physiology. With the advancements in high-throughput sequencing and mathematical modeling, in silico genome-scale metabolic analysis promises to expand our understanding of the metabolic interplay in the community. In an attempt to understand the interactions between microbial species and the phages inside rumen, a genome-scale metabolic modeling approach was utilized by using key members in the rumen microbiome (a bacteroidete, a firmicute, and an archaeon) and the viral phages associated with them. Individual microbial host models were integrated into a community model using multi-level mathematical frameworks. An elaborate and heuristics-based computational procedure was employed to predict previously unknown interactions involving the transfer of fatty acids, vitamins, coenzymes, amino acids, and sugars among the community members. While some of these interactions could be inferred by the available multi-omic datasets, our proposed method provides a systemic understanding of why the interactions occur and how these affect the dynamics in a complex microbial ecosystem. To elucidate the functional role of the virome on the microbiome, local alignment search was used to identify the metabolic functions of the viruses associated with the hosts. The incorporation of these functions demonstrated the role of viral auxiliary metabolic genes in relaxing the metabolic bottlenecks in the microbial hosts and complementing the inter-species interactions. Finally, a comparative statistical analysis of different biologically significant community fitness criteria identified the variation in flux space and robustness of metabolic capacities of the community members. Our elucidation of metabolite exchange among the three members of the rumen microbiome shows how their genomic differences and interactions with the viral strains shape up a highly sophisticated metabolic interplay and explains how such interactions across kingdoms can cause metabolic and compositional shifts in the community and affect the health, nutrition, and pathophysiology of the ruminant animal.

Prevalence of viral photosynthesis genes along a freshwater to saltwater transect in Southeast USA

Bacteriophages encode host-acquired functional genes known as auxiliary metabolic genes (AMGs). Photosynthesis AMGs are commonly found in marine cyanobacteria-infecting Myoviridae and Podoviridae cyanophages, but their ecology remains understudied in freshwater environments. To advance knowledge of this issue, we analysed viral metagenomes collected in the summertime for four years from five lakes and two estuarine locations interconnected by the Chattahoochee River, Southeast USA. Sequences representing ten different AMGs were recovered and found to be prevalent in all sites. Most freshwater AMGs were 10-fold less abundant than estuarine and marine AMGs and were encoded by novel Myoviridae and Podoviridae cyanophage genera. Notably, several of the corresponding viral genomes showed endemism to a specific province along the river. This translated into psbA gene phylogenetic clustering patterns that matched a marine vs. freshwater origin indicating that psbA may serve as a robust classification and source-tracking biomarker. Genomes classified in a novel viral lineage represented by isolate S-EIVl contained psbA, which is unprecedented for this lineage. Collectively, our findings indicated that the acquisition of photosynthesis AMGs is a widespread strategy used by cyanophages in aquatic ecosystems, and further indicated the existence of viral provinces in which certain viral species and/or genotypes are locally abundant.

Rates of Molecular Evolution in a Marine Synechococcus Phage Lineage

Cyanophages are characterized by vast genomic diversity and the formation of stable ecotypes over time. The evolution of phage diversity includes vertical processes, such as mutation, and horizontal processes, such as recombination and gene transfer. Here, we study the contribution of vertical and horizontal processes to short-term evolution of marine cyanophages. Analyzing time series data of Synechococcus-infecting Myoviridae ecotypes spanning up to 17 years, we found a high contribution of recombination relative to mutation (r/m) in all ecotypes. Additionally, we found a molecular clock of substitution and recombination in one ecotype, RIM8. The estimated RIM8 evolutionary rates are 2.2 genome-wide substitutions per year (1.275 × 10⁻⁵ substitutions/site/year) and 29 genome-wide nucleotide alterations due to recombination per year. We found 26 variable protein families, of which only two families have a predicted functional annotation, suggesting that they are auxiliary metabolic genes with bacterial homologs. A comparison of our rate estimates to other phage evolutionary rate estimates in the literature reveals a negative correlation of phage substitution rates with their genome size. A comparison to evolutionary rates in bacterial organisms further shows that phages have high rates of mutation and recombination compared to their bacterial hosts. We conclude that the increased recombination rate in phages likely contributes to their vast genomic diversity.

Nitrogen sourcing during viral infection of marine cyanobacteria

The building blocks of a virus derived from de novo biosynthesis during infection and/or catabolism of preexisting host cell biomass, and the relative contribution of these 2 sources has important consequences for understanding viral biogeochemistry. We determined the uptake of extracellular nitrogen (N) and its biosynthetic incorporation into both virus and host proteins using an isotope-labeling proteomics approach in a model marine cyanobacterium Synechococcus WH8102 infected by a lytic cyanophage S-SM1. By supplying dissolved N as ¹⁵N postinfection, we found that proteins in progeny phage particles were composed of up to 41% extracellularly derived N, while proteins of the infected host cell showed almost no isotope incorporation, demonstrating that de novo amino acid synthesis continues during infection and contributes specifically and substantially to phage replication. The source of N for phage protein synthesis shifted over the course of infection from mostly host derived in the early stages to more medium derived later on. We show that the photosystem II reaction center proteins D1 and D2, which are auxiliary metabolic genes (AMGs) in the S-SM1 genome, are made de novo during infection in an apparently light-dependent manner. We also identified a small set of host proteins that continue to be produced during infection; the majority are homologs of AMGs in S-SM1 or other viruses, suggesting selective continuation of host protein production during infection. The continued acquisition of nutrients by the infected cell and their utilization for phage replication are significant for both evolution and biogeochemical impact of viruses.

Lateral Gene Transfer Shapes the Distribution of RuBisCO among Candidate Phyla Radiation Bacteria and DPANN Archaea

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is considered to be the most abundant enzyme on Earth. Despite this, its full diversity and distribution across the domains of life remain to be determined. Here, we leverage a large set of bacterial, archaeal, and viral genomes recovered from the environment to expand our understanding of existing RuBisCO diversity and the evolutionary processes responsible for its distribution. Specifically, we report a new type of RuBisCO present in Candidate Phyla Radiation (CPR) bacteria that is related to the archaeal Form III enzyme and contains the amino acid residues necessary for carboxylase activity. Genome-level metabolic analyses supported the inference that these RuBisCO function in a CO₂-incorporating pathway that consumes nucleotides. Importantly, some Gottesmanbacteria (CPR) also encode a phosphoribulokinase that may augment carbon metabolism through a partial Calvin–Benson–Bassham cycle. Based on the scattered distribution of RuBisCO and its discordant evolutionary history, we conclude that this enzyme has been extensively laterally transferred across the CPR bacteria and DPANN archaea. We also report RuBisCO-like proteins in phage genomes from diverse environments. These sequences cluster with proteins in the Beckwithbacteria (CPR), implicating phage as a possible mechanism of RuBisCO transfer. Finally, we synthesize our metabolic and evolutionary analyses to suggest that lateral gene transfer of RuBisCO may have facilitated major shifts in carbon metabolism in several important bacterial and archaeal lineages.

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.

The Environment and Cyanophage Diversity: Insights From Environmental Sequencing of DNA Polymerase

Globally distributed and abundant cyanophages in the family Myoviridae have dsDNA genomes with variable gene content, including host-derived auxiliary metabolic genes (AMGs) that potentially can facilitate viral replication. However, it is not well understood how this variation in gene content interacts with environmental variables to shape cyanomyovirus communities. This project correlated the genetic repertoire of cyanomyoviruses with their phyologeny, and investigated cyanomyovirus ecotype distribution as a function of environmental conditions across locations and seasons. Reference cyanomyovirus genomes were compared for their overlap in gene content to infer phyologenetic distances, and these distances were compared to distances calculated based on DNA polymerase (gp43) gene sequences. In turn, gp43 partial gene sequences amplified from natural cyanophage communities were used to describe cyanomyovirus community composition and to assess the relationship between environmental variables. The results showed the following: (1) DNA polymerase gene phylogeny generally correlated with the similarity in gene content among reference cyanomyoviruses, and thus can be used to describe environmental cyanomyovirus communities; (2) spatial and seasonal patterns in cyanomyovirus communities were related to environmental variables; (3) salinity and temperature, combined with nutrient concentration were predictors of cyanomyovirus richness, diversity and community composition. This study shows that environmental variables shape viral communities by drawing on a diverse seed bank of viral genotypes. From these results it is evident that that viral ecotypes with their corresponding genetic repertoires underlie selection pressures. However, the mechanisms involved in selecting for specific viral genotypes remain to be fully understood.

Survey of the bacteriophage phoH gene in wetland sediments in northeast China

PhoH is a host-derived auxiliary metabolic gene that can be used as a new biomarker for surveying phage diversity in marine and paddy waters. However, the applicability of this gene in other environments has not been addressed. In this paper, we surveyed the phoH gene in four wetland sediments in northeast China. DNA was extracted directly from sediments and used for PCR amplification with the degenerate primers vPhoHf and vPhoHr. In total, 44 and 58 phoH sequences were identified as belonging to bacteria and phages, respectively, suggesting that this primer set is not highly specific to the phage phoH gene. A BLASTp search showed that the 58 phage phoH sequences had the highest identity to the known viral sequences, ranging from 48% to 100%. Phylogenetic analysis showed that all phage sequences from wetlands distributed into the previously designated Groups 2, 3, 4 and 6. In addition, two new subgroups, Groups 2c and 4c, which contained sequences exclusively from wetlands, were detected in this study. Nonmetric multidimensional scaling analysis showed that the phage phoH assemblage from a coastal wetland was similar to that in marine environments, while the phage phoH assemblage from a lake wetland was similar to that in paddy waters. These findings indicated that different types of wetlands had distinct phage phoH compositions.

Soil Viruses Are Underexplored Players in Ecosystem Carbon Processing

Rapidly thawing permafrost harbors ∼30 to 50% of global soil carbon, and the fate of this carbon remains unknown. Microorganisms will play a central role in its fate, and their viruses could modulate that impact via induced mortality and metabolic controls. Because of the challenges of recovering viruses from soils, little is known about soil viruses or their role(s) in microbial biogeochemical cycling. Here, we describe 53 viral populations (viral operational taxonomic units [vOTUs]) recovered from seven quantitatively derived (i.e., not multiple-displacement-amplified) viral-particle metagenomes (viromes) along a permafrost thaw gradient at the Stordalen Mire field site in northern Sweden. Only 15% of these vOTUs had genetic similarity to publicly available viruses in the RefSeq database, and ∼30% of the genes could be annotated, supporting the concept of soils as reservoirs of substantial undescribed viral genetic diversity. The vOTUs exhibited distinct ecology, with different distributions along the thaw gradient habitats, and a shift from soil-virus-like assemblages in the dry palsas to aquatic-virus-like assemblages in the inundated fen. Seventeen vOTUs were linked to microbial hosts (in silico), implicating viruses in infecting abundant microbial lineages from Acidobacteria, Verrucomicrobia, and Deltaproteobacteria, including those encoding key biogeochemical functions such as organic matter degradation. Thirty auxiliary metabolic genes (AMGs) were identified and suggested virus-mediated modulation of central carbon metabolism, soil organic matter degradation, polysaccharide binding, and regulation of sporulation. Together, these findings suggest that these soil viruses have distinct ecology, impact host-mediated biogeochemistry, and likely impact ecosystem function in the rapidly changing Arctic. IMPORTANCE This work is part of a 10-year project to examine thawing permafrost peatlands and is the first virome-particle-based approach to characterize viruses in these systems. This method yielded >2-fold-more viral populations (vOTUs) per gigabase of metagenome than vOTUs derived from bulk-soil metagenomes from the same site (J. B. Emerson, S. Roux, J. R. Brum, B. Bolduc, et al., Nat Microbiol 3:870-880, 2018, https://doi.org/10.1038/s41564-018-0190-y). We compared the ecology of the recovered vOTUs along a permafrost thaw gradient and found (i) habitat specificity, (ii) a shift in viral community identity from soil-like to aquatic-like viruses, (iii) infection of dominant microbial hosts, and (iv) carriage of host metabolic genes. These vOTUs can impact ecosystem carbon processing via top-down (inferred from lysing dominant microbial hosts) and bottom-up (inferred from carriage of auxiliary metabolic genes) controls. This work serves as a foundation which future studies can build upon to increase our understanding of the soil virosphere and how viruses affect soil ecosystem services.

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.

Energy limitation of cyanophage development: implications for marine carbon cycling

Marine cyanobacteria are responsible for ~25% of the fixed carbon that enters the ocean biosphere. It is thought that abundant co-occurring viruses play an important role in regulating population dynamics of cyanobacteria and thus the cycling of carbon in the oceans. Despite this, little is known about how viral infections ‘play-out’ in the environment, particularly whether infections are resource or energy limited. Photoautotrophic organisms represent an ideal model to test this since available energy is modulated by the incoming light intensity through photophosphorylation. Therefore, we exploited phototrophy of the environmentally relevant marine cyanobacterium Synechococcus and monitored growth of a cyanobacterial virus (cyanophage). We found that light intensity has a marked effect on cyanophage infection dynamics, but that this is not manifest by a change in DNA synthesis. Instead, cyanophage development appears energy limited for the synthesis of proteins required during late infection. We posit that acquisition of auxiliary metabolic genes (AMGs) involved in light-dependent photosynthetic reactions acts to overcome this limitation. We show that cyanophages actively modulate expression of these AMGs in response to light intensity and provide evidence that such regulation may be facilitated by a novel mechanism involving light-dependent splicing of a group I intron in a photosynthetic AMG. Altogether, our data offers a mechanistic link between diurnal changes in irradiance and observed community level responses in metabolism, i.e., through an irradiance-dependent, viral-induced release of dissolved organic matter (DOM).