Comparative metagenomics of microbial traits within oceanic viral communities

The complete genome of Vibrio sp. 16 unveils two circular chromosomes and a distinctive 46-kb plasmid

The entire 4.6-Mb genome of Vibrio sp. 16, encoding 4,270 genes, best matches with Vibrio rotiferianus. A 46-kb plasmid (pVDT1), alongside two circular chromosomes, showcases parAB/repB partition genes and three toxin/antitoxin systems potentially linked to phage infection.

The role of rhizosphere phages in soil health

While the One Health framework has emphasized the importance of soil microbiomes for plant and human health, one of the most diverse and abundant groups-bacterial viruses, i.e. phages-has been mostly neglected. This perspective reviews the significance of phages for plant health in rhizosphere and explores their ecological and evolutionary impacts on soil ecosystems. We first summarize our current understanding of the diversity and ecological roles of phages in soil microbiomes in terms of nutrient cycling, top-down density regulation, and pathogen suppression. We then consider how phages drive bacterial evolution in soils by promoting horizontal gene transfer, encoding auxiliary metabolic genes that increase host bacterial fitness, and selecting for phage-resistant mutants with altered ecology due to trade-offs with pathogen competitiveness and virulence. Finally, we consider challenges and avenues for phage research in soil ecosystems and how to elucidate the significance of phages for microbial ecology and evolution and soil ecosystem functioning in the future. We conclude that similar to bacteria, phages likely play important roles in connecting different One Health compartments, affecting microbiome diversity and functions in soils. From the applied perspective, phages could offer novel approaches to modulate and optimize microbial and microbe-plant interactions to enhance soil health.

The Impact of Viral Infection on the Chemistries of the Earth's Most Abundant Photosynthesizes: Metabolically Talented Aquatic Cyanobacteria

Cyanobacteria are the most abundant photosynthesizers on earth, and as such, they play a central role in marine metabolite generation, ocean nutrient cycling, and the control of planetary oxygen generation. Cyanobacteriophage infection exerts control on all of these critical processes of the planet, with the phage-ported homologs of genes linked to photosynthesis, catabolism, and secondary metabolism (marine metabolite generation). Here, we analyze the 153 fully sequenced cyanophages from the National Center for Biotechnology Information (NCBI) database and the 45 auxiliary metabolic genes (AMGs) that they deliver into their hosts. Most of these AMGs are homologs of those found within cyanobacteria and play a key role in cyanobacterial metabolism-encoding proteins involved in photosynthesis, central carbon metabolism, phosphate metabolism, methylation, and cellular regulation. A greater understanding of cyanobacteriophage infection will pave the way to a better understanding of carbon fixation and nutrient cycling, as well as provide new tools for synthetic biology and alternative approaches for the use of cyanobacteria in biotechnology and sustainable manufacturing.

Assembly and Annotation of Viral Metagenomes from Short-Read Sequencing Data

Viral metagenomics enables the detection, characterization, and quantification of viral sequences present in shotgun-sequenced datasets of purified virus-like particles and whole metagenomes. Next generation sequencing (Illumina) derived short single or paired-end read runs are a principal platform for metagenomics, and assembly of short reads allows for the identification of distinguishing viral signatures and complex genomic features for taxonomy and functional annotation. Here we describe the identification and characterization of viral genome sequences, bacteriophages, and eukaryotic viruses, from a cohort of human stool samples, using multiple methods. Following the purification of virus-like particles, sequencing, quality refinement, and genome assembly, we begin the protocol with raw short reads deposited in an open-source nucleotide archive. We highlight the use of VIBRANT, an automated computational tool for the characterization of microbial viruses and their viral community function. Finally, we also describe an alternative assembly-free option of mapping reads to established databases of reference genomes and previously characterized metagenome-assembled viral genomes.

Metagenomic and metatranscriptomic insights into sulfate-reducing bacteria in a revegetated acidic mine wasteland

The widespread occurrence of sulfate-reducing microorganisms (SRMs) in temporarily oxic/hypoxic aquatic environments indicates an intriguing possibility that SRMs can prevail in constantly oxic/hypoxic terrestrial sulfate-rich environments. However, little attention has been given to this possibility, leading to an incomplete understanding of microorganisms driving the terrestrial part of the global sulfur (S) cycle. In this study, genome-centric metagenomics and metatranscriptomics were employed to explore the diversity, metabolic potential, and gene expression profile of SRMs in a revegetated acidic mine wasteland under constantly oxic/hypoxic conditions. We recovered 16 medium- to high-quality metagenome-assembled genomes (MAGs) containing reductive dsrAB. Among them, 12 and four MAGs belonged to Acidobacteria and Deltaproteobacteria, respectively, harboring three new SRM genera. Comparative genomic analysis based on seven high-quality MAGs (completeness >90% and contamination <10%; including six acidobacterial and one deltaproteobacterial) and genomes of three additional cultured model species showed that Acidobacteria-related SRMs had more genes encoding glycoside hydrolases, oxygen-tolerant hydrogenases, and cytochrome c oxidases than Deltaproteobacteria-related SRMs. The opposite pattern was observed for genes encoding superoxide reductases and thioredoxin peroxidases. Using VirSorter, viral genome sequences were found in five of the 16 MAGs and in all three cultured model species. These prophages encoded enzymes involved in glycoside hydrolysis and antioxidation in their hosts. Moreover, metatranscriptomic analysis revealed that 15 of the 16 SRMs reported here were active in situ. An acidobacterial MAG containing a prophage dominated the SRM transcripts, expressing a large number of genes involved in its response to oxidative stress and competition for organic matter.

Auxiliary Metabolic Gene Functions in Pelagic and Benthic Viruses of the Baltic Sea

Marine microbial communities are facing various ecosystem fluctuations (e.g., temperature, organic matter concentration, salinity, or redox regimes) and thus have to be highly adaptive. This might be supported by the acquisition of auxiliary metabolic genes (AMGs) originating from virus infections. Marine bacteriophages frequently contain AMGs, which allow them to augment their host’s metabolism or enhance virus fitness. These genes encode proteins for the same metabolic functions as their highly similar host homologs. In the present study, we analyzed the diversity, distribution, and composition of marine viruses, focusing on AMGs to identify their putative ecologic role. We analyzed viruses and assemblies of 212 publicly available metagenomes obtained from sediment and water samples across the Baltic Sea. In general, the virus composition in both compartments differed compositionally. While the predominant viral lifestyle was found to be lytic, lysogeny was more prevalent in sediments than in the pelagic samples. The highest proportion of AMGs was identified in the genomes of Myoviridae. Overall, the most abundantly occurring AMGs are encoded for functions that protect viruses from degradation by their hosts, such as methylases. Additionally, some detected AMGs are known to be involved in photosynthesis, 7-cyano-7-deazaguanine synthesis, and cobalamin biosynthesis among other functions. Several AMGs that were identified in this study were previously detected in a large-scale analysis including metagenomes from various origins, i.e., different marine sites, wastewater, and the human gut. This supports the theory of globally conserved core AMGs that are spread over virus genomes, regardless of host or environment.

