Diversity and distribution of single-stranded DNA phages in the North Atlantic Ocean

Knowledge of marine phages is highly biased toward double-stranded DNA (dsDNA) phages; however, recent metagenomic surveys have also identified single-stranded DNA (ssDNA) phages in the oceans. Here, we describe two complete ssDNA phage genomes that were reconstructed from a viral metagenome from 80 m depth at the Bermuda Atlantic Time-series Study (BATS) site in the northwestern Sargasso Sea and examine their spatial and temporal distributions. Both genomes (SARssφ1 and SARssφ2) exhibited similarity to known phages of the Microviridae family in terms of size, GC content, genome organization and protein sequence. PCR amplification of the replication initiation protein (Rep) gene revealed narrow and distinct depth distributions for the newly described ssDNA phages within the upper 200 m of the water column at the BATS site. Comparison of Rep gene sequences obtained from the BATS site over time revealed changes in the diversity of ssDNA phages over monthly time scales, although some nearly identical sequences were recovered from samples collected 4 years apart. Examination of ssDNA phage diversity along transects through the North Atlantic Ocean revealed a positive correlation between genetic distance and geographic distance between sampling sites. Together, the data suggest fundamental differences between the distribution of these ssDNA phages and the distribution of known marine dsDNA phages, possibly because of differences in host range, host distribution, virion stability, or viral evolution mechanisms and rates. Future work needs to elucidate the host ranges for oceanic ssDNA phages and determine their ecological roles in the marine ecosystem.     

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.     

Narrow-Host-Range Bacteriophages That Infect Rhizobium etli Associate with Distinct Genomic Types

In this work, we isolated and characterized 14 bacteriophages that infect Rhizobium etli. They were obtained from rhizosphere soil of bean plants from agricultural lands in Mexico using an enrichment method. The host range of these phages was narrow but variable within a collection of 48 R. etli strains. We obtained the complete genome sequence of nine phages. Four phages were resistant to several restriction enzymes and in vivo cloning, probably due to nucleotide modifications. The genome size of the sequenced phages varied from 43 kb to 115 kb, with a median size of ∼45 to 50 kb. A large proportion of open reading frames of these phage genomes (65 to 70%) consisted of hypothetical and orphan genes. The remainder encoded proteins needed for phage morphogenesis and DNA synthesis and processing, among other functions, and a minor percentage represented genes of bacterial origin. We classified these phages into four genomic types on the basis of their genomic similarity, gene content, and host range. Since there are no reports of similar sequences, we propose that these bacteriophages correspond to novel species.     

Comparative whole genome analysis of six diagnostic brucellaphages

Whole genome sequencing of six diagnostic brucellaphages, Tbilisi (Tb), Firenze (Fz), Weybridge (Wb), S708, Berkeley (Bk) and R/C, was followed with genomic comparisons including recently described genomes of the Tb phage from Mexico (Tb_M) and Pr phage to elucidate genomic diversity and candidate host range determinants. Comparative whole genome analysis revealed high sequence homogeneity among these brucellaphage genomes and resolved three genetic groups consistent with defined host range phenotypes. Group I was composed of Tb and Fz phages that are predominantly lytic for Brucella abortus and Brucella neotomae; Group II included Bk, R/C, and Pr phages that are lytic mainly for B. abortus, Brucella melitensis and Brucella suis; Group III was composed of Wb and S708 phages that are lytic for B. suis, B. abortus and B. neotomae. We found that the putative phage collar protein is a variable locus with features that may be contributing to the host specificities exhibited by different brucellaphage groups. The presence of several candidate host range determinants is illustrated herein for future dissection of the differential host specificity observed among these phages.     

Identification of a New P335 Subgroup through Molecular Analysis of Lactococcal Phages Q33 and BM13

Lactococcal dairy starter strains are under constant threat from phages in dairy fermentation facilities, especially by members of the so-called 936, P335, and c2 species. Among these three phage groups, members of the P335 species are the most genetically diverse. Here, we present the complete genome sequences of two P335-type phages, Q33 and BM13, isolated in North America and representing a novel lineage within this phage group. The Q33 and BM13 genomes exhibit homology, not only to P335-type, but also to elements of the 936-type phage sequences. The two phage genomes also have close relatedness to phages infecting Enterococcus and Clostridium, a heretofore unknown feature among lactococcal P335 phages. The Q33 and BM13 genomes are organized in functionally related clusters with genes encoding functions such as DNA replication and packaging, morphogenesis, and host cell lysis. Electron micrographic analysis of the two phages highlights the presence of a baseplate more reminiscent of the baseplate of 936 phages than that of the majority of members of the P335 group, with the exception of r1t and LC3.     

