Gene network visualization and quantitative synteny analysis of more than 300 marine T4-like phage scaffolds from the GOS metagenome

Role of Phylogenetic Structure in the Dynamics of Coastal Viral Assemblages

Marine microbes, including viruses, are an essential part of the marine ecosystem, forming the base of the food web and driving biogeochemical cycles. Within this system, the composition of viral assemblages changes markedly with time, and some of these changes are repeatable through time; however, the extent to which these dynamics are reflected within versus among evolutionarily related groups of viruses is largely unexplored. To examine these dynamics, changes in the composition of two groups of ecologically important viruses and communities of their potential hosts were sampled every 2 weeks for 13 months at a coastal site in British Columbia, Canada. We sequenced two marker genes for viruses-the gene encoding the major capsid protein of T4-like phages and their relatives (gp23) and the RNA-dependent RNA polymerase (RdRp) gene of marnavirus-like RNA viruses-as well as marker genes for their bacterial and eukaryotic host communities, the genes encoding 16S rRNA and 18S rRNA. There were strong lagged correlations between viral diversity and community similarity of putative hosts, implying that the viruses influenced the composition of the host communities. The results showed that for both viral assemblages, the dominant clusters of phylogenetically related viruses shifted over time, and this was correlated with environmental changes. Viral clusters contained many ephemeral taxa and few persistent taxa, but within a viral assemblage, the ephemeral and persistent taxa were closely related, implying ecological dynamics within these clusters. Furthermore, these dynamics occurred in both the RNA and DNA viral assemblages surveyed, implying that this structure is common in natural viral assemblages.IMPORTANCE Viruses are major agents of microbial mortality in marine systems, yet little is known about changes in the composition of viral assemblages in relation to those of the microbial communities that they infect. Here, we sampled coastal seawater every 2 weeks for 1 year and used high-throughput sequencing of marker genes to follow changes in the composition of two groups of ecologically important viruses, as well as the communities of bacteria and protists that serve as their respective hosts. Different subsets of genetically related viruses dominated at different times. These results demonstrate that although the genetic composition of viral assemblages is highly dynamic temporally, for the most part the shuffling of genotypes occurs within a few clusters of phylogenetically related viruses. Thus, it appears that even in temperate coastal waters with large seasonal changes, the highly dynamic shuffling of viral genotypes occurs largely within a few subsets of related individuals.

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.

Understanding the enormous diversity of bacteriophages: the tailed phages that infect the bacterial family Enterobacteriaceae

Bacteriophages are the predominant biological entity on the planet. The recent explosion of sequence information has made estimates of their diversity possible. We describe the genomic comparison of 337 fully sequenced tailed phages isolated on 18 genera and 31 species of bacteria in the Enterobacteriaceae. These phages were largely unambiguously grouped into 56 diverse clusters (32 lytic and 24 temperate) that have syntenic similarity over >50% of the genomes within each cluster, but substantially less sequence similarity between clusters. Most clusters naturally break into sets of more closely related subclusters, 78% of which are correlated with their host genera. The largest groups of related phages are superclusters united by genome synteny to lambda (81 phages) and T7 (51 phages). This study forms a robust framework for understanding diversity and evolutionary relationships of existing tailed phages, for relating newly discovered phages and for determining host/phage relationships.

Diversity and genome dynamics of marine cyanophages using metagenomic analyses

Cyanophages are abundant in the oceanic environment and directly impact cyanobacterial distributions, physiological processes and evolution. Two samples collected from coastal Maine in July and September 2009 were enriched for Synechococcus cells using flow cytometry and examined through metagenomic sequencing. Homology-based sequence prediction indicated cyanophages, largely myoviruses, accounted for almost half the reads and provided insights into environmental infection events. T4-phage core-gene phylogenetic reconstruction revealed unique diversity among uncultured cyanophages and reference isolates resulting in identification of a new phylogenetic cluster. Genomic comparison of reference cyanophage strains S-SM2 and Syn1 with putative homologous contigs recovered from metagenomes provided evidence that gene insertion, deletion and recombination have occurred among, and are likely important for diversification of, natural populations. Identification of putative genetic exchange between cyanophage and non-cyanophage viruses, i.e. Micromonas virus and Pelagibacter phage, supports hypotheses related to a significant role for viruses in mediating transfer of genetic material between taxonomically diverse organisms with overlapping ecological niches.

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.

