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.     

Prophage Tracer: precisely tracing prophages in prokaryotic genomes using overlapping split-read alignment

The life cycle of temperate phages includes a lysogenic cycle stage when the phage integrates into the host genome and becomes a prophage. However, the identification of prophages that are highly divergent from known phages remains challenging. In this study, by taking advantage of the lysis-lysogeny switch of temperate phages, we designed Prophage Tracer, a tool for recognizing active prophages in prokaryotic genomes using short-read sequencing data, independent of phage gene similarity searching. Prophage Tracer uses the criterion of overlapping split-read alignment to recognize discriminative reads that contain bacterial (attB) and phage (attP) att sites representing prophage excision signals. Performance testing showed that Prophage Tracer could predict known prophages with precise boundaries, as well as novel prophages. Two novel prophages, dsDNA and ssDNA, encoding highly divergent major capsid proteins, were identified in coral-associated bacteria. Prophage Tracer is a reliable data mining tool for the identification of novel temperate phages and mobile genetic elements. The code for the Prophage Tracer is publicly available at https://github.com/WangLab-SCSIO/Prophage_Tracer.     

The rhizome of Reclinomonas americana, Homo sapiens, Pediculus humanus and Saccharomyces cerevisiae mitochondria

Background

Period 10 dinucleotides are structurally and functionally validated factors that influence the ability of DNA to form nucleosomes, histone core octamers. Robust identification of periodic signals in DNA sequences is therefore required to understand nucleosome organisation in genomes. While various techniques for identifying periodic components in genomic sequences have been proposed or adopted, the requirements for such techniques have not been considered in detail and confirmatory testing for a priori specified periods has not been developed.

Results

We compared the estimation accuracy and suitability for confirmatory testing of autocorrelation, discrete Fourier transform (DFT), integer period discrete Fourier transform (IPDFT) and a previously proposed Hybrid measure. A number of different statistical significance procedures were evaluated but a blockwise bootstrap proved superior. When applied to synthetic data whose period-10 signal had been eroded, or for which the signal was approximately period-10, the Hybrid technique exhibited superior properties during exploratory period estimation. In contrast, confirmatory testing using the blockwise bootstrap procedure identified IPDFT as having the greatest statistical power. These properties were validated on yeast sequences defined from a ChIP-chip study where the Hybrid metric confirmed the expected dominance of period-10 in nucleosome associated DNA but IPDFT identified more significant occurrences of period-10. Application to the whole genomes of yeast and mouse identified ~ 21% and ~ 19% respectively of these genomes as spanned by period-10 nucleosome positioning sequences (NPS).

Conclusions

For estimating the dominant period, we find the Hybrid period estimation method empirically to be the most effective for both eroded and approximate periodicity. The blockwise bootstrap was found to be effective as a significance measure, performing particularly well in the problem of period detection in the presence of eroded periodicity. The autocorrelation method was identified as poorly suited for use with the blockwise bootstrap. Application of our methods to the genomes of two model organisms revealed a striking proportion of the yeast and mouse genomes are spanned by NPS. Despite their markedly different sizes, roughly equivalent proportions (19-21%) of the genomes lie within period-10 spans of the NPS dinucleotides {AA, TT, TA}. The biological significance of these regions remains to be demonstrated. To facilitate this, the genomic coordinates are available as Additional files 1, 2, and 3 in a format suitable for visualisation as tracks on popular genome browsers.

Reviewers

This article was reviewed by Prof Tomas Radivoyevitch, Dr Vsevolod Makeev (nominated by Dr Mikhail Gelfand), and Dr Rob D Knight.     

Analysis of High-Throughput Sequencing and Annotation Strategies for Phage Genomes

Background

Bacterial viruses (phages) play a critical role in shaping microbial populations as they influence both host mortality and horizontal gene transfer. As such, they have a significant impact on local and global ecosystem function and human health. Despite their importance, little is known about the genomic diversity harbored in phages, as methods to capture complete phage genomes have been hampered by the lack of knowledge about the target genomes, and difficulties in generating sufficient quantities of genomic DNA for sequencing. Of the approximately 550 phage genomes currently available in the public domain, fewer than 5% are marine phage.

Methodology/Principal Findings

To advance the study of phage biology through comparative genomic approaches we used marine cyanophage as a model system. We compared DNA preparation methodologies (DNA extraction directly from either phage lysates or CsCl purified phage particles), and sequencing strategies that utilize either Sanger sequencing of a linker amplification shotgun library (LASL) or of a whole genome shotgun library (WGSL), or 454 pyrosequencing methods. We demonstrate that genomic DNA sample preparation directly from a phage lysate, combined with 454 pyrosequencing, is best suited for phage genome sequencing at scale, as this method is capable of capturing complete continuous genomes with high accuracy. In addition, we describe an automated annotation informatics pipeline that delivers high-quality annotation and yields few false positives and negatives in ORF calling.

