The dilemma of phage taxonomy illustrated by comparative genomics of Sfi21-like Siphoviridae in lactic acid bacteria

The complete genome sequences of two dairy phages, Streptococcus thermophilus phage 7201 and Lactobacillus casei phage A2, are reported. Comparative genomics reveals that both phages are members of the recently proposed Sfi21-like genus of Siphoviridae, a widely distributed phage type in low-GC-content gram-positive bacteria. Graded relatedness, the hallmark of evolving biological systems, was observed when different Sfi21-like phages were compared. Across the structural module, the graded relatedness was represented by a high level of DNA sequence similarity or protein sequence similarity, or a shared gene map in the absence of sequence relatedness. This varying range of relatedness was found within Sfi21-like phages from a single species as demonstrated by the different prophages harbored by Lactococcus lactis strain IL1403. A systematic dot plot analysis with 11 complete L. lactis phage genome sequences revealed a clear separation of all temperate phages from two classes of virulent phages. The temperate lactococcal phages share DNA sequence homology in a patchwise fashion over the nonstructural gene cluster. With respect to structural genes, four DNA homology groups could be defined within temperate L. lactis phages. Closely related structural modules for all four DNA homology groups were detected in phages from Streptococcus or Listeria, suggesting that they represent distinct evolutionary lineages that have not uniquely evolved in L. lactis. It seems reasonable to base phage taxonomy on data from comparative genomics. However, the peculiar modular nature of phage evolution creates ambiguities in the definition of phage taxa by comparative genomics. For example, depending on the module on which the classification is based, temperate lactococcal phages can be classified as a single phage species, as four distinct phage species, or as two if not three different phage genera. We propose to base phage taxonomy on comparative genomics of a single structural gene module (head or tail genes). This partially phylogeny-based taxonomical system still mirrors some aspects of the current International Committee on Taxonomy in Virology classification system. In this system the currently sequenced lactococcal phages would be grouped into five genera: c2-, sk1, Sfi11-, r1t-, and Sfi21-like phages.     

Erratum: Causes and Consequences of Bacteriophage Diversification via Genetic Exchanges across Lifestyles and Bacterial Taxa

Bacteriophages (phages) evolve rapidly by acquiring genes from other phages. This results in mosaic genomes. Here, we identify numerous genetic transfers between distantly related phages and aim at understanding their frequency, consequences, and the conditions favoring them. Gene flow tends to occur between phages that are enriched for recombinases, transposases, and nonhomologous end joining, suggesting that both homologous and illegitimate recombination contribute to gene flow. Phage family and host phyla are strong barriers to gene exchange, but phage lifestyle is not. Even if we observe four times more recent transfers between temperate phages than between other pairs, there is extensive gene flow between temperate and virulent phages, and between the latter. These predominantly involve virulent phages with large genomes previously classed as low gene flux, and lead to the preferential transfer of genes encoding functions involved in cell energetics, nucleotide metabolism, DNA packaging and injection, and virion assembly. Such exchanges may contribute to the observed twice larger genomes of virulent phages. We used genetic transfers, which occur upon coinfection of a host, to compare phage host range. We found that virulent phages have broader host ranges and can mediate genetic exchanges between narrow host range temperate phages infecting distant bacterial hosts, thus contributing to gene flow between virulent phages, as well as between temperate phages. This gene flow drastically expands the gene repertoires available for phage and bacterial evolution, including the transfer of functional innovations across taxa.     

Evolutionary history of bacteriophages with double-stranded DNA genomes

Background

Reconstruction of evolutionary history of bacteriophages is a difficult problem because of fast sequence drift and lack of omnipresent genes in phage genomes. Moreover, losses and recombinational exchanges of genes are so pervasive in phages that the plausibility of phylogenetic inference in phage kingdom has been questioned.

Results

We compiled the profiles of presence and absence of 803 orthologous genes in 158 completely sequenced phages with double-stranded DNA genomes and used these gene content vectors to infer the evolutionary history of phages. There were 18 well-supported clades, mostly corresponding to accepted genera, but in some cases appearing to define new taxonomic groups. Conflicts between this phylogeny and trees constructed from sequence alignments of phage proteins were exploited to infer 294 specific acts of intergenome gene transfer.

