First genome sequences of Achromobacter phages reveal new members of the N4 family
© Wittmann et al.; licensee BioMed Central Ltd. 2014
Received: 8 November 2013
Accepted: 21 January 2014
Published: 27 January 2014
Skip to main content
© Wittmann et al.; licensee BioMed Central Ltd. 2014
Received: 8 November 2013
Accepted: 21 January 2014
Published: 27 January 2014
Multi-resistant Achromobacter xylosoxidans has been recognized as an emerging pathogen causing nosocomially acquired infections during the last years. Phages as natural opponents could be an alternative to fight such infections. Bacteriophages against this opportunistic pathogen were isolated in a recent study. This study shows a molecular analysis of two podoviruses and reveals first insights into the genomic structure of Achromobacter phages so far.
Growth curve experiments and adsorption kinetics were performed for both phages. Adsorption and propagation in cells were visualized by electron microscopy. Both phage genomes were sequenced with the PacBio RS II system based on single molecule, real-time (SMRT) technology and annotated with several bioinformatic tools. To further elucidate the evolutionary relationships between the phage genomes, a phylogenomic analysis was conducted using the genome Blast Distance Phylogeny approach (GBDP).
In this study, we present the first detailed analysis of genome sequences of two Achromobacter phages so far. Phages JWAlpha and JWDelta were isolated from two different waste water treatment plants in Germany. Both phages belong to the Podoviridae and contain linear, double-stranded DNA with a length of 72329 bp and 73659 bp, respectively. 92 and 89 putative open reading frames were identified for JWAlpha and JWDelta, respectively, by bioinformatic analysis with several tools. The genomes have nearly the same organization and could be divided into different clusters for transcription, replication, host interaction, head and tail structure and lysis. Detailed annotation via protein comparisons with BLASTP revealed strong similarities to N4-like phages.
Analysis of the genomes of Achromobacter phages JWAlpha and JWDelta and comparisons of different gene clusters with other phages revealed that they might be strongly related to other N4-like phages, especially of the Escherichia group. Although all these phages show a highly conserved genomic structure and partially strong similarities at the amino acid level, some differences could be identified. Those differences, e.g. the existence of specific genes for replication or host interaction in some N4-like phages, seem to be interesting targets for further examination of function and specific mechanisms, which might enlighten the mechanism of phage establishment in the host cell after infection.
During the last years, the number of nosocomially acquired infections caused by multi-resistant Achromobacter strains has rapidly increased and become clinically relevant. Achromobacter xylosoxidans is a motile Gram-negative rod  that apart from being widely distributed in natural environments, mainly in soil  or different water sources , has been recognized as an emerging nosocomial pathogen potentially causing different human infections, including endocarditis [4, 5], bacteremia [6, 7], meningitis , ocular infections [9, 10] or urinary tract infections . Recent studies on antibiotic resistances in Achromobacter strains  showed a high number of resistances against antimicrobial substances from different classification groups, such as penicillins, lincosamides or cephalosporines. Other reports on beta-lactamases [13, 14] and resistances against aminoglycosides  in A. xylosoxidans also stress the urgent need to search for new alternatives to fight this opportunistic pathogen. Bacteriophages as natural opponents of bacteria have become an interesting option for treatment of microbial infections in plants, animals and also humans. However, phages against Achromobacter species were nearly unknown, except for some old publications [16, 17]. Recently, the study by Wittmann et al.  showed that bacteriophages against Achromobacter species can be easily found in the environment as a diverse set of different phages against Achromobacter could be isolated from soil and sewage samples. The majority of these bacteriophages belong to the Siphoviridae with long and flexible tails, but two members of the Podoviridae, phages JWAlpha and JWDelta, could be identified. Both are lytic phages that contain double-stranded DNA genomes, which are very similar, but not identical. Here we present the first genomic data on Achromobacter phages so far and further physiological analysis of both N4-like phages. Escherichia phage N4 was known as a genetic orphan for a long time , having no similarities to any other phages in regard to its conserved genomic structure that includes a large virion-encapsulated RNA polymerase which plays an important role in its transcription mechanism with three different RNA polymerases. Due to the increasing number of sequenced phage genomes, we now have numerous N4-like phages for members of the γ-Proteobacteria [19–21] and α-Proteobacteria . In this study we present the first phages of the N4-like genus from Achromobacter, a member of the β-Proteobacteria.
Incubation at different pHs ranging from pH 4-pH 10 for 1 hour and heating at different temperatures (30–60°C) for 10 min did not have any effect on the titer (data not shown) Storage of lysates over months at 4°C did not result in loss of titer, the phages seem to be stable.
First endonuclease restriction digestions of their DNA gave a first hint that they are not identical . JWAlpha and JWDelta both contain linear, double-stranded DNA genomes with a length of 72329 bp and 73659 bp, respectively, comparable to all N4-like phages known so far. Both genomes reveal terminal redundant ends of 365 bp for JWAlpha and 420 bp for JWDelta. After analysis of the genomic DNA by pulsed field gel electrophoresis only one single band, but no concatemers could be identified (data not shown) which suggests that both genomes possess no cohesive ends as also previously described for phage N4 [25, 26]. Additionally, a time-limited BAL 31 digestion and a following hydrolysis with restriction endonuclease SmaI revealed that the genomes are not circularly permutated since two distinct bands as representatives of the ends of the genomes were truncated (data not shown). The G + C contents of JWAlpha and JWDelta are 54.4% and 54.2%, respectively. In comparison with complete Achromobacter genomes sequenced so far (65-66% GC content) [27, 28], the G + C contents of the phage genomes are significantly lower than those of possible hosts, which also occurs in other N4-like phage genomes . Open reading frames (ORFs) were determined using MyRAST  and ARTEMIS based on three different start codons (ATG, GTG and TTG) and putative Shine Dalgarno sequences upstream the coding regions. Thus, we identified 92 ORFs for JWAlpha (88 ATG, 3 GTG and 1 TTG) and 89 ORFs for JWDelta (86 ATG, 2 GTG and 1 TTG), resulting in a coding percentage of 95.9% and 94.1%, respectively. Analysis of the genome sequences using tRNAscan-SE  did not reveal any genes for tRNAs, which leads to the assumption that both phages are adapted ideally to their host organism in regard to codon usage and do not require additional tRNAs of their own. Comparing the number of encoded tRNAs among all N4-like phages, it could be shown that most N4-like phage harbor between 1–3 genes for different additional tRNAs with codons that the host cell does not provide. Whereas Sulfitobacter phage pCB2047-B (Genbank accession number HQ317387, unpublished) and N4-like phages from Salmonella harbor 10 or more genes for different tRNAs, the genomes of the N4-like phages from the Pseudomonas group did not reveal any tRNA genes. Based on a work of de Paepe et al., Bailly-Bechet et al. suggested that virulent phages have more tRNA genes than temperate phages to ensure optimal translation and therefore might replicate faster [31, 32]. In this context, it would be interesting to compare the number of tRNAs with latency times of all N4-like phages.
E. coli phage N4 has been recognized to be a genetic orphan for a long time, being the only known phage that did not depend on the RNA polymerase (RNAP) of its host to transcribe its early genes. N4-like phages that have been described before harbor at least three genes for RNAPs for the transcription of genes in different stages of their life-cycle . The most striking and highly conserved characteristic of a N4-like phage is the gene for a large RNAP with around 3500 amino acids that is encapsulated in the virion and is injected into the host cell together with the phage DNA for immediate start of early gene transcription [35, 36]. Both Achromobacter phages also harbor three different RNAPs, suggesting the same transcription organization as N4 and its other relatives. The large virion RNAPs consist of 3425 amino acids in both phages, only differing in 10 amino acid residues and containing no cysteine residues as also reported by Kazmierczak et al. . Cysteine residues in a protein might form disulfide bonds and lead to a change of protein folding and tertiary structure, which would make it difficult or even impossible for such a large protein to be injected through the phage tail into the host cell. Like other N4-like phages JWAlpha and JWDelta also carry genes that encode two different RNA polymerase subunits, RNAP1 and RNAP2, which play a role in the transcription of N4 middle genes . Similar to phage N4 the gene for RNAP1 is directly followed by RNAP2, whereas other N4-like phage genomes show insertions of small genes in this area [18, 20, 39]. In contrast to other N4-like phages the gene for RNAP2 is interrupted by a gene similar to gp8 of Celeribacter phage P12053L in both Achromobacter phages, dividing it into two parts and raising the question whether it is still functional.
The so-called early genes of a phage genome often encode proteins that somehow interact with the host, for example to protect the newly-injected phage DNA molecule from degradation in the host cytoplasma and see to its establishment, to shut off the host metabolism or to modify host polymerases for transcription in favor of phage propagation. Phages can shut off their host metabolism in different ways, e.g. by host DNA degradation [40, 41]. Lavigne et al.  recently identified gp13 of Pseudomonas phage LUZ19 as an acetyltransferase that hinders both replication and transcription in Pseudomonas posttranslationally by acetylation. In regard to DNA protection, for example, phages use different strategies to avoid degradation by host nucleases. Phage T4 e.g. protects its free genomic ends from degradation by exonucleases with the help of small DNA-binding proteins , phage T7 possesses an anti-restriction system that inhibits the activity of the bacterial restriction endonuclease by prevention of its binding to the phage DNA [44, 45]. Phage JWAlpha also harbors a gene (gp24) that might play a role in a putative anti-restriction system. Its deduced amino acid sequence reveals a conserved domain with homologies to restriction alleviation protein Lar from the defective prophage Rac  which is able to alleviate restriction and to enhance modification by the restriction and modification system in E. coli at the same time. This is apparently accomplished by altering the activity of a type I methyltransferase to efficiently methylate the phage DNA and protecting it against degradation . Though the two phages have a quite similar genome organization, we could not identify this gene in JWDelta. The deduced amino acid sequence of a putative candidate at the same genomic position as in JWAlpha was comparatively much shorter and did not reveal any homologies to proteins of JWAlpha or Lar-like proteins from the database. Whether JWDelta just lacks this gene or whether the gene product of this putative candidate has a similar function, but a yet unknown structure stays to be examined. In general, the different mechanisms of host takeover after infection are worth closer examination in the future.
Another interesting example for a host interaction system could also be found in both Achromobacter phage genomes. Upstream the replication cluster two genes were assigned that share 99% and 98% identities at the nucleotide level with the homologous gene in the other phage, respectively. Their deduced amino acid sequences reveal homologs to a rIIA-like and a rIIB-like protein from phage N4. This kind of proteins were first described in phage T4 and might play a role in lysis inhibition among other functions that are still not completely clear, but they seem to have an influence on cell energetics, including Mg2+ transport and ATP biosynthesis . Examining all N4-like phage genomes known so far for these genes, we discovered that homologs of both genes are widely distributed among N4-like genomes in different host genera. They could also be identified in N4-like phage genomes from the Roseobacter clade, from E. coli and Pseudomonas, but neither in Erwinia phage vB_EamP_S6 nor in Salmonella or Vibrio phages. Generally, both phage genomes are quite similar in their structure and number of genes. The only differences between them could be found in the area between the early genes for RNA polymerases and the genes for rIIa- and rIIB-like proteins. Apart from the already mentioned gene for a Lar-like protein, JWDelta lacks homologous genes for orf29, orf37 and orf41 from JWAlpha. In comparison with JWDelta, the genome of JWAlpha showed two deletions which resulted in a much shorter orf18 with an unclear function and a missing gap of about 1,3 kb that contains a gene for a HNH endonuclease in JWDelta.
Both phage genomes of this study harbor genes for a complete replication cluster in the same order as N4 and other N4-like phages, including a gene for a nucleoside triphosphate pyrophosphohydrolase (Superfamily cl16941) that might provide the helicase with energy by hydrolysis of NTPs, and a gene for a DNA helicase containing a conserved domain UvrD_C_2 (pfam13538) which resembles a AAA-like structural fold that is often found at the C-terminus of helicases. Moreover, this cluster contains a gene that codes for a DNA polymerase with a 3'-5' exonuclease proofreading domain (35EXOc (smart00474)) and a second conserved domain at the C-terminus, DNA_pol_A (pfam00476). The N4 gp43-like gene for a primase was identified based on its similarities to other primase genes from several phages at the amino acid level and the C-terminal Primase C domain PriCT_1 (pfam08708). At the end of the replication cluster we assigned a gene for a putative single-stranded binding protein (Ssb) that though lacking a conserved domain showed a characteristic glycine- and proline-rich C-terminus. The ssb gene is followed by two genes for putative endonucleases, their deduced amino acid sequences revealed characteristic conserved domains NUMOD4 (pfam07463) and HNH_3 (pfam13392) for putative HNH homing endonucleases. Both phages also possess genes involved in the DNA metabolism, to be specific genes for a deoxycytidine triphosphate deaminase with a conserved dcd (PRK00416) domain, a thymidylate synthase (Thy1 (pfam02511)), a homing endonuclease (HNH_3 (pfam13392)) and a thioredoxin with a conserved TRX family domain (cd02947) that might play a role in the regulation of some enzymes involved in DNA metabolism or replication. In contrast to N4-like phages from the Roseobacter clade that are the only other members that also carry genes for thioredoxin, the Achromobacter phages do not harbor a gene for a ribonucleotide-diphosphate reductase, an enzyme involved in formation of deoxyribonucleotides. In general, the genes of the replication cluster are highly conserved and show strong similarities to genes from other N4-like phages at the amino acid level, e.g. the genes for the helicase, DNA polymerase or Ssb. However, besides these genes, there are other genes (orf52, orf55, orf57 and orf58) with a yet unclear function that reveal similarities only to N4-like phages from Escherichia and Erwinia.
Comparison of homologous structural genes of different N4-like phages
Identified proteins of JWAlpha
Homologous genes in Erwinia phage vB_EamP-S6
Homologous genes in Pseudomonas phage LIT1
Homologous genes in Vibrio phage φJA1
vRNAP, 369 kDa
gp79, 380.6 kDa, 28% aa identities
gp71, 369.7 kDa, 23% aa identities
gp53, 335 kDa, 23% aa identities
Structural protein, gp69 66.2 kDa
gp80, 55.9 kDa, no similarities at the aa level
gp72, 52.5 kDa, no similarities at the aa level
gp54, 47.1 kDa, no similarities at the aa level, not identified
Structural protein, gp70, 17.6 kDa
gp81, 16.3 kDa, 40% aa identities
gp73, 16.6 kDa, 38% aa identities
gp55, 15.3 kDa, 31% aa identities, not identified
Major coat, 44.4 kDa
gp85, 43.7 kDa, 66% aa identities
gp77, 44 kDa, 53% aa identities
gp59, 47 kDa, 44% aa identities
Portal, 85.4 kDa
gp88, 81.3 kDa, 60% aa identities
gp80, 81.7 kDa, 48% aa identities
gp62, 79.2 kDa, 47% aa identities
Tail component, gp82, 87.8 kDa
gp65, 87.3 kDa, no similarities at the aa level, not identified
Tail component, gp83, 26.3 kDa
gp97, 26.7 kDa, 36% aa identities, not identified
gp83, 27.9 kDa, 26% aa identities
gp66, 25.8 kDa, 32% aa identities, not identified
Surface protein, gp20, 33.8 kDa
After assembly of new phage progeny, most phages lyse their bacterial host cell in order to release the new phages into the environment. For that purpose, they often use a dual lysis system that consists of a holin that forms a pore in the cytoplasmic membrane to release the endolysin into the periplasmic space where it enzymatically degrades the peptidoglycan . In general, due to the much higher complexity of the Gram-positive cell wall with a lot of different peptidoglycan subtypes , endolysins of phages with Gram-positive hosts possess a lot of highly specific enzymatic activities. However, most of the phages with Gram-negative hosts use unspecific N-acetylmuramidases for degradation of the peptidoglycan as there is just one peptidoglycan layer to overcome. The lysis clusters of both Achromobacter phages consist of three genes (orf78-orf80 in JWAlpha, orf75-orf77 in JWDelta) and show the same organization as lysis clusters from other phages with Gram-negative hosts. Analysis with TMHMM Server v. 2.0  revealed that the deduced amino acid sequences of orf80 from phage JWAlpha and orf77 from JWDelta both contain two transmembrane regions and show similarities to putative holins of Escherichia phages EC1-UPM and vB_EcoP_G7C, hence both proteins could be assigned as class II holins . The gene for the holin is followed by the endolysin gene, BLASTP analysis of the amino acid sequences identified them both as N-acetylmuramidases. Downstream the endolysin genes, one small ORF for a protein could be assigned that also revealed one transmembrane region and similarities to a Rz-like protein. Besides genes for a holin and an enzymatically active endolysin, phages of gram-negative hosts often harbor genes for proteins Rz and Rz1 that somehow interact with the outer membrane and play a role in the final step of host lysis [54, 55]. Examining all N4-like genomes in regard to their lysis clusters, we could identify five different groups of clusters with apparently genus-specific strategies for host lysis. The largest group of the N4-like phages harbors genes for a class II holin with two transmembrane regions, a N-acetylmuramidase and a Rz protein with one transmembrane region as described above for both Achromobacter phages. Further members of this group are the Escherichia phages N4, EC1-UPM, KBNP21 and vB_EcoP_G7C, Erwinia phage vB_EamP_S6 and Enterobacter phage IME11. The second group comprises the N4-like phages from the Roseobacter clade and Enterobacter phage ECP1 which all encode genes for a lysis protein. Since none of their deduced amino acid sequences revealed homologies to muramidases or to peptidases of the MEROPS database, they were apparently annotated as lysis proteins with yet unclear functions. In all roseophage genomes of this group, this gene is framed by two small orfs with one deduced transmembrane region. Though the Enterobacter phage ECP1 also possesses a gene that encodes a lysis protein, the structure of its lysis cluster is different. The whole cluster is located on the other DNA strand and is not transcribed towards the same direction as the preceding and following genes. In addition to that the gene for the lysis protein is preceded by two genes for proteins with three transmembrane regions each, which is characteristic for class I holins. The N4-like Salmonella phages FSL SP-058 and FSL SP-076 belong to the third group with a different lysis cluster. Like phage ECP1, both Salmonella phages possess two putative genes for class I holins, but in addition to a lysis protein they also harbors a gene for a putative lysozyme and, thus, they are the only N4-like phages with two different enzymes for host lysis. The fourth group consists of three Pseudomonas phages of the N4 family, namely P26, LIT1 and LUZ7. Though two putative holin genes could be identified in all three genomes, no muralytic enzyme could be predicted. The endolysin containing region seems to be completely deleted  which raises the question whether there are other muralytic enzymes with so far unknown mechanisms. Genomes of members of the fifth group, mainly Vibrio phages, even do not reveal any sign of lysis cluster at all which makes them interesting candidates for examining their lysis mechanism and potentially identifying new muralytic enzymes. Since all bacteriophages of the N4-family known so far have gram-negative hosts, it was no surprise that nearly all of them encode genes for rather unspecific muramidases and no specific endopeptidases as often found in phages with gram-positive hosts . Although the genomes of N4-like phages are highly conserved both in sequence identity and structure, it is interesting to observe these differences in their lysis clusters. In contrast to the replication cluster or the cluster for structural genes, a broader diversity could be determined which is worth further examination.
Analysis of the genomes of Achromobacter phages JWAlpha and JWDelta and comparisons of different gene clusters with other phages revealed strong similarities to other N4-like phages, especially of the Escherichia group. Phages of the N4 family infect a lot of different genera with Achromobacter as a new member. Although all these phages show a highly conserved genomic structure and partially strong similarities at the amino acid level, interestingly, some differences could be identified probably depending on their host genus. These differences, e.g. the existence of a gene for a ribonucleotide-diphosphate reductase in N4-like Roseobacter phages or the completely different structure of the replication cluster in N4-like phages of Pseudomonas compared to all other N4-like phages, seem to be interesting targets for further examination of function and specific mechanisms. In particular, further examination of the gene for a Lar-like protein in JWAlpha might enlighten the mechanism of phage establishment in the host cell after infection. A closer look on the mechanisms for host lysis in those phages that revealed no or yet unknown enzymes for cell wall degradation might also be interesting for the identification of new muralytic enzymes.
Genome sequences of JWDelta and JWAlpha were deposited to Genbank under accession numbers KF787094 and KF787095, respectively. Further genome sequences of N4-like phages for comparison could be obtained under following accession numbers: Enterobacteria phage N4 (EF056009), Escherichia phage vB_EcoP_G7C (HQ259105) , Escherichia phage KBNP21 (JX415535) , Enterobacter phage EcP1 (HQ641380), Enterobacter phage IME11 (JX880034) , Escherichia phage EC1-UPM (KC206276) , Erwinia phage vB_EamP-S6 (HQ728266) , Pseudoalteromonas phage pYD6-A, Pseudomonas phage LIT1 (FN422399), Pseudomonas phage LUZ7 (FN422398) , Pseudomonas phage PA26 (JX194238) , Roseovarius sp. 217 phage 1 (FR682616), Roseovarius Plymouth Podovirus 1 (FR719956), Ruegeria phage DSS3φ2 (FJ591093), Sulfitobacter phage EE36φ1 (FJ591094) , Sulfitobacter phage pCB2047-B (HQ317387), Salmonella phage FSL SP-058 (KC139517), Salmonella phage FSL SP-076 (KC139520) , Vibrio phage VBP32 (HQ634196), Vibrio phage VBP47 (HQ634194), Vibrio phage JA-1 (KC438282), Vibrio phage VCO139 (KC438283) .
Achromobacter xylosoxidans strain DSM 11852 was used for phage propagation and adsorption tests. It was cultured in liquid tryptic soy broth TSB (30 g-l, pH 7.5, OXOID) medium or on Petri dishes of TSB medium supplemented with 1.5% agar (w/v) for approximately 16 h (overnight) at 28°C.
Purification of the phages by CsCl gradient centrifugation and of their DNA by phenol extraction was done as previously described .
SMRTbell™ template libraries were prepared according to the instructions from Pacific Biosciences, Menlo Park, CA, USA, following the Procedure & Checklist for “Low-Input 10 kb Template Preparation and Sequencing” using C2-C2 Chemistry.
Briefly, for preparation of 10 kb libraries ~15 μg genomic DNA was sheared in an Eppendorf microcentrifuge 5424 2×3 min at 5500 rpm using g-tubes™ from Covaris, Woburn, MA, USA. Size range was monitored on an Agilent 2100 Bioanalyzer from Agilent Technologies, Santa Clara, USA. DNAs were concentrated, end-repaired and ligated to hairpin adapters applying components from the DNA/Polymerase Binding Kit 2.0 from Pacific Biosciences, Menlo Park, CA, USA. Reactions were carried out according to the manufacturer´s instructions with exception of the ligation, which took place 2 h at 25°C and overnight at 4°C. SMRTbell™ templates were exonuclease treated for removal of incompletely formed reaction products. A mixture of Exonuclease III and Exonuclease VII (Affymetrix, High Wycombe, UK) was utilized. Hereby, an additional exonuclease treatment after the second purification took place, followed by three purification steps as recommended by Pacific Biosciences. Conditions for annealing of sequencing primers and binding of polymerase to purified SMRTbell™ templates were assessed with the Calculator in RS Remote, PacificBiosciences, Menlo Park, CA, USA. SMRT sequencing was carried out on the PacBio RS (PacificBiosciences, Menlo Park, CA, USA) using the DNA Sequencing Reagent 2.0. For each 10 kb libraries one 120 minutes movie was taken.
Phage genome assemblies were performed using the “RS_Preassembler_Allora.1” protocol included in SMRTPortal version 1.3.3. Within that following parameters were applied. For Filtering: Minimum Readlength 100, Minimum Subreadlength 500, Minimum Read Quality 0.8. For Assembly: Minimum Seed Read Length: 3500 for phage JWAlpha, BLASR Options (Advanced) minReadLength 200 maxScore 1000 bestn 24 maxLCPLength 16 nCandidates 24, Trim FASTQ Output true, Use CCS false, Overlap Permissiveness Least permissive, Expected Genome Size 50000, Maximum iterations 10, Detect Chimeras true. In both cases one final contig could be obtained. A quality check of that given consensus sequences regarding overall coverage as well as SNPs was performed using SMRT View and IGV  after application of the “RS_Resequencing.1” protocol with default parameters. The phage genomes were annotated using RAST  with subsequent manual curation in Artemis .
For adsorption assays, exponentially growing cultures of Achromobacter xylosoxidans strain DSM 11852 were infected with phages JWAlpha and JWDelta, respectively, with a multiplicity of infection (moi) of 0.01, for growth curve experiments, cells were infected with a moi of 0.1. Samples for both experiments were taken at different time points and directly filtered (membrane syringe filter 0.45 μm, Sartorius, order number 16555). Filtrates were diluted and spotted on a bacterial lawn of DSM 11852 to determine titer of free phages in the culture.
Treatment of phage samples for SDS-PAGE analysis and SDS-PAGE analysis itself were performed as previously described .
Coomassie stained protein bands were excised from the gel. The slices were washed twice in H2O for 5 min and dehydrated by adding 200 μl 50% acetonitrile for 5 min. Gel pieces were then washed in 200 μl 0.1 M NH4HCO3 for 15 min at room temperature. Washing solution was discarded and slices were dehydrated by adding acetonitrile again. Dehydrated gel pieces were completely dried in Speed Vac for 10 min. For the tryptic digestion of the proteins the gel slices were reswelled by adding 50 μl digestion solution containing trypsin and incubated over night at 37°C. For extraction of peptides 50 μl of acetonitrile were added, followed by incubation under vigorous shaking for 15 min at 37°C. Supernatants were removed and gel slices were treated with 50 μl 5% formic acid under vigorous shaking for 15 min at 37°C. After that 50 μl of acetonitrile were added, followed again by incubation under vigorous shaking for 15 min at 37°C. Supernatants were united with the removed supernatants and concentrated via Speed Vac. 40 μl 32% methanol/0.25% HCOOH were added, followed by 3 min treatment with ultrasound and analyzed by MALDI Ultraflex-TOF/TOF. Results were compared with a database of predicted phage proteins using Mascot. MALDI-TOF analysis was done at the Helmholtz Centre for Infection Research, Braunschweig, Germany.
Thin carbon support films were prepared by sublimation of a carbon thread onto a freshly cleaved mica surface. Phages were negatively stained with 0.5% (w/v) aqueous uranyl acetate, pH 5.0, according to the method of Valentine et al. . Samples were examined in a TEM 910 transmission electron microscope (Carl Zeiss, Oberkochen) at an acceleration voltage of 80 kV.
Bacteria incubated with phage JWAlpha were fixed with 2% glutaraldehyde and 5% formaldehyde, washed with TE-buffer (20 mM TRIS, 1 mM EDTA, pH 6.9) then dehydrated in a graded series of acetone (10, 30, 50, 70, 90, 100%) on ice for 15 min for each step, critical-point dried with liquid CO2 (CPD 30, Bal-Tec, Liechtenstein) and covered with a gold-palladium film by sputter coating (SCD 500, Bal-Tec, Liechtenstein) before being examinated in a field emission scanning electron microscope (Zeiss Merlin) using the Everhart Thornley SE detector and the inlens detector in a 25:75 ratio at an acceleration voltage of 5 kV.
Samples were fixed in medium with 5% formaldehyde after 1 min 2% glutaraldehyde was added and kept on ice for 1 hour, washed with cacodylate buffer, and further fixed with 1% osmium in 0.08 M cacodylate buffer for 1 hour at room temperature. Then, samples were dehydrated with a graded series of acetone (10, 30, 50%) on ice. At the 70% dehydration step samples were incubated overnight in 70% acetone containing 2% uranyl acetate at 4°C. On the next day samples were further dehydrated with 90% and 100% acetone on ice, the 100% aceton step was repeated at room temperature. Samples were embedded in the epoxy resin Spurr according to described procedures . Ultrathin sections were cut with a diamond knife, picked up with butvar-coated grids, counterstained with uranyl acetate (3 min) and lead citrate (15 sec), and examined in a TEM910 transmission electron microscope (Carl Zeiss, Oberkochen) at an acceleration voltage of 80 kV. Images were recorded digitally at calibrated magnifications with a Slow-Scan CCD-Camera (ProScan, 1024×1024, Scheuring, Germany) with ITEM-Software (Olympus Soft Imaging Solutions, Münster, Germany). Contrast and brightness were adjusted with Adobe Photoshop CS3.
The phage dataset comprised of 24 genomes which were compared pairwise on the proteome level using the Genome Blast Distance Phylogeny approach (GBDP) [57–59]. During distance calculation GBDP considers and corrects for several factors that could potentially bias the results such as overlapping HSPs. In general, GBDP was shown to yield robust results even in the presence of a significant amount of paralogous genes, large repeats and partially incomplete genomes. Further, GBDP does not use arbitrary selections of marker genes but rather detects all available matches (i.e., high-scoring segment pairs, HSPs) between pairs of genomes (e.g., via BLAST + ). The resulting set of pairwise distances was subsequently used to infer a phylogenetic tree under the balanced minimum evolution criterion  as shown previously . If distance values could not be obtained due to the lack of HSPs between the respective proteomes, twice the maximum distance observed in the matrix was used as a substitute (as in ). Since GBDP is capable of calculating a set of bootstrap replicates for each pairwise distance , the inferred tree was augmented with branch support values to assess its statistical significance.
We would like to thank Simone Severitt for preparing the library and sequencing and Anja Meier (Helmholtz Centre for Infection Research) for performing the peptide mass fingerprinting analysis. Cordial thanks are addressed to Markus Göker (DSMZ) for helpful comments.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.