CRISPR-Cas: an efficient tool for genome engineering of virulent bacteriophages

Bacteriophages are now widely recognized as major players in a wide variety of ecosystems. Novel genes are often identified in newly isolated phages as well as in environmental metavirome studies. Most of these novel viral genes have unknown functions but appear to be coding for small, non-structural proteins. To understand their biological role, very efficient genetic tools are required to modify them, especially in the genome of virulent phages. We first show that specific point mutations and large deletions can be engineered in the genome of the virulent phage 2972 using the Streptococcus thermophilus CRISPR-Cas Type II-A system as a selective pressure to increase recombination efficiencies. Of significance, all the plaques tested contained recombinant phages with the desired mutation. Furthermore, we show that the CRISPR-Cas engineering system can be used to efficiently introduce a functional methyltransferase gene into a virulent phage genome. Finally, synthetic CRISPR bacteriophage insensitive mutants were constructed by cloning a spacer-repeat unit in a low-copy vector illustrating the possibility to target multiple regions of the phage genome. Taken together, this data shows that the CRISPR-Cas system is an efficient and adaptable tool for editing the otherwise intractable genomes of virulent phages and to better understand phage-host interactions.     

驯化噬菌体提高噬菌体对碳青霉烯类耐药肺炎克雷伯菌的杀菌能力

【目的】本研究旨在通过驯化提高噬菌体的裂解能力并降低其宿主菌耐受性产生的速度,从而提高对重要病原菌-碳青霉烯类耐药肺炎克雷伯菌(carbapenem-resistant Klebsiella pneumoniae, CRKp)的杀菌效果。【方法】以临床CRKp菌株Kp2092为宿主菌,利用双层琼脂平板法从污水中分离噬菌体并分析其裂解谱;对其中的广谱强裂解性噬菌体通过透射电镜观察其形态特征并进行全基因组测序;通过噬菌体-宿主连续培养进行噬菌体驯化,并比较驯化前后噬菌体生物学特性的差异。【结果】分离得到的9株肺炎克雷伯菌噬菌体中,噬菌体P55anc裂解能力强且裂解谱广,透射电镜观察发现其为短尾噬菌体。P55anc基因组全长40 301 bp,包含51个编码序列,其中27个具有已知功能,主要涉及核酸代谢、噬菌体结构蛋白、DNA包装和细胞裂解等。噬菌体P55anc经9 d的驯化后,得到3株驯化噬菌体。驯化后噬菌体杀菌能力增强,主要表现为细菌生长曲线显著下降、噬菌体暴发量增多、裂解谱扩大,且宿主菌对其产生抗性的概率显著降低。与此同时,驯化后的噬菌体在热处理、紫外暴露以及血清等环境下保持较好的稳定性。【结论】利用噬菌体-宿主连续培养的方法可对噬菌体进行驯化和筛选,驯化后的噬菌体杀菌效果更强,且在不同压力处理下的稳定性良好,而细菌产生噬菌体抗性的概率也降低。     

Genomic Diversity of Phages Infecting Probiotic Strains of Lactobacillus paracasei

Strains of the Lactobacillus casei group have been extensively studied because some are used as probiotics in foods. Conversely, their phages have received much less attention. We analyzed the complete genome sequences of five L. paracasei temperate phages: C_L1, C_L2, iLp84, iLp1308, and iA2. Only phage iA2 could not replicate in an indicator strain. The genome lengths ranged from 34,155 bp (iA2) to 39,474 bp (C_L1). Phages iA2 and iLp1308 (34,176 bp) possess the smallest genomes reported, thus far, for phages of the L. casei group. The GC contents of the five phage genomes ranged from 44.8 to 45.6%. As observed with many other phages, their genomes were organized as follows: genes coding for DNA packaging, morphogenesis, lysis, lysogeny, and replication. Phages C_L1, C_L2, and iLp1308 are highly related to each other. Phage iLp84 was also related to these three phages, but the similarities were limited to gene products involved in DNA packaging and structural proteins. Genomic fragments of phages C_L1, C_L2, iLp1308, and iLp84 were found in several genomes of L. casei strains. Prophage iA2 is unrelated to these four phages, but almost all of its genome was found in at least four L. casei strains. Overall, these phages are distinct from previously characterized Lactobacillus phages. Our results highlight the diversity of L. casei phages and indicate frequent DNA exchanges between phages and their hosts.     

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.     