Conclusion

A notoriously reticulate evolutionary history of fast-evolving phages can be reconstructed in considerable detail by quantitative comparative genomics.

Open peer review

This article was reviewed by Eugene Koonin, Nicholas Galtier and Martijn Huynen.     

Evolutionary genomics of a temperate bacteriophage in an obligate intracellular bacteria (Wolbachia)

Genome evolution of bacteria is usually influenced by ecology, such that bacteria with a free-living stage have large genomes and high rates of horizontal gene transfer, while obligate intracellular bacteria have small genomes with typically low amounts of gene exchange. However, recent studies indicate that obligate intracellular species that host-switch frequently harbor agents of horizontal transfer such as mobile elements. For example, the temperate double-stranded DNA bacteriophage WO in Wolbachia persistently transfers between bacterial coinfections in the same host. Here we show that despite the phage's rampant mobility between coinfections, the prophage's genome displays features of constraint related to its intracellular niche. First, there is always at least one intact prophage WO and usually several degenerate, independently-acquired WO prophages in each Wolbachia genome. Second, while the prophage genomes are modular in composition with genes of similar function grouping together, the modules are generally not interchangeable with other unrelated phages and thus do not evolve by the Modular Theory. Third, there is an unusual core genome that strictly consists of head and baseplate genes; other gene modules are frequently deleted. Fourth, the prophage recombinases are diverse and there is no conserved integration sequence. Finally, the molecular evolutionary forces acting on prophage WO are point mutation, intragenic recombination, deletion, and purifying selection. Taken together, these analyses indicate that while lateral transfer of phage WO is pervasive between Wolbachia with occasional new gene uptake, constraints of the intracellular niche obstruct extensive mixture between WO and the global phage population. Although the Modular Theory has long been considered the paradigm of temperate bacteriophage evolution in free-living bacteria, it appears irrelevant in phages of obligate intracellular bacteria.     

Genome characteristics of a novel phage from Bacillus thuringiensis showing high similarity with phage from Bacillus cereus

Bacillus thuringiensis is an important entomopathogenic bacterium belongs to the Bacillus cereus group, which also includes B. anthracis and B. cereus. Several genomes of phages originating from this group had been sequenced, but no genome of Siphoviridae phage from B. thuringiensis has been reported. We recently sequenced and analyzed the genome of a novel phage, BtCS33, from a B. thuringiensis strain, subsp. kurstaki CS33, and compared the gneome of this phage to other phages of the B. cereus group. BtCS33 was the first Siphoviridae phage among the sequenced B. thuringiensis phages. It produced small, turbid plaques on bacterial plates and had a narrow host range. BtCS33 possessed a linear, double-stranded DNA genome of 41,992 bp with 57 putative open reading frames (ORFs). It had a typical genome structure consisting of three modules: the "late" region, the "lysogeny-lysis" region and the "early" region. BtCS33 exhibited high similarity with several phages, B. cereus phage Wβ and some variants of Wβ, in genome organization and the amino acid sequences of structural proteins. There were two ORFs, ORF22 and ORF35, in the genome of BtCS33 that were also found in the genomes of B. cereus phage Wβ and may be involved in regulating sporulation of the host cell. Based on these observations and analysis of phylogenetic trees, we deduced that B. thuringiensis phage BtCS33 and B. cereus phage Wβ may have a common distant ancestor.     

Mycobacteriophage Marvin: a new singleton phage with an unusual genome organization

Mycobacteriophages represent a genetically diverse group of viruses that infect mycobacterial hosts. Although more than 80 genomes have been sequenced, these still poorly represent the likely diversity of the broader population of phages that can infect the host, Mycobacterium smegmatis mc(2)155. We describe here a newly discovered phage, Marvin, which is a singleton phage, having no previously identified close relatives. The 65,100-bp genome contains 107 predicted protein-coding genes arranged in a noncanonical genomic architecture in which a subset of the minor tail protein genes are displaced about 20 kbp from their typical location, situated among nonstructural genes anticipated to be expressed early in lytic growth. Marvin is not temperate, and stable lysogens cannot be recovered from infections, although the presence of a putative xis gene suggests that Marvin could be a relatively recent derivative of a temperate parent. The Marvin genome is replete with novel genes not present in other mycobacteriophage genomes, and although most are of unknown function, the presence of amidoligase and glutamine amidotransferase genes suggests intriguing possibilities for the interactions of Marvin with its mycobacterial hosts.     