Phylogenomic network and comparative genomics reveal a diverged member of the ΦKZ-related group, marine vibrio phage ΦJM-2012

Bacteriophages are the largest reservoir of genetic diversity. Here we describe the novel phage ΦJM-2012. This natural isolate from marine Vibrio cyclitrophicus possesses very few gene contents relevant to other well-studied marine Vibrio phages. To better understand its evolutionary history, we built a mathematical model of pairwise relationships among 1,221 phage genomes, in which the genomes (nodes) are linked by edges representing the normalized number of shared orthologous protein families. This weighted network revealed that ΦJM-2012 was connected to only five members of the Pseudomonas ΦKZ-like phage family in an isolated network, strongly indicating that it belongs to this phage group. However, comparative genomic analyses highlighted an almost complete loss of colinearity with the ΦKZ-related genomes and little conservation of gene order, probably reflecting the action of distinct evolutionary forces on the genome of ΦJM-2012. In this phage, typical conserved core genes, including six RNA polymerase genes, were frequently displaced and the hyperplastic regions were rich in both unique genes and predicted unidirectional promoters with highly correlated orientations. Further, analysis of the ΦJM-2012 genome showed that segments of the conserved N-terminal parts of ΦKZ tail fiber paralogs exhibited evidence of combinatorial assortment, having switched transcriptional orientation, and there was recruitment and/or structural changes among phage endolysins and tail spike protein. Thus, this naturally occurring phage appears to have branched from a common ancestor of the ΦKZ-related groups, showing a distinct genomic architecture and unique genes that most likely reflect adaptation to its chosen host and environment.

The episodic evolution of fibritin: traces of ancient global environmental alterations may remain in the genomes of T4-like phages

The evolutionary adaptation of bacteriophages to their environment is achieved by alterations of their genomes involving a combination of both point mutations and lateral gene transfer. A phylogenetic analysis of a large set of collar fiber protein (fibritin) loci from diverse T4-like phages indicates that nearly all the modular swapping involving the C-terminal domain of this gene occurred in the distant past and has since ceased. In phage T4, this fibritin domain encodes the sequence that mediates both the attachment of the long tail fibers to the virion and also controls, in an environmentally sensitive way, the phage's ability to infect its host bacteria. Subsequent to its distant period of modular exchange, the evolution of fibritin has proceeded primarily by the slow vertical divergence mechanism. We suggest that ancient and sudden changes in the environment forced the T4-like phages to alter fibritin's mode of action or function. The genome's response to such episodes of rapid environmental change could presumably only be achieved quickly enough by employing the modular evolution mechanism. A phylogenetic analysis of the fibritin locus reveals the possible traces of such events within the T4 superfamily's genomes.

The global virome: not as big as we thought?

Viruses likely infect all organisms, serving to unknown extent as genetic vectors in complex networks of organisms. Environmental virologists have revealed that these abundant nanoscale entities are global players with critical roles in every ecosystem investigated. Curiously, novel genes dominate viral genomes and metagenomes, which has led to the suggestion that viruses represent the largest reservoir of unexplored genetic material on Earth with literature estimates, extrapolating from 14 mycobacteriophage genomes, suggesting that two billion phage-encoded ORFs remain to be discovered. Here we examine (meta)genomic data available in the decade since this provocative assertion, and use 'protein clusters' to evaluate whether sampling technologies have advanced to the point that we may be able to sample 'all' of viral diversity in nature.