Genetic Variation of Lactobacillus delbrueckii subsp. lactis Bacteriophages Isolated from Cheese Processing Plants in Finland

The genomes of four Lactobacillus delbrueckii subsp. lactis bacteriophages were characterized by restriction endonuclease mapping, Southern hybridization, and heteroduplex analysis. The phages were isolated from different cheese processing plants in Finland between 1950 and 1972. All four phages had a small isometric head and a long noncontractile tail. Two different types of genome (double-stranded DNA) organization existed among the different phages, the pac type and the cos type, corresponding to alternative types of phage DNA packaging. Three phages belonged to the pac type, and a fourth was a cos-type phage. The pac-type phages were genetically closely related. In the genomes of the pac-type phages, three putative insertion/deletions (0.7 to 0.8 kb, 1.0 kb, and 1.5 kb) and one other region (0.9 kb) containing clustered base substitutions were discovered and localized. At the phenotype level, three main differences were observed among the pac-type phages. These concerned two minor structural proteins and the efficiency of phage DNA packaging. The genomes of the pac-type phages showed only weak homology with that of the cos-type phage. Phage-related DNA, probably a defective prophage, was located in the chromosome of the host strain sensitive to the cos-type phage. This DNA exhibited homology under stringent conditions to the pac-type phages.     

Genome comparison of Pseudomonas aeruginosa large phages

Pseudomonas aeruginosa phage EL is a dsDNA phage related to the giant phiKZ-like Myoviridae. The EL genome sequence comprises 211,215 bp and has 201 predicted open reading frames (ORFs). The EL genome does not share DNA sequence homology with other viruses and micro-organisms sequenced to date. However, one-third of the predicted EL gene products (gps) shares similarity (Blast alignments of 17-55% amino acid identity) with phiKZ proteins. Comparative EL and phiKZ genomics reveals that these giant phages are an example of substantially diverged genetic mosaics. Based on the position of similar EL and phiKZ predicted gene products, five genome regions can be delineated in EL, four of which are relatively conserved between EL and phiKZ. Region IV, a 17.7 kb genome region with 28 predicted ORFs, is unique to EL. Fourteen EL ORFs have been assigned a putative function based on protein similarity. Assigned proteins are involved in DNA replication and nucleotide metabolism (NAD+-dependent DNA ligase, ribonuclease HI, helicase, thymidylate kinase), host lysis and particle structure. EL-gp146 is the first chaperonin GroEL sequence identified in a viral genome. Besides a putative transposase, EL harbours predicted mobile endonucleases related to H-N-H and LAGLIDADG homing endonucleases associated with group I intron and intein intervening sequences.     

Genome sequence of a novel deep-sea vent epsilonproteobacterial phage provides new insight into the co-evolution of Epsilonproteobacteria and their phages

Epsilonproteobacteria are among the predominant primary producers in deep-sea hydrothermal vent ecosystems. However, phages infecting deep-sea vent Epsilonproteobacteria have never been isolated and characterized. Here, we successfully isolated a novel temperate phage, NrS-1, that infected a deep-sea vent chemolithoautotrophic isolate of Epsilonproteobacteria, Nitratiruptor sp. SB155-2, and its entire genome sequence was obtained and analyzed. The NrS-1 genome is linear, circularly permuted, and terminally redundant. The NrS-1 genome is 37,159 bp in length and contains 51 coding sequences. Five major structural proteins including major capsid protein and tape measure protein were identified by SDS-PAGE and mass spectrometry analysis. NrS-1 belongs to the family Siphoviridae, but its sequence and genomic organization are distinct from those of any other previously known Siphoviridae phages. Homologues of genes encoded in the NrS-1 genome were widely distributed among the genomes of diverse Epsilonproteobacteria. The distribution patterns had little relation to the evolutionary traits and ecological and physiological differentiation of the host epsilonproteobacterial species. The widespread occurrence of phage genes in diverse Epsilonproteobacteria supports early co-evolution between temperate phages and Epsilonproteobacteria prior to the divergence of their habitats and physiological adaptation.     