Cluster J Mycobacteriophages: Intron Splicing in Capsid and Tail Genes

Bacteriophages isolated on Mycobacterium smegmatis mc²155 represent many distinct genomes sharing little or no DNA sequence similarity. The genomes are architecturally mosaic and are replete with genes of unknown function. A new group of genomes sharing substantial nucleotide sequences constitute Cluster J. The six mycobacteriophages forming Cluster J are morphologically members of the Siphoviridae, but have unusually long genomes ranging from 106.3 to 117 kbp. Reconstruction of the capsid by cryo-electron microscopy of mycobacteriophage BAKA reveals an icosahedral structure with a triangulation number of 13. All six phages are temperate and homoimmune, and prophage establishment involves integration into a tRNA-Leu gene not previously identified as a mycobacterial attB site for phage integration. The Cluster J genomes provide two examples of intron splicing within the virion structural genes, one in a major capsid subunit gene, and one in a tail gene. These genomes also contain numerous free-standing HNH homing endonuclease, and comparative analysis reveals how these could contribute to genome mosaicism. The unusual Cluster J genomes provide new insights into phage genome architecture, gene function, capsid structure, gene mobility, intron splicing, and evolution.     

Properties of the new D3-like Pseudomonas aeruginosa bacteriophage phiPMG1: Genome structure and prospects for the use in phage therapy

Results of studying the novel virulent phage phiPMG1 active on Pseudomonas aeruginosa are presented. It is shown that phiPMG1 exhibits significant homology and the similarity in the overall structure with the genome of a temperate phage converts D3. Phage phiPMG1 differs from D3 in that it fails to stably lysogenize bacteria and can grow on strains carrying plasmids that cause growth inhibition of phage D3 and some other phages. This significantly diminishes the probability of horizontal gene transfer with phage phiPMG1 and suggests the possible employment of this phage in phage therapy. A comparison of phages phiPMG1 and D3 structures of genomes in demonstrated not only high homology of 65 genes, but also the presence of 16 genes in the phiPMG1 genome that were not included in the in NCBI database. Apparently, the evolution of genomes in phages of this species is mostly associated with migrations into other species of bacteria, and recombinations with phages of other species (for example, F116). A detailed analysis of structure of one region genomes, which significant nonhomology for the three D3-like phages (D3, phiPMG1 and PAJU2), revealed that the phiPMG1 genome possible closest to a hypothetical genome of ancestral phage of this species.     

Marine phage genomics

Marine phages are the most abundant biological entities in the oceans. They play important roles in carbon cycling through marine food webs, gene transfer by transduction and conversion of hosts by lysogeny. The handful of marine phage genomes that have been sequenced to date, along with prophages in marine bacterial genomes, and partial sequencing of uncultivated phages are yielding glimpses of the tremendous diversity and physiological potential of the marine phage community. Common gene modules in diverse phages are providing the information necessary to make evolutionary comparisons. Finally, deciphering phage genomes is providing clues about the adaptive response of phages and their hosts to environmental cues.     

Codon Bias is a Major Factor Explaining Phage Evolution in Translationally Biased Hosts

