Genetic organization of the psbAD region in phages infecting marine Synechococcus strains

Novel phages of Pseudomonas syringae unveil numerous potential auxiliary metabolic genes

Relatively few phages that infect plant pathogens have been isolated and investigated. The Pseudomonas syringae species complex is present in various environments, including plants. It can cause major crop diseases, such as bacterial canker on apricot trees. This study presents a collection of 25 unique phage genomes that infect P. syringae. These phages were isolated from apricot orchards with bacterial canker symptoms after enrichment with 21 strains of P. syringae. This collection comprises mostly virulent phages, with only three being temperate. They belong to 14 genera, 11 of which are newly discovered, and 18 new species, revealing great genetic diversity within this collection. Novel DNA packaging systems have been identified bioinformatically in one of the new phage species, but experimental confirmation is required to define the precise mechanism. Additionally, many phage genomes contain numerous potential auxiliary metabolic genes with diversified putative functions. At least three phages encode genes involved in bacterial tellurite resistance, a toxic metalloid. This suggests that viruses could play a role in bacterial stress tolerance. This research emphasizes the significance of continuing the search for new phages in the agricultural ecosystem to unravel novel ecological diversity and new gene functions. This work contributes to the foundation for future fundamental and applied research on phages infecting phytopathogenic bacteria.

High niche specificity and host genetic diversity of groundwater viruses

Viruses are key members of microbial communities that exert control over host abundance and metabolism, thereby influencing ecosystem processes and biogeochemical cycles. Aquifers are known to host taxonomically diverse microbial life, yet little is known about viruses infecting groundwater microbial communities. Here, we analysed 16 metagenomes from a broad range of groundwater physicochemistries. We recovered 1571 viral genomes that clustered into 468 high-quality viral operational taxonomic units. At least 15% were observed to be transcriptionally active, although lysis was likely constrained by the resource-limited groundwater environment. Most were unclassified (95%), and the remaining 5% were Caudoviricetes. Comparisons with viruses inhabiting other aquifers revealed no shared species, indicating substantial unexplored viral diversity. In silico predictions linked 22.4% of the viruses to microbial host populations, including to ultra-small prokaryotes, such as Patescibacteria and Nanoarchaeota. Many predicted hosts were associated with the biogeochemical cycling of carbon, nitrogen, and sulfur. Metabolic predictions revealed the presence of 205 putative auxiliary metabolic genes, involved in diverse processes associated with the utilization of the host's intracellular resources for biosynthesis and transformation reactions, including those involved in nucleotide sugar, glycan, cofactor, and vitamin metabolism. Viruses, prokaryotes overall, and predicted prokaryotic hosts exhibited narrow spatial distributions, and relative abundance correlations with the same groundwater parameters (e.g. dissolved oxygen, nitrate, and iron), consistent with host control over viral distributions. Results provide insights into underexplored groundwater viruses, and indicate the large extent to which viruses may manipulate microbial communities and biogeochemistry in the terrestrial subsurface.

Diversity of Bathyarchaeia viruses in metagenomes and virus-encoded CRISPR system components

Bathyarchaeia represent a class of archaea common and abundant in sedimentary ecosystems. The virome of Bathyarchaeia so far has not been characterized. Here we report 56 metagenome-assembled genomes of Bathyarchaeia viruses identified in metagenomes from different environments. Gene sharing network and phylogenomic analyses led to the proposal of four virus families, including viruses of the realms Duplodnaviria and Adnaviria, and archaea-specific spindle-shaped viruses. Genomic analyses uncovered diverse CRISPR elements in these viruses. Viruses of the proposed family 'Fuxiviridae' harbor an atypical type IV-B CRISPR-Cas system and a Cas4 protein that might interfere with host immunity. Viruses of the family 'Chiyouviridae' encode a Cas2-like endonuclease and two mini-CRISPR arrays, one with a repeat identical to that in the host CRISPR array, potentially allowing the virus to recruit the host CRISPR adaptation machinery to acquire spacers that could contribute to competition with other mobile genetic elements or to inhibition of host defenses. These findings present an outline of the Bathyarchaeia virome and offer a glimpse into their counter-defense mechanisms.

Characterization and genomic analysis of a novel Synechococcus phage S-H9-2 belonging to Bristolvirus genus isolated from the Yellow Sea

