Reconstruction of novel cyanobacterial siphovirus genomes from Mediterranean metagenomic fosmids

Increased cyanophage infection at the bottom of the euphotic zone, especially in the fall

Picocyanobacteria contribute greatly to offshore primary production with cells extending through the deep euphotic zone. Literature indicates high viral infection of cyanobacteria in ocean transition zones. We postulate that the bottom of the euphotic zone is a transition zone, where communities transition from phototrophic to aphotic processes. We use single-copy core genes to examine cyanophage to cyanobacteria ratios in cellular metagenomes in the subtropical North Atlantic and Pacific. Cyanophage to cyanobacteria terL/rpoB ratios generally increase to >10 in the deep euphotic zone. As light levels decrease in the fall, Prochlorococcus in the deep euphotic zone experience reduced light levels. We find clear differences between spring (Geotraces GA02) and fall (GA03) in the North Atlantic, with terL/rpoB ratios increasing to >40 in the fall. When examining 23 months of the North Pacific Hawaii Ocean Timeseries, the depth of elevated cyanophage to cyanobacteria ratios in cellular metagenomes negatively correlated with surface photosynthetic radiation (PAR), particularly with the change in PAR, which reflected the season. In fall, all picocyanobacteria ecotypes were found at depths enriched with viruses, while in summer, only low light ecotypes were affected. Thus, we find high cyanophage infection both in the deep euphotic zone and during seasonal transitions.

Associations between picocyanobacterial ecotypes and cyanophage host genes across ocean basins and depth

Background

Cyanophages, viruses that infect cyanobacteria, are globally abundant in the ocean's euphotic zone and are a potentially important cause of mortality for marine picocyanobacteria. Viral host genes are thought to increase viral fitness by either increasing numbers of genes for synthesizing nucleotides for virus replication, or by mitigating direct stresses imposed by the environment. The encoding of host genes in viral genomes through horizontal gene transfer is a form of evolution that links viruses, hosts, and the environment. We previously examined depth profiles of the proportion of cyanophage containing various host genes in the Eastern Tropical North Pacific Oxygen Deficient Zone (ODZ) and at the subtropical North Atlantic (BATS). However, cyanophage host genes have not been previously examined in environmental depth profiles across the oceans.

Methodology

We examined geographical and depth distributions of picocyanobacterial ecotypes, cyanophage, and their viral-host genes across ocean basins including the North Atlantic, Mediterranean Sea, North Pacific, South Pacific, and Eastern Tropical North and South Pacific ODZs using phylogenetic metagenomic read placement. We determined the proportion of myo and podo-cyanophage containing a range of host genes by comparing to cyanophage single copy core gene terminase (terL). With this large dataset (22 stations), network analysis identified statistical links between 12 of the 14 cyanophage host genes examined here with their picocyanobacteria host ecotypes.

Results

Picyanobacterial ecotypes, and the composition and proportion of cyanophage host genes, shifted dramatically and predictably with depth. For most of the cyanophage host genes examined here, we found that the composition of host ecotypes predicted the proportion of viral host genes harbored by the cyanophage community. Terminase is too conserved to illuminate the myo-cyanophage community structure. Cyanophage cobS was present in almost all myo-cyanophage and did not vary in proportion with depth. We used the composition of cobS phylotypes to track changes in myo-cyanophage composition.

Conclusions

Picocyanobacteria ecotypes shift with changes in light, temperature, and oxygen and many common cyanophage host genes shift concomitantly. However, cyanophage phosphate transporter gene pstS appeared to instead vary with ocean basin and was most abundant in low phosphate regions. Abundances of cyanophage host genes related to nutrient acquisition may diverge from host ecotype constraints as the same host can live in varying nutrient concentrations. Myo-cyanophage community in the anoxic ODZ had reduced diversity. By comparison to the oxic ocean, we can see which cyanophage host genes are especially abundant (nirA, nirC, and purS) or not abundant (myo psbA) in ODZs, highlighting both the stability of conditions in the ODZ and the importance of nitrite as an N source to ODZ endemic LLV Prochlorococcus.

Exploring Viral Diversity in a Gypsum Karst Lake Ecosystem Using Targeted Single-Cell Genomics