The size and diversity of bacteriophage populations require methodologies to quantitatively study the landscape of phage differences. Statistical approaches are confronted with small genome sizes forbidding significant single-phage analysis, and comparative methods analyzing full phage genomes represent an alternative but they are of difficult interpretation due to lateral gene transfer, which creates a mosaic spectrum of related phage species. Based on a large-scale codon bias analysis of 116 DNA phages hosted by 11 translationally biased bacteria belonging to different phylogenetic families, we observe that phage genomes are almost always under codon selective pressure imposed by translationally biased hosts, and we propose a classification of phages with translationally biased hosts which is based on adaptation patterns. We introduce a computational method for comparing phages sharing homologous proteins, possibly accepted by different hosts. We observe that throughout phages, independently from the host, capsid genes appear to be the most affected by host translational bias. For coliphages, genes involved in virion morphogenesis, host interaction and ssDNA binding are also affected by adaptive pressure. Adaptation affects long and small phages in a significant way. We analyze in more detail the Microviridae phage space to illustrate the potentiality of the approach. The small number of directions in adaptation observed in phages grouped around ϕX174 is discussed in the light of functional bias. The adaptation analysis of the set of Microviridae phages defined around ϕMH2K shows that phage classification based on adaptation does not reflect bacterial phylogeny. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Molecular characterization and comparison of 38 virulent and temperate bacteriophages of Streptococcus lactis

Thirty-three virulent and five temperate phages of Streptococcus lactis and Streptococcus cremoris were differentiated into three groups by DNA homology. A complete lack of DNA homology was demonstrated between the phage groups. Within each group, large parts of the phage genomes were homologous except for a few phages. One group consisted of five temperate and two virulent phages suggesting that virulent phages isolated during abnormal fermentations and temperate phages isolated after induction from lactic streptococcal starter cultures may be related to one another. We observed a good correlation between the grouping of phages by DNA homology and by their protein composition since within the same DNA homology group, the protein composition of a phage exhibited some similarities with that of the other phages of the group. Therefore, the DNA homologies seemed to be located, at least, in the region coding for the structural proteins. By immunoblotting, we confirmed the relatedness between the proteins of the phages belonging to the same DNA homology group. The important host range factor in bacterium-phage interactions appeared to be an unreliable criterion in determining phage taxonomy.     

Characterization of Novel Phages Isolated in Coagulase-Negative Staphylococci Reveals Evolutionary Relationships with Staphylococcus aureus Phages

Despite increasing interest in coagulase-negative staphylococci (CoNS), little information is available about their bacteriophages. We isolated and sequenced three novel temperate Siphoviridae phages (StB12, StB27, and StB20) from the CoNS Staphylococcus hominis and S. capitis species. The genome sizes are around 40 kb, and open reading frames (ORFs) are arranged in functional modules encoding lysogeny, DNA metabolism, morphology, and cell lysis. Bioinformatics analysis allowed us to assign a potential function to half of the predicted proteins. Structural elements were further identified by proteomic analysis of phage particles, and DNA-packaging mechanisms were determined. Interestingly, the three phages show identical integration sites within their host genomes. In addition to this experimental characterization, we propose a novel classification based on the analysis of 85 phage and prophage genomes, including 15 originating from CoNS. Our analysis established 9 distinct clusters and revealed close relationships between S. aureus and CoNS phages. Genes involved in DNA metabolism and lysis and potentially in phage-host interaction appear to be widespread, while structural genes tend to be cluster specific. Our findings support the notion of a possible reciprocal exchange of genes between phages originating from S. aureus and CoNS, which may be of crucial importance for pathogenesis in staphylococci.     

Characterization of the Genome of the Dairy Lactobacillus helveticus Bacteriophage ΦAQ113

The complete genomic sequence of the dairy Lactobacillus helveticus bacteriophage ΦAQ113 was determined. Phage ΦAQ113 is a Myoviridae bacteriophage with an isometric capsid and a contractile tail. The final assembled consensus sequence revealed a linear, circularly permuted, double-stranded DNA genome with a size of 36,566 bp and a G+C content of 37%. Fifty-six open reading frames (ORFs) were predicted, and a putative function was assigned to approximately 90% of them. The ΦAQ113 genome shows functionally related genes clustered together in a genome structure composed of modules for DNA replication/regulation, DNA packaging, head and tail morphogenesis, cell lysis, and lysogeny. The identification of genes involved in the establishment of lysogeny indicates that it may have originated as a temperate phage, even if it was isolated from natural cheese whey starters as a virulent phage, because it is able to propagate in a sensitive host strain. Additionally, we discovered that the ΦAQ113 phage genome is closely related to Lactobacillus gasseri phage KC5a and Lactobacillus johnsonii phage Lj771 genomes. The phylogenetic similarities between L. helveticus phage ΦAQ113 and two phages that belong to gut species confirm a possible common ancestral origin and support the increasing consideration of L. helveticus as a health-promoting organism.     