Classification and quantification of bacteriophage taxa in human gut metagenomes

Bacteriophages have key roles in microbial communities, to a large extent shaping the taxonomic and functional composition of the microbiome, but data on the connections between phage diversity and the composition of communities are scarce. Using taxon-specific marker genes, we identified and monitored 20 viral taxa in 252 human gut metagenomic samples, mostly at the level of genera. On average, five phage taxa were identified in each sample, with up to three of these being highly abundant. The abundances of most phage taxa vary by up to four orders of magnitude between the samples, and several taxa that are highly abundant in some samples are absent in others. Significant correlations exist between the abundances of some phage taxa and human host metadata: for example, ‘Group 936 lactococcal phages'' are more prevalent and abundant in Danish samples than in samples from Spain or the United States of America. Quantification of phages that exist as integrated prophages revealed that the abundance profiles of prophages are highly individual-specific and remain unique to an individual over a 1-year time period, and prediction of prophage lysis across the samples identified hundreds of prophages that are apparently active in the gut and vary across the samples, in terms of presence and lytic state. Finally, a prophage–host network of the human gut was established and includes numerous novel host–phage associations.     

A Virulent Phage Infecting Lactococcus garvieae,with Homology to Lactococcus lactis Phages

A new virulent phage belonging to the Siphoviridae family and able to infect Lactococcus garvieae strains was isolated from compost soil. Phage GE1 has a prolate capsid (56 by 38 nm) and a long noncontractile tail (123 nm). It had a burst size of 139 and a latent period of 31 min. Its host range was limited to only two L. garvieae strains out of 73 tested. Phage GE1 has a double-stranded DNA genome of 24,847 bp containing 48 predicted open reading frames (ORFs). Putative functions could be assigned to only 14 ORFs, and significant matches in public databases were found for only 17 ORFs, indicating that GE1 is a novel phage and its genome contains several new viral genes and encodes several new viral proteins. Of these 17 ORFs, 16 were homologous to deduced proteins of virulent phages infecting the dairy bacterium Lactococcus lactis, including previously characterized prolate-headed phages. Comparative genome analysis confirmed the relatedness of L. garvieae phage GE1 to L. lactis phages c2 (22,172 bp) and Q54 (26,537 bp), although its genome organization was closer to that of phage c2. Phage GE1 did not infect any of the 58 L. lactis strains tested. This study suggests that phages infecting different lactococcal species may have a common ancestor.     

Origins of the Xylella fastidiosa Prophage-Like Regions and Their Impact in Genome Differentiation

Xylella fastidiosa is a Gram negative plant pathogen causing many economically important diseases, and analyses of completely sequenced X. fastidiosa genome strains allowed the identification of many prophage-like elements and possibly phage remnants, accounting for up to 15% of the genome composition. To better evaluate the recent evolution of the X. fastidiosa chromosome backbone among distinct pathovars, the number and location of prophage-like regions on two finished genomes (9a5c and Temecula1), and in two candidate molecules (Ann1 and Dixon) were assessed. Based on comparative best bidirectional hit analyses, the majority (51%) of the predicted genes in the X. fastidiosa prophage-like regions are related to structural phage genes belonging to the Siphoviridae family. Electron micrograph reveals the existence of putative viral particles with similar morphology to lambda phages in the bacterial cell in planta. Moreover, analysis of microarray data indicates that 9a5c strain cultivated under stress conditions presents enhanced expression of phage anti-repressor genes, suggesting switches from lysogenic to lytic cycle of phages under stress-induced situations. Furthermore, virulence-associated proteins and toxins are found within these prophage-like elements, thus suggesting an important role in host adaptation. Finally, clustering analyses of phage integrase genes based on multiple alignment patterns reveal they group in five lineages, all possessing a tyrosine recombinase catalytic domain, and phylogenetically close to other integrases found in phages that are genetic mosaics and able to perform generalized and specialized transduction. Integration sites and tRNA association is also evidenced. In summary, we present comparative and experimental evidence supporting the association and contribution of phage activity on the differentiation of Xylella genomes.     

Novel Virulent and Broad-Host-Range Erwinia amylovora Bacteriophages Reveal a High Degree of Mosaicism and a Relationship to Enterobacteriaceae Phages