Conclusions/Significance

These DNA preparation, sequencing and annotation strategies enable a high-throughput approach to the burgeoning field of phage genomics.     

Protective effect of probiotics on Salmonella infectivity assessed with combined in vitro gut fermentation-cellular models

Background

The alphaproteobacterium Wolbachia pipientis, the most common endosymbiont in eukaryotes, is found predominantly in insects including many Drosophila species. Although Wolbachia is primarily vertically transmitted, analysis of its genome provides evidence for frequent horizontal transfer, extensive recombination and numerous mobile genetic elements. The genome sequence of Wolbachia in Drosophila simulans Riverside (wRi) is available along with the integrated bacteriophages, enabling a detailed examination of phage genes and the role of these genes in the biology of Wolbachia and its host organisms. Wolbachia is widely known for its ability to modify the reproductive patterns of insects. One particular modification, cytoplasmic incompatibility, has previously been shown to be dependent on Wolbachia density and inversely related to the titer of lytic phage. The wRi genome has four phage regions, two WORiBs, one WORiA and one WORiC.

Results

In this study specific primers were designed to distinguish between these four prophage types in wRi, and quantitative PCR was used to measure the titer of bacteriophages in testes, ovaries, embryos and adult flies. In all tissues tested, WORiA and WORiB were not found to be present in excess of their integrated prophages; WORiC, however, was found to be present extrachromosomally. WORiC is undergoing extrachromosomal replication in wRi. The density of phage particles was found to be consistent in individual larvae in a laboratory population. The WORiC genome is organized in conserved blocks of genes and aligns most closely with other known lytic WO phages, WOVitA and WOCauB.

Conclusions

The results presented here suggest that WORiC is the lytic form of WO in D. simulans, is undergoing extrachromosomal replication in wRi, and belongs to a conserved family of phages in Wolbachia.     

DNA repeat arrays in chicken and human genomes and the adaptive evolution of avian genome size

Background

Birds have smaller average genome sizes than other tetrapod classes, and it has been proposed that a relatively low frequency of repeating DNA is one factor in reduction of avian genome sizes.

Results

DNA repeat arrays in the sequenced portion of the chicken (Gallus gallus) autosomes were quantified and compared with those in human autosomes. In the chicken 10.3% of the genome was occupied by DNA repeats, in contrast to 44.9% in human. In the chicken, the percentage of a chromosome occupied by repeats was positively correlated with chromosome length, but even the largest chicken chromosomes had repeat densities much lower than those in human, indicating that avoidance of repeats in the chicken is not confined to minichromosomes. When 294 simple sequence repeat types shared between chicken and human genomes were compared, mean repeat array length and maximum repeat array length were significantly lower in the chicken than in human.

Conclusions

The fact that the chicken simple sequence repeat arrays were consistently smaller than arrays of the same type in human is evidence that the reduction in repeat array length in the chicken has involved numerous independent evolutionary events. This implies that reduction of DNA repeats in birds is the result of adaptive evolution. Reduction of DNA repeats on minichromosomes may be an adaptation to permit chiasma formation and alignment of small chromosomes. However, the fact that repeat array lengths are consistently reduced on the largest chicken chromosomes supports the hypothesis that other selective factors are at work, presumably related to the reduction of cell size and consequent advantages for the energetic demands of flight.     

Characterization of Paenibacillus larvae bacteriophages and their genomic relationships to firmicute bacteriophages

Background

Paenibacillus larvae is a Firmicute bacterium that causes American Foulbrood, a lethal disease in honeybees and is a major source of global agricultural losses. Although P. larvae phages were isolated prior to 2013, no full genome sequences of P. larvae bacteriophages were published or analyzed. This report includes an in-depth analysis of the structure, genomes, and relatedness of P. larvae myoviruses Abouo, Davis, Emery, Jimmer1, Jimmer2, and siphovirus phiIBB_Pl23 to each other and to other known phages.

Results

P. larvae phages Abouo, Davies, Emery, Jimmer1, and Jimmer2 are myoviruses with ~50 kbp genomes. The six P. larvae phages form three distinct groups by dotplot analysis. An annotated linear genome map of these six phages displays important identifiable genes and demonstrates the relationship between phages. Sixty phage assembly or structural protein genes and 133 regulatory or other non-structural protein genes were identifiable among the six P. larvae phages. Jimmer1, Jimmer2, and Davies formed stable lysogens resistant to superinfection by genetically similar phages. The correlation between tape measure protein gene length and phage tail length allowed identification of co-isolated phages Emery and Abouo in electron micrographs. A Phamerator database was assembled with the P. larvae phage genomes and 107 genomes of Firmicute-infecting phages, including 71 Bacillus phages. Phamerator identified conserved domains in 1,501 of 6,181 phamilies (only 24.3%) encoded by genes in the database and revealed that P. larvae phage genomes shared at least one phamily with 72 of the 107 other phages. The phamily relationship of large terminase proteins was used to indicate putative DNA packaging strategies. Analyses from CoreGenes, Phamerator, and electron micrograph measurements indicated Jimmer1, Jimmer2, Abouo and Davies were related to phages phiC2, EJ-1, KC5a, and AQ113, which are small-genome myoviruses that infect Streptococcus, Lactobacillus, and Clostridium, respectively.