Detection of novel recombinases in bacteriophage genomes unveils Rad52, Rad51 and Gp2.5 remote homologs

Homologous recombination is a key in contributing to bacteriophages genome repair, circularization and replication. No less than six kinds of recombinase genes have been reported so far in bacteriophage genomes, two (UvsX and Gp2.5) from virulent, and four (Sak, Redβ, Erf and Sak4) from temperate phages. Using profile–profile comparisons, structure-based modelling and gene-context analyses, we provide new views on the global landscape of recombinases in 465 bacteriophages. We show that Sak, Redβ and Erf belong to a common large superfamily adopting a shortcut Rad52-like fold. Remote homologs of Sak4 are predicted to adopt a shortcut Rad51/RecA fold and are discovered widespread among phage genomes. Unexpectedly, within temperate phages, gene-context analyses also pinpointed the presence of distant Gp2.5 homologs, believed to be restricted to virulent phages. All in all, three major superfamilies of phage recombinases emerged either related to Rad52-like, Rad51-like or Gp2.5-like proteins. For two newly detected recombinases belonging to the Sak4 and Gp2.5 families, we provide experimental evidence of their recombination activity in vivo. Temperate versus virulent lifestyle together with the importance of genome mosaicism is discussed in the light of these novel recombinases. Screening for these recombinases in genomes can be performed at http://biodev.extra.cea.fr/virfam.     

Bioinformatic analysis of the Acinetobacter baumannii phage AB1 genome

As one of the pathogens of hospital-acquired infections, Acinetobacter baumannii poses great challenges to the public health. A. baumannii phage could be an effective way to fight multi-resistant A. baumannii. Here, we completed the whole genome sequencing of the complete genome of A. baumannii phage AB1, which consists of 45,159bp and is a double-stranded DNA molecule with an average GC content of 37.7%. The genome encodes one tRNA gene and 85 open reading frames (ORFs) and the average size of the ORF is 531bp in length. Among 85 ORFs, only 14 have been identified to share significant sequence similarities to the genes with known functions, while 28 are similar in sequence to the genes with function-unknown genes in the database and 43 ORFs are uniquely present in the phage AB1 genome. Fourteen function-assigned genes with putative functions include five phage structure proteins, an RNA polymerase, a big sub-unit and a small sub-unit of a terminase, a methylase and a recombinase and the proteins involved in DNA replication and so on. Multiple sequence alignment was conducted among those homologous proteins and the phylogenetic trees were reconstructed to analyze the evolutionary courses of these essential genes. From comparative genomics analysis, it turned out clearly that the frame of the phage genome mainly consisted of genes from Xanthomonas phages, Burkholderia ambifaria phages and Enterobacteria phages and while it comprises genes of its host A. baumannii only sporadically. The mosaic feature of the phage genome suggested that the horizontal gene transfer occurred among the phage genomes and between the phages and the host bacterium genomes. Analyzing the genome sequences of the phages should lay sound foundation to investigate how phages adapt to the environment and infect their hosts, and even help to facilitate the development of biological agents to deal with pathogenic bacteria.     

The annotated complete DNA sequence of Enterococcus faecalis bacteriophage φEf11 and its comparison with all available phage and predicted prophage genomes

φEf11 is a temperate Siphoviridae bacteriophage isolated by induction from a lysogenic Enterococcus faecalis strain. The φEf11 DNA was completely sequenced and found to be 42,822 bp in length, with a G+C mol% of 34.4%. Genome analysis revealed 65 ORFs, accounting for 92.8% of the DNA content. All except for seven of the ORFs displayed sequence similarities to previously characterized proteins. The genes were arranged in functional modules, organized similar to that of several other phages of low GC Gram-positive bacteria; however, the number and arrangement of lysis-related genes were atypical of these bacteriophages. A 159 bp noncoding region between predicted cI and cro genes is highly similar to the functionally characterized early promoter region of lactococcal temperate phage TP901-1, and possessed a predicted stem-loop structure in between predicted P(L) and P(R) promoters, suggesting a novel mechanism of repression of these two bacteriophages from the λ paradigm. Comparison with all available phage and predicted prophage genomes revealed that the φEf11 genome displays unique features, suggesting that φEf11 may be a novel member of a larger family of temperate prophages that also includes lactococcal phages. Trees based on the blast score ratio grouped this family by tail fiber similarity, suggesting that these trees are useful for identifying phages with similar tail fibers.     