Cyanophages are known to influence the population dynamics and community structure of cyanobacteria and thus play an important role in biogeochemical cycles in aquatic ecosystems. In this study, a novel Synechococcus phage S-H9-2 infecting Synechococcus sp. WH 8102 was isolated from the coastal water of the Yellow Sea. Synechococcus phage S-H9-2 contains a 187,320 bp genome of double-stranded DNA with a G + C content of 40.3%, 202 potential open reading frames (ORFs), and 15 tRNAs. Phylogenetic analysis and nucleotide-based intergenomic similarity suggest that Synechococcus phage S-H9-2 belongs to the Bristolvirus genus under the family Kyanoviridae. Homologs of the S-H9-2 open reading frame can be found in a variety of marine environments, as shown by the results of mapping the genome sequence of S-H9-2 to the Global Ocean Viromes 2.0 dataset. The presence of auxiliary metabolic genes (AMGs) related to photosynthesis, carbon metabolism, and phosphorus assimilation, as well as phylogenetic relationships based on complete genome sequences, reflect the mechanism of phage-host interaction and host-specific strategies for adaptation to environmental conditions. This study enriches the current genomic database of cyanophage and contributed to our understanding of the virus-host interactions and their adaption to the environment.

Cyanophages as an important factor in the early evolution of oxygenic photosynthesis

Cyanophages are viruses that infect cyanobacteria. An interesting feature of many of them is the presence of psbA and psbD, genes that encode D1 and D2 proteins, respectively. The D1 and D2 are core proteins of the photosystem II (PSII) in cyanobacteria, algae and plants and influence the proper function of oxygenic photosynthesis (OP) in all oxyphototrophs on Earth. The frequent occurrence of psbA and psbD in cyanophages raises the question whether these genes coevolved with hosts during the early stages of cyanophage and cyanobacterial evolution, or whether they are direct descendants of genes adopted from the genomes of cyanobacterial hosts. The phylogeny of D1/D2 proteins encoded in the genomes of selected cyanophages and oxyphototrophs was reconstructed. In addition, common ancestral sequences of D1 and D2 proteins were predicted for cyanophages and oxyphototrophs. Based on this, the reconstruction of the 3D structures of D1 and D2 proteins was performed. In addition, the ratio of non-synonymous to synonymous (d_N/d_S) nucleotide substitutions in the coding sequences (CDSs) of psbA and psbD was determined. The results of the predicted spatial structures of the D1 and D2 proteins and purifying selection for the CDSs of psbA and psbD suggest that they belong to the ancient proteins, which may have formed the primordial PSII. It cannot be ruled out that they involved in water oxidation in cyanobacteria-like organisms at early stages of the evolution of life on Earth and coevolved with ancient cyanophages. The data are also discussed in the context of the origin of viruses.

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.

Priorities for ocean microbiome research

Microbial communities have essential roles in ocean ecology and planetary health. Microbes participate in nutrient cycles, remove huge quantities of carbon dioxide from the air and support ocean food webs. The taxonomic and functional diversity of the global ocean microbiome has been revealed by technological advances in sampling, DNA sequencing and bioinformatics. A better understanding of the ocean microbiome could underpin strategies to address environmental and societal challenges, including achievement of multiple Sustainable Development Goals way beyond SDG 14 'life below water'. We propose a set of priorities for understanding and protecting the ocean microbiome, which include delineating interactions between microbiota, sustainably applying resources from oceanic microorganisms and creating policy- and funder-friendly ocean education resources, and discuss how to achieve these ambitious goals.

Phage Genome Annotation: Where to Begin and End

With the renewed interest in phage research, coupled with the rising accessibility to affordable sequencing, ever increasing numbers of phage genomes are being sequenced. Therefore, there is an increased need to assemble and annotate phage genomes. There is a plethora of tools and platforms that allow phage genomes to be assembled and annotated. The choice of tools can often be bewildering for those new to phage genome assembly. Here we provide an overview of the assembly and annotation process from obtaining raw reads to genome submission, with worked examples, providing those new to genome assembly and annotation with a guided pathway to genome submission. We focus on the use of open access tools that can be incorporated into workflows to allow easy repetition of steps, highlighting multiple tools that can be used and common pitfalls that may occur.

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.

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.

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.

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.

Dynamic marine viral infections and major contribution to photosynthetic processes shown by spatiotemporal picoplankton metatranscriptomes

Viruses provide top-down control on microbial communities, yet their direct study in natural environments was hindered by culture limitations. The advance of bioinformatics enables cultivation-independent study of viruses. Many studies assemble new viral genomes and study viral diversity using marker genes from free viruses. Here we use cellular metatranscriptomics to study active community-wide viral infections. Recruitment to viral contigs allows tracking infection dynamics over time and space. Our assemblies represent viral populations, but appear biased towards low diversity viral taxa. Tracking relatives of published T4-like cyanophages and pelagiphages reveals high genomic continuity. We determine potential hosts by matching dynamics of infection with abundance of particular microbial taxa. Finally, we quantify the relative contribution of cyanobacteria and viruses to photosystem-II psbA (reaction center) expression in our study sites. We show sometimes >50% of all cyanobacterial+viral psbA expression is of viral origin, highlighting the contribution of viruses to photosynthesis and oxygen production.

Here, Sieradzki et al. use metatranscriptomics to study active community-wide viral infections at three coastal California sites throughout a year, identify potential viral hosts, and show that viruses can contribute a substantial amount to photosystem-II psbA expression.