Isolation and Characterization of a Novel Cyanophage Encoding Multiple Auxiliary Metabolic Genes

As significant drivers of cyanobacteria mortality, cyanophages have been known to regulate the population dynamics, metabolic activities, and community structure of this most important marine autotrophic picoplankton and, therefore, influence the global primary production and biogeochemical cycle in aquatic ecosystems. In the present study, a lytic Synechococcus phage, namely S-SZBM1, was isolated and identified. Cyanophage S-SZBM1 has a double-stranded DNA genome of 177,834 bp with a G+C content of 43.31% and contains a total of 218 predicted ORFs and six tRNA genes. Phylogenetic analysis and nucleotide-based intergenomic similarity suggested that cyanophage S-SZBM1 belongs to a new genus under the family Kyanoviridae. A variety of auxiliary metabolic genes (AMGs) that have been proved or speculated to relate to photosynthesis, carbon metabolism, nucleotide synthesis and metabolism, cell protection, and other cell metabolism were identified in cyanophage S-SZBM1 genome and may affect host processes during infection. In addition, 24 of 32 predicted structural proteins were identified by a high-throughput proteome analysis which were potentially involved in the assembly processes of virion. The genomic and proteomic analysis features of cyanophage S-SZBM1 offer a valuable insight into the interactions between cyanophages and their hosts during infection.

Viral community analysis in a marine oxygen minimum zone indicates increased potential for viral manipulation of microbial physiological state

Microbial communities in oxygen minimum zones (OMZs) are known to have significant impacts on global biogeochemical cycles, but viral influence on microbial processes in these regions are much less studied. Here we provide baseline ecological patterns using microscopy and viral metagenomics from the Eastern Tropical North Pacific (ETNP) OMZ region that enhance our understanding of viruses in these climate-critical systems. While extracellular viral abundance decreased below the oxycline, viral diversity and lytic infection frequency remained high within the OMZ, demonstrating that viral influences on microbial communities were still substantial without the detectable presence of oxygen. Viral community composition was strongly related to oxygen concentration, with viral populations in low-oxygen portions of the water column being distinct from their surface layer counterparts. However, this divergence was not accompanied by the expected differences in viral-encoded auxiliary metabolic genes (AMGs) relating to nitrogen and sulfur metabolisms that are known to be performed by microbial communities in these low-oxygen and anoxic regions. Instead, several abundant AMGs were identified in the oxycline and OMZ that may modulate host responses to low-oxygen stress. We hypothesize that this is due to selection for viral-encoded genes that influence host survivability rather than modulating host metabolic reactions within the ETNP OMZ. Together, this study shows that viruses are not only diverse throughout the water column in the ETNP, including the OMZ, but their infection of microorganisms has the potential to alter host physiological state within these biogeochemically important regions of the ocean.

Distribution of soil viruses across China and their potential role in phosphorous metabolism

BACKGROUND

Viruses are the most abundant biological entities on the planet and drive biogeochemical cycling on a global scale. Our understanding of biogeography of soil viruses and their ecological functions lags significantly behind that of Bacteria and Fungi. Here, a viromic approach was used to investigate the distribution and ecological functions of viruses from 19 soils across China.

RESULTS

Soil viral community were clustered more significantly by geographical location than type of soil (agricultural and natural). Three clusters of viral communities were identified from North, Southeast and Southwest regions; these clusters differentiated using taxonomic composition and were mainly driven by geographic location and climate factors. A total of 972 viral populations (vOTUs) were detected spanning 23 viral families from the 19 viromes. Phylogenetic analyses of the phoH gene showed a remarkable diversity and the distribution of viral phoH genes was more dependent on the environment. Notably, five proteins involved in phosphorus (P) metabolism-related nucleotide synthesis functions, including dUTPase, MazG, PhoH, Thymidylate synthase complementing protein (Thy1), and Ribonucleoside reductase (RNR), were mainly identified in agricultural soils.

CONCLUSIONS

The present work revealed that soil viral communities were distributed across China according to geographical location and climate factors. In addition, P metabolism genes encoded by these viruses probably drive the synthesis of nucleotides for their own genomes inside bacterial hosts, thereby affecting P cycling in the soil ecosystems.

Hostile Takeover: How Viruses Reprogram Prokaryotic Metabolism

To reproduce, prokaryotic viruses must hijack the cellular machinery of their hosts and redirect it toward the production of viral particles. While takeover of the host replication and protein synthesis apparatus has long been considered an essential feature of infection, recent studies indicate that extensive reprogramming of host primary metabolism is a widespread phenomenon among prokaryotic viruses that is required to fulfill the biosynthetic needs of virion production. In this review we provide an overview of the most significant recent findings regarding virus-induced reprogramming of prokaryotic metabolism and suggest how quantitative systems biology approaches may be used to provide a holistic understanding of metabolic remodeling during lytic viral infection. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Hybrid assembly of an agricultural slurry virome reveals a diverse and stable community with the potential to alter the metabolism and virulence of veterinary pathogens

BACKGROUND

Viruses are the most abundant biological entities on Earth, known to be crucial components of microbial ecosystems. However, there is little information on the viral community within agricultural waste. There are currently ~ 2.7 million dairy cattle in the UK producing 7-8% of their own bodyweight in manure daily, and 28 million tonnes annually. To avoid pollution of UK freshwaters, manure must be stored and spread in accordance with guidelines set by DEFRA. Manures are used as fertiliser, and widely spread over crop fields, yet little is known about their microbial composition. We analysed the virome of agricultural slurry over a 5-month period using short and long-read sequencing.

RESULTS

Hybrid sequencing uncovered more high-quality viral genomes than long or short-reads alone; yielding 7682 vOTUs, 174 of which were complete viral genomes. The slurry virome was highly diverse and dominated by lytic bacteriophage, the majority of which represent novel genera (~ 98%). Despite constant influx and efflux of slurry, the composition and diversity of the slurry virome was extremely stable over time, with 55% of vOTUs detected in all samples over a 5-month period. Functional annotation revealed a diverse and abundant range of auxiliary metabolic genes and novel features present in the community, including the agriculturally relevant virulence factor VapE, which was widely distributed across different phage genera that were predicted to infect several hosts. Furthermore, we identified an abundance of phage-encoded diversity-generating retroelements, which were previously thought to be rare on lytic viral genomes. Additionally, we identified a group of crAssphages, including lineages that were previously thought only to be found in the human gut.

CONCLUSIONS

