Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2008, p. 7552-7560, Vol. 74, No. 24
0099-2240/08/$08.00+0 doi:10.1128/AEM.01709-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Área de Microbiología, Facultad de Medicina, Universidad de Oviedo, Oviedo, Spain,1 Instituto de Productos Lácteos de Asturias (CSIC), Villaviciosa, Spain,2 Centro Comunitario de Sangre y Tejidos and CIBER de Enfermedades Raras (CIBERER), Oviedo, Spain3
Received 24 July 2008/ Accepted 15 October 2008
|
|
|---|
|
|
|---|
The use of bacterial viruses as therapeutic agents was initiated in 1919, 3 years after their discovery, for the treatment of bacillary dysentery and continued up to the decade of the 1940s for the treatment of a variety of infectious diseases such as bubonic plague, cholera, and staphylococcal septicemia (9, 10, 23). These practices were interrupted by the advent of antibiotics, whose wider spectrum of action than that of phages opened the possibility to apply empirical treatments. However, the selection of antibiotic-resistant bacteria and the frequent occurrence of adverse effects associated with their use led to renewed interest in phages as therapeutic agents (5, 20, 27, 28, 33).
In addition, phages are very useful tools both for understanding the biology of their bacterial hosts, especially the regulation of gene expression, and for developing vectors that make their hosts amenable to recombinant DNA manipulation (18, 29).
Three phages active on S. maltophilia have been described previously in some detail. Two of them (
SMA5 and Smp14) are myoviruses with wide host ranges and genomes of about 250 and 160 kb that appear to be resistant to digestion by most restriction endonucleases (3, 4). The other one (
SMA9) has a 6.9-kb genome that harbors seven genes, the largest being an orthologue of the gene encoding the toxin Zot of the filamentous
Lf bacteriophage of Vibrio cholerae (14).
In this communication, the isolation and characterization of several Stenotrophomonas phages are described. In addition, the sequencing of the genome of the first known temperate phage infecting S. maltophilia is reported.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Stenotrophomonas strains used in this work and host ranges of phages S1, S3, and S4
|
Bacteriophage enrichment and isolation from sewage.
Six samples of water, obtained from the aeration and decantation tanks of local sewage plants, were incorporated into preparations of 2x TY liquid medium, which were used (unsterilized) to set up two successive enrichment cultures by inoculation with a mix of the 26 Stenotrophomonas strains (1% of an overnight culture of each strain). After overnight incubation at 30°C with shaking, the cultures were centrifuged and the supernatants were filtered through 0.22-µm-pore-diameter filters and used as a source of phages in a new enrichment step with the same mix of strains (in sterile medium). Supernatants of these second cultures were filtered and tested for the presence of inhibition halos by placing drops onto lawns of each of the 26 bacterial strains. Phages from the lysis zones were suspended in SM buffer (22) and diluted to get single plaques. The procedure was repeated twice, and the purified phages were stored at –70°C in SM buffer with 50% glycerol.
Phage isolation from lysogenic strains.
Exponential cultures of each Stenotrophomonas strain were treated with mitomycin C or with ciprofloxacin at the MIC for the strain and half the MIC (the MICs of mitomycin C ranged from 4 to 16 µg/ml, and those of ciprofloxacin ranged from 4 to 8 µg/ml). Supernatant aliquots of these cultures (5 µl) were placed onto lawns of each of the 26 bacterial strains, and after incubation, the lawns were inspected for the appearance of inhibition halos. Phages were isolated from the lysis zones by following the procedure described above.
Host range, phage purification, and functional studies.
Phage host ranges were determined by placing 5-µl drops of each viral suspension onto lawns of each of the 26 Stenotrophomonas strains. Large-scale phage propagation and purification through CsCl gradients were performed as described in reference 30.
Curves from one-step analyses of growth in 2x TY at a multiplicity of infection of 0.03 were generated as described in reference 15.
The frequency of emergence of bacteriophage-insensitive mutants (BIMs) was determined after the infection of exponential cultures at a multiplicity of infection of 10. The BIM frequency was calculated by dividing the number of surviving colonies by the bacterial titer of a parallel, uninfected culture.
The experiments described above were performed at least in triplicate.
Electron microscopy.
Phage gradient drops were placed onto 0.022-µm-pore-diameter filters floating on SM buffer at room temperature for 30 min (with two changes of buffer) to eliminate CsCl. Aliquots were negatively stained with 2% (wt/vol) uranyl acetate on copper grids covered with a film of carbon and ionized. Grids were observed with a JEOL 1200 EXII microscope stabilized at 100 KV. Images were obtained on Scientia photographic plates (AGFA).
DNA and protein techniques.
Phage DNA extraction from CsCl-purified virion suspensions dialyzed against SM buffer was performed as described in reference 12.
Pulsed-field gel electrophoresis was used to determine the sizes of phage genomes. Phage DNA was loaded into a 1% agarose gel in 0.5x Tris-borate-EDTA buffer (pH 8) and resolved using a CHEF-DR III pulsed-field system (Bio-Rad Laboratories, Hercules, CA) at 6 V/cm for 20 h and 14°C. Linear ramped pulse times of 3 to 63.8 s were employed.
Enzymatic treatments of phage DNA were carried out by following the instructions provided by the suppliers of the enzymes. DNA was visualized after separation in 0.7% agarose gels in Tris-borate-EDTA buffer and staining with ethidium bromide under UV light.
Genotypic identification of the 26 bacterial strains used in this work was performed by multiplex PCR using primers specific for S. maltophilia, the Burkholderia cepacia complex, and Pseudomonas aeruginosa, plus a pair targeting a universal 16S RNA stretch as a positive control (6). The reactions were set up with 25-µl aliquots using 2.5-µl samples of exponential cultures without any processing as DNA sources.
Phage structural proteins were extracted and purified as described in reference 12 and analyzed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) with a Mini-PROTEAN III system (Bio-Rad). Proteins were stained with Coomassie R-250 blue. The bands of interest were excised manually and digested with porcine trypsin (Promega), and the resulting peptides were analyzed by matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass spectrometry, essentially as described previously (13).
Phage genome sequencing and analysis.
Approximately 10 µg of S1 phage DNA was sheared by sonication, and segments were selected on the basis of size (2 to 3 kbp) and cloned into pUC18. Individual clones were sequenced and assembled. Trace assembly was done with the phredPhrap package (11). Generally, the assembled phage genome sequences had at least fivefold coverage in both orientations. Primer walking was sometimes used to close gaps in the sequence.
The cos sites were localized through analysis of the restriction patterns of ligated and unligated S1 DNA samples. To obtain the sequence of the single-stranded cos extensions, the ends of S1 DNA were sequenced directly using primers that pointed outwards. The nucleotide sequences obtained were aligned with that of the circular S1 genome. They overlapped by 12 bases that were considered to constitute the cosN site.
Open reading frames (ORFs) were predicted by the Clone Manager 7 and the ORF Finder (http://ncbi.nlm.nih.gov/gorf.html) software and by visual inspection. The primary nucleotide sequence in all reading frames was scanned for the start codons AUG and GUG, with a threshold of 50 codons. BLASTX and BLASTP (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi) were used to search for proteins homologous to the ORF products. Structural predictions and motif searches were performed with InterProScan (http://www.ebi.ac.uk/InterProScan/), Pfam (http://pfam.sanger.ac.uk/search?tab=searchSequenceBlock), YASPIN (http://zeus.cs.vu.nl/programs/yaspinwww/), and TMHMM Server version 2.0 (http://www.cbs.dtu.dk/services/TMHMM-2.0/). Likely
70 target promoter sequences were identified using BPROM (http://www.softberry.com). Putative terminator sequences were detected with the Terminator function of GCG (version 10.2). The search for putative tRNA-encoding genes was performed with tRNAscan-SE 1.21 (http://selab.janelia.org/tRNAscan-SE/). The prediction of +1 and –1 frameshifting was carried out with FSFinder 2 (http://wilab.inha.ac.kr/fsfinder2/).
Cumulative GC skew analysis with a system from In Silico Biology, Inc. (Yokohama, Japan) was used to locate the origin of replication.
Nucleotide sequence accession number.
The sequence of the phage S1 genome has been deposited in GenBank under accession number EU849489.
|
|
|---|
Two procedures were undertaken for phage isolation: (i) the induction of resident prophages and (ii) the enrichment of sewage samples. In the first procedure, presumptive prophages were induced with mitomycin C or with ciprofloxacin and allowed to form plaques on lawns of each of the 26 bacterial strains. As a result, two phages, named S1 and S2, were isolated from S. maltophilia CECT 4793 and E539, respectively.
Complementarily, water samples taken from a sewage plant were used as potential sources of phages. By this procedure, two successive rounds of enrichment with a mix of all 26 Stenotrophomonas strains were applied prior to the placement of drops of the resulting cultures onto lawns of each strain. This procedure yielded 20 phages, which were named S3 to S22.
Of the two temperate phages isolated after the induction of lysogenic cultures, S1 inhibited the growth of four strains while S2 could be propagated only on its original indicator strain. A narrow host range was also the rule among the phages obtained from sewage, with two exceptions: S3 inhibited the growth of 12 strains, while S4 acted on 18 strains. The host ranges of these phages appeared to be confined to Stenotrophomonas, in that none were able to affect any of a collection of 14 Pseudomonas aeruginosa and 5 Burkholderia cepacia strains isolated from clinical samples (data not shown). The phages S1, S3, and S4 were chosen for further study because they had the widest host ranges within their categories.
Structure of the virions.
S1, S3, and S4 were propagated and purified through two successive CsCl gradients and observed under the electron microscope. All showed an icosahedral capsid and a tail, which allows their inclusion in the order Caudovirales (Fig. 1). S1 and S4 belong to the family Siphoviridae (they present a noncontractile tail), while S3 is a member of the family Myoviridae due to its contractile tail. The diameters of the capsids (means ± standard deviations) of S1, S3, and S4 were, respectively, 61.4 ± 1.35, 82.92 ± 1.25, and 87.5 ± 1.5 nm. The tail of S1, with a length of 129.2 ± 1.3 nm and a width of 9.93 ± 0.66 nm, is remarkable for its extreme flexibility, compared with the tail of S4 (201.87 ± 1.22 nm long and 10.7 ± 0.24 nm wide). The tail of S3 is 136.9 ± 0.59 nm by 21.43 ± 0.0 nm and has a basal plate with protruding spikes.
![]() View larger version (143K): [in a new window] |
FIG. 1. Transmission electron micrographs of the bacteriophage S1, S3, and S4 virions, negatively stained with 2% uranyl acetate.
|
![]() View larger version (46K): [in a new window] |
FIG. 2. SDS-PAGE analysis of the S1 virion proteins. Those identified by mass spectrometry are indicated on the left. The unlabeled bands contained mixes of two or more structural proteins. The migration of the molecular mass markers is indicated on the right. The gel track was heavily overloaded to make evident the relationship between the tape measure and tail fiber proteins.
|
Functional properties of the phages.
Phage S3 showed a much faster development cycle than S1 and S4 (Fig. 3). While its eclipse phase lasted about 30 min, those of the other two phages took 75 min. Its rise period was also shorter, 30 min versus 45 min for S1 and 60 min for S4. The burst sizes were calculated to be about 75 virions/cell for S1, 100 virions/cell for S3, and 80 virions/cell for S4.
![]() View larger version (10K): [in a new window] |
FIG. 3. One-step growth curves for phages S1, S3, and S4 on exponential cultures of different Stenotrophomonas strains (Table 1) incubated in 2x TY at 30°C under agitation. The curves are representative of the data obtained from several experimental replicas.
|
Finally, the appearance of BIMs was tested. The frequency was the lowest for S3, with levels that ranged between 10–5 and 10–6 depending on the host, while in the case of S1, these values fell to levels on the order of 10–3 to 10–4. S4 BIMs emerged at a frequency of approximately 10–4.
The long one-step curve for S4, together with the relatively high frequency with which insensitive mutants were generated, suggested that S4 might be a temperate phage. To test this possibility, colonies of presumptively insensitive mutants were isolated, washed with a chelating agent (10 mM sodium pyrophosphate at 37°C for 30 min) to eliminate free virions (30), and plated onto lawns of a strain susceptible to S4. Colonies that generated an inhibition zone were then subjected to phage induction with mitomycin C, thus confirming the temperate nature of S4. It is then assumed that the high frequency of BIMs observed after exposure to S1 and S4 was due to host lysogenization.
Analysis of the bacteriophage S1 genome.
The genome of S1 was determined to be located in a double-stranded DNA molecule 40,287 bp long (Fig. 4). It has 12-bp single-stranded 5'-end extensions with the sequence 5'-GGTATTCGGACG-3'. The G+C content is 63.7%, slightly lower than the 66% of the two sequenced S. maltophilia genomes (GenBank accession numbers NC011071 and AM743169).
|
View larger version (6K): [in a new window] |
FIG. 4. Physical and genetic map of the bacteriophage S1. The 48 orf genes are correlatively numbered and are indicated by arrows proportional to their length and pointing in the direction of transcription. The functional modules detected are indicated above the scheme, and the products of several putatively or experimentally identified genes are shown below. PAPS, phosphoadenosine phosphosulfate.
|
|
View this table: [in a new window] |
TABLE 2. Characteristics of bacteriophage S1 ORFs and their deduced gene productsa
|
-helices, matches scaffold proteins; and its carboxyl terminus matches major head proteins. MALDI-TOF analysis of the most abundant virion protein yielded peptides that matched gp8 from its 360th residue onwards, possibly indicating that native gp8 became autoprocessed by its protease domain, giving rise to the scaffold and the major capsid polypeptides. The latter would have about 300 amino acids and a mass of about 33 to 35 kDa, in good agreement with the apparent size observed by SDS-PAGE. Unfortunately, the NH2-terminal end of this protein was blocked, so we are at present unable to precisely determine where processing takes place. The lysis cassette would comprise orf11 to orf13, apparently corresponding to the holin, the tandem Rz-Rz1, and the endolysin. The putative holin was unusual since only one transmembrane domain was detected in it. However, its location and size suggest that it may behave as a bona fide holin. The next orf encodes both Rz and Rz1. The Rz protein has a typical signal peptide at its amino terminus, while Rz1 is translated in the +1 frame with respect to Rz and gives rise to a polypeptide rich in proline (accounting for 17 of 106 residues). Finally, gp13 has an endolysin domain and a weak sequence relationship to peptidoglycan binding proteins, which together with its location next to a holin, allows the endolysin function to be ascribed to it.
The tail module appears to extend from orf14 to orf21. gp14 is a minor component of the virion, as are gp15 (major tail protein), gp17 (tape measure), and gp19 (tail fiber) (Fig. 2). As expected for tail measure proteins, the secondary structure of gp17 is almost exclusively
-helix. The protein may be processed during particle formation, since it migrates as a 90-kDa polypeptide, whereas the gene sequence indicates a size of 131.7 kDa. In fact, MALDI-TOF analysis provided only peptides that corresponded to the first 924 amino acids of the protein (orf17 has a coding capacity for a polypeptide of 1,248 residues), suggesting that processing occurs at its C-terminal end. Between the major tail and the tape measure genes lies orf16, which does not have the translational frameshift site frequently encountered in this position in other double-stranded DNA phages (37).
The other two modules that are commonly present in temperate siphoviruses, i.e., the lysis-lysogeny switch and the replication cluster, are not easily recognizable in S1. The only protein showing homology with phage transcription regulation proteins is gp29, whose determinant, in addition, is followed by the putative phage integrase gene. However, gp29 has a motif homologous to the positive transcriptional regulator AlpA (35) rather than to the phage repressor CI, it does not have a clear helix-turn-helix motif, and no cro-like repressor gene is found adjacent to orf29. The product of orf28 may be the phage integrase because it presents the C-terminal catalytic domain of DNA breaking-rejoining enzymes and, more particularly, of Cre integrases.
S1 does not have a replication module, but some genes that may be involved in DNA metabolism are scattered along the genome. To this category belong orf1 and orf23, as well as the putative origin of replication. This cis-acting element was identified between orf44 and orf45, in a 214-bp noncoding region that presents direct and inverted repeats typical of replication origins, through cumulative GC skew analysis. gp1 is predicted to be a topoisomerase-primase that has also been found in phages of Vibrio and Burkholderia species (24, 31). Finally, orf23 encodes a putative N6-adenine methyltransferase.
|
|
|---|
The genomic DNA of S3 and S4 is resistant to digestion by restriction endonucleases. This property may contribute to the wide host ranges of both phages. In fact, the other two Stenotrophomonas phages of the order Caudovirales isolated to date are extremely resistant to DNA restriction as well and infect a wide range of susceptible strains: 42.5 and 70% (of 87 strains) for Smp14 and
SMA5, respectively (3, 4). In contrast, S1, whose DNA is readily digested, is able to infect only 4 of 26 strains (15%). Smp14 DNA presents, in addition to the normal four deoxynucleotides, small proportions of N4-deoxymethylcytidine and two other unusual bases, suggested to be deoxythymidine and deoxyguanosine analogues, which are claimed to be responsible for its resistance to digestion (4). This may be true for S3 and S4 as well. In any case, the genetic determinants that render the S3 and S4 DNA refractory to restriction are probably located in the phages' own genomes rather than being a trait conferred by particular hosts, because the phages remained insensitive to digestion after development on strains whose DNA was a substrate for restriction enzymes.
S3 is proposed as a model phage to be used in therapeutically aimed experiments, not only because of this resistance to restriction, but also because it seems to be virulent (although this property has to be confirmed by genome sequencing), which prevents the risk of lysogen generation; it attacks a wide range of clinical strains; it has a fast developmental cycle; its plaquing frequency is independent of the strain used as a host; and resistant mutants (BIMs) are infrequent.
S1, S2, and S4 are the only temperate phages isolated from Stenotrophomonas to date. Of the two that infected more than one strain, S1 is more amenable to study because its genomic DNA is relatively short compared to that of S4 (40 versus 200 kb) and is susceptible to restriction endonucleases: S1 was therefore chosen as a representative of its group, and its genome sequence was determined. The modular organization observed in most double-stranded DNA bacteriophages occurs only in the left half of the S1 genome, which encodes the morphogenetic functions (DNA packaging, head morphogenesis, cell lysis, and tail morphogenesis). However, no lysis-lysogeny switch or DNA replication modules were obvious, so only informed guesses on the possible implication of a few genes in these functions can be formulated, the tasks of the remaining ORFs being unknown due to the lack of homology of their products to sequences from the databases.
The transcription pattern of S1 presents two groups of ORFs running divergently and only two of them in each row going in the opposite direction. Among them, orf31 may represent a lysogenic conversion gene (moron) because its product resembles GspM, a protein involved in the general secretion system II, frequently used by virulence factors to reach the surrounding medium (21). However, no clear –10 and –35 sequences are found in front of it. Adjacent to orf31 lies the determinant for a homologue of phosphoadenosine phosphosulfate reductase, which has been found previously in other phages that infect gram-negative bacteria, such as
E125 and BcepB1A, which infect Burkholderia (32, 36), and coliphage 186 (26). This protein may confer an advantage on the host in the assimilation of inorganic phosphate through its reduction to organic thiols (31).
The opposite direction of transcription of orf43 and orf44 with respect to that of the surrounding genes may also be considered to be an indication of a moron; however, no function could be ascribed to their products apart from a hypothetical interaction with DNA, suggested by the helix-turn-helix motif of orf44.
The case of orf27 is peculiar; it is located immediately behind the integrase gene, with a gap of only 4 nucleotides. This position, together with the considerable degree of secondary structure implied in its sequence (Fig. 5), suggests that the attP sequence may be included in this gene. On the other hand, the direction of transcription of orf27 is opposite that of the surrounding genes, and overall, orf27 has 75% identity to a hypothetical-protein-encoding orf located in the chromosome of S. maltophilia (Fig. 5). This finding may suggest that the integrity of the cellular gene would be restored after the integration of the phage genome into the bacterial chromosome, thus keeping whatever function was played by its gene product.
![]() View larger version (36K): [in a new window] |
FIG. 5. Comparison of the nucleotide sequences of orf27 and a particular Stenotrophomonas chromosomal orf. The zones of more than 60 nucleotides showing a degree of identity of 85% or higher are shaded. Zones of putative secondary structure are indicated by arrows above the sequences (IR, inverted repeat; DR, direct repeat). Numbers at the ends of the sequences are nucleotide positions.
|
The genome of S1 is a new example of the extraordinary biodiversity shown by bacteriophages: we cannot even guess what may be the roles played by the products of almost half of its genes, even though we know that they must include the determinants for such a widespread and widely studied process as the lysis-lysogeny decision.
The provision of some of the Stenotrophomonas strains used in this work by Jose L. Martinez is gratefully acknowledged. We also thank K. F. Chater for critical reading of the manuscript.
Published ahead of print on 24 October 2008. ![]()
|
|
|---|
SMA5. Appl. Environ. Microbiol. 71:1387-1393.
393-A2, a bacteriophage that infects Lactobacillus casei. Microbiology 140:2585-2590.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»