Previous Article | Next Article ![]()
Journal of Clinical Microbiology, August 2005, p. 4092-4106, Vol. 43, No. 8
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.8.4092-4106.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Institute of Microbiology, ETH Zürich, Wolfgang-Pauli-Str. 10, 8093 Zurich, Switzerland,1 Robert Koch Institut, 38855 Wernigerode, Germany,2 MicroDiscovery GmbH, Marienburger Strasse 1, 10405 Berlin, Germany3
Received 5 January 2005/ Returned for modification 15 March 2005/ Accepted 14 April 2005
|
|
|---|
|
|
|---|
For diagnostic purposes, Salmonella enterica subspecies 1 serovars are differentiated by serotyping (28), which allows reliable identification of organisms that cause infections associated with higher risks (i.e., serovar Typhi). For epidemiological surveys, detection of outbreaks, and uncovering of the chains of transmission, serotyping is insufficient because serovar Enteritidis and serovar Typhimurium often account for >80% of all human salmonelloses. For further differentiation, serotyping is complemented by phage typing as well as molecular typing techniques.
During the past decades serovar Typhimurium strains have been differentiated by phage typing (1). The Anderson phage typing scheme is based on the empirical observation that certain phages can cause lytic infections in one serovar Typhimurium strain but not another. A set of 34 bacteriophages allows more than 200 patterns of resistance or sensitivity, the "phage types" (definitive phage types [DTs]), to be distinguished. Phage typing has proven very useful for distinguishing many serovar Typhimurium strains. It has been invaluable for surveillance of the emergence and spread of epidemic serovar Typhimurium strains (i.e., DT204c and DT104) during the past decades (23, 40, 42-44, 49).
In spite of its success, phage typing has several drawbacks. (i) Phage typing has remained an empirical method, and there is only circumstantial evidence for the molecular mechanisms determining resistance or susceptibility to each Anderson typing phage. (ii) The power of resolution is limited. Several distinct serovar Typhimurium strains can belong to the same phage type. (iii) Phage typing relies on "unique" historical stocks which will not last forever. Replenishment of phage stocks will require propagation in Salmonella strains. Recent research has demonstrated that this can lead to recombination and mixed phage populations. Thus, the resulting phage preparations can show altered lysis patterns (36) and, consequently, can affect the consistency of the Anderson phage typing scheme.
Several molecular typing methods have been devised to differentiate serovar Typhimurium strains. This includes pulsed-field gel electrophoresis (PFGE), ribotyping, random amplification of polymorphic DNA, DNA fingerprinting, amplified length polymorphism analysis, enterobacterial repetitive intragenic consensus sequence typing, determination of outer membrane protein and multilocus enzyme patterns, and IS200 typing. Currently, these methods or combinations thereof are being employed successfully for outbreak investigations and identification of lines of transmission. In the past few years, the DNA microarray technology has been developed and is a potentially powerful tool for strain identification. DNA microarrays harboring probes for each open reading frame (ORF) of serovar Typhimurium strain LT2 (GenBank accession no. NC_003197) and additional probes derived from the serovar Typhi CT18 genome (GenBank accession no. NC_003198) and other Salmonella sequences from public databases have been designed (2, 15, 30). Hybridization of chromosomal DNA from standard Salmonella enterica laboratory collections allowed discrimination of a total of >100 strains of different S. enterica subspecies, various subspecies 1 serovars, and even different strains of the same serovar (29). However, so far no systematic studies have evaluated the suitability of DNA microarrays for the differentiation of serovar Typhimurium strains, e.g., for purposes of epidemiology.
Here, we describe a pilot study in which we assessed DNA microarrays as a tool for serovar Typhimurium strain differentiation. A DNA microarray was designed by using strain-specific probes obtained by differential subtractive hybridization and probes derived from public databases. Repeated hybridization experiments with a set of reference serovar Typhimurium strains demonstrated a high reproducibility of strain differentiation by this technique. Furthermore, we analyzed clinical serovar Typhimurium isolates of phage type DT104 and additional phage types. We found that the DNA microarray is useful for serovar Typhimurium strain differentiation. Implications for future applications and strategies for further optimization of the microarray are discussed.
|
|
|---|
Serovar Typhimurium isolates and serovar Typhi DNA (SARC2) were from the Robert Koch Institut, Wernigerode, Germany (Table 1). Chromosomal DNAs from Escherichia coli strain EDL933 and Yersinia enterocolitica strain 8081 were kindly provided by A. Rakin and S. Schubert (Max von Pettenkofer-Institut, Munich, Germany).
|
View this table: [in a new window] |
TABLE 1. Strains
used in microarray design and hybridization
|
PFGE was performed as described previously (31).
SH and subcloning of the DNA fragments. Subtractive hybridization (SH) was performed by using the Clontech PCR-Select bacterial genome subtraction kit (PT3170-1; Clontech, Palo Alto, Calif.). Briefly, chromosomal serovar Typhimurium DNA (tester DNA) was digested with RsaI (four-base cutter, GTAC), purified, and divided into two aliquots. One aliquot was ligated to adapter 1, and the other was ligated to adapter 2R.
In a third tube, portions of each aliquot were combined (negative control for subtractive hybridization). Two rounds of hybridization with RsaI-digested DNA from another serovar Typhimurium strain (driver DNA) were performed at a temperature (48°C) lower than that recommended by the manufacturer (63°C) in order to increase the fraction of DNA fragments with low or no similarity between the driver and the tester DNAs. Then, a primary PCR and a secondary PCR with primers that anneal to adapters 1 and 2R were performed, and the PCR products were cloned into pCR2.1-Topo (Topo TA cloning kit; Invitrogen, Germany). Plasmids from lac-negative clones growing in Luria-Bertani medium (100 µg/ml ampicillin, 50 µg/ml kanamycin, isopropyl-ß-D-thiogalactopyranoside, 5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside) were analyzed by PCR (with primers T7 [5'-GTAATACGACTCACTATAGGGC-3'] and M13 reverse [5'-GGAAACAGCTATGACCATG-3']; annealing at 52°C and elongation for 1.5 min at 72°C for 35 cycles) and sequenced (AGOWA, Berlin, Germany). As probes for the microarray, we selected 500- to 800-nucleotide (nt)-long inserts (named A2 to G4 and 2SH4 to 3SH93; Table 2) showing <75% identity (blastn;www.ncbi.nlm.nih.gov/BLAST) over >50% of their entire sequence to serovar Typhimurium LT2, E. coli K-12, and Shigella genome sequences (http://www.ncbi.nlm.nih.gov/sutils/genom_table.cgi). Inserts with high similarities to mobile genetic elements (transposons, insertion sequence elements, phages) or pathogenicity islands were also included.
|
View this table: [in a new window] |
TABLE 2. Microarray
probes
|
From each strain to be analyzed, 4 µg of chromosomal DNA was denatured, hybridized to random decamer primers, and labeled by primer extension by using Klenow polymerase and a deoxynucleoside triphosphate mix containing Cy5-dCTP. The labeled DNA was purified (MiniElute; QIAGEN), and the labeling efficiency was verified by gel electrophoresis. Array hybridization was performed in a formamide-based buffer (Sci-Hyb; Scienion) at 42°C for 20 h, followed by washing steps with increasing stringency (Sci-Wash; Scienion). The fluorescence signals (635-nm excitation) obtained with each probe were measured with an Axon 4100 scanner (16-bit; intensity range, 0 to 65,535), and the raw data were corrected for the background and processed with GenePix Pro V5.0 software. This was performed by Scienion.
Data processing. The data were transformed to log scale and normalized by quantile-based methods (37, 50). A two-component normal mixture model (24), p(x; m1, m2, s1, s2, w) = (1 w) · n(x; m1, s1) + w · n2(x; m2, s2), was fitted to the normalized data by a maximum-likelihood method. The discriminant function (12) w · [n2(x)/p(x)] was used to represent the propensity of the sequence for being present (value close to 1) or absent (value close to 0) from the DNA of that strain. Discriminant values for all probes and all samples (strains) were calculated and stored in a signal probability matrix, with each column representing one array measurement and with each row representing the results obtained with one specific probe on the chip. For representation of the probability matrix, a coloring scheme was used: blue indicates state 1 (off), yellow indicates state 2 (on), and dark shading indicates indecisive measurements (p value of about 0.5).
Strain differentiation was studied as follows: pairwise correlation coefficients (16) were calculated from the signal probability values from each pair of hybridization experiments. The pairwise correlation coefficients were stored in a sample correlation matrix. For visualization of the results, the sample correlation matrix was displayed as a color scale ranging from green (weakly correlated) to red (strongly correlated) (adapted from the "sma" package [13]).
Dendrograms were calculated based on the signal probability matrix by using Euclidean distance and the average linkage method (12, 38). All calculations were implemented by using the R-statistics package (22).
|
|
|---|
DNA microarray probes obtained by SH. SH allows the enrichment or retrieval of sequences present in one bacterial genome (designated the "tester") but absent from another (designated the "driver") (see Materials and Methods). These sequences would provide ideal probes for DNA microarray design.
The DNA of
isolates of three serovar Typhimurium strains (DT204 =
STM1690/75, DT104 = 99-00971, and DT208 = 201 IE 888;
Table 1) known to differ
significantly from each other (W. Rabsch, R. Prager, and H.
Tschäpe, unpublished observations) were chosen as testers. DT204
had caused a major epidemic
(23,
42,
44) in the 1970s and
1980s and carries the SopE
prophage, which encodes the type
III effector gene sopE, in its chromosome
(26,
27). Serovar Typhimurium
DT104 (99-00971) is a representative of the antibiotic multiresistant
phage type DT104 strains which emerged in several European countries as
the predominant serovar Typhimurium phage type in human infections
(31,
43,
49). Outbreaks have also
been reported in the United States
(5,
6,
34). Serovar Typhimurium
DT208 (201IE888) is a representative of phage type DT208, which caused
an outbreak in the former Soviet Union but which is also identified in
current epidemiological surveys
(11,
41).
We used DNA from the serovar Typhimurium strains mentioned above and from DT36 (Anderson), a strain highly susceptible to all Anderson typing phages. Tester-driver combinations were DT204 (STM1690/75)-DT104 (99-00971), DT208 (201IE888)-DT204 (STM1690/75), DT104 (99-00971)-DT36 (Anderson), DT204 (STM1690/75)-DT36 (Anderson), and DT208 (201IE888)-DT36 (Anderson). The hybridization steps were performed at 48°C to ensure enrichment of the tester DNA fragments with no significant similarity to the driver genome. DNA fragments were sequenced, and we selected 43 plasmids harboring 500- to 1,600-nt DNA inserts with no or limited similarity to known Salmonella or enterobacterial sequences for probe design (see Materials and Methods). Inserts with high similarities to mobile genetic elements (transposons, insertion sequence elements, phages) or pathogenicity islands were also included. Upon PCR amplification, 38 of these plasmids yielded unique products and were selected for DNA microarray design (Table 2; see Table S1 in the supplemental material). Nine of these probes originated from DT104 (99-00971), 24 originated from DT204 (STM1690/75), and 5 originated from DT208 (201IE888) (Table 2).
DNA microarray probes derived from prophages. The repertoire of prophages residing in the genome is highly diverse between different serovar Typhimurium strains (8, 17, 26). Unique sequences (immunity region, morons) which are present in one prophage but not in most others should be well suited for strain differentiation. PCR products corresponding to 21 unique prophage sequences from serovar Typhimurium strains, serovar Typhi, and serovar Paratyphi A (33) were selected for DNA microarray design (Table 2; see Table S2 in the supplemental material). Three probes were derived from the phage-encoded Shiga-like toxin A subunit (E. coli EDL 933; O157:H7) and cholera toxin A and B subunits (Vibrio cholerae El Tor Nent 720-95). It was unclear whether one of these toxins might also be present in some serovar Typhimurium strains.
DNA microarray probes from TTSSs. Type III secretion systems (TTSSs), their known effector genes, and chaperones are prerequisites for Salmonella virulence (18). The TTSSs and most of the effector genes are highly conserved between all Salmonella strains (25, 30). We designed probes based on sopE2, sopB, sipA, slrP, invB, invC, and prgH as positive controls (Table 2; see Table S2 in the supplemental material).
Effector genes sopE (probe D2 from SH) and avrA are more variable (20, 21, 26, 32) and were also included.
Additional probes for DNA microarray design. Additional probes were based on serovar Typhimurium SL1344 and DT204 (3341/78) DNA fragments retrieved in an earlier study (Table 2, probes pM48, K58, K51, and K9) (K. Ehrbar and W.-D. Hardt, unpublished data) and fimC (serovar Typhimurium LT2). Furthermore, open reading frames STM2240 (serovar Typhimurium LT2) and STY0312, STY1360, and STY1364 (serovar Typhi), which do not show similarity to other sequences in the public database, were selected for probe design (Table 2; see Table S2 in the supplemental material).
In an attempt to control the stringency of microarray hybridization probes E. coli frd (GenBank accession no.DH5_FRB; 90% identity to serovar Typhimurium LT2), E. coli kefc (GenBank accession no. DH5_KEFC; 83% identity to serovar Typhimurium LT2), E. coli agp (GenBank accession no. DH5_AGP; 80% identity to serovar Typhimurium LT2), and Yersinia enterocolitica thrc (GenBank accession no. Y_ENT_THRC; 76% identity to serovar Typhimurium LT2) (Table 2) were designed. The green fluorescent protein (gfp) (10) served as a negative control.
Microarray design. For microarray design, a total of 83 probes (250 to 1,500 nt long) were used, including 38 probes derived from subtractive hybridizations, 4 from earlier genetic screens, and 41 from the sequences in a public database (National Center for Biotechnology Information). Each probe was spotted onto two different locations of the amino slides. Probes E10, pM816, D6, and G3 each were present in four different spots on the microarray. A total of 100 microarrays were produced in two separate batches. Quality control (see Material and Methods) and hybridization experiments with identical DNA samples did not detect differences in the performance of the two different batches (data not shown).
Hybridization with five serovar Typhimurium reference strains. We performed repeated DNA microarray hybridization experiments to analyze the reproducibility and explore the quality of strain differentiation obtained by this method. For this purpose serovar Typhimurium DT36 (Anderson), DT204 (STM1690/75), DT104 (99-00971), DT208 (201IE888), and LT2 were chosen. The chromosomal DNA of one strain was isolated in two independent preparations, and each batch was hybridized with three DNA microarrays (hybridization data not shown). Each of the six hybridizations involved a separate labeling reaction. In this way we could determine whether technical problems associated with poorly reproducible steps of the hybridization procedure (see Materials and Methods) might interfere with DNA microarray-based strain typing.
Hybridization signals were transformed to log scale and normalized (see Materials and Methods). The signals obtained from pairs of identical probes on each array were very similar and showed typical correlation coefficients greater than 0.99 (data not shown). For all further calculations, the two signals were represented by their mean value.
Next, we used the hybridization signals to calculate the probability that a specific sequence is present in or absent from the strain. For this purpose, a two-component normal mixture model was fitted to the hybridization data by a maximum-likelihood method (see Materials and Methods). This yielded a "discriminant" function (see Fig. S1 in the supplemental material) which can use the hybridization signals to calculate for a probe the probability of a match (value close to 1) or no match (value close to 0) with the strain analyzed. The discriminant function was used to transform the 83 hybridization signals into a list of probabilities ("signal probability matrix") for the presence or absence of each probe (data not shown). To study strain differentiation, the signal probability list between each pair of DNA microarray hybridizations was compared. The resulting correlation coefficients were stored in a sample correlation matrix (see Materials and Methods). In Fig. 1, the sample correlation matrix for the analysis of the five representative serovar Typhimurium strains has been displayed in color code: red indicates strongly correlated, and green indicates weakly correlated.
![]() View larger version (60K): [in a new window] |
FIG. 1. Correlation
analysis of the five reference serovar Typhimurium strains. To study
strain differentiation, the signal probability list was compared
between each pair of DNA microarray hybridizations. The resulting
correlation coefficients were stored in a sample correlation matrix
(see Materials and Methods) (data not shown). Red, strong correlation;
green, weak correlation. Each reference strain was hybridized six times
with the microarray. The reference strains are color coded: DT36
(Anderson), marine; DT204 (STM1690/75), blue; DT104 (99-00971), grey;
DT208 (201IE888), pink; and LT2,
yellow.
|
Slight variations were observed within the group of six DNA microarray hybridizations performed with serovar Typhimurium DT36 (Anderson), as indicated by a shift to dark red in the correlation matrix. The same observation was made for serovar Typhimurium LT2. The reasons for this variation are unclear. However, strain differentiation was still clearly feasible.
The weakest correlation was observed between DT204 and the other strains, in particular, strains DT104 (99-00971) and DT208 (201IE888). This is at least partially attributable to the large number of DT204 (STM1690/75)-specific probes (24 in total; Table 2) present on the DNA microarray. In addition, DT204 (STM1690/75) failed to hybridize with many probes specific for the other strains (data not shown).
The hybridization patterns were also compared in a dendrogram. For this purpose, the signal probability matrix was analyzed by using Euclidean distance and the average linkage method (see Materials and Methods). The topology of the resulting dendrogram (Fig. 2) supported the results from the correlation analysis (Fig. 1). The branches from the six hybridizations of each serovar Typhimurium strain clustered closely together. The branch representing DT204 (STM1690/75) diverges early, while the branches for DT104 (99-00971) and DT208 (201IE888) cluster closely together.
![]() View larger version (28K): [in a new window] |
FIG. 2. Dendrogram
for hybridization data of the reference strains serovar Typhimurium
DT36 (Anderson), DT204 (STM1690/75), DT104 (99-00971), DT208
(201IE888), and LT2. The signal probability matrix was analyzed by
using Euclidean distance and average linkage (see Materials and
Methods).
|
Taken together the results from this first series of hybridization experiments provided evidence for the reproducibility of the data and the feasibility of DNA microarray-based differentiation of serovar Typhimurium strains.
Microarray hybridization with various DT104 typed strains. The data presented above showed that the simple prototype DNA microarray could differentiate serovar Typhimurium strains which belong to different phage types. Next, we were interested in finding out whether different serovar Typhimurium isolates of the same phage type would yield similar hybridization results. To analyze this, we have chosen 13 serovar Typhimurium DT104 strains (Tables 3 and 4, strains 02-00822, 02-02682, 02-03881, 03-06256, 97-01178, 97-01193, 97-01864, 97-02103, 97-01174, 02-00429, 02-03122, 96-01186, and 97-05652). Three independent DNA microarray hybridizations (I, II, III) were performed with each strain, as described above. The DT104 (99-00971) reference strain was also hybridized once more.
|
View this table: [in a new window] |
TABLE 3. Strains
studied by microarray
hybridizationa
|
|
View this table: [in a new window] |
TABLE 4. Characteristics
of DT104 serovar Typhimurium strains
|
![]() View larger version (20K): [in a new window] |
FIG. 3. Dendrogram
of 13 serovar Typhimurium DT104 strains (calculated as described in the
legend to Fig. 2). The
five reference strains from Fig.
1 and
2 were included for
comparison (DT36 [Anderson], marine; DT204 [STM1690/75], blue; DT104
[99-00971], grey; DT208 [201IE888], pink; and LT2, yellow). The signal
probability matrix was analyzed by using Euclidean distance and average
linkage (see Materials and Methods). Subgroup 1 strains, 02-00822,
02-02682, 02-03881, 03-06256, 97-01178, 97-01193, 97-01864, 97-02103,
and 97-01174 (similar to reference strain DT104 99-00971); subgroup 2
strains, 02-00429, 02-03122, 96-01186, and
97-05652.
|
The correlation analysis confirmed the finding of two DT104 subgroups (Fig. 4).
![]() View larger version (70K): [in a new window] |
FIG. 4. Correlation
analysis of the 13 serovar Typhimurium DT104 strains. Three
hybridizations were performed per DT104 strain. The signal probability
list was compared between each pair of DNA microarray hybridizations.
Data for the five reference strains were included for comparison (DT36
[Anderson], marine; DT204 [STM1690/75], blue; DT104 [99-00971], grey;
DT208 [201IE888], pink; and LT2, yellow). The resulting correlation
coefficients were stored in a sample correlation matrix (data not
shown) and displayed as red (strong correlation) and green (weak
correlation). The analysis clearly identified two distinct DT104
subgroups.
|
![]() View larger version (50K): [in a new window] |
FIG. 5. PFGE
patterns of 13 serovar Typhimurium strains defined as DT104.
Salmonella serovar Braenderup was used as a reference (lanes
1, 9, and 16). Chromosomal DNA was digested with XbaI and analyzed by
PFGE (see Materials and Methods). Serovar Typhimurium DT104 strains:
lane 2, 96-01186; lane 3, 02-00429; lane 4, 97-05652 lane 5, 02-03122;
lane 6, 97-01147; lane 7, 97-02103; lane 8, 97-01193; lane 10,
02-03881, lane 11, 97-01864; lane 12, 97-01178; lane 13, 02-00822; lane
14, 02-02682, lane 15,
03-06256.
|
Microarray hybridization with various phage-typed strains. So far, our study had focused on the differentiation of strains representing five different phage types. However, the number of serovar Typhimurium strains in nature is much higher. Even the Anderson phage typing scheme is able to differentiate >200 different types of serovar Typhimurium strains. In order to explore the potential and the limitations of the current prototype microarray, we have studied strains representing various phage types. Here, each strain was hybridized only once. Therefore, the data can provide general insights, but they are not statistically significant. Strains belonging to phage types DT204 (00-02379, 00-07509, 00-09603, 00-10022, and 99-07382), DT193 (02-01414, 02-04395, 02-04515), DT120 (03-01606, 03-01663, and 03-01860), DT92 (03-00009, 03-0008), DT12 (03-6841, 03-07846), DT17 (03-01862), DT177 (03-05510), DT186 (03-5115) DTu302 (03-6382), DT7 (03-6600), DT9 (03-01779), DT170 (04-0093), DT2 (04-0110), DT175 (305-70), DT68 (11635), and DT40 (6772/96) were analyzed as described above. For comparison, data for the five reference strains (Fig. 1 and 2) have been included (color coded in blue, grey, pink, marine, and yellow). As expected, all five reference strains were located in distinct branches of the dendrogram (Fig. 6). The vast majority of the other serovar Typhimurium strains clustered in three distinct branches of the dendrogram (termed X, Y, and Z; Fig. 6). These branches were clearly distinct from those of the five reference serovar Typhimurium strains. The hybridization patterns of two strains, DT193 (02-01414) and DT177 (03-05510), differed from those of all of the other strains analyzed. This indicated that the prototype DNA microarray can differentiate at least 11 different types of serovar Typhimurium strains, the five reference strains (Fig. 1); DT104 subgroup 2 (Fig. 3 and 4); clusters X, Y, and Z; and the two strains DT193 (02-01414) and DT177 (03-05510).
![]() View larger version (26K): [in a new window] |
FIG. 6. Dendrogram
of various serovar Typhimurium strains determined by Euclidean distance
and average linkage (see Materials and Methods). The five reference
strains were included for comparison (Fig.
1 and
2; color coded in blue,
grey, pink, marine, and yellow). Strain designations are indicated on
the left. Branch X, DT204 (00-02379, 00-07509, 00-09603, 00-10022,
99-07382) and DT68 (11635); branch Y, DT17 (03-01862), DT2 (04-0110),
DT193 (02-04395), DT12 (03-07846, 03-6841), DT175 (305-70), DT92
(03-00009, 03-0008), DT186 (03-5115), and DT170 (04-0093); branch Z,
DT193 (02-04515), DT9 (03-01779), DT120 (03-01663, 03-01606, 03-01860),
DT7 (03-6600), and DT40
(6772/96).
|
Interestingly, the DT68 (11635)
strain, which carries the SopE
prophage
(26), was also located in
the same branch as the DT204 strains. Similarly, 03-6382 (Dtu302)
clustered closely with the DT104 (99-00971) reference strain. Further
studies are needed to determine whether strains DT68 (11635) and DT204
or strains Dtu302 (03-6382) and DT104 may have a recent common
ancestry.
Cluster Y harbored both DT92 isolates and both DT12 isolates. The three DT120 isolates were located in cluster Z (Fig. 6). Thus, strains of these three phage types might be recognized specifically by the prototype DNA microarray and differentiated from the five reference serovar Typhimurium strains. Further experiments are needed to substantiate this hypothesis.
The three DT193 isolates (02-01414, 02-04395, and 02-04515) were located in separate branches of the dendrogram, implying that this DT group includes heterogeneous serovar Typhimurium strains which can be differentiated by the DNA microarray. Overall, clusters Y and Z harbored quite diverse groups of strains, as judged by phage typing. It is likely that the prototype DNA microarray is simply insufficient to differentiate between these strains. In this case, additional probes would be required in order to improve strain differentiation by future generations of a serovar Typhimurium typing microarray.
Evaluation of the probes chosen for microarray design. For the future design of improved serovar Typhimurium typing arrays, it would be useful to analyze in more detail the performance of each probe on the current DNA microarray. For this purpose, the signal probability matrix (see Materials and Methods) for each probe in each hybridization experiment has been evaluated (see Fig. S2 in the supplemental material). Twelve DNA probes matched every serovar Typhimurium strain analyzed. This included most probes derived from TTSSs as well as probes E6 (transposase), A7 (prophage ST64T), A8 (transposase), A6 (SPI3), and B4 (100% identity to LT2) obtained by subtractive hybridization. Twelve DNA probes did not match any of the serovar Typhimurium strains analyzed. This included the negative control (GFP); prophage-encoded genes for CTXA, CTXB, and EDL-STXA; the serovar Paratyphi B sequences Paratyphi-Contig 8 and Paratyphi-Contig 35; the serovar Typhi sequences Typhi-STY1360, Typhi-STY0312, and Typhi-STY4629; Y. enterocolitica THRC; DH5-AGP; and DH5-KEFC. Some of these probes might be omitted in future. However, a few "positive and negative controls" should be retained. It might be of interest to retain probes GFP (truly negative control), Y_ENT_THRC (76% identity to LT2), DH5_AGP (80% identity to LT2), and DH5_KEFC (83% identity to LT2). In contrast to these "absent" probes, control probe DH5_FRB (90% identity) yielded matches in almost all hybridization experiments. Thus, these probes can provide a useful control for the general stringency of DNA microarray hybridization. In addition, these observations indicate that DNA fragments less than 83% identical between two serovar Typhimurium strains can be employed on future DNA microarrays.
Several probes derived from the serovar Typhimurium LT2 (10 probes) and serovar Typhi CT18 (6 probes) genome sequences contributed to strain differentiation. This pertained especially to probes derived from unique prophage sequences. Thus, additional completed Salmonella genomes will provide promising sources for the design of additional probes.
The vast majority (30 of 38) of probes derived by subtractive hybridization clearly contributed to strain differentiation. Due to the possibility to experimentally enrich DNA fragments that differ between two strains of choice, subtractive hybridization will provide a powerful tool by which additional probes for improved DNA microarray-based strain differentiation can be obtained.
|
|
|---|
Our study included merely 44 serovar Typhimurium strains. This included the most prevalent serovar Typhimurium strains from Germany, but it is just a small fraction of the existing serovar Typhimurium strains with distinct genetic makeups (>>200). Thus, it is conceivable that the prototype DNA microarray can differentiate even more serovar Typhimurium strains than is evident from this study. More hybridization experiments would be required to address this point.
Future generations of the DNA microarray should include additional probes that will allow the accurate discrimination between those strains not differentiated by the current DNA microarray (i.e., subgroup 1 DT104 strains [Fig. 5, lanes 6 to 8 and 10] and strains of clusters Y and Z [Fig. 6]) or strains which are going to emerge in the future.
How many strain-specific probes are required? Optimal strain differentiation was achieved for the reference strain DT204 (STM1690/75; Fig. 1). This strain yielded strong correlations in repeated hybridizations but was clearly distinct from other DT204 strains (Fig. 6) and yielded only a weak correlation with all other reference strains (Fig. 1). The DNA microarray harbored 31 probes that yielded positive hybridization signals with DT204 (STM1690/75) but not with one or more of the other strains. Twenty-one of these "useful" probes were derived from DT204 (STM1690/75) by subtractive hybridization. This number of probes represents an estimate for the upper limit, since all differences are expected to show up in this experimental setup. Sensible strain differentiation can probably be achieved with fewer strain-specific probes. Reference strain DT104 (99-00971) provides an example: only six of the nine probes derived from DT104 (99-00971) by subtractive hybridization yielded positive hybridization signals with DT104 (99-00971) but not with one or more of the other strains. Still, the prototype microarray could differentiate DT104 (99-00971) from all other reference strains and could even distinguish between two different subgroups of DT104 strains (Fig. 4).
How can the prototype DNA microarray be improved? Clearly, additional DNA probes that are present in some serovar Typhimurium strains but absent from others are needed. As discussed above, DNA fragments of several hundred base pairs <83% identical between different strains will be useful for strain typing. It should be noted that (in contrast to many other DNA array studies) it has proven useful to focus on DNA similarity and not on open reading frames. Thus, many probes obtained by subtractive hybridization represent partial open reading frames or noncoding DNA.
In theory, optimal efficiency would be achieved if each of the probes would recognize a different subset of serovar Typhimurium strains. More precisely, an optimal (orthogonal) set of probes would show a vanishing correlation for each pair of probes or samples. Currently, there are two main strategies for acquisition of additional probes: ongoing sequencing projects which include several serovar Typhimurium strains and strains from other S. enterica subspecies 1 serovars will provide rich sources for probe generation. Alternatively, subtractive hybridization can be used. This method has provided many useful probes in this study. It is fairly inexpensive and allows the enrichment of DNA sequences present in one strain but absent from another. Therefore, it is ideally suited for the identification of probes which distinguish between strains which are not clearly distinguished by the prototype DNA microarray (i.e., different DT104 strains of subgroup 1 [Fig. 4 and 5] and different strains of clusters Y and Z [Fig. 6]).
In conclusion, this pilot study has demonstrated that DNA microarrays are useful tools for serovar Typhimurium strain differentiation. Most likely, strains from other S. enterica subspecies 1 serovars can be differentiated in a similar manner. If one considers that the technology for chip production, hybridization, and analysis has become widely available in recent times, DNA microarray strain typing holds many promises for the fast and efficient detection of chains of transmission and for epidemiology in general. However, it should be noted that today's high costs for DNA microarray production and the labeling must drop significantly before this technology will be widely applied in diagnostic and reference laboratories.
This work was funded in part by a grant from the Kompetenznetzwerk "Genomforschung an pathogenen Bakterien" Verbund 3 (BMBF, Germany) to W.-D. Hardt.
Supplemental material for this article may be found at http://jcm.asm.org. ![]()
|
|
|---|
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»