The cattle slurry virome is complex, diverse and dominated by novel genera, many of which are not recovered using long or short-reads alone. Phages were found to encode a wide range of AMGs that are not constrained to particular groups or predicted hosts, including virulence determinants and putative ARGs. The application of agricultural slurry to land may therefore be a driver of bacterial virulence and antimicrobial resistance in the environment. Video abstract.

Treatments of porcine fecal samples affect high-throughput virome sequencing results

The porcine enteric microbiota is currently extensively studied, taking advantage of developments in high-throughput sequencing technologies. However, the viral part of the microbiota, the virome, is being lightly explored, and the impact of the pretreatments used before sequencing the viruses is barely considered. In this study, the impacts of filtration, RNase and DNase treatments on virus reads recovery and diversity after sequencing on a MiSeq platform were assessed on fecal samples individually taken at <3, 5, 12 and 20 weeks from two piglets. None of the four pretreatment series affected the virus read averages or influenced diversity, but the samples with the higher proportion of reads corresponding to an entry in the "nt" database were those receiving the least number of pretreatments. The enzymatic pretreatments affected the detection of the single-stranded RNA viruses of Aichivirus C, porcine astrovirus, Sapovirus and posavirus, which is worrisome, as these viruses can be involved in swine diarrhea. If enzymatic pretreatments are used when sequencing using a high-throughput method, it may impact single-stranded RNA virus recovery, but not the overall virome diversity. Therefore, filtrated samples may be the better option, reducing the amount of bacterial genetic material while preserving the virus reads.

Making sense of virus size and the tradeoffs shaping viral fitness

Viruses span an impressive size range, with genome length varying a thousandfold and virion volume nearly a millionfold. For cellular organisms the scaling of traits with size is a pervasive influence on ecological processes, but whether size plays a central role in viral ecology is unknown. Here, we focus on viruses of aquatic unicellular organisms, which exhibit the greatest known range of virus size. We outline hypotheses within a quantitative framework, and analyse data where available, to consider how size affects the primary components of viral fitness. We argue that larger viruses have fewer offspring per infection and slower contact rates with host cells, but a larger genome tends to increase infection efficiency, broaden host range, and potentially increase attachment success and decrease decay rate. These countervailing selective pressures may explain why a breadth of sizes exist and even coexist when infecting the same host populations. Oligotrophic ecosystems may be enriched in "giant" viruses, because environments with resource-limited phagotrophs at low concentrations may select for broader host range, better control of host metabolism, lower decay rate and a physical size that mimics bacterial prey. Finally, we describe where further research is needed to understand the ecology and evolution of viral size diversity.

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.

Viral elements and their potential influence on microbial processes along the permanently stratified Cariaco Basin redoxcline

Little is known about viruses in oxygen-deficient water columns (ODWCs). In surface ocean waters, viruses are known to act as gene vectors among susceptible hosts. Some of these genes may have metabolic functions and are thus termed auxiliary metabolic genes (AMGs). AMGs introduced to new hosts by viruses can enhance viral replication and/or potentially affect biogeochemical cycles by modulating key microbial pathways. Here we identify 748 viral populations that cluster into 94 genera along a vertical geochemical gradient in the Cariaco Basin, a permanently stratified and euxinic ocean basin. The viral communities in this ODWC appear to be relatively novel as 80 of these viral genera contained no reference viral sequences, likely due to the isolation and unique features of this system. We identify viral elements that encode AMGs implicated in distinctive processes, such as sulfur cycling, acetate fermentation, signal transduction, [Fe–S] formation, and N-glycosylation. These AMG-encoding viruses include two putative Mu-like viruses, and viral-like regions that may constitute degraded prophages that have been modified by transposable elements. Our results provide an insight into the ecological and biogeochemical impact of viruses oxygen-depleted and euxinic habitats.

A network-based integrated framework for predicting virus-prokaryote interactions

Metagenomic sequencing has greatly enhanced the discovery of viral genomic sequences; however, it remains challenging to identify the host(s) of these new viruses. We developed VirHostMatcher-Net, a flexible, network-based, Markov random field framework for predicting virus–prokaryote interactions using multiple, integrated features: CRISPR sequences and alignment-free similarity measures (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$s_2^*$\end{document} and WIsH). Evaluation of this method on a benchmark set of 1462 known virus–prokaryote pairs yielded host prediction accuracy of 59% and 86% at the genus and phylum levels, representing 16–27% and 6–10% improvement, respectively, over previous single-feature prediction approaches. We applied our host prediction tool to crAssphage, a human gut phage, and two metagenomic virus datasets: marine viruses and viral contigs recovered from globally distributed, diverse habitats. Host predictions were frequently consistent with those of previous studies, but more importantly, this new tool made many more confident predictions than previous tools, up to nearly 3-fold more (n > 27 000), greatly expanding the diversity of known virus–host interactions.

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.

Diversities and potential biogeochemical impacts of mangrove soil viruses

BACKGROUND

Mangroves are ecologically and economically important forests of the tropics. As one of the most carbon-rich biomes, mangroves account for 11% of the total input of terrestrial carbon into oceans. Although viruses are considered to significantly influence local and global biogeochemical cycles, little information is available regarding the community structure, genetic diversity and ecological roles of viruses in mangrove ecosystems.

METHODS

Here, we utilised viral metagenomics sequencing and virome-specific bioinformatics tools to study viral communities in six mangrove soil samples collected from different mangrove habitats in Southern China.

RESULTS

Mangrove soil viruses were found to be largely uncharacterised. Phylogenetic analyses of the major viral groups demonstrated extensive diversity and previously unknown viral clades and suggested that global mangrove viral communities possibly comprise evolutionarily close genotypes. Comparative analysis of viral genotypes revealed that mangrove soil viromes are mainly affected by marine waters, with less influence coming from freshwaters. Notably, we identified abundant auxiliary carbohydrate-active enzyme (CAZyme) genes from mangrove viruses, most of which participate in biolysis of complex polysaccharides, which are abundant in mangrove soils and organism debris. Host prediction results showed that viral CAZyme genes are diverse and probably widespread in mangrove soil phages infecting diverse bacteria of different phyla.

CONCLUSIONS

Our results showed that mangrove viruses are diverse and probably directly manipulate carbon cycling by participating in biomass recycling of complex polysaccharides, providing the knowledge essential in revealing the ecological roles of viruses in mangrove ecosystems.

Numerous cultivated and uncultivated viruses encode ribosomal proteins