Are we missing half of the viruses in the ocean?

Viruses are abundant in the ocean and a major driving force in plankton ecology and evolution. It has been assumed that most of the viruses in seawater contain DNA and infect bacteria, but RNA-containing viruses in the ocean, which almost exclusively infect eukaryotes, have never been quantified. We compared the total mass of RNA and DNA in the viral fraction harvested from seawater and using data on the mass of nucleic acid per RNA- or DNA-containing virion, estimated the abundances of each. Our data suggest that the abundance of RNA viruses rivaled or exceeded that of DNA viruses in samples of coastal seawater. The dominant RNA viruses in the samples were marine picorna-like viruses, which have small genomes and are at or below the detection limit of common fluorescence-based counting methods. If our results are typical, this means that counts of viruses and the rate measurements that depend on them, such as viral production, are significantly underestimated by current practices. As these RNA viruses infect eukaryotes, our data imply that protists contribute more to marine viral dynamics than one might expect based on their relatively low abundance. This conclusion is a departure from the prevailing view of viruses in the ocean, but is consistent with earlier theoretical predictions.

Oral T4-like phage cocktail application to healthy adult volunteers from Bangladesh

The genomic diversity of 99 T4-like coliphages was investigated by sequencing an equimolar mixture with Illumina technology and screening them against different databases for horizontal gene transfer and undesired genes. A 9-phage cocktail was given to 15 healthy adults from Bangladesh at a dose of 3×10(9) and 3×10(7) plaque-forming units and placebo respectively. Phages were detected in 64% of the stool samples when subjects were treated with higher titer phage, compared to 30% and 28% with lower-titer phage and placebo, respectively. No Escherichia coli was present in initial stool samples, and no amplification of phage was observed. One percent of the administered oral phage was recovered from the feces. No adverse events were observed by self-report, clinical examination, or from laboratory tests for liver, kidney, and hematology function. No impact of oral phage was seen on the fecal microbiota composition with respect to bacterial 16S rRNA from stool.

An ORFan no more: the bacteriophage T4 39.2 gene product, NwgI, modulates GroEL chaperone function

Bacteriophages are the most abundant biological entities in our biosphere, characterized by their hyperplasticity, mosaic composition, and the many unknown functions (ORFans) encoded by their immense genetic repertoire. These genes are potentially maintained by the bacteriophage to allow efficient propagation on hosts encountered in nature. To test this hypothesis, we devised a selection to identify bacteriophage-encoded gene(s) that modulate the host Escherichia coli GroEL/GroES chaperone machine, which is essential for the folding of certain host and bacteriophage proteins. As a result, we identified the bacteriophage RB69 gene 39.2, of previously unknown function and showed that homologs of 39.2 in bacteriophages T4, RB43, and RB49 similarly modulate GroEL/GroES. Production of wild-type bacteriophage T4 Gp39.2, a 58-amino-acid protein, (a) enables diverse bacteriophages to plaque on the otherwise nonpermissive groES or groEL mutant hosts in an allele-specific manner, (b) suppresses the temperature-sensitive phenotype of both groES and groEL mutants, (c) suppresses the defective UV-induced PolV function (UmuCD) of the groEL44 mutant, and (d) is lethal to the host when overproduced. Finally, as proof of principle that Gp39.2 is essential for bacteriophage growth on certain bacterial hosts, we constructed a T4 39.2 deletion strain and showed that, unlike the isogenic wild-type parent, it is incapable of propagating on certain groEL mutant hosts. We propose a model of how Gp39.2 modulates GroES/GroEL function.

Visualization and quality assessment of de novo genome assemblies

SUMMARY