Little is known about the diversity and distribution of viruses infecting green sulfur bacteria (GSB) thriving in euxinic (sulfuric and anoxic) habitats, including gypsum karst lake ecosystems. In this study, we used targeted cell sorting combined with single-cell sequencing to gain insights into the gene content and genomic potential of viruses infecting sulfur-oxidizing bacteria Chlorobium clathratiforme, obtained from water samples collected during summer stratification in gypsum karst Lake Kirkilai (Lithuania). In total, 82 viral contigs were bioinformatically identified in 62 single amplified genomes (SAGs) of C. clathratiforme. The majority of viral gene and protein sequences showed little to no similarity with phage sequences in public databases, uncovering the vast diversity of previously undescribed GSB viruses. We observed a high level of lysogenization in the C. clathratiforme population, as 87% SAGs contained intact prophages. Among the thirty identified auxiliary metabolic genes (AMGs), two, thiosulfate sulfurtransferase (TST) and thioredoxin-dependent phosphoadenosine phosphosulfate (PAPS) reductase (cysH), were found to be involved in the oxidation of inorganic sulfur compounds, suggesting that viruses can influence the metabolism and cycling of this essential element. Finally, the analysis of CRISPR spacers retrieved from the consensus C. clathratiforme genome imply persistent and active virus–host interactions for several putative phages prevalent among C. clathratiforme SAGs. Overall, this study provides a glimpse into the diversity of phages associated with naturally occurring and highly abundant sulfur-oxidizing bacteria.

The Novel Halovirus Hardycor1, and the Presence of Active (Induced) Proviruses in Four Haloarchaea

The virus Hardycor1 was isolated in 1998 and infects the haloarchaeon Halorubrum coriense. DNA from a frozen stock (HC1) was sequenced and the viral genome found to be 45,142 bp of dsDNA, probably having redundant, circularly permuted termini. The genome showed little similarity (BLASTn) to known viruses. Only twenty-two of the 53 (41%) predicted proteins were significantly similar to sequences in the NCBI nr protein database (E-value ≤ 10^-15). Six caudovirus-like proteins were encoded, including large subunit terminase (TerL), major capsid protein (Mcp) and tape measure protein (Tmp). Hardycor1 was predicted to be a siphovirus (VIRFAM). No close relationship to other viruses was found using phylogenetic tree reconstructions based on TerL and Mcp. Unexpectedly, the sequenced virus stock HC1 also revealed two induced proviruses of the host: a siphovirus (Humcor1) and a pleolipovirus (Humcor2). A re-examination of other similarly sequenced, archival virus stocks revealed induced proviruses of Haloferax volcanii, Haloferax gibbonsii and Haloarcula hispanica, three of which were pleolipoviruses. One provirus (Halfvol2) of Hfx. volcanii showed little similarity (BLASTn) to known viruses and probably represents a novel virus group. The attP sequences of many pleolipoproviruses were found to be embedded in a newly detected coding sequence, split in the provirus state, that spans between genes for integrase and a downstream CxxC-motif protein. This gene might play an important role in regulation of the temperate state.

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.

Assembly-free single-molecule sequencing recovers complete virus genomes from natural microbial communities

Viruses are the most abundant biological entities on Earth and play key roles in host ecology, evolution, and horizontal gene transfer. Despite recent progress in viral metagenomics, the inherent genetic complexity of virus populations still poses technical difficulties for recovering complete virus genomes from natural assemblages. To address these challenges, we developed an assembly-free, single-molecule nanopore sequencing approach, enabling direct recovery of complete virus genome sequences from environmental samples. Our method yielded thousands of full-length, high-quality draft virus genome sequences that were not recovered using standard short-read assembly approaches. Additionally, our analyses discriminated between populations whose genomes had identical direct terminal repeats versus those with circularly permuted repeats at their termini, thus providing new insight into native virus reproduction and genome packaging. Novel DNA sequences were discovered, whose repeat structures, gene contents, and concatemer lengths suggest they are phage-inducible chromosomal islands, which are packaged as concatemers in phage particles, with lengths that match the size ranges of co-occurring phage genomes. Our new virus sequencing strategy can provide previously unavailable information about the genome structures, population biology, and ecology of naturally occurring viruses and viral parasites.

Genome-resolved viral and cellular metagenomes revealed potential key virus-host interactions in a deep freshwater lake

