Previous Article | Next Article 
Applied and Environmental Microbiology, April 2008, p. 2118-2128, Vol. 74, No. 7
0099-2240/08/$08.00+0 doi:10.1128/AEM.02688-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
Molecular Analysis of Shiga Toxin-Producing Escherichia coli Strains Isolated from Hemolytic-Uremic Syndrome Patients and Dairy Samples in France
Nathalie Pradel,1
Yolande Bertin,2
Christine Martin,2 and
Valérie Livrelli3*
Laboratoire de Bactériologie, UFR Médecine, Université d'Auvergne, Clermont-Ferrand,1
Unité de Microbiologie, Institut National de la Recherche Agronomique, St.-Genès-Champanelle,2
Groupe de Recherche Pathogénie Bactérienne Intestinale, Faculté de Pharmacie, Université d'Auvergne Clermont-1, Unité Soutenue par l'INRA, Clermont-Ferrand, France3
Received 28 November 2007/
Accepted 22 January 2008

ABSTRACT
Shiga toxin-producing
Escherichia coli (STEC) has been associated
with food-borne diseases ranging from uncomplicated diarrhea
to hemolytic-uremic syndrome (HUS). While most outbreaks are
associated with
E. coli O157:H7, about half of the sporadic
cases may be due to non-O157:H7 serotypes. To assess the pathogenicity
of STEC isolated from dairy foods in France, 40 strains isolated
from 1,130 raw-milk and cheese samples were compared with 15
STEC strains isolated from patients suffering from severe disease.
The presence of genes encoding Shiga toxins (
stx1,
stx2, and
variants), intimin (
eae and variants), adhesins (
bfp,
efa1),
enterohemolysin (
ehxA), serine protease (
espP), and catalase-peroxidase
(
katP) was determined by PCR and/or hybridization. Plasmid profiling,
ribotyping, and pulsed-field gel electrophoresis (PFGE) were
used to further compare the strains at the molecular level.
A new
stx2 variant,
stx2-CH013, associated with an O91:H10 clinical
isolate was identified. The presence of the
stx2,
eae, and
katP genes, together with a combination of several
stx2 variants,
was clearly associated with human-pathogenic strains. In contrast,
dairy food STEC strains were characterized by a predominance
of
stx1, with a minority of isolates harboring
eae,
espP, and/or
katP. These associations may help to differentiate less virulent
STEC strains from those more likely to cause disease in humans.
Only one dairy O5 isolate had a virulence gene panel identical
to that of an HUS-associated strain. However, the ribotype and
PFGE profiles were not identical. In conclusion, most STEC strains
isolated from dairy products in France showed characteristics
different from those of strains isolated from patients.

