Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2002, p. 1381-1391, Vol. 68, No. 3
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.3.1381-1391.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
BCCM/LMG Bacteria Collection, Laboratory of Microbiology, University of Ghent, B-9000 Ghent,1 Research Group of Industrial Microbiology, Fermentation Technology, and Downstream Processing, Vrije Universiteit Brussel, B-1050 Brussels,7 Department of Immunology, Parasitology, and Ultrastructure, Flemish Interuniversity Institute for Biotechnology, Vrije Universiteit Brussel, B-1640 Sint-Genesius-Rode, Belgium,8 Veneto Agricoltura-Istituto per la Qualità e le Tecnologie Agroalimentari, I-36016 Thiene,2 Dipartimento Scientifico e Tecnologico, Università degli Studi di Verona, I-37134 Verona ,3 DOFATA, University of Catania, I-95124 Catania, Italy,4 Department of Food and Environmental Hygiene, Faculty of Veterinary Medicine, University of Helsinki, SF-00014 Helsinki, Finland,5 Institute for Biotechnology and Molecular Biology, Federal Research Center for Nutrition, D-76131 Karlsruhe, Germany6
Received 9 July 2001/ Accepted 12 December 2001
|
|
|---|
|
|
|---|
Since the late 1970s, E. faecium and E. faecalis have been identified as causes of nosocomial infections, causing an increasing incidence of endocarditis, bacteremia, urinary tract, and neonatal infections (19, 42, 53, 54). Other enterococci have rarely been reported to be pathogenic to humans (15, 33, 35, 39, 44). The ratio of infections due to E. faecalis to those due to all other Enterococcus species was previously approximately 10:1. In recent years, E. faecalis has still dominated as a causative agent in enterococcal bacteremia, but there is a shift which is likely due to the emergence of vancomycin-resistant enterococci, with the predominance of the species E. faecium among this subset of enterococcal isolates (45). It has been demonstrated that enterococci show the ability to take up (27, 28) and transfer antibiotic resistance genes, both vertically (42) and horizontally (36, 48).
On the other hand, these two enterococcal species occur in raw and processed meats, as well as in fermented meat and dairy products, in particular, traditional cheeses in which they may contribute to ripening and product flavor (8, 21, 23, 38, 40, 58). Considering the clinical situation, the safety of E. faecium and E. faecalis strains associated with food fermentations and also their use as probiotics are nowadays questioned (21).
Numerous epidemiological studies involving strain level characterization of enterococci have been done. Many of these studies deal with nosocomial infections and the prevalence of vancomycin-resistant E. faecium in individual or in several hospitals (4, 16, 17, 24, 29, 41, 43) or with the assumed transmission of vancomycin-resistant enterococci from animals to humans (5, 6). However, very few studies have been done to reveal the genomic relationships between vancomycin-resistant and susceptible E. faecium isolates of human, food, and animal origin. Willems et al. (62) and Quednau et al. (51) studied the genomic relationships of antibiotic-resistant E. faecium strains from different sources using amplified fragment length polymorphism (AFLP) analysis and restriction endonuclease analysis (REA) of total chromosomal DNA, respectively. In both studies, host specificity was proposed for strains from, e.g., chicken, pig, and human origin.
In the present study, the intraspecies strain relationships of a large set of susceptible and vancomycin-resistant E. faecium strains from humans, animals, and foods, collected from different European countries, were investigated by different typing methodologies. In total 78 E. faecium isolates were characterized using three genomic typing techniques, random amplified polymorphic DNA (RAPD)-PCR (four different primers), AFLP (two different primer combinations), and pulsed-field gel electrophoresis (PFGE) analysis of SmaI patterns. It is demonstrated that the species E. faecium is composed of two genomic groups which were further divided into, respectively, four and three subclusters. The typing results were evaluated with regard to the origin of the strains, safety aspects such as beta-hemolysis and glycopeptide antibiotic resistance, and bacteriocinogeny.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. E. faecium strains studied
|
RAPD-PCR analysis.
DNA was extracted according to the method of Pitcher et al. (49). The following four random primers were used for DNA amplification: M13 (31), AP4 (4), 1253 (1), and D8635 (1). Amplification, electrophoresis, pattern recognition, and normalization were performed as described previously (2). As in reference 2, two different reaction conditions were applied for primer D8635.
For each strain, the normalized profiles obtained with the four different primers were assembled one after the other into a combined profile using GelCompar version 4.1 software. These combined patterns were imported into the Bionumerics version 1.5 software (Applied Maths) and were analyzed by using the Pearson product moment correlation coefficient and the unweighted pair group with mathematical average clustering algorithm (UPGMA).
AFLP analysis.
Genomic DNA was isolated by using the DNA isolation kit Puregene (Gentra Systems), according to the protocol supplied by the manufacturer. The preparation of templates (using a combination of the restriction enzymes HindIII and MseI) and a preselective PCR amplification was performed as described by Gancheva et al. (25). For selective PCR amplification the following primers were used: M01-ABI (5'GATGAGTCCTGAGTAAA3'), M02-ABI (5'GATGAGTCCTGAGTAAC3'), and the fluorescently labeled primer H01-6FAM (5'GACTGCGTACCAGCTTA3') (Pharmacia Biotech) (selective bases at the 3' end are underlined). PCRs with two primer combinations were performed: H01-6FAM and M01-ABI as well as H01-6FAM and M02-ABI. The selective PCR amplification, sample preparation and separation on a denaturing polyacrylamide gel on an ABI Prism 377 DNA sequencer was as described previously (25).
The data were registered during the electrophoresis run by the ABI Prism 377 data collection software and analyzed with Gene Scan 2.1 software (Applied Biosystems). ABICON software (Applied Maths) was used to import the patterns in GelCompar version 4.2 for normalization. These profiles were transferred to Bionumerics version 1.5 software where patterns generated with both primer combinations were analyzed separately using the Pearson product moment correlation coefficient. The UPGMA cluster analysis was performed using the "composite data set" comparison tool in which the similarity matrices of both individual experiments are averaged.
PFGE analysis.
Cells from 5 ml of an overnight-grown tryptone soy broth (Oxoid) culture were used for DNA isolation. The in situ DNA isolation protocol and digestion of the agarose-embedded DNA were performed as described by Björkroth et al. (7). SmaI restriction endonuclease was selected (Promega) due to its good ability to characterize Enterococcus DNA in PFGE experiments (41, 47). The samples were electrophoresed through 1.0% (wt/vol) agarose gel (SeaKem Gold; FMC BioProducts) in 0.5x TBE (45 mM Tris, 4.5 mM boric acid [pH 8.3], and 1 mM sodium EDTA) at 14°C in a CHEF-DR III system (Bio-Rad). Interpolation ramping from an initial pulse time of 0.1 s to a final 25 s for 18 h at 6 V/cm with a 120° angle was used.
Photographs of the PFGE banding patterns were scanned by Hewlett-Packard ScanJet 4c/T scanner. Numerical analysis of macrorestriction patterns was performed using the GelCompar system version 4.2 (Applied Maths). Based on the banding patterns of internal controls, 1% position tolerance was allowed between the strains. The similarity between all pairs was expressed by a Dice coefficient correlation, and UPGMA clustering was used for the construction of the dendrogram.
DNA-DNA hybridization experiments.
DNA extraction, determination of the G+C content, and DNA-DNA hybridizations were performed as described previously (60).
Phenotypic identification.
The isolates were phenotypically identified with the rapid ID 32 STREP system according to the protocol specified by the manufacturer (bioMérieux).
Hemolysis activity.
Production of hemolysin was determined by streaking the Enterococcus cultures (grown in MRS broth for 18 h) on both human and sheep blood agar plates. Human blood agar was prepared using Columbia blood agar base (BBL) with 5% defibrinated human blood (containing all four blood types) added according to the manufacturer's instructions. Sheep blood agar was commercially available (BBL). Blood agar plates were incubated at 37°C for 24 h, after which plates were examined for hemolysis. Presence and absence of zones of clearing around the colonies were interpreted as beta-hemolysis and gamma-hemolysis activity, respectively.
Detection of glycopeptide resistance.
Susceptibilities to vancomycin and teicoplanin were checked using the broth dilution technique. Strains were inoculated in Trypticase soy broth (Oxoid) with 3 g of yeast extract liter-1 and with either vancomycin or teicoplanin at a concentration ranging from 2 to 64 µg ml-1. After an overnight incubation at 37°C, growth was noted and the MIC (the first dilution at which no growth was observed) was determined. The presence of van genes was checked using the PCR as described by Dutka-Malen et al. (19). Four different primer pairs (vanA, vanB, vanC1, and vanC2) were used in the assay.
Antimicrobial testing and detection of bacteriocin structural genes.
Sixteen-hour cultures of enterococci were checked for inhibitory activity by the agar well diffusion assay towards Clostridium tyrobutyricum LMG 1285T or the agar spot test towards Listeria innocua LMG 11387T and LMG 13568 and Propionibacterium freudenreichii subsp. shermanii LMG 16424T as described in detail by Andrighetto et al. (2). Strains showing a positive reaction in one of the above-mentioned tests were further genotypically investigated by PCR for the presence of structural genes for the bacteriocins enterocin A (3), enterocin B (3), enterocin P (10), and bacteriocin 31 (57). The specific primers used for PCR amplification of these genes are listed in Table 2. PCR amplification was done using total DNA isolated by the rapid alkaline lysis method as described in detail by Dutka-Malen et al. (18). In a first series of PCRs, 2 µl of total DNA was used as a template. Two microliters of the resulting PCR amplifications were used as a template in a second series of nested PCRs. PCR was performed on a DNA thermal cycler (Biometra Thermocycler T3) in a final volume of 50 µl containing 1x PCR buffer (20 mM Tris-HCl [pH 8.4], 50 mM KCl), 1.5 mM MgCl2, a 200 µM concentration of each of the four deoxynucleoside triphosphates, a 0.5 µM concentration of each of the primers, and 1.25 U of Taq DNA polymerase (GibcoBRL). The cycles used were 95°C for 5 min for the first cycle; 95°C for 30 s, 58 or 56°C for 30 s, and 72°C for 30 s for the next 30 cycles; and 72°C for 5 min for the last cycle. PCR products were resolved by electrophoresis on 8% polyacrylamide-Tris-acetate-EDTA gel. The total DNA of E. faecium CTC 492 (encoding the enterocins A and B [3]) and plasmid DNA of E. faecalis FA2-2 carrying the plasmid pYI171 (encoding bacteriocin 31 [57]) was used as a positive control for the corresponding primers. No control was available for enterocin P; therefore, PCR amplicons were sequenced with an automated DNA sequencer (ALF DNA sequencer; Amersham) to confirm specific enterocin P sequences.
|
View this table: [in a new window] |
TABLE 2. Specific terminal and nested primers for the PCR detection of the enterocin structural genes
|
|
|
|---|
RAPD-PCR analysis.
The reproducibility of PCR assays and running conditions estimated by analysis of duplicate DNA extracts of several strains was 0.91. Purified DNA of all E. faecium strains was used as template in five amplification reactions with single primers. To increase the accuracy and discriminative power of the method, the electrophoretic patterns were combined for numerical analysis. Figure 1 (left dendrogram) shows the UPGMA tree. At a correlation level of 0.34, two main groups were distinguished, indicated as group I and II. Each group was further subdivided into, respectively, four (R1 to R4) and three (R5 to R7) subclusters. A photograph of the digitized band patterns obtained with the different PCR assays for representative strains of each of the subclusters is given in Fig. 2. The differences in the band patterns of the two main groups (I and II) were substantial.
![]() View larger version (44K): [in a new window] |
FIG. 1. Dendrograms derived from UPGMA linkage of r values (expressed as percentages) of RAPD-PCR (left side) and AFLP (right side) profiles. Strains in boldface type group in the corresponding subcluster in both analyses.
|
![]() View larger version (44K): [in a new window] |
FIG. 2. Digitized and normalized patterns of selected strains obtained with AFLP fingerprinting with two primer combinations: H01-6FAM and M01-ABI (A) and H01-6FAM and M02-ABI (B) and RAPD-PCR fingerprinting using four primersM13 (C), D8635 (D), AP4 (E), and 1253 (F)and primer D8635 (G) (primer used with two different reaction conditions).
|
PFGE analysis.
The reproducibility between different restriction digests and electrophoretic runs for genomes of individual strains was 0.93. The number of bands per strain varied between 16 and 26. Cluster analysis and visual inspection of the patterns (Fig. 3) showed a high variability among the patterns obtained. High similarity was observed only between a few pairs of isolates (see discussion).
![]() View larger version (45K): [in a new window] |
FIG. 3. Dendrogram derived from UPGMA linkage of r values (expressed as percentages) and normalized banding patterns of PFGE fingerprints.
|
|
View this table: [in a new window] |
TABLE 3. Percentage DNA-DNA binding by reassociation studies and G+C mol% of E. faecium strains
|
|
View this table: [in a new window] |
TABLE 4. Phenotypic characterization of the E. faecium strains studieda
|
Detection of glycopeptide resistance.
Using the broth dilution technique, MICs of vancomycin and teicoplanin were determined. Subsequently, a screening was performed for the presence of van genes and confirmed that all strains for which the MIC of vancomycin was in a range of 32 to >64 µg ml-1 and the MIC of teicoplanin was in a range of 4 to 64 µg ml-1 contain the vanA gene (Table 1). For all other strains, the MIC of both antibiotics was 2 µg ml-1 and no van genes were detected.
Antimicrobial testing.
The presence of one or more enterocin genes (A, B, P, and 31) in strains showing antimicrobial activity towards one of the four indicator strains used (see Materials and Methods) is shown in Table 1. For six strains showing inhibition on plates, none of the above-mentioned genes were detected and for 34 strains one or more of the enterocins A, B, P, and 31 were present.
|
|
|---|
Diversity of the strain collection.
Strains originated from Belgium, Germany, Greece, Ireland, Italy, and The Netherlands, and a few strains were from other European and non-European countries. The study not only included human (clin-ical and nonclinical isolates) and animal (isolates from at least 10 different animals) strains, but also strains originating from food (mainly isolates from cheese) (Table 1). PFGE (using SmaI), which is currently considered to be the "gold standard" for strain differentiation of enterococci (4, 26, 30, 34, 56), was applied on the majority of the strains to verify whether the isolates represent different strains or may be considered clones of the same parent strain. Although some uncertainty remains as to how much similarity isolates should demonstrate in order to be called clonally related (56), we observed a significant diversity among most isolates in order to consider them to originate from different clonal lineages (Fig. 3). Considering the very diverse origin and habitat of the selected cultures, this strain diversity was expected and also confirmed. Numerical and visual comparison of PFGE patterns (Fig. 3) suggests that only the following isolates are possibly clonal lineages: two isolates from human bile in The Netherlands (FAIR-E 134 and FAIR-E 137), four isolates from feta cheese in Greece (FAIR-E 198, FAIR-E 201, FAIR-E 202, and FAIR-E 206), and two isolates from Italian cheeses (FAIR-E 362 and FAIR-E 366).
Delineation of intraspecies genomic groups.
For detecting the intraspecies strain relationships, we applied RAPD-PCR and AFLP, both with multiple primers or primer combinations in order to enlarge the number of bands and to avoid biasing of the numerical analysis due to the presence of too few bands. Although both approaches analyze entire genomes, a crucial difference between the two techniques is that RAPD-PCR is prone to variations in amplification efficiency due to small variations in experimental conditions or the quality of chemicals which may affect the reproducibility of the method, whereas AFLP primers perfectly match their target site and mismatches are not expected to participate in the amplification process (32).
The very high congruency between the cluster analyses obtained by both methods is remarkable (Fig. 1). In both analyses, two main groups, groups I and II, were delineated, comprising 62 and 16 strains, respectively. Within genomic group I, four subclusters were recognized (R1 to R4 and A1 to A4, respectively). A high correlation was found in subclusters R1 and A1 and in R2 and A2, respectively. Most of the strains from the clusters R3 and R4 were found in the AFLP clusters A3 and A4, but not necessarily in the corresponding cluster number. Only six strains assigned to a particular subcluster in the RAPD-PCR dendrogram occupied separate positions in the AFLP analysis. In group II, three identical subclusters with exactly the same set of strains were found in the RAPD-PCR and AFLP cluster analysis (R5 to R7 and A5 to A7, respectively).
No literature data have thus far suggested the presence of two genomic groups in E. faecium, which were also easily recognized after analyses with single primers or primer combinations and not only after combined numerical analyses (Fig. 2; data not shown). Therefore, in a first approach, we investigated whether these two groups represent different genomic species or not. Respectively, four and three strains were selected from different subclusters of both groups (including the type strain), and DNA-DNA hybridizations were performed. A mean DNA homology of 71% was observed between the main groups I and II, whereas in both groups the level of DNA relatedness was significantly higher and showed a mean value of 82% (Table 3). Although there is a clear indication of a lower overall DNA homology between the two groups than within the individual groups, the data confirm that all strains constitute a single genomic species on the basis of the recommendations proposed by Wayne et al. (61). With the rapid ID 32 STREP strip, no phenotypic features were found which clearly differentiate the two genomic groups. Among the many variable features (Table 4), we only noticed different percentages of strains that are positive for particular tests.
Correlation with origin, pathogenicity, and bacteriocinogeny.
Interestingly, comparison of the origin of the strains of the different groups and subclusters revealed that the main groups I and II both comprise strains from (healthy) humans, animals, and food and that no correlation is found with the country of isolation. However, all human clinical strains clustered in group I. Regarding the source of the strains per subcluster in group I, we noticed that all animal isolates and human clinical strains grouped in a single subcluster (R1 or A1), whereas nearly all the food strains grouped in three clusters (R2 to R4 or A2 to A4). Isolates from healthy humans are spread over the different subclusters. In group II, subclusters R5 and R6 contain strains from different countries from food, animal, and human (only R5) origin. Subcluster R7 contains only two Italian food isolates, for which a clonal relationship is suggested (Fig. 3). Considering the low number of strains in group II, conclusions concerning origin in the delineated subclusters are, however, premature.
The genotypic grouping was also compared with possible pathogenic traits of the strains. All resistant strains (containing the vanA gene) clustered in group I, and the majority of them belonged to subcluster R1 and/or A1 (a separate position in one of the analyses was observed for strains FAIR-E 41, FAIR-E 102, FAIR-E 157, and FAIR-E 158). FAIR-E 84 is the only vanA strain that grouped inconsistently by RAPD-PCR and AFLP with regard to its position in subclusters R3 and A1, respectively. All strains were also tested for their hemolytic activity on human and sheep blood. Only five strains showed beta-hemolytic activity on at least one of the blood types and all clustered in genomic group I. One strain, FAIR-E 23, isolated from ostrich in South Africa, was beta-hemolytic on both blood types and grouped in subcluster R1 or occupied a separate position as seen by AFLP analysis. Four strains, all isolated from feta cheese in Greece, were hemolytic only on human blood. They all grouped in subclusters R3 and A4, respectively. The isolation of hemolytic enterococci from cheese may have sanitary implications and also illustrates the notion that virulence factors may occur in enterococci originating from food, as has previously been demonstrated (20, 22).
About half of the strains displayed antagonistic activity towards food-spoiling bacteria or related taxa (C. tyrobutyricum, L. innocua, or P. freudenreichii), and for most of them the presence of one or more of the enterocin genes A, B, P, and 31 was confirmed. All enterocin B-containing strains are found in group I, and with the exception of one strain (FAIR-E 119), all belonged to subcluster R1. Six of the eight strains for which inhibitory activity was seen on plates and for which no known genes were detected belonged to genomic group II. In the latter group, only the two strains of subcluster R7 (A7) contain known but different bacteriocin genes.
When comparing data from the present study with those presented by Willems et al. (62) and Quednau et al. (51), host specificity for human and animal strains (although for most animals only one or a few strains were investigated) is not confirmed by our results. In the study of Willems et al. (62) host specificity of vancomycin-resistant E. faecium isolates was suggested for pigs, chickens, calves, and humans. For the humans, the authors furthermore distinguished between strains from nonhospitalized and hospitalized individuals. Although there is a major difference to our study, in that the majority of the strains are isolated from food and are susceptible to vancomycin, the use of the term "host specificity" in the study of Willems et al. (62) may be questioned, given that farmers and slaughterers (also nonhospitalized individuals) contain the same flora as the animals they are working with. Furthermore, the restricted geographic origin of most of their isolates (mostly from The Netherlands) and presumed epidemiological relatedness among isolates (demonstrated for part of the hospitalized persons; other hosts not studied) are points of discussion. In this context, and as suggested in the discussion of Willems et al. (62), the term "ecovar specificity" would be more appropriate until a more diverse set of isolates is studied. The separate cluster obtained for strains from hospitalized persons nevertheless remains remarkable, although we should bear in mind that these vancomycin-resistant strains may be endemic in the hospital(s) and should not necessarily be seen as part of the natural flora of the patients as suggested in the work of Willems et al. (62).
Host specificity of E. faecium isolates from chicken, pork, and humans is also suggested by Quednau et al. (51), but contrary to the study of Willems et al. (62) no subdivision was obtained according to clinical or healthy-subject strains. The authors delineate multiple clusters for each host, often consisting of only a few isolates from the same and restricted geographic origin. This could be an indication that REA, used in the study, is measuring clonal relationships and not relatedness between epidemiologically unrelated strains. We presume that, as for PFGE analysis, the deeper branches in the REA cluster analysis might not provide valuable information on the genomic relationships among epidemiologically nonrelated enterococcal strains. This cannot, however, be verified as no second typing approach (using an alternative technique) was applied by Quednau et al. (51) on the same set of strains to confirm the REA grouping.
Conclusion.
Data from the present study demonstrate that two main genomic groups, groups I and II, are recognized in the species E. faecium. No host specificity or correlation with the country of isolation is found, but the investigated human clinical strains, antibiotic-resistant strains, and beta-hemolytic strains are found only in genomic group I. In this genomic group, animal and human clinical strains grouped in a single subcluster (including the majority of vancomycin-resistant and enterocin B-producing strains), while nearly all strains originating from food grouped in the other subclusters. Human strains from healthy individuals are found in genomic groups I and II.
Should a separate grouping of more extended sets of potentially pathogenic strains be confirmed in future studies, the observed genomic heterogeneity within the species E. faecium may be useful for the selection of safe strains for starter cultures in food or probiotic preparations.
We thank Karen Lefebvre for her skilled technical assistance.
|
|
|---|
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»