Metagenomics has dramatically expanded the known virosphere, but freshwater viral diversity and their ecological interaction with hosts remain poorly understood. Here, we conducted a metagenomic exploration of planktonic dsDNA prokaryotic viruses by sequencing both virion (<0.22 μm) and cellular (0.22-5.0 μm) fractions collected spatiotemporally from a deep freshwater lake (Lake Biwa, Japan). This simultaneously reconstructed 183 complete (i.e., circular) viral genomes and 57 bacterioplankton metagenome-assembled genomes. Analysis of metagenomic read coverage revealed vertical partitioning of the viral community analogous to the vertically stratified bacterioplankton community. The hypolimnetic community was generally stable during stratification, but occasionally shifted abruptly, presumably due to lysogenic induction. Genes involved in assimilatory sulfate reduction were encoded in 20 (10.9%) viral genomes, including those of dominant viruses, and may aid viral propagation in sulfur-limited freshwater systems. Hosts were predicted for 40 (21.9%) viral genomes, encompassing 10 phyla (or classes of Proteobacteria) including ubiquitous freshwater bacterioplankton lineages (e.g., Ca. Fonsibacter and Ca. Nitrosoarchaeum). Comparison with viral genomes derived from published metagenomes revealed viral phylogeographic connectivity in geographically isolated habitats. Notably, analogous to their hosts, actinobacterial viruses were among the most diverse, ubiquitous and abundant viral groups in freshwater systems, with potential high lytic activity in surface waters.

Phage-centric ecological interactions in aquatic ecosystems revealed through ultra-deep metagenomics

The persistent inertia in the ability to culture environmentally abundant microbes from aquatic ecosystems represents an obstacle in disentangling the complex web of ecological interactions spun by a diverse assortment of participants (pro- and eukaryotes and their viruses). In aquatic microbial communities, the numerically most abundant actors, the viruses, remain the most elusive, and especially in freshwaters their identities and ecology remain unknown. Here, using ultra-deep metagenomic sequencing from pelagic freshwater habitats, we recovered complete genomes of > 2000 phages, including small "miniphages" and large "megaphages" infecting iconic freshwater prokaryotic lineages. For instance, abundant freshwater Actinobacteria support infection by a very broad size range of phages (13-200 Kb). We describe many phages encoding genes that likely afford protection to their host from reactive oxygen species (ROS) in the aquatic environment and in the oxidative burst in protist phagolysosomes (phage-mediated ROS defense). Spatiotemporal abundance analyses of phage genomes revealed evanescence as the primary dynamic in upper water layers, where they displayed short-lived existences. In contrast, persistence was characteristic for the deeper layers where many identical phage genomes were recovered repeatedly. Phage and host abundances corresponded closely, with distinct populations displaying preferential distributions in different seasons and depths, closely mimicking overall stratification and mixis.

Phage-centric ecological interactions in aquatic ecosystems revealed through ultra-deep metagenomics

The persistent inertia in the ability to culture environmentally abundant microbes from aquatic ecosystems represents an obstacle in disentangling the complex web of ecological interactions spun by a diverse assortment of participants (pro- and eukaryotes and their viruses). In aquatic microbial communities, the numerically most abundant actors, the viruses, remain the most elusive, and especially in freshwaters their identities and ecology remain unknown. Here, using ultra-deep metagenomic sequencing from pelagic freshwater habitats, we recovered complete genomes of > 2000 phages, including small "miniphages" and large "megaphages" infecting iconic freshwater prokaryotic lineages. For instance, abundant freshwater Actinobacteria support infection by a very broad size range of phages (13-200 Kb). We describe many phages encoding genes that likely afford protection to their host from reactive oxygen species (ROS) in the aquatic environment and in the oxidative burst in protist phagolysosomes (phage-mediated ROS defense). Spatiotemporal abundance analyses of phage genomes revealed evanescence as the primary dynamic in upper water layers, where they displayed short-lived existences. In contrast, persistence was characteristic for the deeper layers where many identical phage genomes were recovered repeatedly. Phage and host abundances corresponded closely, with distinct populations displaying preferential distributions in different seasons and depths, closely mimicking overall stratification and mixis.

Metagenomic Characterization of the Viral Community of the South Scotia Ridge

Viruses are the most abundant biological entities in aquatic ecosystems and harbor an enormous amount of genetic diversity. Whereas their influence on marine ecosystems is widely acknowledged, current information about their diversity remains limited. We conducted a viral metagenomic analysis of water samples collected during the austral summer of 2016 from the South Scotia Ridge (SSR), near the Antarctic Peninsula. The taxonomic composition and diversity of the viral communities were investigated, and a functional assessment of the sequences was performed. Phylotypic analysis showed that most viruses belonged to the order Caudovirales, especially the family Podoviridae (41.92⁻48.7%), which is similar to the situation in the Pacific Ocean. Functional analysis revealed a relatively high frequency of phage-associated and metabolism genes. Phylogenetic analyses of phage TerL and Capsid_NCLDV (nucleocytoplasmic large DNA viruses) marker genes indicated that many sequences associated with Caudovirales and NCLDV were novel and distinct from known phage genomes. High Phaeocystis globosa virus virophage (Pgvv) signatures were found and complete and partial Pgvv-like were obtained, which influence host⁻virus interactions. Our study expands existing knowledge of viral communities and their diversities from the Antarctic region and provides basic data for further exploring polar microbiomes.

