Previous Article | Next Article ![]()
Journal of Clinical Microbiology, August 2007, p. 2590-2598, Vol. 45, No. 8
0095-1137/07/$08.00+0 doi:10.1128/JCM.00729-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, United Kingdom,1 Health Protection Agency Centre for Infections, 61 Colindale Avenue, London NW9 5EQ, United Kingdom,2 Scottish Salmonella Reference Laboratory, Stobhill Hospital, 133 Balornock Road, Glasgow G21 3UW, Scotland,3 SAC Veterinary Services and Animal Health Group, Drummondhill, Stratherrick Road, Inverness IV2 4JZ, Scotland4
Received 4 April 2007/ Returned for modification 10 May 2007/ Accepted 14 May 2007
|
|
|---|
|
|
|---|
Many methods have been developed to phenotypically distinguish between S. enterica serovar Typhimurium isolates, including antibiotic susceptibility profiling, phage typing (1), pulsed-field gel electrophoresis (PFGE) (7), and plasmid profiling (26) as well as various PCR-based techniques (8, 15, 19). However, since the genome sequences of several S. enterica strains, including different S. enterica serovar Typhimurium strains, are now available, it should be possible to design rational DNA tools based on fully annotated DNA sequences for use in the field to monitor strain diversity. Here, we have used some of the existing classical typing techniques to analyze a matched collection of S. enterica serovar Typhimurium strains isolated from animal and human sources and have extended these approaches to include DNA microarray analysis. Using these techniques, we have been able to identify and map regions of variation on the chromosome of S. enterica serovar Typhimurium that discriminate between isolates circulating in the same geographical region. Using this information, we have designed multiplex PCR assays that are simple to use and that are able to rapidly distinguish between S. enterica serovar Typhimurium isolates in a cost-effective manner. We believe that similar PCR assays, constructed on the basis of regions of variation in other Salmonella serovars, have the potential to improve the local epidemiological analysis of Salmonella outbreaks.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Origin of 30 animal and 31 human S. enterica serovar Typhimurium isolates
|
Plasmid DNA extraction. Plasmid DNA was extracted from each isolate using a modified Kado and Liu protocol (14). Each isolate was cultured in nutrient broth at 37°C overnight, and then 1 ml of culture was centrifuged at 13,000 rpm for 5 min and the supernatant was removed. The pellet was resuspended in 20 µl Tris-EDTA buffer and then mixed with 100 µl 0.05 M Tris-3% sodium dodecyl sulfate (SDS). After the mixture was heated at 55°C for 30 min, 100 µl phenol-chloroform was added and emulsified by brief shaking. After being centrifuged at 13,000 rpm for 5 min, 25 µl of the upper aqueous layer was loaded onto a 0.6% agarose gel containing the Escherichia coli transconjugant 39R861 (27), which harbors plasmids of 7.5, 39, 68, and 160 kb, as a size marker. Plasmid molecular masses were calculated by logarithmic comparison to the mass of this marker. The plasmid preparation was repeated for a number of isolates to assess reproducibility.
Genomic DNA extraction. Bacteria were cultured in LB broth or on LB agar overnight at 37°C. Genomic DNA was extracted using either the cetyltrimethylammonium bromide method (Sigma-Aldrich, Ltd., Dorset, United Kingdom) (3) or a QIAGEN tissue DNA extraction kit according to manufacturer's instructions.
Macrorestriction of genomic DNA visualized by PFGE. One milliliter of the overnight culture was resuspended in 1 ml 0.85% NaCl and centrifuged at 13,000 rpm for 5 min. The pellet was resuspended in 0.85% NaCl at a cell density equivalent to a McFarland standard of 3.0. To prepare the plugs, 500 µl 2% chromosomal agarose (Bio-Rad, Hemel Hempstead, United Kingdom) at 55°C was mixed with an equal volume of cells; 500 µl of the mixture was dispensed into plug molds and cooled at 4°C for 20 min. The plugs were placed in 2 ml lysis solution (500 mM EDTA, 1% [wt/vol] sarkosyl [aqueous]) overnight at 55°C, with 40 µl 1-mg/ml proteinase K. The following day, plugs were washed at 55°C three times with water and four times with ET10 (10 mM Tris-Cl, 10 mM EDTA, pH 8); each wash was 15 min long. Plugs were digested with 1 U of XbaI (Roche) overnight at 37°C. The digested plugs were loaded onto a 1% agarose gel and run on a contour-clamped homogeneous electric field DRII system (Bio-Rad) at 14°C with 2 liters of 0.5x TBE (50 mM Tris, 50 mM boric acid, 0.5 mM EDTA). The switch times ranged from 6 to 72 s during a 44-h run at 5.4 V/cm. XbaI-digested DNA from S. enterica serovar Braenderup H9812 was used as a molecular reference marker every five lanes (22).
Statistical analyses of restriction patterns were performed with BioNumerics software (Applied Maths, St. Martens-Latern, Belgium) using the dice similarity coefficient. A dendrogram was calculated by the unweighted pair group method of averages, with a position tolerance of 1%. Fragments smaller than 30 kb were not included in the final analysis, as recommended by the European guidelines for standardization (20). Different XbaI-PFGE profiles were assigned to different PFGE groups according to a single band difference in restriction pattern (11, 20) rather than according to the criteria of Tenover et al. (24).
Salmonella microarray. Salmonella Microarray Generation III is an extension of the previously described Salmonella Microarray Generation I array constructed at the Wellcome Trust Sanger Institute (2, 4, 25). Generation III includes additional coding sequences from the Salmonella genomes being sequenced at The Wellcome Trust Sanger Institute (http://www.sanger.ac.uk/Projects/Salmonella/). Thus, Microarray Generation III is an essentially nonredundant array that contains features representing the following eight genomes: S. enterica serovar Typhi CT18, S. enterica serovar Typhi Ty2, S. enterica serovar Typhimurium LT2 (ATCC 700220) (17), S. enterica serovar Typhimurium DT104 (NCTC 13348), S. enterica serovar Typhimurium SL1344 (NCTC 13347), S. enterica serovar Enteritidis PT4 (NCTC 13349), S. enterica serovar Gallinarum 287/91 (NCTC 13346), and S. bongori 12419 (ATCC 43975).
Hybridization procedures. DNA was extracted from 23 isolates that had been selected to represent PFGE profiles 1 to 10, and the different phage types, resistance patterns, and plasmid profiles from animal and human origins, where available. Genomic DNA was competitively hybridized to Salmonella Microarray Generation III against S. enterica serovar Typhimurium DT104 (NCTC 13348) genomic DNA as the control. Genomic DNA was sonicated (10 s, level 2; Virsonic sonicator) and then fluorescently labeled with Cy5 (test) or Cy3 (control) using the Bioprime kit (Gibco-BRL). Dye-swap labeling experiments also were performed for each test sample. Labeled DNA was purified using an Autoseq G-50 Amersham column, denatured, and precipitated, and the resulting probes were hybridized to the microarray slide for 16 h at 49°C in a hybridization chamber (Genetix X2530). Three slides per isolate were hybridized. Washing procedures were stringent (two washes at 65°C in 2x SSC-0.1% SDS for 30 min and two washes at 65°C in 0.1x SSC for 30 min [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate]).
Microarray analysis. Hybridization to microarray slides was detected by using a Genepix 4000B scanner (Axon Instruments, Inc.) and quantified by using Genepix Pro software (Axon Instruments, Inc.). Signal intensities were corrected by subtracting the local background values. Normalization was performed across all features on the array before any filtering took place. Data were ordered and labeled according to the gene names of the genomes from which each array feature was derived, and BLASTN analysis, using the array feature for DNA sequence as the query, was used to identify the gene in the S. enterica serovar Typhimurium DT104 (NCTC 13348) genome (loci SDT0001 to SDT4571) that was most similar to the tested gene. Signal ratios were analyzed by three methods to establish the cutoff values for designating genes present, absent, and uncertain: (i) twofold variation, which is the default on GeneSpring microarray analysis software V7.2 (Silicon Genetics); (ii) median value (per strain) ± 3 standard deviations (28); and (iii) the gene-calling program GACK, which is based on log2 ratios, with trinary analysis (16). The sensitivity, specificity, and positive and negative predictive values of each method using the E score and bit score were determined for the hybridization results of S. enterica serovar Typhimurium DT104 (NCTC 13348) compared to those of S. enterica serovar Typhimurium LT2. GACK analysis, using mirroring, produced the highest sensitivity, specificity, and positive and negative predictive values, and it was used to assign the present, absent, or uncertain status to each of the 4,167 loci representing the chromosomal features of S. enterica serovar Typhimurium DT104. Information concerning Salmonella Microarray Generation III has been deposited in ArrayExpress as A-SGRP-4, and experimental data are available under the accession number E-SGRP-6 (http://www.ebi.ac.uk/arrayexpress/).
Multiplex PCR. Oligonucleotides for PCR analyses were designed to explore the main regions of variation highlighted by the microarray data (see Table 3). These regions were Salmonella genomic island 1 (SGI1), which is well described for S. enterica serovar Typhimurium DT104 (5, 6, 18), and five prophage-like elements. Multiplex PCR was performed on the 23 isolates to confirm the microarray hybridization results. PCRs contained template DNA (1 to 4 µl), primer pairs (0.5 to 1 µl), and Invitrogen Platinum PCR Supermix (23 to 46 µl). The following PCR protocol was used for all oligonucleotide pairs, except for pairs P8, P9, and P10: 94°C for 1 min; 30 cycles of 94°C for 30 s (denaturation step), 65°C for 1 min (annealing step), and 72°C for 5 min; and 72°C for 2 min. For oligonucleotide pairs P8, P9, and P10, the annealing step was 61°C for 1 min.
|
View this table: [in a new window] |
TABLE 3. Multiplex PCR for detection of the five prophages harbored by S. enterica serovar Typhimurium DT104 (NCTC 13348) and SGI1
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Results of microbiological analysis of S. enterica serovar Typhimurium isolates
|
Nine isolates of animal origin and seven isolates of human origin were plasmid free (Table 2), sensitive to the antibiotics, and phage typed as DT170 (14), DT193 (1), or nontypeable (1). The remaining 45 isolates harbored between one and five plasmids, ranging in size from 2 kb to approximately 141 kb. PFGE 1 isolates harbored an approximately 100-kb plasmid of the same size as the serovar-specific plasmid. An approximately 100-kb plasmid was also detected in two RNDC (reacts but does not conform to designated types on phage typing) isolates, one DT170 type and one DT193 type. The plasmid content could be classified into 15 different profiles.
DNA microarray analysis. Although we could detect variations in the S. enterica serovar Typhimurium isolates, we were unable to genetically define these regions using the above approaches. Consequently, 23 S. enterica serovar Typhimurium isolates, representative of the different PFGE profiles, antibiotic resistance patterns, and phage sensitivity patterns, were subjected to DNA analysis using a composite Salmonella DNA microarray (Table 2). Not surprisingly, the microarray analysis revealed high levels of similarity between the genomes of the different S. enterica serovar Typhimurium isolates (Fig. 1). However, specific regions were variable, with some differences being present in the majority of isolates while others were restricted to particular strains or groups of strains. This limited variation allowed us to perform a simple cluster analysis of the isolates by using a restricted group of 4,167 genes present in S. enterica serovar Typhimurium DT104 (NTCC 13348), for which the complete genome sequence was available (Fig. 1).
![]() View larger version (44K): [in a new window] |
FIG. 1. Analysis of the microarray data by GACK. Test/reference ratios of the 24 S. enterica serovar Typhimurium isolates analyzed by microarray were assessed for the presence/absence of genes by using GACK software. The input data set was restricted to the 4,167 chromosomal features (SDT loci) expected to be present in one or more of the isolates. The heat map is plotted in the physical order of the SDT loci according to the order of the loci in the S. enterica serovar Typhimurium DT104 (NCTC 13348) genome. The presence of a gene in the test strain that is also present in the control strain (S. enterica serovar Typhimurium DT104 [NCTC 13348]) is shown in yellow, the absence or divergence of the gene from its presence in the control strain is shown in blue, and genes for which their presence was uncertain are in gray. The positions of prophage 1/ST104 and prophages 2 to 5 as well as SGI1 are indicated.
|
![]() View larger version (68K): [in a new window] |
FIG. 2. Schematic diagram of prophages 1 to 5 and SGI1, showing oligonucleotide target sites. The gray schematic diagrams represent the predicted open reading frames in the annotated genome of S. enterica serovar Typhimurium DT104 (NCTC 13348). The arrows show the location of the primer pairs P1 to P13, as listed in Table 3. The corresponding microarray results, processed by the gene-calling program GACK, are shown underneath each prophage or under SGI1. The columns represent the features harbored in each prophage or the genomic island. Yellow represents features that are present in the test strain as well as the control strain (S. enterica serovar Typhimurium DT104 [NCTC 13348]), blue represents genes that are in the control strain but are absent in the test strain, and gray indicates that the presence of the gene in the test strain is uncertain. The rows represent each S. enterica serovar Typhimurium isolate, in the following order: row 1, AWF002008; 2, AWX002821; 3, H20002800; 4, H20021651; 5, AWX000826; 6, AWX000841; 7, H20021958; 8, H20004019; 9, AWF009126; 10, AWF009147; 11, H20021856; 12, AWX003658; 13, AWX004816; 14, AWF007581; 15, H19963477; 16, H19992292; 17, AWX004814; 18, H20004467; 19, H20023824; 20 AWX006747; 21, AWX006748; 22, H20003919; 23, H20023530; 24, DT104 (NCTC 13348); 25, LT2 (ATCC 700220).
|
|
View this table: [in a new window] |
TABLE 4. Multiplex PCR results from 23 S. enterica serovar Typhimurium isolates selected for microarray analysis
|
Prophage 2 was found to be present in all 23 isolates probed with oligonucleotide pair P4, and hence GIFSY-2 phage-related DNA is conserved in this population. In contrast, prophage 3, which was probed for by using oligonucleotide pairs P5, P6, and P7, was present only in the seven PFGE 1 isolates.
Prophage 4 is present in the sequenced DT104 isolate, but interestingly it was not detected in any of the 23 other S. enterica serovar Typhimurium isolates, including other DT104 isolates, by microarray analysis. The presence of prophage 4-related sequences was probed for by using oligonucleotide pair P8, and interestingly a PCR product of the predicted size was produced by using DNA prepared from all but three isolates (H20004467, AWX006747, and H20003919), suggesting that these isolates encode sequences relating to prophage 4.
The PCR results obtained with prophage 5-based probe pairs corresponded to the information gleaned from the microarrays. DNA sequences related to prophage 5, probed by P9 and P10, were present in all the PFGE 1 isolates and the three RDNC isolates belonging to other PFGE types. All other S. enterica serovar Typhimurium isolates generated a product with P9 but not P10, suggesting that they harbored a variant form of prophage 5.
Thus, these sets of primers confirm aspects of the microarray analysis of variation in specific prophage regions of S. enterica serovar Typhimurium and may form the basis of a simple PCR assay for discriminating between different S. enterica serovar Typhimurium strains circulating in a geographical area.
|
|
|---|
Using the recently determined genome sequence of the DT104 isolate NCTC 13348, we were able to identify up to five different regions of the genome encoding prophage-like elements, complementing and extending studies carried out by others using non-genome-wide approaches (12, 13). Prophage 1/ST104, also known as PDT17 (21), previously has been described as a resident prophage in S. enterica serovar Typhimurium DT104 (23). Interestingly, a prophage 1-like element also was detected in the nontypeable S. enterica serovar Typhimurium isolates H20004467 and H19992292 (DT170), although only the former was confirmed by multiplex PCR. These observations are of interest, as prophage 1/ST104 previously has not been described outside of phage type DT104. Prophage 2 (GIFSY-2), which is a well-characterized S. enterica serovar Typhimurium prophage, was present in all 23 isolates and was thus of little value for discriminating between strains. Prophage 3 was present in the seven PFGE 1 isolates and, together with SGI1, could form the basis of a diagnostic test for DT104.
Significantly, although these techniques could distinguish between isolates of S. enterica serovar Typhimurium, organisms derived from humans or animals did not fall into two mutually exclusive groups. This observation serves to illustrate the close genetic relationship between animal and human Salmonella strains circulating in the same geographical region and emphasizes the zoonotic nature of these infections. It seems that very similar strains of S. enterica serovar Typhimurium were causing disease in animals and humans in the same geographical location and in the same time frame. Antibiotic-resistant strains of S. enterica serovar Typhimurium were present in both the human and animal populations. Interestingly, antibiotic-sensitive isolates from animals tended to cluster with sensitive isolates from humans, and antibiotic-resistant isolates also grouped together, irrespective of source.
In conclusion, the examination of the DNA sequences of the different prophages present in S. enterica serovar Typhimurium DT104 (NCTC 13348) allowed us to design PCR primer sets that could be used to interrogate the genomes of different field isolates in multiplex PCR assays. These simple tests were used to rapidly distinguish between different S. enterica serovar Typhimurium isolates by using basic equipment and reagents. Thus, we have developed a rational DNA sequence-based assay for discriminating between S. enterica serovar Typhimurium isolates that can be applied even in a simple laboratory environment. Importantly, we believe that such assays could be refined to form the basis of a simple rapid test to distinguish between S. enterica serovar Typhimurium isolates in a diagnostic or surveillance setting. Further sequence analyses are under way to define the prophage content of multiple S. enterica serovar Typhimurium isolates, with a view to improving the discriminatory power of these assays. This approach easily could be expanded to the design of similar PCR-based assays for additional Salmonella serovars. The advantage of this approach is that it is based on defined DNA sequences physically mapped in the Salmonella genomes, and it targets known hot spots for evolutionary change. Hence, these signatures reasonably might be expected to detect the early emergence of pathogenic variants of Salmonella in the field.
We thank Matloob Qureshi for his help with submission of the microarray data to ArrayExpress.
Published ahead of print on 23 May 2007. ![]()
|
|
|---|
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»