F-CphI represents a new homing endonuclease family using the Endo VII catalytic motif

Background

There are six known families of homing endonucleases, LAGLIDADG, GIY-YIG, HNH, His-Cys box, PD-(D/E)-XK, and EDxHD, which are characterized by their conserved residues. Previously, we discovered a novel homing endonuclease F-CphI encoded by ORF177 of cyanophage S-PM2. F-CphI does not resemble any characterized homing endonucleases. Instead, the C-terminus of F-CphI aligns well with the N-terminal catalytic domain of a Holliday junction DNA resolvase, phage T4 endonuclease VII (Endo VII).

Results

A PSI-BLAST search resulted in a total of 313 Endo VII motif–containing sequences in sequenced genomes. Multiple sequence alignment showed that the catalytically important residues of T4 Endo VII were all well conserved in these proteins. Our site-directed mutagenesis studies further confirmed that the catalytically important residues of T4 Endo VII were also essential for F-CphI activity, and thus F-CphI might use a similar protein fold as Endo VII for DNA cleavage. A phylogenetic tree of the Endo VII motif–containing sequences showed that putative resolvases grouped into one clade while putative homing endonucleases and restriction endonucleases grouped into another clade.

Conclusions

Based on the unique conserved residues, we proposed that F-CphI represents a new homing endonuclease family, which was named the DHHRN family. Our phylogenetic analysis could be used to predict the functions of many previously unknown proteins.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0132-5) contains supplementary material, which is available to authorized users.

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).

Phage hunters: Computational strategies for finding phages in large-scale 'omics datasets

A plethora of tools exist for identifying phage sequences in bacterial genomes, single cell amplified genomes, and host-associated and environmental metagenomes. Yet because the genetics of phages and their hosts are closely intertwined, distinguishing viral from bacterial signal remains an ongoing challenge. Further the size, quantity and fragmentary nature of modern 'omics datasets ushers in a new set of computational challenges. Here, we detail the promises and pitfalls of using currently available gene-centric or k-mer based tools for identifying prophage sequences in genomes and prophage and viral contigs in metagenomes. Each of these methods offers a unique piece of the puzzle to elucidating the intriguing signatures of phage-host coevolution.

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.

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.

VirFinder: a novel k-mer based tool for identifying viral sequences from assembled metagenomic data

BACKGROUND

Identifying viral sequences in mixed metagenomes containing both viral and host contigs is a critical first step in analyzing the viral component of samples. Current tools for distinguishing prokaryotic virus and host contigs primarily use gene-based similarity approaches. Such approaches can significantly limit results especially for short contigs that have few predicted proteins or lack proteins with similarity to previously known viruses.

METHODS

We have developed VirFinder, the first k-mer frequency based, machine learning method for virus contig identification that entirely avoids gene-based similarity searches. VirFinder instead identifies viral sequences based on our empirical observation that viruses and hosts have discernibly different k-mer signatures. VirFinder's performance in correctly identifying viral sequences was tested by training its machine learning model on sequences from host and viral genomes sequenced before 1 January 2014 and evaluating on sequences obtained after 1 January 2014.

RESULTS

VirFinder had significantly better rates of identifying true viral contigs (true positive rates (TPRs)) than VirSorter, the current state-of-the-art gene-based virus classification tool, when evaluated with either contigs subsampled from complete genomes or assembled from a simulated human gut metagenome. For example, for contigs subsampled from complete genomes, VirFinder had 78-, 2.4-, and 1.8-fold higher TPRs than VirSorter for 1, 3, and 5 kb contigs, respectively, at the same false positive rates as VirSorter (0, 0.003, and 0.006, respectively), thus VirFinder works considerably better for small contigs than VirSorter. VirFinder furthermore identified several recently sequenced virus genomes (after 1 January 2014) that VirSorter did not and that have no nucleotide similarity to previously sequenced viruses, demonstrating VirFinder's potential advantage in identifying novel viral sequences. Application of VirFinder to a set of human gut metagenomes from healthy and liver cirrhosis patients reveals higher viral diversity in healthy individuals than cirrhosis patients. We also identified contig bins containing crAssphage-like contigs with higher abundance in healthy patients and a putative Veillonella genus prophage associated with cirrhosis patients.

CONCLUSIONS

This innovative k-mer based tool complements gene-based approaches and will significantly improve prokaryotic viral sequence identification, especially for metagenomic-based studies of viral ecology.

VirFinder: a novel k-mer based tool for identifying viral sequences from assembled metagenomic data

BACKGROUND

Identifying viral sequences in mixed metagenomes containing both viral and host contigs is a critical first step in analyzing the viral component of samples. Current tools for distinguishing prokaryotic virus and host contigs primarily use gene-based similarity approaches. Such approaches can significantly limit results especially for short contigs that have few predicted proteins or lack proteins with similarity to previously known viruses.