Conclusions

This paper represents the first comparison of phage genomes in the Paenibacillus genus and the first organization of P. larvae phages based on sequence and structure. This analysis provides an important contribution to the field of bacteriophage genomics by serving as a foundation on which to build an understanding of the natural predators of P. larvae.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-745) contains supplementary material, which is available to authorized users.     

Methyltransferases acquired by lactococcal 936-type phage provide protection against restriction endonuclease activity

Background

So-called 936-type phages are among the most frequently isolated phages in dairy facilities utilising Lactococcus lactis starter cultures. Despite extensive efforts to control phage proliferation and decades of research, these phages continue to negatively impact cheese production in terms of the final product quality and consequently, monetary return.

Results

Whole genome sequencing and in silico analysis of three 936-type phage genomes identified several putative (orphan) methyltransferase (MTase)-encoding genes located within the packaging and replication regions of the genome. Utilising SMRT sequencing, methylome analysis was performed on all three phages, allowing the identification of adenine modifications consistent with N-6 methyladenine sequence methylation, which in some cases could be attributed to these phage-encoded MTases. Heterologous gene expression revealed that M.Phi145I/M.Phi93I and M.Phi93DAM, encoded by genes located within the packaging module, provide protection against the restriction enzymes HphI and DpnII, respectively, representing the first functional MTases identified in members of 936-type phages.

Conclusions

SMRT sequencing technology enabled the identification of the target motifs of MTases encoded by the genomes of three lytic 936-type phages and these MTases represent the first functional MTases identified in this species of phage. The presence of these MTase-encoding genes on 936-type phage genomes is assumed to represent an adaptive response to circumvent host encoded restriction-modification systems thereby increasing the fitness of the phages in a dynamic dairy environment.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-831) contains supplementary material, which is available to authorized users.     

Whole-genome synthesis and characterization of viable S13-like bacteriophages

Background

Unprecedented progresses in high-throughput DNA sequencing and de novo gene synthesis technologies have allowed us to create living organisms in the absence of natural template.

Methodology/Principal Findings

The sequence of wild-type S13 phage genome was downloaded from GenBank. Two synonymous mutations were introduced into wt-S13 genome to generate m1-S13 genome. Another mutant, m2-S13 genome, was obtained by engineering two nonsynonymous mutations in the capsid protein coding region of wt-S13 genome. A chimeric phage genome was designed by replacing the F capsid protein open reading frame (ORF) from phage S13 with the F capsid protein ORF from phage G4. The whole genomes of all four phages were assembled from a series of chemically synthesized short overlapping oligonucleotides. The linear synthesized genomes were circularized and electroporated into E.coli C, the standard laboratory host of S13 phage. All four phages were recovered and plaques were visualized. The results of sequencing showed the accuracy of these synthetic genomes. The synthetic phages were capable of lysing their bacterial host and tolerating general environmental conditions. While no phenotypic differences among the variant strains were observed when grown in LB medium with CaCl₂, the S13/G4 chimera was found to be much more sensitive to the absence of calcium and to have a lower adsorption rate under calcium free condition.

Conclusions/Significance

The bacteriophage S13 and its variants can be chemically synthesized. The major capsid gene of phage G4 is functional in the phage S13 life cycle. These results support an evolutional hypothesis which has been proposed that a homologous recombination event involving gene F of quite divergent ancestral lineages should be included in the history of the microvirid family.     

Comparative Genome Analysis of Listeria Bacteriophages Reveals Extensive Mosaicism,Programmed Translational Frameshifting,and a Novel Prophage Insertion Site

