Skip to main content

The genome and proteome of Serratia bacteriophage η which forms unstable lysogens



Serratia marcescens phage η is a temperate unclassified member of the Siphoviridae which had been reported as containing hypermodified guanine residues.


The DNA was characterized by enzymatic digestion followed by HPLC analysis of the nucleoside composition, and by DNA sequencing and proteomic analysis. Its ability to form stable lysogens and integrate was also investigated.


Enzymatic digestion and HPLC analysis revealed phage η DNA did not contain modified bases. The genome sequence of this virus, determined using pyrosequencing, is 42,724 nucleotides in length with a mol% GC of 49.9 and is circularly permuted. Sixty-nine putative CDSs were identified of which 19 encode novel proteins. While seven close genetic relatives were identified, they shared sequence similarity with only genes 40 to 69 of the phage η genome, while gp1 to gp39 shared no conserved relationship. The structural proteome, determined by SDS-PAGE and mass spectrometry, revealed seven unique proteins. This phage forms very unstable lysogens with its host S. marcescens.


Serratia marcescens is a Gram-negative opportunistic human pathogen belonging to the family Enterobacteriaceae. In humans it is of particular concern for patients with indwelling medical devices [1]. Outbreaks are almost exclusively caused by contaminated medical solutions or equipment, with S. marcescens even growing in antiseptic solutions [2]. Only recently, community acquired infections have been identified [3, 4], which highlight the growing importance of this pathogen. Treatment of infections is possible with combinations of antibiotics; however the number of clinical strains resistant to multiple antibiotics is on the rise [5, 6], including isolates which produce extended spectrum beta-lactamases [7, 8]. Phage therapy might prove a viable alternative, but historically phages infecting S. marcescens were isolated solely for use in bacterial typing systems [9, 10] or for the transductional mapping of the genome [11, 12] and their biology remains poorly understood [13].