METHODS

We have developed VirFinder, the first k-mer frequency based, machine learning method for virus contig identification that entirely avoids gene-based similarity searches. VirFinder instead identifies viral sequences based on our empirical observation that viruses and hosts have discernibly different k-mer signatures. VirFinder's performance in correctly identifying viral sequences was tested by training its machine learning model on sequences from host and viral genomes sequenced before 1 January 2014 and evaluating on sequences obtained after 1 January 2014.

RESULTS

VirFinder had significantly better rates of identifying true viral contigs (true positive rates (TPRs)) than VirSorter, the current state-of-the-art gene-based virus classification tool, when evaluated with either contigs subsampled from complete genomes or assembled from a simulated human gut metagenome. For example, for contigs subsampled from complete genomes, VirFinder had 78-, 2.4-, and 1.8-fold higher TPRs than VirSorter for 1, 3, and 5 kb contigs, respectively, at the same false positive rates as VirSorter (0, 0.003, and 0.006, respectively), thus VirFinder works considerably better for small contigs than VirSorter. VirFinder furthermore identified several recently sequenced virus genomes (after 1 January 2014) that VirSorter did not and that have no nucleotide similarity to previously sequenced viruses, demonstrating VirFinder's potential advantage in identifying novel viral sequences. Application of VirFinder to a set of human gut metagenomes from healthy and liver cirrhosis patients reveals higher viral diversity in healthy individuals than cirrhosis patients. We also identified contig bins containing crAssphage-like contigs with higher abundance in healthy patients and a putative Veillonella genus prophage associated with cirrhosis patients.

CONCLUSIONS

This innovative k-mer based tool complements gene-based approaches and will significantly improve prokaryotic viral sequence identification, especially for metagenomic-based studies of viral ecology.

Putative archaeal viruses from the mesopelagic ocean

Oceanic viruses that infect bacteria, or phages, are known to modulate host diversity, metabolisms, and biogeochemical cycling, while the viruses that infect marine Archaea remain understudied despite the critical ecosystem roles played by their hosts. Here we introduce "MArVD", for Metagenomic Archaeal Virus Detector, an annotation tool designed to identify putative archaeal virus contigs in metagenomic datasets. MArVD is made publicly available through the online iVirus analytical platform. Benchmarking analysis of MArVD showed it to be >99% accurate and 100% sensitive in identifying the 127 known archaeal viruses among the 12,499 viruses in the VirSorter curated dataset. Application of MArVD to 10 viral metagenomes from two depth profiles in the Eastern Tropical North Pacific (ETNP) oxygen minimum zone revealed 43 new putative archaeal virus genomes and large genome fragments ranging in size from 10 to 31 kb. Network-based classifications, which were consistent with marker gene phylogenies where available, suggested that these putative archaeal virus contigs represented six novel candidate genera. Ecological analyses, via fragment recruitment and ordination, revealed that the diversity and relative abundances of these putative archaeal viruses were correlated with oxygen concentration and temperature along two OMZ-spanning depth profiles, presumably due to structuring of the host Archaea community. Peak viral diversity and abundances were found in surface waters, where Thermoplasmata 16S rRNA genes are prevalent, suggesting these archaea as hosts in the surface habitats. Together these findings provide a baseline for identifying archaeal viruses in sequence datasets, and an initial picture of the ecology of such viruses in non-extreme environments.

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.

In Vitro Structural and Functional Characterization of the Small Heat Shock Proteins (sHSP) of the Cyanophage S-ShM2 and Its Host, Synechococcus sp. WH7803

We previously reported the in silico characterization of Synechococcus sp. phage 18 kDa small heat shock protein (HspSP-ShM2). This small heat shock protein (sHSP) contains a highly conserved core alpha crystalline domain of 92 amino acids and relatively short N- and C-terminal arms, the later containing the classical C-terminal anchoring module motif (L-X-I/L/V). Here we establish the oligomeric profile of HspSP-ShM2 and its structural dynamics under in vitro experimental conditions using size exclusion chromatography (SEC/FPLC), gradient native gels electrophoresis and dynamic light scattering (DLS). Under native conditions, HspSP-ShM2 displays the ability to form large oligomers and shows a polydisperse profile. At higher temperatures, it shows extensive structural dynamics and undergoes conformational changes through an increased of subunit rearrangement and formation of sub-oligomeric species. We also demonstrate its capacity to prevent the aggregation of citrate synthase, malate dehydrogenase and luciferase under heat shock conditions through the formation of stable and soluble hetero-oligomeric complexes (sHSP:substrate). In contrast, the host cyanobacteria Synechococcus sp. WH7803 15 kDa sHSP (HspS-WH7803) aggregates when in the same conditions as HspSP-ShM2. However, its solubility can be maintained in the presence of non-ionic detergent Triton^™X-100 and forms an oligomeric structure estimated to be between dimer and tetramer but exhibits no apparent inducible structural dynamics neither chaperon-like activity in all the assays and molar ratios tested. SEC/FPLC and thermal aggregation prevention assays results indicate no formation of hetero-oligomeric complex or functional interactions between both sHSPs. Taken together these in vitro results portray the phage HspSP-ShM2 as a classical sHSP and suggest that it may be functional at the in vivo level while behaving differently than its host amphitropic sHSP.