The genomes of six Listeria bacteriophages were sequenced and analyzed. Phages A006, A500, B025, P35, and P40 are members of the Siphoviridae and contain double-stranded DNA genomes of between 35.6 kb and 42.7 kb. Phage B054 is a unique myovirus and features a 48.2-kb genome. Phage B025 features 3′ overlapping single-stranded genome ends, whereas the other viruses contain collections of terminally redundant, circularly permuted DNA molecules. Phages P35 and P40 have a broad host range and lack lysogeny functions, correlating with their virulent lifestyle. Phages A500, A006, and B025 integrate into bacterial tRNA genes, whereas B054 targets the 3′ end of translation elongation factor gene tsf. This is the first reported case of phage integration into such an evolutionarily conserved genetic element. Peptide fingerprinting of viral proteins revealed that both A118 and A500 utilize +1 and −1 programmed translational frameshifting for generating major capsid and tail shaft proteins with C termini of different lengths. In both cases, the unusual +1 frameshift at the 3′ ends of the tsh coding sequences is induced by overlapping proline codons and cis-acting shifty stops. Although Listeria phage genomes feature a conserved organization, they also show extensive mosaicism within the genome building blocks. Of particular interest is B025, which harbors a collection of modules and sequences with relatedness not only to other Listeria phages but also to viruses infecting other members of the Firmicutes. In conclusion, our results yield insights into the composition and diversity of Listeria phages and provide new information on their function, genome adaptation, and evolution.The opportunistic pathogen Listeria monocytogenes is ubiquitous in nature and can become endemic in food processing environments, causing contamination of dairy products, meats, vegetables, and processed ready-to-eat food (14). L. monocytogenes is the causative agent of epidemic and sporadic listeriosis. The risk of infection is markedly increased among immunocompromised patients, newborns, pregnant women, and the elderly and is associated with a mortality rate of about 20 to 30% (37).Although all strains of L. monocytogenes are considered potentially pathogenic, epidemiological evidence has shown that certain serovars are more frequently associated with both sporadic cases and larger food-borne outbreaks. However, genetic variation within the virulence genes of wild-type strains appears to be limited and could not be directly linked to differences in pathogenicity (30) or environmental distribution.It is becoming increasingly clear that bacteriophages have an important role in bacterial biology, diversity, and evolution, as indicated by the advances in genome sequencing which revealed a high incidence of phage-related sequences in bacterial genomes. Many phages have been described for the genus Listeria, and lysogeny appears to be widespread (28). Availability of the genome sequences of different L. monocytogenes and L. innocua strains also revealed the presence of several different putative prophages in the genomes of these bacteria, constituting up to 7% of the coding capacities of the genomes (15, 32). Although the genomes of L. monocytogenes were found to be essentially syntenic, a significant portion of sequence variability is apparently based upon phage insertions and subsequent rearrangements. Investigations of prophage contributions to population dynamics in Salmonella suggest that prophages can improve the competitive fitness of the lysogenized host strains (5). This type of selective pressure also results in diversification and generation of new strains by lysogenic conversion. In the case of Listeria, however, the potential influence of prophages on their host strains, such as phenotypic variation or provision of selective benefits, has not been investigated. To gain more insight in bacteriophage-host interactions and the molecular biology and characteristics of Listeria phages, more information on the structure, information content, and variability of different Listeria phage genomes is required.Although a number of Listeria phages have been isolated and described (25, 27, 42), only limited information was available at the sequence level for phages PSA, A118, A511, and P100 (9, 19, 26, 41). As the result of a comprehensive study to determine the diversity of this group of bacterial viruses, we here report the complete nucleotide sequences of a representative set of six different Listeria phages from the Siphoviridae (A006, A500, B025, P35, and P40) and Myoviridae (B054) families in the order Caudovirales. In addition to molecular and in silico analyses, we also determined the physical genome structures and attachment loci of the temperate phages, and we describe integration of the B054 prophage into a highly conserved elongation factor gene. Another interesting finding is the unusual decoding used by some of the phages, which use programmed frameshifting to generate C-terminally modified structural proteins required for assembly of the capsid and tail.     

The Linear Plasmid Prophage Vp58.5 of Vibrio parahaemolyticus Is Closely Related to the Integrating Phage VHML and Constitutes a New Incompatibility Group of Telomere Phages