There are 87 named Serratia bacteriophages listed in Bacteriophage Names 2000 ( belonging to the Siphoviridae[14], Myoviridae[15], and Podoviridae[12] families of tailed phages. Since 2000, only a limited number of additional Serratia phages have been isolated and characterized. Phages PPV [16] and KSP100 [17] belong to family Podoviridae; the former being a tentative member of the T7likevirus genus, while the latter KSP100 is a member of the Phieco32likevirus genus. Only one additional siphovirus, phage SM701, has been described [18]. All other newly characterized phages are all members of the Myoviridae. These include the virulent generalized transducing phage ϕIF3 (phiIF3) [11], ϕMAM1 (phiMAM1) [19] which is related to the newly proposed Viunalikevirus[20], and phages KSP20 and KSP90 [17]. Phylogenetic analysis based upon limited sequence data suggests that KSP20 with its 32 kb genome is part of the Peduovirinae[21], while KSP90 is part of the T4 superfamily.

S. marcescens bacteriophages were chosen as the focus of our research because of the growing importance of the host and its increasing prevalence as a nosocomial pathogen. Bacteriophage η was selected for further study. This virus was isolated from an overnight culture supernatant of S. marcescens ES, in 1960 in Germany [22]. Of considerable interest was von Lohr’s observation that the phage DNA displayed a discrepancy between the mol% GC calculated from melting temperature (Tm) and buoyant density determinations suggesting the presence of a hypermodified base. Preliminary data suggested a hypermodified guanine [22].

Resrults and discussion

Phage morphology and basic biology

Phage η is a B1 siphovirus [23, 24] with an icosahedral head of 60 nm and a tail of 112 × 8 nm which ends in a base plate with three conspicuous spikes (Figure 1). As such, phage η is morphologically identical to the Jersey species of Salmonella phages [23, 24].

Figure 1

Electron micrograph of phage η negatively stained with 2% uranyl acetate. The bar indicates 100 nm.

The phage produces variably sized plaques when grown on S. marcescens using the double-agar overlay technique. The larger size plaque was approximately 3 mm in diameter, and the smaller only 1.5 mm. Both plaque types were turbid, with greater turbidity in the plaque centers. The distribution of turbidity in both plaque types was the same and they arose equally from subcultures of both small and large plaque types.

Phage η forms highly unstable lysogens with its host

Phage-insensitive bacterial strains were selected from turbid plaques and serially propagated to remove any traces of contaminating exogenous phage. Lysogenic strains were distinguished from resistant strains by testing for the ability of chloroform-treated culture supernatants to form plaques on the wild type strain, for resistance to super-infection, and PCR reactions with phage η specific primers (data not shown). Initial subcultures were resistant to superinfection, spontaneous release of phage in absence of an induction agent to a titre of 107 PFU/ml, and 109 PFU/ml when induced with mitomycin C. In addition, phage-specific amplicons were observed with host DNA preparations. Most of the progeny from subsequent subcultures were phage sensitive, though phage production was also still measurable in the bulk culture. This combination of phage production and sensitivity suggests the sub-cultures were reverting to non-lysogenic cells. We choose to dismiss the possibility that this phage displays pseudolysogeny [15, 25, 26] based upon the isolation history of this virus and its inducibility by Mitomycin C; and, favour unstable lysogeny under laboratory culture conditions. Similar observations have been made previously with some Vibrio and Clostridium phages [2729].

Phage η DNA does not contain hypermodified guanine residues

Previous research had suggested the presence of a modified base in the DNA of phage η [22]. This was investigated by digestion of the DNA to nucleosides using a combination of DNase I, phosphodiesterase I and alkaline phosphatase coupled with HPLC separation of the products. We expected to observe a shift in the elution time of one of the four nucleoside peaks or the appearance of a fifth peak in the elution profile [30]. Neither of these were observed (Additional file 1: Figure S1), therefore there was no evidence to support the presence of a modified base. Phage molecular biologists no longer routinely calculate the mol% GC based upon buoyant density or melting temperature [30] and restriction analyses are also becoming a thing of the past. As a consequence there is a strong possibility that DNA modifications may not be discovered since pyrosequencing reads through most modifications (Kropinski, unpublished observations). Third generation sequencing techniques may address this potential problem.

Phage η DNA is circularly permuted and terminally redundant

Determination of the physical location and nature of the genome ends can offer insight into the mechanisms of DNA replication and packaging that a bacteriophage employs [31]. Restriction digestion profiles of phage η DNA matched in silico predictions of a circular DNA molecule, and there was no difference in the profiles generated with ligated and unligated DNA samples (Additional file 2: Figure S2). Furthermore, there was no evidence of sub-molar fragments. Lastly, pulsed-field gel electrophoresis revealed one DNA band with no smearing or concatemeric molecules. These all suggest the impossible - the genome of phage η is circular, necessitating further study.

Exonuclease digestions were carried out to processively truncate the DNA. When carried out before standard restriction digestion, the DNA bands containing the ends of the genome would be expected to become progressively smaller with increased exonuclease digestion time [31, 32]. Time-controlled Bal31 exonuclease treatment followed by complete restriction digestion did not result in the decrease in molecular weight of any of the DNA restriction fragments (Figure 2A), In contrast for the control λ DNA, band mass of the left (band X) and right (band Y) termini were progressively reduced (Figure 2B). Extended Bal31 exonuclease treatment caused a decrease in η DNA concentration over time. As the exonuclease only acts on linear DNA this suggests that the η DNA molecule is linear, as expected, and not a covalently closed circle as evidenced by the restriction mapping. As an additional approach, 16 overlapping primer pairs were designed to span the genome and give products of approximately 4 kb. Each of the primer pairs produced a unique length PCR product which was the same size as predicted. The sum of these data suggests that the dsDNA genome is linear with a high degree of circular permutation.

Figure 2

Determination of the nature of the termini in phage η DNA. A. Bal31 exonuclease treated phage λ DNA with subsequent EcoRI restriction digestion, Lanes are labeled according to the number of minutes the DNA was incubated with the Bal31 before EcoRI digestion with the bands which decrease are boxed. B. phage η DNA was digested with exonuclease and subsequently digested with EcoRI. All lanes labeled L in both 2A and 2B are loaded with Ultra Ranger DNA Ladder ™ with size markers indicated in kb.

Genomic organization and annotation

The genome of phage η is 42,724 bp in length with a mol% GC of 49.9. This value is significantly less than that of the host, S. marcescens at 58 mol% GC, which is somewhat unusual for temperate phages where the overall base composition of the host and phage are usually remarkably similar (Kropinski, unpublished results). A unique 19 bp direct repeat was identified (5′-ATTGCAACTTATTTGTTTA-3′), which was found in four locations in the genome (Additional file 3: Table S1), and is possibly involved in regulation of expression.

The orientation of genome annotation was chosen so that the majority of the genes were on the plus strand. Sixty-nine putative CDSs, identified during the genome annotation, are divided into four operons by four promoters and rho-independent terminators (Figure 3). Two potential morons (genes 45 and 60) also divide the operons, each associated with their own promoter and terminator [33]. A putative oriC was identified using Ori-Finder.

Figure 3

Genetic map of η. The genes are color coded according to function with red: novel, light green: prophage or bacterial related, dark green: phage related, black: lysis cassette, orange: recombination, brown: replication, yellow: packaging related, dark blue; confirmed structural, light blue: hypothetical, white: morons, grey: possible DNA binding and Kil. Promoters (P), terminators (T), direct repeats (R) and the ori are also included. Promoters are indicated with black arrows on stalks, while rho-independent terminators are indicated . The GenomeVx software package was used to construct this figure.

The 69 proteins were grouped into six categories based on the number of homologs and the degree of amino acid conservation. Thirteen predicted proteins (90-100% identical sequence coverage) were very highly conserved to bacterial proteins with no homologs present in the viral database (data not shown). Nineteen predicted proteins were entirely novel. There were seven structural proteins confirmed using mass spectrometry (Additional file 4: Table S2). The fourth category contained 25 proteins which were conserved hypothetical proteins. Though based on genome organization of phage η and related phages, a structural or morphogenesis role could be hypothesized for nine of these proteins. Of the proteins with no predicted function, four DNA-binding proteins were tentatively identified, none of them related to previously characterized phage regulatory proteins. Functional predictions could be made for the remaining 13 gene products. These functions were associated mainly with the lytic replication cycle, with protein functions involved in DNA replication, DNA packaging, structural proteins and host cell lysis. No integrase was identified during the annotation, although recombination related proteins and putative DNA-binding proteins could be identified (Additional file 3: Table S1). Some notable genes will be discussed in greater detail.

Phage integration

There are two lysogenic pathways open to temperate phages - integration into the host genome or maintenance as a plasmid. The former, as exhibited by coliphage lambda, involves recombination between homologous integration (att) sites on the phage and host genome mediated by integrase. Integration can also involve, as it does with the transposable phage such as Mu, a transposase, leading to nonhomologous recombination, and random insertion into the host genome. Phage Eta lacks both integrase and transposase homologs. The genome of this virus shares a 33 bp sequence (GTGGAGGTGCGTGATGCCAGCAAA TGAACTGAA) located between residues 15009 and 15041 within the complete genome sequences of Serratia sp. AS13 (CP002775), Serratia sp. AS12 (CP002774) and Serratia plymuthica AS9 (CP002773). In each case, the homologous sequence is located within a prophage related to enterobacterial phage cdtI identified using PHAST[33]. In the case of strain AS13, this site was within the coding sequence of hypothetical protein SerAS13_3448. Whether this serves as an integration site is not known; however, it is long enough that it could function in homologous recombination.

The alternative strategy is displayed by coliphage N15 [34], Yersinia phage PY54 [35], Klebsiella phage ΦKO2 [36], and Halomonas aquamarina phage ΦHAP-1 [37]. These phages encode protelomerase, along with partitioning proteins ParA and ParB. Coliphage P1 also possesses ParAB homologs. A ParB homolog was identified in Eta (gene 11; Additional file 3: Table S1), though we can find no evidence for a ParS sequence nor a ParA homolog. It is possible that the prophage is maintained as a plasmid-type element, and is segregated using a combination of the phage ParB homolog (centromere-binding protein) and host factors.

DNA packaging

Packaging of phage DNA into the precapsid during phage replication is carried out by the phage-encoded terminase holoenzyme, which is typically a heterodimer composed of large and small subunit. The large subunit of the terminase (gene 40) of phage η was well conserved with previously identified subunits (Additional file 3: Table S1), unlike the small subunit (gene 39), which had an entirely novel N-terminus. Basic features and domain organization in terminase enzyme complexes are shared in the tailed phages [38]. The large terminase subunit responsible for the ATPase and nuclease packaging-associated enzymatic functions, with the small terminase is involved in sequence specific substrate recognition and holoenzyme formation [39].

Structural proteins

Seven distinct bands in the SDS-PAGE gel of phage η were identified using high resolution UPLC LTQ-FT mass spectrometry (Additional file 4: Table S2). Reliable identification was obtained by mass measurements of the peptides and high protein sequence coverage of the proteolytic fragments observed by tandem mass spectrometric (MS/MS) analyses. In the case of tailspike protein, the N-terminus of the protein was modified by acetylation (Additional file 4: Table S2). In addition to these high abundance proteins, an additional nine proteins could be annotated as likely as structural or morphogenesis proteins (Additional file 3: Table S1). It is possible that these proteins are either present at low abundance which could not be detected on the SDS-PAGE gel. While otherwise a standard lysis cassette, that of phage η is unique in that it is located in the middle of the capsid proteins with confirmed structural proteins both upstream and downstream (Figure 3). Typically, the structural proteins are co-regulated in one operon, with the lysis cassette upstream, which allows for a simplified regulation.

Genomic relatives and classification

The structural cassette of phage η is the only segment of the genome which showed any degree of conserved relationship with phage genome sequences available in the NCBI database. Discontiguous megablast analysis revealed that phage η shared limited similarity (<22% sequence coverage) with Escherichia phages K1H [40] and K1G [40], Salmonella phages SETP3 [41], vB_SenS-Ent1 [42], SE2 [43], wksl3 [44], and SS3e [45]. The homologous regions were restricted to 22-43kb on the η genome, with the average shared sequence identity above 80% (Additional file 3: Table S1).


Serratia phage η was originally characterized as a temperate siphovirus possessing hypermodified guanine residues in its DNA. The putative presence of hypermodified bases in a temperate phage creates a conceptual problem since expression of modifying enzymes would probably modify the host genome with catastrophic results. Our chromatographic results do not support the existence of modified bases in the DNA of this virus. Furthermore, the genomic results suggest that this phage is not classically temperate, since it lacks a demonstrable repressor and integrase homologs, but is capable of establishing an unstable relationship with its host which can best be described as unstable lysogeny. The DNA genome of phage η is terminally redundant and highly circularly permuted. The latter observation has led other researchers [19, 4649] to the mistaken assumption that their phage genomes are circular, which has never been observed among the Caudovirales, and indeed would pose problems in DNA packaging.

Comparative genomic analysis reveals that Serratia phage η is a unique virus which, except for the structural gene module, is unrelated to any of the current population of phages which have the complete genomes sequenced. Interestingly, the genome of its structural homolog, Salmonella phage Jersey, has recently been sequenced as part of a project to characterize the “Jerseylikevirus” genus (H. Anany, personal communication), and the only relationship between these viruses is within the morphogenesis cassette.

Materials and methods

Bacteria and bacteriophages

S. marcescens CV/rc3 (HER 1311) and phage η were obtained from the Centre de Référence pour les Virus Bactériens Félix d’Hérelle (Université Laval, QC, Canada). The host was grown in Difco Luria-Bertani broth (LB; Fisher Scientific, Toronto, ON, Canada) at 30°C with shaking at 180 rpm or on LB agar plates (LB with 1.5% [wt/vol] agar) at 30°C. The bacteriophage titer was determined using the double agar overlay method [48]. Phage particles were partially purified by PEG precipitation and continuous CsCl purification [49]. Purified phage was dialyzed using a Pierce 10 K Slide-A-Lyzer cassette (Thermo Fisher Scientific, Rockford, IL, USA) against 10 mM Tris-0.1 mM EDTA (pH 8.0) (TE) buffer to remove CsCl. Purified phage was stored at 4°C.

Electron microscopy

Purified phages were sedimented by centrifugation at 25,000 g for 60 min, using a Beckman Coulter J-E (Palo Alto, CA) centrifuge and a JA19.1 fixed angle rotor. The pellet was washed twice under the same conditions in neutral ammonium acetate (0.1 M). Phages were deposited on copper grids with carbon-coated Formvar films, stained with 2% uranyl acetate (pH 4.5) and examined using a Philips EM 300 electron microscope.

DNA isolation, digestion and HPLC analysis: Phage

DNA was extracted by standard phenol chloroform extraction and ethanol precipitation [50]. The purified phage η genomic DNA was digested to the nucleoside level in preparation for HPLC analysis according to the protocol described by [30] using DNase I, phosphodiesterase I and alkaline phosphatase (Worthington Biochemical Corp. Vassar, NJ, USA). The number of units of phosphodiesterase I were increased from 0.01 units in the protocol to 0.1 units. The phosphodiesterase I and alkaline phosphatase digestions were done in a water bath at 25°C. The digested DNA was filtered using an Ultrafree–MC membrane with a pore size of 0.1 μm (Millipore Corp, Billerica, MA, USA) to remove particles in preparation for loading onto the HLPC column. Standard nucleosides were purchased from Sigma-Aldrich (Sigma-Aldrich, Oakville, Canada).

The protocol developed was based on the method outlined by [30]. The analysis was performed on an Agilent 1100 series LC/MSD_SL Trap system. Samples (100 μl) were injected into the LC/MSD system through an Agilent 1100 series thermostated auto-sampler. Separations were carried out on a 5 μm Agilent ZORBAX SB C-18 column (4.6 × 150 mm) at a flow rate of 0.75 ml/min. Two mobile phase solvents were used; Solvent A was 0.1% formic acid in deionized water; Solvent B was 100% HPLC grade acetonitrile. Dimethyl sulfoxide (DMSO) was used to dissolve the concentrated nucleosides. HPLC method: 98/2 A/B for five minutes then increasing to 95/5 A/B at ten minutes, 80/20 A/B at 12 minutes then returning to 100% A.

DNA sequencing and annotation

Phage DNA was digested, ligated and treated with Bal31 exonuclease according to manufacturer’s specifications (NEB) and resolved using a 1% agarose gel. Genome end determination was carried out following protocols published by Loessner et al., [32] and Casjens et al., [51]. Pulsed-field gel analysis was carried out following the protocol of Lingohr et al., [52].

The genome sequence was determined using 454 sequencing technology at the McGill University and Genome Quebec Innovation Centre (Montreal, QC, Canada) with 49.1 fold coverage. The genome was annotated using Kodon (Applied Maths, Austin, TX, USA) with the proteins screened for homologs using the BLASTP feature of Geneious Pro 6.2 (Biomatters Ld., Auckland, New Zealand). Their molecular mass and pI were calculated using the Batch MW and pI Finder server at Potential protein motifs were identified using the NCBI Batch CD-Search tool at with an E-value threshold of 0.0001. Transmembrane domains were screened for using TMHMM [53] and Phobius [54]. In certain cases HHpred [55, 56] was employed to analyze the protein sequences. The annotated sequence was deposited with GenBank under accession number KC460990 (GI:511624446). RNA secondary structure prediction was used to support the identification of rho-independent terminators [57]. The oriC was identified using Ori-Finder [58]. The gene diagram was constructed using Genome Vx [59].

Structural protein analysis

Phage particles which had been purified through CsCl twice were used for SDS-PAGE structural protein analysis. Purified whole phages were denatured at 100°C for 5 min. in Laemmli Sample Buffer supplemented with 100 mM 2-mercaptoethanol. The structural proteins were resolved in a Bio-Rad Ready 5-15% Tris-HCl gel, in Running Buffer (25 mM Tris, 192 mM glycine, 0.1% SDS, pH 8.3 (Bio-Rad)) at 150 V for an hour. Visible bands were excised, and then digested in-gel by sequencing-grade trypsin. The cleaved peptide fragments of the proteins were extracted using acentonitrile/0.1% trifluoroacetic acid (TFA) (v/v, 60:40), then dried by a Savant vacuum centrifuge (Thermo Fisher Scientific, Nepean, Ontario). In the subsequent mass spectrometric analyses, the peptides were reconstituted in 8 μl of 0.2% formic acid (FA) and identified using a nanoAcquity ultra-performance liquid chromatograph (UPLC, Waters, Milford, MA) – linear ion-trap Fourier transform ion cyclotron resonance (LTQ-FT ICR, Thermo Fisher, San Jose, CA) mass spectrometer. The sample was normally trapped by a reverse-phase (RP) Symmetry C18 column (180 μm i.d. × 20 mm length, 5 μm) at 5 μl/min of solvent A (0.1% FA) for 3 minutes, and then separated through a C18 analytical column (100 μm i.d. × 100 mm, 1.7 μm, BEH 130) at 400 nl/min for 65 minutes. UPLC gradient was set up to a linear gradient from 5% to 30% solvent B (0.1% FA in acetonitrile), followed by to 85% solvent B over 10 min. for peptide elution. FT-MS scans were acquired with high resolution (100,000) MS from m/z 300 to 2000, and MS/MS measurements in linear ion-trap mode by data-dependent scans of the top eight intense precursor ions at multiply charged states of 2+, 3+ and 4+. Dynamic exclusion was set to a period of time at 180 s.

Protein identification was performed using Mascot Server (version 2.3.0, Matrix Science, London, UK), and UPLC MS/MS raw data were searched against the in-house database of the protein sequence derived from the annotated genome of phage η. The search parameters were restricted to tryptic peptides for a maximum of 2 missed cleavages. Cysteine carbamidomethylation was designated as a fixed modification, and deamidation of asparagine and glutamine, methionine oxidation were considered as variable modifications. Mass tolerances were set up to 10 ppm for the FT-MS ions and 1 Da for ion trap MS/MS fragment ions. Peptide assignments were filtered by an ion score cut off at 20, and the identified MS/MS spectra were also verified manually.



Basic local alignment search tool


Gene product


High pressure liquid chromatography


Not applicable


Coding sequence


Plaque forming unit, a measure of the number of viable viral particles


Denaturing polyacrylamide gel electrophoresis


Transmembrane prediction using hidden Markov models.


  1. 1.

    Wolcott RD, Gontcharova V, Sun Y, Zischakau A, Dowd SE: Bacterial diversity in surgical site infections: not just aerobic cocci any more. J Wound Care 2009, 18: 317-323.

    PubMed  CAS  Article  Google Scholar 

  2. 2.

    Bosi C, Davin-Regli A, Charrel R, Rocca B, Monnet D, Bollet C: Serratia marcescens nosocomial outbreak due to contamination of hexetidine solution. J Hosp Infect 1996, 33: 217-224. 10.1016/S0195-6701(96)90005-5

    PubMed  CAS  Article  Google Scholar 

  3. 3.

    Perez C, Fujii Y, Fauls M, Hummel J, Breitschwerdt E: Fatal aortic endocarditis associated with community-acquired Serratia marcescens infection in a dog. J Am Anim Hosp Assoc 2011, 47: 133-137. 10.5326/JAAHA-MS-5616

    PubMed  Article  Google Scholar 

  4. 4.

    Pinna A, Usai D, Sechi LA, Carta A, Zanetti S: Detection of virulence factors in Serratia strains isolated from contact lens-associated corneal ulcers. Acta Ophthalmol 2009, 89: 382-387.

    PubMed  Article  Google Scholar 

  5. 5.

    Dowzicky MJ, Park CH: Update on antimicrobial susceptibility rates among gram-negative and gram-positive organisms in the United States: results from the Tigecycline Evaluation and Surveillance Trial (TEST) 2005 to 2007. Clin Ther 2008, 30: 2040-2050. 10.1016/j.clinthera.2008.11.006

    PubMed  CAS  Article  Google Scholar 

  6. 6.

    Samonis G, Vouloumanou EK, Christofaki M, Dimopoulou D, Maraki S, Triantafyllou E, Kofteridis DP, Falagas ME: Serratia infections in a general hospital: characteristics and outcomes. Eur J Clin Microbiol Infect Dis 2011, 30: 653-660. 10.1007/s10096-010-1135-4

    PubMed  CAS  Article  Google Scholar 

  7. 7.

    Crivaro V, Bagattini M, Salza MF, Raimondi F, Rossano F, Triassi M, Zarrilli R: Risk factors for extended-spectrum beta-lactamase-producing Serratia marcescens and Klebsiella pneumoniae acquisition in a neonatal intensive care unit. J Hosp Infect 2007, 67: 135-141. 10.1016/j.jhin.2007.07.026

    PubMed  CAS  Article  Google Scholar 

  8. 8.

    Rizvi M, Fatima N, Shukla I, Malik A: Epidemiology of extended spectrum beta-lactamases in Serratia and Citrobacter species in North India. Indian J Pathol Microbiol 2010, 53: 193-194. 10.4103/0377-4929.59237

    PubMed  Article  Google Scholar 

  9. 9.

    Hamilton RL, Brown WJ: Bacteriophage typing of clinically isolated Serratia marcescens . Appl Microbiol 1972, 24: 899-906.

    PubMed  CAS  PubMed Central  Google Scholar 

  10. 10.

    Wassermann MM, Seligmann E: Serratia marcescens bacteriophages. J Bacteriol 1953, 66: 119-120.

    PubMed  CAS  PubMed Central  Google Scholar 

  11. 11.

    Petty NK, Foulds IJ, Pradel E, Ewbank JJ, Salmond GP: A generalized transducing phage (phiIF3) for the genomically sequenced Serratia marcescens strain Db11: a tool for functional genomics of an opportunistic human pathogen. Microbiology 2006, 152: 1701-1708. 10.1099/mic.0.28712-0

    PubMed  CAS  Article  Google Scholar 

  12. 12.

    Regue M, Fabregat C, Vinas M: A generalized transducing bacteriophage for Serratia marcescens . Res Microbiol 1991, 142: 23-27. 10.1016/0923-2508(91)90093-P

    PubMed  CAS  Article  Google Scholar 

  13. 13.

    Calendar R: Introduction. In Bacteriophages. Edited by: Calendar R. Interscience publishers Inc: New York; 2006:1-6.

    Google Scholar 

  14. 14.

    Casjens S, Winn-Stapley DA, Gilcrease EB, Morona R, Kuhlewein C, Chua JE, Manning PA, Inwood W, Clark AJ: The chromosome of Shigella flexneri bacteriophage Sf6: complete nucleotide sequence, genetic mosaicism, and DNA packaging. J Mol Biol 2004, 339: 379-394. 10.1016/j.jmb.2004.03.068

    PubMed  CAS  Article  Google Scholar 

  15. 15.

    Abedon ST: Disambiguating bacteriophage pseudolysogeny: an historical analysis of lysogeny, pseudolysogeny, and the phage carrier state. In Contemporary trends in bacteriophage research. Edited by: Adams HT. New York: Nova Science Publishers, Inc.; 2009:285-307.

    Google Scholar 

  16. 16.

    Friman VP, Hiltunen T, Jalasvuori M, Lindstedt C, Laanto E, Ormala AM, Laakso J, Mappes J, Bamford JK: High temperature and bacteriophages can indirectly select for bacterial pathogenicity in environmental reservoirs. PloS One 2011, 6: e17651. 10.1371/journal.pone.0017651

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  17. 17.

    Matsushita K, Uchiyama J, Kato S, Ujihara T, Hoshiba H, Sugihara S, Muraoka A, Wakiguchi H, Matsuzaki S: Morphological and genetic analysis of three bacteriophages of Serratia marcescens isolated from environmental water. FEMS Microbiol Lett 2009, 291: 201-208. 10.1111/j.1574-6968.2008.01455.x

    PubMed  CAS  Article  Google Scholar 

  18. 18.

    Yu L, Wen Z, Yang W, Li N, Wang J, Lu J, Li J: Isolation and characterization of Serratia marcescens phage. Wei sheng wu xue bao 2008, 48: 498-502.

    PubMed  CAS  Google Scholar 

  19. 19.

    Matilla MA, Salmond GP: Complete Genome Sequence of Serratia plymuthica Bacteriophage ϕMAM1. J Virol 2012, 86: 13872-13873. 10.1128/JVI.02702-12

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  20. 20.

    Adriaenssens EM, Ackermann HW, Anany H, Blasdel B, Connerton IF, Goulding D, Griffiths MW, Hooton SP, Kutter EM, Kropinski AM, Lee JH, Maes M, Pickard D, Ryu S, Sepehrizadeh Z, Shahrbabak SS, Toribio AL, Lavigne R: A suggested new bacteriophage genus: “Viunalikevirus”. Arch Virol 2012, 157: 2035-2046. 10.1007/s00705-012-1360-5

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  21. 21.

    Lavigne R, Darius P, Summer EJ, Seto D, Mahadevan P, Nilsson AS, Ackermann HW, Kropinski AM: Classification of Myoviridae bacteriophages using protein sequence similarity. BMC Microbiol 2009, 9: 224. 10.1186/1471-2180-9-224

    PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Pons FW: Untersuchung der DNS einiger Serratiastämme und deren Phagen. Biochemische Zeit 1966, 346: 26-40.

    CAS  Google Scholar 

  23. 23.

    Ackermann HW, Gershman M: Morphology of phages of a general Salmonella typing set. Res Virol 1992, 143: 303-310.

    PubMed  CAS  Article  Google Scholar 

  24. 24.

    Bradley DE: The morphology and physiology of bacteriophages as revealed by the electron microscope. J R Microsc Soc 1965, 84: 257.

    PubMed  CAS  Google Scholar 

  25. 25.

    Łoś M, Węgrzyn G: Pseudolysogeny. Adv Virus Res 2012, 82: 339-349.

    PubMed  Article  Google Scholar 

  26. 26.

    Miller RV, Day MJ: Contribution of lysogeny, pseudolysogeny and starvation to phage ecology. In Bacteriophage ecology: population growth, evolution, and impact of bacterial viruses. advances in molecular and cellular microbiology series. Edited by: Abedon ST. Cambridge, UK: Cambridge University Press; 2008:114-146.

    Google Scholar 

  27. 27.

    Khemayan K, Pasharawipas T, Puiprom O, Sriurairatana S, Suthienkul O, Flegel TW: Unstable lysogeny and pseudolysogeny in Vibrio harveyi siphovirus-like phage 1. Appl Environ Microbiol 2006, 72: 1355-1363. 10.1128/AEM.72.2.1355-1363.2006

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  28. 28.

    Sakaguchi Y, Hayashi T, Kurokawa K, Nakayama K, Oshima K, Fujinaga Y, Ohnishi M, Ohtsubo E, Hattori M, Oguma K: The genome sequence of Clostridium botulinum type C neurotoxin-converting phage and the molecular mechanisms of unstable lysogeny. Proc Natl Acad Sci U S A 2005, 102: 17472-17477. 10.1073/pnas.0505503102

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  29. 29.

    Pasharawipas T, Thaikua S, Sriurairatana S, Ruangpan L, Direkbusarakum S, Manopvisetcharean J, Flegel TW: Partial characterization of a novel bacteriophage of Vibrio harveyi isolated from shrimp culture ponds in Thailand. Virus Res 2005, 114: 63-69. 10.1016/j.virusres.2005.05.012

    PubMed  CAS  Article  Google Scholar 

  30. 30.

    Manderville RAM, Kropinski AM: Approaches to the compositional analysis of DNA. In Bacteriophages: methods and protocols, volume 2: molecular and applied aspects. Edited by: Clokie MRJ, Kropinski AM. Totowa, NJ: Humana Press; 2009:11-18.

    Google Scholar 

  31. 31.

    Casjens SR, Gilcrease EB: Determining DNA packaging strategy by analysis of the termini of the chromosomes in tailed-bacteriophage virions. In Bacteriophages: methods and protocols, volume 2: molecular and applied aspects. Edited by: Clokie MRJ, Kropinski AM. Totowa, NJ: Humana Press; 2009:91-111.

    Google Scholar 

  32. 32.

    Loessner MJ, Inman RB, Lauer P, Calendar R: Complete nucleotide sequence, molecular analysis and genome structure of bacteriophage A118 of Listeria monocytogenes: implications for phage evolution. Mol Microbiol 2000, 35: 324-340. 10.1046/j.1365-2958.2000.01720.x

    PubMed  CAS  Article  Google Scholar 

  33. 33.

    Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS: PHAST: a fast phage search tool. Nucleic Acids Res 2011,39(Web Server issue):W347-W352.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  34. 34.

    Ravin NV: N15: The linear plasmid prophage. In The Bacteriophages. Edited by: Calendar R. New York, NY: Oxford University Press; 2006:448-456.

    Google Scholar 

  35. 35.

    Hammerl JA, Klein I, Appel B, Hertwig S: Interplay between the temperate phages PY54 and N15, linear plasmid prophages with covalently closed ends. J Bacteriol 2007, 189: 8366-8370. 10.1128/JB.01066-07

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  36. 36.

    Casjens SR, Gilcrease EB, Huang WM, Bunny KL, Pedulla ML, Ford ME, Houtz JM, Hatfull GF, Hendrix RW: The pKO2 linear plasmid prophage of Klebsiella oxytoca . J Bacteriol 2004, 186: 1818-1832. 10.1128/JB.186.6.1818-1832.2004

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  37. 37.

    Mobberley JM, Authement RN, Segall AM, Paul JH: The temperate marine phage ΦHAP-1 of Halomonas aquamarina possesses a linear plasmid-like prophage genome. J Virol 2008, 82: 6618-6630. 10.1128/JVI.00140-08

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  38. 38.

    Zhang X, Studier FW: Multiple roles of T7 RNA polymerase and T7 lysozyme during bacteriophage T7 infection. J Mol Biol 2004, 340: 707-730. 10.1016/j.jmb.2004.05.006

    PubMed  CAS  Article  Google Scholar 

  39. 39.

    Kanamaru S, Kondabagil K, Rossmann MG, Rao VB: The functional domains of bacteriophage T4 terminase. J Biol Chem 2004, 279: 40795-40801. 10.1074/jbc.M403647200

    PubMed  CAS  Article  Google Scholar 

  40. 40.

    Bull JJ, Vimr ER, Molineux IJ: A tale of tails: Sialidase is key to success in a model of phage therapy against K1-capsulated Escherichia coli . Virology 2010, 398: 79-86. 10.1016/j.virol.2009.11.040

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  41. 41.

    De Lappe N, Doran G, O’Connor J, O’Hare C, Cormican M: Characterization of bacteriophages used in the Salmonella enterica serovar Enteritidis phage-typing Scheme. J Med Microbiol 2009, 58: 86-93. 10.1099/jmm.0.000034-0

    PubMed  CAS  Article  Google Scholar 

  42. 42.

    Turner D, Hezwani M, Nelson S, Salisbury V, Reynolds D: Characterization of the Salmonella bacteriophage vB_SenS-Ent1. J Gen Virol 2012, 93: 2046-2056. 10.1099/vir.0.043331-0

    PubMed  CAS  Article  Google Scholar 

  43. 43.

    Tiwari BR, Kim S, Kim J: Complete genomic sequence of Salmonella enterica serovar Enteritidis phage SE2. J Virol 2012, 86: 7712. 10.1128/JVI.00999-12

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  44. 44.

    Kang HW, Kim JW, Jung TS, Woo GJ: wksl3, a New biocontrol agent for Salmonella enterica serovars enteritidis and typhimurium in foods: characterization, application, sequence analysis, and oral acute toxicity study. Appl Environ Microbiol 2013, 79: 1956-1968. 10.1128/AEM.02793-12

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  45. 45.

    Kim SH, Park JH, Lee BK, Kwon HJ, Shin JH, Kim J, Kim S: Complete genome sequence of Salmonella bacteriophage SS3e. J Virol 2012, 86: 10253-10254. 10.1128/JVI.01550-12

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  46. 46.

    Salifu SP, Valero-Rello A, Campbell SA, Inglis NF, Scortti M, Foley S, Vázquez-Boland JA: Genome and proteome analysis of phage E3 infecting the soil-borne actinomycete Rhodococcus equi. Environ Microbiol Rep 2013, 5: 170-178. 10.1111/1758-2229.12028

    PubMed  CAS  Article  Google Scholar 

  47. 47.

    Fan H, Fan H, An X, Huang Y, Zhang Z, Mi Z, Tong Y: Complete genome sequence of IME11, a new N4-like bacteriophage. J Virol 2012, 86: 13861. 10.1128/JVI.02684-12

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  48. 48.

    Jeon J, Kim JW, Yong D, Lee K, Chong Y: Complete genome sequence of the podoviral bacteriophage YMC/09/02/B1251 ABA BP, which causes the lysis of an OXA-23-producing carbapenem-resistant Acinetobacter baumannii isolate from a septic patient. J Virol 2012, 86: 12437-12438. 10.1128/JVI.02132-12

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  49. 49.

    Lee YD, Park JH: Complete genome of temperate phage ENT39118 from Cronobacter sakazakii . J Virol 2012, 86: 5400-5401. 10.1128/JVI.00345-12

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  50. 50.

    Sambrook J, Russell DW: Molecular cloning - a laboratory manual, vol 1-3. 3rd edition. New York: Cold Spring Harbor Laboratory Press; 2001.

    Google Scholar 

  51. 51.

    Adams MH: Methods of study of bacterial viruses. In Bacteriophages. New York: Interscience publishers, Inc; 1959:443-457.

    Google Scholar 

  52. 52.

    Lingohr EJ, Frost S, Johnson RP: Determination of bacteriophage genome size by pulsed-field gel electrophoresis. In Bacteriophages: methods and protocols, volume 2: molecular and applied aspects. Edited by: Clokie MRJ, Kropinski AM. Totowa, NJ: Humana Press; 2009:19-26.

    Google Scholar 

  53. 53.

    Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol 2001, 305: 567-580. 10.1006/jmbi.2000.4315

    PubMed  CAS  Article  Google Scholar 

  54. 54.

    Käll L, Krogh A, Sonnhammer EL: A combined transmembrane topology and signal peptide prediction method. J Mol Biol 2004, 338: 1027-1036. 10.1016/j.jmb.2004.03.016

    PubMed  Article  Google Scholar 

  55. 55.

    Hildebrand A, Remmert M, Biegert A, Söding J: Fast and accurate automatic structure prediction with HHpred. Proteins 2009, 77: 128-132. 10.1002/prot.22499

    PubMed  CAS  Article  Google Scholar 

  56. 56.

    Sonnhammer EL, von Heijne G, Krogh A: A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol 1998, 6: 175-182.

    PubMed  CAS  Google Scholar 

  57. 57.

    Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acids Res 2003, 31: 3406-3415. 10.1093/nar/gkg595

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  58. 58.

    Goa F, Zhang CT: Ori-finder: a web-based system for finding oriCs in unannotated bacterial genomes. BMC Bioinformatics 2008, 9: 79. 10.1186/1471-2105-9-79

    Article  Google Scholar 

  59. 59.

    Conant GC, Wolfe KH: GenomeVX: simple web-based creation of editable circular chromosome maps. Bioinformatics 2008, 6: 861.

    Article  Google Scholar 

Download references


AMK was supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada. This funding body did not have any role in the design of the experiments, in the collection, analysis, and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.

Author information



Corresponding author

Correspondence to Andrew M Kropinski.

Additional information

Competing interests

The authors have no competing interests to disclose.

Authors’ contributions

AMK and PJK designed and guided the project; JMD did most of the writing of the paper and the laboratory work and informatic analysis; RAM directed the HPLC nucleoside analysis; HWA did the electron microscopic examination of the phage; YMS preformed the mass spectrometry, all authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Figure S1: HPLC analysis Digested samples were separated using concentration gradients of H2O with 0.1% formic acid (Buffer A) and acetonitrile (Buffer B) which varied with time and are listed in Additional file 3 and in Figure 3. A mixture of the four unmodified nucleosides was used as a standard and the chromatogram produced by their separation was overlaid on the separation achieved with the experimental sample DNA. Peak elution times remained consistent and no new peaks were evident, supporting the conclusion that there is no modification present in the DNA. The small peak visible at 14 minutes was also fluorescent (fluorescence data not shown), and it was concluded to be a result of contamination of the column and unrelated to this experiment, as the concentration of acetonitrile was extremely high which would result in the elution of any contaminants. Absorbance was measured at 254 nm. (TIFF 1 MB)

Additional file 2: Figure S2: Restriction digestion of phage η DNA. Lane L contains Norgen UltraRanger DNA Ladder™, with size markers indicated in kb. The even numbered lanes contained un-ligated η DNA while ligated DNA was used as the substrate in the odd numbered lanes. Lanes 1 and 2 - digested with NdeI, 3 and 4 with BglI, and lanes 5 and 6 with EcoRI. (TIFF 779 KB)

Additional file 3: Table S1: List of CDSs, putative promoters and terminators found in phage η. An E-value threshold of 0.0001 was used for all homology analysis. (DOC 200 KB)

Additional file 4: Table S2: An SDS-PAGE analysis of the purified structural phage proteins on a 12% SDS-PAGE separation gel alongside a PageRuler™ prestained protein ladder (Fermentas) is presented. Visible protein bands were excised from the gel, and labeled as shown in the figure. Subsequently the nature of the protein in the visible band was identified by UPLC-LTQFT -MS/MS analysis. For every detected protein the protein name, the predicted molecular size (Da), the maximum number of unique spectra and sequence coverage (%) is listed. (DOC 868 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Denyes, J.M., Krell, P.J., Manderville, R.A. et al. The genome and proteome of Serratia bacteriophage η which forms unstable lysogens. Virol J 11, 6 (2014).

Download citation


  • Serratia marcescens
  • Phage evolution
  • Genome
  • Proteome
  • Bioinformatics
  • Lysogeny
  • Unstable lysogeny
  • Modified nucleosides
  • Siphoviridae