Previous Article | Next Article ![]()
Journal of Bacteriology, March 2005, p. 1966-1973, Vol. 187, No. 6
0021-9193/05/$08.00+0 doi:10.1128/JB.187.6.1966-1973.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
,
Nadine Fornelos,
and
Jacques Mahillon*
Université Catholique de Louvain, Louvain-la-Neuve, Belgium
Received 5 October 2004/ Accepted 8 December 2004
|
|
|---|
NS11) hosts. Several reports have shown that tectiviruses infecting gram-negative bacteria are closely related, whereas no information is currently available on the genetic relationship among those infecting gram-positive bacteria. The present study reports the sequence of GIL16, a new isolate originating from Bacillus thuringiensis, and a genetic comparison of this isolate with the tectiviral bacteriophages Bam35 and GIL01, which originated from B. thuringiensis serovars Alesti and Israelensis, respectively. In contrast to PRD1 and its relatives, these are temperate bacteriophages existing as autonomous linear prophages within the host cell. Mutations in a particular motif in both the GIL01 and GIL16 phages are also shown to correlate with a switch to the lytic cycle. Interestingly, both bacterial viruses displayed narrow, yet slightly different, host spectrums. We also explore the hypothesis that pBClin15, a linear plasmid hosted by the Bacillus cereus reference strain ATCC 14579, is also a prophage. Sequencing of its inverted repeats at both extremities and a comparison with GIL01 and GIL16 emphasize its relationship to the Tectiviridae. |
|
|---|
In addition to GIL01, three other tectiviruses that infect gram-positive bacteria have been identified so far. The bacteriophages AP50 and
NS11, isolated from Bacillus anthracis (19) and Bacillus acidocaldarius (23, 24), respectively, have only been morphologically characterized. Bam35 was isolated from B. thuringiensis serovar Alesti in 1978 (1) and was recently sequenced (21), revealing that it differs from GIL01 by 11 bp. A previous distribution analysis of similar prophages among the B. cereus group had unveiled the presence of a GIL01-related molecule in B. thuringiensis strain B16. PCR and restriction analyses performed with this element, named GIL16, indicated that it differs from GIL01 (28). In order to further analyze this isolate, we sequenced GIL16, and we report its sequence in this study. A genome comparison effectively showed a high degree of conservation with GIL01. We also demonstrate here that both phages have an extremely narrow host range and that GIL16 displays a slightly different spectrum than that of GIL01.
Furthermore, recent genome sequencing of the Bacillus cereus reference strain ATCC 14579 revealed the sequence of pBClin15 (11), a 15.1-kb linear plasmid that was previously detected in total DNA preparations from this bacterium (7). Although no experimental evidence is currently available regarding its nature, its genetic organization, similar to that of GIL01 and GIL16, strongly supports the hypothesis that it is a prophage. Moreover, we report the existence of unresolved inverted terminal repeats (ITRs), the common end structure of linear plasmids and bacteriophages.
Considering that all gram-negative tectiviruses identified so far share a very high level of sequence identity, it seemed interesting to investigate whether this close relationship could also be found among members that infect gram-positive bacteria. The present study reports a comparison at the genetic level of three tectiviral bacteriophages, GIL16, GIL01, and Bam35, and of the pBClin15 plasmid.
|
|
|---|
Phage propagation. Phage stocks were obtained by centrifuging lysogenic cell cultures at 2,500 x g for 15 min and then filtering the supernatants (0.45-µm-pore-size filter). Exponentially growing cell cultures (200 µl) were infected with 200 µl of a series of diluted (1:10 to 1:100) phage stocks, and the mixtures were incubated at room temperature for 30 min. The phage-bacterium mixtures were plated with molten top agar (0.7%) onto nutrient or Trypticase soy broth plates and incubated at 30°C overnight. Single turbid plaques were used for PCR screening. Clear plaques were picked, and pure lines were obtained by at least three consecutive single plaque purification steps. Clear plaque (cp) variants were used to make lytic stocks by collecting the soft agar (0.3%) from semiconfluent plates. Cell debris was removed by centrifugation (2,500 x g, 15 min). The supernatants were collected for use in phage DNA extractions and PCRs.
DNA manipulation and sequencing strategy. A clear plaque mutant designated GIL16c was isolated from the centers of plaques of B. thuringiensis serovar Israelensis GBJ002 infected with GIL16. After three runs of single-colony isolations, the phage GIL16c DNA was extracted as previously described for phage GIL01 (28). Two PCR products, of 500 and 2,500 bp, obtained by the use of primer pairs specific for GIL01, were first cloned into the positive selection vector pCR4-TOPO (Invitrogen) and then were sequenced. The remaining sequences of GIL16c were obtained by direct primer walking sequencing of the phage DNA. pBClin15 was isolated by pulsed-field gel electrophoresis (PFGE) of strain ATCC 14579 chromosomal DNA and was extracted by use of a QIAquick gel extraction kit (QIAGEN). Subsequent DNA purification enabled direct sequencing of the extremities by use of the pBClin15-specific primers pBClin-L (5'-CATGCTATGTCATGTTTTGAC) and pBClin-R (5'-CGTTAGAGAAGTTGACAGGTG). DNA purifications were performed with a QIAquick PCR purification kit (QIAGEN).
PFGE of genomic DNA. Chromosomal DNAs were prepared according to the method of Léonard et al. (15), with the exception that the LMP (GIBCO-BRL) agarose-embedded DNA plugs were not digested. PFGE was performed with the CHEF-DR II system (Bio-Rad) at 14°C and 4.5 V cm1 in a 1.0% agarose (Sigma) gel in 0.5x TBE (45 mM Tris-borate, 1 mM EDTA). Electrophoresis was performed for 20 h, with the pulse divided into two phases, from 16 to 19 s for 10 h and from 50 to 55 s for 10 h. After being stained with an ethidium bromide solution (1 µg ml1), the gel was analyzed with the BioCaptMW, version 99.053, system.
Electron microscopy analysis. Following exposition to UV light (254 nm for 10 s) and incubation at 30°C for 4 h, bacterial debris was cleared by centrifugation. The supernatant, containing viral particles, was filtered through a 0.45-µm-pore-size filter and then precipitated with a polyethylene glycol 6000 (10%) and NaCl (5 mM) solution. The phage pellet was then resuspended in Tris-HCl (pH 7.5), deposited on carbon-coated Formvar grids, stained with 1% potassium phosphotungstate (pH 7.2), and analyzed under an electron microscope (JEOL) at 80 kV.
Curing experiments. The curing of strain ATCC 14579 was achieved by heat treatment. The strain was initially grown at 30°C in LB medium for a few cycles. This culture was then shifted to 44°C by use of a 1% (vol/vol) inoculum. Growth at this temperature was repeated every 4 h for 2 days with fresh medium, and upon transfer, samples were taken, diluted, plated onto LB medium, and incubated at 30°C for 12 h. Single colonies were isolated and subjected to molecular and plasmid profile analyses by PCR and PFGE, respectively.
PCR amplification. Phage plaques, or bacterial lawns for the negative control, were recovered from titration plates and diluted in 100 µl of 0.9% NaCl. PCRs were then performed with 1 µl of these samples by use of the GIL01-specific primers Lex1 (5'-GGATCCATGTTGACGCCAAGGG) and R2GIL01 (5'-AAGCTTCAGTCATCCTTCTTCCC) or by use of Lex1 and the GIL16-specific primer R2GIL16 (5'-AAGCTTCTAGTCCTTTTCCGCATTTTC). Cell lines that were cured of pBClin15 were screened by the use of pBClin1 (5'-GTCAAAACATGACATAGCATGC) and pBClin2 (5'-CTACCGCGATTTCGTTTACC). The cereolysin O gene (1,146 bp) was amplified with the primer pair Hem4.1 (5'-ACGTCACCAGTNGATATWTC) and Hem4.2 (5'-TCTCCACCATTCCCAWGCAAG).
Computational analyses. DNA and protein sequences were analyzed with the EMBOSS package at the Belgian EMBL Node and with the Accelrys DS, version 1.5, gene program. Possible homologies to known proteins were searched with PSI-BLAST (http://www.ncbi.nlm.nih.gov/BLAST/) and ACLAME (http://aclame.ulb.ac.be). The Bam35 genomic sequence was obtained with accession number AY257527, and pBClin15 original and recently updated sequences were obtained with the accession numbers NC_004721 and AE016878, respectively.
Nucleotide sequence accession number. The GIL16c nucleotide sequence has been deposited under GenBank accession number AY701338.
|
|
|---|
GIL16 is able to produce small turbid plaques on B. thuringiensis GBJ002, which is also the native host for GIL01 (28). In order to avoid any contamination by other bacteriophages potentially present in strain B16, we isolated B. thuringiensis GBJ002 lysogens for GIL16 and confirmed the presence of the pGIL16 prophage by PCR. The release of viral particles by GBJ002 lysogens was then analyzed by electron microscopy analysis. As shown in Fig. 1, GIL16 viral particles displayed a similar morphology to those of other identified tectiviral phages, such as PRD1 (3), AP50 (20), and Bam35 (1, 27). These tectiviruses are characterized by an icosahedral capsid with a diameter of about 50 to 60 nm and are occasionally associated with a tail-like structure serving as a DNA ejection device during infection. The latter could clearly be distinguished on the micrograph (Fig. 1) and had the same size as that previously described for other tectiviruses (1, 10, 20, 21).
![]() View larger version (81K): [in a new window] |
FIG. 1. Electron micrograph of GIL16 particles with a tail-like structure. Bar, 50 nm.
|
Thr or Val
Leu). Considering that all of these modifications do not significantly alter the genetic organization of either bacteriophage, GIL01 may be considered a subspecies of Bam35. |
View this table: [in a new window] |
TABLE 1. Nucleotide differences between the GIL01 and Bam35 genomes
|
|
View larger version (20K): [in a new window] |
FIG. 2. ITRs at the GIL16 extremities and comparison with those of GIL01. Conserved residues within the GIL16 ITRs are indicated by full vertical lines. Conserved residues between the GIL01 and GIL16 left and right ends are shown by dashed vertical lines.
|
|
View this table: [in a new window] |
TABLE 2. Comparison of GIL16c, GIL01, and pBClin15 predicted gene products
|
|
View larger version (15K): [in a new window] |
FIG. 3. Genetic map of the GIL16c and GIL01 bacteriophages and of pBClin15. Predicted ORFs are depicted by block arrows and ORFs sharing similarities with known proteins are shown in dark gray. DNA pol, DNA polymerase; LexA, LexA-like repressor; Pack, DNA-packaging protein. Similar ORFs are connected by vertical lines. ORFs shared by two genomes are shown in pale gray, while unique ORFs are shown in black.
|
GIL16c, GIL01, and pBClin15 phage-borne predicted functions. So far, among the tectiviruses that infect gram-positive bacteria, only the endolysin-degrading activities of GIL01 have been experimentally confirmed by zymogram analysis (29). GIL16c ORF27 shares 95.6% identity with the Mur1 endolysin encoded by GIL01 (ORF26) and 74% identity with the pBClin15 (ORF23) putative endolysin. The high degree of similarity among these three hydrolases suggests that they possess comparable lytic activities. It is noteworthy that most of the variations observed between the pBClin15 predicted endolysin and those encoded by GIL01 and GIL16 are located at the C-terminal end. Since that domain has been shown to be implicated in cell wall binding (14, 17), the observed differences may reflect the target specificity of these endolysins.
GIL16c ORF5 showed similarity to several B-type DNA polymerases. Its deduced gene product is highly conserved in GIL01 (97.3% identity), whereas only 47.5% identity was observed with its counterpart in pBClin15. All of the conserved motifs of this protein family are found in the three proteins. The B-type family includes several DNA polymerases encoded by linear plasmids carrying genes for terminal proteins at their extremities. Linear molecules, such as
29 and PRD1, are replicated via a protein-priming mechanism in which the terminal proteins are used as primers to initiate replication. The presence of genes at the 5' extremities of the GIL01 genome was experimentally confirmed (28), supporting the idea that the two related linear molecules GIL16 and pBClin15 may also possess a similar end structure.
One interesting BLAST result was the DNA-packaging protein P9 of PRD1, the prototype member of the Tectiviridae. GIL16c ORF13, GIL01 ORF14, and pBClin15 ORF12 had 26, 21, and 27% identities, respectively, with this protein. Protein P9 is thought to be involved in PRD1 genome encapsidation by providing the required energy (9, 18). The hypothesis that these predicted proteins correspond to DNA-packaging proteins was reinforced by the presence of an ATP-binding consensus motif (28) in all three ORFs (Fig. 4).
![]() View larger version (78K): [in a new window] |
FIG. 4. Amino acid comparison between the putative DNA-packaging proteins of GIL16 (ORF12), GIL01 (ORF13), and pBClin15 (ORF12) and the P9 protein of PRD1. The ATP-binding motif (GXXGXGKXXXXXXXL) is framed, and conserved residues are indicated by asterisks. All residues of GIL16 that are present in at least two other ORFs are indicated by gray boxes.
|
Occasionally, clear plaques can be observed among the turbid plaques formed by GIL01 and GIL16 on GBJ002 lawns. Clear plaques are produced by lysogeny-defective phages, designated clear plaque (cp) mutants. These phages are characterized by mutations either in the repressor or in the operator sequence recognized by the repressor, leading to permanent expression of the lytic genes and the production of clear plaques. The GIL16 sequence reported here corresponds to that of a cp mutant (GIL16c). The region that is supposed to code for the putative repressor as well as the downstream predicted gene (ORF7) was sequenced from a wild-type (wt) particle and several cp mutants in order to identify putative variations. No mutation was observed in ORF6, which encodes the putative repressor, but two types of mutations were identified in a specific region of ORF7. As shown in Fig. 5, wt GIL16 harbors an 11-bp motif that is repeated twice with just one mismatch. In the cp1 mutant, corresponding to GIL16c, one of these 11-bp motifs was deleted, while in the cp2 mutant, the mutation was a repetition of a 28-bp motif that is present once in the wt. Interestingly, both mutations occurred at nearly the same position in ORF7. The corresponding regions in wt GIL01 and in two cp mutants have also been investigated, giving similar results. The GIL01 cp1 and cp2 mutants harbored a repetition or a deletion of the 11-bp motif identical to that implicated in the GIL16 mutation (Fig. 5). It is intriguing that all of the isolated GIL16 and GIL01 cp clones had an insertion or deletion of a particular motif in ORF7, while no mutation was ever observed in ORF6. These motifs may be part of repressor binding sites controlling the lytic genes. Since it is common within phage genomes to find genes encoding DNA-binding proteins situated close to the site of action of the gene product (12), the possibility of ORF6 and/or ORF7 being involved in the regulation of the lysogenic cycle should be explored. Future research will consist of establishing the precise role of these ORFs as well as identifying potential repressor binding sites across the phage genomes.
![]() View larger version (36K): [in a new window] |
FIG. 5. Mutations observed in the ORF6-ORF7 region of GIL16 and GIL01 clear plaque (cp) mutants compared to the wt sequence. Motifs affected by the mutations are indicated by white letters over a black background (mismatches are indicated by gray boxes), and mutations observed in each cp mutant are indicated by black letters over a gray background.
|
![]() View larger version (96K): [in a new window] |
FIG. 6. PFGE patterns of undigested genomic DNAs from B. cereus and B. thuringiensis strains used in this study. The molecular size markers are the lambda 48.5-kb size marker (Sigma) (left) and a Gene Ruler mix (Fermentas) (right).
|
|
View this table: [in a new window] |
TABLE 3. Host range of GIL01, GIL16, and their corresponding clear plaque mutants
|
The recent genome sequencing of Bam35c (21) revealed 99% identity with GIL01, suggesting that these two particles are virtually the same. Nevertheless, the present study revealed that GIL01 does not infect HER1047, a strain which is sensitive to Bam35. It would therefore be pertinent to reexamine the Bam35 host range alongside with GIL01 and GIL16 to rule out the possibility that HER1047 has not mutated since 1978. Surprisingly, GIL16 was able to infect HER1047 and thus displayed a slightly different host spectrum than that of GIL01. While PRD1 and its relatives have a broad host spectrum due to the wide distribution of their receptor-encoding plasmid, what determines the narrow range of gram-positive bacterium-infecting tectiviruses remains to be seen.
pBClin15 is delineated by imperfect ITRs. Sequencing of the B. cereus reference strain ATCC 14579 led to the identification of a 15.1-kb linear plasmid, pBClin15 (11). This molecule has a size and genetic organization comparable to those of GIL01, suggesting that it might also be a prophage. Yet the most intriguing distinction of this plasmid relative to tectiviruses remained the absence of ITRs. The hypothesis that these inverted repeats might exist at both extremities prompted us to perform sequencing of both ends. pBClin15 was isolated as described in Materials and Methods, and subsequent runoff sequencing at the extremities gave the whole sequence of this linear molecule. Sequences of 100 and 74 bp were missing from the left and right extremities, respectively, bringing the actual size of pBClin15 to 15,274 bp instead of 15,100 bp as originally reported (11). The pBClin15 sequence has since been corrected (AE016878). The alignment of the resulting sequences revealed the existence of imperfect 74-bp ITRs that share 64.6% identity with each other (Fig. 7). In contrast to those of GIL01 and GIL16, the pBClin15 inverted repeats include several gaps and mismatches, and they only share 51 to 58% identity with those of GIL01 and GIL16 (data not shown).
|
View larger version (7K): [in a new window] |
FIG. 7. Inverted terminal repeats of pBClin15. Conserved nucleotides are indicated by vertical lines.
|
wiThis work was supported by grants from the National Fund for Scientific Research (FNRS, Belgium). C.V. and N.F. hold research fellowships from FRIA (Fonds pour la Formation à la Recherche dans l'Industrie et l'Agriculture).
C.V. and N.F. contributed equally to this work. ![]()
Present address: Center for Molecular Genetics and Division of Biological Sciences, University of CaliforniaSan Diego, La Jolla, CA 92093-0322. ![]()
|
|
|---|
ius, and D. H. Bamford. 2002. Sequential model of phage PRD1 DNA delivery: active involvement of the viral membrane. Mol. Microbiol. 46:1199-1209.[CrossRef][Medline]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»