Vibrio parahaemolyticus O3:K6 pandemic strains recovered in Chile frequently possess a 42-kb plasmid which is the prophage of a myovirus. We studied the prototype phage VP58.5 and show that it does not integrate into the host cell chromosome but replicates as a linear plasmid (Vp58.5) with covalently closed ends (telomeres). The Vp58.5 replicon coexists with other plasmid prophages (N15, PY54, and ΦKO2) in the same cell and thus belongs to a new incompatibility group of telomere phages. We determined the complete nucleotide sequence (42,612 nucleotides) of the VP58.5 phage DNA and compared it with that of the plasmid prophage. The two molecules share the same nucleotide sequence but are 35% circularly permuted to each other. In contrast to the hairpin ends of the plasmid, VP58.5 phage DNA contains 5′-protruding ends. The VP58.5 sequence is 92% identical to the sequence of phage VHML, which was reported to integrate into the host chromosome. However, the gene order and termini of the phage DNAs are different. The VHML genome exhibits the same gene order as does the Vp58.5 plasmid. VHML phage DNA has been reported to contain terminal inverted repeats. This repetitive sequence is similar to the telomere resolution site (telRL) of VP58.5 which, after processing by the phage protelomerase, forms the hairpin ends of the Vp58.5 prophage. It is discussed why these closely related phages may be so different in terms of their genome ends and their lifestyle.Most temperate bacteriophages integrate into the host chromosome during lysogeny. However, there are some phages (telomere phages) whose prophages are linear plasmids with covalently closed ends. Members of this group of phages are the siphoviruses N15, PY54, and ΦKO2 isolated from Escherichia coli, Yersinia enterocolitica, and Klebsiella oxytoca, respectively, and the recently described myoviruses ΦHAP-1 of Halomonas aquamarina and VP882 of Vibrio parahaemolyticus (6, 20, 23, 26, 37). Despite their different origins (enterobacteria versus marine bacterium) and morphologies, all known telomere phages share similar genome organizations and some protein similarities. The linear DNA of each phage is a circular permutation of the respective linear plasmid prophage. For the generation of the terminal hairpins of the linear plasmid, the protelomerase (Tel) is essential (8). This enzyme has cleaving/joining activity; its target is a large palindromic DNA sequence called the telomere resolution site (telRL) located upstream of tel on the phage genome. After cleaving telRL by staggered cuts, the resulting self-complementary single-stranded DNA overhangs fold back and are rejoined by the protelomerase (9). Besides tel, all telomere phages possess the gene repA, encoding a multifunctional replication protein. repA of N15 and PY54 was shown to harbor the prophage replication origin and to function as a circular minimal replicon (35, 42). Compatibility studies demonstrated that the N15 and ΦKO2 plasmids belong to the same incompatibility group, whereas the PY54 plasmid is able to coexist with these two prophages in doubly lysogenic E. coli and Y. enterocolitica hosts (19).There are some reports on the presence of tel and repA in prophages (VP882, VHML, and Vp58.5) of marine Vibrio strains (28, 41). V. parahaemolyticus phage VP882 is a close relative of the Halomonas phage ΦHAP-1 (26). VHML was isolated from a toxin-producing Vibrio harveyi strain, pathogenic for some crustaceans and fish (30). Similarly to ΦHAP-1 and VP882, VHML has a myovirus-like morphology. The phage contains genes for products similar to Tel and RepA, suggesting that its prophage is a linear plasmid with terminal hairpins. However, it was surmised that VHML integrates into the Vibrio chromosome (28, 29). Phage VP58.5 was isolated from a V. parahaemolyticus strain belonging to the serovar O3:K6 pandemic clonal complex (41). During the last several years, this clone has been associated with many seafood-borne diarrhea outbreaks in Southeast Asia and South America, particularly Chile (5, 12, 13, 15). Up to 33% of the Chilean isolates harbored a 42-kb plasmid which was shown to be the prophage of a myovirus inducible by mitomycin C. VP58.5 is the prototype of these phages.In this work we demonstrate that VP58.5 is closely related to the V. harveyi phage VHML but that its prophage is a linear plasmid with covalently closed ends. The Vp58.5 prophage belongs to a new incompatibility group of telomere phages.     

Evolution of Lactococcus lactis Phages within a Cheese Factory