Marine and giant viruses as indicators of a marine microbial community in a riverine system

Viral communities are important for ecosystem function as they are involved in critical biogeochemical cycles and controlling host abundance. This study investigates riverine viral communities around a small rural town that influences local water inputs. Myoviridae, Siphoviridae, Phycodnaviridae, Mimiviridae, Herpesviridae, and Podoviridae were the most abundant families. Viral species upstream and downstream of the town were similar, with Synechoccocus phage, salinus, Prochlorococcus phage, Mimivirus A, and Human herpes 6A virus most abundant, contributing to 4.9-38.2% of average abundance within the metagenomic profiles, with Synechococcus and Prochlorococcus present in metagenomes as the expected hosts for the phage. Overall, the majority of abundant viral species were or were most similar to those of marine origin. At over 60 km to the river mouth, the presence of marine communities provides some support for the Baas-Becking hypothesis "everything is everywhere, but, the environment selects." We conclude marine microbial species may occur more frequently in freshwater systems than previously assumed, and hence may play important roles in some freshwater ecosystems within tens to a hundred kilometers from the sea.

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.

Genome analysis and gene nblA identification of Microcystis aeruginosa myovirus (MaMV-DC) reveal the evidence for horizontal gene transfer events between cyanomyovirus and host

The genome sequence, genetic characterization and nblA gene function of Microcystis aeruginosa myovirus isolated from Lake Dianchi in China (MaMV-DC) have been analysed. The genome DNA is 169 223 bp long, with 170 predicted protein-coding genes (001L–170L) and a tRNA gene. About one-sixth of these genes have homologues in the host cyanobacteria M. aeruginosa. The genome carries a gene homologous to host nblA, which encodes a protein involved in the degradation of cyanobacterial phycobilisome. Its expression during MaMV-DC infection was confirmed by reverse transcriptase PCR and Western blot detection and abundant expression was companied by the significant decline of phycocyanin content and massive release of progeny MaMV-DC. In addition, expressing MaMV-DC nblA reduced the phycocyanin peak and the phycocyanin to chlorophyll ratio in model cyanobacteria. These results confirm that horizontal gene transfer events have occurred between cyanobacterial host and cyanomyovirus and suggest that MaMV-DC carrying host-derived genes (such as 005L, that codes for NblA) is responsible for more efficient expression of cyanophage genes and release of progeny cyanophage. This study provides novel insight into the horizontal gene transfer in cyanophage and the interactions between cyanophage and their host.

Assessing Illumina technology for the high-throughput sequencing of bacteriophage genomes

Bacteriophages are the most abundant biological entities on the planet, playing crucial roles in the shaping of bacterial populations. Phages have smaller genomes than their bacterial hosts, yet there are currently fewer fully sequenced phage than bacterial genomes. We assessed the suitability of Illumina technology for high-throughput sequencing and subsequent assembly of phage genomes. In silico datasets reveal that 30× coverage is sufficient to correctly assemble the complete genome of ~98.5% of known phages, with experimental data confirming that the majority of phage genomes can be assembled at 30× coverage. Furthermore, in silico data demonstrate it is possible to co-sequence multiple phages from different hosts, without introducing assembly errors.

PhagePhisher: a pipeline for the discovery of covert viral sequences in complex genomic datasets

