Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2005, p. 1604-1611, Vol. 43, No. 4
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.4.1604-1611.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Infectious Diseases Laboratories, Institute of Medical and Veterinary Science, Adelaide, South Australia, Australia
Received 21 September 2004/ Returned for modification 18 November 2004/ Accepted 6 December 2004
|
|
|---|
|
|
|---|
There are currently a number of potential problems with both classical and most molecular-based methods for typing Salmonella. Phage typing requires a well-maintained phage library, precise methodology, and experience in interpretation of results (8, 30). Interpretation of plaque patterns is often subjective and prone to variation between laboratories, especially between phage types with similar reaction patterns. For example, serovar Typhimurium definitive type 12 (DT12) and DT108 are phenotypically similar phage types with regard to lysis patterns of serovar Typhimurium and are separated on the Anderson typing panel by the sensitivity of DT108 with bacteriophages A21, A27, and A35. In many Australian isolates of DT108, there is no observed reaction with phage A21. Furthermore, there are cases where there are very weak reactions with phages A27 and A35.
Perceived problems with phage typing are further compounded, as a significant number of isolates can be untypable by this method. Often, one phage type may predominate in a geographical area, making discrimination between outbreak and nonoutbreak isolates difficult (21, 30). Molecular methods such as PFGE, although less subjective than phage typing, at times do not provide clear discrimination between isolates. Furthermore, some isolates are not typable by PFGE, for example, Salmonella enterica serovar Kentucky strains and Salmonella enterica subsp. salamae serovar Sofia (8, 22).
Multilocus sequence typing (MLST) is a relatively new high-resolution typing system developed for both evolutionary and epidemiological studies of microorganisms (6, 8, 10). MLST involves the nucleotide sequencing of approximately 400-bp regions of at least seven genes, usually housekeeping genes. Although multilocus enzyme electrophoresis (MLEE) is similar to MLST in that both systems target enzymes or their genes, MLST has advantages over MLEE. MLST will detect single-nucleotide polymorphisms, especially third-base redundancy within a codon, whereas MLEE relies on nonconservative changes in the amino acid sequence of a target enzyme that result in a change of its electrophoretic mobility on a starch gel (28). For example, Kotetishvili et al. (8) found that the ratio of nonsynonymous (amino acid substitution) to synonymous base changes (no change in amino acids) in three Salmonella housekeeping genes was >1. This meant that there was a low incidence of amino acid substitutions that would have affected enzyme mobility in MLEE studies.
Sequencing of Salmonella genes has generally been restricted to comparing sequences of housekeeping genes between the major groups of Salmonella or serovars of S. enterica subsp. enterica, including serovars Typhimurium and Enteritidis. Sequence variation can be demonstrated for geographically and epidemiologically unrelated isolates (5, 8, 16, 34). Consequently, MLST is a useful tool for studying evolution and global epidemiology of the salmonellae. However, for investigations of outbreaks, the fine discrimination of closely related isolates based on MLST of housekeeping genes may be not be possible, due to the high identity of the housekeeping gene sequences.
The presence of genomic sequences derived from temperate phages in Salmonella may offer an alternative target for MLST of Salmonella (25, 26). Although a typical prophage genome only constitutes approximately 1% of the bacterial chromosome, these loci are subject to mutational events at a higher rate than housekeeping genes. Furthermore, recombination between phages has resulted in an array of mosaic genetic structures of many bacteriophages, both within individual genes and in the overall genome organization (3, 17, 24, 32). Indeed, it has been stated that prophage genes account for most of the genetic diversity among closely related Salmonella strains (7). The presence of immunity exclusion genes such as the cI repressor gene of the immC region of some phages may confer a selective advantage for the host cell, rendering the cell immune to superinfection by lytic phage (20). This suggests that prophage genes, both in their presence or absence in a host cell and in their nucleotide sequence, vary considerably within a serovar. Consequently, heterogenicity of prophage genes is likely to be greater than that of host housekeeping genes, making them superior targets for MLST.
Previous studies by our laboratory have identified two bacteriophages (ST64T and ST64B) induced from a strain of serovar Typhimurium DT64. These two temperate phages were separated from each other by a cesium chloride gradient. The two bacteriophages were fully sequenced, and the sequences were deposited in GenBank (www.ncbi.nlm.nih.gov) under accession numbers AY052766 (ST64T) and AY055382 (ST64B) (14, 15). Southern hybridization studies have suggested that these bacteriophages are present in a number of serovar Typhimurium phage types as well as other S. enterica subsp. enterica serovars, including Enteritidis, Virchow, Heidelberg, and Hadar (13, 31). Based on the widespread nature of these phages, we used this sequence data to design a number of primer sets to analyze sequence variations of prophage loci. As well as using locus sequences from bacteriophages ST64B and ST64T, we also examined a range of similar loci from the well-characterized bacteriophage P22 (18, 32). A number of Salmonella housekeeping genes were also assessed to compare their level of sequence identity to the level of sequence identity of the prophage loci. We used this data to discriminate DT12 and DT108 strains of serovar Typhimurium phenotypically that were closely related by lysis pattern, in addition to DT126-related food poisoning by MLST and by amplification of phage loci alone.
|
|
|---|
Thirty-seven DT126 isolates were also included in the study. This group comprised 13 isolates from a 2001 restaurant outbreak in New South Wales, Australia. Another 10 DT126 isolates were obtained during an outbreak in 2003. Four of these isolates were designated DT126 var, due to a variation in reaction to the Anderson typing panel. A further 14 DT126 isolates were epidemiologically unrelated.
All serovar Typhimurium isolates used in this study were provided by the Australian Salmonella Reference Centre, Institute of Medical and Veterinary Science, Adelaide, South Australia. Serotyping had previously been undertaken using the Kaufmann-White scheme, and bacteriophage typing was performed using the Anderson scheme of 31 bacteriophages (1), both by the Australian Salmonella Reference Centre.
Amplification of phage loci. PCR of Salmonella was performed directly from cell lysates. Isolates were grown overnight in Luria-Bertani broth at 37°C with gentle shaking. A 30-µl PCR mixture was prepared as follows: 3.0 µl of 10x MgCl2-free buffer, 1.0 µM each forward and reverse primer, 3.0 µl of 200 µM each deoxynucleoside triphosphate, 1.5 mM MgCl2, 1.0 U of Taq polymerase, and 2.0 µl of overnight cell culture. The volume was made up to 30 µl with H2O. Buffer, MgCl2, and Taq polymerase were supplied by Roche. Primers were designed based on published sequences (Table 1.). All primers were supplied by Geneworks, Adelaide, South Australia. The deoxynucleoside triphosphates were supplied by Amersham Biosciences, Piscataway, N.J..
|
View this table: [in a new window] |
TABLE 1. Primers used in this study
|
The amplification product was detected and prepared for sequencing by running 5.0 µl of PCR product on a 2.0% agarose gel (Progen, Darra, Queensland) in 1.0x TBE buffer (23) at 5.0 V/cm with pUC19 digested with HpaII (Biotech, Belmont, Western Australia) as a marker. Bands were visualized with UV light after being stained with ethidium bromide. The remaining 25-µl PCR product of positive samples was prepared for sequencing by passage through a QIAquick PCR purification column (QIAGEN, Hilden, Germany) and collection in 30 µl of elution buffer per the manufacturer's instructions. Amplicons were stored at 20°C prior to sequencing.
MLST. Sequencing was performed in both directions with Big Dye Terminator, version 3-1 (Applied Biosystems, Foster City, Calif.). A 20-µl reaction mixture comprised 4.0 µl of BD3-1 master mix, 1.5 µl of 3.0 mM either forward or reverse primer, 12.5 µl of H2O, and 2.0 µl of template. Sequencing was performed with a Corbett Research PC-960G gradient thermal cycler and comprised 25 cycles with the following parameters: 96°C (30 s), 50°C (15 s), and 60°C (4.0 min). The sequence product was precipitated, washed with 75% isopropanol, and then dried. Sequencing was performed with an Applied Biosystems 3700 DNA analyzer.
PFGE. The protocol for PFGE followed that of Maslow et al. (11). Briefly, cells grown overnight in brain heart infusion broth (Oxoid, Basingstoke, United Kingdom) were embedded in agarose and lysed by incubation of the plug in 4 ml of lysis buffer supplemented with 4 mg of lysozyme (Roche, Mannheim, Germany) and 80 µg of RNase/ml. Plugs were then digested with proteinase K and washed, and the DNA was digested overnight with XbaI restriction endonuclease (New England BioLabs, Beverley, Mass.). The plugs were placed into the wells of a 1% agarose gel prepared with PFGE-grade agarose (Bio-Rad Laboratories, Hercules, Calif.) in 0.5x TBE buffer. Staphylococcus aureus strain NCTC 8325 digested with SmaI was used as a molecular size marker (29). The PFGE was run on a Bio-Rad CHEF-DR III system for 19 h at 6.0 V cm1 at 4°C, with an initial switch time of 2.0 s and a final switch time of 50 s. After being run, the gel was stained in ethidium bromide, destained in water, and photographed under UV light.
Data analysis. Isolates were initially separated based on PCR results. Each primer set was given a number, and data were entered into a spreadsheet consisting of all positive results for each strain. The spreadsheet was imported into GelCompar IV (Applied Maths, Kortrijk, Belgium) with appropriate formatting. A dendrogram of isolate PCR profiles was generated by a Dice coefficient and clustering by unweighted-pair group method using average linkages (UPGMA) (Fig. 1).
![]() View larger version (29K): [in a new window] |
FIG. 1. Dendrogram of PCR profiles for separation of serovar Typhimurium isolates based on MAPLT. The dendrogram was generated by Dice coefficient with clustering by UPGMA, based on the presence or absence of amplified product. A total of six clusters of isolates with identical PCR profiles were generated, as well as a total of 21 isolates with unique profiles. Boldface numbers in the dendrogram refer to cluster numbers. Cluster 1 comprises the eight non-DT126 isolates, and cluster 5 comprises the four DT126 var isolates. Abbreviations: N.S.W., New South Wales; N.T., Northern Territory; Qld., Queensland; S.A., South Australia; Vic., Victoria.
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. The incidence of loci, number of alleles, and their distribution in 73 S. enterica serovar Typhimurium isolates
|
MLST analysis. All PCR products for all relevant primer sets were sequenced with both forward and reverse primers, and the number of different alleles for each primer set was determined (Table 2). Most bacteriophage-derived primer sets produced at least two different alleles when sequenced, except those primer sets where only a small number of strains generating a PCR product were observed. Some primer sets produced only one allele even when a significant number of isolates produced PCR product, for example, SB28ST64B and the mntST64T gene. In a few cases where significant numbers of isolates produced product, only one or two isolates contained an allele distinct from the majority of isolates. For example, a second SB26ST64B allele was present in only one isolate, 01-9-001.
The housekeeping genes fhuA, glnA, and sucA amplified in all 73 isolates. In contrast, 46 and 15 isolates yielded product for manB and tonB, respectively. More than one allele was present only in glnA and manB. Isolate 01-126-114 possessed a manB allele different from that of the other 45 that were positive with manB primers. Six DT126 and four DT126 var isolates from the 2003 outbreak had a characteristic glnA allele.
MLST of phage loci was undertaken to further discriminate between isolates that clustered as depicted in Fig. 1. Separation of the eight isolates in cluster 1 produced three subgroups (Table 3). Separation was based on a single allele in each case, with the two-outbreak strains possessing a SB37ST64B allele different from that of the other six isolates. Isolate 02-12-002 could be separated from the other seven isolates in this cluster based upon the immCST64B c1 5' sequence.
|
View this table: [in a new window] |
TABLE 3. The separation of MAPLT-clustered non-DT126 isolates by sequencing
|
MLST of the DT126 isolates indicated a lower degree of variation than observed above (Table 4). No allelic variation was observed for the 13 isolates from the 2001 restaurant outbreak or the 10 2003 outbreak isolates for prophage loci. However, separation of the 2003 isolates could be achieved based upon glnA. Of the remaining DT126 isolates only one isolate, 01-126-123, exhibited different prophage alleles from the other DT126 isolates.
|
View this table: [in a new window] |
TABLE 4. Separation of DT126 and DT126 var isolates by MAPLT and sequencing
|
![]() View larger version (72K): [in a new window] |
FIG. 2. Separation of serovar Typhimurium isolates based on PFGE. Dendrogram generated by Dice coefficient with clustering by UPGMA. Five different profiles were generated. Profiles 2 and 3 had >90% similarity, suggesting isolates with these profiles are genetically closely related. Serovar Typhimurium DT126 isolates were either profile 4 or 5; the remaining serovar Typhimurium isolates were found in all five profiles. Marker sizes are in kilobases x 103.
|
|
|
|---|
Amplification of prophage loci provided the first level of typing of the serovar Typhimurium isolates based on the presence or absence of PCR product. We term this method multiple amplification of phage locus typing (MAPLT). PCR of prophage loci in serovar Typhimurium generated 27 distinct profiles (Fig. 1), compared to only 5 profiles generated by PFGE (Fig. 2). The MAPLT profiles comprised either single isolates or clusters of isolates with identical PCR profiles within each cluster. Some MAPLT profiles were closely related with the only product from a single primer set separating them; for example, clusters 1 and were separated by the absence of the SB26ST64B gene amplification in isolates of cluster 2.
DT126 isolates (as well as the single DT9 and DT135 isolates) were clearly different from the other phage types analyzed. There were six MAPLT profiles for DT126, in comparison to PFGE, which only generated two profiles, even though the isolates were from diverse sources. These results demonstrate the limitations of PFGE and the potential of MAPLT as an epidemiological tool when typing closely related strains.
It was also observed that most DT126 isolates contained only ST64B sequences. The apparent lack of other phage-related sequences in this phage type reduced the ability of MAPLT to discriminate between isolates. This observation may partially explain the susceptibility of this phage type to almost all (except A6 and A8) of the bacteriophages in the Anderson typing panel, since no P22- and ST64T-like immunity genes were amplified. In contrast, DT12 and DT108 contain a mosaic of sequences that are related to phages other than ST64B, including genes related to immunity. It is significant that DT12 and DT108 are susceptible to only three and six members of the Anderson panel, respectively. It is this lack of diversity of susceptibility of certain phage types of serovar Typhimurium to the typing panel that can prove problematical in classical phage typing. Typing is then dependent on the more subjective assessment of lysis intensity. The likely presence of a mosaic of phage loci in these DTs makes MAPLT a suitable method for fine discrimination between isolates, as a high number of suitable targets will be present for PCR.
MLST was undertaken with all PCR products, to separate clustered isolates as illustrated in Fig. 1 and to determine which prophage loci exhibit the greatest sequence variability. This information is important, since nucleotide sequence may sometimes be required for further discrimination of isolates. It was observed that the more unrelated isolates by MAPLT showed greater variability of prophage sequences (data not shown). Interestingly, DT170 and DT108 in cluster 3 could not be separated further (Table 3); these DTs can be difficult to separate using the Anderson typing panel with only variable reactions, with phages A14 and A21 providing separation of the two phage types. It could be postulated that these three isolates do indeed represent the same phage type.
Limited separation of the DT126 isolates was achieved by sequencing of prophage loci (Table 4). In addition, PFGE and MAPLT suggest that DT126 was more clonal than the other DTs tested. In 2001, DT126 became the most frequently isolated phage type in Australia and was associated with a number of food-associated outbreaks (2). After 2001, the frequency of DT126 isolation declined to pre-2001 numbers. It is possible that the DT126 isolates isolated during and after 2001 were from a common clone from within Australia or an imported clone. This clone found its way into the food chain, resulting in the increased number of outbreak incidences.
The sequences derived from bacteriophage ST64B have been found to be highly variable in the non-DT126 isolates tested in both amplification and nucleotide sequence analyses (Table 2) (Fig. 3). Genomic sequences of ST64B are widespread in serovars of S. enterica subsp. enterica, including serovar Typhimurium (13, 31). The apparent widespread incidence of ST64B genomic material in serovar Typhimurium observed with MAPLT and sequence data suggests that loci from ST64B will make excellent targets for isolate discrimination within serovar Typhimurium. It is also likely that these sequences could be of use in other Salmonella serovars (C. P. Tucker and M. W. Heuzenroeder, unpublished data). In contrast, the incidence of ST64T and P22 loci was generally restricted to the non-DT126 isolates (Tables 2 and 4). Although discrimination of some isolates was enhanced by targeting ST64T and P22 loci, it should be noted that these loci (a) occur at lower frequency or (b) have limited allelic variation. For example, mntST64T and sieBP22, both of which occur in high frequency in non-DT126 isolates, were found to have little sequence variation (Table 2). While these loci provided less discrimination between isolates tested in this study than those derived from ST64B, they may be potentially useful for isolate discrimination in other DTs or serovars that are yet to be tested.
![]() View larger version (12K): [in a new window] |
FIG. 3. Sequence types (ST) within the immC region of phage ST64B. The base position is the position within the sequenced region; hence, position 1 is position 28047 in GenBank accession number AY055382, which includes ST1. Eight different sequence types including ST1 were detected in the 73 serovar Typhimurium isolates. Where there was no nucleotide difference to ST1, this is indicated by a . The filled boxes indicate regions covered by primers BIM1 to BIM4.
|
Like PFGE, MAPLT and subsequent sequencing of loci did not separate the DT12 and DT108 isolates into groups based on phage type. MAPLT provided a greater level of separation of isolates than PFGE, regardless of phage type. It can be difficult to discriminate between DT108 and DT170, since they react with six and eight members of the typing panel, respectively, and share five reactions, with three being variable or weak. These results suggest that DT170 isolates may be distinguished from the majority of DT108 isolates by MAPLT when phage typing delivers ambiguous results. Further testing with more DT170 isolates will be required to confirm this initial observation.
Ten DT126 isolates were obtained from an outbreak in a restaurant in 2003. MAPLT and sequencing analyses indicated that these isolates were similar to all the other DT126 isolates tested, with the exception of the glnA sequence. When phage typed, four of the isolates did not react with phages A12 and A13, which normally produce semiconfluent lysis on DT126. These four isolates were designated DT126 var, based upon this difference. MAPLT separated DT126 var from the other 2003 outbreak isolates. It could be postulated that DT126 var is the result of phage type conversion or some other genetic event by a mobile element rendering the DT126 var isolates resistant to A12 and A13 typing phages. This demonstrates that when multiple phage-derived primers are used, MAPLT has the power to discern underlying genetic relationships between strains, even when a phage type conversion or acquisition of a mobile element occurs during an outbreak.
In conclusion, this pilot study has shown the potential of utilizing prophage loci for MLST of phenotypically closely related phage types by lysis patterns of serovar Typhimurium. It is clear from the data that phage type does not necessarily reflect a close genetic relationship between isolates. The results show that the genetic variation of prophage genomes makes them excellent targets for PCR amplification and sequencing, compared to housekeeping genes. The variation of frequency of occurrence of many of the loci selected indicates that discrimination of isolates can often be achieved by MAPLT alone, without the need to determine the nucleotide sequence. As this is the first report of the utilization of prophage loci for typing, further work is required to obtain a minimal number of primer sets that will give the maximum level of isolate discrimination. This will mean the examination of genes from other Salmonella phages as a source of new primer sets, for example, A12 and A13, which are genetically distinct from the P22- and ES18-related phages in the Anderson typing panel (27).
This project was undertaken with the generous assistance of the Rural Industries Research and Development Corporation (Chicken Meat Program).
|
|
|---|
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»