We have sequenced the double-stranded DNA genomes of six lactococcal phages (SL4, CB13, CB14, CB19, CB20, and GR7) from the 936 group that were isolated over a 9-year period from whey samples obtained from a Canadian cheese factory. These six phages infected the same two industrial Lactococcus lactis strains out of 30 tested. The CB14 and GR7 genomes were found to be 100% identical even though they were isolated 14 months apart, indicating that a phage can survive in a cheese plant for more than a year. The other four genomes were related but notably different. The length of the genomes varied from 28,144 to 32,182 bp, and they coded for 51 to 55 open reading frames. All five genomes possessed a 3′ overhang cos site that was 11 nucleotides long. Several structural proteins were also identified by nano-high-performance liquid chromatography-tandem mass spectrometry, confirming bioinformatic analyses. Comparative analyses suggested that the most recently isolated phages (CB19 and CB20) were derived, in part, from older phage isolates (CB13 and CB14/GR7). The organization of the five distinct genomes was similar to the previously sequenced lactococcal phage genomes of the 936 group, and from these sequences, a core genome was determined for lactococcal phages of the 936 group.The manufacture of cheeses requires the inoculation of carefully selected bacterial cultures, known as starter cultures, at concentrations of at least 10⁷ live bacteria per ml of heat-treated milk. The purpose of this process is to control the fermentation and to obtain high-quality fermented products (29). Starter cultures are a combination of lactic acid bacteria (LAB), of which one of the most important species is Lactococcus lactis. L. lactis is a low-GC gram-positive bacterium used to metabolize lactose into lactic acid during the production of several cheese varieties. Because large amounts of lactococcal cells are cultivated each day in large-scale fermentation vats and because these cells are susceptible to bacteriophage infection, it is not surprising that most cheese factories have experienced problems with phage contamination (13). Even a single phage infecting a starter strain is enough to begin a chain reaction that can eventually inhibit bacterial growth and cause production delays, taste and texture variations, and even complete fermentation failures (1, 29).Phage infections are unpredictable in food fermentations. Their presence and persistence in a dairy factory can be explained in many ways. First, raw milk can introduce new phages into an industrial plant (25). Madera et al. (22) also reported that newly isolated lactococcal phages were more resistant to pasteurization. Whey, a liquid by-product of cheese manufacturing, is another reservoir that can spread phages in a factory environment (25). Airborne phage dissemination may also be important since concentrations of up to 10⁶ PFU/m³ have been observed close to a functional whey separation tank (32).For decades, the dairy industry has been working to curtail the propagation of virulent phages using a variety of practical strategies, including, among others, sanitation, optimized factory design, air filtration units, rotation of bacterial strains, and the use of phage resistance systems (13). Yet new virulent phages emerge on a regular basis. Indeed, large-scale industrial milk fermentation processes can be slowed down by virulent phages of the Caudovirales order. Members of three lactococcal phage groups, namely, 936, c2, and P335, are mostly found in dairy plants. The 936-like phages are by far the most predominant worldwide (3, 18, 22, 27).Phages of the 936 group have a double-stranded DNA genome and possess a long noncontractile tail connected to a capsid with icosahedral symmetry characteristic of the Siphoviridae family. Currently, six complete phage genomes of the lactococcal 936 group are available in public databases, including sk1 (6), bIL170 (10), jj50, 712, P008 (23), and bIBB29 (16). Their comparative analysis revealed a conserved gene organization despite being isolated from different countries. Most of the differences have been observed in the early gene module, where insertions, deletions, and point mutations likely occurred (16, 23). Moreover, it is assumed that these phages can also exchange DNA through recombination with other bacterial viruses present in the same ecosystem.Because new members of this lactococcal phage group are regularly isolated, a better understanding of their evolution is warranted to better control them. A cheese factory is a particular man-made niche where rapidly growing bacterial strains encounter ubiquitous phages. Such active environments provide ample opportunities for phage evolution, especially to dodge phage resistance mechanisms that may be present in host cells. Nonetheless, the evolutionary dynamics that shape the diversity of lactococcal phage populations are still not well understood.In this study, we analyzed the genome and structural proteome of six 936-group phages (SL4, CB13, CB14, CB19, CB20, and GR7) that infected the same L. lactis strains and were isolated over a 9-year period from a cheese factory.     

Comparative Genomic Analysis of Ten Streptococcus pneumoniae Temperate Bacteriophages