Analysis of six prophages in Lactococcus lactis IL1403: different genetic structure of temperate and virulent phage populations

We report the genetic organisation of six prophages present in the genome of Lactococcus lactis IL1403. The three larger prophages (36–42 kb), belong to the already described P335 group of temperate phages, whereas the three smaller ones (13–15 kb) are most probably satellites relying on helper phage(s) for multiplication. These data give a new insight into the genetic structure of lactococcal phage populations. P335 temperate phages have variable genomes, sharing homology over only 10–33% of their length. In contrast, virulent phages have highly similar genomes sharing homology over >90% of their length. Further analysis of genetic structure in all known groups of phages active on other bacterial hosts such as Escherichia coli, Bacillus subtilis, Mycobacterium and Streptococcus thermophilus confirmed the existence of two types of genetic structure related to the phage way of life. This might reflect different intensities of horizontal DNA exchange: low among purely virulent phages and high among temperate phages and their lytic homologues. We suggest that the constraints on genetic exchange among purely virulent phages reflect their optimal genetic organisation, adapted to a more specialised and extreme form of parasitism than temperate/lytic phages.     

Ecogenomics and genome landscapes of marine Pseudoalteromonas phage H105/1

Marine phages have an astounding global abundance and ecological impact. However, little knowledge is derived from phage genomes, as most of the open reading frames in their small genomes are unknown, novel proteins. To infer potential functional and ecological relevance of sequenced marine Pseudoalteromonas phage H105/1, two strategies were used. First, similarity searches were extended to include six viral and bacterial metagenomes paired with their respective environmental contextual data. This approach revealed ‘ecogenomic'' patterns of Pseudoalteromonas phage H105/1, such as its estuarine origin. Second, intrinsic genome signatures (phylogenetic, codon adaptation and tetranucleotide (tetra) frequencies) were evaluated on a resolved intra-genomic level to shed light on the evolution of phage functional modules. On the basis of differential codon adaptation of Phage H105/1 proteins to the sequenced Pseudoalteromonas spp., regions of the phage genome with the most ‘host''-adapted proteins also have the strongest bacterial tetra signature, whereas the least ‘host''-adapted proteins have the strongest phage tetra signature. Such a pattern may reflect the evolutionary history of the respective phage proteins and functional modules. Finally, analysis of the structural proteome identified seven proteins that make up the mature virion, four of which were previously unknown. This integrated approach combines both novel and classical strategies and serves as a model to elucidate ecological inferences and evolutionary relationships from phage genomes that typically abound with unknown gene content.     

Complete genome sequence and analysis of the<Emphasis Type="Italic">Streptomyces aureofaciens</Emphasis> phage μ1/6

The entire double-stranded DNA genome of the Streptomyces aureofaciens phage mu1/6 was sequenced and analyzed. Its size is 38.194 kbp with an overall molar G+C content of 71.2 %. Fifty-two potential open reading frames (orfs) were identified, divided into two oppositely transcribed regions. In the left arm of the mu1/6 genome, an identified putative integrase and possible regulation proteins were identified. The rightwards transcribed region contains genes organized into apparently four functional units responsible for: (i) replication, (ii) DNA packaging and head assembly, (iii) tail morphogenesis, and (iv) lysis. Putative functions were assigned to twelve orfs based on bioinformatic analysis or experimental substantiation. Comparative analysis with three complete genomes of streptomycete phages revealed resemblance with respect to the organization of their genes into functional modules. Closer relationship was observed only between mu1/6 and S. venezuelae phage VWB.