A New Freshwater Cyanosiphovirus Harboring Integrase

Pelagic cyanobacteria are key players in the functioning of aquatic ecosystems, and their viruses (cyanophages) potentially affect the abundance and composition of cyanobacterial communities. Yet, there are few well-described freshwater cyanophages relative to their marine counterparts, and in general, few cyanosiphoviruses (family Siphoviridae) have been characterized, limiting our understanding of the biology and the ecology of this prominent group of viruses. Here, we characterize S-LBS1, a freshwater siphovirus lytic to a phycoerythrin-rich Synechococcus isolate (Strain TCC793). S-LBS1 has a narrow host range, a burst size of ∼400 and a relatively long infecting step before cell lysis occurs. It has a dsDNA 34,641 bp genome with putative genes for structure, DNA packing, lysis, replication, host interactions, DNA repair and metabolism. S-LBS1 is similar in genome size, genome architecture, and gene content, to previously described marine siphoviruses also infecting PE-rich Synechococcus, e.g., S-CBS1 and S-CBS3. However, unlike other Synechococcus phages, S-LBS1 encodes an integrase, suggesting its ability to establish lysogenic relationships with its host. Sequence recruitment from viral metagenomic data showed that S-LBS1-like viruses are diversely present in a wide range of aquatic environments, emphasizing their potential importance in controlling and structuring Synechococcus populations. A comparative analysis with 16 available sequenced cyanosiphoviruses reveals the absence of core genes within the genomes, suggesting high degree of genetic variability in siphoviruses infecting cyanobacteria. It is likely that cyanosiphoviruses have evolved as distinct evolutionary lineages and that adaptive co-evolution occurred between these viruses and their hosts (i.e., Synechococcus, Prochlorococcus, Nodularia, and Acaryochloris), constituting an important driving force for such phage diversification.

Pseudomonas PB1-Like Phages: Whole Genomes from Metagenomes Offer Insight into an Abundant Group of Bacteriophages

Despite the abundance, ubiquity and impact of environmental viruses, their inherent genomic plasticity and extreme diversity pose significant challenges for the examination of bacteriophages on Earth. Viral metagenomic studies have offered insight into broader aspects of phage ecology and repeatedly uncover genes to which we are currently unable to assign function. A combined effort of phage isolation and metagenomic survey of Chicago&rsquo;s nearshore waters of Lake Michigan revealed the presence of Pbunaviruses, relatives of the Pseudomonas phage PB1. This prompted our expansive investigation of PB1-like phages. Genomic signatures of PB1-like phages and Pbunaviruses were identified, permitting the unambiguous distinction between the presence/absence of these phages in soils, freshwater and wastewater samples, as well as publicly available viral metagenomic datasets. This bioinformatic analysis led to the de novo assembly of nine novel PB1-like phage genomes from a metagenomic survey of samples collected from Lake Michigan. While this study finds that Pbunaviruses are abundant in various environments of Northern Illinois, genomic variation also exists to a considerable extent within individual communities.

Diversity in a Polymicrobial Community Revealed by Analysis of Viromes, Endolysins and CRISPR Spacers