Streptococcus pneumoniae is an important human pathogen that often carries temperate bacteriophages. As part of a program to characterize the genetic makeup of prophages associated with clinical strains and to assess the potential roles that they play in the biology and pathogenesis in their host, we performed comparative genomic analysis of 10 temperate pneumococcal phages. All of the genomes are organized into five major gene clusters: lysogeny, replication, packaging, morphogenesis, and lysis clusters. All of the phage particles observed showed a Siphoviridae morphology. The only genes that are well conserved in all the genomes studied are those involved in the integration and the lysis of the host in addition to two genes, of unknown function, within the replication module. We observed that a high percentage of the open reading frames contained no similarities to any sequences catalogued in public databases; however, genes that were homologous to known phage virulence genes, including the pblB gene of Streptococcus mitis and the vapE gene of Dichelobacter nodosus, were also identified. Interestingly, bioinformatic tools showed the presence of a toxin-antitoxin system in the phage φSpn_6, and this represents the first time that an addition system in a pneumophage has been identified. Collectively, the temperate pneumophages contain a diverse set of genes with various levels of similarity among them.Streptococcus pneumoniae (the pneumococcus) is an important human pathogen and a major etiological agent of pneumonia, bacteremia, and meningitis in adults and of otitis media in children. The casualties due to the pneumococcus are estimated to be over 1.6 million deaths per year, and most of these deaths are of young children in developing countries (40). S. pneumoniae is also a human commensal that resides in the upper respiratory tract, and it is asymptomatically carried in the nasopharynx of up to 60% of the normal population (48).Bacteriophages of S. pneumoniae (pneumophages) were first identified in 1975 from samples isolated from throat swabs of healthy children by two independent groups (46, 65). Since then, pneumophages have been identified from different sources and a variety of locations (44). The abundance of temperate bacteriophages in S. pneumoniae has been reported in different studies in the past (6, 54). Up to 76% of clinical isolates have been showed to contain prophages (or prophage remnants) when studied with a DNA probe specific for the major autolysin gene, lytA, which hybridizes with many of the endolysin genes of temperate pneumococcal phages (54). Hybridization analyses have identified highly similar prophages among pneumococcal clinical isolates even of different capsular serotypes, a result which indicates the widespread distribution of these mobile genetic elements among virulent strains (26).Only three S. pneumoniae bacteriophage genomes have been characterized in detail, and their sequences have been determined. Dp-1 and Cp-1 are lytic bacteriophages, whereas MM1 is a temperate pneumophage (45, 50, 52). Genes coding for virulence factors such as toxins or secreted enzymes have been associated with the presence of prophages in both gram-negative (67) and gram-positive bacteria, such as Streptococcus pyogenes (7) and Staphylococcus aureus (23). Because a considerable number of toxin genes are located in prophages, phage dynamics are of apparent importance for bacterial pathogenesis. Unfortunately, the role of temperate bacteriophages in the virulence of S. pneumoniae remains mostly unknown.Recently, the availability of relatively inexpensive next-generation sequencing technologies has permitted the complete genomic analysis of dozens of genomes of pneumococcal clinical isolates. In this report, we present a comparative genomic analysis of 10 pneumophages identified in the genomes of newly sequenced S. pneumoniae strains. The proteome of these phages has been predicted and annotated by comparative sequence analyses by using the available databases at the National Center for Biotechnological Information website (http://www.ncbi.nlm.nih.gov/). This systematic characterization of pneumophage genomes provides for a substantial increase in our knowledge of the global proteome and the overall genetic diversity of this important human pathogen. The comparative analysis of multiple temperate bacteriophages from a single species offers a unique opportunity to study one of the mechanisms of lateral gene transfer that drive prokaryotic genetic diversity.     

Core and accessory genome architecture in a group of Pseudomonas aeruginosa Mu-like phages

Background

Bacteriophages that infect the opportunistic pathogen Pseudomonas aeruginosa have been classified into several groups. One of them, which includes temperate phage particles with icosahedral heads and long flexible tails, bears genomes whose architecture and replication mechanism, but not their nucleotide sequences, are like those of coliphage Mu. By comparing the genomic sequences of this group of P. aeruginosa phages one could draw conclusions about their ontogeny and evolution.

Results

Two newly isolated Mu-like phages of P. aeruginosa are described and their genomes sequenced and compared with those available in the public data banks. The genome sequences of the two phages are similar to each other and to those of a group of P. aeruginosa transposable phages. Comparing twelve of these genomes revealed a common genomic architecture in the group. Each phage genome had numerous genes with homologues in all the other genomes and a set of variable genes specific for each genome. The first group, which comprised most of the genes with assigned functions, was named “core genome”, and the second group, containing mostly short ORFs without assigned functions was called “accessory genome”. Like in other phage groups, variable genes are confined to specific regions in the genome.

Conclusion

Based on the known and inferred functions for some of the variable genes of the phages analyzed here, they appear to confer selective advantages for the phage survival under particular host conditions. We speculate that phages have developed a mechanism for horizontally acquiring genes to incorporate them at specific loci in the genome that help phage adaptation to the selective pressures imposed by the host.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1146) contains supplementary material, which is available to authorized users.     

A novel termini analysis theory using HTS data alone for the identification of Enterococcus phage EF4-like genome termini

Background

Enterococcus faecalis and Enterococcus faecium are typical enterococcal bacterial pathogens. Antibiotic resistance means that the identification of novel E. faecalis and E. faecium phages against antibiotic-resistant Enterococcus have an important impact on public health. In this study, the E. faecalis phage IME-EF4, E. faecium phage IME-EFm1, and both their hosts were antibiotic resistant. To characterize the genome termini of these two phages, a termini analysis theory was developed to provide a wealth of terminal sequence information directly, using only high-throughput sequencing (HTS) read frequency statistics.

Results

The complete genome sequences of phages IME-EF4 and IME-EFm1 were determined, and our termini analysis theory was used to determine the genome termini of these two phages. Results showed 9 bp 3′ protruding cohesive ends in both IME-EF4 and IME-EFm1 genomes by analyzing frequencies of HTS reads. For the positive strands of their genomes, the 9 nt 3′ protruding cohesive ends are 5′-TCATCACCG-3′ (IME-EF4) and 5′-GGGTCAGCG-3′ (IME-EFm1). Further experiments confirmed these results. These experiments included mega-primer polymerase chain reaction sequencing, terminal run-off sequencing, and adaptor ligation followed by run-off sequencing.

Conclusion