Obtaining meaningful viral information from large sequencing datasets presents unique challenges distinct from prokaryotic and eukaryotic sequencing efforts. The difficulties surrounding this issue can be ascribed in part to the genomic plasticity of viruses themselves as well as the scarcity of existing information in genomic databases. The open-source software PhagePhisher (http://www.putonti-lab.com/phagephisher) has been designed as a simple pipeline to extract relevant information from complex and mixed datasets, and will improve the examination of bacteriophages, viruses, and virally related sequences, in a range of environments. Key aspects of the software include speed and ease of use; PhagePhisher can be used with limited operator knowledge of bioinformatics on a standard workstation. As a proof-of-concept, PhagePhisher was successfully implemented with bacteria–virus mixed samples of varying complexity. Furthermore, viral signals within microbial metagenomic datasets were easily and quickly identified by PhagePhisher, including those from prophages as well as lysogenic phages, an important and often neglected aspect of examining phage populations in the environment. PhagePhisher resolves viral-related sequences which may be obscured by or imbedded in bacterial genomes.

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.

Shedding new light on viral photosynthesis

Viruses infecting the environmentally important marine cyanobacteria Prochlorococcus and Synechococcus encode 'auxiliary metabolic genes' (AMGs) involved in the light and dark reactions of photosynthesis. Here, we discuss progress on the inventory of such AMGs in the ever-increasing number of viral genome sequences as well as in metagenomic datasets. We contextualise these gene acquisitions with reference to a hypothesised fitness gain to the phage. We also report new evidence with regard to the sequence and predicted structural properties of viral petE genes encoding the soluble electron carrier plastocyanin. Viral copies of PetE exhibit extensive modifications to the N-terminal signal peptide and possess several novel residues in a region responsible for interaction with redox partners. We also highlight potential knowledge gaps in this field and discuss future opportunities to discover novel phage-host interactions involved in the photosynthetic process.

Using signature genes as tools to assess environmental viral ecology and diversity

Viruses (including bacteriophages) are the most abundant biological entities on the planet. As such, they are thought to have a major impact on all aspects of microbial community structure and function. Despite this critical role in ecosystem processes, the study of virus/phage diversity has lagged far behind parallel studies of the bacterial and eukaryotic kingdoms, largely due to the absence of any universal phylogenetic marker. Here we review the development and use of signature genes to investigate viral diversity, as a viable strategy for data sets of specific virus groups. Genes that have been used include those encoding structural proteins, such as portal protein, major capsid protein, and tail sheath protein, auxiliary metabolism genes, such as psbA, psbB,and phoH, and several polymerase genes. These marker genes have been used in combination with PCR-based fingerprinting and/or sequencing strategies to investigate spatial, temporal, and seasonal variations and diversity in a wide range of habitats.

Phage Phenomics: Physiological Approaches to Characterize Novel Viral Proteins

Current investigations into phage-host interactions are dependent on extrapolating knowledge from (meta)genomes. Interestingly, 60 - 95% of all phage sequences share no homology to current annotated proteins. As a result, a large proportion of phage genes are annotated as hypothetical. This reality heavily affects the annotation of both structural and auxiliary metabolic genes. Here we present phenomic methods designed to capture the physiological response(s) of a selected host during expression of one of these unknown phage genes. Multi-phenotype Assay Plates (MAPs) are used to monitor the diversity of host substrate utilization and subsequent biomass formation, while metabolomics provides bi-product analysis by monitoring metabolite abundance and diversity. Both tools are used simultaneously to provide a phenotypic profile associated with expression of a single putative phage open reading frame (ORF). Representative results for both methods are compared, highlighting the phenotypic profile differences of a host carrying either putative structural or metabolic phage genes. In addition, the visualization techniques and high throughput computational pipelines that facilitated experimental analysis are presented.

Polar freshwater cyanophage S-EIV1 represents a new widespread evolutionary lineage of phages

Cyanobacteria are often the dominant phototrophs in polar freshwater communities; yet, the phages that infect them remain unknown. Here, we present a genomic and morphological characterization of cyanophage S-EIV1 that was isolated from freshwaters on Ellesmere Island (Nunavut, High Arctic Canada), and which infects the polar Synechococcus sp., strain PCCC-A2c. S-EIV1 represents a newly discovered evolutionary lineage of bacteriophages whose representatives are widespread in aquatic systems. Among the 130 predicted open reading frames (ORFs) there is no recognizable similarity to genes that encode structural proteins other than the large terminase subunit and a distant viral morphogenesis protein, indicating that the genes encoding the structural proteins of S-EIV1 are distinct from other viruses. As well, only 19 predicted coding sequences on the 79 178 bp circularly permuted genome have homology with genes encoding proteins of known function. Although S-EIV1 is divergent from other sequenced phage isolates, it shares synteny with phage genes captured on a fosmid from the deep-chlorophyll maximum in the Mediterranean Sea, as well as with an incision element in the genome of Anabaena variabilis (ATCC 29413). Sequence recruitment of metagenomic data indicates that S-EIV1-like viruses are cosmopolitan and abundant in a wide range of aquatic systems, suggesting they have an important ecological role.

Genomic characteristics and environmental distributions of the uncultivated Far-T4 phages

Viral metagenomics (viromics) is a tremendous tool to reveal viral taxonomic and functional diversity across ecosystems ranging from the human gut to the world's oceans. As with microbes however, there appear vast swaths of “dark matter” yet to be documented for viruses, even among relatively well-studied viral types. Here, we use viromics to explore the “Far-T4 phages” sequence space, a neighbor clade from the well-studied T4-like phages that was first detected through PCR study in seawater and subsequently identified in freshwater lakes through 454-sequenced viromes. To advance the description of these viruses beyond this single marker gene, we explore Far-T4 genome fragments assembled from two deeply-sequenced freshwater viromes. Single gene phylogenetic trees confirm that the Far-T4 phages are divergent from the T4-like phages, genome fragments reveal largely collinear genome organizations, and both data led to the delineation of five Far-T4 clades. Three-dimensional models of major capsid proteins are consistent with a T4-like structure, and highlight a highly conserved core flanked by variable insertions. Finally, we contextualize these now better characterized Far-T4 phages by re-analyzing 196 previously published viromes. These suggest that Far-T4 are common in freshwater and seawater as only four of 82 aquatic viromes lacked Far-T4-like sequences. Variability in representation across the five newly identified clades suggests clade-specific niche differentiation may be occurring across the different biomes, though the underlying mechanism remains unidentified. While complete genome assembly from complex communities and the lack of host linkage information still bottleneck virus discovery through viromes, these findings exemplify the power of metagenomics approaches to assess the diversity, evolutionary history, and genomic characteristics of novel uncultivated phages.

An inexpensive, accurate, and precise wet-mount method for enumerating aquatic viruses

Viruses affect biogeochemical cycling, microbial mortality, gene flow, and metabolic functions in diverse environments through infection and lysis of microorganisms. Fundamental to quantitatively investigating these roles is the determination of viral abundance in both field and laboratory samples. One current, widely used method to accomplish this with aquatic samples is the "filter mount" method, in which samples are filtered onto costly 0.02-μm-pore-size ceramic filters for enumeration of viruses by epifluorescence microscopy. Here we describe a cost-effective (ca. 500-fold-lower materials cost) alternative virus enumeration method in which fluorescently stained samples are wet mounted directly onto slides, after optional chemical flocculation of viruses in samples with viral concentrations of <5×10(7) viruses ml(-1). The concentration of viruses in the sample is then determined from the ratio of viruses to a known concentration of added microsphere beads via epifluorescence microscopy. Virus concentrations obtained by using this wet-mount method, with and without chemical flocculation, were significantly correlated with, and had precision equivalent to, those obtained by the filter mount method across concentrations ranging from 2.17×10(6) to 1.37×10(8) viruses ml(-1) when tested by using cultivated viral isolates and natural samples from marine and freshwater environments. In summary, the wet-mount method is significantly less expensive than the filter mount method and is appropriate for rapid, precise, and accurate enumeration of aquatic viruses over a wide range of viral concentrations (≥1×10(6) viruses ml(-1)) encountered in field and laboratory samples.

Genomic and proteomic characterization of two novel siphovirus infecting the sedentary facultative epibiont cyanobacterium Acaryochloris marina

Acaryochloris marina is a symbiotic species of cyanobacteria that is capable of utilizing far-red light. We report the characterization of the phages A-HIS1 and A-HIS2, capable of infecting Acaryochloris. Morphological characterization of these phages places them in the family Siphoviridae. However, molecular characterization reveals that they do not show genetic similarity with any known siphoviruses. While the phages do show synteny between each other, the nucleotide identity between the phages is low at 45-67%, suggesting they diverged from each other some time ago. The greatest number of genes shared with another phage (a myovirus infecting marine Synechococcus) was four. Unlike most other cyanophages and in common with the Siphoviridae infecting Synechococcus, no photosynthesis-related genes were found in the genome. CRISPR (clustered regularly interspaced short palindromic repeats) spacers from the host Acaryochloris had partial matches to sequences found within the phages, which is the first time CRISPRs have been reported in a cyanobacterial/cyanophage system. The phages also encode a homologue of the proteobacterial RNase T. The potential function of RNase T in the mark-up or digestion of crRNA hints at a novel mechanism for evading the host CRISPR system.

Diversity of viral photosystem-I psaA genes

Marine photosynthesis is one of the major contributors to the global carbon cycle and the world's oxygen supply. This process is largely driven by cyanobacteria, namely Synechococcus and Prochlorococcus. Genes encoding photosystem-II (PSII) reaction center proteins are found in many cyanophage genomes, and are expressed during the infection of their hosts. On the basis of metagenomics, cyanophage photosystem-I (PSI) gene cassettes were recently discovered with two gene arrangements psaJF→C→A→B→K→E→D and psaD→C→A→B. It was suggested that the horizontal transfer of PSII and PSI genes is increasing phage fitness. To better understand their diversity, we designed degenerate primers to cover a wide diversity of organisms, and using PCR we targeted the psaC→A arrangement, which is unique to cyanophages cassettes. We examined viral concentrates from four islands in the Pacific Ocean and found samples containing the psaC→A arrangement. Analyses of the amplified viral psaA gene revealed six subgroups varying in their level of similarity and %G+C content, suggesting that the diversity of cyanophage PSI genes is greater than originally thought.

Phages in nature

Bacteriophages or phages are the most abundant organisms in the biosphere and they are a ubiquitous feature of prokaryotic existence. A bacteriophage is a virus which infects a bacterium. Archaea are also infected by viruses, whether these should be referred to as 'phages' is debatable, but they are included as such in the scope this article. Phages have been of interest to scientists as tools to understand fundamental molecular biology, as vectors of horizontal gene transfer and drivers of bacterial evolution, as sources of diagnostic and genetic tools and as novel therapeutic agents. Unraveling the biology of phages and their relationship with their hosts is key to understanding microbial systems and their exploitation. In this article we describe the roles of phages in different host systems and show how modeling, microscopy, isolation, genomic and metagenomic based approaches have come together to provide unparalleled insights into these small but vital constituents of the microbial world.

Supersize me: Cronobacter sakazakii phage GAP32

Cronobacter sakazakii is a Gram-negative pathogen found in milk-based formulae that causes infant meningitis. Bacteriophages have been proposed to control bacterial pathogens; however, comprehensive knowledge about a phage is required to ensure its safety before clinical application. We have characterized C. sakazakii phage vB_CsaM_GAP32 (GAP32), which possesses the second largest sequenced phage genome (358,663bp). A total of 571 genes including 545 protein coding sequences and 26 tRNAs were identified, thus more genes than in the smallest bacterium, Mycoplasma genitalium G37. BLASTP and HHpred searches, together with proteomic analyses reveal that only 23.9% of the putative proteins have defined functions. Some of the unique features of this phage include: a chromosome condensation protein, two copies of the large subunit terminase, a predicted signal-arrest-release lysin; and an RpoD-like protein, which is possibly involved in the switch from immediate early to delayed early transcription. Its closest relatives are all extremely large myoviruses, namely coliphage PBECO4 and Klebsiella phage vB_KleM-RaK2, with whom it shares approximately 44% homologous proteins. Since the homologs are not evenly distributed, we propose that these three phages belong to a new subfamily.

Comparative metagenomic analyses reveal viral-induced shifts of host metabolism towards nucleotide biosynthesis

BACKGROUND

Viral genomes often contain metabolic genes that were acquired from host genomes (auxiliary genes). It is assumed that these genes are fixed in viral genomes as a result of a selective force, favoring viruses that acquire specific metabolic functions. While many individual auxiliary genes were observed in viral genomes and metagenomes, there is great importance in investigating the abundance of auxiliary genes and metabolic functions in the marine environment towards a better understanding of their role in promoting viral reproduction.

RESULTS

In this study, we searched for enriched viral auxiliary genes and mapped them to metabolic pathways. To initially identify enriched auxiliary genes, we analyzed metagenomic microbial reads from the Global Ocean Survey (GOS) dataset that were characterized as viral, as well as marine virome and microbiome datasets from the Line Islands. Viral-enriched genes were mapped to a "global metabolism network" that comprises all KEGG metabolic pathways. Our analysis of the viral-enriched pathways revealed that purine and pyrimidine metabolism pathways are among the most enriched pathways. Moreover, many other viral-enriched metabolic pathways were found to be closely associated with the purine and pyrimidine metabolism pathways. Furthermore, we observed that sequential reactions are promoted in pathways having a high proportion of enriched genes. In addition, these enriched genes were found to be of modular nature, participating in several pathways.

CONCLUSIONS

Our naïve metagenomic analyses strongly support the well-established notion that viral auxiliary genes promote viral replication via both degradation of host DNA and RNA as well as a shift of the host metabolism towards nucleotide biosynthesis, clearly indicating that comparative metagenomics can be used to understand different environments and systems without prior knowledge of pathways involved.

The evolutionary divergence of psbA gene in Synechococcus and their myoviruses in the East China Sea

Marine Synechococcus is a principal component of the picophytoplankton and makes an important contribution to primary productivity in the ocean. Synechophages, infecting Synechococcus, are believed to have significant influences on the distribution and abundance of their hosts. Extensive previous ecological studies on cyanobacteria and viruses have been carried out in the East China Sea (ECS). Here we investigate the diversity and divergence of Synechococcus and their myoviruses (Synechomyoviruses) based on their shared photosynthesis psbA gene. Synechococcus is dominated by subclades 5.1A I, 5.1A II and 5.1A IV in the ECS, and clades I and II are the dominant groups in the Synechomyoviruses. As two phylogenetically independent clades, there is much higher diversity of the Synechomyoviruses than Synechococcus. Obvious partitioning characteristics of GC and GC3 (the GC content at the third codon position) contents are obtained among different picophytoplankton populations and their phages. The GC3 content causes the psbA gene in Synechococcus to have a higher GC content, while the opposite is true in the Synechomyoviruses. Analyzing more than one-time difference of the codon usage frequency of psbA sequences, the third position nucleotides of preferred codons for Synechococcus are all G and C, while most Synechomyoviral sequences (72.7%) have A and T at the third position of their preferred codons. This work shed light on the ecology and evolution of phage-host interactions in the environment.