Recent technological progress has greatly facilitated de novo genome sequencing. However, de novo assemblies consist in many pieces of contiguous sequence (contigs) arranged in thousands of scaffolds instead of small numbers of chromosomes. Confirming and improving the quality of such assemblies is critical for subsequent analysis. We present a method to evaluate genome scaffolding by aligning independently obtained transcriptome sequences to the genome and visually summarizing the alignments using the Cytoscape software. Applying this method to the genome of the red fire ant Solenopsis invicta allowed us to identify inconsistencies in 7%, confirm contig order in 20% and extend 16% of scaffolds.

CONTACT

oksana.ribagrognuz@unil.ch; yannick.wurm@unil.ch

AVAILABILITY

Scripts that generate tables for visualization in Cytoscape from FASTA sequence and scaffolding information files are publicly available at https://github.com/ksanao/TGNet.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The gp38 adhesins of the T4 superfamily: a complex modular determinant of the phage's host specificity

The tail fiber adhesins are the primary determinants of host range in the T4-type bacteriophages. Among the indispensable virion components, the sequences of the long tail fiber genes and their associated adhesins are among the most variable. The predominant form of the adhesin in the T4-type phages is not even the version of the gene encoded by T4, the archetype of the superfamily, but rather a small unrelated protein (gp38) encoded by closely related phages such as T2 and T6. This gp38 adhesin has a modular design: its N-terminal attachment domain binds at the tip of the tail fiber, whereas the C-terminal specificity domain determines its host receptor affinity. This specificity domain has a series of four hypervariable segments (HVSs) that are separated by a set of highly conserved glycine-rich motifs (GRMs) that apparently form the domain’s conserved structural core. The role of gp38’s various components was examined by a comparative analysis of a large series of gp38 adhesins from T-even superfamily phages with differing host specificities. A deletion analysis revealed that the individual HVSs and GRMs are essential to the T6 adhesin’s function and suggests that these different components all act in synergy to mediate adsorption. The evolutionary advantages of the modular design of the adhesin involving both conserved structural elements and multiple independent and easily interchanged specificity determinants are discussed.

Campylobacter jejuni group III phage CP81 contains many T4-like genes without belonging to the T4-type phage group: implications for the evolution of T4 phages

CP81 is a virulent Campylobacter group III phage whose linear genome comprises 132,454 bp. At the nucleotide level, CP81 differs from other phages. However, a number of its structural and replication/recombination proteins revealed a relationship to the group II Campylobacter phages CP220/CPt10 and to T4-type phages. Unlike the T4-related phages, the CP81 genome does not contain conserved replication and virion modules. Instead, the respective genes are scattered throughout the phage genome. Moreover, most genes for metabolic enzymes of CP220/CPt10 are lacking in CP81. On the other hand, the CP81 genome contains nine similar genes for homing endonucleases which may be involved in the attrition of the conserved gene order for the virion core genes of T4-type phages. The phage apparently possesses an unusual modification of C or G bases. Efficient cleavage of its DNA was only achieved with restriction enzymes recognizing pure A/T sites. Uncommonly, phenol extraction leads to a significant loss of CP81 DNA from the aqueous layer, a property not yet described for other phages belonging to the T4 superfamily.

Comparative metagenomics of microbial traits within oceanic viral communities

Viral genomes often contain genes recently acquired from microbes. In some cases (for example, psbA) the proteins encoded by these genes have been shown to be important for viral replication. In this study, using a unique search strategy on the Global Ocean Survey (GOS) metagenomes in combination with marine virome and microbiome pyrosequencing-based datasets, we characterize previously undetected microbial metabolic capabilities concealed within the genomes of uncultured marine viral communities. A total of 34 microbial gene families were detected on 452 viral GOS scaffolds. The majority of auxiliary metabolic genes found on these scaffolds have never been reported in phages. Host genes detected in viruses were mainly divided between genes encoding for different energy metabolism pathways, such as electron transport and newly identified photosystem genes, or translation and post-translation mechanism related. Our findings suggest previously undetected ways, in which marine phages adapt to their hosts and improve their fitness, including translation and post-translation level control over the host rather than the already known transcription level control.