Using this termini analysis theory, the termini of two newly isolated antibiotic-resistant Enterococcus phages, IME-EF4 and IME-EFm1, were identified as the byproduct of HTS. Molecular biology experiments confirmed the identification. Because it does not require time-consuming wet lab termini analysis experiments, the termini analysis theory is a fast and easy means of identifying phage DNA genome termini using HTS read frequency statistics alone. It may aid understanding of phage DNA packaging.     

Genome Organization and Characterization of the Virulent Lactococcal Phage 1358 and Its Similarities to Listeria Phages

Virulent phage 1358 is the reference member of a rare group of phages infecting Lactococcus lactis. Electron microscopy revealed a typical icosahedral capsid connected to one of the smallest noncontractile tails found in a lactococcal phage of the Siphoviridae family. Microbiological characterization identified a burst size of 72 virions released per infected host cell and a latent period of 90 min. The host range of phage 1358 was limited to 3 out of the 60 lactococcal strains tested. Moreover, this phage was insensitive to four Abi systems (AbiK, AbiQ, AbiT, and AbiV). The genome of phage 1358 consisted of a linear, double-stranded, 36,892-bp DNA molecule containing 43 open reading frames (ORFs). At least 14 ORFs coded for structural proteins, as identified by SDS-PAGE coupled to liquid chromatography-tandem mass spectrometry (LC-MS/MS) analyses. The genomic organization was similar to those of other siphophages. All genes were on the same coding strand and in the same orientation. This lactococcal phage was unique, however, in its 51.4% GC content, much higher than those of other phages infecting this low-GC Gram-positive host. A bias for GC-rich codons was also observed. Comparative analyses showed that several phage 1358 structural proteins shared similarity with two Listeria monocytogenes phages, P35 and P40. The possible origin and evolution of lactococcal phage 1358 is discussed.The first sequenced genome of a phage infecting Lactococcus lactis (bIL67) was reported in 1994 (57). Its genomic characterization was performed with the prospect of a better understanding of lactococcal phage biology. L. lactis is a Gram-positive bacterium added to milk to produce an array of fermented dairy products. In this human-made environment, substantial amounts of lactococcal cells are cultivated on a daily basis in large fermentation vats, and these added cells randomly encounter virulent phages present in heat-treated but nonsterile milk. Moreover, it is widely acknowledged that the increased use of the same bacterial strains within existing dairy facilities inevitably leads to milk fermentation failures due to the multiplication of virulent phages. This biotechnological problem reduces yields and lowers the quality of fermented products (51).Over 700 lactococcal phage isolates have been reported in the literature (3). To date, more than 25 complete genome sequences of lactococcal phages are publicly available in the NCBI database, and the sequencing of others is under way. These numbers indicate that Lactococcus phages are among the most studied of the bacterial viruses. All lactococcal phages belong to the order Caudovirales and are included within two families according to their tail morphology: the Siphoviridae (long noncontractile tail [most lactococcal phages]) and the Podoviridae (short noncontractile tail [few lactococcal phages]) (14). Currently, phages infecting L. lactis strains have been divided into 10 genetically distinct groups (14). The complete genomic sequence is available for at least one representative of 8 of the groups.Early sequencing efforts concentrated on the genomes of lactococcal phages belonging to the 936, c2, and P335 groups (Siphoviridae), because members of these groups were regularly isolated in dairy plants (8, 36, 50). PCR-based methods were also devised to rapidly classify these phages (41). These Siphoviridae phages pose a significant risk to the dairy industry, and their characterization is important for developing adapted antiphage strategies to limit their propagation and evolution.In recent years, representatives of the less recognized lactococcal phage groups have been characterized, including phages Q54 (22), KSY1 (13), 1706 (23), asccφ28 of the P034 group (39), and P087 (63). Their molecular characterizations were aimed at understanding why some phage groups (936, c2, and P335) predominate while the others have remained marginal, at best. However, it was recently reported that P034-like phages may be emerging in certain regions (52). Genomic and microbiological analyses indicated that members of these rare phage groups were likely the result of recombination between different lactococcal phages and phages infecting other Gram-positive bacteria, and they may not be fit to multiply rapidly in milk. For example, lactococcal phage 1706 shares similarities with Ruminococcus and Clostridium prophages (23). Similarly, L. lactis phage P087 structural proteins share identity with gene products found in a prophage in the Enterococcus faecalis genome (63). It was also shown previously that lactococcal phage asccφ28 was related to Streptococcus pneumoniae phage Cp-1 and Bacillus subtilis φ29-like phages (39). It was suggested that phages 1706, asccφ28, and P087 acquired a receptor-binding protein complex from another lactococcal phage that enabled them to infect a L. lactis host.Here, we report the complete genome sequence and analysis of phage 1358, a virulent representative of the 9th lactococcal phage group.