Viruses modulate ecosystems by directly altering host metabolisms through auxiliary metabolic genes. However, viral genomes are not known to encode the core components of translation machinery, such as ribosomal proteins (RPs). Here, using reference genomes and global-scale viral metagenomic datasets, we identify 14 different RPs across viral genomes arising from cultivated viral isolates and metagenome-assembled viruses. Viruses tend to encode dynamic RPs, easily exchangeable between ribosomes, suggesting these proteins can replace cellular versions in host ribosomes. Functional assays confirm that the two most common virus-encoded RPs, bS21 and bL12, are incorporated into 70S ribosomes when expressed in Escherichia coli. Ecological distribution of virus-encoded RPs suggests some level of ecosystem adaptations as aquatic viruses and viruses of animal-associated bacteria are enriched for different subsets of RPs. Finally, RP genes are under purifying selection and thus likely retained an important function after being horizontally transferred into virus genomes.

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.

A giant virus infecting green algae encodes key fermentation genes

The family Mimiviridae contains uncommonly large viruses, many of which were isolated using a free-living amoeba as a host. Although the genomes of these and other mimivirids that infect marine heterokont and haptophyte protists have now been sequenced, there has yet to be a genomic investigation of a mimivirid that infects a member of the Viridiplantae lineage (green algae and land plants). Here we characterize the 668-kilobase complete genome of TetV-1, a mimivirid that infects the cosmopolitan green alga Tetraselmis (Chlorodendrophyceae). The analysis revealed genes not previously seen in viruses, such as the mannitol metabolism enzyme mannitol 1-phosphate dehydrogenase, the saccharide degradation enzyme alpha-galactosidase, and the key fermentation genes pyruvate formate-lyase and pyruvate formate-lyase activating enzyme. The TetV genome is the largest sequenced to date for a virus that infects a photosynthetic organism, and its genes reveal unprecedented mechanisms by which viruses manipulate their host's metabolism.

Dietary energy drives the dynamic response of bovine rumen viral communities

BACKGROUND

Rumen microbes play a greater role in host energy acquisition than that of gut-associated microbes in monogastric animals. Although genome-enabled advancements are providing access to the vast diversity of uncultivated microbes, our understanding of variables shaping rumen microbial communities is in its infancy. Viruses have been shown to impact microbial populations through a myriad of processes, including cell lysis and reprogramming of host metabolism. However, little is known about the processes shaping the distribution of rumen viruses or how viruses may modulate microbial-driven processes in the rumen. To this end, we investigated how rumen bacterial and viral community structure and function responded in five steers fed four randomized dietary treatments in a crossover design.

RESULTS

Total digestible nutrients (TDN), a measure of dietary energy, best explained the variation in bacterial and viral communities. Additional ecological drivers of viral communities included dietary zinc content and microbial functional diversity. Using partial least squares regression, we demonstrate significant associations between the abundances of 267 viral populations and variables driving the variation in rumen viral communities. While rumen viruses were dynamic, 14 near ubiquitous viral populations were identified, suggesting the presence of a core rumen virome largely comprised of novel viruses. Moreover, analysis of virally encoded auxiliary metabolic genes (AMGs) indicates rumen viruses have glycosidic hydrolases to potentially augment the breakdown of complex carbohydrates to increase energy production. Other AMGs identified have a role in redirecting carbon to the pentose phosphate pathway and one carbon pools by folate to boost viral replication.

CONCLUSIONS

We demonstrate that rumen bacteria and viruses have differing responses and ecological drivers to dietary perturbation. Our results show that rumen viruses have implications for understanding the structuring of the previously identified core rumen microbiota and impacting microbial metabolism through a vast array of AMGs. AMGs in the rumen appear to have consequences for microbial metabolism that are largely in congruence with the current paradigm established in marine systems. This study provides a foundation for future hypotheses regarding the dynamics of viral-mediated processes in the rumen.

Peptide deformylases from Vibrio parahaemolyticus phage and bacteria display similar deformylase activity and inhibitor binding clefts

Unexpected peptide deformylase (PDF) genes were recently retrieved in numerous marine phage genomes. While various hypotheses dealing with the occurrence of these intriguing sequences have been made, no further characterization and functional studies have been described thus far. In this study, we characterize the bacteriophage Vp16 PDF enzyme, as representative member of the newly identified C-terminally truncated viral PDFs. We show here that conditions classically used for bacterial PDFs lead to an enzyme exhibiting weak activity. Nonetheless, our integrated biophysical and biochemical approaches reveal specific effects of pH and metals on Vp16 PDF stability and activity. A novel purification protocol taking in account these data allowed strong improvement of Vp16 PDF specific activity to values similar to those of bacterial PDFs. We next show that Vp16 PDF is as sensitive to the natural inhibitor compound of PDFs, actinonin, as bacterial PDFs. Comparison of the 3D structures of Vp16 and E. coli PDFs bound to actinonin also reveals that both PDFs display identical substrate binding mode. We conclude that bacteriophage Vp16 PDF protein has functional peptide deformylase activity and we suggest that encoded phage PDFs might be important for viral fitness.

Cyanophage-encoded lipid desaturases: oceanic distribution, diversity and function

Cyanobacteria are among the most abundant photosynthetic organisms in the oceans; viruses infecting cyanobacteria (cyanophages) can alter cyanobacterial populations, and therefore affect the local food web and global biochemical cycles. These phages carry auxiliary metabolic genes (AMGs), which rewire various metabolic pathways in the infected host cell, resulting in increased phage fitness. Coping with stress resulting from photodamage appears to be a central necessity of cyanophages, yet the overall mechanism is poorly understood. Here we report a novel, widespread cyanophage AMG, encoding a fatty acid desaturase (FAD), found in two genotypes with distinct geographical distribution. FADs are capable of modulating the fluidity of the host’s membrane, a fundamental stress response in living cells. We show that both viral FAD (vFAD) families are Δ9 lipid desaturases, catalyzing the desaturation at carbon 9 in C16 fatty acid chains. In addition, we present a comprehensive fatty acid profiling for marine cyanobacteria, which suggests a unique desaturation pathway of medium- to long-chain fatty acids no longer than C16, in accordance with the vFAD activity. Our findings suggest that cyanophages are capable of fiddling with the infected host’s membranes, possibly leading to increased photoprotection and potentially enhancing viral-encoded photosynthetic proteins, resulting in a new viral metabolic network.

Metagenomic Analysis of Dairy Bacteriophages: Extraction Method and Pilot Study on Whey Samples Derived from Using Undefined and Defined Mesophilic Starter Cultures