The polymicrobial biofilm communities in Mushroom and Octopus Spring in Yellowstone National Park (YNP) are well characterized, yet little is known about the phage populations. Dominant species, Synechococcus sp. JA-2-3B'a(2–13), Synechococcus sp. JA-3-3Ab, Chloroflexus sp. Y-400-fl, and Roseiflexus sp. RS-1, contain multiple CRISPR-Cas arrays, suggesting complex interactions with phage predators. To analyze phage populations from Octopus Spring biofilms, we sequenced a viral enriched fraction. To assemble and analyze phage metagenomic data, we developed a custom module, VIRITAS, implemented within the MetAMOS framework. This module bins contigs into groups based on tetranucleotide frequencies and CRISPR spacer-protospacer matching and ORF calling. Using this pipeline we were able to assemble phage sequences into contigs and bin them into three clusters that corroborated with their potential host range. The virome contained 52,348 predicted ORFs; some were clearly phage-like; 9319 ORFs had a recognizable Pfam domain while the rest were hypothetical. Of the recognized domains with CRISPR spacer matches, was the phage endolysin used by lytic phage to disrupt cells. Analysis of the endolysins present in the thermophilic cyanophage contigs revealed a subset of characterized endolysins as well as a Glyco_hydro_108 (PF05838) domain not previously associated with sequenced cyanophages. A search for CRISPR spacer matches to all identified phage endolysins demonstrated that a majority of endolysin domains were targets. This strategy provides a general way to link host and phage as endolysins are known to be widely distributed in bacteriophage. Endolysins can also provide information about host cell wall composition and have the additional potential to be used as targets for novel therapeutics.

Metagenomic recovery of phage genomes of uncultured freshwater actinobacteria

Low-GC Actinobacteria are among the most abundant and widespread microbes in freshwaters and have largely resisted all cultivation efforts. Consequently, their phages have remained totally unknown. In this work, we have used deep metagenomic sequencing to assemble eight complete genomes of the first tailed phages that infect freshwater Actinobacteria. Their genomes encode the actinobacterial-specific transcription factor whiB, frequently found in mycobacteriophages and also in phages infecting marine pelagic Actinobacteria. Its presence suggests a common and widespread strategy of modulation of host transcriptional machinery upon infection via this transcriptional switch. We present evidence that some whiB-carrying phages infect the acI lineage of Actinobacteria. At least one of them encodes the ADP-ribosylating component of the widespread bacterial AB toxins family (for example, clostridial toxin). We posit that the presence of this toxin reflects a 'trojan horse' strategy, providing protection at the population level to the abundant host microbes against eukaryotic predators.

Analysis of defence systems and a conjugative IncP-1 plasmid in the marine polyaromatic hydrocarbons-degrading bacterium Cycloclasticus sp. 78-ME

Marine prokaryotes have evolved a broad repertoire of defence systems to protect their genomes from lateral gene transfer including innate or acquired immune systems and infection-induced programmed cell suicide and dormancy. Here we report on the analysis of multiple defence systems present in the genome of the strain Cycloclasticus sp. 78-ME isolated from petroleum deposits of the tanker 'Amoco Milford Haven'. Cycloclasticus are ubiquitous bacteria globally important in polyaromatic hydrocarbons degradation in marine environments. Two 'defence islands' were identified in 78-ME genome: the first harbouring CRISPR-Cas with toxin-antitoxin system, while the second was composed by an array of genes for toxin-antitoxin and restriction-modification proteins. Among all identified spacers of CRISPR-Cas system only seven spacers match sequences of phages and plasmids. Furthermore, a conjugative plasmid p7ME01, which belongs to a new IncP-1θ ancestral archetype without any accessory mobile elements was found in 78-ME. Our results provide the context to the co-occurrence of diverse defence mechanisms in the genome of Cycloclasticus sp. 78-ME, which protect the genome of this highly specialized PAH-degrader. This study contributes to the further understanding of complex networks established in petroleum-based microbial communities.

Genomes of Abundant and Widespread Viruses from the Deep Ocean

The deep sea is a massive, largely oligotrophic ecosystem, stretched over nearly 65% of the planet’s surface. Deep-sea planktonic communities are almost completely dependent upon organic carbon sinking from the productive surface, forming a vital component of global biogeochemical cycles. However, despite their importance, viruses from the deep ocean remain largely unknown. Here, we describe the first complete genomes of deep-sea viruses assembled from metagenomic fosmid libraries. “Candidatus Pelagibacter” (SAR11) phage HTVC010P and Puniceispirillum phage HMO-2011 are considered the most abundant cultured marine viruses known to date. Remarkably, some of the viruses described here recruited as many reads from deep waters as these viruses do in the photic zone, and, considering the gigantic scale of the bathypelagic habitat, these genomes provide information about what could be some of the most abundant viruses in the world at large. Their role in the viral shunt in the global ocean could be very significant. Despite the challenges encountered in inferring the identity of their hosts, we identified one virus predicted to infect members of the globally distributed SAR11 cluster. We also identified a number of putative proviruses from diverse taxa, including deltaproteobacteria, bacteroidetes, SAR11, and gammaproteobacteria. Moreover, our findings also indicate that lysogeny is the preferred mode of existence for deep-sea viruses inhabiting an energy-limited environment, in sharp contrast to the predominantly lytic lifestyle of their photic-zone counterparts. Some of the viruses show a widespread distribution, supporting the tenet “everything is everywhere” for the deep-ocean virome.