A diverse set of 24 novel phages infecting the fire blight pathogen Erwinia amylovora was isolated from fruit production environments in Switzerland. Based on initial screening, four phages (L1, M7, S6, and Y2) with broad host ranges were selected for detailed characterization and genome sequencing. Phage L1 is a member of the Podoviridae, with a 39.3-kbp genome featuring invariable genome ends with direct terminal repeats. Phage S6, another podovirus, was also found to possess direct terminal repeats but has a larger genome (74.7 kbp), and the virus particle exhibits a complex tail fiber structure. Phages M7 and Y2 both belong to the Myoviridae family and feature long, contractile tails and genomes of 84.7 kbp (M7) and 56.6 kbp (Y2), respectively, with direct terminal repeats. The architecture of all four phage genomes is typical for tailed phages, i.e., organized into function-specific gene clusters. All four phages completely lack genes or functions associated with lysogeny control, which correlates well with their broad host ranges and indicates strictly lytic (virulent) lifestyles without the possibility for host lysogenization. Comparative genomics revealed that M7 is similar to E. amylovora virus ΦEa21-4, whereas L1, S6, and Y2 are unrelated to any other E. amylovora phage. Instead, they feature similarities to enterobacterial viruses T7, N4, and ΦEcoM-GJ1. In a series of laboratory experiments, we provide proof of concept that specific two-phage cocktails offer the potential for biocontrol of the pathogen.     

Plasticity of the gene functions for DNA replication in the T4-like phages

We have completely sequenced and annotated the genomes of several relatives of the bacteriophage T4, including three coliphages (RB43, RB49 and RB69), three Aeromonas salmonicida phages (44RR2.8t, 25 and 31) and one Aeromonas hydrophila phage (Aeh1). In addition, we have partially sequenced and annotated the T4-like genomes of coliphage RB16 (a close relative of RB43), A. salmonicida phage 65, Acinetobacter johnsonii phage 133 and Vibrio natriegens phage nt-1. Each of these phage genomes exhibited a unique sequence that distinguished it from its relatives, although there were examples of genomes that are very similar to each other. As a group the phages compared here diverge from one another by several criteria, including (a) host range, (b) genome size in the range between approximately 160 kb and approximately 250 kb, (c) content and genetic organization of their T4-like genes for DNA metabolism, (d) mutational drift of the predicted T4-like gene products and their regulatory sites and (e) content of open-reading frames that have no counterparts in T4 or other known organisms (novel ORFs). We have observed a number of DNA rearrangements of the T4 genome type, some exhibiting proximity to putative homing endonuclease genes. Also, we cite and discuss examples of sequence divergence in the predicted sites for protein-protein and protein-nucleic acid interactions of homologues of the T4 DNA replication proteins, with emphasis on the diversity in sequence, molecular form and regulation of the phage-encoded DNA polymerase, gp43. Five of the sequenced phage genomes are predicted to encode split forms of this polymerase. Our studies suggest that the modular construction and plasticity of the T4 genome type and several of its replication proteins may offer resilience to mutation, including DNA rearrangements, and facilitate the adaptation of T4-like phages to different bacterial hosts in nature.     

Evolutionary Relationships among Actinophages and a Putative Adaptation for Growth in Streptomyces spp.

The genome sequences of eight Streptomyces phages are presented, four of which were isolated for this study. Phages R4, TG1, ϕHau3, and SV1 were isolated previously and have been exploited as tools for understanding and genetically manipulating Streptomyces spp. We also extracted five apparently intact prophages from recent Streptomyces spp. genome projects and, together with six phage genomes in the database, we analyzed all 19 Streptomyces phage genomes with a view to understanding their relationships to each other and to other actinophages, particularly the mycobacteriophages. Fifteen of the Streptomyces phages group into four clusters of related genomes. Although the R4-like phages do not share nucleotide sequence similarity with other phages, they clearly have common ancestry with cluster A mycobacteriophages, sharing many protein homologues, common gene syntenies, and similar repressor-stoperator regulatory systems. The R4-like phage ϕHau3 and the prophage StrepC.1 (from Streptomyces sp. strain C) appear to have hijacked a unique adaptation of the streptomycetes, i.e., use of the rare UUA codon, to control translation of the essential phage protein, the terminase. The Streptomyces venezuelae generalized transducing phage SV1 was used to predict the presence of other generalized transducing phages for different Streptomyces species.