Comparison of the genome for phylogenetically related bacteriophages ϕKZ and EL of <Emphasis Type="Italic">Pseudomonas aeruginosa</Emphasis>: Evolutionary aspects and minimal genome size

Bacteriophages of the family Myoviridae represent one of the most widespread domains of the biosphere substantially affecting the ecological balance of microorganisms. Interestingly, sequence analysis of genomic DNAs of large bacteriophages revealed many genes coding for proteins with unknown functions. A new approach is proposed to improve the functional identification of genes. This approach is based on comparing the genome sequence for phylogenetically and morphologically related phages showing no considerable homology at the level of genomic DNA. It is assumed that gene functions essential for the development of phages of a given family are conserved and that the corresponding genes code for similar orthologous proteins even when lacking sequence homology. The genome was sequenced and compared for two Pseudomonas aeruginosa giant bacteriophages, KZ and EL, which belong to a group of KZ-related phages. A substantial difference in genome organization was observed, suggesting specific features of phage evolution. In addition, the problem of the minimal genome of the superfamily is discussed on the basis of the difference in size and structure between the KZ and EL genomes.__________Translated from Genetika, Vol. 41, No. 4, 2005, pp. 455–465.Original Russian Text Copyright © 2005 by Krylov, Pleteneva, Lavigne, Hertveldt, Volckaert, Sernova, Georgopoulos, Korchevskii, Kurochkina, Mesyanzhinov.     

Comparison of DNA sizes in a group of giant Pseudomonas aeruginosa phages by the PFGE method

Genome sizes of Pseudomonas aeruginosa phages phiKZ and EL earlier determined by sequence analysis were shown to correspond to sizes of their DNAs assessed by pulse-electrophoresis (PFGE). Putative "redundant" genes in phiKZ phage genome are supposed to control functions promoting vigorous growth of the phage belonging to this species, compared to phages of EL species.     

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.     

Comparison of new giant bacteriophages OBP and Lu11 of soil pseudomonads with bacteriophages of phiKZ-supergroup of Pseudomonas aeruginosa

Study of two recently isolated giant bacteriophages Lu11 and OBP that are active on Pseudomonas putida var. Manila and Pseudomonas fluorescens, respectively, demonstrated their similarity in morphology, genome size, and size of phage particles, with giant bacteriophages of Pseudomonas aeruginosa assigned to the supergroup of phiKZ-like phages of the family Myoviridae designated in this manner according to the best studied phage phiKZ that belongs to the species of this group widely distributed in nature. Comparison of major polypeptide sizes of mature particles suggests the similarity of certain proteins in the phages examined. In OBP particles visualized with an electron microscope, an "inner body" was detected, which points to the specific DNA package intrinsic to phages of phiKZ group. In the meantime, phages Lul11 and OBP do not exhibit resemblance among themselves or with any of earlier described phiKZ-like phages in respect to other traits; particularly, they have no detectable DNA homology. Note that phage Lu11 of P. putida var. Manila exhibits very slight homology with phage Lin68 of the family of P. aeruginosa phiKZ-like phages detected only in blot hybridization. This suggests the possible involvement of these phages in interspecies recombination ("gene shuffling") between phages of various bacterial species. Results of partial sequencing of phage genomes confirmed the phylogenetic relatedness of phage OBP to phages of the phiKZ-supergroup, whereas phage Lu11 most probably belongs to a novel species that is not a member of supergroup phiKZ composition. The results of the study are discussed in terms of the evolution of these phages.     

Phenogenetic characterization of a group of giant Phi KZ-like bacteriophages of Pseudomonas aeruginosa

A comparative study was made of a group of Pseudomonas aeruginosa virulent giant DNA bacteriophages similar to phage phi KZ in several genetic and phenotypic properties (particle size, particle morphology, genome size, appearance of negative colonies, high productivity, broad spectrum of lytic activity, ability to overcome the suppressing effect of plasmids, absence of several DNA restriction sites, capability of general transduction, pseudolysogeny). We have recently sequenced the phage phi KZ genome (288,334 bp) [J. Mol. Biol., 2002, vol. 317, pp. 1-19]. By DNA homology, the phages were assigned to three species (represented by phage phi KZ, Lin68, and EL, respectively) and two new genera (phi KZ and EL). Restriction enzyme analysis revealed the mosaic genome structure in four phages of the phi KZ species (phi KZ, Lin21, NN, and PTB80) and two phages of the EL species (EL and RU). Comparisons with respect to phage particle size, number of structural proteins, and the N-terminal sequences of the major capsid protein confirmed the phylogenetic relatedness of the phages belonging to the phi KZ genus. The origin and evolution of the phi KZ-like phages are discussed. Analysis of protein sequences encoded by the phage phi KZ genome made it possible to assume wide migration of the phi KZ-like phages (wandering phages) among various prokaryotes and possibly eukaryotes. Since the phage phi KZ genome codes for potentially toxic proteins, caution must be exercised in the employment of large bacteriophages in phage therapy.     

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.     

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.     