The deep sea is among the largest known habitats and a critical cog in biogeochemical cycling but remains underexplored in its microbiology. Even more than is the case for its prokaryotic community, our knowledge of its viral component has remained limited by the paucity of information provided by studies dependent upon short sequence fragments. In this work, we attempt to fill this existing gap by using a combination of classical fosmid libraries with next-generation sequencing and assembly to recover long viral genomic fragments. We have sequenced ca. 6,000 fosmids from two metagenomics libraries made from prokaryotic biomass from the deep Mediterranean Sea and recovered twenty-eight complete viral genomes, all of them novel and quite distinct from all previously described viral genomes. They are preferentially found in deeper waters and are widely distributed all over the oceans. To our knowledge, this is the first report on complete and cosmopolitan viral genomes from the bathypelagic habitat.

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.

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.

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.

Combining genomic sequencing methods to explore viral diversity and reveal potential virus-host interactions

Viral diversity and virus-host interactions in oxygen-starved regions of the ocean, also known as oxygen minimum zones (OMZs), remain relatively unexplored. Microbial community metabolism in OMZs alters nutrient and energy flow through marine food webs, resulting in biological nitrogen loss and greenhouse gas production. Thus, viruses infecting OMZ microbes have the potential to modulate community metabolism with resulting feedback on ecosystem function. Here, we describe viral communities inhabiting oxic surface (10 m) and oxygen-starved basin (200 m) waters of Saanich Inlet, a seasonally anoxic fjord on the coast of Vancouver Island, British Columbia using viral metagenomics and complete viral fosmid sequencing on samples collected between April 2007 and April 2010. Of 6459 open reading frames (ORFs) predicted across all 34 viral fosmids, 77.6% (n = 5010) had no homology to reference viral genomes. These fosmids recruited a higher proportion of viral metagenomic sequences from Saanich Inlet than from nearby northeastern subarctic Pacific Ocean (Line P) waters, indicating differences in the viral communities between coastal and open ocean locations. While functional annotations of fosmid ORFs were limited, recruitment to NCBI's non-redundant “nr” database and publicly available single-cell genomes identified putative viruses infecting marine thaumarchaeal and SUP05 proteobacteria to provide potential host linkages with relevance to coupled biogeochemical cycling processes in OMZ waters. Taken together, these results highlight the power of coupled analyses of multiple sequence data types, such as viral metagenomic and fosmid sequence data with prokaryotic single cell genomes, to chart viral diversity, elucidate genomic and ecological contexts for previously unclassifiable viral sequences, and identify novel host interactions in natural and engineered ecosystems.

Marine cyanophages demonstrate biogeographic patterns throughout the global ocean

Myoviruses and podoviruses that infect cyanobacteria are the two major groups of marine cyanophages, but little is known of how their phylogenetic lineages are distributed in different habitats. In this study, we analyzed the phylogenetic relationships of cyanopodoviruses and cyanomyoviruses based on the existing genomes. The 28 cyanomyoviruses were classified into four clusters (I to IV), and 19 of the 20 cyanopodoviruses were classified into two clusters, MPP-A and MPP-B, with four subclusters within cluster MPP-B. These genomes were used to recruit cyanophage-like fragments from microbial and viral metagenomes to estimate the relative abundances of these cyanophage lineages. Our results showed that cyanopodoviruses and cyanomyoviruses are both abundant in various marine environments and that clusters MPP-B, II and III appear to be the most dominant lineages. Cyanopodoviruses and cluster I and IV cyanomyoviruses exhibited habitat-related variability in their relative levels of abundance, while cluster II and III cyanomyoviruses appeared to be consistently dominant in various habitats. Multivariate analyses showed that reads that mapped to Synechococcus phages and Prochlorococcus phages had distinct distribution patterns that were significantly correlated to those of Synechococcus and Prochlorococcus, respectively. The Mantel test also revealed a strong correlation between the community compositions of cyanophages and picocyanobacteria. Given that cyanomyoviruses tend to have a broad host range and some can cross-infect Synechococcus and Prochlorococcus, while cyanopodoviruses are commonly host specific, the observation that their community compositions both correlated significantly with that of picocyanobacteria was unexpected. Although cyanomyoviruses and cyanopodoviruses differ in host specificity, their biogeographic distributions are likely both constrained by the picocyanobacterial community.