Despite being potentially highly useful for characterizing the biodiversity of phages, metagenomic studies are currently not available for dairy bacteriophages, partly due to the lack of a standard procedure for phage extraction. We optimized an extraction method that allows the removal of the bulk protein from whey and milk samples with losses of less than 50% of spiked phages. The protocol was applied to extract phages from whey in order to test the notion that members of Lactococcus lactis 936 (now Sk1virus), P335, c2 (now C2virus) and Leuconostoc phage groups are the most frequently encountered in the dairy environment. The relative abundance and diversity of phages in eight and four whey mixtures from dairies using undefined mesophilic mixed-strain cultures containing Lactococcus lactis subsp. lactis biovar diacetylactis and Leuconostoc species (i.e., DL starter cultures) and defined cultures, respectively, were assessed. Results obtained from transmission electron microscopy and high-throughput sequence analyses revealed the dominance of Lc. lactis 936 phages (order Caudovirales, family Siphoviridae) in dairies using undefined DL starter cultures and Lc. lactis c2 phages (order Caudovirales, family Siphoviridae) in dairies using defined cultures. The 936 and Leuconostoc phages demonstrated limited diversity. Possible coinduction of temperate P335 prophages and satellite phages in one of the whey mixtures was also observed.IMPORTANCE The method optimized in this study could provide an important basis for understanding the dynamics of the phage community (abundance, development, diversity, evolution, etc.) in dairies with different sizes, locations, and production strategies. It may also enable the discovery of previously unknown phages, which is crucial for the development of rapid molecular biology-based methods for phage burden surveillance systems. The dominance of only a few phage groups in the dairy environment signifies the depth of knowledge gained over the past decades, which served as the basis for designing current phage control strategies. The presence of a correlation between phages and the type of starter cultures being used in dairies might help to improve the selection and/or design of suitable, custom, and cost-efficient phage control strategies.

The C-terminal residue of phage Vp16 PDF, the smallest peptide deformylase, acts as an offset element locking the active conformation

Prokaryotic proteins must be deformylated before the removal of their first methionine. Peptide deformylase (PDF) is indispensable and guarantees this mechanism. Recent metagenomics studies revealed new idiosyncratic PDF forms as the most abundant family of viral sequences. Little is known regarding these viral PDFs, including the capacity of the corresponding encoded proteins to ensure deformylase activity. We provide here the first evidence that viral PDFs, including the shortest PDF identified to date, Vp16 PDF, display deformylase activity in vivo, despite the absence of the key ribosome-interacting C-terminal region. Moreover, characterization of phage Vp16 PDF underscores unexpected structural and molecular features with the C-terminal Isoleucine residue significantly contributing to deformylase activity both in vitro and in vivo. This residue fully compensates for the absence of the usual long C-domain. Taken together, these data elucidate an unexpected mechanism of enzyme natural evolution and adaptation within viral sequences.

From virus isolation to metagenome generation for investigating viral diversity in deep-sea sediments

Viruses are the most abundant and, likely, one of the most diverse biological components in the oceans. By infecting their hosts, they play key roles in biogeochemical cycles and ecosystem functioning at a global scale. The ocean interior hosts most of the microbial life, and, despite deep-sea sediments represent the main repository of this component and the largest biome on Earth, viral diversity in these ecosystems remains almost completely unknown. We compared a physical-chemical procedure and a previously published sediment washing-based procedure for isolating viruses from benthic deep-sea ecosystems to generate viromes through high-throughput sequencing. The procedure based on a physical-chemical dislodgment of viral particles from the sediments, followed by vacuum filtration was much more efficient allowing us to recover >85% of the extractable viruses. By using this procedure, a high fraction of viral DNA was recovered and new viromes from different benthic deep-sea sites were generated. Such viromes were diversified in terms of both viral families and putative functions. Overall, the results presented here provide new insights for evaluating benthic deep-sea viral diversity through metagenomic analyses, and reveal that deep-sea sediments are a hot spot of novel viral genotypes and functions.

Biological action in Read-Write genome evolution

Many of the most important evolutionary variations that generated phenotypic adaptations and originated novel taxa resulted from complex cellular activities affecting genome content and expression. These activities included (i) the symbiogenetic cell merger that produced the mitochondrion-bearing ancestor of all extant eukaryotes, (ii) symbiogenetic cell mergers that produced chloroplast-bearing ancestors of photosynthetic eukaryotes, and (iii) interspecific hybridizations and genome doublings that generated new species and adaptive radiations of higher plants and animals. Adaptive variations also involved horizontal DNA transfers and natural genetic engineering by mobile DNA elements to rewire regulatory networks, such as those essential to viviparous reproduction in mammals. In the most highly evolved multicellular organisms, biological complexity scales with 'non-coding' DNA content rather than with protein-coding capacity in the genome. Coincidentally, 'non-coding' RNAs rich in repetitive mobile DNA sequences function as key regulators of complex adaptive phenotypes, such as stem cell pluripotency. The intersections of cell fusion activities, horizontal DNA transfers and natural genetic engineering of Read-Write genomes provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

A myovirus encoding both photosystem I and II proteins enhances cyclic electron flow in infected Prochlorococcus cells

Cyanobacteria are important contributors to primary production in the open oceans. Over the past decade, various photosynthesis-related genes have been found in viruses that infect cyanobacteria (cyanophages). Although photosystem II (PSII) genes are common in both cultured cyanophages and environmental samples ^1-4 , viral photosystem I (vPSI) genes have so far only been detected in environmental samples ^5,6 . Here, we have used a targeted strategy to isolate a cyanophage from the tropical Pacific Ocean that carries a PSI gene cassette with seven distinct PSI genes (psaJF, C, A, B, K, E, D) as well as two PSII genes (psbA, D). This cyanophage, P-TIM68, belongs to the T4-like myoviruses, has a prolate capsid, a long contractile tail and infects Prochlorococcus sp. strain MIT9515. Phage photosynthesis genes from both photosystems are expressed during infection, and the resultant proteins are incorporated into membranes of the infected host. Moreover, photosynthetic capacity in the cell is maintained throughout the infection cycle with enhancement of cyclic electron flow around PSI. Analysis of metagenomic data from the Tara Oceans expedition ⁷ shows that phages carrying PSI gene cassettes are abundant in the tropical Pacific Ocean, composing up to 28% of T4-like cyanomyophages. They are also present in the tropical Indian and Atlantic Oceans. P-TIM68 populations, specifically, compose on average 22% of the PSI-gene-cassette carrying phages. Our results suggest that cyanophages carrying PSI and PSII genes are likely to maintain and even manipulate photosynthesis during infection of their Prochlorococcus hosts in the tropical oceans.

Marine viruses discovered via metagenomics shed light on viral strategies throughout the oceans

Marine viruses are key drivers of host diversity, population dynamics and biogeochemical cycling and contribute to the daily flux of billions of tons of organic matter. Despite recent advancements in metagenomics, much of their biodiversity remains uncharacterized. Here we report a data set of 27,346 marine virome contigs that includes 44 complete genomes. These outnumber all currently known phage genomes in marine habitats and include members of previously uncharacterized lineages. We designed a new method for host prediction based on co-occurrence associations that reveals these viruses infect dominant members of the marine microbiome such as Prochlorococcus and Pelagibacter. A negative association between host abundance and the virus-to-host ratio supports the recently proposed Piggyback-the-Winner model of reduced phage lysis at higher host densities. An analysis of the abundance patterns of viruses throughout the oceans revealed how marine viral communities adapt to various seasonal, temperature and photic regimes according to targeted hosts and the diversity of auxiliary metabolic genes.