Isolation of new Pseudomonas tolaasii bacteriophages and genomic investigation of the lytic phage BF7

Sixteen lytic bacteriophages that infect Pseudomonas tolaasii LMG 2342(T) were isolated from smashed sporocarps of oyster mushroom (Pleurotus ostreatus) showing necrotic symptoms. On the basis of the host range investigation of the phages, they have wide infection abilities against the genus Pseudomonas, mainly in the case of phages Bf3, Bf7, Bf10, and Bf15. Molecular investigations have revealed that they all have dsDNA genomes about 40?kbp in size. Identical restriction patterns resulting from restriction enzyme analysis suggest that the isolates probably belong to the same phage species. However, there was a difference between these phage isolates in their infecting abilities. Phage isolate Bf7 was investigated and characterized more deeply. Morphological characterization of Bf7 by transmission electron microscopy (TEM) has shown that it has a short, noncontractile tail, an icosahedral phage head, and the size is about 60?nm in diameter, suggesting that it belongs to the Podoviridae family. Complete genome sequence analysis of the Bf7 phage isolate revealed a 40?058?bp genome, 58.4% G+C content, 46 open reading frames encoding different proteins showing homology to proteins of the bacteriophage Caulobacter crescentus φCd1 from the Podoviridae family. On the basis of these results and comparative genomic studies, we classified the Bf7 phage to the subfamily of Autographivirinae, φKMV-like phages.     

Use of lambda gt11 to isolate genes for two pseudorabies virus glycoproteins with homology to herpes simplex virus and varicella-zoster virus glycoproteins.

A library of pseudorabies virus (PRV) DNA fragments was constructed in the expression cloning vector lambda gt11. The library was screened with antisera which reacted with mixtures of PRV proteins to isolate recombinant bacteriophages expressing PRV proteins. By the nature of the lambda gt11 vector, the cloned proteins were expressed in Escherichia coli as beta-galactosidase fusion proteins. The fusion proteins from 35 of these phages were purified and injected into mice to raise antisera. The antisera were screened by several different assays, including immunoprecipitation of [14C]glucosamine-labeled PRV proteins. This method identified phages expressing three different PRV glycoproteins: the secreted glycoprotein, gX; gI; and a glycoprotein that had not been previously identified, which we designate gp63. The gp63 and gI genes map adjacent to each other in the small unique region of the PRV genome. The DNA sequence was determined for the region of the genome encoding gp63 and gI. It was found that gp63 has a region of homology with a herpes simplex virus type 1 (HSV-1) protein, encoded by US7, and also with varicella-zoster virus (VZV) gpIV. The gI protein sequence has a region of homology with HSV-1 gE and VZV gpI. It is concluded that PRV, HSV, and VZV all have a cluster of homologous glycoprotein genes in the small unique components of their genomes and that the organization of these genes is conserved.     

Complete genome sequences of three Erwinia amylovora phages isolated in north america and a bacteriophage induced from an Erwinia tasmaniensis strain

Fire blight, a plant disease of economic importance caused by Erwinia amylovora, may be controlled by the application of bacteriophages. Here, we provide the complete genome sequences and the annotation of three E. amylovora-specific phages isolated in North America and genomic information about a bacteriophage induced by mitomycin C treatment of an Erwinia tasmaniensis strain that is antagonistic for E. amylovora. The American phages resemble two already-described viral genomes, whereas the E. tasmaniensis phage displays a singular genomic sequence in BLAST searches.     

Comparison of genomes of new gigantic Pseudomonas aeruginosa phages from native populations from different regions

To study the genome diversity of bacteriophages from geographically distant natural populations, new giant phi KZ-like Pseudomonas aeruginosa phages isolated in two different regions were compared with earlier known phages of three species (phi KZ, Lin68, EL). A broad spectrum of lytic activity was demonstrated for all phi KZ-like phages. Phages of the phi KZ species proved to be common in natural populations of various regions, while IL- and Lin68-related phages were extremely rare. Most phi KZ-related phages had unique DNA restriction patterns, but the differences between these were only minor, and the genomes did not contain nonhomologous fragments. The spectrum of capsid polypeptides proved to be conserved in each species, and was proposed as a character necessary and sufficient for express classification of phages with an accuracy of species. Phages isolated in different geographical regions showed no substantial difference. Some phages only slightly differing in DNA restriction pattern from phi KZ may be used to study the origin of phi KZ genes coding for orthologs of proteins of unrelated species (other phages, pathogenic bacteria, eukaryotes).     

Nucleotide sequence and genetic organization of the genome of the N-specific filamentous bacteriophage IKe. Comparison with the genome of the F-specific filamentous phages M13, fd and f1