Novel group of podovirus infecting the marine bacterium Alteromonas macleodii

Four novel, closely related podoviruses, which displayed lytic activity against the gamma-proteobacterium Alteromonas macleodii, have been isolated and sequenced. Alterophages AltAD45-P1 to P4 were obtained from water recovered near a fish farm in the Mediterranean Sea. Their morphology indicates that they belong to the Podoviridae. Their linear and dsDNA genomes are 100–104 kb in size, remarkably larger than any other described podovirus. The four AltAD45-phages share 99% nucleotide sequence identity over 97% of their ORFs, although an insertion was found in AltAD45-P1 and P2 and some regions were slightly more divergent. Despite the high overall sequence similarity among these four phages, the group with the insertion and the group without it, have different host ranges against the A. macleodii strains tested. The AltAD45-P1 to P4 phages have genes for DNA replication and transcription as well as structural genes, which are similar to the N4-like Podoviridae genus that is widespread in proteobacteria. However, in terms of their genomic structure, AltAD45-P1 to P4 differ from that of the N4-like phages. Some distinguishing features include the lack of a large virion encapsidated RNA polymerase gene, very well conserved among all the previously described N4-like phages, a single-stranded DNA binding protein and different tail protein genes. We conclude that the AltAD45 phages characterized in this study constitute a new genus within the Podoviridae.

Deciphering diversity and ecological function from marine metagenomes

Metagenomic sequencing now represents a common, powerful approach for investigating diversity and functional relationships in marine ecosystems. High-throughput datasets generated from random fragments of environmental DNA can provide a less biased view of organismal abundance (versus PCR-based amplicon sequencing) and enable novel exploration of microbial genomes by recovering genome assemblies from uncultured species, identifying ecological functions, and reconstructing metabolic pathways. This review highlights the current state of knowledge in marine metagenomics, focusing on biological insights gained from recent environmental studies and detailing commonly employed methods for data collection and analysis.

Tales from a thousand and one phages

The sequencing of marine metagenomic fosmids led to the discovery of several new complete phage genomes. Among the 21 major sequence groups, 10 totally novel groups of marine phages could be identified. Some of these represent the first phages infecting large marine prokaryotic phyla, such as the Verrucomicrobia and the recently described Ca. Actinomarinales. Coming from a single deep photic zone sample the diversity of phages found is astonishing, and the comparison with a metavirome from the same location indicates that only 2% of the real diversity was recovered. In addition to this large macro-diversity, rich micro-diversity was also found, affecting host-recognition modules, mirroring the variation of cell surface components in their host marine microbes.

Variations in abundance, genome size, morphology, and functional role of the virioplankton in Lakes Annecy and Bourget over a 1-year period

We sampled the surface waters (2-50 m) of two deep peri-alpine lakes over a 1-year period in order to examine (1) the abundance, vertical distribution, genome size, and morphology structures of the virioplankton; (2) the virus-mediated bacterial mortality; and (3) the specific genome size range of double-stranded DNA (dsDNA) phytoplankton viruses. Virus-like particle (VLP) concentrations varied between 4.16 × 10(7) (January) and 2.08 × 10(8) part mL(-1) (May) in Lake Bourget and between 2.7 × 10(7) (June) and 8.39 × 10(7) part mL(-1) (November) in Lake Annecy. Our flow cytometry analysis revealed at least three viral groups (referred to as virus-like particles 1, 2, and 3) that exhibited distinctive dynamics suggestive of different host types. Phage-induced bacterial mortality varied between 6.1% (June) and 33.2% (October) in Lake Bourget and between 7.4% (June) and 52.6% (November) in Lake Annecy, suggesting that viral lysis may be a key cause of mortality of the bacterioplankton. Virioplankton genome size ranged from 27 to 486 kb in Lake Bourget, while it reached 620 kb in Lake Annecy for which larger genome sizes were recorded. Our analysis of pulsed field gel electrophoresis bands using different PCR primers targeting both cyanophages and algal viruses showed that (1) dsDNA viruses infecting phytoplankton may range from 65 to 486 kb, and (2) both cyanophage and algal "diversity" were higher in Lake Annecy. Lakes Annecy and Bourget also differed regarding the proportions of both viral families (with the dominance of myoviruses vs. podoviruses) and infected bacterial morphotypes (short rods vs. elongated rods), in each of these lakes, respectively. Overall, our results reveal that (1) viruses displayed distinct temporal and vertical distribution, dynamics, community structure in terms of genome size and morphology, and viral activity in the two lakes; (2) the Myoviridae seemed to be the main cause of bacterial mortality in both lakes and this group seemed to be related to VLP2; and (3) phytoplankton viruses may have a broader range of genome size than previously thought. This study adds to growing evidence that viruses are diverse and play a significant role in freshwater microbial dynamics and more globally lake functioning. It highlights the importance of further considering this biological compartment for a better understanding of plankton ecology in peri-alpine lakes.