In silico approach to designing rational metagenomic libraries for functional studies

Background

With the development of Next Generation Sequencing technologies, the number of predicted proteins from entire (meta-) genomes has risen exponentially. While for some of these sequences protein functions can be inferred from homology, an experimental characterization is still a requirement for the determination of protein function. However, functional characterization of proteins cannot keep pace with our capabilities to generate more and more sequence data.

Results

Here, we present an approach to reduce the number of proteins from entire (meta-) genomes to a reasonably small number for further experimental characterization without loss of important information. About 6.1 million predicted proteins from the Global Ocean Sampling Expedition Metagenome project were distributed into classes based either on homology to existing hidden markov models (HMMs) of known families, or de novo by assessment of pairwise similarity. 5.1 million of these proteins could be classified in this way, yielding 18,437 families. For 4,129 protein families, which did not match existing HMMs from databases, we could create novel HMMs. For each family, we then selected a representative protein, which showed the closest homology to all other proteins in this family. We then selected representatives of four families based on their homology to known and well-characterized lipases. From these four synthesized genes, we could obtain the novel esterase/lipase GOS54, validating our approach.

Conclusions

Using an in silico approach, we were able improve the success rate of functional screening and make entire (meta-) genomes amenable for biochemical characterization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1668-y) contains supplementary material, which is available to authorized users.

The use of informativity in the development of robust viromics-based examinations

Metagenomics-based studies have provided insight into many of the complex microbial communities responsible for maintaining life on this planet. Sequencing efforts often uncover novel genetic content; this is most evident for phage communities, in which upwards of 90% of all sequences exhibit no similarity to any sequence in current data repositories. For the small fraction that can be identified, the top BLAST hit is generally posited as being representative of a viral taxon present in the sample of origin. Homology-based classification, however, can be misleading as sequence repositories capture but a small fraction of phage diversity. Furthermore, lateral gene transfer is pervasive within phage communities. As such, the presence of a particular gene may not be indicative of the presence of a particular viral species. Rather, it is just that: an indication of the presence of a specific gene. To circumvent this limitation, we have developed a new method for the analysis of viral metagenomic datasets. BLAST hits are weighted, integrating the sequence identity and length of alignments as well as a taxonomic signal, such that each gene is evaluated with respect to its information content. Through this quantifiable metric, predictions of viral community structure can be made with confidence. As a proof-of-concept, the approach presented here was implemented and applied to seven freshwater viral metagenomes. While providing a robust method for evaluating viral metagenomic data, the tool is versatile and can easily be customized to investigations of any environment or biome.

Environmental Viral Genomes Shed New Light on Virus-Host Interactions in the Ocean

Viruses are diverse and play significant ecological roles in marine ecosystems. However, our knowledge of genome-level diversity in viruses is biased toward those isolated from few culturable hosts. Here, we determined 1,352 nonredundant complete viral genomes from marine environments. Lifting the uncertainty that clouds short incomplete sequences, whole-genome-wide analysis suggests that these environmental genomes represent hundreds of putative novel viral genera. Predicted hosts include dominant groups of marine bacteria and archaea with no isolated viruses to date. Some of the viral genomes encode many functionally related enzymes, suggesting a strong selection pressure on these marine viruses to control cellular metabolisms by accumulating genes.

Metagenomics has revealed the existence of numerous uncharacterized viral lineages, which are referred to as viral “dark matter.” However, our knowledge regarding viral genomes is biased toward culturable viruses. In this study, we analyzed 1,600 (1,352 nonredundant) complete double-stranded DNA viral genomes (10 to 211 kb) assembled from 52 marine viromes. Together with 244 previously reported uncultured viral genomes, a genome-wide comparison delineated 617 genus-level operational taxonomic units (OTUs) for these environmental viral genomes (EVGs). Of these, 600 OTUs contained no representatives from known viruses, thus putatively corresponding to novel viral genera. Predicted hosts of the EVGs included major groups of marine prokaryotes, such as marine group II Euryarchaeota and SAR86, from which no viruses have been isolated to date, as well as Flavobacteriaceae and SAR116. Our analysis indicates that marine cyanophages are already well represented in genome databases and that one of the EVGs likely represents a new cyanophage lineage. Several EVGs encode many enzymes that appear to function for an efficient utilization of iron-sulfur clusters or to enhance host survival. This suggests that there is a selection pressure on these marine viruses to accumulate genes for specific viral propagation strategies. Finally, we revealed that EVGs contribute to a 4-fold increase in the recruitment of photic-zone viromes compared with the use of current reference viral genomes.

IMPORTANCE Viruses are diverse and play significant ecological roles in marine ecosystems. However, our knowledge of genome-level diversity in viruses is biased toward those isolated from few culturable hosts. Here, we determined 1,352 nonredundant complete viral genomes from marine environments. Lifting the uncertainty that clouds short incomplete sequences, whole-genome-wide analysis suggests that these environmental genomes represent hundreds of putative novel viral genera. Predicted hosts include dominant groups of marine bacteria and archaea with no isolated viruses to date. Some of the viral genomes encode many functionally related enzymes, suggesting a strong selection pressure on these marine viruses to control cellular metabolisms by accumulating genes.

Gene Expression Patterns during Light and Dark Infection of Prochlorococcus by Cyanophage

Cyanophage infecting the marine cyanobacteria Prochlorococcus and Synechococcus require light and host photosystem activity for optimal reproduction. Many cyanophages encode multiple photosynthetic electron transport (PET) proteins, which are presumed to maintain electron flow and produce ATP and NADPH for nucleotide biosynthesis and phage genome replication. However, evidence suggests phage augment NADPH production via the pentose phosphate pathway (PPP), thus calling into question the need for NADPH production by PET. Genes implicated in cyclic PET have since been identified in cyanophage genomes. It remains an open question which mode of PET, cyclic or linear, predominates in infected cyanobacteria, and thus whether the balance is towards producing ATP or NADPH. We sequenced transcriptomes of a cyanophage (P-HM2) and its host (Prochlorococcus MED4) throughout infection in the light or in the dark, and analyzed these data in the context of phage replication and metabolite measurements. Infection was robust in the light, but phage were not produced in the dark. Host gene transcripts encoding high-light inducible proteins and two terminal oxidases (plastoquinol terminal oxidase and cytochrome c oxidase)-implicated in protecting the photosynthetic membrane from light stress-were the most enriched in light but not dark infection. Among the most diminished transcripts in both light and dark infection was ferredoxin-NADP+ reductase (FNR), which uses the electron acceptor NADP+ to generate NADPH in linear photosynthesis. The phage gene for CP12, which putatively inhibits the Calvin cycle enzyme that receives NADPH from FNR, was highly expressed in light infection. Therefore, both PET production of NADPH and its consumption by carbon fixation are putatively repressed during phage infection in light. Transcriptomic evidence is thus consistent with cyclic photophosphorylation using oxygen as the terminal electron acceptor as the dominant mode of PET under infection, with ATP from PET and NADPH from the PPP producing the energy and reducing equivalents for phage nucleotide biosynthesis and replication.

Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses

Ocean microbes drive biogeochemical cycling on a global scale. However, this cycling is constrained by viruses that affect community composition, metabolic activity, and evolutionary trajectories. Owing to challenges with the sampling and cultivation of viruses, genome-level viral diversity remains poorly described and grossly understudied, with less than 1% of observed surface-ocean viruses known. Here we assemble complete genomes and large genomic fragments from both surface- and deep-ocean viruses sampled during the Tara Oceans and Malaspina research expeditions, and analyse the resulting 'global ocean virome' dataset to present a global map of abundant, double-stranded DNA viruses complete with genomic and ecological contexts. A total of 15,222 epipelagic and mesopelagic viral populations were identified, comprising 867 viral clusters (defined as approximately genus-level groups). This roughly triples the number of known ocean viral populations and doubles the number of candidate bacterial and archaeal virus genera, providing a near-complete sampling of epipelagic communities at both the population and viral-cluster level. We found that 38 of the 867 viral clusters were locally or globally abundant, together accounting for nearly half of the viral populations in any global ocean virome sample. While two-thirds of these clusters represent newly described viruses lacking any cultivated representative, most could be computationally linked to dominant, ecologically relevant microbial hosts. Moreover, we identified 243 viral-encoded auxiliary metabolic genes, of which only 95 were previously known. Deeper analyses of four of these auxiliary metabolic genes (dsrC, soxYZ, P-II (also known as glnB) and amoC) revealed that abundant viruses may directly manipulate sulfur and nitrogen cycling throughout the epipelagic ocean. This viral catalog and functional analyses provide a necessary foundation for the meaningful integration of viruses into ecosystem models where they act as key players in nutrient cycling and trophic networks.

Metabolic Genes within Cyanophage Genomes: Implications for Diversity and Evolution

Cyanophages, a group of viruses specifically infecting cyanobacteria, are genetically diverse and extensively abundant in water environments. As a result of selective pressure, cyanophages often acquire a range of metabolic genes from host genomes. The host-derived genes make a significant contribution to the ecological success of cyanophages. In this review, we summarize the host-derived metabolic genes, as well as their origin and roles in cyanophage evolution and important host metabolic pathways, such as the light-dependent reactions of photosynthesis, the pentose phosphate pathway, nutrient acquisition and nucleotide biosynthesis. We also discuss the suitability of the host-derived metabolic genes as potential diagnostic markers for the detection of genetic diversity of cyanophages in natural environments.

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

Viruses of marine cyanobacteria frequently contain auxiliary metabolic genes (AMGs) that augment host metabolism during infection, but little is known about their adaptive significance. We analyzed the distribution and genomic context of 33 AMGs across 60 cyanomyovirus genomes. Similarity in AMG content among cyanomyoviruses was only weakly correlated with phylogenetic relatedness; however, AMG content was generally conserved within the same operational taxonomic unit (OTU). A virus' AMG repertoire was also correlated with its isolation host and environment (coastal versus open ocean). A new analytical method based on shared co-linear blocks revealed that variation in the genomic location of an AMG was negatively correlated with its frequency across the genomes. We propose that rare AMGs are more frequently gained or lost as a result of fluctuating selection pressures, whereas common AMGs are associated with stable selection pressures. Finally, we describe a unique cyanomyovirus (S-CAM7) that lacks many AMGs including the photosynthesis gene psbA.

iVirus: facilitating new insights in viral ecology with software and community data sets imbedded in a cyberinfrastructure

Microbes affect nutrient and energy transformations throughout the world's ecosystems, yet they do so under viral constraints. In complex communities, viral metagenome (virome) sequencing is transforming our ability to quantify viral diversity and impacts. Although some bottlenecks, for example, few reference genomes and nonquantitative viromics, have been overcome, the void of centralized data sets and specialized tools now prevents viromics from being broadly applied to answer fundamental ecological questions. Here we present iVirus, a community resource that leverages the CyVerse cyberinfrastructure to provide access to viromic tools and data sets. The iVirus Data Commons contains both raw and processed data from 1866 samples and 73 projects derived from global ocean expeditions, as well as existing and legacy public repositories. Through the CyVerse Discovery Environment, users can interrogate these data sets using existing analytical tools (software applications known as ‘Apps') for assembly, open reading frame prediction and annotation, as well as several new Apps specifically developed for analyzing viromes. Because Apps are web based and powered by CyVerse supercomputing resources, they enable scalable analyses for a broad user base. Finally, a use-case scenario documents how to apply these advances toward new data. This growing iVirus resource should help researchers utilize viromics as yet another tool to elucidate viral roles in nature.

Marine phage genomics: the tip of the iceberg

Marine viruses are the most abundant biological entity in the oceans, the majority of which infect bacteria and are known as bacteriophages. Yet, the bulk of bacteriophages form part of the vast uncultured dark matter of the microbial biosphere. In spite of the paucity of cultured marine bacteriophages, it is known that marine bacteriophages have major impacts on microbial population structure and the biogeochemical cycling of key elements. Despite the ecological relevance of marine bacteriophages, there are relatively few isolates with complete genome sequences. This minireview focuses on knowledge gathered from these genomes put in the context of viral metagenomic data and highlights key advances in the field, particularly focusing on genome structure and auxiliary metabolic genes.

Only a tiny fraction of marine phages have been discovered, yet are known to have important roles in the ocean.

Freshwater Metaviromics and Bacteriophages: A Current Assessment of the State of the Art in Relation to Bioinformatic Challenges

Advances in bioinformatics and sequencing technologies have allowed for the analysis of complex microbial communities at an unprecedented rate. While much focus is often placed on the cellular members of these communities, viruses play a pivotal role, particularly bacteria-infecting viruses (bacteriophages); phages mediate global biogeochemical processes and drive microbial evolution through bacterial grazing and horizontal gene transfer. Despite their importance and ubiquity in nature, very little is known about the diversity and structure of viral communities. Though the need for culture-based methods for viral identification has been somewhat circumvented through metagenomic techniques, the analysis of metaviromic data is marred with many unique issues. In this review, we examine the current bioinformatic approaches for metavirome analyses and the inherent challenges facing the field as illustrated by the ongoing efforts in the exploration of freshwater phage populations.

Nothing in Evolution Makes Sense Except in the Light of Genomics: Read-Write Genome Evolution as an Active Biological Process

