Previous Article | Next Article ![]()
Journal of Bacteriology, February 2007, p. 1482-1487, Vol. 189, No. 4
0021-9193/07/$08.00+0 doi:10.1128/JB.01111-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Département de biochimie et de microbiologie, Faculté des sciences et de génie, Groupe de recherche en écologie buccale, Faculté de médecine dentaire, Félix d'Hérelle Reference Center for Bacterial Viruses, Université Laval, Québec, Canada G1K 7P4
Received 25 July 2006/ Accepted 2 October 2006
|
|
|---|
|
|
|---|
The evolution of phages infecting the low-G+C-content gram-positive bacterium Lactococcus lactis is the subject of ongoing studies because of the economic value of the host strains in fermented dairy products, as well as the frequent emergence of new virulent phages that are responsible for delays in milk fermentation. Lactococcal phages have been reclassified recently into 10 genetically distinct groups of dsDNA and tail-containing phages (15). However, members of only three L. lactis phage groups (936, c2, and P335) are regularly isolated. While virulent members of the 936 and c2 groups are rather homogeneous, there is considerable genetic heterogeneity in members of the P335 group, which contains both temperate and lytic phages (15).
One effective way to control lactococcal phages in dairy processes is through the use, in rotation, of L. lactis strains harboring phage defense mechanisms (8, 26). These mechanisms are divided based on their general modes of action, as follows: inhibition of phage adsorption, DNA ejection blocking, restriction/modification systems, and abortive infection mechanisms (Abi). The latter systems block phage multiplication and cause premature cell death upon phage infection. With the constant use of L. lactis strains carrying Abi systems, new phages resistant to Abi systems have emerged, but these phages remain largely uncharacterized (4, 5, 16, 17, 19, 27).
In this study, the evolution of the lytic lactococcal phage ul36 (P335 group) was studied (25). Phage ul36 was sequentially propagated on a prophage-containing L. lactis strain (SMQ-86) harboring either AbiK (4-7, 20, 21) or AbiT (3). Phage mutants resistant to the Abi systems were isolated and characterized by genome sequencing and electron microscopy analyses.
Isolation and characterization of phage mutants. Lactococcal phage mutants resistant to Abi systems were isolated in the laboratory from GM17 lysates (27) as shown in Fig. 1A. The wild-type lytic phage ul36 (GenBank accession number AF349457) was first propagated on L. lactis SMQ-86 containing the AbiK system. A phage mutant resistant to AbiK was isolated at a frequency of 107 and designated ul36.k1. This phage was then propagated on the same L. lactis host carrying AbiT. Another phage mutant (ul36.k1t1) was obtained at a frequency of 108, and this mutant was not sensitive to either Abi system. In parallel, a similar experiment was conducted by first propagating ul36 on an AbiT-containing host strain and then propagating it on an AbiK+ strain. Two phage mutants, ul36.t1 and ul36.t1k1, were isolated at frequencies of 108 and 107, respectively.
![]() View larger version (32K): [in a new window] |
FIG. 1. Electron micrographs (A) and EcoRI restriction profiles (B) of lactococcal phages analyzed in this study. Lane M contained markers. The method used to isolate the mutant phages is also shown.
|
![]() View larger version (26K): [in a new window] |
FIG. 2. Alignment of the genetic maps of wild-type phage u136, mutant phages u136.k1, u136.k1t1, u136.t1, and u136.t1k1, prophage smq86, and host strain L. lactis SMQ-86 (4). Identical ORFs are the same color. Cross-hatched regions are unique among the phage genomes shown.
|
Phage ul36.k1t1. The AbiK- and AbiT-resistant phage mutant ul36.k1t1 has a genome that contains 35,594 bp and has an overall G+C content of 35.7% (GenBank accession number DB394807). Comparative analysis indicated that only 27% of the genome of phage ul36.k1t1 is present in wild-type phage ul36, while 34% is present in ul36.k1. However, 23,430 bp of the phage ul36.k1t1 genome was not found in the genomes of ul36.k1 and ul36. The divergent region contains all the genes involved in phage morphogenesis and cell lysis, which may explain the differences in tail length and in the structure of the baseplate.
Phage ul36.t1. The genome of the AbiT-resistant ul36.t1 phage contains 35,992 bp and has an overall G+C content of 36.0% (GenBank accession number DB394808), and it is 806 bp shorter than the genome of wild-type phage ul36. Sequence analysis showed that 36% of the ul36.t1 genome is present in the genome of ul36, confirming the relationship of these phages (Fig. 2). The remaining 23,179 bp (64%) of the ul36.t1 genome is different from the phage ul36 genome, indicating that there was significant genomic modification. Again, the morphogenesis and lysis modules were divergent in the two phages, which is in agreement with the morphological observations (Fig. 1A). Interestingly, the genome of phage ul36.t1 exhibits 78% identity with the genome of phage ul36.k1t1, indicating that these two phages acquired the same DNA. The difference between ul36.k1t1 and ul36.t1 was mainly in the area containing the genes coding for proteins involved in DNA replication and transcription.
Phage ul36.t1k1. The genome of the AbiT- and AbiK-resistant phage mutant ul36.t1k1 contains 34,897 bp and has an overall G+C content of 35.8% (GenBank accession number DB394809). This phage has the shortest genome of the five lytic phages analyzed here. Comparative analysis demonstrated that the ul36.t1k1 genome exhibits 84% identity with the genome of ul36.t1 (Fig. 2). The main difference was again in the genes coding for proteins involved in DNA replication and transcription. Surprisingly, ul36.t1k1 contains only 21% of the original core DNA from ul36.
Taken altogether, these results showed that the four derivatives of phage ul36 picked up new DNA. The acquisition of such large DNA segments is rather exceptional, and it occurred without affecting the functional capacities of the resulting phage hybrids. This clearly indicates that there are exchangeable and compatible modules. Thus, some of the gene products are highly fit to work together, which would be expected for proteins necessary for the phage structure. Therefore, the functional new DNA must have come from another phage. This prompted us to sequence the complete genome of the only known inducible prophage of L. lactis SMQ-86,
smq86.
Prophage
smq86.
Following induction with mitomycin C (5 µg/ml), the prophage was purified for electron microscopy observation, and its DNA was isolated for sequencing. Prophage
smq86 had the same morphological features as the lytic phages ul36.t1, ul36.k1t1, and ul36.t1k1 (Fig. 1A). The genome of prophage
smq86 contains 33,641 bp and has an overall G+C content of 36.0% (GenBank accession number DB394810). Thus, the genome of
smq86 is smaller than the genomes of five lytic phages. Fifty-one open reading frames (ORFs) consisting of 40 or more codons were identified, and functions were attributed to 17 ORFs based on homology with proteins with putative functions or conserved domains (Table 1). The genome of
smq86 is organized into different functional modules, as observed in other phages (10, 14, 25) (Fig. 2). All ORFs are on the same strand, except for the four genes (orf1 to orf4) of the lysogeny module. Many gene products exhibit homology with counterparts found in other phages belonging to the P335 group, prompting classification of
smq86 in this group (15). However, the deduced proteins involved in packaging and capsid morphogenesis exhibit no similarity to known proteins of L. lactis phages. Comparative nucleotide sequence analysis confirmed that most of the DNA acquired by the new lytic phages obtained in this study came from
smq86.
|
View this table: [in a new window] |
TABLE 1. Features of prophage smq86 ORFs and putative functions of their products
|
smq86, ul36.t1, and ul36.t1k1.
All of the new DNA picked up by ul36.t1k1 came from prophage
smq86, except for a 100-bp region that is repeated three times from coordinates 8216 to 8315, 8316 to 8415, and 8416 to 8515 in ul36.t1k1. The same observation was made for ul36.t1 as a short 1,165-nucleotide insertion in orf46, possibly coding for the receptor-binding protein (Fig. 2), was not found in the inducible prophage. But this 1.1-kb region of ul36.t1 exhibits homology with the rbp gene of other L. lactis phages, suggesting that a second recombination event took place with another prophage (most likely defective) of L. lactis SMQ-86. It is known that rbp genes are the site of frequent DNA shuffling that favors the generation of phage variants with altered host ranges (18, 22, 28).
Comparisons between
smq86, ul36.k1, and ul36.k1t1.
As indicated previously, the exchanged DNA (2.5 kb) in phage ul36.k1 has been investigated (4). Phage ul36.k1t1 contains this 2.5-kb segment but also acquired a large DNA fragment from prophage
smq86 (Fig. 2). However, another region (2.4 kb) covering the rbp gene (orf46) and the lysis genes was unique to ul36.k1t1 (Fig. 2). Again, this segment is likely the result of another recombination event with a resident prophage in L. lactis SMQ-86.
Comparison between
smq86 and ul36.
Pairwise comparison of the entire genomes of the lytic phage ul36 and the inducible prophage
smq86 revealed a very limited number of identical nucleotides that could serve as substrates for homologous recombination. Most matching DNA regions were between coordinates 11,000 and 14,000 in ul36. This finding supports the hypothesis that one of the recombination events in phages ul36.k1t1 and ul36.t1 occurred by homologous recombination in this region. Otherwise, it appears that the other recombination events were illegitimate.
Influence of Abi systems and strain rotation on phage evolution. Viruses are known to mutate during amplification cycles, and the presence of various selective pressures, such as host diversity and barriers, influences the nature of the infectious population. Lactococcal phages typically face such shifting selective pressures in the dairy environment. It has been reported previously that two general types of lactococcal phage mutants resulting from the selective pressure of Abi systems can be isolated (13). The first class of Abi-resistant phage mutants carry only point mutations (5, 13), while the second class is the result of recombination with phage-related sequences present in the host chromosome (4, 19, 27). All the mutant phages isolated in this study belonged to the second class of mutants. To different extents, they all acquired new DNA by homologous or illegitimate recombination, mostly from an inducible prophage. The percentage of DNA exchanged in the genomes ranged from 6.8% (ul36.k1) to a remarkable 79% (ul36.t1k1). The frequency at which the ul36 derivatives were obtained (107 to 108) suggests that they were already present in the high-titer lysates used to challenge the Abi-containing strains. In fact, these lysates likely contained a population consisting mostly of phages with the ul36 genetic make-up but also containing some functional and nonfunctional derivatives carrying genetic variations. In the presence of a specific selective pressure, the fittest organisms rapidly multiply, while the other organisms remain at low levels or eventually are eliminated. This study provided biological evidence that Abi systems and prophage DNA can significantly influence the genetic make-up of lytic phages. It also demonstrated that the genome plasticity of lactococcal phages allows them to rapidly adapt to new environments. In fact, the practice of rotating isogenic strains carrying different antiphage systems appears to significantly contribute to the emergence of new lytic phage variants, at least for P335-like phages. Interestingly, sequential encounters with a specific phage resistance mechanism also determine the direction of phage evolution.
Phage-bacterium coevolution. Since prophage elements are found in many bacterial genomes (11), the probability that phages will infect a cell and recombine with a prophage genome(s) by homologous or illegitimate recombination is rather high. This genomic reshuffling inexorably leads to the emergence of new phages and quite possibly to new bacterial strains as well. Considering the astonishing genomic rearrangement observed in some of the phage ul36 derivatives and the presence of numerous phage defense systems in bacteria, this phage-bacterium coevolution may be more significant than previously envisaged.
This study was funded by strategic grants from the Natural Sciences and Engineering Research Council (NSERC) of Canada.
Published ahead of print on 13 October 2006. ![]()
|
|
|---|
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»