Expanding the marine virosphere using metagenomics

Viruses infecting prokaryotic cells (phages) are the most abundant entities of the biosphere and contain a largely uncharted wealth of genomic diversity. They play a critical role in the biology of their hosts and in ecosystem functioning at large. The classical approaches studying phages require isolation from a pure culture of the host. Direct sequencing approaches have been hampered by the small amounts of phage DNA present in most natural habitats and the difficulty in applying meta-omic approaches, such as annotation of small reads and assembly. Serendipitously, it has been discovered that cellular metagenomes of highly productive ocean waters (the deep chlorophyll maximum) contain significant amounts of viral DNA derived from cells undergoing the lytic cycle. We have taken advantage of this phenomenon to retrieve metagenomic fosmids containing viral DNA from a Mediterranean deep chlorophyll maximum sample. This method allowed description of complete genomes of 208 new marine phages. The diversity of these genomes was remarkable, contributing 21 genomic groups of tailed bacteriophages of which 10 are completely new. Sequence based methods have allowed host assignment to many of them. These predicted hosts represent a wide variety of important marine prokaryotic microbes like members of SAR11 and SAR116 clades, Cyanobacteria and also the newly described low GC Actinobacteria. A metavirome constructed from the same habitat showed that many of the new phage genomes were abundantly represented. Furthermore, other available metaviromes also indicated that some of the new phages are globally distributed in low to medium latitude ocean waters. The availability of many genomes from the same sample allows a direct approach to viral population genomics confirming the remarkable mosaicism of phage genomes.

Prevalence of viral photosynthetic and capsid protein genes from cyanophages in two large and deep perialpine lakes

Cyanophages are important components of aquatic ecosystems, but their genetic diversity has been little investigated in freshwaters. A yearlong survey was conducted in surface waters of the two largest natural perialpine lakes in France (Lake Annecy and Lake Bourget) to investigate part of this cyanophage diversity through the analysis of both structural (e.g., g20) and functional (e.g., psbA) genes. We found that these cyanophage signature genes were prevalent throughout the year but that the community compositions of g20 cyanomyoviruses were significantly different between the two lakes. In contrast, psbA-containing cyanophages seemed to be more similar between the two ecosystems. We also found that a large proportion of g20 sequences grouped with cyanomyophage isolates. psbA sequences, belonging to phages of Synechococcus spp., were characterized by distinct triplet motifs (with a novel viral triplet motif, EFE). Thus, our results show that cyanophages (i) are a diverse viral community in alpine lakes and (ii) are clearly distinct from some other freshwater and marine environments, suggesting the influence of unique biogeographic factors.

Assembly of a marine viral metagenome after physical fractionation

Metagenomic analyses of marine viruses generate an overview of viral genes present in a sample, but the percentage of the resulting sequence fragments that can be reassembled is low and the phenotype of the virus from which a given sequence derives is usually unknown. In this study, we employed physical fractionation to characterize the morphological and genomic traits of a subset of uncultivated viruses from a natural marine assemblage. Viruses from Kāne‘ohe Bay, Hawai‘i were fractionated by equilibrium buoyant density centrifugation in a cesium chloride (CsCl) gradient, and one fraction from the CsCl gradient was then further fractionated by strong anion-exchange chromatography. One of the fractions resulting from this two-dimensional separation appeared to be dominated by only a few virus types based on genome sizes and morphology. Sequences generated from a shotgun clone library of the viruses in this fraction were assembled into significantly more numerous contigs than have been generated with previous metagenomic investigations of whole DNA viral assemblages with comparable sequencing effort. Analysis of the longer contigs (up to 6.5 kb) assembled from our metagenome allowed us to assess gene arrangement in this subset of marine viruses. Our results demonstrate the potential for physical fractionation to facilitate sequence assembly from viral metagenomes and permit linking of morphological and genomic data for uncultivated viruses.