The 21st century genomics-based analysis of evolutionary variation reveals a number of novel features impossible to predict when Dobzhansky and other evolutionary biologists formulated the neo-Darwinian Modern Synthesis in the middle of the last century. These include three distinct realms of cell evolution; symbiogenetic fusions forming eukaryotic cells with multiple genome compartments; horizontal organelle, virus and DNA transfers; functional organization of proteins as systems of interacting domains subject to rapid evolution by exon shuffling and exonization; distributed genome networks integrated by mobile repetitive regulatory signals; and regulation of multicellular development by non-coding lncRNAs containing repetitive sequence components. Rather than single gene traits, all phenotypes involve coordinated activity by multiple interacting cell molecules. Genomes contain abundant and functional repetitive components in addition to the unique coding sequences envisaged in the early days of molecular biology. Combinatorial coding, plus the biochemical abilities cells possess to rearrange DNA molecules, constitute a powerful toolbox for adaptive genome rewriting. That is, cells possess "Read-Write Genomes" they alter by numerous biochemical processes capable of rapidly restructuring cellular DNA molecules. Rather than viewing genome evolution as a series of accidental modifications, we can now study it as a complex biological process of active self-modification.

New biosynthetic pathway for pink pigments from uncultured oceanic viruses

The pink open-chain tetrapyrrole pigment phycoerythrobilin (PEB) is employed by marine cyanobacteria, red algae and cryptophytes as a light-harvesting chromophore in phycobiliproteins. Genes encoding biosynthesis proteins for PEB have also been discovered in cyanophages, viruses that infect cyanobacteria, and mimic host pigment biosynthesis with the exception of PebS which combines the enzymatic activities of two host enzymes. In this study, we have identified novel members of the PEB biosynthetic enzyme families, heme oxygenases and ferredoxin-dependent bilin reductases. Encoding genes were found in metagenomic datasets and could be traced back to bacteriophage but not cyanophage origin. While the heme oxygenase exhibited standard activity, a new bilin reductase with highest homology to the teal pigment producing enzyme PcyA revealed PEB biosynthetic activity. Although PcyX possesses PebS-like activity both enzymes share only 9% sequence identity and likely catalyze the reaction via two independent mechanisms. Our data point towards the presence of phycobilin biosynthetic genes in phages that probably infect alphaproteobacteria and, therefore, further support a role of phycobilins outside oxygenic phototrophs.

High coverage metabolomics analysis reveals phage-specific alterations to Pseudomonas aeruginosa physiology during infection

Phage-mediated metabolic changes in bacteria are hypothesized to markedly alter global nutrient and biogeochemical cycles. Despite their theoretic importance, experimental data on the net metabolic impact of phage infection on the bacterial metabolism remains scarce. In this study, we tracked the dynamics of intracellular metabolites using untargeted high coverage metabolomics in Pseudomonas aeruginosa cells infected with lytic bacteriophages from six distinct phage genera. Analysis of the metabolomics data indicates an active interference in the host metabolism. In general, phages elicit an increase in pyrimidine and nucleotide sugar metabolism. Furthermore, clear phage-specific and infection stage-specific responses are observed, ranging from extreme metabolite depletion (for example, phage YuA) to complete reorganization of the metabolism (for example, phage phiKZ). As expected, pathways targeted by the phage-encoded auxiliary metabolic genes (AMGs) were enriched among the metabolites changing during infection. The effect on pyrimidine metabolism of phages encoding AMGs capable of host genome degradation (for example, YuA and LUZ19) was distinct from those lacking nuclease-encoding genes (for example, phiKZ), which demonstrates the link between the encoded set of AMGs of a phage and its impact on host physiology. However, a large fraction of the profound effect on host metabolism could not be attributed to the phage-encoded AMGs. We suggest a potentially crucial role for small, 'non-enzymatic' peptides in metabolism take-over and hypothesize on potential biotechnical applications for such peptides. The highly phage-specific nature of the metabolic impact emphasizes the potential importance of the 'phage diversity' parameter when studying metabolic interactions in complex communities.

Computational approaches to predict bacteriophage-host relationships

Metagenomics has changed the face of virus discovery by enabling the accurate identification of viral genome sequences without requiring isolation of the viruses. As a result, metagenomic virus discovery leaves the first and most fundamental question about any novel virus unanswered: What host does the virus infect? The diversity of the global virosphere and the volumes of data obtained in metagenomic sequencing projects demand computational tools for virus–host prediction. We focus on bacteriophages (phages, viruses that infect bacteria), the most abundant and diverse group of viruses found in environmental metagenomes. By analyzing 820 phages with annotated hosts, we review and assess the predictive power of in silico phage–host signals. Sequence homology approaches are the most effective at identifying known phage–host pairs. Compositional and abundance-based methods contain significant signal for phage–host classification, providing opportunities for analyzing the unknowns in viral metagenomes. Together, these computational approaches further our knowledge of the interactions between phages and their hosts. Importantly, we find that all reviewed signals significantly link phages to their hosts, illustrating how current knowledge and insights about the interaction mechanisms and ecology of coevolving phages and bacteria can be exploited to predict phage–host relationships, with potential relevance for medical and industrial applications.

New viruses infecting bacteria are increasingly being discovered in many environments through sequence-based explorations. To understand their role in microbial ecosystems, computational tools are indispensable to prioritize and guide experimental efforts. This review assesses and discusses a range of bioinformatic approaches to predict bacteriophage–host relationships when all that is known is their genome sequence.

Closing the gaps on the viral photosystem-I psaDCAB gene organization

Marine photosynthesis is largely driven by cyanobacteria, namely Synechococcus and Prochlorococcus. Genes encoding for photosystem (PS) I and II reaction centre proteins are found in cyanophages and are believed to increase their fitness. Two viral PSI gene arrangements are known, psaJF→C→A→B→K→E→D and psaD→C→A→B. The shared genes between these gene cassettes and their encoded proteins are distinguished by %G + C and protein sequence respectively. The data on the psaD→C→A→B gene organization were reported from only two partial gene cassettes coming from Global Ocean Sampling stations in the Pacific and Indian oceans. Now we have extended our search to 370 marine stations from six metagenomic projects. Genes corresponding to both PSI gene arrangements were detected in the Pacific, Indian and Atlantic oceans, confined to a strip along the equator (30°N and 30°S). In addition, we found that the predicted structure of the viral PsaA protein from the psaD→C→A→B organization contains a lumenal loop conserved in PsaA proteins from Synechococcus, but is completely absent in viral PsaA proteins from the psaJF→C→A→B→K→E→D gene organization and most Prochlorococcus strains. This may indicate a co-evolutionary scenario where cyanophages containing either of these gene organizations infect cyanobacterial ecotypes biogeographically restricted to the 30°N and 30°S equatorial strip.