The nucleotide sequence and genetic organization of the genome of the N-specific filamentous single-stranded DNA phage IKe has been established and compared with that of the F-specific filamentous phages M13, fd and f1 (Ff). The IKe DNA sequence comprises 6883 nucleotides, which is 476 (475) nucleotides more than the nucleotide sequence of the Ff genome. The data indicate that IKe and Ff have evolved from a common ancestor (overall homology approx. 55%) and that their genomes contain ten homologous genes, the order of which is identical. Similar to Ff, the major coat protein and the gene III-encoded pilot protein of IKe are synthesized via precursor molecules. The extent of homology between the genes of IKe and Ff differs significantly from one gene to another. Genes that code for viral capsid proteins are less homologous than genes whose products are involved in the processes of DNA replication and phage morphogenesis. During evolution, large nucleotide sequence rearrangements have occurred in the gene (gene III) whose product is needed for the attachment of the virion to the conjugative pili of the host cell, suggesting that these rearrangements have led to phages with different host specificities. Extensive nucleotide sequence homology was noted between the structural elements involved in DNA replication and phage morphogenesis, indicating that the mechanisms involved in DNA replication and morphogenesis are highly conserved.     

Study of morphology and genome structure of Pseudomonas putida bacteriophages for their classification

A group of 27 bacteriophages specific for Pseudomonas putida strains PpG1 and PpN has been isolated. The phages were characterized and compared with the previously described virulent (pf 16, af, tf and PMW) and temperate (PP56 and PP71) phages. The new phages belong to B1 and C1 morphotypes, according to Ackerman's classification. Phage DNAs were digested with several endonucleases; the molecular weights and homology of the DNAs were determined. All phages of P. putida isolated up to now were distributed into 10 species (groups), on the basis of particle morphology, genome size and the results of homology studies. Recombination processes are believed to participate in formation of phages belonging to certain species.     

Two-dimensional gel electrophoretic analysis of the HindIII 1.8-kb repetitive-sequence family in the human genome

The structure and organization of a human repetitive DNA family containing the HindIII 1.8-kb repetitive sequence were studied, using two-dimensional (2D) gel electrophoresis. The HindIII 1.8-kb sequence proved to be part of a repetitive sequence about 5 kb long and interspersed on the genome. The long repetitive sequence family contained several subgroups, as based on polymorphism of the restriction site. Recombinant phages containing the long repetitive sequence were isolated from the human genomic DNA library. Heteroduplex and restriction analysis showed that the structure of the repetitive sequence carried by the phages was close to that expected from 2D gel electrophoretic analysis. The 2D gel electrophoretic analysis was shown to be a reliable and useful approach for surveying and mass analysis of repetitive sequence families.     

Diversity of temperate bacteriophages induced in bovine and ovine Mannheimia haemolytica isolates and identification of a new P2-like phage

The diversity of temperate bacteriophages was examined in 32 Mannheimia haemolytica, six Mannheimia glucosida and four Pasteurella trehalosi isolates. Phage particles were induced and identified by electron microscopy in 24 (75%) M. haemolytica isolates, but in only one (17%) M. glucosida and one (25%) P. trehalosi isolate. The M. haemolytica phages were relatively diverse as seven Siphoviridae, 15 Myoviridae and two Podoviridae-like phages were identified; the Myoviridae-type phages also exhibited structural variation of their tails. The bacteriophages induced in M. glucosida and P. trehalosi were of the Myoviridae type. Restriction endonuclease (RE) analysis identified nine distinct RE types among the M. haemolytica bacteriophages, providing further evidence of their relative diversity. A limited number of phages caused plaques on indicator strains and the phages exhibited a narrow host range. A subgroup of 11 bovine serotype A1 and A6 isolates contained Myoviridae-type phages of the same RE type (type A), but these differed in their abilities to infect and form plaques on the same panel of indicator strains. A P2-like phage (phiPHL213.1), representative of the RE type A phages, was identified from the incomplete M. haemolytica genome sequence. The phiPHL213.1 genome contains previously unidentified genes and represents a new member of the P2 phage family.     

Analysis of Arabidopsis genome sequence reveals a large new gene family in plants

A detailed analysis of the currently available Arabidopsis thaliana genomic sequence has revealed the presence of a large number of open reading frames with homology to the stigmatic self-incompatibility (S) genes of Papaver rhoeas. The products of these potential genes are all predicted to be relatively small, basic, secreted proteins with similar predicted secondary structures. We have named these potential genes SPH (S-protein homologues). Their presence appears to have been largely missed by the prediction methods currently used on the genomic sequence. Equivalent homologues could not be detected in the human, microbial, Drosophila or C. elegans genomic databases, suggesting a function specific to plants. Preliminary RT-PCR analysis indicates that at least two members of the family (SPH1, SPH8) are expressed, with expression being greatest in floral tissues. The gene family may total more than 100 members, and its discovery not only illustrates the importance of the genome sequencing efforts, but also indicates the extent of information which remains hidden after the initial trawl for potential genes.