INTRODUCTION
Shiga toxin-producing
Escherichia coli (STEC) causes human diseases,
ranging from uncomplicated diarrhea to hemorrhagic colitis (HC)
and life-threatening complications such as the hemolytic-uremic
syndrome (HUS). STEC food-borne infections appear worldwide,
with bovine feces being the main source of food contamination.
While most outbreaks are associated with
E. coli O157:H7, about
half of the sporadic cases may be due to non-O157:H7 serotypes.
Several studies have shown a high prevalence of STEC, belonging
to a wide range of serotypes, in animals and food products (
5,
34,
37). However, only a limited number of serotypes have been
associated with human disease. In France, the first identified
nationwide outbreak, consisting of 69 cases related to O157:H7
in ground meat, occurred in 2005 (A. Mailles, P. Mariani-Kurkdjian,
C. Vernozy-Rozand, F. Grimont, N. Pihier, E. Bingen, B. Horen,
V. Doireau, B. Llanas, E. Espié, and V. Vaillant, poster
presented at the 6th International Symposium on Shiga Toxin
(Verocytotoxin)-Producing
Escherichia coli Infections, Melbourne,
Australia, 2006,
http://www.invs.sante.fr/surveillance/shu/poster_melbourne_o157.pdf).
The major characteristic of STEC linked to virulence is the production of one or more Shiga toxins (Stx1, Stx2, and variants of Stx2) (31). Epidemiological studies, together with in vivo and in vitro experiments, have shown that Stx2 is the most common virulence factor associated with severe human disease. Three stx1 (stx1-933J, stx1c, stx1d) and at least 12 stx2 gene variants have been described, including stx2-EDL933, stx2c (stx2vh-a, stx2vh-b), stx2d (stx2d-Ount, stx2d-OX3a), stx2e, stx2f, and stx2-NV206 (3, 10, 24, 33, 44, 45). Some of the variants, which are associated with STEC strains isolated from specific hosts such as sheep (stx2d) and pig (stx2e), are probably less pathogenic for humans (18, 38). The type of variant could thus reflect both the origin and the clonality of STEC as well as its pathogenicity (18, 31).
Several other determinants have been involved in STEC virulence. The locus of enterocyte effacement (LEE) contains genes encoding proteins responsible for the "attaching and effacing" lesions in epithelial cells. One of these proteins, intimin, an outer membrane protein encoded by eae, is involved in the binding of bacteria to target cells. Several intimin types that may indicate the phylogenetic origin of STEC and determine the host tropism have been identified (1, 29). The LEE constitutes a pathogenicity island, a region of the bacterial chromosome transmitted by horizontal transfer, usually inserted adjacent to either the selC or the pheU tRNA locus (22, 27, 43). Large plasmids of STEC encode determinants that are thought to be additional virulence factors. These include enterohemolysin, which is encoded by the ehxA gene and acts as a pore-forming cytolysin on eukaryotic cells, an extracellular serine protease (EspP), and a catalase-peroxidase (KatP) (8, 9, 41). Phylogenetic studies have indicated that E. coli is composed of four main phylogenetic groups, A, B1, B2, and D, with STEC and enterohemorrhagic E. coli strains falling into phylogenetic groups A, B1, and D (14, 20). According to their incidence and association with HUS and outbreaks, STEC strains have also been classified into five seropathotypes (A to E) by Karmali et al. (25).
One of the more common sources of human infection is the consumption of undercooked ground beef, and many outbreaks described in the United States have been linked to this source. Dairy products have also been involved in outbreaks, although less frequently (39). In France, several clusters of HUS cases linked to dairy products have been identified and investigated (15), and there have been two recent reports of O157:H7- and O26-associated outbreaks (16; E. Espié, P. Mariani-Kurkdjian, F. Grimont, N. Pihier, V. Vaillant, S. Francart, H. de Valk, and C. Vernozy-Rozand, poster presented at the 6th International Symposium on Shiga Toxin (Verocytotoxin)-Producing Escherichia coli Infections, Melbourne, Australia, 2006, http://www.invs.sante.fr/surveillance/shu/poster_melbourne_o26.pdf). In the present study, we performed a detailed analysis of the molecular characteristics of STEC strains isolated from HUS patients and from dairy samples during prospective studies. Our aims were (i) to identify the characteristics of the strains isolated from HUS patients, (ii) to analyze the factors previously described as potentially associated with STEC virulence in a collection of dairy strains, (iii) to compare strains at the molecular level using powerful techniques such as ribotyping and pulsed-field gel electrophoresis (PFGE), and (iv) to assess the virulence of strains isolated from dairy products.

MATERIALS AND METHODS
Bacterial strains and culture conditions.
A total of 55 STEC isolates collected in France were examined
in the study, of which 40 were isolated from 1,130 dairy samples,
which included 205 raw-milk (166 from cows and 39 from ewes),
828 cheese (781 from cows, 32 from ewes, and 15 from goats),
and 97 "dairy environment" (87 from cows and 10 from ewes) samples
taken during previous studies (
17,
37). The samples were randomly
collected over a 2-year period from several regions in France.
They were purchased from local retail stores, provided by the
producers, or obtained from the Laboratoires Départementaux
d'Analyses Vétérinaires et Biologiques. The different
types of cheeses and milk were not epidemiologically related
and were representative of French production as a whole. A total
of 15 STEC samples were isolated from patients in central France
admitted to hospitals for HUS or HC during the same period.
The strains are representative of STEC isolated from patients
in a French hospital and include strains of serogroups O157,
O103, and O91. Serotype determination of HUS/HC-associated strains
was performed by Flemming Scheutz at the International
Escherichia and
Klebsiella Reference Centre (WHO) in Copenhagen, Denmark.
For the dairy STEC sample provided by Fach et al. (
17), only
the O group had been determined.
E. coli EDL933 (ATCC 43895)
serotype O157:H7 was used as an enterohemorrhagic
E. coli reference
strain, and
E. coli K-12 was used as a negative control for
PCR experiments.
Detection of genes associated with virulence.
This detection was done by PCR using the primers listed in Table 1. The DNA to be amplified was released from whole organisms by boiling the samples. The PCR cycles included denaturation for 60 s at 94°C, primer annealing for 60 s at 45 to 68°C (primer-specific annealing temperatures are shown in Table 1), and an extension for 60 or 120 s (according to the length of the fragment) at 72°C (30 cycles). The amplification reactions were performed on a Perkin-Elmer GeneAmp 2400 PCR system. The presence of three genes (espP, ehxA, katP) was also determined by colony blot hybridization, using probes labeled with [
-32P]dCTP, as described previously (36). The probes used for the detection of target genes were obtained by PCR from strain EDL933, using the primer pairs espA-espB for espP, RH35-RH37 for ehxA, and wkatb-wkatf for katP (Table 1). The probe sizes were 1,830 bp, 321 bp, and 2,125 bp, respectively.
PCRs for LEE insertions.
The insertion of the LEE adjacent to
selC was tested using the
four primer pairs described by McDaniel et al. (
27): K260-K261
(intact
selC, 527 bp), K255-K260 (right junction of the LEE
inserted into
selC, 418 bp), K295-K296 (left junction of the
LEE inserted into
selC, 405 bp), and K260-K295 (intact
selC,
2,173 bp). The primer pair K913-K914 was used to amplify the
pheU gene (intact
pheU, 303 bp) (
43).
stx1 subtyping.
The three stx1 gene variants (stx1-933J, stx1c, and stx1d; GenBank accession numbers M19473, Z36901, and AY170851, respectively) were amplified using the GannonF-GannonR primer pair. PCR products were then subjected to restriction endonuclease digestion with BglI, HaeI, and RsaI (Roche Applied Science). In silico analysis revealed that the 603-bp amplified product from STEC strains possessing stx1-933J (strain EDL933), stx1c (strain DG131/3), and stx1d (strain MHI813) generated fragments of 215 and 387 bp with BglI, 220 and 382 bp with RsaI, and 415 and 187 bp with HaeI, respectively. The presence of the stx1c variant was confirmed by PCR amplification using the VT1AvarF-VT1AvarR and lin-up-AOX3 primer pairs (26).
Detection of stx2 variants by PCR.
Restriction fragment length polymorphism (RFLP)-PCR systems were used to distinguish stx2 subtypes, as previously described (3, 33, 38, 44). The VT2c-VT2d primer pair was used to detect the stx2-EDL933, stx2vh-a, stx2vh-b, and stx2-NV206 genes. The 285-bp PCR product obtained was then subjected to restriction endonuclease digestion with HaeIII, RsaI, and NciI (Boehringer Mannheim, Meylan, France), as described by Piérard et al. (33). To identify STEC strains harboring the stx2-NV206 gene, the HpaI enzyme was also used (Fig. 1) (3). The VT2cm-VTf primer pair was used to specifically detect the stx2d genes; discrimination between stx2d-Ount and stx2d-OX3a was performed with the VT2e and VT2f primers, the 348-bp amplicon being digested with HaeIII and PvuII. The stx2e gene variant was detected using the VT2ea-VT2eb primer pair. DNA from the E. coli strains EDL933 (O157:H7 stx2-EDL933), B2F1 (O91:H21 stx2vh-a, stx2vh-b), NV206 (stx2-NV206), Fac9 (stx2e), an OX3:H11 (stx2d) strain, and K12C600 (negative control) was included in each PCR run. Digestion products were separated on 2% agarose gels, and the standard-size 100-bp DNA ladder marker (Promega, France) was used as the referee.
Phylogenetic-group determination.
The main phylogenetic groups (A, B1, B2, and D) of the
E. coli strains were determined by triplex PCR amplification as described
by Clermont et al. (
14).
Preparation of plasmid DNA.
Plasmid DNA was prepared from 1.5-ml overnight cultures of bacteria grown in Luria broth (Difco Laboratories, Detroit, MI) by the alkaline lysis method. Plasmid DNA patterns were obtained by electrophoresis on 0.7% agarose gels in Tris-acetate-EDTA buffer, stained with ethidium bromide, and observed under UV light.
Ribotyping.
Total genomic DNA of STEC isolates was prepared from 10-ml overnight cultures in Mueller-Hinton broth (Biokar Diagnostics, Beauvais, France) by the method described by Picard-Pasquier et al. (32). About 3 µg of genomic DNA was digested independently with two restriction endonucleases, EcoRI and HindIII (Boehringer Mannheim), according to the manufacturer's instructions. Digested DNA was separated, transferred to Hybond N+ nylon membranes (Amersham Pharmacia Biotech, Orsay, France), and probed with an rrnB7 probe as previously described (36). The RFLP Extension Ladder System (Life Technologies, Cergy Pontoise, France) was used as the molecular weight marker. The Dice similarity coefficient and the Diversity Database software (Bio-Rad) were used to compare profiles. The hierarchic unweighted pair group arithmetic average algorithm was used for cluster analysis, and dendrographic trees were constructed.
PFGE analysis.
Genomic DNA was prepared according to the protocol described by Böhm and Karch (7) using low-melting-point agarose (Gibco BRL). For restriction endonuclease digestion, 50 U of XbaI (Life Technologies, Cergy Pontoise, France) was used at 37°C for 18 h. PFGE was performed in 1% agarose (Invitrogen) gels using a contour-clamped homogeneous electric field DRIII apparatus (Bio-Rad, Ivry, France) in 0.5x Tris-borate-EDTA buffer at 14°C at 6 V/cm. The pulse time was increased from 5 to 50 s over 22 h. A lambda ladder (Bio-Rad) was used as the size marker.
Sequencing of the new stx2 variant.
The PCR product for sequencing was prepared by using the QIAquick PCR purification kit (Qiagen S.A.). DNA sequencing was carried out on an automated sequencer by the fluorescent-dye termination method, and the sequence was edited using the manufacturer's software (GENOME Express, Grenoble, France). BLASTN and BLASTX sequence homology analyses were performed using the National Center for Biotechnology Information BLAST network service (2).
Statistical analysis.
The data were analyzed with Epi-Info version 6.02 software by the
2 test, except for the variable needing the two-tailed Fisher exact test. A P value of <0.05 was considered statistically significant.
Nucleotide sequence accession number.
The sequence of the new stx2 variant has been deposited in the GenBank database under accession no. EU184879.

RESULTS
Detection of genes associated with virulence and significant associations.
Forty STEC strains from 1,130 dairy products or dairy environments
and 15 pathogenic STEC strains were analyzed. Data regarding
the virulence factor-encoding genes
stx1,
stx2,
eae,
efa1,
ehxA,
espP,
katP, and
astA and their distribution according to the
origin of the STEC strains (dairy products versus HUS patients)
are given in Table
2. The characteristics of each of the 55
STEC isolates tested are shown in Table
3. The
stx1 gene was
found, alone or in association with
stx2, in 25 out of 40 dairy
STEC strains (62%) and in 4 out of 15 HUS-associated strains
(27%) (
P = 0.018). In contrast, the
stx2 gene was present in
20 (50%) dairy isolates and in 12 (80%) HUS-associated strains
and was therefore significantly associated with pathogenic strains
(
P = 0.045). The
stx1-
stx2 association was not frequent and
was almost equally distributed in both categories (six dairy
strains [15%] and one HUS patient isolate [7%]). Only 5 dairy
strains out of 40 (12%) were positive for
eae; they belonged
to the same serogroup (O5) and shared the same characteristics
(Table
3). Interestingly, 8 of the 15 HUS strains (53%) harbored
the
eae gene, indicating a strong association with pathogenic
isolates (
P = 0.003). As expected, similar results were obtained
with the
efa1 gene (
28), which was always associated with
eae,
except in the CH071 strain (Table
3). The enterohemolysin-encoding
gene
ehxA, found in 15 dairy STEC isolates (37%) and in 9 HUS
patient isolates (60%), was not clearly associated with any
category. The
espP gene was found more frequently in strains
from patients (53%) than in dairy isolates (22%) (
P = 0.032).
PCR and Southern blot hybridization demonstrated that the
katP gene was rarely present in dairy STEC strains (only one strain,
PL475, harbored this gene) but was found in 40% of the HUS isolates
(
P = 0.001). The enteroaggregative heat-stable enterotoxin-encoding
gene
astA was present in nine dairy STEC strains but was not
found in any patient isolates. Finally, the
bfp gene, encoding
the pili involved in the adhesion of enteropathogenic
E. coli strains to intestinal cells, was absent from all the strains
tested, which indicates that all tested isolates belong to the
STEC group of diarrheogenic
E. coli organisms.
Analysis of the stx1 gene variants.
Three main
stx1 variants,
stx1-933J (identified in the EDL933
reference strain),
stx1-OX3 (also called
stx1c), and
stx1d,
have been described in the literature (
10,
26). Using the lin-up-lin-OX3
primers, we showed that 9 strains (all isolated from dairy products)
possessed
stx1c and 20 strains possessed
stx1-933J (the 4
stx1-positive
HUS isolates and 16 out of 25
stx1-positive dairy strains).
No
stx1d variant was detected. No more than one variant was
observed.
Analysis of the stx2 gene variants.
To further assess the pathogenicity of the STEC isolates, the stx2 gene subtypes were identified by PCR and RFLP-PCR. Four restriction enzymes were used in this study: HaeIII, RsaI, NciI, and HpaI. Representative patterns are shown in Fig. 1, and the results of this analysis are summarized in Table 3. The 32 stx2-positive isolates were found to contain stx2-EDL933, stx2vh-a, stx2vh-b, stx2-NV206, or stx2d, but none was found to possess the stx2e gene. The strains contained either one stx2 gene type or a combination of two or three. The CH013 strain showed an atypical RsaI profile (Fig. 1). Cloning and sequencing experiments revealed that it contained two variants: stx2-EDL933 and a new stx2 variant, stx2-CH013 (GenBank accession number EU184879), whose sequence showed only 92% identity with the sequence of stx2-NV206. The distribution of the stx2 variants according to the origin of STEC (dairy products versus HUS patients) is given in Table 4. Of the 12 stx2-positive HUS strains, 9 (75%) harbored the stx2-EDL933 variant, indicating a strong association with pathogenic isolates (P = 0.0003). The stx2vh-a variant was more frequent in the HUS isolates (P = 0.018), and the stx2vh-b and stx2-NV206 genes were more frequent among dairy strains but not to a significant degree. The stx2d variant, which was recently associated with ovine isolates (38), was absent from the HUS strains and was found in only two dairy isolates. Interestingly, both of the cheese samples from which the last-named strains were isolated were made from ewe's milk. Finally, the stx2e variant was absent from all the strains tested, which was not surprising, since the variant is clearly associated with STEC strains isolated from pigs. Six out of 12 (50%) HUS patient strains harbored several stx2 variants (Table 3), while only 4 dairy isolates harbored a combination of two genes (stx2vh-b and stx2-NV206).
Analysis of the LEE.
The
eae gene variants were identified from the 13
eae-positive
isolates (Table
3). β Intimin was present in the five dairy
strains of serogroup O5 (ECA34 to ECA43). The strains isolated
from patients were more heterogeneous: two harbored the β
variant (the O157:H26 isolate CH071 and the O5:H
– isolate
CH123), three harbored the

variant (including the two strains
of serotype O103:H2, CH087 and CH089), and three possessed the

variant (all of serotype O157:H7). The sites of insertion of
the LEE were investigated using PCR systems previously described
(
27,
43), and the results of this analysis are shown in Table
5. In the five dairy STEC isolates belonging to the O5 serogroup
(ECA34 to ECA43), both the
selC and
pheU sites were not intact.
However, the signals indicating the insertion of the LEE into
selC (i.e., a 418-bp fragment with K255 and K260 and a 405-bp
fragment with K295 and K296) were missing. Therefore, if the
LEE is inserted into
selC or
pheU, it is "modified" in the region
where the primers were chosen. Identical results were obtained
for two isolates from HUS patients, CH123 (which belongs to
serotype O5:H
– and carries a β intimin) and CH023
(O
+:H
–,

intimin). The LEE was inserted into
selC in the
four pathogenic isolates belonging to the O157 serogroup (CH071,
CH075, CHVi-1, CH1898). For three of them, the LEE was modified
either in the left junction (CH071) or in the right junction
(CHVi-1, CH1898). An intact
selC region was observed in three
isolates. For two of them, strains CH087 and CH089, which belong
to the O103:H2 serotype and carry the

intimin, the LEE might
be inserted into
pheU.
Plasmid content analysis.
The presence or absence of plasmids and their number and size
were investigated in the 55 STEC strains. The plasmid patterns
were heterogeneous: 11 isolates did not carry any visible plasmids,
and 44 strains carried one to three plasmids that ranged from
approximately 3 to more than 90 kb in size (data not shown).
Some isolates belonging to the same serotypes presented distinct
profiles, further evidence of the high heterogeneity of STEC
isolates with regard to mobile elements such as plasmids. The
numbers of plasmids according to the origins of the STEC isolates
(dairy product versus HUS patients) are shown in Table
6. No
significant association was observed.
View this table:
[in this window]
[in a new window]
|
TABLE 6. Distribution of the plasmids observed after extraction and migration in a 0.8% agarose gel according to the origins of the STEC isolates
|
Analysis of the phylogenetic group.
Of the 40 STEC dairy isolates, 16 belonged to the A group, 20
to the B1 group, and 4 to the D group (Table
3). Most of the
strains isolated from patients belonged to the B1 group (12
strains), and the three strains of serotype O157:H7 belonged
to the D group, but none fell into the A group, confirming that
the most virulent strains belong either to the D or to the B1
groups (
20). Interestingly, a large number of isolates from
dairy products belong to the A group (Fisher exact test,
P =
0.002), which is associated with less virulent STEC strains.
Ribotyping profiles of the STEC isolates.
The 55 STEC isolates were analyzed by ribotyping using the rrnB7 probe. A total of 23 different patterns were obtained with HindIII-digested genomic DNA, but two isolates (NV63 and PL915) could not be typed. When the EcoRI enzyme was used, 28 different patterns were observed, but six isolates (CH016, PL27, 1429, PL325, PL915, ECA64) could not be typed. The HindIII ribotyping patterns showed between 7 and 10 fragments ranging from 1 to 26 kb in length. Figure 2 shows a dendrographic tree constructed from these HindIII profiles. The six O6 dairy isolates exhibited the same pattern, which was identical to those of O110 strains PL475 and PL650, O113 strain NV81, and O136 strain PL652. The four O157:H7 strains (including the reference strain EDL933) presented identical patterns, which were the same as those of the O3 and O15 isolates. Five O5 isolates (ECA34, ECA36, ECA37, ECA38, ECA43) out of seven shared the same profile. The O5 strain isolated from a patient (CH123) presented a slightly different pattern, and the last O5 dairy isolate, 2206, had a distinct ribotype profile, which indicates a more distant relationship. Four nontypeable isolates (14, 15, 40, PL324) out of five had identical ribotype profiles, as did two O77 isolates, three O91 isolates, three O76 isolates, two O103 isolates, a group of two isolates (one O117 and one O+ isolate), and a group of three isolates including one O174 strain and two O8 strains. The STEC strains isolated from patients were scattered along the tree and showed patterns similar to those of dairy STEC strains (Fig. 2). However, only six of the latter showed profiles identical to those of clinical isolates, and only one of them (ECA95, an O91 strain) also exhibited an EcoRI profile identical to that of clinical isolates (CH013 and CH085, both of which are O91:H10 strains) (data not shown). This may indicate that the great majority of dairy STEC strains are distantly related to pathogenic isolates.
PFGE analysis.
To further determine relatedness among human and dairy isolates
showing similar ribotypes, we analyzed their XbaI patterns by
PFGE. Ten to 19 fragments, ranging in size from about 50 kb
to 580 kb, were observed. Strains showing identical ribotypes
but belonging to different serotypes had only distantly related
PFGE patterns, confirming the discriminatory power of PFGE.
Most of the isolates were not genetically related, since their
PFGE patterns showed more than six band differences. Of the
strains with identical ribotypes, four O6 strains (PL599, -600,
-602, -603) had identical PFGE patterns, suggesting that they
belong to the same clone (named hereinafter the PL599 clone).
The PFGE pattern of the O6 strain PL479 showed more than six
band differences from the latter. Interestingly, this strain
harbored only the
stx2-NV206 variant, whereas the PL599 clone
harbored a combination of two genes (
stx2vh-b and
stx2-NV206).
The PFGE pattern of PL230, the last O6 isolate, was completely
different from the others, indicating a more distant relatedness.
The five O5 isolates ECA34, ECA36, ECA37, ECA38, and ECA43 showed
highly related restriction patterns, with just one additional
band for ECA34 and two additional bands for ECA43. The profiles
of the two other O5 isolates (including the HUS isolate CH123)
were different and more distantly related. The profiles of the
four O157:H7 isolates were related (11 bands in common) but
differed by 3 to 5 bands. Among the O91 isolates, the dairy
strain ECA95 had a profile close to that of the clinical O91:H10
strain CH013, but these two profiles differed by more than six
bands.

DISCUSSION
In this study, we performed a detailed analysis of the molecular
epidemiology of STEC, using a collection of 55 STEC isolates
from sporadic cases of HUS or HC and from milk or other dairy
products from France. The two groups of strains had clearly
different characteristics. Regarding virulence genes, STEC strains
associated with severe human disease were statistically associated
with the
stx2 gene (particularly with the
stx2-EDL933 and
stx2vh-a variants and frequently with several
stx2 gene variants), but
they were also associated with
eae,
espP, and
katP. In contrast,
dairy food STEC strains were characterized by a predominance
of
stx1, with a minority of isolates harboring
eae,
espP, and/or
katP. Recent studies of strains isolated from cheese and other
dairy products evidenced the difficulties encountered in defining
virulence markers of STEC (
4,
11,
40,
46). Our findings supply
further information. Although the
ehxA gene had previously been
suspected to be a marker for virulence (
41), possibly because
of a preferential association of
ehxA with
eae-positive isolates
(
6,
37), we did not find any association of
ehxA with pathogenic
isolates. In contrast, our data revealed a very strong association
of
katP (and to a lesser extent
espP) with pathogenic strains.
KatP is a catalase-peroxidase which is thought to enhance the
resistance of strains to oxidative stress, but its role in virulence
remains unclear (
8). EspP is able to degrade coagulation factor
V, which could increase hemorrhage into the gastrointestinal
tract (
9). In our study, the
katP gene was found only in
eae-positive
strains, with one exception (the dairy O110 isolate PL475),
and similar results were found when analyzing bovine STEC strains
(
37; unpublished data). In STEC strains of serogroup O103, O26,
O111, and O157, the
katP gene was also always present in
eae-positive
isolates (
42). This raises the question of the preferential
association of plasmid-borne genes with chromosomal genes such
as
eae.
Among bovine STEC strains isolated in central France, the most frequent stx2 subtypes were stx2vh-b (39.5% of the isolates), stx2-EDL933 (39%), and stx2vh-a (25.5%), followed by stx2-NV206 (14.5%) and stx2d (8.5%) (3). In dairy STEC strains, the stx2vh-b variant predominated (60%) over stx2-NV206 (30%), stx2-EDL933 (10%), stx2vh-a (10%), and stx2d (10%). In contrast, the stx2-EDL933 and stx2vh-a variants were overrepresented in strains isolated from patients (75% and 50%, respectively), suggesting a major role in the genesis of severe disease in humans. The percentages of stx2d-positive strains were similar in dairy STEC strains and in bovine strains. The stx2d variants, which were recently associated with STEC strains isolated from sheep (38), are rarely found in STEC strains involved in severe disease, suggesting that stx2d production is not a characteristic of human-pathogenic strains. This observation is in correlation with the results obtained by Vernozy-Rozand et al. (46). Finally, we identified a new stx2 variant, stx2-CH013 (GenBank accession number EU184879), in an O91:H10 HUS isolate and observed a combination of several stx2 variants in 27% of the clinical isolates. We previously showed that HUS-associated strains were more cytotoxic than their bovine counterparts belonging to the same serotype (37). The number of stx2 gene variants could affect the level of cytotoxicity and thus the level of pathogenicity of clinical isolates.
The most frequent serogroups isolated from cattle during a previous study were OX3/O174, O113, O22, O91, O172, O6, OX178, O171, O46, and O74. Strains belonging to these serogroups accounted for 70% of the total number of bovine isolates but for only 32% of the dairy STEC isolates (37). This, together with data regarding virulence factors and stx2 variants, confirms that the STEC strains isolated from cheese samples and from infected patients represent only particular subsets of strains isolated from cattle. It is conceivable that STEC strains isolated from cheese possess special properties enabling them to survive during the cheese-making process and that the properties needed for full virulence in humans are different.
One of the molecular techniques used to assess relationships between the STEC strains was the determination of the LEE insertion site. Because of numerous modifications in the regions surrounding the tRNA genes, which can be considered "hot spots" of insertion, the PCR systems described to date are not efficient enough to determine into which site the LEE is inserted. Although these PCR systems provide important information on clonal relationships between isolates, new tools are needed to better define the locus into which the LEE is inserted. Ribotyping was useful for establishing relationships among STEC strains belonging to a wide range of serogroups, but it was not able to discriminate between STEC isolates of identical serotypes. Plasmid profile analysis has been a useful tool in molecular epidemiology for the interstrain differentiation of STEC isolates belonging to the same serotype (36). In the present study, this method provided confirmation of a common clonal origin of STEC strains of identical serotypes having identical virulence patterns, but no significant association was found between the presence and/or number of the plasmids and the pathogenicity of the strains.
The molecular techniques allowing analysis of clonal relationships (i.e., subtyping of intimin, determination of the LEE insertion site, plasmid profiling, ribotyping, and PFGE) revealed that two groups of four O6 and five O5 isolates corresponded to two clones. The O5 "ECA34 clone" exhibited characteristics very close to those of a strain isolated from a child suffering from HC in 2000 (CH123); both belong to the O5 serogroup, harbor the β intimin variant, and have identical profiles for LEE insertion (Table 5). However, the HindIII ribotypes and PFGE patterns were different, indicating that the O5 dairy clone was probably not the origin of this clinical strain. A dairy strain belonging to the O91 serogroup (ECA95) showed HindIII and EcoRI ribotypes identical to those of two STEC strains belonging to the O91:H10 serotype, isolated from HUS patients in March 1997 (CH013) and in December 1998 (CH085). The pulsotypes of the clinical strains CH013 and CH085 differed by more than six bands from that of ECA95. These strains harbored none of the virulence factors described previously except the stx2 gene. Interestingly, the dairy isolate harbored only the stx2vh-a variant, whereas the clinical isolates harbored two stx2 variants, which may account for the differences in the pulsotypes and possibly in virulence power.
In conclusion, our results indicate that STEC strains isolated from dairy food samples in France possessed major characteristics different from those of strains isolated from patients. STEC strains involved in severe human disease were statistically associated with the presence of several stx2 gene variants, mainly stx2-EDL933 and stx2vh-a, but also with eae, espP, and katP. In contrast, dairy food STEC strains were characterized by a predominance of stx1 and belonged to phylogenetic A group. These findings may be helpful for the selective detection of strains that present a risk to public health.

ACKNOWLEDGMENTS
We thank Patrick Fach for providing the dairy STEC strains (isolates
from a previous program supported in part by ARILAIT Recherches)
and Chantal Rich and Christiane Forestier for the isolation
of the STEC strains from HUS/HC patients (supported by a Programme
Hospitalier de Recherche Clinique National). We thank Christophe
De Champs for invaluable help in statistical analysis and Muriel
Millet and Nathalie Banchet for technical assistance.
This study was supported in part by ARILAIT Recherches and by the Ministère de l'Enseignement Supérieur et de la Recherche (EA2348).

FOOTNOTES
* Corresponding author. Mailing address: Groupe de Recherche Pathogénie Bactérienne Intestinale, Faculté de Pharmacie, Université d'Auvergne, 28 Place Henri-Dunant, 63 001 Clermont-Ferrand, France. Phone: (33) 473 178 375. Fax: (33) 473 178 371. E-mail:
Valerie.Livrelli{at}u-clermont1.fr 
Published ahead of print on 1 February 2008. 

REFERENCES
1 - Adu-Bobie, J., G. Frankel, C. Bain, A. G. Goncalves, L. R. Trabulsi, G. Douce, S. Knutton, and G. Dougan. 1998. Detection of intimins
, β,
, and
, four intimin derivatives expressed by attaching and effacing microbial pathogens. J. Clin. Microbiol. 36:662-668.[Abstract/Free Full Text] 2 - Altschul, S. F., T. L. Madden, A. A. Schäffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25:3389-3402.[Abstract/Free Full Text]
3 - Bertin, Y., K. Boukhors, N. Pradel, V. Livrelli, and C. Martin. 2001. Stx2 subtyping of Shiga toxin-producing Escherichia coli isolated from cattle in France: detection of a new Stx2 subtype and correlation with additional virulence factors. J. Clin. Microbiol. 39:3060-3065.[Abstract/Free Full Text]
4 - Beutin, L., A. Miko, G. Krause, K. Pries, S. Haby, K. Steege, and N. Albrecht. 2007. Identification of human-pathogenic strains of Shiga toxin-producing Escherichia coli from food by a combination of serotyping and molecular typing of Shiga toxin genes. Appl. Environ. Microbiol. 73:4769-4775.[Abstract/Free Full Text]
5 - Blanco, M., J. E. Blanco, J. Blanco, A. Mora, C. Prado, M. P. Alonso, M. Mourino, C. Madrid, C. Balsalobre, and A. Juarez. 1997. Distribution and characterization of faecal verotoxin-producing Escherichia coli (VTEC) isolated from healthy cattle. Vet. Microbiol. 54:309-319.[CrossRef][Medline]
6 - Boerlin, P., S. A. McEwen, F. Boerlin-Petzold, J. B. Wilson, R. P. Johnson, and C. L. Gyles. 1999. Associations between virulence factors of Shiga toxin-producing Escherichia coli and disease in humans. J. Clin. Microbiol. 37:497-503.[Abstract/Free Full Text]
7 - Böhm, H., and H. Karch. 1992. DNA fingerprinting of Escherichia coli O157:H7 strains by pulsed-field gel electrophoresis. J. Clin. Microbiol. 30:2169-2172.[Abstract/Free Full Text]
8 - Brunder, W., H. Schmidt, and H. Karch. 1996. KatP, a novel catalase-peroxidase encoded by the large plasmid of enterohaemorrhagic Escherichia coli O157:H7. Microbiology 142:3305-3315.[Abstract/Free Full Text]
9 - Brunder, W., H. Schmidt, and H. Karch. 1997. EspP, a novel extracellular serine protease of enterohaemorrhagic Escherichia coli O157:H7 cleaves human coagulation factor V. Mol. Microbiol. 24:767-778.[CrossRef][Medline]
10 - Bürk, C., R. Dietrich, G. G. Açar, M. Moravek, M. Bülte, and E. Märtlbauer. 2003. Identification and characterization of a new variant of Shiga toxin 1 in Escherichia coli ONT:H19 of bovine origin. J. Clin. Microbiol. 41:2106-2112.[Abstract/Free Full Text]
11 - Caro, I., and M. R. Garcia-Armesto. 2007. Occurrence of Shiga toxin-producing Escherichia coli in a Spanish raw ewe's milk cheese. Int. J. Food Microbiol. 116:410-413.[CrossRef][Medline]
12 - Cebula, T. A., W. L. Payne, and P. Feng. 1995. Simultaneous identification of strains of Escherichia coli serotype O157:H7 and their Shiga-like toxin type by mismatch amplification mutation assay-multiplex PCR. J. Clin. Microbiol. 33:248-250.[Abstract]
13 - China, B., F. Goffaux, V. Pirson, and J. Mainil. 1999. Comparison of eae, tir, espA and espB genes of bovine and human attaching and effacing Escherichia coli by multiplex polymerase chain reaction. FEMS Microbiol. Lett. 178:177-182.[CrossRef][Medline]
14 - Clermont, O., S. Bonacorsi, and E. Bingen. 2000. Rapid and simple determination of the Escherichia coli phylogenetic group. Appl. Environ. Microbiol. 66:4555-4558.[Abstract/Free Full Text]
15 - Decludt, B., P. Bouvet, P. Mariani-Kurkdjian, F. Grimont, P. A. Grimont, B. Hubert, and C. Loirat. 2000. Haemolytic uraemic syndrome and Shiga toxin-producing Escherichia coli infection in children in France. Epidemiol. Infect. 124:215-220.[CrossRef][Medline]
16 - Espié, E., V. Vaillant, P. Mariani-Kurkdjian, F. Grimont, R. Martin-Schaller, H. de Valk, and C. Vernozy-Rozand. 2006. Escherichia coli O157 outbreak associated with fresh unpasteurized goats' cheese. Epidemiol. Infect. 134:143-146.[Medline]
17 - Fach, P., S. Perelle, F. Dilasser, and J. Grout. 2001. Comparison between a PCR-ELISA test and the vero cell assay for detecting Shiga toxin-producing Escherichia coli in dairy products and characterization of virulence traits of the isolated strains. J. Appl. Microbiol. 90:809-818.[CrossRef][Medline]
18 - Friedrich, A. W., M. Bielaszewska, W. L. Zhang, M. Pulz, T. Kuczius, A. Ammon, and H. Karch. 2002. Escherichia coli harboring Shiga toxin 2 gene variants: frequency and association with clinical symptoms. J. Infect. Dis. 185:74-84.[CrossRef][Medline]
19 - Gannon, V. P., R. K. King, J. Y. Kim, and E. J. Golsteyn Thomas. 1992. Rapid and sensitive method for detection of Shiga-like toxin-producing Escherichia coli in ground beef using the polymerase chain reaction. Appl. Environ. Microbiol. 58:3809-3815.[Abstract/Free Full Text]
20 - Girardeau, J. P., A. Dalmasso, Y. Bertin, C. Ducrot, S. Bord, V. Livrelli, C. Vernozy-Rozand, and C. Martin. 2005. Association of virulence genotype with phylogenetic background in comparison to different seropathotypes of Shiga toxin-producing Escherichia coli isolates. J. Clin. Microbiol. 43:6098-6107.[Abstract/Free Full Text]
21 - Gunzburg, S. T., N. G. Tornieporth, and L. W. Riley. 1995. Identification of enteropathogenic Escherichia coli by PCR-based detection of the bundle-forming pilus gene. J. Clin. Microbiol. 33:1375-1377.[Abstract]
22 - Hacker, J., and J. B. Kaper. 2000. Pathogenicity islands and the evolution of microbes. Annu. Rev. Microbiol. 54:641-679.[CrossRef][Medline]
23 - Hall, R. H., and J. G. Xu. 1995. December A new and distinctive DNA sequence of E. coli O157:H7 and its use for the rapid, sensitive and specific detection of O157:H7 and other enterohemorrhagic E. coli. International patent WO 95/34682.
24 - Johnson, W. M., D. R. Pollard, H. Lior, S. D. Tyler, and K. R. Rozee. 1990. Differentiation of genes coding for Escherichia coli verotoxin 2 and the verotoxin associated with porcine edema disease (VTe) by the polymerase chain reaction. J. Clin. Microbiol. 28:2351-2353.[Abstract/Free Full Text]
25 - Karmali, M. A., M. Mascarenhas, S. Shen, K. Ziebell, S. Johnson, R. Reid-Smith, J. Isaac-Renton, C. Clark, K. Rahn, and J. B. Kaper. 2003. Association of genomic O island 122 of Escherichia coli EDL 933 with verocytotoxin-producing Escherichia coli seropathotypes that are linked to epidemic and/or serious disease. J. Clin. Microbiol. 41:4930-4940.[Abstract/Free Full Text]
26 - Koch, C., S. Hertwig, R. Lurz, B. Appel, and L. Beutin. 2001. Isolation of a lysogenic bacteriophage carrying the stx1(OX3) gene, which is closely associated with Shiga toxin-producing Escherichia coli strains from sheep and humans. J. Clin. Microbiol. 39:3992-3998.[Abstract/Free Full Text]
27 - McDaniel, T. K., K. G. Jarvis, M. S. Donnenberg, and J. B. Kaper. 1995. A genetic locus of enterocyte effacement conserved among diverse enterobacterial pathogens. Proc. Natl. Acad. Sci. USA 92:1664-1668.[Abstract/Free Full Text]
28 - Nicholls, L., T. H. Grant, and R. M. Robins-Browne. 2000. Identification of a novel genetic locus that is required for in vitro adhesion of a clinical isolate of enterohaemorrhagic Escherichia coli to epithelial cells. Mol. Microbiol. 35:275-288.[CrossRef][Medline]
29 - Oswald, E., H. Schmidt, S. Morabito, H. Karch, O. Marchès, and A. Caprioli. 2000. Typing of intimin genes in human and animal enterohemorrhagic and enteropathogenic Escherichia coli: characterization of a new intimin variant. Infect. Immun. 68:64-71.[Abstract/Free Full Text]
30 - Paton, A. W., and J. C. Paton. 1998. Detection and characterization of Shiga toxigenic Escherichia coli by using multiplex PCR assays for stx1, stx2, eaeA, enterohemorrhagic E. coli hlyA, rfbO111, and rfbO157. J. Clin. Microbiol. 36:598-602.[Abstract/Free Full Text]
31 - Paton, J. C., and A. W. Paton. 1998. Pathogenesis and diagnosis of Shiga toxin-producing Escherichia coli infections. Clin. Microbiol. Rev. 11:450-479.[Abstract/Free Full Text]
32 - Picard-Pasquier, N., M. Ouagued, B. Picard, P. Goullet, and R. Krishnamoorthy. 1989. A simple, sensitive method of analyzing bacterial ribosomal DNA polymorphism. Electrophoresis 10:186-189.[CrossRef][Medline]
33 - Piérard, D., G. Muyldermans, L. Moriau, D. Stevens, and S. Lauwers. 1998. Identification of new verocytotoxin type 2 variant B-subunit genes in human and animal Escherichia coli isolates. J. Clin. Microbiol. 36:3317-3322.[Abstract/Free Full Text]
34 - Piérard, D., L. Van Damme, L. Moriau, D. Stevens, and S. Lauwers. 1997. Virulence factors of verocytotoxin-producing Escherichia coli isolated from raw meats. Appl. Environ. Microbiol. 63:4585-4587.[Abstract]
35 - Pollard, D. R., W. M. Johnson, H. Lior, S. D. Tyler, and K. R. Rozee. 1990. Differentiation of Shiga toxin and verocytotoxin type 1 genes by polymerase chain reaction. J. Infect. Dis. 162:1195-1198.[Medline]
36 - Pradel, N., K. Boukhors, Y. Bertin, C. Forestier, C. Martin, and V. Livrelli. 2001. Heterogeneity of Shiga toxin-producing Escherichia coli strains isolated from hemolytic-uremic syndrome patients, cattle, and food samples in central France. Appl. Environ. Microbiol. 67:2460-2468.[Abstract/Free Full Text]
37 - Pradel, N., V. Livrelli, C. De Champs, J.-B. Palcoux, A. Reynaud, F. Scheutz, J. Sirot, B. Joly, and C. Forestier. 2000. Prevalence and characterization of Shiga toxin-producing Escherichia coli isolated from cattle, food, and children during a one-year prospective study in France. J. Clin. Microbiol. 38:1023-1031.[Abstract/Free Full Text]
38 - Ramachandran, V., M. A. Hornitzky, K. A. Bettelheim, M. J. Walker, and S. P. Djordjevic. 2001. The common ovine Shiga toxin 2-containing Escherichia coli serotypes and human isolates of the same serotypes possess a Stx2d toxin type. J. Clin. Microbiol. 39:1932-1937.[Abstract/Free Full Text]
39 - Rangel, J. M., P. H. Sparling, C. Crowe, P. M. Griffin, and D. L. Swerdlow. 2005. Epidemiology of Escherichia coli O157:H7 outbreaks, United States, 1982-2002. Emerg. Infect. Dis. 11:603-609.[Medline]
40 - Rey, J., S. Sánchez, J. E. Blanco, J. Hermoso de Mendoza, M. Hermoso de Mendoza, A. García, C. Gil, N. Tejero, R. Rubio, and J. M. Alonso. 2006. Prevalence, serotypes and virulence genes of Shiga toxin-producing Escherichia coli isolated from ovine and caprine milk and other dairy products in Spain. Int. J. Food Microbiol. 107:212-217.[CrossRef][Medline]
41 - Schmidt, H., L. Beutin, and H. Karch. 1995. Molecular analysis of the plasmid-encoded hemolysin of Escherichia coli O157:H7 strain EDL 933. Infect. Immun. 63:1055-1061.[Abstract]
42 - Schmidt, H., C. Geitz, P. I. Tarr, M. Frosch, and H. Karch. 1999. Non-O157:H7 pathogenic Shiga toxin-producing Escherichia coli: phenotypic and genetic profiling of virulence traits and evidence for clonality. J. Infect. Dis. 179:115-123.[CrossRef][Medline]
43 - Sperandio, V., J. B. Kaper, M. R. Bortolini, B. C. Neves, R. Keller, and L. R. Trabulsi. 1998. Characterization of the locus of enterocyte effacement (LEE) in different enteropathogenic Escherichia coli (EPEC) and Shiga-toxin producing Escherichia coli (STEC) serotypes. FEMS Microbiol. Lett. 164:133-139.[CrossRef][Medline]
44 - Tyler, S. D., W. M. Johnson, H. Lior, G. Wang, and K. R. Rozee. 1991. Identification of verotoxin type 2 variant B subunit genes in Escherichia coli by the polymerase chain reaction and restriction fragment length polymorphism analysis. J. Clin. Microbiol. 29:1339-1343.[Abstract/Free Full Text]
45 - Unkmeir, A., and H. Schmidt. 2000. Structural analysis of phage-borne stx genes and their flanking sequences in Shiga toxin-producing Escherichia coli and Shigella dysenteriae type 1 strains. Infect. Immun. 68:4856-4864.[Abstract/Free Full Text]
46 - Vernozy-Rozand, C., M. P. Montet, M. Berardin, C. Bavai, and L. Beutin. 2005. Isolation and characterization of Shiga toxin-producing Escherichia coli strains from raw milk cheeses in France. Lett. Appl. Microbiol. 41:235-241.[CrossRef][Medline]
47 - Yamamoto, T., N. Wakisaka, F. Sato, and A. Kato. 1997. Comparison of the nucleotide sequence of enteroaggregative Escherichia coli heat-stable enterotoxin 1 genes among diarrhea-associated Escherichia coli. FEMS Microbiol. Lett. 147:89-95.[CrossRef][Medline]
Applied and Environmental Microbiology, April 2008, p. 2118-2128, Vol. 74, No. 7
0099-2240/08/$08.00+0 doi:10.1128/AEM.02688-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
This article has been cited by other articles:
-
Bielaszewska, M., Stoewe, F., Fruth, A., Zhang, W., Prager, R., Brockmeyer, J., Mellmann, A., Karch, H., Friedrich, A. W.
(2009). Shiga Toxin, Cytolethal Distending Toxin, and Hemolysin Repertoires in Clinical Escherichia coli O91 Isolates. J. Clin. Microbiol.
47: 2061-2066
[Abstract]
[Full Text]
-
Bielaszewska, M., Prager, R., Vandivinit, L., Musken, A., Mellmann, A., Holt, N. J., Tarr, P. I., Karch, H., Zhang, W.
(2009). Detection and Characterization of the Fimbrial sfp Cluster in Enterohemorrhagic Escherichia coli O165:H25/NM Isolates from Humans and Cattle. Appl. Environ. Microbiol.
75: 64-71
[Abstract]
[Full Text]
-
Laing, C., Pegg, C., Yawney, D., Ziebell, K., Steele, M., Johnson, R., Thomas, J. E., Taboada, E. N., Zhang, Y., Gannon, V. P. J.
(2008). Rapid Determination of Escherichia coli O157:H7 Lineage Types and Molecular Subtypes by Using Comparative Genomic Fingerprinting. Appl. Environ. Microbiol.
74: 6606-6615
[Abstract]
[Full Text]