Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2007, p. 2624-2630, Vol. 73, No. 8
0099-2240/07/$08.00+0 doi:10.1128/AEM.02465-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Division of Bacteriology, Chiba Prefectural Institute of Public Health, Chiba,1 Laboratory of Veterinary Public Health, College of Bioresource Science, Nihon University, Kanagawa,2 Kanagawa Prefectural Institute of Public Health, Kanagawa,3 Department of Microbiology and Molecular Genetics, Graduate School of Pharmaceutical Science, Chiba University, Chiba, Japan4
Received 22 October 2006/ Accepted 30 January 2007
|
|
|---|
|
|
|---|
Rodents are public health hazards since they can be reservoirs for human salmonellosis (9, 12, 26, 30-32). We have previously reported the prevalence of Salmonella enterica in Rattus norvegicus and Rattus rattus rats captured at two buildings located in Yokohama City, Japan, during 1997 and 1998 (25). In that study, 60 of 339 (17.7%) house rats carried S. enterica, and of these 60 rats, 58 (96.7%) carried serovar Typhimurium and two carried serovar Enteritidis, pointing to the importance of rats as reservoirs of S. enterica serovar Typhimurium.
While there are some reports about salmonellae in rats (9, 12, 25, 26, 30-32), no information has been published regarding the incidence of DT104. The lack of information about DT104 is partially due to the difficulty of performing definitive phage typing. Although phage typing is the "gold standard" for DT104 identification, PCR has been used recently to identify DT104-specific DNA sequences (11, 28) and could facilitate studies of the DT104 prevalence in animal reservoirs. However, more comparative data on DT104 identification by PCR are needed.
A major problem for the investigation of DT104 prevalence in animal reservoirs is the experimental difficulty of determining genetic relationships among DT104 isolates. For example, molecular epidemiological analysis using pulsed-field gel electrophoresis (PFGE) analysis cannot readily determine genetic relationships (1, 4, 23). Recently, variable-number tandem repeat (VNTR) typing was reported to be a new molecular epidemiological method for differentiating pathogenic bacteria (21). Furthermore, VNTR typing has been used to differentiate DT104 strains having the same PFGE patterns (20).
In this study, we have investigated both the prevalence and the genetic properties of DT104 strains among serovar Typhimurium isolates from house rats by studies of antimicrobial resistance, detection of class I integrons, and antimicrobial resistance genes. Molecular epidemiological analysis was done using PFGE and VNTR typing.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Antimicrobial resistance patterns of serovar Typhimurium DT104 isolates from house rats in building J
|
|
View this table: [in a new window] |
TABLE 2. Antimicrobial resistance patterns of serovar Typhimurium DT104 isolates from house rats in building YS
|
Characterization of DT104 strains.
Antimicrobial sensitivities of the DT104 strains were determined by the Kirby-Bauer single disk method (2). Thirteen different Sensi-Discs (Becton Dickinson, Franklin Lakes, NJ) were used, and the manufacturer's instructions were followed for testing. Each disk contained one of the following: 10 µg ampicillin, 30 µg cephalothin, 30 µg cefotaxime, 30 µg kanamycin, 10 µg gentamicin, 10 µg streptomycin, 30 µg tetracycline, 30 µg chloramphenicol, 5 µg ciprofloxacin, 250 µg sulfisoxazole, 30 µg nalidixic acid, 50 µg fosfomycin, and 5 µg trimethoprim.
The presence of class I integrons and antibiotic resistance genes was investigated by PCR as described by Lai-King et al. (18), using primers for the 5' and 3' conserved segments of class I integrons and pse1, aadA2, sulI, pasppflo-like, and tet(G) genes. For each PCR, 2 ng DNA was amplified using FastStart Taq DNA polymerase (Roche) with a GeneAmp 9700 PCR system (Applied Biosystems). Each PCR involved one cycle at 94°C for 5 min; 30 cycles with the reaction conditions for each cycle being 94°C for 30 s, 58°C for 30 s, and 72°C for 2 min; and one cycle at 72°C for 7 min. Class I integron amplicon images were saved as tagged image format files and analyzed using BioNumerics version 4 software (Applied Maths, Sint-Martens-Latem, Belgium) to calculate amplicon size.
PFGE analysis.
PFGE analysis was performed as previously described (40) with some modifications. Briefly, the plug containing bacteria from an overnight culture was made with Seakem gold agarose (Cambrex, Rockland, ME) using a sample plug caster (Bio-Rad, Hercules, CA). The plug was treated for 18 h at 50°C with a solution of 1 mg proteinase K (Roche) per ml. After incubation, the plug was treated twice for 20 min each with Tris-EDTA (TE) buffer containing 4 mM Pefabloc (Roche) at 50°C and then washed twice on ice for 20 min each with TE buffer. The plug was digested for 18 h at 37°C with either XbaI (Roche) or BlnI (Roche). Salmonella enterica serovar Braenderup H9812 was used as the molecular weight standard (14). The standard was set at every fourth lane for the normalization of PFGE patterns (40).
The cluster analysis method used in these studies has been previously described (40). Briefly, PFGE patterns were imported into BioNumerics software (Applied Maths). Similarities among the PFGE patterns were calculated using the Pearson product-moment correlation coefficient (Pearson correlation) with a 0.5% optimization parameter. Dendrograms were made using the unweighted pair group method using arithmetic averages (UPGMA). Variations in grouping DT104 strains were investigated by changing the cutoff value from 95% to values at which all strains were grouped. A cluster was defined as a group of strains with higher similarity values than the cutoff value for cluster analysis of XbaI-digested PFGE patterns and BlnI-digested PFGE patterns.
VNTR typing.
DNA was isolated from DT104 strains using InstaGene Matrix (Bio-Rad). Two-nanogram DNA samples were amplified using FastStart Taq DNA polymerase (Roche). Primers for five VNTR loci, described by Lindstedt et al. (20), were used. A GeneAmp 9700 PCR system (Applied Biosystems) was used for PCR. Each PCR involved one cycle at 94°C for 5 min; 30 cycles with the reaction conditions for each cycle being 94°C for 30 s, 58°C for 30 s, and 72°C for 1 min; and one cycle at 72°C for 7 min. PCR products were partially purified using EXOSAP-it (USB, Cleveland, OH), and 10-fold dilution samples were made using 0.1x TE buffer. The samples were analyzed by capillary electrophoresis using an ABI-310 Genetic Analyzer (Applied Biosystems) with POP4 polymer (Applied Biosystems) and GeneMark 1000 Fluorescent DNA Ladder Rox (Northernbiotech, Weston, WI) as the internal control. Electrophoresis was carried out at 60°C for 40 min. The sizes of the DNA fragments in the samples were measured using GeneScan version 3.12 (Applied Biosystems), and the number of repeats in each VNTR locus was calculated as described by Lindstedt et al. (20). The number of repeats was rounded to two decimal places (41).
The number of repeats in each sample was imported into BioNumerics software (Applied Maths), and similarities among strains were calculated using the Pearson correlation. A dendrogram was made using the UPGMA method with the "use square root" option (41). The percent similarities based on the "use square root" option were calculated as follows:
![]() |
Variations in grouping DT104 strains were investigated following the same methods described above under "PFGE analysis." A cluster was defined as a group of strains with higher similarity values than the cutoff value for cluster analysis of VNTR types.
Correlation between the three clustering methods.
The congruence among the results of three cluster analyses of XbaI- and BlnI-digested PFGE patterns and VNTR data was measured using BioNumerics software (Applied Maths). The pairwise comparison of DT104 strains using one cluster analysis was compared to pairwise comparison of the isolates using another cluster analysis. The values of the similarity matrices of the two typing methods were plotted on an xy graph. Each point in the graph was the corresponding similarity value for the two typing methods. Such a plot shows the degree of correlation between the two methods. The congruence values among the three cluster analyses were calculated using BioNumerics software (Applied Maths) with Kendall's tau correlation under the option of 0% of "minimum similarity" and 100% of "maximum similarity" and without the option of "included self matches."
Of the three cluster analysis methods used in these studies (i.e., XbaI, BlnI, and VNTR cluster analysis), the method showing the most consistency with epidemiological information was used as the first step in a two-step cluster analysis. The clusters made by this first method were further analyzed using each of the two other cluster analysis methods to produce a two-step cluster analysis.
Statistical analysis.
Statistical analysis was performed using the Statcel2 (OMS Inc., Saitama, Japan) add-on package for Microsoft Excel. Chi-square analysis was used to compare numbers of rats carrying DT104 at each of the two study sites, as well as the numbers of R. rattus and R. norvegicus house rats carrying DT104 in each building. A P level of <0.05 was considered to be significant for all analyses.
|
|
|---|
Characteristics of DT104.
The 52 DT104 strains showed antimicrobial resistance against at least two drugs. Three antibiotic resistance patterns were observed: resistance to ampicillin, chloramphenicol, streptomycin, sulfisoxazole, and tetracycline (ACSSuT); resistance to ampicillin and sulfisoxazole (ASu); and resistance to streptomycin and sulfisoxazole (SSu). The ACSSuT resistance pattern was found only among DT104 strains isolated from house rats in building J. Two DT104 strains isolated from the rats in building J showed the SSu resistance pattern (Table 1). All DT104 strains isolated from house rats in building YS had only the ASu resistance pattern (Table 2).
PCR analysis using primers to the 5' and 3' conserved segments of class I integrons showed that all DT104 strains with the ACSSuT resistance pattern produced both the 1.0- and 1.2-kbp amplicons, and by the use of primers for specific antibiotic resistance genes, pse1, pasppflo-like, aadA2, sulI, and tet(G) were identified. All the strains with the SSu resistance pattern produced the 1.0-kbp amplicon of the class I integron with aadA2 and sulI, and all the strains with the ASu resistance pattern produced the 1.2-kbp amplicon of the class I integron with pse1 and sulI.
Molecular epidemiological analysis of DT104 strains.
Cluster analysis of XbaI-digested PFGE patterns (XbaI cluster analysis) showed that most or all DT104 strains from the two buildings formed one large cluster at all cutoff values (Table 3; see the XbaI dendrogram in the supplemental material). At the 95% cutoff value, the large cluster (cluster 95-X1) was made up of 41 of the 52 DT104 strains: 21 of the 26 building J strains (80.8%) and 20 of the 26 building YS strains (76.9%). At the 90% cutoff value, a large cluster (90-X1) was formed by 22 building J strains and 22 building YS strains. At the 85%, 80%, and 75% cutoff values, only a large cluster was formed, containing 49, 51, and 52 strains, respectively.
|
View this table: [in a new window] |
TABLE 3. Relationship between cutoff value and clustering of DT104 strains by use of three different cluster analysis methods
|
Cluster analysis of the VNTR data (VNTR cluster analysis) found one major and one minor cluster when the cutoff value was 95 or 90%. At the 95% cutoff value, one major cluster (95-V1) was formed by 24 building J strains and 15 building YS strains and one minor cluster (95-V2) was formed by eight building YS strains (Table 3; see the VNTR dendrogram in the supplemental material). At the 90% cutoff value, one major cluster (90-V1) was formed by all 26 building J strains and 15 building YS strains and one minor cluster (90-V2) was formed by eight building YS strains. At the 85%, 80%, and 75% cutoff values, 49 strains formed a large cluster at all values.
Correlation between the three clustering methods.
The congruence between the results of the XbaI and BlnI cluster analyses was low. The corresponding similarity values for these analyses are plotted in Fig. 1A, and the congruence value between the two analyses was 23.5 ± 1.83%. The congruence between the results of the BlnI and VNTR cluster analyses was also low. The corresponding similarity values for the two analyses are plotted in Fig. 1B, and the congruence value between the two analysis was 10.7 ± 1.83%.
![]() View larger version (12K): [in a new window] |
FIG. 1. Correlation of XbaI, BlnI, and VNTR cluster analyses. (A) Corresponding similarities of BlnI and XbaI cluster analysis. (B) Corresponding similarities of BlnI and VNTR cluster analysis.
|
|
View this table: [in a new window] |
TABLE 4. Relationship between cutoff value and clustering of DT104 strains in the clusters made by BlnI cluster analysis
|
|
|
|---|
The prevalence of DT104 strains in house rats, demonstrated in this study, supports the hypothesis of global expansion of DT104 by clonal dissemination (5, 29). Generally, the increase in antibiotic-resistant bacteria is due to selective pressure by antibiotics (7, 36). Since it is doubtful that such selective pressure could occur among house rats inhabiting buildings, dissemination by clonal expansion is a more likely conclusion. The hypothesis of clonal dissemination is supported by the molecular epidemiological results reported here. DT104 strains with similar PFGE patterns after XbaI digestion must have arisen from the same genetic clone (1, 4, 23). Most of the DT104 strains of this study had the same PFGE patterns following XbaI digestion (1, 16) and class I integrons with the same ACSSuT resistance pattern (18).
The global expansion of DT104 by clonal dissemination among house rats raises the difficult problem of tracing the local prevalence of the organism by comparison of PFGE patterns with XbaI digestion. Even though PFGE is the "gold standard" method for determining genetic relatedness of bacteria in molecular epidemiological studies (35), DT104 is highly clonal and XbaI digests of DT104 DNAs showed very similar PFGE patterns (1, 16). Therefore, other tools are necessary to investigate the molecular epidemiology of DT104 in house rats. Several molecular epidemiological methods have been investigated to differentiate DT104 strains (4, 19, 23). Plasmid profiling has been reported to be the most effective method for differentiating DT104 strains, compared to PFGE analysis, IS200 typing, and randomly amplified polymorphic DNA typing (1, 23). In contrast, plasmid profiling was not suitable for identifying DT104 clones in long-term studies because the plasmid studied was unstable (19). Since our study isolated DT104 strains over a 1-year period, we did not use plasmid profiling for molecular epidemiological analysis.
Another problem with the molecular epidemiological method is the question of the relationship between the similarity of PFGE patterns and the genetic relatedness of bacterial strains. Tenover et al. (35) recommended that the criteria for a genetic relationship could be based on differences in the numbers of restriction digest bands among PFGE patterns due to single-nucleotide polymorphisms at the restriction sites. However, the change of a PFGE pattern was due not primarily to the single-nucleotide polymorphism but to large-scale inversion (15, 17). Therefore, computer-based image analysis of PFGE patterns and cluster analysis of these data should be used to identify genetic relationships between bacterial strains. Cluster analysis can sort data into groups and reveal associations that are not otherwise evident. Although a number of studies have previously investigated cluster analysis of PFGE patterns, there were problems due to overlapping restriction bands (6) and the need to use a band matching coefficient, such as the Dice coefficient, for similarity calculations (40). In this study, the Pearson correlation, a densitometric curve-based coefficient, was used. Use of the Pearson correlation removes the problem of overlapping restriction bands, and the effectiveness of this approach has been confirmed by cluster analysis of PFGE patterns of other enteric bacteria (40).
In this study, XbaI cluster analysis could not differentiate DT104 strains. Most DT104 strains isolated from the two buildings studied were grouped in one large cluster. In contrast, analysis of the class I integron in these strains showed that the DT104 strains isolated from the two buildings contained different integrons, suggesting that different DT104 clones are prevalent in the house rats in the two buildings. The absence of DT104 strains with the ACSSuT or SSu resistance pattern in building YS suggests that DT104 strains with the ASu resistance pattern in building YS may have originated from ACSSuT strains by a single crossover at intI1 (3) before DT104 spread in this building. These results indicate that PFGE analyses with additional restriction enzymes will be required to trace the local prevalence of DT104 strains in house rats.
BlnI cluster analysis is more suitable than XbaI cluster analysis for tracing the local prevalence of DT104. PFGE analysis of BlnI digests has been reported to be more effective for differentiating DT104 strains than XbaI digestion (1, 23). In this study, BlnI cluster analysis showed that most DT104 strains from each building formed a single cluster per building at the 95% cutoff value. The BlnI cluster analysis results were in agreement with the integron analysis results for the DT104 strains in this study. All DT104 strains from building YS had the ASu resistance pattern, encoded by the same resistance genes in the same-size integron, while most strains from building J had the ACSSuT resistance pattern, encoded by the same resistance genes in the two same-size integrons. The BlnI cluster analysis and integron data and the circulation of salmonellae through habitats via carrier animals (13) suggest that different DT104 strains were prevalent in the house rat colonies of each building in this study. However, DT104 strains isolated from each of the two buildings showed similar BlnI-digested PFGE patterns throughout this 1-year study, indicating that it is impossible to follow the transmission of DT104 among house rats by this method.
Any single method of molecular analysis used in this study could not reveal the transmission of DT104 among house rats at a single location. DT104 strains derived from building YS were obviously divided into two distinct clusters using VNTR cluster analysis, but information about the appropriate dendrogram cutoff value is lacking. This result indicates that VNTR cluster analysis is useful for further differentiation of DT104 strains to follow the transmission of DT104 among house rats. However, strain differentiation among rats from the two buildings was not possible with only VNTR cluster analysis. Although the differentiating power of VNTR cluster analysis for DT104 strains will be improved by adding more VNTR loci (21), there is no information about the usefulness of VNTR loci other than those used in this study. One possible solution is to combine several analytical methods. The BioNumerics software can make composite data based on the results of several types of analytical methods. However, since similarities of PFGE patterns calculated by band matching coefficients are required to produce composite data (BioNumerics version 4 manual; Applied Maths), such a composite analysis could not be done in this study.
A two-step cluster analysis might partially trace the transmission of DT104 among house rats at a single location, with a VNTR cluster analysis of BlnI cluster analysis data. In this study, the congruence values of XbaI and VNTR cluster analysis results for BlnI cluster analysis results were low. These results indicated that clusters identified by BlnI cluster analysis could be further divided by XbaI or VNTR cluster analysis. Of the three clusters identified by BlnI cluster analysis, two were divided further into smaller clusters and unclustered strains by VNTR cluster analysis, while XbaI cluster analysis did not divide any of the BlnI clusters at any cutoff values. At the present time, it is unclear which cutoff values should be applied for VNTR cluster analysis. To determine the relatedness of bacterial strains, epidemiological information should be combined with genetic relatedness data (35). In this study, detailed epidemiological contact information among house rats was lacking. The efficiency of the two-step cluster analysis results for DT104 should be confirmed using epidemiological information.
It is unknown how many of the reported human cases of serovar Typhimurium in Japan are caused by DT104. Phage type is not always examined for human salmonellosis in Japan. However, DT104 has been expanding internationally (10). The DT104 strains isolated from human cases in Japan had the same genetic profile as the original DT104 clone (16), indicating that DT104 has already spread in Japan. Moreover, there have been reports that DT104 had also spread among livestock in Japan (8, 34), indicating a possible risk of carrier livestock contaminating food (13). Despite numerous reports of human DT104 infection in foreign countries (10), only one case was reported during a mass outbreak in Japan (33). During 1997-1998, one mass outbreak and three sporadic cases of serovar Typhimurium were reported in Yokohama City, but there were no phage type data (38, 39). Some cases of DT104 might be identified if serovar Typhimurium isolates from outbreaks were subjected to phage typing. PCR identification would be a useful tool to investigate whether serovar Typhimurium strains are DT104, although one strain showed contradictory results by PCR and phage typing. The positive reaction of the strain by Pritchett's methods suggested that the strain might be U302 (28). However, phage typing in this study was done to identify "DT104" or "not DT104," rather than to identify U302. In the study by Hermans et al. (11), no U302 strains were investigated, so it is still unclear whether PCR identification can identify DT104 only or a DT104 complex including U302.
In conclusion, the house rat plays an important role in the epidemiology of DT104, since DT104 can be maintained in the house rat colonies for a long time. BlnI cluster analysis is useful for investigating the molecular epidemiology of DT104, and two-step cluster analysis using VNTR typing may reveal further molecular epidemiology of DT104.
This work was supported by grants for Academic Frontier Project "Surveillance and Control for Zoonosis" from the Ministry of Education, Culture, Sports, Science and Technology of Japan.
Published ahead of print on 16 February 2007. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»